西安石油大学电子工程学院：《自动控制理论 Modern Control System》精品课程教学资源（电子讲稿）第二章 控制系统的数学模型（1/2）

第二章控制系统的数学模型 2.1引言 控制系统是由控制对象、执行机构、放大器、检测(测量)装置和控制器等组成。 数学模型:在分析和设计系统时,了解这些砖砖瓦瓦的工作原理及运动过程是很重要的, 更重要的是深入研究它们的动态特性,正确列写出它们的数学表达式 我们把描述系统或元件的动态特性的数学表达式叫做系统或元件的数学模型。 深入了解元件及系统的动态特性,淮确建立它们的数学模型一称建模,只有得到较为准确的 数学建模,才能设计出性能良好的控制系统。 动态特性控制系统所采用的元件种类繁多,虽然各自服从的规律,但它们有一共同点 即任何系统或元件总有物质或能量流入,同时又有某些物质或能量流出,系统通常又是有贮 存物质或能量的能力,贮存量的多少用状态变量来表示。状态变量是反应系统流入量或流出 量之间平衡的物理量,由于外部供给系统的物质或能量的速率是有限的,不可能是无穷大 因此,系统的状态变量有一个状态变到另一个状态不可能瞬间完成,而要经过一段时间。这 样,状态变量的变化就有一个过程,这就是动态过程。例如,电路中电容上的电压是一个状 态变量,它由一个值变到另一个值不可能瞬间完成。具有一定惯量的物体的转速是一个状态 变量,转速的变化也是一个过渡过程,具有一定质量的物体的温度是一个状态变量,它由温 度T0变到T,同样有一个动态过程;又如容器中液位也是一个状态变量,液位的变化也要 定的时间 物理模型任何元件或系统实际上都是很复杂的,难以对它作出精确、全面的描述,必 须进行简化或理想化。简化后的元件或系统为该元件或系统的物理模型。简化是有条件的, 要根据问题的性质和求解的精确要求,来确定出合理的物理模型 建立控制系统数学模型的方法有 分析法一对系统各部分的运动机理进行分析,物理规律、化学规律 建立系统数学模型的几个步骤 建立物理模型。 列写原始方程。利用适当的物理定律一如牛顿定律、基尔霍夫电流和电压定律、能 量守恒定律等) 选定系统的输入量、输出量及状态变量(仅在建立状态模型时要求),消去中间变 量,建立适当的输入输出模型或状态空间模型。 实验法一人为施加某种测试信号,记录基本输出响应。 输入(已知) 输出(已知) 黑匣子 数学模型的逼近 数学模型有多种表现形式

14 第二章 控制系统的数学模型 2.1 引言 控制系统是由控制对象、执行机构、放大器、检测（测量）装置和控制器等组成。 数学模型： 在分析和设计系统时，了解这些砖砖瓦瓦的工作原理及运动过程是很重要的， 更重要的是深入研究它们的动态特性，正确列写出它们的数学表达式。 我们把描述系统或元件的动态特性的数学表达式叫做系统或元件的数学模型。 深入了解元件及系统的动态特性，准确建立它们的数学模型－称建模，只有得到较为准确的 数学建模，才能设计出性能良好的控制系统。 动态特性 控制系统所采用的元件种类繁多，虽然各自服从的规律，但它们有一共同点： 即任何系统或元件总有物质或能量流入，同时又有某些物质或能量流出，系统通常又是有贮 存物质或能量的能力，贮存量的多少用状态变量来表示。状态变量是反应系统流入量或流出 量之间平衡的物理量，由于外部供给系统的物质或能量的速率是有限的，不可能是无穷大， 因此，系统的状态变量有一个状态变到另一个状态不可能瞬间完成，而要经过一段时间。这 样，状态变量的变化就有一个过程，这就是动态过程。例如，电路中电容上的电压是一个状 态变量，它由一个值变到另一个值不可能瞬间完成。具有一定惯量的物体的转速是一个状态 变量，转速的变化也是一个过渡过程，具有一定质量的物体的温度是一个状态变量，它由温 度T0变到T，同样有一个动态过程；又如容器中液位也是一个状态变量，液位的变化也要一 定的时间。 物理模型 任何元件或系统实际上都是很复杂的，难以对它作出精确、全面的描述，必 须进行简化或理想化。简化后的元件或系统为该元件或系统的物理模型。简化是有条件的， 要根据问题的性质和求解的精确要求，来确定出合理的物理模型。 建立控制系统数学模型的方法有 分析法－对系统各部分的运动机理进行分析，物理规律、化学规律 建立系统数学模型的几个步骤： • 建立物理模型。 • 列写原始方程。利用适当的物理定律—如牛顿定律、基尔霍夫电流和电压定律、能 量守恒定律等） • 选定系统的输入量、输出量及状态变量（仅在建立状态模型时要求），消去中间变 量，建立适当的输入输出模型或状态空间模型。 实验法－人为施加某种测试信号，记录基本输出响应。 黑匣子 输入（已知） 输出（已知） 数学模型的逼近 数学模型有多种表现形式

