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

电网电能质量控制作业姓名：李明善 学号：20071901088 u(t)=sin(wt)+0.3sin(3wt)+0.1sin(5wt)+0.08sin(7wt)对该电压方程用离散化傅里叶级数编程，求各次谐波含量 N=512; % 采样点数NT=0.02; % 采样时间Tt=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名称 “电压波形”%hold on;%plot(t,zeros(1,length(t),k) % 绘y=0图线%hold off;for k=1:N for n = 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的名称“振幅-频率图”for k=1:N for n = 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); 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); % x轴变量 “n”%ylabel(Cn); % y轴变量“cn”%title(各次谐波) % 图像3名称“各次谐波”for n =2:20 disp(n);disp(次谐波为);disp(C(n,N); % 输出谐波含量endthd=0;for n =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.001810次谐波为： 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-0041818次谐波为：6.8207e-00419次谐波为：6.4228e-00420次谐波为：6.07