实验4傅立叶变换及其性质..docVIP

傅立叶变换及其性质 一、涉及的MATLAB函数 fourier 功能：实现信号的傅立叶变换。 调用格式： F＝fourier(f)：是符号函数f的傅立叶变换，默认返回函数F是关于的函数。 F＝fourier(f，v)：是符号函数f的傅立叶变换，默认返回函数F是关于v的函数。 F＝fourier(f，u，v)：是关于u的函数f的傅立叶变换，返回函数F是关于v的函数。 ifourier 功能：实现信号的逆傅立叶变换。 调用格式： f＝ifourier(F)：是符号函数F的傅立叶逆变换，默认的独立变换为，默认返回是关于x的函数。 f＝ifourier(F，u)：返回函数F是关于u的函数，而不是默认的x的函数。 f＝ifourier(F，v，u)：是对关于v的函数F的傅立叶逆变换，返回关于u的函数f。 二、实验内容 1、利用fourier函数计算常见信号的傅立叶变换 (1) 已知连续时间信号，通过程序完成信号的傅立叶变换。 MATLAB程序： syms t; f= exp(-2*abs(t)); % 等价于f=sym(‘exp(-2*abs(t))’) Fw=fourier(f); ezplot(Fw); 信号的傅立叶变换如下图所示。 图4.1 信号f(t)的傅立叶变换 （2）试画出单边指数信号的波形及其幅频特性曲线： MATLAB程序如下（a=3）： syms t f=2/3*exp(-3*t)*heaviside(t); -- heaviside 是Symbolic Math Toolbox一个阶跃函数 Fw=fourier(f); subplot(211); ezplot(f); subplot(212); w=(-2:.01:2)*pi; Fw1=subs(Fw, w); %将Fw中在符号变量w用数值来取代 Plot(w,abs(Fw1);grid on 信号的波形及其幅频特性曲线如下图4.2所示。 图4.2信号的波形及其幅频特性曲线图 思考：改变单边指数信号的的大小，观察的变化对频谱的影响？ 2、利用ifourer计算常见信号的傅立叶逆变换 (1) 已知，求信号的逆傅立叶变换。 syms w Fw=1/(1+(w^2)); ft=ifourier(Fw,t) 结果如下： ft =1/2*exp(-t)*Heaviside(t)+1/2*exp(t)*Heaviside(-t) 其中Heaviside(t)是阶跃函数。 已知，求信号的逆傅立叶变换。 syms t w Fw=heaviside(w+10)-heaviside(w-10); ft=ifourier(Fw,t); 结果如下： ft =1/t*sin(10*t)/pi 3、在MATLAB 中验证傅立叶变换的性质 (1) 傅立叶变换的时移特性 分别绘出信号与信号的频谱图，并观察信号时移对信号频谱的影响。 MATLAB程序： %验证傅立叶变换的时移性质 clf; syms t ft=1/2*exp(-2*t)*heaviside(t); subplot(221) ezplot(ft,[-0.5 2]);grid;xlabel(t); title(信号f(t)的波形); ft1=1/2*exp(-2*(t-0.2))*heaviside(t-0.2); subplot(222) ezplot(ft1,[-0.5 2]);grid;xlabel(t); title(信号f(t-0.2)的波形); w=(-4:0.01:4)*pi; Fw=fourier(ft); subplot(223) plot(w,abs(subs(Fw)));grid; xlabel(\omega); title(信号f(t)的幅度频谱); w=(-4:0.01:4)*pi; Fw1=fourier(ft1); subplot(224) plot(w,abs(subs(Fw1)));grid; xlabel(\omega); title(信号f(t-0.2)的幅度频谱); 傅立叶变换的时移特性如下图4.3所示。 图4.3 傅立叶变换的时移特性 (2) 傅立叶变换的频移特性 信号为门函数，试绘出信号频谱图，并与原信号频谱图进行比较。 MATLAB程序： %验证傅立叶变换的频移性质 clear clf; ft=inline(stepfun(t,-1)-stepfun(t,1)); % stepfun函数可以用heaviside 函数代替 t=-3:0.01:3; subplot(221) plot(t,ft(t));grid;xlabel(t); title(信号f(t)的波形); axis([-3 3 -1.5 1.5]); w0=10; ft