(西北工业大学信号系统上机实验3.docVIP

上机实验3 连续LTI系统的频域分析 实验目的 掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法; 了解傅立叶变换的频移特性及其应用 掌握函数fourier和函数ifourier的调用格式及作用; 掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。 实验内容 .验证性实验 编程实现信号的傅立叶变换和傅立叶逆变换 傅立叶变换。 已知连续时间信号f(t)=,通过程序完成f(t)的傅立叶变换。 MATLAB程序： syms t; f=fourier(exp(-2*abs(t))); ezplot(f); 信号f(t)的傅立叶变换如图。 试画出f(t)=U(t)的波形及其幅频特性曲线。 MATLAB程序： syms t v w f f=2/3*exp(-3*t)*sym(Heaviside(t)); F=fourier(f); subplot(2,1,1); ezplot(f); subplot(2,1,2); ezplot(abs(F)); 信号f(t)=U(t)的波形及其幅频特性曲线如图 傅立叶逆变换 已知f(jw)=,求信号F（jw）的逆傅立叶变换。 MATLAB程序： syms t w ifourier(1/(1+w^2),t) 结果如下： 傅立叶变换数值计算 已知门函数f(t)=g2(t)=U(t+1)-U(t-1),试采用数值计算方法确定信号的傅立叶变换F（jw）。 MATLAB程序： R=0.02; t=-2:R:2; f=stepfun(t,-1)-stepfun(t,1); W1=2*pi*5; N=500; k=0:N; W=k*W1/N; F=f*exp(-j*t*W)*R; F=real(F); W=[-fliplr(W),W(2:501)]; F=[fliplr(F),F(2:501)]; subplot(2,1,1);plot(t,f); xlabel(t);ylabel(f(t));axis([-2,2,-1.5,2]); title(f(t)=U(t+1)-U(t-1)); subplot(2,1,2);plot(W,F); xlabel(w);ylabel(F(w));title(f(t)的傅立叶变换); 信号的傅立叶变换如图。 连续函数的傅立叶变换 MATLAB程序： clf; dt=2*pi/8; w=linspace(-2*pi,2*pi,2000)/dt; k=-2:2;f=ones(1,5);F=f*exp(-j*k*w); f1=abs(F);plot(w,f1);grid; 连续函数的傅立叶变换如图。 程序设计实验 试确定下列信号的傅立叶变换的数学表达式。 f(t)=U(t+1)- U(t-1) f(t)=U(t) f(t)=U(t) 试画出信号f(t)=U(t), f(t-4)以及信号f(t)的频谱图。 f(t)=U(t) f(t-4) f(t)