test.m

经常用于电力系统中的同步发电机两相短路的分析。

clc;
clear;
A=load('d:\Program Files\MATLAB71\matlab\work\NewFolder1\SXDL.txt');
XD=A(1,1);XAD=A(1,2);RXF=A(1,3);
RXD=A(1,4);RF=A(1,5);RD=A(1,6);
RXFD=A(1,7);XQ=A(1,8);XAQ=A(1,9);
RXH=A(1,10);RXQ=A(1,11);RH=A(1,12);
RQ=A(1,13);RR=A(1,14);DLTTY=A(1,15);
TMAX=A(1,16);U0=A(1,17);SITA0=A(1,18);
OMG=1.0;T=0.0;TY=0.0;
DLTT=DLTTY*314.15927;
DLTT2=DLTT/2.0;
num=TMAX/DLTTY;   %计算数据长度
I=zeros(6,num+1); %创建电流数据(Id,Iq,If,ID,IH,IQ形式)
UD=0.0;
UQ=U0;
% UQ=0;
I(3,1)=U0/XAD;
UF=I(3,1)*RF;    %15~18 故障边界条件（决定故障类型）
C=CYZFC(XD,XAD,RXF,RXD,RF,RD,RXFD,XQ,XAQ,RXH,RXQ,RH,RQ,RR,DLTT2,OMG);  %形成C因子表
D=zeros(6,2);
D(1,1)=DLTT2;D(2,2)=DLTT2; %形成D因子表
AIA=zeros(1,num+1);
AIB=zeros(1,num+1);
AIC=zeros(1,num+1);  %创建电流数据(Ia,Ib,Ic形式)
AIA(1,1)=I(1,1)*cos(SITA0)-I(2,1)*sin(SITA0);
AIB(1,1)=I(1,1)*cos(SITA0-2.0*pi/3.0)-I(2,1)*sin(SITA0-2.0*pi/3.0);
AIC(1,1)=I(1,1)*cos(SITA0+2.0*pi/3.0)-I(2,1)*sin(SITA0+2.0*pi/3.0);
for i=2:num+1
    T=T+DLTT;
    E=EYZFC(XD,XAD,RXF,RXD,RF,RD,RXFD,XQ,XAQ,RXH,RXQ,RH,RQ,RR,DLTT2,OMG,I(:,i-1),UD,UQ,UF);
   
    AG=C*D;BG=C*E; 
    str=sprintf('E=%2d \t\n',BG);disp(str);
    AG11=AG(1,1);AG12=AG(1,2);AG21=AG(2,1);AG22=AG(2,2);
    BG1=BG(1);BG2=BG(2);
    G11=cos(SITA0+T-2.0*pi/3.0)-cos(SITA0+T+2.0*pi/3.0);
    G12=-sin(SITA0+T-2.0*pi/3.0)+sin(SITA0+T+2.0*pi/3.0);
    G21=AG11*cos(SITA0+T)-AG21*sin(SITA0+T);
    G22=AG12*cos(SITA0+T)-AG22*sin(SITA0+T);
    F=-BG1*cos(SITA0+T)+BG2*sin(SITA0+T);
    UD=-G12*F/(G11*G22-G12*G21);
    UQ=G11*F/(G11*G22-G12*G21);
    I(:,i)=C*D(:,1)*UD+C*D(:,2)*UQ+BG;
    str=sprintf('%2d \t',I(:,i));disp(str);
    AIA(1,i)=I(1,i)*cos(SITA0+T)-I(2,i)*sin(SITA0+T);
    AIB(1,i)=I(1,i)*cos(SITA0+T-2.0*pi/3.0)-I(2,i)*sin(SITA0+T-2.0*pi/3.0);
    AIC(1,i)=I(1,i)*cos(SITA0+T+2.0*pi/3.0)-I(2,i)*sin(SITA0+T+2.0*pi/3.0);
end
% figure(1)
% plot(I(3,:))
% figure(2)
% plot(AIB)
% figure(3)
% plot(AIA)