信号与系统实验报告实验三 连续时间LTI系统的频域分析.docx

信号与系统实验报告实验三 连续时间LTI系统的频域分析.docx

《信号与系统实验报告实验三 连续时间LTI系统的频域分析.docx》由会员分享，可在线阅读，更多相关《信号与系统实验报告实验三 连续时间LTI系统的频域分析.docx（14页珍藏版）》请在冰豆网上搜索。

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

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

一、实验目的

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

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

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

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

基本要求:

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

二、实验原理及方法

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

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

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

上图中x（t）、y（t）分别为系统的时域激励信号与响应信号,h（t）就是系统的单位冲激响应,它们三者之间的关系为:

由傅里叶变换的时域卷积定理可得到:

3、1

或者:

3、2

为系统的频域数学模型,它实际上就就是系统的单位冲激响应h（t）的傅里叶变换。

即

3、3

由于H（jω）实际上就是系统单位冲激响应h（t）的傅里叶变换,如果h（t）就是收敛的,或者说就是绝对可积（Absolutlyintegrabel）的话,那么H（jω）一定存在,而且H（jω）通常就是复数,因此,也可以表示成复数的不同表达形式。

在研究系统的频率响应时,更多的就是把它表示成极坐标形式:

3、4

上式中,

称为幅度频率相应（Magnituderesponse）,反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,

称为相位特性（Phaseresponse）,反映信号经过系统后,信号各频率分量在相位上发生变换的情况。

与

都就是频率ω的函数。

对于一个系统,其频率响应为H（jω）,其幅度响应与相位响应分别为

与

如果作用于系统的信号为

则其响应信号为

3、5

若输入信号为正弦信号,即x（t）=sin（ω0t）,则系统响应为

3、6

可见,系统对某一频率分量的影响表现为两个方面,一就是信号的幅度要被

加权,二就是信号的相位要被

移相。

由于

与

都就是频率ω的函数,所以,系统对不同频率的频率分量造成的幅度与相位上的影响就是不同的。

2LTI系统的群延时

从信号频谱的观点瞧,信号就是由无穷多个不同频率的正弦信号的加权与（Weightedsum）所组成。

正如刚才所述,信号经过LTI系统传输与处理时,系统将会对信号中的所有频率分量造成幅度与相位上的不同影响。

从相位上来瞧,系统对各个频率分量造成一定的相位移（Phaseshifting）,相位移实际上就就是延时（Timedelay）。

群延时（Groupdelay）的概念能够较好地反映系统对不同频率分量造成的延时。

LTI系统的群延时定义为:

3、7

群延时的物理意义:

群延时描述的就是信号中某一频率分量经过线性时不变系统传输处理后产生的响应信号在时间上造成的延时的时间。

如果系统的相位频率响应特性就是线性的,则群延时为常数,也就就是说,该系统对于所有的频率分量造成的延时时间都就是一样的,因而,系统不会对信号产生相位失真（Phasedistortion）。

反之,若系统的相位频率响应特性不就是线性的,则该系统对于不同频率的频率分量造成的延时时间就是不同的,因此,当信号经过系统后,必将产生相位失真。

3用MATLAB计算系统频率响应

在本实验中,表示系统的方法仍然就是用系统函数分子与分母多项式系数行向量来表示。

实验中用到的MATLAB函数如下:

[H,w]=freqs（b,a）:

b,a分别为连续时间LTI系统的微分方程右边的与左边的系数向量（Coefficientsvector）,返回的频率响应在各频率点的样点值（复数）存放在H中,系统默认的样点数目为200点;

Hm=abs（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系统,下面的微分方程描述其输入输出之间的关系

编写的MATLAB范例程序,绘制系统的幅度响应特性、相位响应特性、频率响应的实部与频率响应的虚部。

程序如下:

%Program3_1

%ThisProgramisusedtocomputeanddrawtheplotsofthefrequencyresponse

%ofacontinuous-timesystem

b=[1];%Thecoefficientvectoroftherightsideofthedifferentialequation

a=[132];%Thecoefficientvectoroftheleftsideofthedifferentialequation

[H,w]=freqs（b,a）;%ComputethefrequencyresponseH

Hm=abs（H）;%ComputethemagnituderesponseHm

phai=angle（H）;%Computethephaseresponsephai

Hr=real（H）;%Computetherealpartofthefrequencyresponse

Hi=imag（H）;%Computetheimaginarypartofthefrequencyresponse

subplot（221）

plot（w,Hm）,gridon,title（'Magnituderesponse'）,xlabel（'Frequencyinrad/sec'）

subplot（223）

plot（w,phai）,gridon,title（'Phaseresponse'）,xlabel（'Frequencyinrad/sec'）

subplot（222）

plot（w,Hr）,gridon,title（'Realpartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

subplot（224）

plot（w,Hi）,gridon,title（'Imaginarypartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

三、实验内容及步骤

实验前,必须首先阅读本实验原理,了解所给的MATLAB相关函数,读懂所给出的全部范例程序。

实验开始时,先在计算机上运行这些范例程序,观察所得到的信号的波形图。

并结合范例程序所完成的工作,进一步分析程序中各个语句的作用,从而真正理解这些程序。

实验前,一定要针对下面的实验项目做好相应的实验准备工作,包括事先编写好相应的实验程序等事项。

给定三个连续时间LTI系统,它们的微分方程分别为

系统1:

Eq、3、1

系统2:

Eq、3、2

系统3:

Eq、3、3

Q3-1修改程序Program3_1,并以Q3_1存盘,使之能够能够接受键盘方式输入的微分方程系数向量。

并利用该程序计算并绘制由微分方程Eq、3、1、Eq、3、2与Eq、3、3描述的系统的幅度响应特性、相位响应特性、频率响应的实部与频率响应的虚部曲线图。

抄写程序Q3_1如下:

%Q3_1

b=input（'请输入右边向量系数'）;%Thecoefficientvectoroftherightsideofthedifferentialequation

a=input（'请输入左边向量系数'）;%Thecoefficientvectoroftheleftsideofthedifferentialequation

[H,w]=freqs（b,a）;%ComputethefrequencyresponseH

Hm=abs（H）;%ComputethemagnituderesponseHm

phai=angle（H）;%Computethephaseresponsephai

Hr=real（H）;%Computetherealpartofthefrequencyresponse

Hi=imag（H）;%Computetheimaginarypartofthefrequencyresponse

subplot（221）

plot（w,Hm）,gridon,title（'Magnituderesponse'）,xlabel（'Frequencyinrad/sec'）

subplot（223）

plot（w,phai）,gridon,title（'Phaseresponse'）,xlabel（'Frequencyinrad/sec'）

subplot（222）

plot（w,Hr）,gridon,title（'Realpartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

subplot（224）

plot（w,Hi）,gridon,title（'Imaginarypartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

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

从系统1的幅度频率响应曲线瞧,系统1就是低通、高通、全通、带通还就是带阻滤波器？

答:

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

从系统2的幅度频率响应曲线瞧,系统2低通、高通、全通、带通还就是带阻滤波器？

答:

执行程序Q3_1,绘制的系统3的频率响应特性曲线如下:

从系统3的幅度频率响应曲线瞧,系统3就是低通、高通、全通、带通还就是带阻滤波器？

答:

这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系？

请根据傅里叶变换的性质说明为什么会具有这些对称关系？

答:

Q3-2编写程序Q3_2,使之能够能够接受键盘方式输入的输入信号x（t）的数学