发布时间 : 星期三 文章信号与系统实验报告实验三 连续时间LTI系统的频域分析更新完毕开始阅读838c03dc00d276a20029bd64783e0912a3167c4d
信号与系统实验报告实验三 连续时间LTI系统的频域分析
实验三 连续时间LTI系统的频域分析
一、实验目的
1、掌握系统频率响应特性的概念及其物理意义;
2、掌握系统频率响应特性的计算方法与特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;
3、学习与掌握幅度特性、相位特性以及群延时的物理意义; 4、掌握用MATLAB语言进行系统频响特性分析的方法。
基本要求:掌握LTI连续与离散时间系统的频域数学模型与频域数学模型的MATLAB描述方法,深刻理解LTI系统的频率响应特性的物理意义,理解滤波与滤波器的概念,掌握利用MATLAB计算与绘制LTI系统频率响应特性曲线中的编程。
二、实验原理及方法
1 连续时间LTI系统的频率响应
所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response),就是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况与响应的相位随频率的变化情况两个方面。
x(t)X(j?)LTI系统h(t)H(j?)y(t)Y(j?)连续时间LTI系统的时域及频域分析图 上图中x(t)、y(t)分别为系统的时域激励信号与响应信号,h(t)就是系统的单位冲激响应,它们三者之间的关系为:y(t)?x(t)*h(t),由傅里叶变换的时域卷积定理可得到:
Y(j?)?X(j?)H(j?)
或者: H(j?)?3、1
Y(j?) 3、2
X(j?)H(j?)为系统的频域数学模型,它实际上就就是系统的单位冲激响应h(t)的傅里叶变换。即
? H(j?)????j?th(t)edt 3、3 ?信号与系统实验报告实验三 连续时间LTI系统的频域分析
由于H(j?)实际上就是系统单位冲激响应h(t)的傅里叶变换,如果h(t)就是收敛的,或者说就是绝对可积(Absolutly integrabel)的话,那么H(j?)一定存在,而且H(j?)通常就是复数,因此,也可以表示成复数的不同表达形式。在研究系统的频率响应时,更多的就是把它表示成极坐标形式:
H(j?)?H(j?)ej?(?) 3、4
上式中,H(j?)称为幅度频率相应(Magnitude response),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,?(?)称为相位特性(Phase response),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。H(j?)与?(?)都就是频率?的函数。
对于一个系统,其频率响应为H(j?),其幅度响应与相位响应分别为H(j?)与?(?),如果作用于系统的信号为x(t)?ej?0t,则其响应信号为
y(t)?H(j?0)ej?0t ?H(j?0)ej?(?0)ej?0t?H(j?0)ej(?0t??(?0)) 3、
5
若输入信号为正弦信号,即x(t) = sin(?0t),则系统响应为
y(t)?H(j?0)sin(?0t)?|H(j?0)|sin(?0t??(?0)) 3、6
可见,系统对某一频率分量的影响表现为两个方面,一就是信号的幅度要被H(j?)加权,二就是信号的相位要被?(?)移相。
由于H(j?)与?(?)都就是频率?的函数,所以,系统对不同频率的频率分量造成的幅度与相位上的影响就是不同的。
2 LTI系统的群延时
从信号频谱的观点瞧,信号就是由无穷多个不同频率的正弦信号的加权与(Weighted sum)所组成。正如刚才所述,信号经过LTI系统传输与处理时,系统将会对信号中的所有频率分量造成幅度与相位上的不同影响。从相位上来瞧,系统对各个频率分量造成一定的相位移(Phase shifting),相位移实际上就就是延时(Time delay)。群延时(Group delay)的概念能够较好地反映系统对不同频率分量造成的延时。
信号与系统实验报告实验三 连续时间LTI系统的频域分析
LTI系统的群延时定义为:
?(?)??7
d?(?) 3、d?群延时的物理意义:群延时描述的就是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。
如果系统的相位频率响应特性就是线性的,则群延时为常数,也就就是说,该系统对于所有的频率分量造成的延时时间都就是一样的,因而,系统不会对信号产生相位失真(Phase distortion)。反之,若系统的相位频率响应特性不就是线性的,则该系统对于不同频率的频率分量造成的延时时间就是不同的,因此,当信号经过系统后,必将产生相位失真。
3 用MATLAB计算系统频率响应
在本实验中,表示系统的方法仍然就是用系统函数分子与分母多项式系数行向量来表示。实验中用到的MATLAB函数如下:
[H,w] = freqs(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):计算系统的相位频率响应所对应的群延时。 计算频率响应的函数freqs()的另一种形式就是:
H = freqs(b,a,w):在指定的频率范围内计算系统的频率响应特性。在使用这种形式的
freqs/freqz函数时,要在前面先指定频率变量w的范围。
例如在语句H = freqs(b,a,w)之前加上语句:w = 0:2*pi/256:2*pi。
下面举例说明如何利用上述函数计算并绘制系统频率响应特性曲线的编程方法。 假设给定一个连续时间LTI系统,下面的微分方程描述其输入输出之间的关系
d2y(t)dy(t)?3?2y(t)?x(t)
dt2dt编写的MATLAB范例程序,绘制系统的幅度响应特性、相位响应特性、频率响应的实部与
频率响应的虚部。程序如下:
% Program3_1
% This Program is used to compute and draw the plots of the frequency response % of a continuous-time system
b = [1]; % The coefficient vector of the right side of the differential equation
信号与系统实验报告实验三 连续时间LTI系统的频域分析
a = [1 3 2]; % 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')
三、实验内容及步骤
实验前,必须首先阅读本实验原理,了解所给的MATLAB相关函数,读懂所给出的全部范例程序。实验开始时,先在计算机上运行这些范例程序,观察所得到的信号的波形图。并结合范例程序所完成的工作,进一步分析程序中各个语句的作用,从而真正理解这些程序。
实验前,一定要针对下面的实验项目做好相应的实验准备工作,包括事先编写好相应的实验程序等事项。
给定三个连续时间LTI系统,它们的微分方程分别为
d2y(t)dy(t)dx(t)?1?25y(t)?系统1: Eq、3、1
dtdtdt2系统2: 系统3:
dy(t)dx(t)?y(t)??x(t) Eq、3、2 dtdtd6y(t)d5y(t)d4y(t)d3y(t)d2y(t)dy(t)?10?48?148?306?401?262y(t)?262x(t)65432dtdtdtdtdtdt Eq、3、3
Q3-1 修改程序Program3_1,并以Q3_1存盘,使之能够能够接受键盘方式输入的微分方程系数
向量。并利用该程序计算并绘制由微分方程Eq、3、1、Eq、3、2与Eq、3、3描述的系统的幅度响应特性、相位响应特性、频率响应的实部与频率响应的虚部曲线图。