节点导纳矩阵程序

n=5;%节点数为5
L=5;%支路数为5
B1=[125.931.50;1531101.05;230.8230.954545;34651000;45651000];
%每行是一条支路
%第一列是支路的一个端点I
%第二列是支路的另一个端点J
%第三列是支路的电阻R
%第四列是支路的电抗X
%第五列是支路的对地电纳或变压器变比k
%规定k*始终在I边Y=zeros（5）;
G=zeros（5）;
B=zeros（5）;
for m1=1:L
    I=B1（m11）;
    J=B1（m12）;
    R=B1（m13）;
    X=B1（m14）;
    k=B1（m15）;
    if J==0            %变压器励磁支路
            G（II）= G（II）+R/（R^2+X^2）;
            B（II）=  B（II）-X/（R^2+X^2）;
        end
    if I>0            %线路支路
            G（IJ）=-R/（R^2+X^2）;
            B（IJ）=X/（R^2+X^2）;
            G（JI）=G（IJ）;
            B（