资源描述：
《上机实验3 连续lti系统的频域分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、上机实验3连续LTI系统的频域分析一、实验目的（1）掌握连续时间信号傅立叶变换和傅里叶逆变换的实现方法，以及傅里叶变换的时移特性，傅立叶变换的频移特性的实现方法。（2）了解傅立叶变换的特点及应用；（3）掌握函数fourier和函数ifourier的调用格式和作用；（4）掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。二、实验原理1.系统的频率特性连续的LTI系统的频率特性又称为频率响应特性，是指系统在正弦信号激励下稳态响应随激励信号频率的变化而变化的情况，又称系统函数H(w)。对于一个零状态的线性
2、系统，如图2.3-1所示Y(w)X(w)H(w)图2.3-1LTI系统框图其系统函数H(w)=Y(w)/X(w)式中，X(w)为系统信号的傅里叶变换，Y(w)为系统在零状态条件下输出响应信号的傅里叶变换。系统函数H(w)反映了系统内在的的固有的特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。H(w)是w的复函数，可以表示为：H(w)=
3、H(w)
4、e^jψ(w)其中，
5、H(w)
6、随w的变化而变化的称为系统的幅频特性；ψ(w)随w变化的规律称为系统的相频特性。
7、频率特性不仅可以用函数表达式表示，还可以用随频率f或者w变化的曲线来描述。当频率特性曲线采用对数坐标表示时，又称为波特图。2.连续时间信号的傅里叶变换的数值计算方法·算法理论依据：F(jw)=当f(t)为限时信号时，或可近似看做限时信号时，上式的n可认为是有限的，记为N则可得F(k)=式中:编程中要注意正确生成信号f(t)的N个样本f(Nt)的向量及向量三、涉及的matlab函数·fourier函数功能：实现信号f(t)的傅里叶变换。调用格式：F=fourier(f):是符号函数f的傅里叶变换，默认返回函
8、数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。四、实
9、验内容与方法1.验证性实验1）编程实现信号的傅里叶变换和傅里叶逆变换（1）傅里叶变换。已知连续时间信号f(t)=e^-2
10、t
11、,通过程序完成f(t)的傅里叶变换。MATLAB程序：symst;f=fourier(exp(-2*abs(t)));ezplot(f);信号f(t)的傅里叶变换图如下：试画出f(t)=2/3*e^-3t*U(t)的波形及其幅频特性曲线。MATLAB程序：symstvwf>>f=2/3*exp(-3*t)*sym('heaviside(t)');>>F=fourier(f);>>s
12、ubplot(2,1,1);>>ezplot(f);>>subplot(2,1,2);>>ezplot(abs(F));信号f(t)的波形及其幅频特性性曲线如上图所示:(2)傅立叶逆变换已知f(jw)=,求信号F（jw）的逆傅立叶变换。MATLAB程序：symstwifourier(1/(1+w^2),t)结果如下：(3)傅立叶变换数值计算已知门函数f(t)=g2(t)=U(t+1)-U(t-1),试采用数值计算方法确定信号的傅立叶变换F（jw）。MATLAB程序：R=0.02;t=-2:R:2;f=st
13、epfun(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
14、);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)傅里叶变化的时移性
