Test ensemble stability

ensembleStabilityTest(ensemble, eigThreshold, strucIdx)

Takes in a model ensemble, calculates the jacobian and respective eigenvalues for each model. If any eigenvalue’s real part is higher than eighThreshold, the model is considered unstable.

USAGE:

stabilityRes = ensembleStabilityTest(ensemble, eigThreshold, strucIdx)

INPUT:

ensemble (struct): model ensemble, see buildEnsemble for fields description eigThreshold (double): threshold for positive eigenvalues’ real part

OPTIONAL INPUT:

strucIdx (int): number of the model structure considered

OUTPUT:

stabilityRes (struct): stability analysis results

• posEig (cell) : positive eigenvalues for unstable models

• unstableModels(int vector) : list of unstable models