时1微分方程一输入量和状态变量都是连续的。集总参数 偏微分方程,分布参数 2差分方程一离散系统 复域1传递函数 2结构图一信号流图 三、频率一频率特性,波特l 连续系统一微分方程1线性微分方程线性系统 2常系数线性微分方程线性定常系统线性时变系统 3偏微分方程一分布参数系统 4非线性微分方程 离散系统一差分方程 2.2控制系统的时域数学模型 2.2.1线性元件的微分方程 例2-1图2-1为由一RC组成的四端无源网络。试列写以U1(t)为输入量,U2(t)为输出 量的网络微分方程 RI R2 U1 C1 U2 图2-1RC组成的四端网络 解:设回路电流i1、i2 根据克希霍夫定律,列写方程 U=RA+UC Ue=cja-i2)dr UG=r,i2+u dt ④ 由④、⑤得 du dU 由②导出

15 一、时 1微分方程－输入量和状态变量都是连续的。集总参数 偏微分方程，分布参数 2差分方程－离散系统 二、复域 1传递函数 2 结构图－信号流图 三、频率－频率特性，波特图 连续系统－微分方程 1线性微分方程 线性系统 时 2常系数线性微分方程 线性定常系统 线性时变系统 3偏微分方程－分布参数系统 域 4非线性微分方程 离散系统－差分方程 2.2 控制系统的时域数学模型 2.2.1线性元件的微分方程 例2-1 图2-1为由一RC组成的四端无源网络。试列写以U1（t）为输入量，U2(t)为输出 量的网络微分方程。 U1 R1 R2 C1 C2 U2 图2-1 RC组成的四端网络 解： 设回路电流i1、i2 根据克希霍夫定律，列写方程 1 1 1 Uc1 U  R i  ①   i  i dt C Uc ( ) 1 1 2 1 1 ② c1 2 2 Uc2 U  R i  ③   i dt C Uc 2 2 2 1 ④ U2  Uc2 ⑤ 由④、⑤得 dt dU C dt dU i C c 2 2 2 2  2  由②导出

