资源描述：
《实验二-LTI系统的频域分析.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验2：——连续LTI系统的频域特性及频域分析实验性质：提高性实验级别：必做开课单位：机械电子工程学院学时：2一、实验目的1、掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。2、掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。二、实验设备计算机，MATLAB软件三、实验原理1、连续时间LTI系统的频率响应所谓频率特性，也称为频率响应特性，简称频率响应（Frequencyresponse），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。
2、图1连续时间LTI系统的时域及频域分析图图1中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：，由傅里叶变换的时域卷积定理可得到：(1)或者：(2)为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即(3)由于H(jw)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutlyintegrabel）的话，那么H(jw)一定存在，而且H(jw)通常是复数，因此，也可以表示成复数的不同表达形式
3、。在研究系统的频率响应时，更多的是把它表示成极坐标形式：(4)上式中，称为幅度频率响应（Magnituderesponse），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，称为相位特性（Phaseresponse），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。和都是频率w的函数。对于一个系统，其频率响应为H(jw)，其幅度响应和相位响应分别为和，如果作用于系统的信号为，则其响应信号为2、用MATLAB计算系统频率响应在本实验中，表示系统的方法仍然是用频响函数分子和分母多项式系数行向量来表示
4、。实验中用到的MATLAB函数如下：[H,w]=freqs(b,a)：b,a分别为连续时间LTI系统的微分方程右边的和左边的系数向量（Coefficientsvector），返回的频率响应在各频率点的样点值（复数）存放在H中，系统默认的样点数目为200点；计算频率响应的函数freqs()的另一种形式是：H=freqs(b,a,w)：在指定的频率范围内计算系统的频率响应特性。在使用这种形式的freqs/freqz函数时，要在前面先指定频率变量w的范围。例如在语句H=freqs(b,a,w)之前加上语句：w=0:2*p
5、i/256:2*pi。Hm=abs(H)：求模数，即进行运算，求得系统的幅度频率响应，返回值存于Hm之中。real(H)：求H的实部；imag(H)：求H的虚部；phi=angle(H)：求相位频率响应特性；tao=grpdelay(b,a,w)：计算系统的相位频率响应所对应的群延时。下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。假设给定一个连续时间LTI系统，下面的微分方程描述其输入输出之间的关系编写的MATLAB范例程序，绘制系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚
6、部。程序如下：%ThisProgramisusedtocomputeanddrawtheplotsofthefrequencyresponse%ofacontinuous-timesystemb=[1];%Thecoefficientvectoroftherightsideofthedifferentialequationa=[132];%Thecoefficientvectoroftheleftsideofthedifferentialequation[H,w]=freqs(b,a);%Computethefre
7、quencyresponseHHm=abs(H);%ComputethemagnituderesponseHmphai=angle(H);%ComputethephaseresponsephaiHr=real(H);%ComputetherealpartofthefrequencyresponseHi=imag(H);%Computetheimaginarypartofthefrequencyresponsesubplot(221)plot(w,Hm),gridon,title('Magnituderesponse
8、'),xlabel('Frequencyinrad/sec')subplot(223)plot(w,phai),gridon,title('Phaseresponse'),xlabel('Frequencyinrad/sec')subplot(222)plot(w,Hr),gridon,title('Realpartoffrequencyresponse')
