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

精品文档 精品文档 0 0欢迎。下载 目录 TOC \o 1-5 \h \z \o Current Document 技术要求 1 \o Current Document 基本原理 1 \o Current Document 建立模型描述 2 \o Current Document 模块功能分析或源程序代码 2 \o Current Document 调试过程及结论 4 5.1 调试过程 . 4 \o Current Document 5.2 结论 5 \o Current Document 心得体会 7 \o Current Document 参考文献 . 8 精品文档 精品文档 PAGE PAGE #欢迎下载 方波分解为多次正弦波之和的设计 1技术要求 已知某一周期性方波，用 matlab仿真软件演示谐波合成情况，讨论参数对分解和合 成波形的影响。 2基本原理 根据三角傅里叶级数: cos2 (n 3 )dt =1+Tsin2(n 3 )dt = cos2 (n 3 )dt = 1+T sin2(n 3 )dt = * t1+T 1 cos(mwt) cos(n 3 ) dt = f+T sin(m3t) sin(n 3 ) dt = 0,m 工 n (1) (2) t1+T sin( m3t) cos(n 3 ) dt = 0 m、 n为任何整数 (3) 对于任何一个周期为T的周期信号f (t),都可以求出它在上述各函数中的分量，从 而将此函数在区间(t1 , t1+T)内表示为 f(t) = 0° + a1 cos( 3 ) + a1 cos(2 3 ) + ? + an cos(n 3 t + ? + b1 sin( 3)+6 sin(2 3 ) + TOC \o 1-5 \h \z + bn sin(n 3 + ?=乎 + ?冰1 [an cos(n 3 ) + bnSin(n 3 t)] (4) 这就是函数f(t)在上述区间内的三角傅里叶级数表达式 方波函数表达式： { f(t) = 1 OvtvT/2 (5) { f(t) T/2^ t T () 先把这个函数展开为三角级数，为此就要求出分量系数 a和bo T 2 T 2 一 T 0 = a ￡ f(t)dt = -[ ￡ dt - f dt] = 0 (6) 1 1 2 Tn a= 2 彳 T n a= 2 彳f(t)cos (n 3 t)dt = 2[彳cos(n 3 ) dt - /cos(n 3 )] = 0 2 bn = 2 / f(t)sin (n 3 t)dt= 2[ #sin(n 3 ) dt- / sin(n 3 ) dt] = { 2 当n为奇数 当n为偶数 (7) (8) 因此，该非周期性方波在区间(0, T)内可以表示为 TOC \o 1-5 \h \z 4 1 1 (t ) f= [si n( 3) + 3si n(3 3 ) + 5sin(5 3 ) + ? ] (9) n 3 5 3建立模型描述 利用Matlab画出方波波形根据公式（9）令3=1 - 1分别作出各次谐波的波形 —与方波比较 U*将各波形图叠加作图 分别作出不同次谐波叠加后的波形 得出结论：方波可以分解为多次正弦波之和 4模块功能分析或源程序代码 figure⑴ ts=0.0001; t=0:ts:4*pi;%t 的取值范围为0到4*pi，间隔0.0001取一个点 f=1/ts; N=le ngth(t); y1=square(0.32*pi*t);% 产生一个方波 plot(t,y1); title(产生一个方波);pause hold on figure(2) subplot 321 y2=4/pi*sin(t);% 频率为 1 (f=1/2*pi )的正弦基波 % plot(t,y2); title( 正弦基波) 精品文档 精品文档 PAGE PAGE #欢迎。下载 subplot 322 y3=4/pi*(sin(3*t)/3);% 三次谐波 plot(t,y3); title( 三次谐波 ) subplot 323 y4=4/pi*(sin(5*t)/5);% 五次谐波 plot(t,y4); title( 五次谐波 ) subplot 324 y5=4/pi*(sin(7*t)/7);% 七次谐波 plot(t,y5); title( 七次谐波 ) subplot 325 y6=4/pi*(sin(9*t)/9);% 九次谐波 plot(t,y6); title( 九次谐波 );pause figure(3) subplot 221 y7=4/pi*(sin