yzero.m

一个用matlab编写的用于故障诊断的源程序

function Yz=Yzero(nb,ng,nl,nt)

Yz=zeros(nb);
%incorporation of elements corresponding to line impedances
fid=fopen('linedata.txt','r');
bb=textread('linedata.txt');
fclose(fid);

realtap=bb(:,5);
imagtap=bb(:,6);
x=bb(:,1);
y=bb(:,2);
zp=bb(:,4);
for i=1:length(x)
    tap_effect11=(1/(zp(i))/((abs(realtap(i)+imagtap(i)*j))^2));
    tap_effect12=(-1/(zp(i))/((realtap(i)+imagtap(i)*j)'));
    tap_effect21=(-1/(zp(i))/((realtap(i)+imagtap(i)*j)));
    tap_effect22=(1/(zp(i)));
    Yz(x(i),x(i))=Yz(x(i),x(i))+tap_effect11;
    Yz(y(i),y(i))=Yz(y(i),y(i))+tap_effect22;
    Yz(x(i),y(i))=Yz(x(i),y(i))+tap_effect12;
    Yz(y(i),x(i))=Yz(y(i),x(i))+tap_effect21;
end

%incorporation of elements corresponding to transformers
fid=fopen('transformerdata.txt','r');
bb=textread('transformerdata.txt');
fclose(fid);
x=bb(:,1);
y=bb(:,3);
hvx=bb(:,2);
lvx=bb(:,4);
z=bb(:,6);
for i=1:length(x)
    if hvx(i)==0 & lvx(i)==0 %represents the case of grounded star at both ends
        Yz(x(i),x(i))=Yz(x(i),x(i))+1/z(i);
        Yz(y(i),y(i))=Yz(y(i),y(i))+1/z(i);
        Yz(x(i),y(i))=Yz(x(i),y(i))+(-1)*1/z(i);
        Yz(y(i),x(i))=Yz(y(i),x(i))+(-1)*1/z(i);
    end
    
    if hvx(i)==0 & lvx(i)==1 %represents the case of grounded star and ungrounded star
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==0 & lvx(i)==2 %represents the case of grounded star and delta
        Yz(x(i),x(i))=Yz(x(i),x(i))+1/z(i);
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==1 & lvx(i)==0 %represents the case of ungrounded star and grounded star
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==1 & lvx(i)==1 %represents the case of ungrounded star at both ends
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==1 & lvx(i)==2 %represents the case of ungrounded star and delta
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==2 & lvx(i)==0 %represents the case of delta and grounded star
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+1/z(i);
    end
    
    if hvx(i)==2 & lvx(i)==1 %represents the case of delta and ungrounded star
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
    
    if hvx(i)==2 & lvx(i)==2 %represents the case of delta at both ends
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
        Yz(y(i),y(i))=Yz(y(i),y(i))+0;
        Yz(x(i),y(i))=Yz(x(i),y(i))+0;
        Yz(y(i),x(i))=Yz(y(i),x(i))+0;
    end
end

%incorporation of generator elements
fid=fopen('gendata.txt','r');
bb=textread('gendata.txt');
fclose(fid);

x=bb(:,1);
z=bb(:,3);
gnd=bb(:,4);
for i=1:length(x)
    if gnd(i)==1000
        Yz(x(i),x(i))=Yz(x(i),x(i))+0;
    end
    if gnd(i)~=1000
        Yz(x(i),x(i))=Yz(x(i),x(i))+1/(z(i)+3*gnd(i));
    end
end