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

《实验三 连续时间LTI系统的频域分析报告》由会员分享，可在线阅读，更多相关《实验三 连续时间LTI系统的频域分析报告（11页珍藏版）》请在金锄头文库上搜索。

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)，由傅里叶变换的时域卷积定理可得到：3.1或者：H(jo)二Y (jo)X (jo)3.2H (jo)为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即H (jo ) = J h(t)e- jotdt3.3由于Hjo)实际上是系统单位

3、冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者 说是绝对可积(Absolutly integrabel)的话，那么H(jo) 定存在，而且Hjo)通常是 复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它 表示成极坐标形式：H (jo) = |H (jo) |ej9(o)3.4上式中，|H(jo)|称为幅度频率相应(Magnitude response)，反映信号经过系统之后，信 号各频率分量的幅度发生变化的情况，9(o)称为相位特性(Phase response)，反映信号经过系统后，信号各频率分量在相位上发生变换的情况。|H(j)|和甲()都是频率的

4、函 数。对于一个系统，其频率响应为H(j),其幅度响应和相位响应分别为|H ( jo )和甲()，如果作用于系统的信号为x(t) = ejo0，则其响应信号为二 |H (jo ) |eje (o0) e jo0t二 |H (jo ) |ej(o0t+e (o0)3.5若输入信号为正弦信号，即x(t) = sin(o#)，则系统响应为y(t)二 H(jo )sin(o t) =1 H(jo ) | sin(o t + Q(o )0 0 0 0 03.6可见，系统对某一频率分量的影响表现为两个方面,一是信号的幅度要被|H(j)加权, 二是信号的相位要被9(o)移相。由于|H(jo)|和Q(o)都是

5、频率o的函数，所以，系统对不同频率的频率分量造成的幅 度和相位上的影响是不同的。2 LTI 系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的 加权和(Weighted sum)所组成。正如刚才所述，信号经过LTI系统传输与处理时，系统将会对信号中的所有 频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相 位移(Phase shifting),相位移实际上就是延时(Time delay)。群延时(Group delay)的 概念能够较好地反映系统对不同频率分量造成的延时。LTI系统的群延时定义为:T (o)=旳(o)do3.7群延时的物理意义:群

6、延时描述的是信号中某一频率分量经过线性时不变系统传输处理 后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说，该系统对于所有 的频率分量造成的延时时间都是一样的，因而，系统不会对信号产生 相位失真( Phase distortion)o反之，若系统的相位频率响应特性不是线性的，则该系统对于不同频率的频率 分量造成的延时时间是不同的，因此，当信号经过系统后，必将产生相位失真。3 用 MATLAB 计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。 实验中用到的MATLAB函数如下：H,w = freqs

7、(b,a)： b,a分别为连续时间LTI系统的微分方程右边的和左边的系数向 量(Coefficients vector)，返回的频率响应在各频率点的样点值(复数)存放在H中，系统默认的样点数目为 200点；Hm = abs(H)：求模数，即进行Hm = |H|运算，求得系统的幅度频率响应，返回值存于 Hm 之中。real(H)：求H的实部；imag(H)：求H的虚部；phi = atan(-imag(H)./(real(H)+eps)：求相位频率相应特性，atan()用来计算反正切值； 或者phi = angle(H)：求相位频率相应特性；tao = grpdelay(num,den,w)：计

8、算系统的相位频率响应所对应的群延时。 计算频率响应的函数freqs()的另一种形式是：H = freqs(b,a,w)：在指定的频率范围内计算系统的频率响应特性。在使用这种形式 的freqs/freqz函数时，要在前面先指定频率变量w的范围。例如在语句H = freqs(b,a,w)之前加上语句：w = 0:2*pi/256:2*pi。 下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。 假设给定一个连续时间LTI系统，下面的微分方程描述其输入输出之间的关系如 + 3 皱 + 2 = x(t)dt2dt编写的 MATLAB 范例程序，绘制系统的幅度响应特性、相位响应特性、频

9、率响应的实 部和频率响应的虚部。程序如下：% Program3_1% This Program is used to compute and draw the plots of the frequency response% of a continuous-time systemb = 1;% The coefficient vector of the right side of the differential equationa = 1 3 2;% The coefficient vector of the left side of the differential equationH,w

10、 = freqs(b,a); Hm = abs(H); phai = angle(H); Hr = real(H); Hi = imag(H); subplot(221)% Compute the frequency response H% Compute the magnitude response Hm% Compute the phase response phai% Compute the real part of the frequency response% Compute the imaginary part of the frequency responseplot(w,Hm)

11、, 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) 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(Imagi

12、nary part of frequency response), xlabel(Frequency in rad/sec)三、实验内容及步骤实验前，必须首先阅读本实验原理，了解所给的MATLAB相关函数，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并 结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的 实验程序等事项。给定三个连续时间LTI系统，它们的微分方程分别为系统1系统 2系统 3d2 y (t)dt 2+1 dy(t) dt+ 2

13、5 y (t)=畔 + y (t)=畔-x(t)dtdtdx(t)dtEq.3.1Eq.3.2竺 +10 空 + 48 空 +148 血 + 306 化2 + 401 如 + 262y (t) = 262x(t) dt6dt5dt4dt3dt2dtEq.3.3Q3-1修改程序Program3_l,并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅 度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程序 Q3_1 如下：clear,close all;a = input(微分方

14、程左边的系数：)；b = input(微分方程右边的系数：)；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) subplot(222) plot(w

15、,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_l,绘制的系统1的频率响应特性曲线如下:Real part of frequency responsePhase responseImaginary part of frequency response从系统1的幅度频率响应曲线看,系统1是低通、高通、全通、带通还是带阻滤波器？ 答：带通滤波器。执行程序Q3_l,绘制的系统2的频率响应特性曲线如下:Magnitude responseReal part of frequency responsePhase responseImaginary part of frequency response