马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
%本人系电压稳定菜鸟,对电压稳定分析中基础知识左特征根向量右特征根向量,模态,参与因子的概念通过一个具体的矩阵进行了MATLAB演算 %希望对和我一样的初学者有点直观的认识
% P, c4 t5 V4 ]0 \A=[1 0 1;0 1 1; 0 0 2];+ r0 n( S: \! F8 J
[X,Y]=eig(A);. f" |1 b; |6 k9 s
[m,n]=eig(A');
4 C0 p* H9 U Z* T W4 N1 @%右特征根, ~( D/ Z/ v% R7 ]8 X
x10=1;
' \& R6 r$ M1 G( Yx1=[1 0 0]';! _5 f% X* @: @; }+ e6 H
x20=1;& R0 H' }$ m0 { {
x2=[0 1 0]';
" Z9 ]1 L( F/ x2 Hx30=2;
) W5 U. a$ c5 |/ C' xx3=[0.57735 0.57735 0.57735]';
' t6 @' j; x3 U, G' K1 b Y. v%左特征根8 S% a1 s7 p2 l/ V8 B$ M5 d' |
e9 F# @6 z5 k! k7 xy10=1;7 M" d( @& U3 y% U% u; T. U
y1=[0.70711 0 -0.70711]';6 p9 s! k$ U6 T: x; r
y20=1;
. ]. P6 ~. a5 `- n% Ry2=[0 0.70711 -0.70711]';
, d0 O1 ^. n" hy30=2;+ ]8 K! ]7 R2 m6 N/ ^0 u
y3=[0 0 1]';
$ W9 R/ P( p1 K4 C
" Q( J- x' X% R i5 h; r" N1 LMM=[x1,x2,x3];
( y! p3 a, P A$ oNN=[y1,y2,y3];
+ E$ N( v' A2 }/ s: yaa=[1 0 0;0 1 0;0 0 2];8 L, i2 W0 ?; ?0 v/ Y V+ g
%归一化处理! U3 F6 K1 j$ @/ `9 P
k1=y1'*x1;
; I' g' _, J7 Z `! p5 d- J" _- kk2=y2'*x2; }4 w5 a' [6 | _
k3=y3'*x3;& @) ]$ y3 t' J3 f
y1p=y1/k1;
\+ g& N) g8 f0 j, J; x0 K$ Qy2p=y2/k2;
3 ?% h# N. x5 h* B' [) sy3p=y3/k3;
; s6 s+ N7 {3 Z) Z, H1 R+ fNNP=[y1p,y2p,y3p]';
9 }! Y' c# k6 ]. J* x* r$ z5 V2 Z4 X0 N$ u) ], P; v$ d
%此时矩阵A=MM(右特征向量矩阵)*aa(对角阵)*NNP(左特征向量矩阵)
' b; j: B9 v+ M$ p# b%A=(x1,x2,x3)*[lemda1 0 0;0 lemda2 0;0 0 lemda3]*(y1p';y2p';y3p');
! t5 {: S( T2 ~7 b2 B$ S3 P- \testA=1*x1*y1p'+1*x2*y2p'+1/2*x3*y3p'; %将一个矩阵的逆写成各个特征值倒数乘以相对应的右特征列向量和左特征行向量相乘的和。6 g& c( t# v9 X/ y! |
dtest=testA-inv(A);, \, K2 ^' n, K1 J: f) n: w
, f8 r6 [0 e6 ^. _$ Z
%模态分析法分析电压稳定时 deta V= inv(J)* deta Q 而雅克比矩阵逆可以写成 左特征矩阵NNP *deta V=
2 }3 x" K8 J* A%inv([lemda1 0 0;0 lemda2 0;0 0 lemda3])*左特征矩阵NNP*deta Q 而将左特征向量矩阵NNP*detaV
& u2 n( L: e6 K% e/ t%称为模态电压变化向量,而将左特征向量NNP*detaQ 称作是模态无功变化向量 假设以上矩阵A 就是雅克比矩阵,即存在 deta) r& s8 T. R6 p1 A
%V=MM*inv(aa)*NNP* deta Q 即可以验证 NNP*det V=inv(aa)*NNP*deta Q' `, |7 [/ N* ]5 V. T
%即电压模态变化向量=inv(aa)*无功变化向量( Q1 E6 u6 K# W3 W! o
%特征值越小,说明电压稳定性越差,该模态电压降崩溃,本矩阵A特征值为1,1,2,则第三个模态下相比前连个模态电压更稳定。
) `: o' c0 Q/ Y% o5 E2 B
) N4 {' s0 C6 F) V) K%假设 Deta Q=[0 0 1]';( j( N% X2 K: ?; Q( P1 O
Q=[0 0 1]';! K- y4 E- b- Y! E
cc=inv(A)*Q
4 v& x$ k6 q' n/ a1 Jdd=1*y1p(3)*x1+1*y2p(3)*x2+0.5*y3p(3)*x3;. V7 F0 t0 d' I/ L; V0 c- n4 v
cctest=cc-dd;
+ g, X/ U' I0 F) t9 I; I, \%而dUk/dQk:/ p4 a( W; v- H8 P! Y) v
Qa=[0 0 1]';
0 R$ _& P: \+ \+ L5 l! {cca=inv(A)*Qa
' h U7 D0 E7 idda=1*y1p(3)*x1(3)+1*y2p(3)*x2(3)+0.5*y3p(3)*x3(3);' J" a$ p' E* l% \! b/ b% E
ccatest=cca(3)-dda;
' V: F$ ~8 @4 r) r A/ E- @& w! V" ^
%y1p(3)*x1(3) 即为参与度因子 |