连续时间系统的变换域分析系统函数的零极点分布决定时域特性文档格式.docx
- 文档编号:22465125
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:16
- 大小:97.82KB
连续时间系统的变换域分析系统函数的零极点分布决定时域特性文档格式.docx
《连续时间系统的变换域分析系统函数的零极点分布决定时域特性文档格式.docx》由会员分享,可在线阅读,更多相关《连续时间系统的变换域分析系统函数的零极点分布决定时域特性文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
3、检查编程、运行结果、答辩
4、上交课程设计报告
指导教师:
2014年6月日
专业负责人:
学院教学副院长:
1课程设计目的
1.掌握系统函数的零极点分布与系统冲激响应时域特性之间的关系。
2.学习MATLAB软件知识及应用。
3.利用MATLAB编程,完成相应的信号分析和处理。
2实验原理
拉普拉斯变换将时域函数f(t)变换为s域函数F(s);
反之,拉普拉斯逆变换将F(s)变换为相应的f(t)。
由于f(t)与F(s)之间存在一定的对应关系,故可以从函数F(s)的典型形式透视出f(t)的内在性质。
当F(s)为有理函数时,其分子多项式和分母多项式皆可分解为因子形式,各项因子指明了F(s)零点和极点的位置,显然,从这些零点和极点的分布情况,便可确定原函数的性质。
设连续系统的系统函数为
,冲激响应为
,则
显然,
必然包含了
的本质特性。
对于集中参数的LTI连续系统,其系统函数可表示为关于s的两个多项式之比,即
其中
为
的M个零点,
的N个极点。
3实现过程
3.1MATLAB简介
MALAB译于矩阵实验室(MATrixLABoratory),是用来提供通往LINPACK和EISPACK矩阵软件包接口的。
后来,它渐渐发展成了通用科技计算、图视交互系统和程序语言。
MATLAB的基本数据单位是矩阵。
它的指令表达与数学、工程中常用的习惯形式十分相似。
比如,矩阵方程Ax=b,在MATLAB中被写成A*x=b。
而若要通过A,b求x,那么只要写x=A\b即可,完全不需要对矩阵的乘法和求逆进行编程。
因此,用MATLAB解算问题要比用C、Fortran等语言简捷得多。
MATLAB发展到现在,已经成为一个系列产品:
MATLAB“主包”和各种可选的toolbox“工具包”。
主包中有数百个核心内部函数。
迄今所有的三十几个工具包又可分为两类:
功能性工具包和学科性工具包。
功能性工具包主要用来扩充MATLAB的符号计算功能、图视建模仿真功能、文字处理功能以及硬件实时交互功能。
这种功能性工具包用于多种学科。
而学科性工具包是专业性比较强的,如控制工具包(ControlToolbox)、信号处理工具包(SignalProcessingToolbox)、通信工具包(CommunicationToolbox)等都属此类。
开放性也许是MATLAB最重要、最受人欢迎的特点。
除内部函数外,所有MATLAB主包文件和各工具包文件都是可读可改的源文件,用户可通过对源文件的修改或加入自己编写文件去构成新的专用工具包。
MATLAB已经受了用户的多年考验。
在欧美发达国家,MATLAB已经成为应用线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具;
成为攻读学位的大学生、硕士生、博士生必须掌握的基本技能。
在设计研究单位和工业部门,MATLAB被广泛地用于研究和解决各种具体工程问题。
利用Matlab7.0编程完成习题设计
在熟悉了MATLAB7.0的基本界面之后,可以通过简单的编程与相关函数的调用,实现一些常用时间信号的可视化操作。
例如:
编程实现矩形波的仿真。
程序如下,直接在命令窗口键入如下程序:
%rectpuls
t=-4:
0.001:
4;
T=2;
ft=rectpuls(t,T);
plot(t,ft)
axis([-4,4,-0.5,1.5])
仿真图形如下:
3.2系统函数极点分布情况
3.2.1极点为单实根
若系统函数的N个极点是单极点,则可将
进行部分分式展开为:
3.2.2极点为共轭复根
若系统函数的N个极点包含共轭复数极点,则可将
3.2.3极点为重根
若系统函数的N个极点有多重极点,则可将
3.2.4用MATLAB绘制系统函数的零极点分布图
绘制系统函数零极点分布图的函数如下:
function[p,q]=sjdt(A,B)
%绘制连续系统零极点图程序
%A:
系统函数分母多项式系数向量
%B:
系统函数分子多项式系数向量
%p:
函数返回的系统函数极点位置行向量
%q:
函数返回的系统函数零点位置行向量
p=roots(A);
%求系统极点
q=roots(B);
%求系统零点
p=p'
;
%将极点列向量转置为行向量
q=q'
%将零点列向量转置为行向量
x=max(abs([pq]));
%确定纵坐标范围
x=x+0.1;
y=x;
%确定横坐标范围
clf
holdon
axis([-xx-yy]);
%确定坐标轴显示范围
axis('
square'
)
plot([-xx],[00])%画横坐标轴
plot([00],[-yy])%画纵坐标轴
plot(real(p),imag(p),'
x'
)%画极点
plot(real(q),imag(q),'
o'
)%画零点
title('
连续系统零极点图'
)%标注标题
text(0.2,x-0.2,'
虚轴'
text(y-0.2,0.2,'
实轴'
3.3系统函数的零极点分布与冲击响应时域特性的关系
3.3.1用MATLAB绘制冲击响应的时域函数
函数impulse()将绘出由向量a和b表示的连续系统在指定时间范围内的冲击响应h(t)的时域波形图,并能求出指定时间范围内冲击响应的数值解。
Impulse()函数有如下几种调用格式:
(1)impulse(b,a)
该调用格式以默认方式绘出向量a和b定义的连续系统的冲击响应的时域波形图。
(2)impulse(b,a,t)
该调用格式以默认方式绘出向量a和b定义的连续系统在0~t时间范围冲击响应的时域波形图。
(3)impulse(b,a,t1:
p:
t2)
该调用格式以默认方式绘出向量a和b定义的连续系统在t1~t2时间范围内,且以时间间隔p均匀抽样的冲击响应的时域波形。
(4)y=impulse(b,a,t1:
该调用格式并不绘出系统冲击响应的波形,而是求出向量a和b定义的连续系统在t1~t2时间范围内以时间间隔p抽样的系统冲击响应的数值解。
在MATLAB中分别输入以下各程序,分别绘制出极点在S域的不同分布其对应的时域函数波形。
a.极点为单实根,且在原点处,即:
Oσ
绘制冲激响应时域波形的MATLAB命令如下:
a=[10];
b=[1];
impulse(b,a)
运行结果如图5-1所示:
图5-1时域函数波形图
b.极点为单实根,且落于原点左侧,即:
aOσ
令a=10
a=[110];
运行结果如图5-2所示:
图5-2时域函数波形图
c.极点为单实根,且落于原点右侧,即:
Oaσ
a=[1-10];
运行结果如图5-3所示:
图5-3时域函数波形图
d.极点为一阶共轭复根,且在虚轴上,即:
Oσ
令β=5
a=[1025];
impulse(b,a,10)
运行结果如图5-4所示:
图5-4时域函数波形图
e.极点为一阶共轭复根,且在虚轴左侧,即:
令α=0.2,β=5
a=[10.425.04];
运行结果如图5-5所示:
图5-5时域函数波形图
f.极点为一阶共轭复根,且在虚轴右侧,即:
令α=-0.2,β=5
a=[1-0.425.04];
运行结果如图5-6所示:
图5-6时域函数波形图
g.极点为二阶重根,且在原点处,即:
a=[100];
运行结果如图5-7所示:
图5-7时域函数波形图
h.极点为二阶重根,且在原点左侧,即:
a=[120100];
运行结果如图5-8所示:
图5-8时域函数波形图
i.极点为二阶重根,且在原点右侧,即:
a=[1-20100];
运行结果如图5-9所示:
图5-9时域函数波形图
j.极点二阶共轭复根,且在虚轴上,即:
Oσ
a=[10500625];
运行结果如图5-10所示:
图5-10时域函数波形图
k.极点二阶共轭复根,且在虚轴左侧,即:
令α=2,β=5
a=[1874232841];
运行结果如图5-11所示:
图5-11时域函数波形图
l.极点二阶共轭复根,且在虚轴右侧,即:
令α=-2,β=5
a=[1-874-232841];
运行结果如图5-12所示:
图5-12时域函数波形图
3.3.2极点的类型决定时间函数的时间连续形式
由图5-1~图5-12得:
当系统函数的极点为一阶极点时,若极点是实根且位于原点处,时域函数为阶跃函数;
若极点是实根,但不在原点处,时域函数为指数函数。
若极点是共轭复根且位于原点处,时域函数为正弦函数;
若极点是共轭复根,但不在原点处,时域函数为指数函数与正弦函数的乘积。
当系统函数的极点为二阶极点时,若极点是实根且位于原点处,时域函数为一次函数;
若极点是实根,但不在原点处,时域函数为一次函数与指数函数的乘积。
若极点是共轭复根且位于原点处,时域函数为一次函数与正弦函数的乘积;
若极点是共轭复根,但不在原点处,时域函数为一次函数、指数函数和正弦函数的乘积。
不难推论得,若H(s)具有多重极点,那么,部分分式展开式各项所对应的时间函数可能具有等与指数函数相乘的形式,t的幂次由极点阶次决定。
3.3.3极点在S平面的位置决定时间函数的波形特点
由图5-1~图5-12可以看出:
若H(s)极点落于左半平面,则h(t)波形为衰减形式;
若H(s)极点落于右半平面,则h(t)增长;
落于虚轴上的一阶极点对应的h(t)成等幅震荡或阶跃;
而虚轴上的二阶极点将使h(t)呈增长形势。
3.3.4零点分布与时域函数的对应关系
以上分析了H(s)极点分布与时域函数的对应关系。
至于H(s)零点分布的情况则只影响到时域函数的幅度和相位;
s平面中零点变动对于t平面波形的形式没有影响。
其H(s)的零极点分布以及h(t)波形如图6-1,6-2所示:
图6-1H(s)的零极点分布图
图6-2h(t)波形图
由图6-2知,h(t)波形呈衰减震荡形式。
假定保持极点不变,只移动零点的位置,使零点有-0.5变为0,则H(s)变为:
则改变后的H(s)的零极点分布以及h(t)波形如图7-1,7-2所示:
图7-1改变后的H(s)的零极点分布图
图7-2改变后的h(t)波形图
比较可得,保持极点不变,只移动零点的位置,那么h(t)波形将仍呈衰减震荡形式,震荡频率也不改变,只是幅度和相位有变化。
4设计体会
通过本次课程设计,让我更加深刻地理解和掌握了系统函数零极点的分布与冲击响应时域特性之间的关系,也使我了解了MATLAB软件,并学会了如何应用MATLAB软件进行数字信号的处理和分析。
本课程设计利用MATLAB软件成功绘制出了系统函数的零极点分布图。
根据系统函数零极点的不同分布情况,利用MATLAB软件绘制出了冲击响应的时域函数波形。
通过分析和比较,最终得出:
系统函数的极点分布与时域特性的关系,即极点的类型决定时间函数的时间连续形式,极点在S平面的位置决定时间函数的波形特点;
系统函数的零点分布与时域特性的关系,即系统函数零点分布的情况只影响到时域函数的幅度和相位。
通过设计研究,我对这些理论知识有了全新的认识,因为通过自学,可以找出自己想要理解的内容。
我也发现,在设计之前,我们必须把理论知识打扎实,否则很难对题目进行整体的把握。
这次的课程设计也提醒着我要在理论知识上花点功夫,为以后毕业论文撰写打下基础,我也希望自己以后把这些知识运用到实践中去。
在设计之前,我阅读了一些书籍来获得一定的感性知识,然后才有自己的想法和观点,这次的题目并不难,最大的困难是步骤多,考验的是耐心。
通过设计研究,我还了解了运用MATLAB设计的基本步骤,也了解了MATLAB方便,简洁,高效的优点,掌握此软件定能为以后的学习打下了基础。
在研究由系统函数零、极点分布决定时域特性的题目时,由于自己的粗心,导致坐标标错,辛苦画出的5张图全部得重新画,浪费大量的时间和精力。
粗心的结果很残酷,在以后的学习生活中,我会记住这次的教训。
在设计中,我也遇到了些许困难,例如,刚开始设计题目的意思理解错误,走入了一个很复杂的歪路,导致的结果是不知道如何分析波形,停滞在那不知如何,经过同学的指点,才正确理解了题目的意思。
这个教训也在告诉这我,不能钻牛角尖,抓紧某一点不放手,应该开拓思路,可以与同学讨论,也可以查阅相关书籍,从而得到灵感。
此外,本次课程设计报告是按照毕业设计报告的要求来的,因此我抓住此次机会,对word的使用作了进一步的熟悉,学会了公式产生,目录的制作以及文档的排版等,如果在毕业设计时再去熟悉这些小细节,势必会花费较多的时间,所以在平时应该注意这些细节。
总之,这次实训受益颇多,在以后的学习中,我会好好利用这次的学习经验,结合实际在专业领域内进行更深入的学习。
5.参考文献
[1]郑君里.信号与系统引论.第1版.北京:
高等教育出版社,2009:
191-214
[2]梁虹,梁洁,陈跃斌.信号与系统分析及MATLAB实现.第1版.北京:
电子工业出版社,2002:
240-244
[3]陈泽,占海明.详解MATLAN在科学计算中的应用.第1版.北京:
电子工业出版社,2011:
87-93
[4]龚纯,王正林.MATLAB语言常用算法程序集.第2版.北京:
103-121
[5]郑君里(杨为理).信号与系统[M](第二版)上册北京:
高等教育出版社,2009.
[6]张志涌(杨祖樱).MATLAB教程[M]北京:
北京航空航天大学出版社,2009.
[7]徐金明(张孟喜、丁涛).MATLAB实用教程[M]北京:
清华大学出版社,2009.
[8]党宏社.信号与系统实验(MATLAB版)[M]西安:
西安电子科技大学出版社,2007.
[9]刘会灯(朱飞).信号与系统:
MATLAB综合实验[M]北京:
高等教育出版社,2008.
[10]张小虹.《信号与系统》.西安.西安电子科技大学出版社.2008年
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 连续 时间 系统 变换 分析 函数 极点 分布 决定 时域 特性