matlab对频谱进行周期重复,《实验1信号的频谱图.docx

《实验1信号的频谱图

实验一信号的频谱图实验目的1. 掌握周期信号的傅里叶级数展开2. 掌握周期信号的有限项傅里叶级数逼近3. 掌握周期信号的频谱分析4. 掌握连续非周期信号的傅立叶变换5. 掌握傅立叶变换的性质相关知识1 周期信号的傅里叶级数设周期信号，其周期为T，角频率为，该信号可展开为三角形式的傅里叶级数，即为：其中，正弦项与余弦项的系数和成为傅里叶系数，根据函数的正交性，得(2)其中，。积分区间通常取为或。若将(2)式中同频率项合并，可改写为(3)从物理概念上来说，(3)中的即是信号的直流分量；式中的第二项称为信号的基波或者基波分量，它的角频率与原周期信号相同；式中第三项称为信号的二次谐波，他的频率是基波频率的二倍；以此类推。一般而言称为信号的n次谐波；n比较大的分量统称为信号的高次谐波。我们还常用到复指数形式的傅里叶。设周期信号，其周期为T，角频率为，该信号复指数形式的傅里叶级数为其中，称为复指数形式傅里叶级数系数。利用MATLAB可以直观地观察和分析周期信号傅里叶级数及其收敛性。【例1-1】周期方波信号如图所示，画出该信号的傅里叶级数，利用MATLAB编程实现其各次谐波的叠加。图1-1 周期方波信号解：从理论上分析可知，已知周期方波信号的傅里叶级数展开为取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:N n=1:2: n_max(k);b=4./(pi*n);x=b*sin(omega*n'*t);figure;plot(t,y,'b');hold on;plot(t,x,'r');hold off;xlabel('t'),ylabel('部分和的波形');axis([-1 1 -1.5 1.5]);grid on;title(['最大谐波数=',num2str(n_max(k))])end运行后的各项部分和的波形如图图1-2 周期方波信号的有限项傅里叶级数逼近很多项的时候，部分和的波形和周期方波信号的波形很接近，但在信号的跳变点附近，却总是存在一个过冲，这就是所谓的Gibbs现象。2 周期信号的频谱分析周期信号通过傅里叶级数分解可展开成一些列相互正交的正弦信号或复指数信号分量的加权和。在三角形是傅里叶级数中，各分量的形式为；在指数形式的傅里叶级数中，各分量的形式为。对实信号而言，和成对出现。对不同的周期信号，它们各个分量的数目、角频率、幅度或、相位或不同。傅里叶系数的幅度或随角频率的变化关系绘制成图形，称为信号的幅度频谱，简称幅度谱。相位或随角频率的变化关系绘制成图形，称为信号的相位频谱，简称相位谱。幅度谱和相位谱统称为信号的频谱。信号的频谱是信号的另一种表示，它提供了从另一个角度来观察和分析信号的途径。利用MATLAB命令可对周期信号的频谱及其特点进行观察验证和分析。【例1-2】已知周期矩形脉冲如图所示，设脉冲幅度为A=1，宽度为，重复周期为T(角频率)。将其展开为复指数形式傅里叶级数，研究周期矩形脉冲的宽度和周期变化时，对其频谱的影响。图1-3 周期矩形脉冲信号解：根据傅里叶级数理论可知，周期矩形脉冲信号的傅里叶系数为各谱线之间的间隔为。图画出了；和三种情况下傅里叶系数。为了能在同一时间段对比，第二种情况由于周期T不一样，所以谱线之间的间隔也不一样，因此对横坐标进行了调整，使它与第一种和第三种情况一致。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);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);subplot(313)stem(n*w3,fn),grid on;title('tao=2,T=10');图1-4 周期矩形脉冲信号的傅里叶系数从图中可以看出，脉冲宽度越大，信号的频谱带宽越小；而周期越小，谱线之间间隔越大，验证了傅里叶级数理论。【练习】已知周期三角信号如图所示，试求出