電網(wǎng)電能質(zhì)量控制作業(yè)

1、電網(wǎng)電能質(zhì)量控制作業(yè)對(duì)該電壓方程用姓名：李明善 學(xué)號(hào)： 20071901088u(t)=sin(wt)+0.3sin(3wt)+0.1sin(5wt)+0.08sin(7wt)離散化傅里葉級(jí)數(shù)編程，求各次諧波含量% 采樣點(diǎn)數(shù) N% 采樣時(shí)間 T% 給出 N 個(gè)采樣時(shí)間 ti( i=1:N)% w 賦值% 求各采樣點(diǎn)樣本值 x% 繪制電壓波形圖像N=512;T=0.02;t=linspace(0,T,N);w=100*pi;U=sin(w*t)+0.3*sin(3*w*t)+0.1*sin(5*w*t)+0.08*sin(7*w*t);%plot(t,U)%grid% x 軸變量 “ t”% y

2、 軸變量 “U(t) ”% 圖像 1 名稱(chēng) “電壓波形”% 繪 y=0 圖線%xlabel(t);%ylabel(U(t);%title( 電壓波形 );%hold on;%plot(t,zeros(1,length(t),k)%hold off;for k=1:Nfor n = 1:Nx(n) = U(n)*exp(-j*2*pi*k/N*n);endX(k) = sum(x)/N*2;% 繪制諧波含量圖% x 軸變量 “ Frequency ”% y 軸變量“ F(K)| ”% 圖像 2 的名稱(chēng)“振幅 - 頻率圖”end %bar(abs(X(1:10)%grid %xlabel(Freq

3、uency); %ylabel(|F(K)|);%title( 振幅頻率圖 );for k=1:Nfor n = 1:Na(n) = U(n)*cos(2*pi*k/N*n);b(n) = U(n)*sin(2*pi*k/N*n);A(k,n)=sum(a)/N*2;B(k,n)=sum(b)/N*2;C(k, n)=sqrt(a1(k, n)9+b1(k, n) A2);endend%plot(1:N,C(1,1:N),1:N,C(3,1:N),1:N,C(5,1:N),1:N,C(7,1:N)%grid%xlabel(n);%ylabel(Cn);%title( 各次諧波 )x 軸變量“n

4、 ”y 軸變量“ cn ” 圖像 3 名稱(chēng)“各次諧波輸出諧波含量for n =2:20disp(n);disp( 次諧波為 );disp(C(n,N);endthd=0;for n =2:20thd=thd+C(n,N)*C(n,N);% 求畸變率enddisp( 總畸變率 );disp(sqrt(thd);2 次諧波為： 9.3130e-0043 次諧波為： 0.29874 次諧波為： 0.00195 次諧波為： 0.09846 次諧波為： 0.00157 次諧波為： 0.07818 次諧波為： 0.00269 次諧波為： 0.001810 次諧波為： 0.001511 次諧波為： 0.001312 次諧波為： 0.001113 次諧波為： 0.001014 次諧波為： 9.1501e-00415 次諧波為： 8.4160e-00416 次諧波為： 7.8004e-004 次諧波為17 次諧波為：7.2752e-0041818 次諧波為：6.8207e-00419 次諧波為：6.4228e-00420 次諧波為：6.0713e-004總畸變率：0.3241抿幅亠頻率團(tuán)09 -0.8 -070B5 0