系统函数的零极点分布决定时域特性
- 格式:doc
- 大小:397.72 KB
- 文档页数:26
信号与系统课程设计报告--傅里叶变换的对称性和时移特性课程设计任务书2沈阳理工大学摘要本文研究的是傅里叶变换的对称性和时移特性,傅里叶变换的性质有:对称性、线性(叠加性)、奇偶虚实性、尺度变换特性、时移特性、频移特性、微分特性、积分特性、卷积特性(时域和频域);从信号与系统的角度出发,给出了激励信号的具体模型;应用Matlab软件进行仿真,将研究的信号转化成具体的函数形式,在Matlab得到最终变换结果。
使用傅里叶变换的方法、卷积的求解方法以及函数的微分等方法研究题目。
关键词: 傅里叶变换;对称性;时移特性;Matlab3沈阳理工大学目录1、Matlab介绍........................... 错误!未定义书签。
2.利用Matlab实现信号的频域分析—傅里叶变换的对称性与时移特性设计 (5)2.1.傅里叶变换的定义及其相关性质 (5)2.2.傅里叶变换的对称性验证编程设计及实现 (7)2.3.傅里叶变换的时移特性验证编程设计及实现 (11)3.总结 (13)4.参考文献 (13)4沈阳理工大学1、Matlab介绍MATLAB作为一种功能强大的工程软件,其重要功能包括数值处理、程序设计、可视化显示、图形用户界面和与外部软件的融合应用等方面。
MATLAB软件由美国Math Works公司于1984年推出,经过不断的发展和完善,如今己成为覆盖多个学科的国际公认的最优秀的数值计算仿真软件。
MATLAB具备强大的数值计算能力,许多复杂的计算问题只需短短几行代码就可在MATLAB中实现。
作为一个跨平台的软件,MATLAB已推出Unix、Windows、Linux和Mac等十多种操作系统下的版本,大大方便了在不同操作系统平台下的研究工作。
MATLAB软件具有很强的开放性和适应性。
在保持内核不变的情况下,MATLAB 可以针对不同的应用学科推出相应的工具箱(toolbox),目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱,极大地方便了不同学科的研究工作。
实验一离散时间LTI系统的时域分析与Z域分析一、实验目的1、掌握用MATLAB求解离散时间系统的零状态响应、单位脉冲响应和单位阶跃响应;2、掌握离散时间系统系统函数零极点的计算方法和零极点图的绘制方法,并能根据零极点图分析系统的稳定性。
二、实验原理1、离散时间系统的时域分析(1)离散时间系统的零状态响应离散时间LTI系统可用线性常系数差分方程来描述,即MATLAB中函数filter可对式(1-1)的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter的语句格式为:y=filter(b,a,x)其中,x为输入的离散序列;y为输出的离散序列;y的长度与x的长度一样;b与a分别为差分方程右端与左端的系数向量。
(2)离散时间系统的单位脉冲响应系统的单位脉冲响应定义为系统在 (n)激励下系统的零状态响应,用h(n)表示。
MATLAB求解单位脉冲响有两种方法:一种是利用函数filter;另一种是利用函数impz。
impz函数的常用语句格式为impz(b,a,n),其中b和a的定义见filter,n表示脉冲响应输出的序列个数。
(3)离散时间系统的单位阶跃响应系统的单位阶跃响应定义为系统在ε(n)激励下系统的零状态响应。
MATLAB求解单位脉冲响应有两种方法:一种是利用函数filter,另一种是利用函数stepz。
stepz函数的常用语句格式为stepz(b,a,N)其中,b和a的定义见filter,N表示脉冲响应输出的序列个数。
2、离散时间系统的Z域分析(1)系统函数的零极点分析离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比,即如果系统函数H(z)的有理函数表示式为那么,在MATLAB中系统函数的零极点就可通过函数roots得到,也可借助函数tf2zp得到。
roots的语法格式为:Z=roots(b)%计算零点b=[b1b2…bmbm+1]P=roots(a)%计算极点a=[a1a2…anan+1]tf2zp的语句格式为[Z,P,K]=tf2zp(b,a)其中,b与a分别表示H(z)的分子与分母多项式的系数向量。
目录1.引言 (2)2.虚拟仪器开发软件Labview入门 (3)2.1 Labview简介 (3)2.2 利用Labview编程完成习题设计 (3)3.利用LabVIEW实现系统函数的零极点分布决定时域特性的设计 (20)3.1系统函数的零极点分布决定时域特性的基本原理 (20)3.2系统函数的零极点分布决定时域特性的编程设计及实现 (22)3.3运行结果及分析 (23)4. 总结 (25)5.参考文献 (25)1.引言冲激响应h(t)与系统函数H(s) 从时域和变换域两方面表征了同一系统的本性。
在s 域分析中,借助系统函数在s平面零点与极点分布的研究,可以简明、直观地给出系统响应的许多规律。
系统的时域、频域特性集中地以其系统函数的零、极点分布表现出来。
主要优点:可以预言系统的时域特性;便于划分系统的各个分量(自由/强迫,瞬态/稳态);可以用来说明系统的正弦稳态特性。
2.虚拟仪器开发软件Labview入门2.1 Labview简介LabVIEW是一种程序开发环境,由美国国家仪器(NI)公司研制开发的,类似于C 和BASIC开发环境,但是LabVIEW与其他计算机语言的显著区别是:其他计算机语言都是采用基于文本的语言产生代码,而LabVIEW使用的是图形化编辑语言G编写程序,产生的程序是框图的形式。
LabVIEW(Laboratory Virtual Instrument Engineering Workbench)是一种用图标代替文本行创建应用程序的图形化编程语言。
传统文本编程语言根据语句和指令的先后顺序决定程序执行顺序,LabVIEW 则采用数据流编程方式,程序框图中节点之间的数据流向决定VI及函数的执行顺序。
VI指虚拟仪器,是LabVIEW]的程序模块。
LabVIEW 提供很多外观与传统仪器(如示波器、万用表)类似的控件,可用来方便地创建用户界面。
用户界面在LabVIEW中被称为前面板。
使用图标和连线,可以通过编程对前面板上的对象进行控制。
实验八由系统函数零、极点分布决定时域、频域特性实验目的用MA TLAB模拟系统函数零、极点分布决定时域、频域特性的过程实验原理计算零、极点可以用roots函数,若参数为传统函数H(s)的分子多项式系数b,则得到零点;若为分母多项式系数a,则得到极点。
MA TLAB还提供了zero(sys)和pole(sys)函数直接计算零、极点,其中sys表示系统传递函数。
另外,【p,z】=pzmap(sys)函数也具有计算极点p和零点z的功能。
不带返回值的pzmap(sys)则绘制出系统的零、极点分布图。
零极点和传递函数的多项式系数一样,可以作为LTI系统的描述方法。
MATLAB提供了(b,a)=zp2tf(z,p,k)和【z,p,k】=tf2zp(b,a)两个函数用于在上述两种描述方法之间进行转换,其中k为用零、极点表示传函数时的系统增益。
实验过程1.当F(s)极点(一阶)落于图所示s屏幕图中的各方框所处位置时,画出对应的f(t)波形填入方框中。
解:clear all, close all, clc;t = [0:.1:40]';figure, id = 1;for omega = .5:-.25:0for sigma = -.06:.03:.06p = sigma + j*omega;if omega ~= 0p = [p;p'];end[b a] = zp2tf([],p,1);subplot(3,5,id);impulse(b,a,t);set(gca,'YLim',[-20,20]);endend输出结果如下图所示,可见随着极点从虚轴左侧移动到右侧,其冲激响应有衰减变为发散;随着极点由实轴向上、下两侧移动,冲激响应由单调变化转为振荡变化,且振荡周期逐渐减小。
MATLAB知识点(9)-虚数单位i,j以及常量恢复和数学公式中经常以i、j作为下标一样,包括Fortran和C语言在内的打部分机器语言教科书都习惯的把这两个字符作为循环变量使用,但在MATLAB中i和j被默认为定义为表示虚数单位的常量。
2020年硕士研究生统一入学考试《信号与系统》第一部分考试说明一、考试性质信号与系统是信息与通信工程(一级学科)、电子与通信工程领域硕士研究生入学考试的专业基础课。
考试对象为参加东北大学计算机科学与工程学院2020年全国硕士研究生入学考试的准考考生。
二、考试形式与试卷结构(一)答卷方式:闭卷,笔试(二)答题时间:180分钟(三)考试题型及比例基本概念解释15%选择填空 25%计算题30%应用题 30%(四)参考书目郑君里、应启珩、杨为理《信号与系统》 (第三版) 上、下册高等教育出版社 2011.03第二部分考查要点(一)基本概念1、信号的描述、分类和典型示例2、阶跃信号与冲激信号3、信号的分解4、系统的模型及分类(二)连续时间系统的时域分析1、微分方程式的建立与求解2、零输入响应和零状态响应3、冲激响应与阶跃响应4、卷积5、卷积的性质(三)连续时间系统的频域分析1、周期信号的傅立叶分解2、典型周期信号的傅立叶级数3、傅立叶变换4、典型非周期信号的傅立叶变换5、冲激函数和阶跃函数傅立叶变换6、傅立叶变换的基本性质7、卷积特性8、周期信号的傅立叶变换9、抽样信号的傅立叶变换10、抽样定理(四)拉普拉斯变换、连续时间系统的S域分析1、拉普拉斯变换的定义、收敛域2、拉普拉斯变换的性质3、拉普拉斯反变换4、用拉普拉斯变换分析电路、S域元件模型5、系统函数H(s)6、由系统函数的零极点分布决定时域特性7、由系统函数零极点分布决定濒响特性8、连续系统稳定性9、双边拉氏变换10、拉氏变换与付氏变换的关系(五)付里叶变换应用于通信系统1、利用系统函数H(jw)求响应2、无失真传输3、理想低通滤波器4、调制与解调(六)离散系统时域分析1、离散时间信号2、离散时间系统数学模型-差分方程3、常系数线性差分方程的求解4、离散时间系统的单位样值响应5、离散时间系统卷积和(七)Z变换、离散时间系统的Z域分析1、Z变换的定义、典型序列的Z变换2、Z变换的收敛域3、反Z变换4、Z变换的基本性质5、Z变换与拉拉氏变换的关系6、利用Z变换解差分方程7、离散系统函数H(Z)样题:。
目录1.引言 (2)2.虚拟仪器开发软件Labview入门 (3)2.1 Labview简介 (3)2.2 利用Labview编程完成习题设计 (3)3.利用LabVIEW实现系统函数的零极点分布决定时域特性的设计 (20)3.1系统函数的零极点分布决定时域特性的基本原理 (20)3.2系统函数的零极点分布决定时域特性的编程设计及实现 (22)3.3运行结果及分析 (23)4. 总结 (25)5.参考文献 (25)1.引言冲激响应h(t)与系统函数H(s) 从时域和变换域两方面表征了同一系统的本性。
在s 域分析中,借助系统函数在s平面零点与极点分布的研究,可以简明、直观地给出系统响应的许多规律。
系统的时域、频域特性集中地以其系统函数的零、极点分布表现出来。
主要优点:可以预言系统的时域特性;便于划分系统的各个分量(自由/强迫,瞬态/稳态);可以用来说明系统的正弦稳态特性。
2.虚拟仪器开发软件Labview入门2.1 Labview简介LabVIEW是一种程序开发环境,由美国国家仪器(NI)公司研制开发的,类似于C 和BASIC开发环境,但是LabVIEW与其他计算机语言的显著区别是:其他计算机语言都是采用基于文本的语言产生代码,而LabVIEW使用的是图形化编辑语言G编写程序,产生的程序是框图的形式。
LabVIEW(Laboratory Virtual Instrument Engineering Workbench)是一种用图标代替文本行创建应用程序的图形化编程语言。
传统文本编程语言根据语句和指令的先后顺序决定程序执行顺序,LabVIEW 则采用数据流编程方式,程序框图中节点之间的数据流向决定VI及函数的执行顺序。
VI指虚拟仪器,是LabVIEW]的程序模块。
LabVIEW 提供很多外观与传统仪器(如示波器、万用表)类似的控件,可用来方便地创建用户界面。
用户界面在LabVIEW中被称为前面板。
使用图标和连线,可以通过编程对前面板上的对象进行控制。
摘要本文详细分析了系统函数零极点的分布与冲击响应时域特性之间的关系。
首先论述了如何通过MATLAB软件绘制出系统函数的零极点分布图。
然后根据系统函数极点的不同分布情况,通过MATLAB软件绘制出冲击响应的时域函数,通过对图像的观察和比较,得出了极点的类型决定时间函数的时间连续形式,极点在S平面的位置决定时间函数的波形特点。
最后,在极点相同,但零点不同的情况下,通过比较时域函数的波形,得出零点分布与时域函数的对应关系,即零点分布的情况只影响到时域函数的幅度和相位。
关键词:系统函数的零极点;时域特性;MATLAB软件目录1课程设计目的 (1)2实验原理 (1)3实现过程 (1)3.1MATLAB简介 (1)3.2系统函数极点分布情况 (2)3.2.1极点为单实根 (2)3.2.2极点为共轭复根 (2)3.2.3极点为重根 (2)3.2.4用MATLAB绘制系统函数的零极点分布图 (2)3.3系统函数的零极点分布与冲击响应时域特性的关系 (6)3.3.1用MATLAB绘制冲击响应的时域函数 (6)3.3.2极点的类型决定时间函数的时间连续形式 (19)3.3.3极点在S平面的位置决定时间函数的波形特点 (19)3.3.4零点分布与时域函数的对应关系 (19)4设计体会 (23)5参考文献 (24)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)零点和极点的位置,显然,从这些零点和极点的分布情况,便可确定原函数的性质。
设连续系统的系统函数为)(s H ,冲激响应为)(t h ,则⎰+∞-=0)()(dt e t h s H st显然,)(s H 必然包含了)(t h 的本质特性。
对于集中参数的LTI 连续系统,其系统函数可表示为关于s 的两个多项式之比,即其中),,2,1(M j q j =为)(s H 的M 个零点,),,2,1(N i p i =为)(s H 的N 个极点。
3 实现过程3.1 MATLAB 简介MALAB 译于矩阵实验室(MATrix LABoratory ),是用来提供通往 LINPACK 和EISPACK 矩阵软件包接口的。
后来,它渐渐发展成了通用科技计算、图视交互系统和程序语言。
MATLAB 的基本数据单位是矩阵。
它的指令表达与数学、工程中常用的习惯形式十分相似。
比如,矩阵方程Ax=b ,在MATLAB 中被写成A*x=b 。
而若要通过A ,b 求x ,那么只要写x =A \b 即可,完全不需要对矩阵的乘法和求逆进行编程。
因此,用MATLAB 解算问题要比用C 、Fortran 等语言简捷得多。
MATLAB 发展到现在,已经成为一个系列产品:MATLAB “主包”和各种可选的toolbox “工具包”。
主包中有数百个核心内部函数。
迄今所有的三十几个工具包又可分为两类:功能性工具包和学科性工具包。
功能性工具包主要用来扩充MATLAB 的符号计∏∏11)-()-()()()(Ni i Mj j p sq sCs A s B s H ====算功能、图视建模仿真功能、文字处理功能以及硬件实时交互功能。
这种功能性工具包用于多种学科。
而学科性工具包是专业性比较强的,如控制工具包(Control Toolbox )、信号处理工具包(SignalProcessing Toolbox )、通信工具包(Communication Toolbox )等都属此类。
开放性也许是MATLAB 最重要、最受人欢迎的特点。
除内部函数外,所有MATLAB 主包文件和各工具包文件都是可读可改的源文件,用户可通过对源文件的修改或加入自己编写文件去构成新的专用工具包。
MATLA B 已经受了用户的多年考验。
在欧美发达国家,MATLAB 已经成为应用线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具;成为攻读学位的大学生、硕士生、博士生必须掌握的基本技能。
在设计研究单位和工业部门,MATLAB 被广泛地用于研究和解决各种具体工程问题。
3.2 系统函数极点分布情况3.2.1 极点为单实根若系统函数的N 个极点是单极点,则可将)(s H 进行部分分式展开为:3.2.2 极点为共轭复根若系统函数的N 个极点包含共轭复数极点,则可将)(s H 进行部分分式展开为:3.2.3 极点为重根若系统函数的N 个极点有多重极点,则可将)(s H 进行部分分式展开为:3.2.4 用MATLAB 绘制系统函数的零极点分布图绘制系统函数零极点分布图的函数如下:∑1-)(Ni iip s k s H ==))()(()()(βαβαj s j s s D s A s H +++=)()()()1(s D s A s H p s k =function [p,q]=sjdt(A,B)%绘制连续系统零极点图程序%A:系统函数分母多项式系数向量%B:系统函数分子多项式系数向量%p:函数返回的系统函数极点位置行向量%q:函数返回的系统函数零点位置行向量p=roots(A); %求系统极点q=roots(B); %求系统零点p=p';%将极点列向量转置为行向量q=q';%将零点列向量转置为行向量x=max(abs([p q])); %确定纵坐标范围x=x+0.1;y=x;%确定横坐标范围clfhold onaxis([-x x -y y]);%确定坐标轴显示范围axis('square')plot([-x x],[0 0]) %画横坐标轴plot([0 0],[-y y]) %画纵坐标轴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,'实轴')(1)首先打开MATLAB软件,如图1所示。
图1 MATLAB界面(2)点击File-->New-->Function,创建M文件,用于建立一个求解系统函数零极点的函数。
(3)在编辑窗口中输入绘制系统函数零极点分布图的函数,如图2所示。
图2 M文件编辑窗口(4)点击Debug-->Save File and Run设置其保存路径为MATLAB的工作路径,如图3所示。
图3 设置保存路径(5)输入在MATLAB的命令窗口中输入计算系统零极点的具体程序:例如某连续系统的系统函数为:则输入的程序为:a=[10 8 6 4 2];b=[5 4 3];[p,q]=sjdt(a,b)运行结果如图4所示。
图4 绘制系统零极点分布图的运行结果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:p:t2)该调用格式并不绘出系统冲击响应的波形,而是求出向量a和b定义的连续系统在t1~t2时间范围内以时间间隔p抽样的系统冲击响应的数值解。
在MATLAB中分别输入以下各程序,分别绘制出极点在S域的不同分布其对应的时域函数波形。
a.极点为单实根,且在原点处,即绘制冲激响应时域波形的MATLAB命令如下:a=[1 0];b=[1];impulse(b,a)运行结果如图5-1所示。
图5-1 时域函数波形图b. 极点为单实根,且落于原点左侧,即令a=10绘制冲激响应时域波形的MATLAB命令如下:a=[1 10];b=[1];impulse(b,a)运行结果如图5-2所示。
图5-2 时域函数波形图c. 极点为单实根,且落于原点右侧,即令a=10绘制冲激响应时域波形的MATLAB命令如下:a=[1 -10];b=[1];impulse(b,a)运行结果如图5-3所示。
图5-3 时域函数波形图d.极点为一阶共轭复根,且在虚轴上,即令β=5绘制冲激响应时域波形的MATLAB命令如下:a=[1 0 25];b=[1];impulse(b,a,10)运行结果如图5-4所示。
图5-4 时域函数波形图e. 极点为一阶共轭复根,且在虚轴左侧,即令α=0.2,β=5绘制冲激响应时域波形的MATLAB命令如下:a=[1 0.4 25.04];b=[1];impulse(b,a,10)运行结果如图5-5所示。
图5-5 时域函数波形图f. 极点为一阶共轭复根,且在虚轴右侧,即令α=-0.2,β=5绘制冲激响应时域波形的MATLAB命令如下:a=[1 -0.4 25.04];b=[1];impulse(b,a,10)运行结果如图5-6所示。
图5-6 时域函数波形图g. 极点为二阶重根,且在原点处:绘制冲激响应时域波形的MATLAB命令如下:a=[1 0 0];b=[1];impulse(b,a)运行结果如图5-7所示。
图5-7 时域函数波形图h. 极点为二阶重根,且在原点左侧:令a=10绘制冲激响应时域波形的MATLAB命令如下:a=[1 20 100];b=[1];impulse(b,a)运行结果如图5-8所示。
图5-8 时域函数波形图i. 极点为二阶重根,且在原点右侧:令a=10绘制冲激响应时域波形的MATLAB命令如下:a=[1 -20 100];b=[1];impulse(b,a)运行结果如图5-9所示。
图5-9 时域函数波形图j. 极点二阶共轭复根,且在虚轴上:令β=5绘制冲激响应时域波形的MATLAB命令如下:a=[1 0 50 0 625];b=[1];impulse(b,a)运行结果如图5-10所示。