实验三 连续时间LTI系统的频域分析汇总文档格式.docx

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

2、换的情况。和都是频率 的函数。对于一个系统，其频率响应为H（j ），其幅度响应和相位响应分别为和，如果作用于系统的信号为，则其响应信号为 3.5若输入信号为正弦信号，即x（t） = sin（ 0t），则系统响应为 3.6可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被加权，二是信号的相位要被移相。由于和都是频率 的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。2 LTI系统的群延时从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted sum）所组成。正如刚才所述，信号经过LTI系统传输与处理时，系统将会对信号中的所有频率分量造

3、成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phase shifting），相位移实际上就是延时（Time delay）。群延时（Group delay）的概念能够较好地反映系统对不同频率分量造成的延时。LTI系统的群延时定义为： 3.7群延时的物理意义：群延时描述的是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。如果系统的相位频率响应特性是线性的，则群延时为常数，也就是说，该系统对于所有的频率分量造成的延时时间都是一样的，因而，系统不会对信号产生相位失真（Phase distortion）。反之，若系统的相位频率响应特性

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

5、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。 下面举例说明如何利

6、用上述函数计算并绘制系统频率响应特性曲线的编程方法。假设给定一个连续时间LTI系统，下面的微分方程描述其输入输出之间的关系编写的MATLAB范例程序，绘制系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部。程序如下：% 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 differen

7、tial 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 = angle（H）; % Compute the phase response phaiHr = real（H）; % Compute the real part of the frequency r

8、esponseHi = 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（Phase responsesubplot（222）plot（w,Hr）, grid on, title（Real part of frequency response）, xlabel（su

9、bplot（224）plot（w,Hi）, grid on, title（Imaginary part of frequency response三、实验内容及步骤实验前，必须首先阅读本实验原理，了解所给的MATLAB相关函数，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。给定三个连续时间LTI系统，它们的微分方程分别为系统1： Eq.3.1系统2： Eq.3.2系统3： E

10、q.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（微分方程右边的系数：plot（w,Hm）, grid on, title（）,plot（w,phai）,plot（w,Hr）,plot（w,Hi）, 执行程序Q3_1，绘制的系统1的频率响应特性曲线如下：从系统1的幅

11、度频率响应曲线看，系统1是低通、高通、全通、带通还是带阻滤波器？答：带通滤波器。执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：从系统2的幅度频率响应曲线看，系统2低通、高通、全通、带通还是带阻滤波器？低通滤波器。执行程序Q3_1，绘制的系统3的频率响应特性曲线如下：从系统3的幅度频率响应曲线看，系统3是低通、高通、全通、带通还是带阻滤波器？带阻滤波器。这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系？请根据傅里叶变换的性质说明为什么会具有这些对称关系？Q3-2 编写程序Q3_2，使之能够能够接受键盘方式输入的输入信号x（t）的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号y（t）的幅度频谱，系统的单位冲激响应h（t），并按照下面的图Q3-2的布局，绘制出各个信号的时域和频域图形。图Q3-2你编写的程序Q3_2抄写如下：t = 0:0.01:40;T =0.01;dw =0.1;w=-4*pi:dw:4*pi;x = input（表达式的输入信号x（t）：subplot（323）;impulse（b,a,40）;axis（0 40 -0.2 1）;grid on,系统单位冲击响应h（t）subplot（321）, plot（t,x）输入信号x（t）:t/s