傅 立 叶 变 换 五 大 性 质 的 matlab 实 现 20092426 2012-5-10 xx 远 整 理 一 . 傅 立 叶变换的 时移 性 质 若)()(Ftf, 则00)(0)()()(tjtjeFeFttf 结论: )(tf延 时( 或 超 前 )0t 后 , 其 对应的 幅 度 谱保 持 不 变, 但 相 位谱中 一 切 频率 分 量 的 相 位 均 滞后 ( 或 超 前 )0t。 例1 ( 1) 用matlab 画)(21)(2tetft 及 频谱( 幅 度 谱及 相 位谱) ( 2) 用 matlab 画)5.0(1tff及 频谱( 幅 度 谱及 相 位 谱)。 ( 1) 程 序 : N=256;t=linspace(-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=linspace(-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=angle(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)相 位 谱') ...