电压稳定模态分析基础MATLAB实例

Treestech

%本人系电压稳定菜鸟，对电压稳定分析中基础知识左特征根向量右特征根向量，模态，参与因子的概念通过一个具体的矩阵进行了MATLAB演算 %希望对和我一样的初学者有点直观的认识
  b0 K- n4 F" l# f2 [A=[1 0 1;0 1 1; 0 0 2];
[X,Y]=eig(A);
[m,n]=eig(A');
& H4 \& g& g. A* Y%右特征根
" P$ v% w% k& j" N) ix10=1;
- p) c' B! n. \: `  cx1=[1 0 0]';
. x+ F1 W! u, ~7 X) px20=1;
x2=[0 1 0]';
9 Q0 r2 o6 z& T& xx30=2;
% }2 `4 ^$ c% I4 Z9 j8 K( Rx3=[0.57735 0.57735 0.57735]';
' g4 N+ D5 |2 F+ v2 l6 |. Z" |; ]%左特征根
% x! ^$ v9 Y' J9 e! U/ X1 ~
& x1 b, s9 E+ V6 x& M8 qy10=1;
y1=[0.70711 0 -0.70711]';
: \# v$ K3 |5 Q: Q; w  ^, Jy20=1;
y2=[0 0.70711 -0.70711]';
/ Z& S1 p9 c7 d& iy30=2;
y3=[0 0 1]';

0 F5 h3 `) k7 nMM=[x1,x2,x3];
NN=[y1,y2,y3];
aa=[1 0 0;0 1 0;0 0 2];
5 t) R9 ~/ m. Q& Q& Q! m%归一化处理
: ^$ W0 t: X  y1 R" K. jk1=y1'*x1;
k2=y2'*x2;
6 I+ n' V$ J! R0 w/ f& q3 T4 T, q0 Ck3=y3'*x3;
% w3 D3 x  C# D. g( M8 Ny1p=y1/k1;
y2p=y2/k2;
- g1 V- U- a% |1 uy3p=y3/k3;
8 Y" w# |$ F: N7 `. }2 q$ f- gNNP=[y1p,y2p,y3p]';
7 c. r0 i0 [4 o# ?8 t! _5 y) ^
9 @9 k7 D- N/ m) i: w%此时矩阵A=MM(右特征向量矩阵)*aa（对角阵）*NNP（左特征向量矩阵）
%A=(x1,x2,x3)*[lemda1 0 0;0 lemda2 0;0 0 lemda3]*(y1p';y2p';y3p');
5 {9 g/ v& F& r- RtestA=1*x1*y1p'+1*x2*y2p'+1/2*x3*y3p';  %将一个矩阵的逆写成各个特征值倒数乘以相对应的右特征列向量和左特征行向量相乘的和。
dtest=testA-inv(A);
( y6 J% W8 b  [$ ^& ~) J7 C
%模态分析法分析电压稳定时  deta V= inv(J)* deta Q 而雅克比矩阵逆可以写成 左特征矩阵NNP *deta V=
: h$ P' ]0 |5 F5 n. w7 P2 _# P%inv([lemda1 0 0;0 lemda2 0;0 0 lemda3])*左特征矩阵NNP*deta Q 而将左特征向量矩阵NNP*detaV
/ ^: _  b& D& r% ?%称为模态电压变化向量，而将左特征向量NNP*detaQ 称作是模态无功变化向量  假设以上矩阵A 就是雅克比矩阵，即存在 deta
%V=MM*inv(aa)*NNP* deta Q 即可以验证 NNP*det V=inv(aa)*NNP*deta Q
: q  D! S3 T/ Q& i0 ^%即电压模态变化向量=inv(aa)*无功变化向量
. X- t4 a3 M  T1 C%特征值越小，说明电压稳定性越差，该模态电压降崩溃，本矩阵A特征值为1，1，2，则第三个模态下相比前连个模态电压更稳定。
3 {1 P( r3 I: c9 H; @# G
%假设 Deta Q=[0 0 1]';
, F: Y+ t0 K7 U  ~1 D* OQ=[0 0 1]';
; H( }. p$ \1 P, l, P* P: qcc=inv(A)*Q
dd=1*y1p(3)*x1+1*y2p(3)*x2+0.5*y3p(3)*x3;
* T0 H" W7 J3 w0 }cctest=cc-dd;
) r; e" J8 v, V+ {, D1 \7 g%而dUk/dQk:
6 k, W2 F0 V' i- V. {* h5 [3 F& rQa=[0 0 1]';
cca=inv(A)*Qa
dda=1*y1p(3)*x1(3)+1*y2p(3)*x2(3)+0.5*y3p(3)*x3(3);
ccatest=cca(3)-dda;

%y1p(3)*x1(3) 即为参与度因子

redplum

感谢楼主的资料

liliang0603

谢谢，很有用

11216788

正在研究特征值计算，学习了。

redplum

这程序很好嘛

周沭铭

谢谢  楼主分享

gczhu

学习了

interzsm

感谢楼主的资料分享

redplum

不用费流量啊阿，好

电网工人