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

实验三 连续时间LTI系统的频域分一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用MATLAB语言进行系统频响特性分析的方法。基本要求：掌握LTI连续和离散时间系统的频域数学模型和频域数学模型的MATLAB描述方法，深刻理解LTIMATLAB计算和绘制LTI系统频率响应特性曲线中的编程。二、实验原理及方法LTI系统的频率响应（Frequencyrespons，是指系统在位随频率的变化情况两个方面。x(x(t)X(j)LTI系统h(t)H(j)y(t)Y(j)连续时间LTI系统的时域及频域分析图上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)y(t)x(t)*h(t，由傅里叶变换的时域卷积定理可得到：Y(j)X(j)H(j)

3.1或者： H(j)

Y(j)X(j) 3.2H(j)为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即H(j)

h(t)edt 3.3PAGEPAGE13由于H(j)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutlyintegrabel）的话，那么H(j)一定存在，而且H(j)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式：H(j)H(j)ej() 3.4上式中，H(称为幅度频率相（Magnituderesponse，反映信号经过系统之后，信号频率分量的幅度发生变化的情况，)称为（Phaseresponse，反映信号经过系统后，信号各频率分量在相位上发生变换的情况。H(j)和()都是频的函数。对于一个系统，其频率响应为H(j)，其幅度响应和相位响应分别为H(j)和()，x(t)e，则其响应信号为y(t)H(j

)ej

H(

)ej(

)ejt H(

)ej(t())3.5

0 0 0 0 00 0 03.6

若输入信号为正弦信号，即x(t)=sin(0t)，则系统响应为y(t)H(j)sin(t)|H(j)|sin(t())0 0 0 0 0可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被H(j)加权，二是信号的相位要被()移相。由于H(j)和()都是频的函数所以系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。LTI系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted所组成。正如刚才所述，信号经过LTI系统传输与处理时，系统将会对信号中的所有频率分量造（Phaseshiftin（Timedela群延时Groupdela映系统对不同频率分量造成的延时。LTI系统的群延时定义为：()d() 3.7d群延时的物理意义：群延时描述的是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说，该系统对于所有的频率分量造成的延时时间都是一样的Phasedistortio。反之，若系统的相位频率响应特性不是线性的，则该系统对于不同频率的频率分量造成的延时时间是不同的，因此，当信号经过系统后，必将产生相位失真。MATLAB计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。实验中用到的MATLAB函数如下：[H,w]=freqs(b,a)：b,a分别为连续时间LTI系统的微分方程右边的和左边的系数向量（Coefficientsvecto，返回的频率响应在各频率点的样点值（复数）存放在H样点数目200点；Hm=abs(H)：求模数，即进行HmH运算，求得系统的幅度频率响应，返回值存于Hm之中。real(H)：求H的实部；imag(H)：求Hphi=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系统，下面的微分方程描述其输入输出之间的关系d2y(t)dt2

3dy(t)dt

2y(t)x(t)编写的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')subplot(222)plot(w,Hr),gridon, title('Realpartoffrequencyresponse'),xlabel('Frequencyinrad/sec')subplot(224)plot(w,Hi),gridon, title('Imaginarypartoffrequencyresponse'),xlabel('Frequencyinrad/sec')三、实验内容及步骤MATLAB例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。给定三个连续时间LTI系统，它们的微分方程分别为系统1：

d2y(t)

1dy(t)

25y(t)

dx(t)

Eq.3.1dt2 dt dt系统2：

dy(t)

y(t)

dx(t)

x(t) Eq.3.2dt dt系统3：d6y(t)

10d5y(t)

48d4y(t)

148d3y(t)

306d2y(t)401dy(t)

262y(t)262x(t)dt6

dt5

dt4

dt3

dt2

Eq.3.3Q3-1修改程序Program3_1，并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程Q3_1如下：% Q3_1b=请输入右边向量系'); %Thecoefficientvectoroftherightsideofthedifferentialequationa=请输入左边向量系数'); %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')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的频率响应特性曲线如下：Magnituderesponse Realpartoffrequency0.4 0.30.3 0.20.2 0.10.1 0-0.10 5 10 0 5 10-0.10 5 10 0 5 10Frequencyinrad/secFrequencyinrad/secPhaseresponse Imaginarypartoffrequency0 0-1-0.1-2-0.2-34 -0.30 5 10 0 5 104 -0.30 5 10 0 5 10Frequencyinrad/secFrequencyinrad/sec从11、、、答：执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：Magnituderesponse Realpartoffrequency1 10.5101-0.5-10 5 10 0 5 10-10 5 10 0 5 10Frequencyinrad/secFrequencyinrad/secPhaseresponse Imaginarypartoffrequency4 132 0.5100 5 10 0 5 1000 5 10 0 5 10Frequencyinrad/secFrequencyinrad/sec从22、、、答：执行程序Q3_1，绘制的系统3的频率响应特性曲线如下：Magnituderesponse1 0.50.5

