周期信号的分解与合成.doc

《周期信号的分解与合成.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验一周期信号的分解与合成一、实验目的1．用同时分析法观测50Hz非正弦周期信号的频谱。2．观测基波和其谐波的合成。二、实验原理1．一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示，其中与非正弦具有相同频率的成分称为基波或一次谐波，其它成分则根据其频率为基波频率的2、3、4、...、n等倍数分别称二次、三次、四次、...、n次谐波，其幅度将随谐波次数的增加而减小，直至无穷小。2．不同频率的谐波可以合成一个非正弦周期波，反过来，一个非正弦周期波也可以分解为无限个不同频率的谐波成分。3．一个非正弦周期函数可用傅里叶级

2、数来表示，级数各项系数之间的关系可用一各个频谱来表示，不同的非正弦周期函数具有不同的频谱图，各种不同波形及其傅氏级数表达式见表1-1表1-1各种不同波形的傅里叶级数表达式（下）1．方波2．三角波3．半波4．全波5．矩形波三、预习要求在做实验前必须认真复习教材中关于周期性信号傅利叶级数分解的有关内容。四、实验内容1.50HZ方波信号的频谱。2.周期矩形脉冲的频谱；脉冲宽度为1；周期为4；则基波角频率为0.5pi3.使用不同频率的谐波合成方波信号；注意观察随着谐波数的增加合成的波形发生的变化。4.使用不同频率的谐波合成矩形脉冲

3、信号；注意观察随着谐波数的增加合成的波形。五、思考题1．什么样的周期性函数没有直流分量和余弦项？附：1.50HZ方波信号的频谱。>>w1=;%基波角频率>>n=0:1:30;>>bn=;%三角级数中系数bn，参考书p122>>stem(n*w1,bn),gridon>>xlabel('omega(rad/s)'),ylabel('bn')>>title('方波信号频谱分析图')2.周期矩形脉冲的频谱；脉冲宽度为1；周期为4；则基波角频率为0.5pitao=;w1=;n=-15:1:15;fn=;%矩形脉冲级数系数fn，参

4、考书p130，用matlab自带函数sincstem(n,fn),gridonxlabel('n');ylabel('Fn');title('周期矩形脉冲的频谱图');3.%使用不同频率的谐波合成方波信号；注意观察随着谐波数的增加合成的波形%发生的变化。t=-1:0.001:1;omega=2*pi;y=square(2*pi*t,50);plot(t,y);gridonxlabel('t');ylabel('周期方波信号');axis([-11-1.51.5]);n_max=[1351147];N=length(n_ma

5、x);fork=1:Nn=1:2:n_max(k);b=4./(pi*n);x=b*sin(omega*n'*t);figure;plot(t,y)holdon;plot(t,x);holdoff;xlabel('t');ylabel('部分和的波形');axis([-11-1.51.5]),gridontitle(['最大谐波数=',num2str(n_max(k))])end4.%使用不同频率的谐波合成矩形脉冲信号；注意观察随着谐波数的增加合成的波形%发生的变化。t=-2:0.001:2;omega=pi;y=rect

6、puls(t)+rectpuls(t-2)+rectpuls(t+2);axis([-2.22.201.2]);plot(t,y);gridonxlabel('t');ylabel('周期矩形脉冲信号');axis([-2.22.2-1.51.5]);n_max=[1351147];N=length(n_max);fork=1:Nn=1:1:n_max(k);a=2*sin(n*omega/2)./(n*pi);x=a*cos(omega*n'*t);x=x+1/2;figure;plot(t,y)holdon;plot(

7、t,x);holdoff;xlabel('t');ylabel('部分和的波形');axis([02-1.51.5]),gridontitle(['最大谐波数=',num2str(n_max(k))])end