方波分解为多次正弦波之和.docxVIP

目录 TOC \o 1-5 \h \z \o Current Document 1技术要求 1 \o Current Document 2基本原理 1 \o Current Document 3建立模型描述 2 4模块功能分析或源程序代码 2 \o Current Document 5调试过程及结论 4 \o Current Document 5.1调试过程 4 \o Current Document 5.2结论 5 \o Current Document 6心得体会 7 \o Current Document 7参考文献 8 方波分解为多次正弦波之和的设计 1技术要求 已知某一周期性方波，用mat lab仿真软件演示谐波合成情况，讨论参数对分解和合 成波形的影响。 2基本原理 根据三角傅里叶级数： J：； + Fcos (na)t)dt = J；； + rsin2(na)t)dt = f (1) ?f；； +「cos (mo)t)cos (ncot)dt = J；； + sin (mcot)sin (ncot)dt = 0,ni 工 n (2) J*：； * rsin (mo)t)cos (n(ot)dt = 0 m、n为任何整数 (3) 对于任何一个周期为T的周期信号f (t),都可以求出它在上述各函数中的分量，从 而将此函数在区间(tl, tl+T)内表示为 f(t) ao =y + a〔cos (a)t) + atcos (2a)t) + ... + ancos (niot) + ... + bjSin (o)t) + b2sin (2a)t) + ... + bnsin (no)) + ... = j + 宵二 Jancos (nu)t) + bnsin (na)t)] (4) 这就是函数f(t)在上述区间内的三角傅里叶级数表达式。 方波函数表达式： (f(t) = 1 0 tT/2 (f(t) T/2tT 先把这个函数展开为三角级数，为此就要求出分量系数和b。 (7) (7) T an=続f(t)cos (ncot)dt = | 几cos (ncot)dt - Jtcos (no) an =辛J：f(t)sin (na)t)dt = | J^sin (na)t)dt - jTsin (na)t)dt = 4 —当n为奇数 0当n为偶数 (8) 因此，该非周期性方波在区间(0, T)内可以表示为 4 1 1 f(t) = —[sin(3t) + 3sin (3(ot) + §sin (Scot) + …] (9) 3建立模型描述 figure (1) ts二0. 0001; t=0:ts:4*pi;%t的取值范围为0到4*pi,间隔0. 0001取一个点 f二1/ts; N=length(t); yl=square (0. 32*pi*t) ;%产生一个方波 plot (t, yl); titlef产生一个方波);pause hold on figure(2) subplot 321 y2=4/pi*sin(t);%频率为 1 (f=l/2*pi)的正弦基波% plot(t, y2); titleC正弦基波) subplot 322 y3=4/pi*(sin(3*t)/3) ;%三次谐波 plot(t,y3); titleC三次谐波) subplot 323 y4=4/pi*(sin(5*t)/5)俺五次谐波 plot(t, y4); titleC五次谐波) subplot 324 y5=4/pi*(sin(7*t)/7) ;%七次谐波 plot(t,y5); titleC七次谐波) subplot 325 y6=4/pi*(sin(9*t)/9)型九次谐波 plot(t,y6); title C 九次谐波):pause figure(3) subplot 221 y7=4/pi*(sin(t) +sin(3*t)/3);%将前两次谐波叠加 plot(t, y7); titleC前两次谐波叠加) subplot 222 y8=4/pi* (sin(t)+sin(3*t)/3+sin(5*t)/5) ;%将前三次谐波叠加 plot(t, y8); titleC前三次谐波叠加) subplot 223 y9=4/pi*(sin (t) +sin (3*t) /3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/t) ;%将 前五次谐波叠 加 plot(t, y9); titleC前五次谐波叠加) 弘重新定义y,把各次波形数据存为一个三维数组 y二zeros(10, max(size(t)));x二zeros(size(t)); for k=l:2:15 x=x+4/pi*(