ypositive.m

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

function Yp=Ypositive(nb,ng,nl,nt)

Yp=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(:,3);
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)));
    Yp(x(i),x(i))=Yp(x(i),x(i))+tap_effect11;
    Yp(y(i),y(i))=Yp(y(i),y(i))+tap_effect22;
    Yp(x(i),y(i))=Yp(x(i),y(i))+tap_effect12;
    Yp(y(i),x(i))=Yp(y(i),x(i))+tap_effect21;
end

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

x=bb(:,1);
y=bb(:,2);
for i=1:length(x)
    Yp(x(i),x(i))=Yp(x(i),x(i))+(1/y(i));
end

%incorporation of transformer elements
fid=fopen('transformerdata.txt','r');
bb=textread('transformerdata.txt');
fclose(fid);
x=bb(:,1);
y=bb(:,3);
z=bb(:,5);
for i=1:length(x)
    Yp(x(i),x(i))=Yp(x(i),x(i))+1/z(i);
    Yp(y(i),y(i))=Yp(y(i),y(i))+1/z(i);
    Yp(x(i),y(i))=Yp(x(i),y(i))+(-1)*1/z(i);
    Yp(y(i),x(i))=Yp(y(i),x(i))+(-1)*1/z(i);
end