傅里叶变换(FT)与其性质.doc

实验四 傅里叶变换(FT)及其性质 一、实验目的 1、学会运用Matlab求连续时间信号的傅里叶 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 二．实验内容 clear all;close all;clc ft=sym(sin(2*pi*(t-1))/(pi*(t-1))); fw=fourier(ft) ftt=sym((sin(pi*t)/(pi*t))^2); fww=fourier(ftt) subplot(2,2,1) ezplot(abs(fw)),grid on title(幅度谱a) phase=atan(imag(fw)/real(fw)); subplot(2,2,2) ezplot(phase),grid on title(相位谱A) subplot(2,2,3) ezplot(abs(fww)),grid on title(幅度谱b) phasee=atan(imag(fww)/real(fww)); subplot(2,2,4) ezplot(phasee),grid on title(相位谱B) 2、clear all;close all;clc syms t fw1=sym(10/(3+iw)-4/(5+iw)); ft1=ifourier(fw1,t) fw2=sym(exp(-4*w*w)); ft2=ifourier(fw2,t) subplot(2,1,1) ezplot(t,ft1,[-15,15]),grid on title(时域信号a) subplot(2,1,2) ezplot(t,ft2,[-10,10]),grid on title(时域信号b) [注意：(1)写代码时j(i] 3、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT，并画出其频谱图。 clear all;close all;clc; ft=sym((t+2)*(heaviside(t-2)-heaviside(t-1))+heaviside(t-1)-heaviside(t+1)+(-t+2)*(heaviside(t-2)-heaviside(t-1))); fw=fourier(ft) subplot(2,1,1); ezplot(abs(fw));grid on; title(幅度谱); phase=atan(imag(fw)/real(fw)); subplot(2,1,2); ezplot(phase);grid on; title(相位谱); fw = i/w*(-4*exp(-2*i*w)+3*exp(-i*w)+exp(i*w))  4、已知门函数自身卷积为三角波信号，试用Matlab命令验证FT的时域卷积定理。 4. clear all;close all;clc; syms t ft=sym(heaviside(t-1)-heaviside(t+1)); fw=fourier(ft) f0=ft*ft; f1=fourier(f0) f2=fw*fw subplot(2,1,1); ezplot(f1);grid on; title(时域卷积); subplot(2,1,2); ezplot(f2);grid on; title(频域相乘); fw = -2/w*sin(w) f1 = 2/w*sin(w) f2 = 4/w^2*sin(w)^2 实验思考： 感觉第三个题目挺麻烦的，得多看例题。Matlab的功能很强大，但是太多的东西都忘记了，课下得好好地看Matlab。