程序2_电能质量_电能分析_电能_电力系统功率_powerquality_

close all;clear all;clc A=xlsread('C:\Users\lenovo\Desktop\电能质量\数据\数据\德昌数据.xlsx'); %-------------------------------scott谐波电流-----------------------------% for i=1:28800 Uamo(i)=A(i,1);Uaag(i)=A(i,4); Ua(i)=Uamo(i)*(cos(Uaag(i))+1i*sin(Uaag(i))); Ubmo(i)=A(i,13);Ubag(i)=A(i,16); Ub(i)=Ubmo(i)*(cos(Ubag(i))+1i*sin(Ubag(i))); Iamo(i)=A(i,7);Ia3mo(i)=A(i,8);Ia5mo(i)=A(i,9); Iaag(i)=A(i,10);Ia3ag(i)=A(i,11);Ia5ag(i)=A(i,12); Iaag1(i)=Iaag(i)-Uaag(i)-pi/2;Ia3ag1(i)=Ia3ag(i)-Uaag(i)-3*pi/2;Ia5ag1(i)=Ia5ag(i)-Uaag(i)-5*pi/2; Ia(i)=Iamo(i)*exp(Iaag1(i)*1j);Ia3(i)=Ia3mo(i)*exp(Ia3ag1(i)*1j);Ia5(i)=Ia5mo(i)*exp(Ia5ag1(i)*1j); Ibmo(i)=A(i,19);Ib3mo(i)=A(i,20);Ib5mo(i)=A(i,21); Ibag(i)=A(i,22);Ib3ag(i)=A(i,23);Ib5ag(i)=A(i,24); Ibag1(i)=Ibag(i)-Ubag(i);Ib3ag1(i)=Ib3ag(i)-Ubag(i);Ib5ag1(i)=Ib5ag(i)-Ubag(i); Ib(i)=Ibmo(i)*exp(Ibag1(i)*1j);Ib3(i)=Ib3mo(i)*exp(Ib3ag1(i)*1j);Ib5(i)=Ib5mo(i)*exp(Ib5ag1(i)*1j); IA1(i)=-0.25/sqrt(3).*Ib(i)+0.25.*Ia(i); IA13(i)=-0.25/sqrt(3).*Ib3(i)+0.25.*Ia3(i); IA15(i)=-0.25/sqrt(3).*Ib5(i)+0.25.*Ia5(i); IB1(i)=-0.25/sqrt(3).*Ib(i)-0.25.*Ia(i); IB13(i)=-0.25/sqrt(3).*Ib3(i)-0.25.*Ia3(i); IB15(i)=-0.25/sqrt(3).*Ib5(i)-0.25.*Ia5(i); IC1(i)=0.5/sqrt(3)*Ib(i); IC13(i)=0.5/sqrt(3).*Ib3(i); IC15(i)=0.5/sqrt(3).*Ib5(i); IA(i)=abs(IA1(i));IA3(i)=abs(IA13(i));IA5(i)=abs(IA15(i)); IB(i)=abs(IB1(i));IB3(i)=abs(IB13(i));IB5(i)=abs(IB15(i)); IC(i)=abs(IC1(i));IC3(i)=abs(IC13(i));IC5(i)=abs(IC15(i)); T(i)=3*i; UN=110; Sd1=300;Sd2=500;Sd3=750; Z1=UN^2/Sd1;Z2=UN^2/Sd2;Z3=UN^2/Sd3; Uamo(i)=A(i,1);Uaag(i)=A(i,4); Ua(i)=Uamo(i)*(cos(Uaag(i))+1i*sin(Uaag(i))); Ubmo(i)=A(i,13);Ubag(i)=A(i,16); Ub(i)=Ubmo(i)*(cos(Ubag(i))+1i*sin(Ubag(i))); UA1(i)=-2./sqrt(3).*Ub(i)+2.*Ua(i);UA(i)=abs(UA1(i)); UB1(i)=-2./sqrt(3).*Ub(i)-2.*Ua(i);UB(i)=abs(UB1(i)); UC1(i)=4./sqrt(3).*Ub(i);UC(i)=abs(UC1(i)); %-----------------------------------------------------------------------------% UA31(i)=3*Z1*IA3(i)/1000;UB31(i)=3*Z1*IB3(i)/1000;UC31(i)=3*Z1*IC3(i)/1000; UA32(i)=3*Z2*IA3(i)/1000;UB32(i)=3*Z2*IB3(i)/1000;UC32(i)=3*Z2*IC3(i)/1000; UA33(i)=3*Z3*IA3(i)/1000;UB33(i)=3*Z3*IB3(i)/1000;UC33(i)=3*Z3*IC3(i)/1000; %-----------------------------------------------------------------------------% UA51(i)=5*Z1*IA5(i)/1000;UB51(i)=5*Z1*IB5(i)/1000;UC51(i)=5*Z1*IC5(i)/1000; UA52(i)=5*Z2*IA5(i)/1000;UB52(i)=5*Z2*IB5(i)/1000;UC52(i)=5*Z2*IC5(i)/1000; UA53(i)=5*Z3*IA5(i)/1000;UB53(i)=5*Z3*IB5(i)/1000;UC53(i)=5*Z3*IC5(i)/1000; %-----------------------------------Sd1=300-----------------------------------% THDA1(i)=sqrt((UA31(i)./UA(i)).^2+(UA51(i)./UA(i)).^2); THDB1(i)=sqrt((UB31(i)./UB(i)).^2+(UB51(i)./UB(i)).^2); THDC1(i)=sqrt((UC31(i)./UC(i)).^2+(UC51(i)./UC(i)).^2); HRUA31(i)=UA31(i)./UA(i);HRUA51(i)=UA51(i)./UA(i); HRUB31(i)=UB31(i)./UB(i);HRUB51(i)=UB51(i)./UB(i); HRUC31(i)=UC31(i)./UC(i);HRUC51(i)=UC51(i)./UC(i); %-----------------------------------Sd1=500-----------------------------------% THDA2(i)=sqrt((UA32(i)./UA(i)).^2+(UA52(i)./UA(i)).^2); THDB2(i)=sqrt((UB32(i)./UB(i)).^2+(UB52(i)./UB(i)).^2); THDC2(i)=sqrt((UC32(i)./UC(i)).^2+(UC52(i)./UC(i)).^2); HRUA32(i)=UA32(i)./UA(i);HRUA52(i)=UA52(i)./UA(i); HRUB32(i)=UB32(i)./UB(i);HRUB52(i)=UB52(i)./UB(i); HRUC32(i)=UC32(i)./UC(i);HRUC52(i)=UC52(i)./UC(i); %-----------------------------------Sd1=750-----------------------------------% THDA3(i)=sqrt((UA33(i)./UA(i)).^2+(UA53(i)./UA(i)).^2); THDB3(i)=sqrt((UB33(i)./UB(i)).^2+(UB53(i)./UB(i)).^2); THDC3(i)=sqrt((UC33(i)./UC(i)).^2+(UC53(i)./UC(i)).^2); HRUA33(i)=UA33(i)./UA(i);HRUA53(i)=UA53(i)./UA(i); HRUB33(i)=UB33(i)./UB(i);HRUB53(i)=UB53(i)./UB(i); HRUC33(i)=UC33(i)./UC(i);HRUC53(i)=UC53(i)./UC(i); end figure plot(T,THDA1);grid on; title('A相电压畸变率(Sd=300MVA)'); figure plot(T,THDA2);grid on; title('A相电压畸变率(Sd=500MVA)'); figure plot(T,THDA3);grid on; title('A相电压畸变率(Sd=750MVA)'); figure plot(T,THDB1);grid on; title('B相电压畸变率(Sd=300MVA)'); figure plot(T,THDB2);grid on; title('B相电压畸变率(Sd=500MVA)'); figure plot(T,THDB3);grid on; title('B相电压畸变率(Sd=750MVA)'); figure plot(T,THDC1);grid on; title('C相电压畸变率(Sd=300MVA)'); figure plot(T,THDC2);grid on; title('C相电压畸变率(Sd=500MVA)'); figure plot(T,THDC3);grid on; title('C相电压畸变率(Sd=750MVA)');