信号与系统实验报告实验三 连续时间LTI系统的频域分析报告Word文件下载.docx
- 文档编号:16015218
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:14
- 大小:164.56KB
信号与系统实验报告实验三 连续时间LTI系统的频域分析报告Word文件下载.docx
《信号与系统实验报告实验三 连续时间LTI系统的频域分析报告Word文件下载.docx》由会员分享,可在线阅读,更多相关《信号与系统实验报告实验三 连续时间LTI系统的频域分析报告Word文件下载.docx(14页珍藏版)》请在冰豆网上搜索。
,由傅里叶变换的时域卷积定理可得到:
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'
subplot(222)
plot(w,Hr),gridon,title('
Realpartoffrequencyresponse'
),
xlabel('
subplot(224)
plot(w,Hi),gridon,title('
Imaginarypartoffrequencyresponse'
三、实验容及步骤
实验前,必须首先阅读本实验原理,了解所给的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('
请输入右边向量系数'
);
a=input('
请输入左边向量系数'
执行程序Q3_1,绘制的系统1的频率响应特性曲线如下:
从系统1的幅度频率响应曲线看,系统1是低通、高通、全通、带通还是带阻滤波器?
答:
执行程序Q3_1,绘制的系统2的频率响应特性曲线如下:
从系统2的幅度频率响应曲线看,系统2低通、高通、全通、带通还是带阻滤波器?
执行程序Q3_1,绘制的系统3的频率响应特性曲线如下:
从系统3的幅度频率响应曲线看,系统3是低通、高通、全通、带通还是带阻滤波器?
这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系?
请根据傅里叶变换的性质说明为什么会具有这些对称关系?
Q3-2编写程序Q3_2,使之能够能够接受键盘方式输入的输入信号x(t)的数学表达式,系统微分方程的系数向量,计算输入信号的幅度频谱,系统的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号与系统实验报告实验三 连续时间LTI系统的频域分析报告 信号 系统 实验 报告 连续 时间 LTI 分析