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

1、.实验三连续时间 LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用 MATLAB 语言进行系统频响特性分析的方法。基本要求 ：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB描述方法 ，深刻理 LTI系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用 MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。二、实验原理及方法1 连续时间 LTI 系统的频率

2、响应所谓 频率特性 ，也称为 频率响应特性 ，简称 频率响应 （Frequency response），是指系统在 正弦信号 激励下的稳态响应随频率变化的情况 ，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。x(t) 、 y(t) 分别为系统的时域激励信号和响应信号， h(t) 是系统的单位冲激响应，它们三者之间的关系为： y(t )x(t ) * h(t ) ，由傅里叶变换的时域卷积定理可得到：Y( j )X ( j ) H ( j )3.1或者 ：H ( jY( j)3.2)X ( jH ( j) 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t) 的傅里叶

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

4、 称为 相位特性 （ Phase response ），反映信号经过系统后 ，信号各频率分量在相位上发生变换的情况。 H ( j ) 和( ) 都是频率的函数 。对于一个系统 ，其频率响应为H( j )， 其幅度响应和相位响应分别为( ) ，如果作用于系统的信号为 x(t )e j 0 t，则其响应信号为H ( j) 和y(t)H ( j0 )e j 0tH ( j0 ) e j ( 0 ) e j 0tH ( j0 ) e j ( 0t(0 )3.5若输入信号为正弦信号 ，即 x(t) = sin(0 t) ，则系统响应为y(t ) H ( j 0 ) sin( 0t )| H ( j 0

5、) | sin( 0t( 0 )3.6可见 ，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被H ( j) 加权，二是信号的相位要被() 移相 。由于 H ( j) 和() 都是频率的函数 ，所以 ，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。.专业 .专注.2 LTI 系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和 （ Weightedsum ）所组成 。 正如刚才所述 ，信号经过 LTI 系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移 （ Phaseshifti

6、ng）， 相位移实际上就是延时 （ Timedelay ） 。 群延时 （ Groupdelay ）的概念能够较好地反映系统对不同频率分量造成的延时。LTI 系统的群延时定义为：()3.7d()d群延时的物理意义：群延时描述的是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说 ，该系统对于所有的频率分量造成的延时时间都是一样的，因而 ，系统不会对信号产生相位失真 （ Phasedistortion）。 反之 ， 若系统的相位频率响应特性不是线性的，则该系统对于不同频率的频率分量造成的延时时间是不同

7、的，因此，当信号经过系统后，必将产生相位失真。3 用 MATLAB 计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。实验中用到的MATLAB 函数如下 ：H,w=freqs(b,a) ： b,a 分别为连续时间LTI 系统的微分方程右边的和左边的系数向量 （ Coefficientsvector ） ，返回的频率响应在各频率点的样点值（复数 ）存放在H中，系统默认的 样点数目 为 200 点；Hm = abs(H)：求模数 ，即进行 HmH 运算，求得系统的幅度频率响应，返回.专业 .专注.值存于 Hm 之中 。real(H) ：求 H 的实部 ；i

8、mag(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*p

9、i/256:2*pi。下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。假设给定一个连续时间LTI 系统，下面的微分方程描述其输入输出之间的关系d2 y(t )dy(t)()dt 223 dtx t编写的 MATLAB 范例程序 ，绘制系统的幅度响应特性、相位响应特性 、频率响应的实部和频率响应的虚部。程序如下 ：% Program3_1% This Program is used to compute and draw the plots of the frequencyresponse% of a continuous-time systemb =1;% The coe

10、fficientvectorof therightsideofthe differentialequation.专业 .专注.a = 1 3 2;% The coefficientvectoroftheleftside ofthedifferentialequationH,w = freqs(b,a);% Compute the frequency response HHm = abs(H);% Compute the magnitude response Hmphai = angle(H);% Compute the phase response phaiHr = real(H);% Com

11、pute the real part of the frequency responseHi = imag(H);% Compute the imaginary part of the frequencyresponsesubplot(221)plot(w,Hm),gridon,title('Magnituderesponse'),xlabel('Frequencyinrad/sec')subplot(223)plot(w,phai), grid on,title('Phase response'),xlabel('Frequency i

12、n rad/sec')subplot(222)plot(w,Hr), grid on,title('Real part of frequency response'),xlabel('Frequency in rad/sec')subplot(224)plot(w,Hi), grid on,title('Imaginary part of frequency response'),xlabel('Frequency in rad/sec')三、实验内容及步骤实验前 ，必须首先阅读本实验原理，了解所给的MATLAB 相关函数

13、 ，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序， 观察所得到的信号的波形.专业 .专注.图。 并结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序 。实验前 ，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。给定三个连续时间LTI 系统 ，它们的微分方程分别为系统 1：d 2 y(t )1dy (t )25 y(t )dx (t)Eq.3.1dt 2dtdt系统 2：dy(t)y(t)dx(t )x(t )Eq.3.2dtdt系统 3：d 6 y(t )10 d 5 y(t )48 d 4 y(t)148

14、d 3 y(t)306 d 2 y(t)401 dy(t)262 y(t ) 262 x(t)dt 6dt 5dt 4dt 3dt 2dtEq.3.3Q3-1修改程序Program3_1 ，并以 Q3_1 存盘，使之能够能够接受键盘方式输入的微分方程系数向量 。 并利用该程序计算并绘制由微分方程Eq.3.1 、 Eq.3.2 和 Eq.3.3 描述的系统的幅度响应特性 、相位响应特性、频率响应的实部和频率响应的虚部曲线图 。抄写程序 Q3_1 如下：clear,close all;a = input('微分方程左边的系数： ');b = input('微分方程右边的系数

15、： ');H,w = freqs(b,a);Hm = abs(H);.专业 .专注.phai = angle(H);Hr = real(H);Hi = imag(H);subplot(221)plot(w,Hm),grid on,title('Magnitude response'),xlabel('Frequency in rad/sec')subplot(223)plot(w,phai),grid on,title('Phase response'),xlabel('Frequency in rad/sec')subp

16、lot(222)plot(w,Hr),grid on,title('Real part of frequency response'),xlabel('Frequency in rad/sec')subplot(224)plot(w,Hi),grid on,title('Imaginary part of frequency response'),.专业 .专注.xlabel('Frequency in rad/sec')执行程序 Q3_1 ，绘制的系统1 的频率响应特性曲线如下：Magnitude responseReal pa

17、rt of frequency response110.50.50510051000Frequency in rad/secFrequency in rad/secPhase responseImaginary part of frequency response20.5100-1-2510-0.551000Frequency in rad/secFrequency in rad/sec从系统 1 的幅度频率响应曲线看 ，系统 1是低通、高通、全通、带通还是带阻滤波器 ？答：带通滤波器 。.专业 .专注.执行程序Q3_1 ，绘制的系统 2的频率响应特性曲线如下：Magnitude respon

18、seReal part of frequency response1110.501-0.51510-151000Frequency in rad/secFrequency in rad/secPhase responseImaginary part of frequency response41320.510510051000Frequency in rad/secFrequency in rad/sec从系统 2 的幅度频率响应曲线看，系统 2低通、高通、全通、带通还是带阻滤波器 ？答：低通滤波器 。执行程序Q3_1 ，绘制的系统 3的频率响应特性曲线如下：.专业 .专注.Magnitude

19、 responseReal part of frequency response110.50.50-0.50510-151000Frequency in rad/secFrequency in rad/secPhase responseImaginary part of frequency response4120.500-2-0.5-4510-151000Frequency in rad/secFrequency in rad/sec从系统 3 的幅度频率响应曲线看，系统 3 是低通、高通、全通 、带通 还是带阻 滤波器 ？答：带阻滤波器 。这三个系统的幅度频率响应、相位频率相应、频率响应的

20、实部以及频率响应的虚部分别具有何种对称关系？请根据傅里叶变换的性质说明为什么会具有这些对称关系？答：Q3-2 编写程序Q3_2，使之能够能够接受键盘方式输入的输入信号x(t) 的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号 y(t) 的幅度频谱 ，系统的单位冲激响应h(t) ，并按照下面的图Q3-2 的布局 ，绘制出各个信号的时域和频域图形。.专业 .专注.图 Q3-2你编写的程序Q3_2 抄写如下 ：clear,close all;t = 0:0.01:40;T =0.01; dw =0.1;w=-4*pi:dw:4*pi;a = input(

21、'微分方程左边的系数： ');b = input('微分方程右边的系数： ');x = input('表达式的输入信号x（ t）： ');subplot(323);impulse(b,a,40);.专业 .专注.axis(0 40 -0.2 1);grid on,title(' 系统单位冲击响应h(t)')subplot(321),plot(t,x)title(' 输入信号x(t):');xlabel('t/s');subplot(325),y=lsim(b,a,x,t);plot(t,y)titl

22、e(' 输出信号y(t)');xlabel('t/s');X=x*exp(-j*t'*w)*T;X1=abs(X);subplot(322);plot(w,X1),axis(-4*pi 4*pi 0 20);grid on,title(' 输入信号x(t)的幅度频谱 ').专业 .专注.xlabel(' 频率弧度 / 秒 ')Y=y'*exp(-j*t'*w)*T;Y1=abs(Y);subplot(326)plot(w,Y1),axis(-4*pi 4*pi 0 20);grid on,title(

23、9; 输出信号y(t) 的幅度频谱 ');xlabel(' 频率弧度 / 秒 ')H,w = freqs(b,a);Hm = abs(H);phai = angle(H);subplot(324)plot(w,Hm),grid on,title(' 系统的幅度频率响应:'),xlabel(' 频率弧度 / 秒 ').专业 .专注.执行程序 Q3_2 ，输入信号x(t) = sin(t) + sin(8t)，输入由 描述的系统 。 得到的图形如下 ：此处粘帖执行程序Q3_2 所得到的图形输 入 信 号 x(t):输 入 信 号 x(t) 的 幅 度 频 谱220115010-15-2102030400-505100-10t/s频 率弧度/秒系统的幅度频率响应:系统单位冲击响应 h(t)11ed0.50.5utilpmA005100102030400Time (sec)频 率弧度/秒输 出 信 号 y(t