程序二：dq變換法程序

fid=fopen('dxjd.txt','rt');
u=fscanf(fid,'%f',[10,2000]);
Va=u(8,:)                              %a相電壓
t=0:0.0001:0.1999;
%figure(1);
%plot(t,Va);
%title('110kV輸電線路A相短路故障負(fù)荷側(cè)A相電壓');xlabel('t(秒)'),ylabel('伏');
status=fclose(fid); %k(i)=i;                                                                                                             
wt=100*pi*t
％@@@@@@@@@@@@構(gòu)造b相和c相@@@@@@@@@@@@@@@@@
for i=1:2000
    if (i>33)                                 
        j=i-33;                               %33點相當(dāng)于60度
    else
        j=i;
    end
    Vc(i)=-Va(j);                    
    Vb(i)=-Vc(i)-Va(i);
end
    a=2*pi/3.0;
for  i=1:2000
      t=0:0.0001:0.1999;
      wt=100*pi*t;
        dsys(i)=sin(wt(i))*Va(i)+sin(wt(i)-a)*Vb(i)+sin(wt(i)+a)*Vc(i);      %d軸
        qsys(i)=-cos(wt(i))*Va(i)-cos(wt(i)-a)*Vb(i)-cos(wt(i)+a)*Vc(i);     %q軸
        Vd(i)=0.81649658*dsys(i);                                   %d軸有效值   
        Vq(i)=0.81649658*qsys(i);                                   %q軸有效值
        Vsag(i)=sqrt(Vd(i)*Vd(i)+Vq(i)*Vq(i))/1.7320508;                %瞬時有效值
        da(i)=asin(-Vq(i)/Vsag(i)/1.7320508)*180/pi-90;                  %相位跳變
        dwt(i)=asin(-Vq(i)/Vsag(i)/1.7320508);  
        %df(i)=dwt(i)/(2*pi*t(i));
end
       figure(1);
       subplot(3,1,1); plot(t,Va);%title('a相電壓');
       ylabel('伏'),%xlabel('t（秒）');
        %subplot(2,2,1);plot(t,Vd);title('d軸有效值');
        %ylabel('伏'),xlabel('t(秒)');
        %subplot(2,2,2);plot(t,Vq);title('q軸有效值');
        %ylabel('伏'),xlabel('t(秒');
       subplot(3,1,2);plot(t,Vsag);%title('dq變換所得的a相電壓有效值');
       ylabel('伏'),%xlabel('t（秒）');
        subplot(3,1,3);plot(t,da);%title('相位跳變'); 
        ylabel('da'),xlabel('t（秒）');