实验三 连续时间LTI系统的频域分析Word文档格式.docx
- 文档编号:19258399
- 上传时间:2023-01-04
- 格式:DOCX
- 页数:14
- 大小:158.71KB
实验三 连续时间LTI系统的频域分析Word文档格式.docx
《实验三 连续时间LTI系统的频域分析Word文档格式.docx》由会员分享,可在线阅读,更多相关《实验三 连续时间LTI系统的频域分析Word文档格式.docx(14页珍藏版)》请在冰豆网上搜索。
在研究系统的频率响应时,更多的是把它表示成极坐标形式:
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如下:
clear,
closeall;
a=input('
微分方程左边的系数:
'
);
b=input('
微分方程右边的系数:
plot(w,Hm),
gridon,
title('
),
plot(w,phai),
plot(w,Hr),
plot(w,Hi),
执行程序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抄写如下:
t=0:
0.01:
40;
T=0.01;
dw=0.1;
w=-4*pi:
dw:
4*pi;
x=input('
表达式的输入信号x(t):
subplot(323);
impulse(b,a,40);
axis([040-0.21]);
gridon,
系统单位冲击响应h(t)'
subplot(321),
plot(t,x)
输入信号x(t):
t/s'
subplot(325),
y=lsim(b,a,x,t);
plot(t,y)
输出信号y(t)'
X=x*exp(-j*t'
*w)*T;
X1=abs(X);
subplot(322);
plot(w,X1),
axis([-4*pi4*pi020]);
输入信号x(t)的幅度频谱'
频率弧度/秒'
Y=y'
*exp(-j*t'
Y1=abs(Y);
subplot(326)
plot(w,Y1),
输出信号y(t)的幅度频谱'
subplot(324)
系统的幅度频率响应:
执行程序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)的幅度频谱图如下:
结合所学的有关滤波的知识,根据上面所得到的信号的时域和频域图形,请从时域和频域两个方面解释滤波的概念。
滤波在时域的理解即是输入信号与单位冲激响应的卷积,也即是说,当前时刻的滤波结果,是对当前及之前若干个输入信号加上不同系数的权值,然后再累加的结果。
权系数的不同,代表了不同的滤波特性。
滤波在频域的理解即是首先对原信号f(x,y)进行傅里叶变换,将信号从时域变换到频域,获得频域信号F(u,v),然后用滤波器H(u,v)和F(u,v)相乘,改变原信号的频谱成分,对改变后的频域信号进行傅里叶逆变换,将信号从频域变换回时域,获得滤波后的图像。
四、实验体会
本实验完成时间:
年月日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验三 连续时间LTI系统的频域分析 实验 连续 时间 LTI 系统 分析