用MATLAB形成节点导纳矩阵(电力系统稳态分析)-Toy模板网

这篇具有很好参考价值的文章主要介绍了用MATLAB形成节点导纳矩阵(电力系统稳态分析)。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

前言

本文通过一个简单电力网络图来说明计算机编写形成节点导纳矩阵的具体方法。

一、看懂电力网络图及确定原始数据

1.1 所求电力网络图

matlab计算导纳矩阵,潮流计算,电力系统分析,matlab,矩阵,学习

1.2确定形成节点导纳矩阵的原始数据

（1）节点数：n=5；
（2）支路数：nl=5；
（3）支路参数矩阵B：包括六个数据[i,j,z,b,t,it],i,j为支路两端节点号,z为支路的阻抗，b为线路电纳,t为变比,it为高低压侧标志（高为1，低为0）。
这里有5条支路故有5组数据

B=[1,2,0.03i,0,1.05,0;2,3,0.08i+0.3i,0.5i,1,1;2,5,0.1+0.35i,0,1,1;3,5,0.04+0.25i,0.5i,1,1;3,4,0.015i,0,1.05,0]

（4）节点对地阻抗矩阵X（由节点号与接地阻抗构成）。图中节点均未接地，没有对地阻抗故为零

X=[1 0;2 0;3 0;4 0;5 0;]

前面的参数为节点编号，后面为对地阻抗。

二、M程序的编写

n=input('n=');
nl=input('nl=');
B=input('B=');
X=input('X=');
Y=zeros(n);
for i=1:n
    if X(i,2)~=0;
        p=X(i,1);
        Y(p,p)=1./X(i,2);
    end
end
for i=1:nl
    if B(i,6)==0
        p=B(i,1);q=B(i,2);
    else p=B(i,2);q=B(i,1);
    end
    Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));
    Y(q,p)=Y(p,q);
    Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)*2)+B(i,4)./2;
    Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;
end
disp('Y=');
disp(Y)