叶变换+频率响应+抽样定理仿真.docxVIP

第6章周期信号的傅里叶级数及频谱分析6.1 实验目的学会运用MATLAB分析傅里叶级数展开，深入理解傅里叶级数的物理含义学会运用MATLAB分析周期信号的频谱特性6.2 实验原理及实例分析6.2.1 周期信号的傅里叶级数设周期信号,其周期为T，角频率为,则该信号可展开围三角形式的傅里叶级数，即（6-1）其中，各正弦项与余弦项的系数、成为傅里叶系数，更具函数的正交性，得（6-2）其中，n=1,2……。积分区间通常取为（0，T）或（，），将上式同频率项合并，可以改写为（6-3）由此，可以得出傅里叶级数中各系数间的关系为：（6-4）从物理概念上来说，式6-3中的即是信号的直流分量；式中第二项成为信号的基波或基波分量，它的角频率与元周期信号相同；式中第三项成为信号的二次谐波，它的频率是基波频率的二倍，以此类推。一般而言，成为信号的n次谐波，n比较法的那些分量统称为高次谐波。我们还成用到负指数形式的傅里叶级数，设周期信号，其周期为T，角频率为，该信号负指数形式的傅里叶级数为其中， n=0,……成为复指数形式傅里叶级数系数。利用MATLAB可直观地观察和分析周期洗脑傅里叶级数及其收敛性。【实例6-1】周期方波信号如图6-1所示，试求出该信号的傅里叶基数，利用MATLAB编程实现各次谐波的叠加，并验证其收敛性。解：从理论分析可知，已知周期方波信号的傅里叶级数展开式为+)取A=1，T=1,可分别求出1、3、5、11、47项傅里叶级数求和的结果，其MATLAB源程序为 t=-1:0.001:1; omega=2*pi; y=square(2*pi*t,50); plot(t,y),grid on xlabel(t),ylabel(周期方波信号) axis([-1 1 -1.5 1.5]) n_max=[1 3 5 11 47]; N=length(n_max); for k=1:Nn=1:2:n_max(k);b=4./(pi*n);x=b*sin(omega*n*t);figure;plot(t,y);hold on;plot(t,x);hold off;xlabel(t),ylabel(部分和的波形)axis([-1 1 -1.5 1.5]),grid ontitle([最大谐波数=,num2str(n_max(k))])end程序运行后，画出各项部分和的波形如图6-2所示：图6-2 周期方波信号的有限项傅里叶级数逼近从图可以看出，随着傅里叶级数项数的增多，部分和与周期方波型号的误差越来与小。在N=47项的时候，部分和的波箱与周期方波信号的波形很接近，但在信号的跳变点附近，却总是存在一个过冲，这就是所谓的Gibbs现象。6.2.2周期信号的频谱分析周期信号通过傅里叶级数分解可展开成一系列互相正交的正弦信号或复指数信号分量的加权和。在三角形式傅里叶级数中，各分量的形式为;在指数形式的傅里叶级数中，各量的形式为。对实信号而言，与承兑出现。对不同的周期信号，它们各个分量的数目、角频率、幅度或相位不同。傅里叶系数的幅度或随角频率变化关系绘制成图形，成为信号的幅度频谱，简称幅度谱。相位随角频率变化关系绘制成图形，称为信号的相位频谱，简称相位谱。幅度谱和相位朋友统称为信号的频谱。信号的频谱是信号的另一种表示，它提供了从另一个角度来观察和分析信号的途径。利用MATLAB命令可对周期信号的频谱及其特点进行观察验证分析。6.3 编程练习2.试用MATLAB分析图6-5中周期三角信号的频谱。当周期三角信号的周期好三角信号的宽度变化时，试观察分析其频谱的变化。解：根据傅里叶级数理论可知，周期三角信号的傅里叶系数为各谱线之间的间隔为图6-6画出了，和三种情况下傅里叶系数。为了能在同一时间段对比，第2种情况由于周期T不一样，所以谱线之间的间隔也不一样，因此，对横坐标做了调整，使它与第1、3种情况一致。MATLAB源程序为 n=-30:30;tao=1;T=10;w1=2*pi/T; x=n*tao/T;fn=tao*sinc(x); subplot(311) stem(n*w1,fn),grid on title(tao=1,T=10) tao=1;T=5;w2=2*pi/T; x=n*tao/T;fn=tao*(sinc(x)).^2; m=round(30*w1/w2); n1=-m:m; fn=fn(30-m+1:30+m+1); subplot(312) stem(n1*w2,fn),grid on title(tao=1,T=5) tao=2;T=10;w3=2*pi/T; x=n*tao/T;fn=tao*(sinc(x)).^2; subplot(313) stem(n*w3,fn),grid ontitle(tao=2,T=10)图6-6 周期三角信