Realpartoffrequencyresponse00 5Frequencyinrad/sec

-0.5-110

5 10Frequencyinrad/secPhaseresponseImaginarypartoffrequencyresponse4120.500-2-0.5-40 5 Frequencyinrad/sec

-10 5 10Frequencyinrad/sec从33是、、、答：这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系？请根据傅里叶变换的性质说明为什么会具有这些对称关系？答：Q3-2编写程序Q3_2，使之能够能够接受键盘方式输入的输入信号x(t)的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号y(t)的幅度频谱，系统的单位冲激响应h(t)，并按照下面的图Q3-2的布局，绘制出各个信号的时域和频域图形。你编写的程序Q3_2抄写如下：%Q3_2

图Q3-2b=input('Typeintherightcoefficientvectorofdifferentialequation ：'); %Thecoefficientvectoroftherightsideofthedifferentialequationa=input('Typeintheleftcoefficientvectorofdifferentialequation：'); %Thecoefficientoftheleftsideofthedifferentialequationw=-10:0.01:10;H=freqs(b,a,w); %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');subplot(222);plot(w,Hr);gridon, title('Realpartoffrequencyresponse');xlabel('Frequencyinrad/sec');subplot(224);plot(w,Hi);gridon, title('Imaginarypartoffrequencyresponse');xlabel('Frequencyinrad/sec');执行程序Q3_2，输入信号x(t)=sin(t)+sin(8t)，输入由Eq.3.3描述的系统。得到的图形如下：此处粘帖执行程序Q3_2所得到的图形Magnituderesponse15105

21.510.5

Realpartoffrequencyresponse0-10 -5

5

0-10 -5

5 10Frequencyinrad/secPhaseresponseFrequencyinrad/secImaginarypartoffrequencyresponse2101500-1-5-2-10 -5

5

-10-10 -5

5 10Frequencyinrad/sec Frequencyinrad/sec请手工绘制出信号x(t)=sin(t)+sin(8t)的幅度频谱图如下：你手工绘制的信号x(t)=sin(t)+sin(8t)的幅度频谱图与执行程序Q3_2得到的x(t)=sin(t)+sin(8t)的幅度频谱图是否相同？如不同，是何原因造成的？答：执行程序Q3_2得到的x(t)=sin(t)+sin(8t)的幅度频谱图实际上是另外一个信号x1(t)的幅度谱，这个信号的时域数学表达式为 x1(t)=请利用傅里叶变换的相关性质计算并绘制信号x1(t)的幅度频谱图。计算过程：手工绘制的x1(t)的幅度频谱图如下：结合所学的有关滤波的知识，根据上面所得到的信号的时域和频域图形，请从时域和频域两个方面解释滤波的概念。答：Q3-3编写程序Q3_3系统的群延时曲线图。抄写程序Q3_3如下：%Q3_3b1=input('Typeintherightcoefficientvectorofdifferentialequation：'); %Thecoefficientvectoroftherightsideofthedifferentialequationa1=input('Typeintheleftcoefficientvectorofdifferential%Thecoefficientoftheleftsideofthedifferentialequationb2=input('Typeintherightcoefficientvectorofdifferentialequation：'); %Thecoefficientvectoroftherightsideofthedifferentialequationa2=input('Typeintheleftcoefficientvectorofdifferential%Thecoefficientoftheleftsideofthedifferentialequationb3=input('Typeintherightcoefficientvectorofdifferentialequation：'); %Thecoefficientvectoroftherightsideofthedifferentialequationa3=input('Typeintheleftcoefficientvectorofdifferential%Thecoefficientoftheleftsideofthedifferentialequationw=-10:0.01:10;H1=freqs(b1,a1,w); %Computethefrequencyresponsephi1=angle(H1);H2=freqs(b2,a2,w); %Computethefrequencyresponsephi2=angle(H2);H3=freqs(b3,a3,w); %Computethefrequencyresponsephi3=angle(H3);tao1=grpdelay(b1,a1,w);tao2=grpdelay(b2,a2,w);tao3=grpdelay(b3,a3,w);subplot(321);plot(w,phi1);gridon, title('Phaseresponseofsubplot(323);plot(w,phi2);gridon, title('Phaseresponseofs