Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am solving a system of equations derived from power system applications. It consists of index-1 differential and algebraic equations in the form: $$\dot{x}=f(x,y) \\ 0=g(x,y)$$

To get the eigenvalues of the system I eliminate to get the reduced jacobian matrix $$\Delta \dot{x}=(f_x-f_y*g_y^{-1}*g_x)\Delta x$$

From the eigenvalue analysis I get the eigenvalues that were expected by theory. But, I also get a pair of positive eigenvalues. $0.3\pm3.5i$. The participation factors point to a PI controller and by tuning the Kp gain I can eliminate them (move to the left plane).

My question is: Although I have these positive eigenvalues, they don't show up in time domain simulation (solving the system using trapezoidal method) over discretized time horizon. I make pulse disturbances and the system appears to be stable while I would expect to observe instability because of the positive eigenvalues.

Could this be explained or should I assume that there is something wrong with the simulation software or the eigenvalue analysis?

Thanks in advance!

share|cite|improve this question

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.