连续时间LTI系统的频域分析

连续时间LTI 系统的频域分析

一、实验目的

1、掌握系统频率响应特性的概念及其物理意义；

2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；

3、学习和掌握幅度特性、相位特性以及群延时的物理意义；

4、掌握用MATLAB 语言进行系统频响特性分析的方法。

基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。

二、实验原理及方法

1 连续时间LTI 系统的频率响应

所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。

上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响。它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到：

)()()(ωωωj H j X j Y =

3.1

或者： )

()

()(ωωωj X j Y j H =

3.2

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

即

⎰

∞

∞

--=

dt e t h j H t

j ωω)()( 3.3 由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数

，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把

它表示成极坐标形式：

)()()(ωϕωωj e j H j H = 3.4

上式中，)j (ωH 称为幅度频率相应（Magnitude response ），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(ωϕ称为相位特性（Phase response ），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)(ωj H 和)(ωϕ都是频率ω的函数。

对于一个系统，其频率响应为H(j ω)，其幅度响应和相位响应分别为)(ωj H 和

)(ωϕ，如果作用于系统的信号为t j e t x 0

)(ω=，则其响应信号为

t

j e j H t y 0)()(0ωω=

t j j e e j H 00)(0)(ωωϕω=))((000)(ωϕωω+=t j e j H

3.5

若输入信号为正弦信号，即x(t) = sin(ω0t)，则系统响应为

))(sin(|)(|)sin()()(00000ωϕωωωω+==t j H t j H t y 3.6

可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)(ωj H 加权，二是信号的相位要被)(ωϕ移相。

由于)(ωj H 和)(ωϕ都是频率ω的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。

三、实验内容及步骤

Q3-1 修改程序Program3_1，并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。 抄写程序Q3_1如下： %Program Q3_1

%This program is used to compute and draw the the plots of the frequency response %Of a continuous-time system clc,clear,close all;

b = input('Type in the right coefficient vector of differential equation b'); %The coefficient vector of the right side of the differential equation a = input('Type in the left coefficient vector of differential equation a'); %The coefficient vector of the left side of the differential equation

[H,w]= freqs(b,a);%Compute the frequency response H

Hm= abs(H);%Compute the magnitude response Hm

phai= angle(H);%Compute the phase response phai

Hr= real(H);%Compute the real part of the frequency response

Hi= imag(H);%Compute the imaginary part of the frequency response

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,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')

%Eq.3.1 b=[1 0];a=[1 1 25];

%Eq.3.2 b=[1 -1];a=[1 1];

%Eq.3.3 b=[262];a=[1 10 48 148 306 401 262];

执行程序Q3_1，绘制的系统1的频率响应特性曲线如下：

从系统1的幅度频率响应曲线看，系统1是低通、高通、全通、带通还是带阻滤波器？答：带通滤波器。

执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：