傅立叶变换五大性质的matlab实现.docx

傅立叶变换五大性质的matlab实现2012-5-10 xx 远整理 一. 傅立叶变换的时移性质 若 f(t) F()，则 f(t t0) F ( )e j t0 |F( )ej ( ) t0 结论：f (t)延时(或超前)t0后，其对应的幅度谱保持不变，但相 位谱中一切频率分量的相位均滞后(或超前) t0。 例 1 ( 1)用 matlab 画 f (t) (2)用 matlab 画 f 1 1 -e 2 f(t 2t (t)及频谱(幅度谱及相位谱) 0.5)及频谱(幅度谱及相位谱)。 (1)程序： N=256;t=li nspace(-2,2,N); % 进行时间分割，在【-2，2】 内均匀产生N点，分割成N-1段 f=1/2*exp(-2*t).*heaviside(t); % 建立信号f(t)，这里点乘 ‘ .* 不能用*，点乘是对应元 素相乘，‘ * 是矩阵相乘。 dt=4/(N-1); % 时间长度为4,均匀分割成N-1 段，相邻两时间点的间隔为dt M=401; w=li nspace(-2*pi,2*pi,M); %进行频率分割，在 [-2*pi,2*pi] 内均匀产生M点， 分割成M-1段 F=f*exp(-j*t*w)*dt; % 求信号f(t)的傅立叶变换 F1=abs(F);P1=a ngle(F); % 求幅度谱和相位谱 subplot(3,1,1); plot(t,f);grid on xlabel(t);ylabel(f(t)); title(f(t)) subplot(3,1,2); plot(w,F1);grid on xlabel(w);ylabel(abs(F(w))); subplot(3,1,3); plot(w,P1);grid on xlabel(w);ylabel(angle(F(w))); (2)程序： N=256; t=linspace(-2,2,N); f=1/2*exp(-2*t).*heaviside(t); % 建立时间信号 f(t) f1=1/2*exp(-2*(t-0.5)).*heaviside(t-0.5); %建立 时间信 号 f(t-0.3) dt=4/(N-1); M=401; w=linspace(-2*pi,2*pi,M); F=f*exp(-j*t*w)*dt; % 求信号 f(t) 的傅立叶变换 F1=f1*exp(-j*t*w)*dt; % 求信号 f(t-0.5) 的傅立叶变换 subplot(3,1,1); plot(t,f,t,f1,r),grid on xlabel(t);ylabel(f), title(f(t),f(t-0.5)) subplot(3,1,2); plot(w,abs(F),w,abs(F1),r),grid on xlabel(w); ylabel( f(t) 和 f(t-0.5) 幅度谱 ) ； subplot(3,1,3); plot(w,angle(F),w,angle(F1),r),grid on xlabel(w); ylabel( f(t) 和 f(t-0.5) 相位谱 ) 二．傅立叶变换的频移性质 若 f (t) F( )，则 f(t) ej 0t F( 0) 结论：将信号 f ( t )乘以因子 ej 0t ，对应于将频谱函数沿轴 右移 0； 将信号 f ( t )乘以因子 ej 0t ，对应于将频谱函数沿轴 右移 0。 例 2 已知 f (t) (t 1) (t 1)， 且 f1(t) f(t) e20jt, f2(t) f (t) e20jt，求： ( 1)用 matlab 在同一个图中画它们的幅度谱； (2)用 matlab 在同 一个图中画它们的幅度谱的实部；验证傅立叶变换的频移特性 程序： N=256;M=500; t=linspace(-2,2,N); w=linspace(-10*pi,10*pi,M); % 在 [-10*pi,10*pi] 内进行频率分 dt=4/(N-1); f=heaviside(t+1)-heaviside(t-1); f1=f.*exp(j*20*t); f2=f.*exp(-j*20*t); % 这里必须用‘ .* F=f*exp(-j*t*w)*dt; % 求 f(t) 的傅立叶变换 F1=f1*exp(-j*t*w)*dt; F2=f2*exp(-j*t*w)*dt; subplot(2,1,1); plot(w,real(F),w,real(F1),r,w,real(F2),g), grid on xlabel(w);ylabel(real(F(w))); title( 信号傅立叶变换的实部 ) subplot(2,1,2);