连续LTI系的频域分析.doc

上机实验3 连续LTI系统的频域分析 实验目的 掌握连续时间信号傅立叶变换和傅里叶逆变换的实现方法，以及傅里叶变换的时移特性，傅立叶变换的频移特性的实现方法。 了解傅立叶变换的特点及应用； 掌握函数fourier和函数ifourier的调用格式和作用； 掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。 二、实验原理 1.系统的频率特性 连续的LTI系统的频率特性又称为频率响应特性，是指系统在正弦信号激励下稳态响应随激励信号频率的变化而变化的情况，又称系统函数H(w)。对于一个零状态的线性系统，如图2.3-1所示 Y(w)X(w) H(w) Y(w) X(w) H(w) 图 2.3-1 LTI 系统框图 其系统函数H(w)=Y(w)/X(w)式中，X(w)为系统信号的傅里叶变换，Y(w)为系统在零状态条件下输出响应信号的傅里叶变换。 系统函数H(w)反映了系统内在的的固有的特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。H(w)是w的复函数，可以表示为： H(w)=|H(w)|e^jψ(w) 其中，|H(w)|随w的变化而变化的称为系统的幅频特性；ψ(w)随w变化的规律称为系统的相频特性。频率特性不仅可以用函数表达式表示，还可以用随频率f或者w变化的曲线来描述。当频率特性曲线采用对数坐标表示时，又称为波特图。 2. 连续时间信号的傅里叶变换的数值计算方法 算法理论依据： F(w)=dt=当f(t)为限时信号时，或可近似看做限时信号时，上式的n可认为是有限的，记为N 则可得 F(k)=tt 式中=2π/(Nt)*k 编程中要注意正确生成信号f(t)的N个样本f(Nt)的向量及向量 三、涉及的matlab函数 fourier函数 功能：实现信号f(t)的傅里叶变换。 调用格式： F=fourier(f):是符号函数f的傅里叶变换，默认返回函数F是关于w的函数； F=fourier(f,v)：是符号函数f的傅里叶变换，返回函数F是关于v的函数。 F=fourier(f,u,v)：是关于u的函数的f的傅里叶变换，返回函数F是关于v的函数。 ifourier函数 功能：实现信号F(jw)的傅里叶逆变换。 F=ifourier(F)：是函数F的傅里叶逆变换，默认返回函数F是关于x的函数； F=ifourier(F,v):返回函数f是v的函数，而不是关于x的函数； F=ifourier(F,v,u)：是对关于v的函数F进行傅里叶逆变换，返回关于u的函数f。 四、实验内容与方法 1.验证性实验 （1）傅里叶变换。 已知连续时间信号f(t)=e^-2|t|,通过程序完成f(t)的傅里叶变换。 MATLAB程序： syms t; f=fourier(exp(-2*abs(t))); ezplot(f); 信号f(t)的傅里叶变换图如下： （2）试画出f(t)=2/3*e^-3t*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)的波形及其幅频特性性曲线如图: (3)傅里叶变化的时移性 分别绘出信号f(t)=1/2*e^-2t*U(t)与信号f(t-1)的频谱图，并观察信号时移对信号频谱的影响。 MATLAB的程序： r=0.02;t=-5:r:5;N=200;Y=2*pi;k=-N:N;w=k*Y/N; f1=1/2*exp(-2*t).*stepfun(t,0); 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(相位（度）); 傅里叶变换的时移特性如下图： (4)傅里叶变化的频移特性 信号f(t)=g2(t)为门函数，试绘出f1(t)=f(t)*e^-j10以及信号f2(t)=f(t)*e^j10t的频谱图，并与原信号频谱图进行比较。 MATLAB程序： . xlabel(w);ylab