电压稳定模态分析基础MATLAB实例
%本人系电压稳定菜鸟，对电压稳定分析中基础知识左特征根向量右特征根向量，模态，参与因子的概念通过一个具体的矩阵进行了MATLAB演算 %希望对和我一样的初学者有点直观的认识A=;
=eig(A);
=eig(A');
%右特征根
x10=1;
x1=';
x20=1;
x2=';
x30=2;
x3=';
%左特征根
y10=1;
y1=';
y20=1;
y2=';
y30=2;
y3=';
MM=;
NN=;
aa=;
%归一化处理
k1=y1'*x1;
k2=y2'*x2;
k3=y3'*x3;
y1p=y1/k1;
y2p=y2/k2;
y3p=y3/k3;
NNP=';
%此时矩阵A=MM(右特征向量矩阵)*aa（对角阵）*NNP（左特征向量矩阵）
%A=(x1,x2,x3)**(y1p';y2p';y3p');
testA=1*x1*y1p'+1*x2*y2p'+1/2*x3*y3p';%将一个矩阵的逆写成各个特征值倒数乘以相对应的右特征列向量和左特征行向量相乘的和。
dtest=testA-inv(A);
%模态分析法分析电压稳定时deta V= inv(J)* deta Q 而雅克比矩阵逆可以写成 左特征矩阵NNP *deta V=
%inv()*左特征矩阵NNP*deta Q 而将左特征向量矩阵NNP*detaV
%称为模态电压变化向量，而将左特征向量NNP*detaQ 称作是模态无功变化向量假设以上矩阵A 就是雅克比矩阵，即存在 deta
%V=MM*inv(aa)*NNP* deta Q 即可以验证 NNP*det V=inv(aa)*NNP*deta Q
%即电压模态变化向量=inv(aa)*无功变化向量
%特征值越小，说明电压稳定性越差，该模态电压降崩溃，本矩阵A特征值为1，1，2，则第三个模态下相比前连个模态电压更稳定。
%假设 Deta Q=';
Q=';
cc=inv(A)*Q
dd=1*y1p(3)*x1+1*y2p(3)*x2+0.5*y3p(3)*x3;
cctest=cc-dd;
%而dUk/dQk:
Qa=';
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) 即为参与度因子 感谢楼主的资料 谢谢，很有用 正在研究特征值计算，学习了。 这程序很好嘛 谢谢楼主分享 学习了 感谢楼主的资料分享 不用费流量啊阿，好 ::handshake::
页:
[1]
2
