《周期信号的合成与分解以及在matlab上的实现.docx》由会员分享，可在线阅读，更多相关《周期信号的合成与分解以及在matlab上的实现.docx（5页珍藏版）》请在三一办公上搜索。
1、周期信号的合成与分解以及在matlab上的实现广东工业大学实验报告 信息工程学院 成绩评定_ 学号： 姓名： 教师签名_ 预习情况 操作情况 考勤情况 数据处理情况 实验4题目：周期信号的合成与分解 第7周 星期3第 8、9 节 一、实验目的 1.熟悉信号的合成、分解原理，加深对傅立叶级数的理解； 2.了解和认识吉布斯现象。 一、 实验原理 信号可以分解为一个直流分量和许多不同频率的正弦分量之和。主要表现各频率的正弦分量在信号所占比重大小的不同。 根据周期信号的傅立叶级数展开式可知，任何非正弦周期信号，只要满足狄里赫利条件都可以分解为一直流分量和由基波及各次谐波分量的叠加。 同样，由基波及各次
2、谐波分量也可以叠加出一个周期方波信号。至于叠加出来的信号与原信号的误差，则取决于傅立叶级数的项数。根据傅立叶的原理，任何周期信号都可以分解为用一组三角函数sin，cos的组合表示。 合成波形所包含的谐波分量越多，除间断点附近外，它越接近于原方波信号，在间短点附近，随着所含谐波次数增高，合成波形的尖峰越靠近间断点，但尖峰幅度并未明显减小，即当合成波形包含的谐波次数n，在间断点附近仍有9%的偏差，这种现象称为吉布斯现象。 二、 验证性实验 1、信号的合成、分解方法一： clear all clc x0=-pi:0.01:pi; sum=0.0; for(n0=1:1:5) p0=4*sin(n0*
3、pi/2)*cos(n0*x0)/(n0*pi); sum=sum+p0; end plot(x0,sum,b) text(2.1,-1.2,n=5) hold on; x1=-pi:0.01:pi; sum=0.0; for(n1=1:1:20) p1=4*sin(n1*pi/2)*cos(n1*x1)/(n1*pi); sum=sum+p1; end plot(x1,sum,r) text(3.1,-1.2,n=20) hold on; x2=-pi:0.01:pi; sum=0.0; for(n2=1:1:150) p2=4*sin(n2*pi/2)*cos(n2*x2)/(n2*pi)
4、; sum=sum+p2; end plot(x2,sum,y) text(1.6,1.3,n=150) hold on; x3=-pi:0.01:pi; sum=0.0; for(n3=1) p3=4*sin(n3*pi/2)*cos(n3*x3)/(n3*pi); sum=sum+p3; end plot(x3,sum,g) text(2,-0.5,n=1) hold on; x4=-pi:0.01:pi/2; y1=-1; plot(x4,y1,g) hold on x5=-pi/2:0.01:pi/2; y2=2; plot(x5,y2,g) hold on x6=pi/2:0.01:
5、pi; y3=-1; plot(x6,y3,g) hold on y4=-1:0.01:1; x7=-pi/2; plot(x7,y4,g) hold on x8=pi/2; plot(x8,y4,g); xlabel(x); ylabel(y); 2、信号的合成、分解方法二： clear all clc T=2; w0=2*pi/T; fw=3; t=-fw:0.001:fw; N=input(Number Of Harmonics:); a0=0; An=a0*ones(1,length(t); for n=1:2:N An=An+4/n/pi*sin(n*w0*t); end plot(
6、t,An); ylim(-1.5 1.5); hold on plot(t,zeros(1,length(t); hold off clear all clc T=2; w0=2*pi/T; fw=3; t=-fw:0.001:fw; N=input(Number Of Harmonics:); F0=0; Fn=F0*ones(1,length(t); for n=1:2:N Fn=Fn+2/pi/j*(exp(j*n*w0*t)-exp(-j*n*w0*t)/n); end plot(t,Fn); ylim(-1.5 1.5); hold on plot(t,zeros(1,length(
7、t); hold off 三、 程序设计实验： clear all clc T=0.02; w0=2*pi/T; fw=0.03; t=-fw:0.001:fw; N=9; a0=0; An=a0*ones(1,length(t); for n=1:2:N; An=An+3*(4/n/pi*sin(n*w0*t); end plot(t,An); ylim(-5 5); hold on plot(t,zeros(1,length(t); hold off 五、思考与讨论 设计一个三角波合成实验，写出实验步骤和程序。 clear all clc T=2; w0=2*pi/T; fw=3; t=-fw:0.001:fw; N=100; a0=0; An=a0*ones(1,length(t); for n=1:2:N; An=An+3*(4/n2/pi2*cos(n*w0*t); end plot(t,An); ylim(-5 5); hold on plot(t,zeros(1,length(t); hold off 543210-1-2-3-4-5-3-2-10123
