文档介绍：用程序实现节点导纳矩阵的形成:
主程序:
%节点导纳矩阵的形成
%网络的节点数为n
n=input('网络的节点数');
A=input('节点相关性矩阵');%1表示两节点相连,0表示两节点不连
Y=zeros(n,n);
for i=n:-1:1
for j=i-1:-1:1
if (A(i,j))
i
j
Y(i,j)=input('i,j支路之间的导纳值');
Y(j,i)=Y(i,j);
end
end
Y(i,i)=sum(Y(i,:));
end
Y
%对节点导纳矩阵的修改
%0表示不修改;
%1表示给原有网络增加节点,且节点处的阻抗为Zij;
%2表示在节点i,j之间增加一条阻抗为Zij支路;
%3表示在原有网络i,j之间切除一条阻抗为Zij支路;
%4表示原有网络i,j之间支路阻抗由Zij变成Z1ij
m=1;
for m=1:100
m=input('导纳矩阵修改的情况分类');
switch(m)
case(0)
break
case(1)
i=input('在第i点增加阻抗');
Zij=input('增加节点的阻抗');
Y(n+1,:)=zeros(1,n);
Y(:,n+1)=zeros(n+1,1);
Y(n+1,n+1)= Y(i,i);
Y(i,n+1)=-1/Zij;
Y(n+1,i)= Y(i,n+1);
Y(i,i)= Y(i,i)+ Y(i,i);
Y
case(2)
i=input('请输入节点i');
j=input('请输入节点j');
Zij=input('请输入在节点i,j间增加的阻抗');
Y(i,i)= Y(i,i)+1/Zij;
Y(j,j)=Y(j,j)+1/Zij;
Y(i,j)=Y(j,j)-1/Zij;
Y(j,i)=Y(i,j);
Y
case(3)
i=input('请输入节点i');
j=input('请输入节点j');
Zij=input('请输入在节点i,j间增加的阻抗');
Y(i,i)= Y(i,i)-1/Zij;
Y(j,j)=Y(j,j)-1/Zij;
Y(i,j)=Y(j,j)+1/Zij;
Y(j,i)=Y(i,j);
Y
case(4)
i=input('请输入节点i');
j=input('请输入节点j');
Z1ij=input('请输入在节点i,j间变换的阻抗');
Y(i,i)= Y(i,i)-1/Y(i,j)+1/Z1ij;
Y(j,j)=Y(j,j)-1/Y(i,j)+1/Z1ij;
Y(i,j)=Y(j,j)+1/Y(i,j)-1/Z1ij;
Y(j,i)=Y(i,j);
Y
otherwise
disp('输入数据有误,请重新输入数据');
end
end
运行结果:
网络的节点数3
节点相关性矩阵[1,1,1;1,1,0;1,0,1]
i =
3
j =
1
i,j支路之间的导纳值5
i =
2
j =
1
i,j支路之间的导纳值2
Y =
7 2