信号与系统实验报告实验三 连续时间LTI系统的频域分析文档格式.docx
- 文档编号:20540179
- 上传时间:2023-01-23
- 格式:DOCX
- 页数:14
- 大小:164.60KB
信号与系统实验报告实验三 连续时间LTI系统的频域分析文档格式.docx
《信号与系统实验报告实验三 连续时间LTI系统的频域分析文档格式.docx》由会员分享,可在线阅读,更多相关《信号与系统实验报告实验三 连续时间LTI系统的频域分析文档格式.docx(14页珍藏版)》请在冰豆网上搜索。
由于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)的数学表达式,系统微分方程的系数向量,计算输入信号的幅度频谱,系统的幅度频率响应,系统输出信号y(t)的幅度频谱,系统的单位冲激响应h(t),并按照下面的图Q3-2的布局,绘制出各个信号的时域与频域图形。
图Q3-2
您编写的程序Q3_2抄写如下:
%Q3_2
Typeintherightcoefficientvectorofdifferentialequation:
'
Typeintheleftcoefficientvectorofdifferentialequation:
w=-10:
0、01:
10;
H=freqs(b,a,w);
subplot(221);
plot(w,Hm);
gridon,title('
subplot(223);
plot(w,phai);
subplot(222);
plot(w,Hr);
subplot(224);
plot(w,Hi);
执行程序Q3_2,输入信号x(t)=sin(t)+sin(8t),输入由Eq、3、3描述的系统。
得到的图形如下:
此处粘帖执行程序Q3_2所得到的图形
请手工绘制出信号x(t)=sin(t)+sin(8t)的幅度频谱图如下:
您手工绘制的信号x(t)=sin(t)+sin(8t)的幅度频谱图与执行程序Q3_2得到的x(t)=sin(t)+sin(8t)的幅度频谱图就是否相同?
如不同,就是何原因造成的?
执行程序Q3_2得到的x(t)=sin(t)+sin(8t)的幅度频谱图实际上就是另外一个信号x1(t)的幅度频谱,这个信号的时域数学表达式为x1(t)=
请利用傅里叶变换的相关性质计算并绘制信号x1(t)的幅度频谱图。
计算过程:
手工绘制的x1(t)的幅度频谱图如下:
结合所学的有关滤波的知识,根据上面所得到的信号的时域与频域图形,请从时域与频域两个方面解释滤波的概念。
Q3-3编写程序Q3_3,能够接受从键盘输入的系统微分方程系数向量,并分别绘制所给三个系统的群延时曲线图。
抄写程序Q3_3如下:
%Q3_3
b1=input('
a1=input('
b2=input('
a2=input('
b3=input('
a3=input('
H1=freqs(b1,a1,w);
phi1=angle(H1);
H2=freqs(b2,a2,w);
phi2=angle(H2);
H3=freqs(b3,a3,w);
phi3=angle(H3);
tao1=grpdelay(b1,a1,w);
tao2=grpdelay(b2,a2,w);
tao3=grpdelay(b3,a3,w);
subplot(321);
plot(w,phi1);
Phaseresponseofnum1'
subplot(323);
plot(w,phi2);
Phaseresponseofnum2'
subplot(325);
plot(w,phi3);
Phaseresponseofnum3'
subplot(322);
plot(w,tao1);
Groupdelayofnum1'
subplot(324);
plot(w,tao2);
Groupdelayofnum2'
subplot(326);
plot(w,tao3);
Groupdelayofnum3'
系统Eq、3、1的群延时曲线图系统Eq、3、3的群延时曲线图
根据上面的群延时曲线图可以瞧出,对系统Eq、3、1,当频率为5弧度/秒时,群延时为秒,当频率为10弧度/秒时,群延时为秒,如何解释这两个群延时时间?
根据上面的群延时曲线图,说明这两个系统就是否会造成对信号的相位失真?
为什么?
从系统Eq、3、3的群延时曲线图中可以瞧出,当信号的频率为1弧度/秒时,系统Eq、3、3对这一频率的信号的延时就是秒。
所以,执行程序Q3_2时,当作用于系统Eq、3、3的输入信号为x(t)=sin(t)+sin(8t)时,其输出信号y(t)的数学表达式为:
四、实验报告要求
1、按要求完整书写您所编写的全部MATLAB程序
2、详细记录实验过程中的有关信号波形图(存于自带的U盘中),图形要有明确的标题。
全部的MATLAB图形应该用打印机打印,然后贴在本实验报告中的相应位置,禁止复印件。
3、实事求就是地回答相关问题,严禁抄袭。
本实验完成时间:
年月日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号与系统实验报告实验三 连续时间LTI系统的频域分析 信号 系统 实验 报告 连续 时间 LTI 分析