实验十七傅立叶变换和傅立叶反变换的数值计算汇编.ppt

* 实验十七 傅里叶变换和傅里叶反变换的数值计算 一、实验目的 1. 掌握实现傅里叶变换和反变换的原理； 2. 会用数值计算方法求给定信号的傅里叶变换和反变换，并绘制相应图形。 二、实验原理 Matlab提供了能直接求解傅里叶变换及反变换的函数fourier( )和ifourier( )。 例如键入命令： syms t fourier(exp(-2*abs(t))) %求 则得到 ans=4/(4+w^2) %F(ω)=4/(4+ ω2) 用数值计算来求信号的傅里叶变换的方法 若信号f(t)是时限的，或当t的绝对值大于某个给定值时，f(t)取值已近似趋向于零，则f(t)可近似看作时限信号，则n的取值个数就是有限的。 设n从n1~n2共取N个值，则 对角频率进行离散化，取离散间隔为Δω 例17-1 画出图17-1所示 f(t)=te-tu(t)的幅度频谱和相位频谱。 R=0.01; t=0:R:10; %忽略掉f(t)在t10之后的非零取值 f=t.*exp(-t).*u(t); W1=2*pi; K=500; k=-K:K; w=k*W1/K; %生成角频率的矢量 F=f*exp(-j*t*w)*R; %计算F(ω) absF=abs(F); angF=angle(F); %幅度谱和相位谱 subplot(1,2,1), plot(w, absF) axis([-2*pi, 2*pi, 0, max(abs(F))+0.1]), grid on set(gca, xtick, [-1 0 1]), set(gca, ytick, [1/2 1]) subplot(1,2,2), plot(w, angF); axis([-2*pi, 2*pi, -pi, pi]), grid on set(gca, xtick, [-1 0 1]), set(gca, ytick, [ -pi -pi/2 0 pi/2 pi]) t 1 3 1/e 0 例17-2 画出F(ω)=2Sa(ω)的傅里叶反变换f(t)的波形。 R=0.02; w=-10*pi:R:10*pi; %忽略掉界外非零取值 F=2*sin(w)./w; N=100; n=-N: N; t=n*2/N; %生成t在区间[-2,2]上的矢量 f=F*exp(j*w*t)*R/(2*pi); %计算f(t) f=real(f); subplot(1,2,1), plot(w, F) axis([-3.5*pi, 3.5*pi, -0.6, 2.2]), grid on set(gca, xtick, [-3*pi -pi 0 pi 3*pi]), set(gca, ytick, [0 2]) subplot(1,2,2), plot(t, f) axis([-2,2,-0.2,1.2]), grid on set(gca, xtick, [-1 0 1]), set(gca, ytick, [0 1]) 说明：按傅里叶变换的定义， F(ω)=2Sa(ω)的傅里叶反变换应为 f(t)=u(t+1)-u(t-1)，这与图17-3画出的波形有所不同。 原因：其原因在于本题程序只取了 F(ω)在-10π ≤ ω ≤ 10π区间内的部分作反变换，这相当于将 F(ω)通过一个截止频率为 ωc= 10π的低通滤波器，由此导致反变换结果 f(t)偏离矩形脉冲而呈现图17-3所示的形状，这就是著名的吉布斯现象。 加大低通滤波器的截止频率ωc以允许更多的高频分量通过可使 f(t)更接近矩形脉冲 。 三. 实验内容 1. 利用Matlab画出教材习题3-26所示梯形脉冲信号在t2=2t1情况下的频谱图，并与手工绘制草图比较。 t2/2 t1/2 -t1/2 -t2/2 0 E *