du i1=C1=,+i2=C1 将i1、i2代入①、③,则得 U1=R1+R2i2+U2 dU du. du R C2=,2)+R2 r,lCi d (R2i2+U,)+c dU+R,G dU,+Ux dt d2U du. R, CRC, dt R,C R,C d+,c dl dUe d-U du R,( d2+(RC+RC2+R2C2)2+U,=U dt 这就是RC组成的四端网络的数学模型,是一个二阶线性微分方程 例2-2试证明图2-2(a)、(b)所示的机、电系统是相似系统(即两系统具有相同的数学模型)。 解:对机械网络:输入为Xr,输出为Ⅺc,根据力平衡,可列出其运动方程式 KI(X-X+B,(Xr-Xe)=K2x+B (B,+B2)X+(K+K2)X=B, X+Kx 对电气网络(b),列写电路方程如下: idt +ri+ CU,=C R,I (RI+R2)1+Uc+U=U 利用②、③、④求出 Ur-(1+Uc RI+ R2-(1+-Rl 代入① 将①两边微分得 (R1+R2)U+(+=RU+U 力-电压相似 机械 电气 B1阻尼R1电阻 K1弹性系数C1 K2弹性系数C2

16 dt dU C dt dU i C dt dU i C c c 2 2 1 2 1 1 1  1    将i1、i2代入①、③，则得 1 1 2 2 Uc2 U  R  R i  2 2 2 2 2 2 1 1 1 ( ) U dt dU R C dt dU C dt dU R C c     2 2 2 2 2 1 1 2 2 2 2 [ ( ) ] U dt dU R C dt dU R i U C dt d  R C     2 2 2 2 2 1 2 2 2 1 1 2 2 1 1 2 2 U dt dU R C dt dU R C dt dU R C dt d U  R C R C     2 1 2 2 1 1 1 2 2 2 2 2 1 2 1 2 ( ) U U dt dU R C R C R C dt d U R R C C      这就是RC组成的四端网络的数学模型，是一个二阶线性微分方程。 例2-2 试证明图2-2(a)、(b)所示的机、电系统是相似系统(即两系统具有相同的数学模型)。 解： 对机械网络：输入为Xr，输出为Xc，根据力平衡，可列出其运动方程式 c 2 c 2 r c K1 (Xr - Xc ) B1 (X - X ) K X B X       r r (B1  B2 )Xc  (K1  K2 )Xc  B1 X  K1X   对电气网络(b)，列写电路方程如下：       Ur idt C idt R i C R i 1 1 2 2 1 1 ① Uc2 Uc1 C1Uc1  C2Uc2 ② c 1 c1 U  R i U ③ 1 2 Uc1 Uc2 Ur (R  R )i    ④ 利用②、③、④求出 ) 1 2 1 1 2 (1 ) 2 1 (1 R C C R R Uc C C Ur i       代入① 将①两边微分得 c c r Ur C U R U C C R R U 1 1 1 2 1 2 1 ) 1 1 (  )  (      力-电压相似 机械 电气 B1 阻尼 R1 电阻 B2 R2 K1 弹性系数 1 1 C K2 弹性系数 2 1 C

可见,机系统(a)和电系统(b)具有相同的数学模型,故这些物理系统为相似系统。(即 电系统为即系统的等效网络) 相似系统揭示了不同物理现象之间的相似关系。 为我们利用简单易实现的系统(如电的系统)去研究机械系统. 因为一般来说,电的或电子的系统更容易,通过试验进行研究 例3图2-3所示为电枢控制直流电动机的微分方程,要求取电枢电压Ua(t)(v)为输入量, 电动机转速ωm(t)(rad/s)为输出量,列写微分方程。图中Ra(9)、La(H)分别是电枢电 路的电阻和电感,Mc(N·M)是折合到电动机轴上的总负载转距。激磁磁通为常值 ifa 1a W Ea +(SM-2- Jm, m 图2-3电枢控制直流电动机原理图 解:电枢控制直流电动机的工作实质是将输入的电能转换为机械能,也就是由输入的电枢 电压Ua(t)在电枢回路中产生电枢电流ia(t),再由电流ia(t)与激磁磁通相互作用产生电 磁转距Mm(t),从而拖动负载运动。因此,直流电动机的运动方程可由以下三部分组成。 电枢回路电压平衡方程 di() Ea是电枢反电势,它是当电枢旋转时产生的反电势,其大小与激磁磁通及转速成正 比,方向与电枢电压Ua(t)相反,即 Ea=Ce om (t) Ce一反电势系数(v/rad/s)② 电磁转距方程:Mm(t)=Cmin(t) Cm电动机转距系数(N·M/A)是电动机转距系数 Mm(t)-是由电枢电流产生的电磁转距(N·M) 电动机轴上的转距平衡方程: fm@m(o=M(o-M(o Jm一转动惯量(电动机和负载折合到电动机轴上的)kg·ms2 fm-电动机和负载折合到电动机轴上的粘性摩擦系数(N·m/rad/s) ③、④求出ia(t),代入①同时②亦代入①得: LJ.h+(Lnfm+8人on +(r,f+CmC)o(1) =CU(O-L dM(o) -RM(() 在工程应用中,由于电枢电路电感La较小,通常忽略不计,因而⑤可简化为

17 可见，机系统（a）和电系统（b）具有相同的数学模型，故这些物理系统为相似系统。（即 电系统为即系统的等效网络） 相似系统揭示了不同物理现象之间的相似关系。 为我们利用简单易实现的系统（如电的系统）去研究机械系统...... 因为一般来说，电的或电子的系统更容易，通过试验进行研究。 例3 图2-3 所示为电枢控制直流电动机的微分方程，要求取电枢电压Ua(t)（v）为输入量， 电动机转速ωm（t）（rad/s）为输出量，列写微分方程。图中Ra(Ω)、La(H)分别是电枢电 路的电阻和电感，Mc(N·M)是折合到电动机轴上的总负载转距。激磁磁通为常值。 图2-3 电枢控制直流电动机原理图 SM 负 载 ＋ - ＋ - La Ra Ea Wm Jm,fm Ua if ia 解： 电枢控制直流电动机的工作实质是将输入的电能转换为机械能，也就是由输入的电枢 电压Ua(t)在电枢回路中产生电枢电流ia(t)，再由电流ia（t）与激磁磁通相互作用产生电 磁转距Mm(t)，从而拖动负载运动。因此，直流电动机的运动方程可由以下三部分组成。 ·电枢回路电压平衡方程： a a a a a a R i t E dt di t U t  L  ( )  ( ) ( ) ① Ea是电枢反电势，它是当电枢旋转时产生的反电势，其大小与激磁磁通及转速成正 比，方向与电枢电压Ua(t)相反，即 Ea=Ceωm（t） Ce－反电势系数（v/rad/s） ② ·电磁转距方程： Mm(t)=Cmia(t) ③ Cm-电动机转距系数（N·M/A）是电动机转距系数。 Mm(t)-是由电枢电流产生的电磁转距（N·M） ·电动机轴上的转距平衡方程： ( ) ( ) ( ) ( ) f t M t M t dt d t J m m m c m m      ④ Jm－转动惯量（电动机和负载折合到电动机轴上的） kg·m· 2 s fm-电动机和负载折合到电动机轴上的粘性摩擦系数（N·m/rad/s） ③、④求出ia(t)，代入①同时②亦代入①得： ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 R M t dt dM t C U t L R f C C t dt d t L f R J dt d t L J a c c m a a a m m e m m a m a m m a m           ⑤ 在工程应用中，由于电枢电路电感La较小，通常忽略不计，因而⑤可简化为

T don+mn()=K2(0)-k2M2( T RJ 式中 电动机机电时间常数(s) R Rf+CC Rafn +Cmc 电动机传递系数 如果电枢电阻Ra和电动机的转动惯量Jm都很小而忽略不计时⑥还可进一步简化为 C,O(0=U( 电动机的转速Wm(t)与电枢电压Ua(t)成正比,于是电动机可作为测速发电机使用。 系统最基本的数学模型是它的微分方程式。建立微分方程的步骤如下: ①确定系统的输入量和输出量 ②将系统划分为若干环节,从输入端开始,按信号传递的顺序,依据各变量所遵循的物理学 定律,列出各环节的线性化原始方程 ③消去中间变量,写出仅包含输入、输出变量的微分方程式。 2.2.2线性微分方程的求解 初条 列出方程—一求解方程一求解微分方程 输入量 2.2.3非线性元件微分方程的线性化 具有连续变化的非线性函数的线性化,可用切线法或小偏差法。在一个小范围内,将非线 性特性用一断直线来代替。(分段定常系统) 个变量的非线性函数y=f(x) 在x0处连续可微,则可将它在该点附件用台劳级数展开 y=f(x)=f(x0)+f(x0)x-x0)+f(x0)(x-x0)2+ 增量较小时略去其高次幂项,则有 y-y0=f(x)-f(x0)=f(x0)(x-x0) k △y=k△x 比例系数,函数在X0点切线的斜率 两个变量的非线性函数 f(x1,x2),同样可在某工作点(x10,x20)附件用台劳级数展开为 y=f(xx2)=f(x10,x20)+/9(x0,x20)(x-x10)+9(x10,x20(x2-x20) ax 2 1(0.x20)(x1-x10)2+29(0.x20)(x-x0(x-x20 axle axlax2 a2f(x10,x20) 2 略去二级以上导数项,并令△y=y-f(x10,x20)

18 ( ) ( ) ( ) ( ) 1 2 t K U t K M t dt d t T m a c m m     ⑥ 式中 a m m e a m m R f C C R J T   电动机机电时间常数（s） a m m e m R f C C C K  1  m m e a Raf C C R K  2  电动机传递系数 如果电枢电阻Ra和电动机的转动惯量Jm都很小而忽略不计时 ⑥还可进一步简化为 C (t) U (t) e m  a ⑦ 电动机的转速Wm(t)与电枢电压Ua(t)成正比，于是 电动机可作为测速发电机使用。 系统最基本的数学模型是它的微分方程式。建立微分方程的步骤如下： ①确定系统的输入量和输出量 ②将系统划分为若干环节，从输入端开始，按信号传递的顺序，依据各变量所遵循的物理学 定律，列出各环节的线性化原始方程。 ③消去中间变量，写出仅包含输入、输出变量的微分方程式。 2.2.2 线性微分方程的求解 列出方程 求解方程 求解微分方程 初条 输入量 2.2.3 非线性元件微分方程的线性化 具有连续变化的非线性函数的线性化，可用切线法或小偏差法。在一个小范围内，将非线 性特性用一断直线来代替。（分段定常系统） ·一个变量的非线性函数 y=f(x) 在x0处连续可微，则可将它在该点附件用台劳级数展开         ' '' 2 ( 0)( 0) 2! 1 y f (x) f (x0) f (x0)(x x0) f x x x 增量较小时略去其高次幂项，则有 0 ( ) ( 0) ( 0)( 0) ' y  y  f x  f x  f x x  x Δy k Δx Δy=kΔx 比例系数，函数在x0点切线的斜率 ·两个变量的非线性函数 y=f(x1,x2)，同样可在某工作点（x10,x20）附件用台劳级数展开为                           ( 1 20) ] 2 ( 10, 20) ( 10)( 20) 1 2 ( 10, 20) ( 1 10) 2 1 ( 10, 20) [ 2! 1 ( 2 20)] 2 ( 10, 20) ( 1 10) 1 ( 10, 20) ( 1, 2) ( 10, 20) [ 2 2 2 2 2 2 x x x f x x x x x x x x f x x x x x f x x x x x f x x x x x f x x y f x x f x x 略去二级以上导数项，并令Δy＝y-f(x10,x20) Δx1=x-x10 Δx2=x-x20

(x10.x20)Ax1+9(x0.x20 Ax2=K1△xl+K2A axle 这种小偏差线性化方法对于控制系统大多数工作状态是可行的,平衡点附近,偏差一般不会 很大,都是“小偏差点”。 例2-4试把非线性方程z=xy在区域5≤x≤7、10≤y≤12上线性化。求用线性化方程来计算 当x=5,y=10时z值所产生的误差 解:由于研究的区域为5≤x≤7、10≤y≤12,故选择工作点x0=6,y0=11。于是z0=x0y0 =6×11=66 求在点x0=6,y0=11,z0=66附近非线性方程的线性化表达式。将非线性方程在点 0,y0,z0处展开成泰勒级数,并忽略其高阶项,则有 z-20=a(x-x0)+b(y-y0) 式中G/m=n=y0=11 b= exo 0=6 因此,线性化方程式为: z-66=11(x-6)+6(y-11) z=11x+6y-66 当x=5,y=10时,z的精确值为z=xy=5×10=50 由线性化方程求得的z值为z=11x+6y=55+60-66=49 因此,误差为50-49=1,表示成百分数502% 第4讲 数学工具一拉普拉斯变换与反变换 (1)拉氏变换定义设函数f(t)满足①t0时,f(t)分段连续 f(eldt <oo 则f(t)的拉氏变换存在,其表达式记作 F(s)=LIf(]= f(oe"dt (2)拉氏变换基本定理 线性定理La1f()+a22()=a1F1(s)+a2F2(s) 位移定理ef()=F(s+a) 延迟定理L(-可)=eF(s) lim f(o=lim SF(s) 终值定理 lim f(o=lim sF(S) 初值定理 df(o) L[,]=sF(s)-f(0) 微分定理 a1=sFs)-y(0)-

19 2 1 1 2 2 2 ( 10, 20) 1 10 ( 10, 20) x K x K x x f x x x x f x x y              这种小偏差线性化方法对于控制系统大多数工作状态是可行的，平衡点附近，偏差一般不会 很大，都是“小偏差点”。 例2-4 试把非线性方程 z=xy 在区域5≤x≤7、10≤y≤12上线性化。求用线性化方程来计算 当x=5，y=10时z值所产生的误差。 解： 由于研究的区域为5≤x≤7、10≤y≤12，故选择工作点x0=6，y0=11。于是z0=x0y0 =6×11=66. 求在点x0=6，y0=11，z0=66附近非线性方程的线性化表达式。将非线性方程在点 x0,y0,z0处展开成泰勒级数，并忽略其高阶项，则有 z-z0=a(x-x0)+b(y-y0) 式中 0 11 0 0        y x z a y y x x 0 6 0 0        x y z b y y x x 因此，线性化方程式为： z-66=11(x-6)+6(y-11) z=11x+6y-66 当x=5,y=10时，z的精确值为z=xy=5×10=50 由线性化方程求得的z值为z=11x+6y=55+60-66=49 因此，误差为50-49=1，表示成百分数 2% 50 1  第4讲 数学工具－拉普拉斯变换与反变换 ⑴ 拉氏变换定义 设函数f(t)满足 ①t0时，f(t)分段连续     f t e dt st 0 ( ) 则f(t)的拉氏变换存在，其表达式记作 F s L f t f t e dt st     0 ( ) [ ( )] ( ) ⑵拉氏变换基本定理 ·线性定理 [ ( ) ( )] ( ) ( ) 1 1 2 2 1 1 2 2 L a f t  a f t  a F s  a F s ·位移定理 L[e f (t)] F(s a) at    ·延迟定理 L[ f (t )] e F(s) s     ·终值定理 lim ( ) lim ( ) 0 f t sF s t  s   ·初值定理 lim ( ) lim ( ) 0 f t sF s t  s   ·微分定理 ] ( ) (0) ( ) [ sF s f dt df t L   ] ( ) (0) (0) ( ) [ 2 ' 2 2 s F s sf f dt d f t L   

叫(0=2(s)_f-(0) ·积分定理 S 可0m=(s)-o)-f°(o) (3)拉氏反变换 F(s)化成下列因式分解形式 F(s)=B6)2=k(s+Xs+2)(+=n) A( P1S+P2)…(S+P a.F(s)中具有不同的极点时,可展开为 F(s)=-+ S+ p, B(S) ar=L-(S+ Pr) --Pk b.F(s)含有共扼复数极点时,可展开为 s+a F(s) l3 (s+P1)(s+P2) B(s) [a, S+a2Jss-p- A(s) (s+P1)(s+P2) =-P C.F(s)含有多重极点时,可展开为 b = p1)(s+p1) (s+P1)(s+p+) (S+pu) b=/2() +P1)] d b(s) (s+P1)]} ds A(s) dB(s) j! dsA(s) (s+P1)]} b 1d- B(s) (r-1)!ds=A(s) (s+P1)]}= 其余各个极点的留数确定方法与上同

20 ·积分定理     s f s F s L f t dt ( ) (0) [ ( ) ] 1 s f s f s F s L f t dt ( ) (0) (0) [ ( ) ] 2 2 1 2       ⑶ 拉氏反变换 F(s)化成下列因式分解形式： ( )( ) ( ) ( )( ) ( ) ( ) ( ) ( ) 1 2 1 2 nm s p s p s p k s z s z s z A s B s F s           a. F(s)中具有不同的极点时，可展开为 n n s p a s p a s p a F s         2 2 1 1 ( ) k k k s p s p A s B s a    ( )] ( ) ( ) [ b. F(s)含有共扼复数极点时，可展开为 n n s p a s p a s p s p a s a F s           3 3 1 2 1 2 ( )( ) ( ) 1 1 ( )( )] ( ) ( ) [ ] [ 1 2 s p 1 2 s p s p s p A s B s a s a       c. F(s)含有多重极点时，可展开为 ( ) ( ) ( ) ( ) ( ) ( ) 1 1 1 1 1 1 1 1 n n r r r r r r s p a s p a s p b s p b s p b F s                   1 ( ) ] ( ) ( ) [ 1 s p r r s p A s B s b    1 1 1 ( ) ]} ( ) ( ) { [ s p r r s p A s B s ds d b     1 1 ( ) ]} ( ) ( ) { [ ! 1 s p r j j r j s p A s B s ds d j b     1 1 1 1 1 ( ) ]} ( ) ( ) { [ ( 1)! 1 s p r r r s p A s B s ds d r b       其余各个极点的留数确定方法与上同

F(s) f(t) 1(t) t s+a) r"(n=1,2,3…) t"e(n=1,2,3…) (S+a) (b-a) s+a)+o s+a e coS or (s+a)2+ s-(s+a) n(on√1-52n +20.s+ 第4讲补 补充电路的运算分析法一拉普拉斯变换分析法 (1)引言 拉普拉斯的由来即拉普拉斯法的作用与地位 定义 大量的实例求O()、1(t)、t、 sin ot、e等函数的拉氏变换 性质 ①线性②相似(尺度)③时延④位移(平域平移) ⑤微分⑥积分⑦初值⑧终值

21 第4讲补 补充 电路的运算分析法－拉普拉斯变换分析法 ⑴ 引言 拉普拉斯的由来即拉普拉斯法的作用与地位 一、 定义 大量的实例 求 (t)、1(t)、t、sint 、 at e  等函数的拉氏变换。 二、 性质 ①线性 ②相似（尺度） ③时延 ④位移（平域平移） ⑤微分 ⑥积分 ⑦初值 ⑧终值 F(s) f(t) 1  (t) 1(t) s 1 t 2 1 s at e  s  a 1 at te  2 ( ) 1 s  a sint 2 2   s  cost 2 2 s  s t (n  1,2,3) n 1 ! n s n (  1,2,3)  t e n n at 1 ( ) !   n s a n ( ) ( ) 1 at bt e e b a     ( )( ) 1 s  a s  b e t at sin  2 2 ( )   s  a  e t at cos  2 2 (  )   s a s a ( 1 ) 1 2 at at e a    ( ) 1 2 s s  a sin( 1 ) 1 2 2 e t n n t n         2 2 2 2 n n n s  s    

三、拉普拉斯逆变换 展开定理 F(s)= N(S) >m①D(s)=0有n个不相等的单根 Ak =[(s-SkF(S)Ises ②D(s)=0有实根的情况m阶重根。 Am d[(s-sl)F(s) (m-1)!ds ③D(s)=0有复根的情况 5.2符号法(拉普拉斯变换法) 、元件的运算模型 1、电阻U(s)=RI(s) u=R 2、电感U(s)=Ls(s)-L0.)u(s) 3、电 I(S)=CSU(S)-CU(O) 电路定律的运算形式 时域2()=0分2(s)=0 ∑l()=04>∑U(s)=0 线性定理 向量法中各种分析方法①阻抗的串并联 ②回路电流法 均可移植到运算法中 ③节点电压法 ④叠加定理和戴维南定理 运算法 应用拉氏变换求解线性动态电路的方法称为运算法 从原理上说,分析线性动态电路,应先列出电路的微分方程,然后应用拉氏变换将微 分方程转化为代数方程求解。但在工程分析中,往往直接利用运算电路图写出运算形 式的代数方程,从而简化求解过程 三步: ①计算t=0时电路中储能元件的初始值,画出换路后的运算电路 ②列出电路运算形式的方程,求出响应的象函数 ③应用展开定理求响应的原函数 例接通单位冲激电压源,试求Uc(t) 解:①画出运算电路 ②应用节点电压法 US) SCI= R R2 R U2(S)= C R,+ R R2 RR,C U(1) R,C 据电学基本定律可列出下列方程组:

22 三、 拉普拉斯逆变换 展开定理 ( ) ( ) ( ) D s N s F s  m>m ① D(s)=0 有n个不相等的单根 k k k s s A s s F s    [( ( )] ) ② D(s)＝0 有实根的情况 m阶重根。 1 1 ( 1)! [( 1) ( )]      m m m m ds d s s F s Am ③ D(s)=0 有复根的情况 5.2 符号法 （拉普拉斯变换法） 一、 元件的运算模型 1、电阻 U(s)=RI(s) u Ri R  2、电感 U(s) LsI(s)- Li(0 )  - U(s) dt di u l l  3、电容 I(s) CsU(s) - CU(0 )  - dt du i C c c  二、 电路定律的运算形式 时域 i(t)  0   I(s)  0 u(t)  0  U (s)  0 线性定理 向量法中各种分析方法 ① 阻抗的串并联 ② 回路电流法 均可移植到运算法中 ③ 节点电压法 ④ 叠加定理和戴维南定理 三、 运算法 应用拉氏变换求解线性动态电路的方法称为运算法。 从原理上说，分析线性动态电路，应先列出电路的微分方程，然后应用拉氏变换将微 分方程转化为代数方程求解。但在工程分析中，往往直接利用运算电路图写出运算形 式的代数方程，从而简化求解过程。 三步： ①计算t=0时电路中储能元件的初始值，画出换路后的运算电路。 ②列出电路运算形式的方程，求出响应的象函数。 ③应用展开定理求响应的原函数。 例 接通单位冲激电压源，试求Uc(t)。 解：①画出运算电路 ②应用节点电压法 1 2 1 1 ] 1 1 ( )[ R SC R R U S c    R R C R R S R C R R SC R U S c 1 2 1 2 1 1 2 1 1 1 1 1 ( )       t R R C R R c e R C U t 1 2 1 2 1 1 ( )    据电学基本定律可列出下列方程组：

Ri(o)+i(odt=U,(n) U0(1)=[(n)d→()= dU。(t) 消去中间变量i(t),得: dU0() Ro d+00()=U1() 在初始条件为零的情况下,进行拉氏变换得 U0() U,(O RCS+1 可见,输入输出象函数之比只与本电路结构参数R、C有关。它可用来表征电路本身的 特性,称为传递函数。由上述电路得到的这一概念可推广到一般情况。 第4讲 2.3控制系统的复域数学模型 2.3.1传递函数 是在用拉氏变换求解线性常微分方程的过程中引申出来的概念。 微分方程是在时域中描述系统动态性能的数学模型,在给定外作用和初使条件下,解 微分方程可以得到系统的输出响应。系统结构和参数变化是分析较麻烦。 用拉氏变化法求结微分方程时,可以得到控制系统在复数域的数学模型一传递函数。 定义:线性定常系统的传递函数,定义为零初始条件下,系统输出量的拉氏变换与输 入量的拉氏变换之比 设线性定常系统由下述n阶线性常微分方程描述 c(o+a al dr-c(()+.+a-c(0)+a, (0 d +b dim-7r()++bm-i dt r()+mr(o) 式中c()是系统输出量,r(1)是系统输入量,a1(=12.3…,m)和b(=12…m)是 与系统结构和参数有关的常系数 设r(t)和c(t)及其各阶系数在t=0是的值均为零,即零初始条件,则对上式中各项分别 求拉氏变换,并令c(s)=L[c(t)],R(s)=[r(t)],可得s的代数方程为 [aos"+a1s"+…+an-1S+an]C(s)=[bs"+b1s"+…+bn1S+am]R(s) 于是,由定义得系统传递函数为 C(s)_Dos M(s) G(s)= R(s) aos"+a S+a, N(s) 式中M(s)=bs"+bs"+…+bn1s+bm N(s)=aos"+a,s"+.+a-S+a X(s) U(s) 例5求例2机械系统与电路系统的传递函数X()和U,(S) (B1+B2)X2+(K1+K2)X=B1X+K1X (B1+B2)sXx(s)+(K1+K2)X2(s)=B1sX,(s)+K1X,(s)

23   ( )  ( ) 1 ( ) i t dt U t C Ri t i     dt dU t i t dt i t C C U t o ( ) ( ) ( ) 1 ( ) 0 消去中间变量i(t)，得： ( ) ( ) ( ) 0 0 U t U t dt dU t RC   i 在初始条件为零的情况下，进行拉氏变换得： 1 1 ( ) ( ) 0   U t RCS U t i 可见，输入输出象函数之比只与本电路结构参数R、C有关。它可用来表征电路本身的 特性，称为传递函数。由上述电路得到的这一概念可推广到一般情况。 第4讲 2.3 控制系统的复域数学模型 2.3.1 传递函数 是在用拉氏变换求解线性常微分方程的过程中引申出来的概念。 微分方程是在时域中描述系统动态性能的数学模型，在给定外作用和初使条件下，解 微分方程可以得到系统的输出响应。系统结构和参数变化是分析较麻烦。 用拉氏变化法求结微分方程时，可以得到控制系统在复数域的数学模型－传递函数。 定义：线性定常系统的传递函数，定义为零初始条件下，系统输出量的拉氏变换与输 入量的拉氏变换之比。 设线性定常系统由下述n阶线性常微分方程描述： ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 1 1 0 1 1 1 1 0 1 r t b r t dt d r t b dt d r t b dt d b c t a c t dt d c t a dt d c t a dt d a m m m m m m n n n n n n                  式中c(t)是系统输出量，r(t)是系统输入量， a (i 1,2,3, , n) i   和 b ( j 1,2, ,m) j   是 与系统结构和参数有关的常系数。 设r(t)和c(t)及其各阶系数在t=0是的值均为零，即零初始条件，则对上式中各项分别 求拉氏变换，并令c(s)＝L[c(t)]，R(s)=L[r(t)]，可得s的代数方程为： [ ] ( ) [ ] ( ) 1 1 1 0 1 1 0 1 a s a s a s a C s b s b s b s a R s m m m m n n n n                于是，由定义得系统传递函数为： ( ) ( ) ( ) ( ) ( ) 1 1 0 1 1 1 0 1 N s M s a s a s a s a b s b s b s b R s C s G s n n n n m m m m                  式中 m m m m M s  b s  b s    b s  b   1 1 0 1 ( ) n n n n N s  a s  a s    a s  a   1 1 0 1 ( ) 例5 求例2机械系统与电路系统的传递函数 ( ) ( ) X s X s rc 和 ( ) ( ) U s U s rc B1 B2 Xc K1 K2 Xc B1 Xc K1X r (  )  (  )     ( ) ( ) ( ) ( ) ( ) ( ) 1 2 1 2 1 1 B B sX s K K X s B sX s K X s  c   c  r  r