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

Treestech · 发表于 2012-12-7 21:07:25

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

%本人系电压稳定菜鸟，对电压稳定分析中基础知识左特征根向量右特征根向量，模态，参与因子的概念通过一个具体的矩阵进行了MATLAB演算 %希望对和我一样的初学者有点直观的认识
A=[1 0 1;0 1 1; 0 0 2];
[X,Y]=eig(A);
[m,n]=eig(A');
%右特征根
x10=1;
x1=[1 0 0]';
x20=1;
x2=[0 1 0]';
x30=2;
x3=[0.57735 0.57735 0.57735]';
%左特征根

y10=1;
y1=[0.70711 0 -0.70711]';
y20=1;
y2=[0 0.70711 -0.70711]';
y30=2;
y3=[0 0 1]';

MM=[x1,x2,x3];
NN=[y1,y2,y3];
aa=[1 0 0;0 1 0;0 0 2];
%归一化处理
k1=y1'*x1;
k2=y2'*x2;
k3=y3'*x3;
y1p=y1/k1;
y2p=y2/k2;
y3p=y3/k3;
NNP=[y1p,y2p,y3p]';

%此时矩阵A=MM(右特征向量矩阵)*aa（对角阵）*NNP（左特征向量矩阵）
%A=(x1,x2,x3)*[lemda1 0 0;0 lemda2 0;0 0 lemda3]*(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([lemda1 0 0;0 lemda2 0;0 0 lemda3])*左特征矩阵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=[0 0 1]';
Q=[0 0 1]';
cc=inv(A)*Q
dd=1*y1p(3)*x1+1*y2p(3)*x2+0.5*y3p(3)*x3;
cctest=cc-dd;
%而dUk/dQk:
Qa=[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 · 发表于 2012-12-8 00:37:29

感谢楼主的资料

liliang0603 · 发表于 2012-12-12 16:33:29

谢谢，很有用

11216788 · 发表于 2013-12-25 18:35:14

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

redplum · 发表于 2013-12-26 09:00:09

这程序很好嘛

周沭铭 · 发表于 2013-12-26 09:07:23

谢谢楼主分享

gczhu · 发表于 2014-1-3 10:52:26

学习了

interzsm · 发表于 2014-1-25 18:01:58

感谢楼主的资料分享

redplum · 发表于 2014-1-30 01:45:34

不用费流量啊阿，好

电网工人 · 发表于 2014-1-30 08:51:11

账号		自动登录	找回密码
密码			立即加入

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

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

评分

相关帖子