节点导纳矩阵程序代码

程序代码

n=input('请输入节点数：n=');

L=input('请输入支路数：L=');

W=input('请输入支路参数标幺値：W=');

m=input('请输入发电机数：m=');

M=input('请输入与发电机连接的节点号及其阻抗标幺値：M='); F=input('请输入短路节点号：F=');

Y=zeros(n);

for i=1:m

s=M(i,1);

Y(s,s)=1./M(i,2);

end

for i=1:L

if W(i,6)~=0

s=W(i,2);g=W(i,1);

else

s=W(i,1);g=W(i,2);

end

Y(s,g)=Y(s,g)-1./(W(i,3)*W(i,5));

Y(g,s)=Y(s,g);

Y(g,g)=Y(g,g)+1./(W(i,3)*(W(i,5)^2))+W(i,4); Y(s,s)=Y(s,s)+1./W(i,3)+W(i,4);

end

disp(Y);

Yb=Y;

Yb(F,:)=[];

Yb(:,F)=[];

disp(Yb);

Ya=Y;

Ya(3,3)=Y(3,3)+1/(W(4,3)+i*W(4,4))-W(4,5);

Ya(4,4)=Y(4,4)+1/(W(4,3)+i*W(4,4))-W(4,5);

Ya(3,4)=0;

Ya(4,3)=0;

disp(Ya)

运行结果