LTI的频域分析matlab,matlab与信号实验连续lti系统的频域分析

matlab与信号实验连续lti系统的频域分析

上机实验 3 连续 LTI 系统的频域分析一．实验目的 (1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法;(2).了解傅立叶变换的频移特性及其应用;(3).掌握函数 fourier 和函数 ifourier 的调用格式及作用;(4).掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。2．实验原理1.系统的频率特性连续 LTI 系统的频率特性又称频率响应特性，是指系统在正弦信号激励下的稳态响应随频率变化的情况，又称系统函数 。对于一个零状态H的线性系统，如图 2.3-1 所示。其系统函数 定义为 式中， 为系统激励信号H)()(jXYj的傅里叶变换， 为系统在零状态条件下输出响应的傅里叶变换。Y系统函数 反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。 是 的复函数，可以表示为 = 。其中，HHej随 变化的规律称为系统的幅频特性： 随 变化的规律称为系统的相频特性。频率特性不仅可用函数表达式表示，还可以随频率 f 变化的曲线来表示。当频率特性曲线采用对数坐标时，又称为波特图。)(jX)(j )(jY图 2.3-1 LTI 系统框图2.连续时间信号傅里叶变换的数值计算方法算法理论依据： (2.2-eenjtj fdfjF)(lim)(1)当 为时限信号时，或和近似的看做时限信号时，式(2-1)中的 n 取)(tf值可认作是有限的，设为 N，则可得 ，ejNnkfkF100<=k<=N (2.2-2) 式(2.2-2 )中 。编程中需要注意的是：要正k2确生成信号 的 N 个样本 的向量及向量 。)(tf)(nf njk3．涉及的 MATLAB 函数1、.fourier 函数功能：实现信号 的傅里叶变换。)(tf调用格式:F=fourier(f):符号函数 f 的傅里叶变换，默认返回函数 F 是关于 的函数。F=fourier(f,v)：是符号函数 f 的傅里叶变换，返回函数 F 是关于 v 的函数。F=fourier(f,u,v):是关于 u 的函数 f 的傅里叶变换，返回函数 F 式关于 v的函数。2、.ifourier 函数功能:实现信号 的傅里叶变换)(jF调用格式:F=ifourier(F):是函数 F 的傅里叶反变换，默认的独立变量 ，默认返回是关于 x 的函数。F=ifourier(F，v)：返回函数 f 是 u 的函数，而不是默认的 x 的函数。F=ifourier(F,v,u)：是对关于 v 的函数 F 进行傅里叶逆变换，返回关于 u的函数 f。四、实验内容与方法1. .验证性实验1) 编程实现信号的傅立叶变换和傅立叶逆变换(1) 傅立叶变换。已知连续时间信号 f(t)= ,通过程序完成 f(t)的傅立叶变换。et|2MATLAB 程序：syms t; f=fourier(exp(-2*abs(t))); ezplot(f);信号 f(t)的傅立叶变换如图。试画出 f(t)= U(t)的波形及其幅频特性曲线。et32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)的波形及其幅频特性曲线如图。et32(2) 傅立叶逆变换已知 f(jw)= ,求信号 F(jw)的逆傅立叶变换。w21MATLAB 程序：syms t wifourier(1/(1+w^2),t)结果如下：ans= *exp(-t)*U(t)+ *exp(t)*U(t)2121(3) 傅立叶变换数值计算已知门函数 f(t)=g2(t)=U(t+1)-U(t-1),试采用数值计算方法确定信号的傅立叶变换 。jFMATLAB 程序：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=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)];Subpolt(2,1,1);plot(t,f);Xlabel(‘t’);ylabel(‘f(t)’);axis([-2,2,-o.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)傅里叶变换‘);信号的傅里叶变换如图所示。(4) 连续函数的傅立叶变换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;连续函数的傅立叶变换如图。2. 程序设计实验(1) 试确定下列信号的傅立叶变换的数学表达式。(a)f(t)=U(t+1)-U(t-1)MATLAB 程序：syms t ωfourier(-heaviside(t-1)+heaviside(t+1))结果如下：ans=(1/exp(ω*j))*(- pi*dirac(-ω) + j/ω) - exp(ω*j)*(- pi*dirac(- ω) + j/ω)(b)f(t)=e^-3t*U(t)MATLAB 程序：fourier(exp(-3*t)*heaviside(t));结果如下：ans=1/(3 + w*i)(c)f(t)=e^-t*U(t)MATLAB 程序：f=fourier(exp(-1*t)*heaviside(t));结果如下：ans=1/(1 + w*i)(d)f(t)=σ”*U(t)MATLAB 程序：syms tfourier(diff(diff(dirac(t)))*heavisid