matlab求节点导纳矩阵,Matlab个节点导纳矩阵形成的程序问题

下面这个是我改过的 应该更对点 不过还是没有什么显示

clear

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

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

B1=input('请输入支路信息：B1=');

%每行是一条支路

%第一列是支路的一个端点I

%第二列是支路的另一个端点J

%第三列是支路的电阻R

%第四列是支路的电抗X

%第五列是支路的对地电纳或变压器变比k

Y=zeros(N1);

G=zeros(N1);

B=zeros(N1);

for m1=1:N1

I=B1(m1,1);

J=B1(m1,2);

R=B1(m1,3);

X=B1(m1,4);

k=B1(m1,5);

if I*J==0

if I==0

G(J,J)=R;

B(J,J)=X;

else

G(I,I)=R;

B(I,I)=X;

end

else if I*J>0

B(I,I)=k;

B(J,J)=k;

k=1;

else

if I<0

t=I;

I=J;

J=t;

end

J=-J;

if k>0

B(I,I)=k;

B(J,J)=k;

k=1;

else

if k<0

k=-1/k;

end

G(I,J)=-(1/k)*R/(R^2+X^2);

B(I,J)=(1/k)*R/(R^2+X^2);

G(I,J)=G(J,