《实验三 连续时间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)是系统的单位冲激响应，它们三者之间的关系为： ，由傅里叶变换的时域卷积定理可得到：)(*)(thxty3.1)(jHXjY或者： 3.2)()(jj为系统的频域数学模型，它实际上就是系统的单位冲激响应 h(t)的傅里叶变换。即 )(jH3.3dtethjHj)(由于 H(j)实际上是系统单位冲激响应 h(t)的傅里叶变
3、换，如果 h(t)是收敛的，或者说是绝对可积（Absolutly integrabel）的话，那么 H(j)一定存在，而且 H(j)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式：3.4)()(jeHj上式中， 称为幅度频率相应 （Magnitude response） ，反映信号经过系统之后，信)j(H号各频率分量的幅度发生变化的情况， 称为相位特性（Phase response） ，反映信号)(经过系统后，信号各频率分量在相位上发生变换的情况。 和 都是频率的函)(jH(数。对于一个系统，其频率响应为 H(j)，其幅度响应和相位响应分
4、别为 和)(j，如果作用于系统的信号为 ，则其响应信号为)(tjetx0(tjHty0)() tjjeH0)(0( )(00)(tjeH3.5若输入信号为正弦信号，即 x(t) = sin(0t)，则系统响应为)(sin|(|sin()() 000 tjjty3.6可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被 加)(jH权，二是信号的相位要被 移相。)(由于 和 都是频率 的函数，所以，系统对不同频率的频率分量造成的幅)(jH度和相位上的影响是不同的。2 LTI 系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted sum）所组成。正
5、如刚才所述，信号经过 LTI 系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phase shifting） ，相位移实际上就是延时（Time delay） 。群延时（Group delay）的概念能够较好地反映系统对不同频率分量造成的延时。LTI 系统的群延时定义为：d)()(3.7群延时的物理意义：群延时描述的是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说，该系统对于所有的频率分量造成的延时时间都是一样的，因而，
6、系统不会对信号产生相位失真（Phase distortion） 。反之，若系统的相位频率响应特性不是线性的，则该系统对于不同频率的频率分量造成的延时时间是不同的，因此，当信号经过系统后，必将产生相位失真。3 用 MATLAB 计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。实验中用到的 MATLAB 函数如下：H,w = freqs(b,a)：b,a 分别为连续时间 LTI 系统的微分方程右边的和左边的系数向量（Coefficients vector） ，返回的频率响应在各频率点的样点值（复数）存放在 H 中，系统默认的样点数目为 200 点；Hm
7、= abs(H)：求模数，即进行 运算，求得系统的幅度频率响应，返回Hm值存于 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 函数时，要
8、在前面先指定频率变量 w 的范围。例如在语句 H = freqs(b,a,w)之前加上语句：w = 0:2*pi/256:2*pi。 下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。假设给定一个连续时间 LTI 系统，下面的微分方程描述其输入输出之间的关系)(2)(3)(2txdtyt编写的 MATLAB 范例程序，绘制系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部。程序如下：% Program3_1% This Program is used to compute and draw the plots of the frequency response
9、 % 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 = freqs(b,a); % Compute the frequency response HHm = abs(H); % Compute the magnitude response Hmphai = angl
10、e(H); % Compute the phase response phaiHr = real(H); % Compute the real part of the frequency responseHi = imag(H); % Compute the imaginary part of the frequency responsesubplot(221)plot(w,Hm), grid on, title(Magnitude response), xlabel(Frequency in rad/sec)subplot(223)plot(w,phai), grid on, title(P
11、hase 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(Imaginary part of frequency response), xlabel(Frequency in rad/sec)三、实验内容及步骤实验前，必须首先阅读本实验原理，了解所给的 MATLAB 相关函数，读懂所给出的
12、全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。给定三个连续时间 LTI 系统，它们的微分方程分别为系统 1： Eq.3.1dtxydtty)(25)(1)(2系统 2： Eq.3.2)()(txt系统 3：)(26)()(401)(306)(148)()(10)( 2456 txtydttydtydtytydty Eq.3.3Q3-1 修改程序 Program3_1，并以 Q3_1 存盘
13、，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程 Eq.3.1、Eq.3.2 和 Eq.3.3 描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程序 Q3_1 如下：clear,close all;a = input(微分方程左边的系数：);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(Magnitud
14、e 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(Imaginary part of frequency response
15、), xlabel(Frequency in rad/sec)执行程序 Q3_1，绘制的系统 1 的频率响应特性曲线如下：0 5 1000.51 Magnitude responseFrequency in rad/sec0 5 10-2-1012 Phase responseFrequency in rad/sec0 5 1000.51 Real part of frequency responseFrequency in rad/sec0 5 10-0.500.5Imaginary part of frequency responseFrequency in rad/sec从系统 1 的幅度频率响应曲线看，系统 1 是低通、高通、全通、带通还是带阻滤波器？答：带通滤波器。执行程序 Q3_1，绘制的系统 2 的频率响应特性曲线如下：0 5 101111 Magnitude responseFrequency in rad/sec0 5 1001234 Phase responseFrequ
