dq_transform.txt

dq变换法程序,用于电能质量的监测和有功电流或无功电流的计算等。

程序二：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相短路故障负荷侧A相电压');xlabel('t(秒)'),ylabel('伏');
status=fclose(fid); %k(i)=i;                                                                                                             
wt=100*pi*t
％@@@@@@@@@@@@构造b相和c相@@@@@@@@@@@@@@@@@
for i=1:2000
    if (i>33)                                 
        j=i-33;                               %33点相当于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（秒）');