信号与系统实验报告实验三---连续时间LTI系统的频域分析

实验三连续时间LTI系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用MATLAB语言进行系统频响特性分析的方法。基本要求：掌握LTI连续和离散时间系统的频域数学模型和频域数学模型的MATLAB描述方法，深刻理解LTI系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB计算和绘制LTI系统频率响应特性曲线中的编程。二、实验原理及方法1连续时间LTI系统的频率响应所谓频率特性，也称为频率响应特性，简称频率响应（Frequencyresponse），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(thtxty，由傅里叶变换的时域卷积定理可得到：)()()(jHjXjY3.1或者：)()()(jXjYjH3.2)(jH为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即dtethjHtj)()(3.3由于H(j)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutlyintegrabel）的话，那么H(j)一定存在，而且H(j)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极连续时间LTI系统的时域及频域分析图系统LTI)(th)(jH)(ty)(jX)(jY)(tx2坐标形式：)()()(jejHjH3.4上式中，)j(H称为幅度频率相应（Magnituderesponse），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(称为相位特性（Phaseresponse），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)(jH和)(都是频率的函数。对于一个系统，其频率响应为H(j)，其幅度响应和相位响应分别为)(jH和)(，如果作用于系统的信号为tjetx0)(，则其响应信号为tjejHty0)()(0tjjeejH00)(0)())((000)(tjejH3.5若输入信号为正弦信号，即x(t)=sin(0t)，则系统响应为))(sin(|)(|)sin()()(00000tjHtjHty3.6可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)(jH加权，二是信号的相位要被)(移相。由于)(jH和)(都是频率的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。2LTI系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weightedsum）所组成。正如刚才所述，信号经过LTI系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phaseshifting），相位移实际上就是延时（Timedelay）。群延时（Groupdelay）的概念能够较好地反映系统对不同频率分量造成的延时。LTI系统的群延时定义为：dd)()(3.7群延时的物理意义：群延时描述的是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说，该系统对于所有的频率分量造成的延时时间都是一样的，因而，系统不会对信号产生相位失真（Phasedistortion）。反之，若系统的相位频率响应特性不是线性的，则该系统对于不同频率的频率分量造成的延时3时间是不同的，因此，当信号经过系统后，必将产生相位失真。3用MATLAB计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。实验中用到的MATLAB函数如下：[H,w]=freqs(b,a)：b,a分别为连续时间LTI系统的微分方程右边的和左边的系数向量（Coefficientsvector），返回的频率响应在各频率点的样点值（复数）存放在H中，系统默认的样点数目为200点；Hm=abs(H)：求模数，即进行HHm运算，求得系统的幅度频率响应，返回值存于Hm之中。real(H)：求H的实部；imag(H)：求H的虚部；phi=atan(-imag(H)./(real(H)+eps))：求相位频率相应特性，atan()用来计算反正切值；或者phi=angle(H)：求相位频率相应特性；tao=grpdelay(num,den,w)：计算系统的相位频率响应所对应的群延时。计算频率响应的函数freqs()的另一种形式是：H=freqs(b,a,w)：在指定的频率范围内计算系统的频率响应特性。在使用这种形式的freqs/freqz函数时，要在前面先指定频率变量w的范围。例如在语句H=freqs(b,a,w)之前加上语句：w=0:2*pi/256:2*pi。下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。假设给定一个连续时间LTI系统，下面的微分方程描述其输入输出之间的关系)()(2)(3)(22txtydttdydttyd编写的MATLAB范例程序，绘制系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部。程序如下：%Program3_1%ThisProgramisusedtocomputeanddrawtheplotsofthefrequencyresponse%ofacontinuous-timesystemb=[1];%Thecoefficientvectoroftherightsideofthedifferentialequationa=[132];%Thecoefficientvectoroftheleftsideofthedifferentialequation[H,w]=freqs(b,a);%ComputethefrequencyresponseHHm=abs(H);%ComputethemagnituderesponseHmphai=angle(H);%ComputethephaseresponsephaiHr=real(H);%ComputetherealpartofthefrequencyresponseHi=imag(H);%Computetheimaginarypartofthefrequencyresponsesubplot(221)plot(w,Hm),gridon,title('Magnituderesponse'),xlabel('Frequencyinrad/sec')subplot(223)plot(w,phai),gridon,title('Phaseresponse'),xlabel('Frequencyinrad/sec')4subplot(222)plot(w,Hr),gridon,title('Realpartoffrequencyresponse'),xlabel('Frequencyinrad/sec')subplot(224)plot(w,Hi),gridon,title('Imaginarypartoffrequencyresponse'),xlabel('Frequencyinrad/sec')三、实验内容及步骤实验前，必须首先阅读本实验原理，了解所给的MATLAB相关函数，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。给定三个连续时间LTI系统，它们的微分方程分别为系统1：dttdxtydttdydttyd)()(25)(1)(22Eq.3.1系统2：)()()()(txdttdxtydttdyEq.3.2系统3：)(262)(262)(401)(306)(148)(48)(10)(2233445566txtydttdydttyddttyddttyddttyddttydEq.3.3Q3-1修改程序Program3_1，并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程序Q3_1如下：%Q3_1b=input('请输入右边向量系数');%Thecoefficientvectoroftherightsideofthedifferentialequationa=input('请输入左边向量系数');%Thecoefficientvectoroftheleftsideofthedifferentialequation[H,w]=freqs(b,a);%ComputethefrequencyresponseHHm=abs(H);%ComputethemagnituderesponseHmphai=angle(H);%ComputethephaseresponsephaiHr=real(H);%ComputetherealpartofthefrequencyresponseHi=imag(H);%Computetheimaginarypartofthefrequencyresponsesubplot(221)5plot(w,Hm),gridon,title('Magnituderesponse'),xlabel('Frequencyinrad/sec')subplot(223)plot(w,phai),gridon,title('Phaseresponse'),xlabel('Frequencyinrad/sec')subplot(222)plot(w,Hr),gridon,title('Realpartoffrequencyresponse'),xlabel('Frequencyinrad/sec')subplot(224)plot(w,Hi),gridon,title('Imaginarypartoffrequencyresponse'),xlabel('Frequencyinrad/sec')执行程序Q3_1，绘制的系统1的频率响应特性曲线如下：从系统1的幅度频率响应曲线看，系统1是低通、高通、全通、带通还是带阻滤波器？051000.10.20.30.4MagnituderesponseFrequencyinrad/sec0510-4-3-2-10PhaseresponseFrequencyinrad/sec0510-0.100.10.20.3RealpartoffrequencyresponseFrequencyinrad/sec0510-0.3-0.2-0.10ImaginarypartoffrequencyresponseFrequencyinrad/sec6答：执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：从系统2的幅度频率响应曲线看，系统2低通、高通、全通、带通还是带阻滤波器？答：执行程序Q3_1，绘制的系统3的频率响应特性曲线如下：05101111MagnituderesponseFrequencyinrad/sec051001234PhaseresponseFrequencyinrad/sec0510-1-0.500.51Realpar