1、ft1+Tcos2(nu)t)dt = +Tsin2(nu)t)dt =扌 (1)tl+T tl+T21 cos(mcot) cos(noot) dt =sinfmcot) sin(nu)t) dt = O, mHn (2)J;t;+Tsin(m3t)cos(n3t)dt = 0 m、n 为任何整数 (3)对于任何一个周期为T的周期信号f (t),都可以求出它在上述各函数中的分量,从 而将此函数在区间(tl, tl+T)内表示为f(t)=乎 + d cos(a)t) + ax cos(2at) + - + an cos(noot) + bx sin(oot) + b2 sin(2cot) +.
2、+ bn sin(nu) + =+ Z=i3n cos(n(x)t) + bnsin (nu)t) (4)这就是函数f(t)在上述区间内的三角傅里叶级数表达式。方波函数表达式:先把这个函数展开为三角级数,为此就要求出分量系数a和b。(6)ao =訂:F(t)dt =辛fjdt - g dt = 0Qn =訂0 f(t)cos (n3t)dt =jcos(na)t) dt J? cos(nu)t)j = 0(7)当n为奇数当n为偶数=| 历sin(n3t) dt g sin(nu)t) dtj = nTT(8)因此,该非周期性方波在区间(0, T)内可以表示为(9)4 11f(t) = sin(
3、a)t) + -sin(3a)t) + - sin (Scot) + 3 S3建立模型描述4模块功能分析或源程序代码figure ts=0.0001;t=0:ts:4*pi;%t的取值范圉为0到4切,间隔0.0001取一个点f=l/ts;N=lengtli(t);yl=square(0.32*pi*t);% 产生一个方波plot(t,yl);title/产生一个方波ipausehold onfigure subplot 321y2=4/pi*sin(t);%频率为 1 (f“/2*pi)的正弦基波plot(t,y2);title/正弦基波)subplot 322y3=4/p 计(sin(3*t
4、)/3);% 三次谐波plot(t,y3);title(三次谐波Jsubplot 323y4=4/pi*(sin(5*t)/5);% 五次谐波plot(t,y4);titleC五次谐波Jsubplot 324y5=4/pi*(sin(7*t)/7);% 七次谐波plot(t,y5);title/七次谐波Jsubplot 325y6=4/pi * (sin(9 * t)/9); % 九次谐波plot(t,y6);title(九次谐波1) ;pausefigure(3)subplot 221y7=4/pi*(sin(t)+sin(3*t)/3);% 将前两次谐波叠加plot(t,y7);title
5、C前两次谐波叠加】subplot 222y8=4/pi*(sin(t)4-sin(3*t)/3+sin(5*t)/5);% 将前三次谐波叠加plot(t,y8);titled前三次谐波叠加Jsubplot 223y9=4/p 计(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:15x=x+4/pi*(sin(k*t)/k);y
6、(k+l)/2,:)=x;end%将各波形叠合绘出pause,figure(4),plot(t,y(l :9,:), titleC各波形叠合5调试过程及结论5.1调试过程在调试过程中,编写方波时由于当时没有注意方波扫描的精确度,将ts的数值选取的 过于大了而导致方波的波形出现了大的失真,经过反复儿次的调试后终于得出了一个比较 理想的方波图形。由于在设计之初是要将方波与后来要形成的多次正弦波贰加后的波形放 到一起,所以对丁方波周期的选取以及方波振幅的确定也尝试了好多次。最终才将这两个 数值确定下来。本事课程设计的日的是要证明方波能够分解为多次正眩波之和,所以程序 的设计目的在于体现分解的过程以及
7、最后n次谐波叠加后所形成的波形是否为方波。根据 上述公式9得出的方波表达式,令3“后分别画出y=4/pi*sin(t)的正弦基波、 y=4/pi*(siii(3*t)/3)的三次谐波、y=4/pi*(sin(5*t)/5)的五次谐波、y=4/pi*(siii(7*t)/7)的七次谐 波以及y=4/pi*(sin(9*t)/9)W九次谐波,为了能够充分的显示风波分解为多次正弦波之和的 过程,又分別作出了前两、三、五次谐波叠加的波形。最后为了能够完美的证明方波分解 为多次正弦波之和,利用了一个for循环语句画出了经过前n次谐波廉加后的波形(其中 令n二15),并且与之前程序最开始画出的方波波形相互
8、叠加。由于经过三角级数傅里叶分解得到的方波表达式是一个非周期性函数,故只有当无限 次谐波蒂加后才能形成正真的方波,所以本次的实验是存在一定的误差的,本想在程序中 作出其误差的图形表示以及计算出不同次数的方波卷加形成的方波所造成的误差,但由于知识以及能力有限,经过多次尝试无果后放弃了计算误差这部分项目。造成了一定的遗憾。5.2结论经过此次的实验,充分的证明了方波能够分解为多次正弦波z和,叠加次数越多误差越小越接近方波真实的波形。所得波形如下所示:产生一个方波图5-1方波波形图5-3叠加后的波形图5-4前n次廉加后的波形6心得体会通过此次的课程设计可以知道利用Matlab所得出的结果与利用理论知识
9、通过计算得 出的结果是一样的,方波,正弦波,以及正弦波廉加后的波形都能够利用程序编写出來的, 并且通过利用Matlab所得出的波形能够很清晰的表达出各个波形的物理意义。Mat lab是 一款非常强大非常实用的软件,利用Mat lab我们不仅能够证明许多数学上的问题,而且 许多我们通过理论知识所无法计算的问题都能够解决,比如说我们这次的课程设计,通过 人工计算我们能够得出方波的公式,但是这个公式却是一个非周期性的,我们无法准确的 算出方波的波形。但是利用Mat lab却是能够做到。通过利用Matlab进行仿真,仅仅用一 个循环语句就能够让我们的波形无限的接近方波的真实波形,而且我们还能利用Mat
10、 lab 计算出不同次正弦波卷加后所产生的误差并且通过图形模式表示出来以及我们还可以通 过建立一个三维坐标模型将方波分解为多次正弦波之和的趋势表示出來,使我们更加山观 的了解方波是如何分解为多次正弦波之和的。这次的课程设计让我真的很难忘,经过两个多星期的努力,终丁顺利完成了课程设计。 开始做课程设计不知道从何入手,困难很多,经过査阅资料,和同学讨论,终于了解了许 多。在做课程设计的过程中,我学会了很多,最主要的是我对Mat lab的运用更加熟练了, 同时对信号这门课的知识乂弄懂了不少,增加了对信号课程的学习兴趣。课程设计是每个 大学生必须面临的一项综合素质的考验,如果说在我们的学习阶段是一个知
11、识的积累过程, 那么现在的课程设计就是对过去所学的知识的综合应用,是对理论进行深化和重新认识的 实践活动。在这期间,我们有艰辛的付出,当然也有丰收的喜悦。首先,学习能力和解决 问题的信心都得到了提高。通过这次课程设计,我不仅对理论有了更深一步的认识,还培 养了自学能力和解决问题的能力,更重要的是,培养了克服困难的勇气和信心。7参考文献1吴大正.信号与线性系统分析(第四版).北京:高等教育出版社.2010刘卫国.MATLAB程序设计教程(第二版).北京:中国水利水电出版社.20103梁虹,梁洁,陈跃斌.信号与线性系统分析及MATLAB实现北京:电子工业出版社.20024韩利你王华.MATLAB电子仿真与应用M.北京:国防工业出版社,2001.5陈怀,高西全.MATLAB及在电子信息课程中的应用M.北京:电子工业出版社,2003.
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1