实验7-8连续LTI时间系统的频域及S域分析.docVIP

河北工业大学通信091 杨泽祺 092295 实验七 连续 LTI 时间系统的频域分析 一·实验目的： 了解傅里叶变换、傅里叶逆变换的实现方式，以及傅立叶变换的时移特性、傅立叶 变换的频移特性的实现方式。 二．实验内容： (1) 了解如何采用函数 fourier()对时域信号进行傅里叶变换，以及如何采用函数 ifourier()对频域信号进行傅里叶反变换。 (2) 了解傅里叶变换的数值计算方法，以及信号频谱图。 (3) 了解傅里叶变换的时移特性、频移特性的实现方法。 （（1傅里叶变换 使用函数 fourier()可以实现信号 f(t)的傅里叶变换，调用格式如下所示： F=fourier(f) 是符号函数 f的傅里叶变换，默认返回函数 F是关于ω的函数。 F=fourier(f,v) 是符号函数 f的傅里叶变换，返回函数 F是关于 v 的函数。 F=fourier(f,u,v) 是关于 u的函数 f的傅里叶变换，返回函数 F是关于 v 的函数 （（2 傅里叶逆变换 f=ifourier(F) 是函数 F的傅里叶逆变换，默认的独立变量为ω，默认返回是关于 x 的函数。 f=ifourier(F,u) 返回函数f是 u 的函数，而不是默认的 x 的函数。 f=ifourier(F,v,u) 是对关于 v 的函数F进行傅里叶逆变换，返回关于 u 的函数 f。 ）） 三、实验结果 1：试求下列信号的傅里叶变换的数学表达式。 （1） f (t) = ε (t +1) ?ε (t ?1) syms t f f=heaviside(t+1)-heaviside(t-1); F=fourier(f) 结果如下： F = exp(i*w)*(pi*dirac(w) - i/w) - (pi*dirac(w) - i/w)/exp(i*w)（2） f (t) = e^(?3t)*ε (t) syms t f f=exp(-3*t)*heaviside(t); F=fourier(f) 结果如下： F = 1/(3+i*w) （3） f (t) = e^(?t)*ε (t) syms t f f=exp(-t)*heaviside(t); F=fourier(f) 结果如下： F = 1/(1+i*w)’’ (t) syms t f f=diff(diff(diff(heaviside(t)))); F=fourier(f) 结果如下： F = -w^2 2：试画出信号f (t) = e^(?3t)*ε (t) 的频谱图，并画出信号f (t ? 4) 以及信号f (t)e^(? j4t)的频谱图。f (t) = e^(?3t)*ε (t) 的频谱图 r=0.02; t=-5:r:5; N=200; W=2*pi; k=-N:N; w=k*W/N; f1=exp(-3*t).*eaviside(t); F=r*f1*exp(-j*t*w); F1=abs(F); P1=angle(F); subplot(3,1,1) plot(t,f1) grid xlabel(t) ylabel(f(t)) title(f(t)) subplot(3,1,2) plot(w,F1); xlabel(w) grid ylabel(F(jw)); subplot(3,1,3) plot(w,P1*180/pi); grid xlabel(w); ylabel(度) 2、f (t ? 4)的频谱图 r=0.02; t=-5:r:5; N=200; W=2*pi; k=-N:N; w=k*W/N; f1=exp(-3*(t-4)).*heaviside(t-4); F=r*f1*exp(-j*t*w); F1=abs(F); P1=angle(F); subplot(3,1,1) plot(t,f1) grid on xlabel(t) ylabel(f(t)) title(f(t-4)) subplot(3,1,2) plot(w,F1); xlabel(w) grid on ylabel(F(jw)的模); subplot(3,1,3) plot(w,P1*180/pi); grid xlabel(w); ylabel(相位（度）) f (t)e^(? j4t)的频谱图r=0.02; t=-2:r:2 N=500; W=2*pi*5; k=-N:N; w=k*W/N; f=exp(-3*t).*heaviside(t); f1=f.*exp(j*4*t); F=