MATLAB与信号实验 连续LTI系统的频域分析.doc

上机实验3 连续LTI系统的频域分析 一．实验目的 (1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法; (2).了解傅立叶变换的频移特性及其应用; (3).掌握函数fourier和函数ifourier的调用格式及作用; (4).掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。 实验原理 系统的频率特性 连续LTI系统的频率特性又称频率响应特性，是指系统在正弦信号激励下的稳态响应随频率变化的情况，又称系统函数。对于一个零状态的线性系统，如图2.3-1所示。 其系统函数定义为式中，为系统激励信号的傅里叶变换，为系统在零状态条件下输出响应的傅里叶变换。 系统函数反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。是的复函数，可以表示为=。其中，随变化的规律称为系统的幅频特性：随变化的规律称为系统的相频特性。 频率特性不仅可用函数表达式表示，还可以随频率f变化的曲线来表示。 当频率特性曲线采用对数坐标时，又称为波特图。 连续时间信号傅里叶变换的数值计算方法 算法理论依据：（2.2-1） 当为时限信号时，或和近似的看做时限信号时，式（2-1）中的n取值可认作是有限的，设为N，则可得，0=k=N (2.2-2) 式（2.2-2）中。编程中需要注意的是：要正确生成信号的N个样本的向量及向量。 涉及的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，v)：返回函数f是u的函数，而不是默认的x的函数。 F=ifourier(F,v,u）：是对关于v的函数F进行傅里叶逆变换，返回关于u的函数f。 四、实验内容与方法 .验证性实验 编程实现信号的傅立叶变换和傅立叶逆变换 傅立叶变换。 已知连续时间信号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) 结果如下： ans=*exp(-t)*U(t)+*exp(t)*U(t) 傅立叶变换数值计算 已知门函数f(t)=g2(t)=U(t+1)-U(t-1),试采用数值计算方法确定信号的傅立叶变换。 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程序： R=0.02;t=-2:R:2; F=stepfun(t,-1)-stepfun(t,1); W1=