Hello everyone I have been using MATPOWER for quite some time now. My work has to do with voltage stability in power systems and hence I find myself frequently using the AC power flow in MATPOWER. Recently I have been trying to perform Modal Analysis for voltage stability monitoring. Now Modal Analysis requires me to calculate the Eigen values and eigen vectors of the reduced Jacobian or full Jacobian matrix. Now in my analysis I increase the load at all the PQ buses of a network incrementally while performing power flow (using runpf) at each load increase. Now I need to extract the full jacobian everytime I run the power flow. Once I have the jacobian I can then go ahead and study the eigen values and eigen vectors.
Now my question is how can i extract the full Jacobian every time I run the power flow (using runpf). Here is a small part of my code that I have been using: results=runpf(this_case,mpopt); J=full(makeJac(this_case)); Doing this gives me a jacobian matrix but I am running this in a loop where I update the case file (this_case) by changing the loads on all the buses and running the power flow again. So for every new case file that I am using, the jacobian matrix should be different but in my case it doesn't change. I would appreciate any help in this regard. Thanks -- *Aaqib A Peerzada* *Graduate Student * *Department of Electrical and Computer Engineering* *Wisenbaker Engineering Research Center* *Texas A&M University* *College Station, Texas* *United States.* *Email: [email protected] <[email protected]>*
