电网电能质量控制作业.doc

电网电能质量控制作业.doc

《电网电能质量控制作业.doc》由会员分享，可在线阅读，更多相关《电网电能质量控制作业.doc（4页珍藏版）》请在冰豆网上搜索。

电网电能质量控制作业

姓名：

李明善学号：

20071901088

u（t）=sin（wt）+0.3sin（3wt）+0.1sin（5wt）+0.08sin（7wt）对该电压方程用离散化傅里叶级数编程，求各次谐波含量

N=512;%采样点数N

T=0.02;%采样时间T

t=linspace（0,T,N）;%给出N个采样时间ti（i=1:

N）

w=100*pi;%w赋值

U=sin（w*t）+0.3*sin（3*w*t）+0.1*sin（5*w*t）+0.08*sin（7*w*t）;%求各采样点样本值x

%plot（t,U）%绘制电压波形图像

%grid

%xlabel（'t'）;%x轴变量“t”

%ylabel（'U（t）'）;%y轴变量“U（t）”

%title（'电压波形'）;%图像1名称“电压波形”

%holdon;

%plot（t,zeros（1,length（t））,'k'）%绘y=0图线

%holdoff;

fork=1:

N

forn=1:

N

x（n）=U（n）*exp（-j*2*pi*k/N*n）;

end

X（k）=sum（x）/N*2;

end

%bar（abs（X（1:

10）））%绘制谐波含量图

%grid

%xlabel（'Frequency'）;%x轴变量“Frequency”

%ylabel（'|F（K）|'）;%y轴变量“F（K）|”

%title（'振幅—频率图'）;%图像2的名称“振幅-频率图”

fork=1:

N

forn=1:

N

a（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）^2+b1（k,n）^2）;

end

end

%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'）;%x轴变量“n”

%ylabel（'Cn'）;%y轴变量“cn”

%title（'各次谐波'）%图像3名称“各次谐波”

forn=2:

20

disp（n）;disp（'次谐波为'）;disp（C（n,N））;%输出谐波含量

end

thd=0;

forn=2:

20

thd=thd+C（n,N）*C（n,N）;

end

disp（'总畸变率'）;disp（sqrt（thd））;%求畸变率

2次谐波为：

9.3130e-004

3次谐波为：

0.2987

4次谐波为：

0.0019

5次谐波为：

0.0984

6次谐波为：

0.0015

7次谐波为：

0.0781

8次谐波为：

0.0026

9次谐波为：

0.0018

10次谐波为：

0.0015

11次谐波为：

0.0013

12次谐波为：

0.0011

13次谐波为：

0.0010

14次谐波为：

9.1501e-004

15次谐波为：

8.4160e-004

16次谐波为：

7.8004e-004次谐波为

17次谐波为：

7.2752e-00418

18次谐波为：

6.8207e-004

19次谐波为：

6.4228e-004

20次谐波为：

6.0713e-004

总畸变率：

0.3241