计算机控制系统4教案.docx
- 文档编号:2009652
- 上传时间:2023-05-02
- 格式:DOCX
- 页数:46
- 大小:157.48KB
计算机控制系统4教案.docx
《计算机控制系统4教案.docx》由会员分享,可在线阅读,更多相关《计算机控制系统4教案.docx(46页珍藏版)》请在冰点文库上搜索。
第4章计算机控制系统的数学描述及脉冲传递函数
本章将系统地讲述有关计算机控制系统的数学描述问题,并根据系统的数学描述对系统进行动态分析。
线性离散系统的数学描述形式和线性连续系统的数学描述形式是相对应的,通常有差分方程,脉冲传递函数
(又称Z传递函数),单位脉冲响应序列(又称权序列),离散状态空间表达式等四种数学描述形式。
它们分别与连续系统的四种数学描述形式相对应。
离散状态空间表达式在第9章进行研究。
系统分析和系统设计时所采用的方法有关,一般是,不同的方法采用不同形式的数学描述。
4.1离散系统
离散时间系统(简称离散系统),简单地说就是其输入和输出信号均为离散信号的物理系统。
在数学上,离散系统可以抽象为一种由系统的离散输入信号x(k)到系统的离散输出y(k),(k=0,±1,±2,L)的数学变换或映射。
若将这种变换或映射以符号T[×]表示,则离散系统可表示为
y(k)=T[x(k)]
-32-
(4.1)
离散系统可用框图表示如图4.1所示。
其中x(k)和
x(k)
输入
y(k)
输出
离散系统
图4.1离散系统框图
y(k)分别表示系统的输入和输出在kT时刻的数值。
1.线性离散系统:
如果离散系统的输入信号到输出信号的变换关系满足比例叠加原理,即当输入信号为x(k)=ax1(k)+bx2(k)时,其中a,b为任意常数,系统相应的输出信号可表示为
y(k)=T[x(k)]=aT[x1(k)]+bT[x2(k)]
则该系统就称为线性离散系统。
若不满足比例叠加原理,就是非线性离散系统。
2.时不变离散系统:
是指由输入信号到输出信号之间的变换关系不随时间变化而变化的离散系统,即时不变离散系统应满足如下关系,若y(k)=T[x(k)],那么当系统输入信号为x(k-n)时,则相应的输出信号为
y(k-n)=T[x(k-n)],n=0,±1,±2,L
时不变离散系统又称为定常离散系统。
3.线性时不变离散系统:
是指系统的输入信号到输出信号之间的变换关系既满足比例叠加原理,同时其变换关系又不随时间变化而变化的离散系统。
工程中大多数计算机控制系统可以近似为线性时不变离散系统来处理。
所以本书以后的论述仅限于线性时不变离散系统。
4.2差分方程
线性时不变离散系统的基本数学描述是常系数线性差分方程。
差分方程有前向差分方程与后向差分方程之分。
为了方便,这里提到差分方程,若无特别说明,均指线性常系数差分方程;系统是指线性时不变离散系统。
4.2.1线性常系数差分方程
x(k)
y(k)
输入
输出
离散系统
图4.1离散系统框图
(图4.1重绘与此)
设有一单输入、单输出的线性时不变离散系统,如图
4.1所示。
显然,在某一采样时刻的系统输出值y(k)不仅与
该时刻的输入值x(k)有关,而且与过去时刻的输入值x(k-1),
x(k-2),L有关,还和过去时刻的输出值y(k-1),y(k-2),L有关。
这种关系可以描述为
y(k)+a1y(k-1)+a2y(k-2)+L+any(k-n)
=b0x(k)+b1x(k-1)+b2x(k-2)+L+bmx(k-m)
(4.2)
式中,a1,a2,Lan,b0,b1,b2,Lbm均为实常数,n为方程的阶次。
因此,式(4.2)称为n阶后向非齐次差分
方程。
对于n阶差分方程,an
¹0,其余系数a,La
n-1都有
1
可能为零。
若an
=0,就相当于方程的阶次降为n-1阶。
若b0=0则相应离散系统有一拍(即一个采样周期T)的延迟,即系统在k时刻输出y(k)只与k以前各时刻的输入x(k-i),i=1,2,L,m有关,而与当前时刻的输入值x(k)无
关。
若b0
=b1
=L=bl
=0,则相应离散系统存在l+1拍
延迟,即系统当前时刻的输出y(k)只与(k-l)以前时刻的输入(k-l-1),L,(k-l-m)有关。
与方程(4.2)类似,非齐次n阶前向差分方程基本形式为
y(k+n)+a1y(k+n-1)+L+any(k)
=b0x(k+m)+b1x(k+m-1传+L+bmx(k)
(4.3)
式中,a1,a2,Lan,b0,b1,b2,Lbm均为实常数。
对于有因果关系的物理系统,方程中总是m£n。
若m>n,将式(4.3)的两边右移n拍,即
y(k)+a1y(k-1)+L+any(k-n)
=b0x(k+m-n)+b1x(k+m-n-1传+L+bmx(k-n)
上式右边第一项为b0x(k+m-n),令m-n=d
,即有
b0x(k+d)项,说明当前的输出y(k)与未来输入x(k+d)有关,即不是因果关系。
表明方程描述的离散系统输出信号超前于输入信号,即输入信号尚未作用于系统,其对应的输出信号就已出现,或者说系统当前时刻的输出y(k)与未来时刻输入值
x(k+d),d>0有关。
这种情况在现实的物理系统是不可能出现的。
当m 工程上差分方程都是采用其标准形式如方程(4.2)和 (4.3)形式,至于前向差分方程和后向差分方程,并无本质区别,前向差分方程多用于描述非零初始值的离散系统,而后向差分多用于描述全零初始值的离散系统。 若不考虑 系统初始值,就系统输入与输出关系而言两者完全等价,可以相互转换。 4.2.2差分方程求解 差分方程求解,就是在系统初始值(即系统输入、输出的初始值)和输入序列已知的条件下, 求解差分方程描述的系统在任何时刻的输出序列值。 差分方程解的形式与微分方程解相似,非齐次差分方程的解是由通解加特解组成的。 通解表示方程描述的离散系统在输入为零情况下 (即无外界作用)由系统非零初始值所引起的自由运动,它反映系统本身所固有的动态特性;特解表示方程描述的离散系统在外界输入作用下所产生的强迫运动,它既与系统本身的动态特性有关,又与外界输入作用有关,但与系统初始值无关。 求解线性时不变差分方程有三种基本方法,即经典解法、计算机迭代编程法以及Z变换法。 1.差分方程的经典解法 与(4.2)式相应的齐次方程为 y(k)+a1y(k-1)+a2y(k-2)+L+any(k-n)=0 (4.4) 设满足齐次方程(4.4)的通解具有Aak的形式,并将其代入(4.4)式,得 1 Aak+a Aak-1+a Aak-2+L+a Aak-n=0 2 n Aak(1+aa-1+aa-2+L+aa-n)=0 1 2 n 因为Aak¹0,故得 1+aa-1+aa-2+L+aa-n=0 1 2 n (4.5) 用an乘(4.5)式的两边,得 an+aan-1+aan-2 +L+a =0 1 2 n (4.6) 式(4.6)称为齐次方程(4.4)的特征方程。 若特征方程有两两相异的特征根a1,a2,L,an,则 y(k)=Cak +Cak +L+Cak (4.7) 11 2 2 n n 式(4.7)称为齐次方程(4.4)的通解。 式中系数 C1,C2,LCn由初始值求出。 若特征方程有a1的m重特征根,那么差分方程的通解为 y(k)=(Ckm-1+Ckm-2+L+C )ak+C ak+L+Cak-m+1 1 2 (4.8) m 1 m+12 n n 若ai为复数或虚数时,总是成对出现,每对复数特征值所对应的自由运动分量呈现衰减振荡或发散振荡,每对虚特征值所对应的自由运动分量呈现等幅振荡。 当差分方程包含输入作用时,称该方程为非齐次方程。 非齐次方程的特解与微分方程的特解类似,特解的形式要经过试探才能确定。 表4.1列出了非齐次差分方程常见的特解形式。 表4.1非齐次差分方程常见的特解形式 输入量x(k) 输出量y(k) km pkm+pkm-1+L+p 0 1 m ak a不是差分方程的任何特征 根 pak a是差分方程 的特征根之 一 相异根 pkak+pak 1 2 m-1次重根 pkm-1ak+pkm-2ak+L+pak 1 2 m 例4.1求解差分方程 y(k)+2y(k-1)=2k-1 (4.9) 的解y(k)。 初始值y(0)=1。 解该差分方程的特征方程为 所以,特征根为 a+2=0 a=-2 因此差分方程的通解为 C(-2)k 设差分方程的特解为 p1k+ p2代入(4.9)式,得 p1k+ p2+2[p1(k-1)+ p2]=2k-1 3p1k-2p1+3p2 (4.10) =2k-1 比较(4.10)式两边系数,得 2 p1=3, p =1 2 9 则,差分方程的特解为 2k+1 3 9 因此,差分方程的全解为 y(k)=C(-2)k+2k+1 3 9 利用初始值y(0)=1代入上式, 1=C+1 9 则 C=8将C值代入上式,可得差分方程的全 9 解为 y(k)=8(-2)k+2k+1 k=0,1,2,L 9 3 9 2.利用计算机编程解差分方程 高阶差分方程不论前向或后向差分方程,都是一种递推算法,任何差分方程都可以用递推算法求解。 现对一般 n阶前向差分方程递推求解予以说明,为便于计算,将n阶 前向差分方程(4.3)改写为 y(k+n)=-a1y(k+n-1)-a2y(k+n-2)-L-any(k) +b0x(k+m)+L+bmx(k) n m =-åaiy(k+n-i)+åbix(k+m-i) (4.11) i=1 i=0 只要知道输出序列初始值y(0),y (1),L,y(n-1)和任何时刻的输入序列x(i),i=0,1,2,L,那么系统任何时刻的输出序列y(k),k³n,都可以由式(4.11)逐步递推计算出来。 例4.2求下列差分方程的解y(k)。 y(k)+y(k-1)=x(k)-x(k-1), k³0 式中 ì 1, í x(k)= î0传 k传传 k奇数 且 y(-1)=x(-1)=0。 解①令k=0,1,2,L一步一步迭代解差分方程。 y(0)=-y(-1)+x(0)-x(-1)=0+1-0=1y (1)=-y(0)+x (1)-x(0)=-1+0-1=-2y (2)=-y (1)+x (2)-x (1)=2+1-0=3y(3)=-y (2)+x(3)-x (2)=-3+0-1=-4 y(4)=-y(3)+x(4)-x(3)=4+1+0=5 M ②利用MATLAB语言编程 y (1)=0,x (1)=0; ;MATLAB语言数组编号由1开始,这里y (1)=0,x (1)=0表示实际y(-1)=0,x(-1)=0。 FOR i=2: n ;n为要计算y(k)的拍数,根据需要 设置。 x(i)=1-x(i-1); ;输入 1, í ì x(k)= î0传 k传传 k奇数 y(i)=-y(i-1)+x(i)-x(i-1); ;用迭代式计算y(k)的 值。 end 以上程序计算出y(k)的值,k³0(k=i-2)。 y(0)=1y (1)=-2y (2)=3y(4)=-4y(5)=5 M 差分方程的迭代算法,虽然计算简明,不需要更多的数学知识,但它只能计算出有限个序列值,在一般情况下,得不到方程解的解析表达式,即系统输出序列的一般项表达式。 3.Z变换法 用Z变换方法解差分方程同用拉氏变换求解微分方程类似,其步骤如下: ①利用Z变换线性性质和位移定理对差分方程两边分别进行Z变换,将差分方程变为以Z为变量的代数方程; ②代入系统初始值,通过同类项合并、整理,得到输出Z变换Y(z)的表达式; ③对已知的输入序列进行Z变换,并将其Z变换代入 Y(z)的表达式中,使Y(z)成为确定的Z的函数; ④对Y(z)进行Z反变换,求得相应的输出序列y(k)的表达式。 例4.3试用Z变换方法求解例4.2差分方程。 y(k)+y(k-1)=x(k)-x(k-1),k³0 (4.12) 差分方程输入 1, í ì x(k)= î0传 k传传 k奇数 ,且y(-1)=x(-1)=0。 解 利用Z变换的实平移定理对式(4.12)两边求Z 变换,得 Y(z)+z-1Y(z)=X(z)-z-1X(z) 即 Y(z)=z-1X(z)z+1 因为输入 1, í ì x(k)= î0传 -2 k传传 k奇数 -4 1 z2 所以 X(z)=1+z z-1 +z +L=1-z-2 z2 z2 =z2-1 因而 Y(z)= z+1×z2-1= (z+1)2 用反演积分法求Y(z)的反变换y(k), Y(z)z k-1 = zk+1 (z+1)2 上式含有p=-1二重极点,所以 y(k)=limd(zk+1)=lim[(k+1)zk]=(k+1)(-1)k k=0,1,2,L z®-1dz z®-1 y(0)=1y (1)=-2y (2)=3y(3)=-4 M 与例4.2所求结果是一样的。 上例中,没有给出初始值,因而对后向差分方程作Z变换。 如果给出输出的初始值,则应对前向差分方程求Z变换。 例4.4 试用Z变换方法求解例4.2差分方程,本例给出初始值。 y(k)+y(k-1)=x(k)-x(k-1),k³0 (4.13) 差分方程输入 1, í ì x(k)= î0传 k传传 k奇数 ,且y(0)=2。 解 因为本例中只给出了一个初始值y(0)=2,因此将 (4.13)式前移一步,写成前向差分方程的形式,即 y(k+1)+y(k)=x(k+1)-x(k) (4.14) 对(4.14)式两边作Z变换,得 zY(z)-zy(0)+Y(z)=zX(z)-zx(0)-X(z) 代入初始值y(0)=2,x(0)=1,整理得 Y(z)= (4.15) z-1X(z)+ z+1 z z+1 上式右边第一项为差分方程的特解Z变换,表示系统在外界输入作用下的强迫运动;右边第二项为差分方程的通解 Z变换,表示系统由初始值引起的自由运动。 显然,上式右边两项的分母多项式就是差分方程(4.14)的特征多项式。 由例4.3可知,X(z)= z2 z2-1 ,代入式(4.15),得 Y(z)= (4.16) z-1× z+1 z2 + z2-1 z = z+1 z(2z+1)(z+1)2 用反演积分法求上式的Z反变换 Y(z)z k-1 =(2z+1)zk (z+1)2 上式右边表达式含有z=-1二重极点,因此 dé 2zk(2z+1)ù y(k)=lim ê(z+1) 2 ú z®-1dzë y(k)=lim(k+2)zkz®-1 (z+1) û =(k+2)(-1)k k=0,1,2,L (4.17) y(0)=2 y (1)=-3 M 由差分方程的解(4.17)表达式可以求出y(0)=2,即为本例给出的初始值。 4.3脉冲传递函数(Z传递函数) 我们仿效连续系统理论的思路,利用Z变换引出离散系统脉冲传递函数的概念。 脉冲传递函数可以简称为Z传递函数,为后面进一步研究离散系统复域分析与设计方法奠定基础。 4.3.1脉冲传递函数定义 在连续系统理论中,传递函数定义为: 连续系统在初始静止状态下,即系统初始值为零,系统输出信号的拉氏变换Y(s)与对应的输入信号拉氏变换X(s)之比,即 G(s)= Y(s) X(s) (4.18) 与此类似,线性离散控制系统中,在初始值为零的条件下,一个系统(或环节)输出序列Z变换Y(z)与输入序列Z变换X(z)之比,定义为该系统(或环节)的脉冲传递函数。 (4.19) G(z)= Y(z) X(z) 由脉冲传递函数的定义式(4.19)可知,离散系统输出信号的Z变换可以表示为系统的脉冲传递函数与输入信号Z变换的乘积,即 Y(z)=G(z)X(z) X(z) Y(z) 输入 输出 G(z) 图4.2离散系统方框图 (4.20) 离散系统也可以采用方框图直观表示形式,如图4.2所示。 脉冲传递函数是离散系统动态特性的一种数学描述形式,它只与系统本身特性有关,而与外部输入形式无关。 此外,脉冲传递函数和连续系统传递函数一样,仅适用于线性、时不变系统,而不适用于非线性和时变系统。 离散系统的脉冲传递函数可以由描述离散系统的差分方程通过Z变换求出或离散系统的单位脉冲响应序列通过 Z变换求出。 4.3.2脉冲传递函数的求取 脉冲传递函数是离散系统在Z域的描述形式,差分方程则是离散系统特性在时间域的描述形式,若不考虑系统初始值的作用,两者是相互对应的,可以互相转换。 若给定系统的差分方程,只要令系统初始值为零,对差分方程等号两边分别作Z变换,通过整理,便可获得相应的脉冲传递函数。 1、由离散系统差分方程求脉冲传递函数 设离散系统的差分方程为 y(k+n)+a1y(k+n-1)+L+any(k) =b0x(k+m)+b1x(k+m-1)+L+bmx(k) (4.21) 令系统初始值均为零,即 y(i)=0,i=0,1,2,L,n-1;x(i)=0,i=0,1,2,L,m-1,利用前向实 平移定理,对上式两边作Z变换,得 znY(z)+azn-1Y(z)+L+aY(z) 1 n =bzmX(z)+bzm-1X(z)+L+bX(z) 0 1 m (4.22) 整理后,便得到输出Z变换与输入Z变换之比,即脉冲传递函数为 Y(z) bzm+bzm-1+L+b G(z)= X(z) =0 1 m zn+azn-1+L+a 1 n (4.23) 对于后向差分方程 (4.24) y(k)+a1y(k-1)+L+any(k-n) =b0x(k)+b1x(k-1)+L+bmx(k-m) 由于后向差分方程我们已经约定只描述初始值为零的 系统,所以不必考虑初始值,直接利用实平移延迟定理对方程中各项作Z变换,得 Y(z)+az-1Y(z)+L+az-nY(z) 1 n =bX(z)+bz-1X(z)+L+bz-mX(z) 0 1 m (4.25) 整理后,得相应的脉冲传递函数为 Y(z) b+bz-1+L+bz-m G(z)= =0 1 m 1 n (4.26) X(z) 1+az-1+L+az-n 由上式可以看出,脉冲传递函数有两种形式: 一种是 (4.23)式的形式,为复变量Z的有理形式;另一种则是 (4.26)式的形式,为复变量z-1的有理分式形式。 这两种形式是等价的,可以相互转换。 脉冲传递函数也可以转换为相应的差分方程,假设给定脉冲传递函数G(z)如(4.23)式,则首先化成如下方程 znY(z)+azn-1Y(z)+L+aY(z) 1 n =bzmX(z)+bzm-1X(z)+L+bX(z) 0 1 m 再利用零初始条件下输出、输入序列与其Z变换之间的对 应关系,进一步得到相应的差分方程为 y(k+n)+a1y(k+n-1)+L+any(k) =b0x(k+m)+b1x(k+m-1)+L+bmx(k) 若给定的脉冲传递函数G(z)如式(4.26)形式,即是 z-1的有理分式形式,通过类似操作,同样可以转化为 (4.24
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机控制 系统 教案