1欢迎各位来到《自动控制原理》课堂!2第2讲程向红控制系统的数学模型时域数学模型3第二章控制系统的数学模型2.1引言2.2时域数学模型2.3频域数学模型下讲2.4信号流图与梅逊公式下讲4数学模型的几种表示方式数学模型时域模型频域模型方框图和信号流图状态空间模型52.1引言描述系统或元件的动态特性的数学表达式叫做系统或元件的数学模型深入了解元件及系统的动态特性,准确建立它们的数学模型-称建模物理模型任何元件或系统实际上都是很复杂的,难以对它作出精确、全面的描述,必须进行简化或理想化。简化后的元件或系统为该元件或系统的物理模型。简化是有条件的,要根据问题的性质和求解的精确要求,来确定出合理的物理模型。6电子放大器看成理想的线性放大环节。通讯卫星看成质点。7建立控制系统数学模型的方法有:分析法-对系统各部分的运动机理进行分析,物理规律、化学规律。实验法-人为施加某种测试信号,记录基本输出响应。8分析法建立系统数学模型的几个步骤:建立物理模型。列写原始方程。利用适当的物理定律—如牛顿定律、基尔霍夫电流和电压定律、能量守恒定律等)选定系统的输入量、输出量及状态变量(仅在建立状态模型时要求),消去中间变量,建立适当的输入输出模型或状态空间模型。9实验法-基于系统辨识的建模方法黑匣子输入(已知)输出(已知)已知知识和辨识目的实验设计--选择实验条件模型阶次--适合于应用的适当的阶次参数估计--最小二乘法模型验证—将实际输出与模型的计算输出进行比较,系统模型需保证两个输出之间在选定意义上的接近102.2控制系统的时域数学模型2.2.1线性元件的微分方程U1R1R2U2C1C2图2-1RC组成的四端网络例2-1图2-1为由一RC组成的四端无源网络。试列写以U1(t)为输入量,U2(t)为输出量的网络微分方程。111111cUiRUdtiiCUc)(121112221ccUiRUdtiCUc222122cUU解:设回路电流i1、i2,根据克希霍夫定律,列写方程如下:①②③④⑤12dtdUCdtdUCic22222dtdUCdtdUCidtdUCicc2211211122211cUiRRU222222111)(UdtdUCRdtdUCdtdUCRc由④、⑤得由②导出将i1、i2代入①、③,则得1322222222211])([UdtdUCRdtdUCUiRdtdCR22222212112222211UdtdUCRdtdUCRdtdUCRdtUdCRCR1222221112222121)(UUdtdUCRCRCRdtUdCCRR这就是RC组成的四端网络的数学模型,是一个二阶线性微分方程。14试证明图2-2(a)、(b)所示的机、电系统是相似系统(即两系统具有相同的数学模型)。图2-2机电相似系统B1B2K1K2XrXc(a)机械系统R2C2R1C1UrUc(b)电气系统例2-215对电气网络(b),列写电路方程如下:c2c2cr1cr1XBXK)X-X(B)X-(XKrrcXKBXKKBB1121c21X)(X)(rUidtCiRidtCiR112211c22c11UCUCc11ciURUrc2c121UUU)i(RR解:对机械网络:输入为Xr,输出为Xc,根据力平衡,可列出其运动方程式③②④161)211(21)211(RCCRRUcCCUrirrccUCURUCCURR1121211)11()(利用②、③、④求出代入①将①两边微分得比较两个金色的公式,可得出如下机-电相似系统17力-电压相似机系统(a)和电系统(b)具有相同的数学模型,故这些物理系统为相似系统。(即电系统为即系统的等效网络)相似系统揭示了不同物理现象之间的相似关系。为我们利用简单易实现的系统(如电的系统)去研究机械系统......提供了方便。因为一般来说,电的或电子的系统更容易,通过试验进行研究。机械电阻R1电阻R2弹性系数K1弹性系数K2电气阻尼B1阻尼B21/C11/C218图2-3所示为电枢控制直流电动机的微分方程,要求取电枢电压Ua(t)(v)为输入量,电动机转速ωm(t)(rad/s)为输出量,列写微分方程。图中Ra(Ω)、La(H)分别是电枢电路的电阻和电感,Mc(N·M)是折合到电动机轴上的总负载转距。激磁磁通为常值。图2-3电枢控制直流电动机原理图SM负载+-+-LaRaEaWmJm,fmUaifia例2-319解:电枢控制直流电动机的工作实质是将输入的电能转换为机械能,也就是由输入的电枢电压Ua(t)在电枢回路中产生电枢电流ia(t),再由电流ia(t)与激磁磁通相互作用产生电磁转距Mm(t),从而拖动负载运动。因此,直流电动机的运动方程可由以下三部分组成。电枢回路电压平衡方程电磁转距方程电动机轴上的转距平衡方程20aaaaaaEtiRdttdiLtU)()()(①Ea是电枢反电势,它是当电枢旋转时产生的反电势,其大小与激磁磁通及转速成正比,方向与电枢电压Ua(t)相反,即Ea=Ceωm(t)②Ce-反电势系数(v/rad/s)•电枢回路电压平衡方程:21)()(tiCtMammmC)(tMm)()()()(tMtMtfdttdJcmmmmm③④-电动机转距系数(N·m/A)是电动机转距系数-是由电枢电流产生的电磁转距(N·m)•电动机轴上的转距平衡方程:fm-电动机和负载折合到电动机轴上的粘性摩擦系数(N·m/rad/s)Jm-转动惯量(电动机和负载折合到电动机轴上的)kg·m·•电磁转距方程:22)()()()()()()()(2tMRdttdMLtUCtCCfRdttdJRfLdttdJLcacaammemmammamamma)()()()(21tMKtUKtdttdTcammmemmamamCCfRJRTemmaCCRafRK2emmamCCfRCK1电动机机电时间常数(s)⑤⑥在工程应用中,由于电枢电路电感La较小,通常忽略不计,因而⑤可简化为③、④求出ia(t),代入①同时②亦代入①得:23如果电枢电阻Ra和电动机的转动惯量Jm都很小而忽略不计时⑥还可进一步简化为)(tm)()(tUtCame⑦)(tUa系统最基本的数学模型是它的微分方程式。建立微分方程的步骤如下:①确定系统的输入量和输出量②将系统划分为若干环节,从输入端开始,按信号传递的顺序,依据各变量所遵循的物理学定律,列出各环节的线性化原始方程。③消去中间变量,写出仅包含输入、输出变量的微分方程式。电动机的转速与电枢电压成正比,于是电动机可作为测速发电机使用。242.2.2线性微分方程的求解2.2.3非线性元件微分方程的线性化具有连续变化的非线性函数的线性化,可用切线法或小偏差法。在一个小范围内,将非线性特性用一断直线来代替。(分段定常系统)一个变量的非线性函数y=f(x)列出方程求解方程求解微分方程初条输入量在x0处连续可微,则可将它在该点附件用台劳级数展开2''')0)(0(!21)0)(0()0()(xxxfxxxfxfxfy25增量较小时略去其高次幂项,则有令Δy=kΔxk比例系数,函数在x0点切线的斜率两个变量的非线性函数y=f(x1,x2),同样可在某工作点(x10,x20)附近用台劳级数展开为0xxx)0)(0()0()(0'xxxfxfxfyy])201(2)20,10()20)(10(21)20,10(2)101(1)20,10([!21)]202(2)20,10()101(1)20,10([)20,10()2,1(222222xxxxxfxxxxxxxxfxxxxxfxxxxxfxxxxxfxxfxxfy0yyy26略去二级以上导数项,并令Δy=y-f(x10,x20)这种小偏差线性化方法对于控制系统大多数工作状态是可行的,平衡点附近,偏差一般不会很大,都是“小偏差点”。20221011xxxxxx2211222010112010),(),(xKxKxxxxfxxxxfy2710≤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处展开成泰勒级数,并忽略其高阶项,则有)()(000yybxxazz11000yxzayyxx6000xyzbyyxx因此,线性化方程式为: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,表示成百分数%2501例2-4试把非线性方程z=xy在区域5≤x≤7、28数学工具-拉普拉斯变换与反变换⑴拉氏变换定义设函数f(t)满足①t0时f(t)=0②t0时,f(t)分段连续则f(t)的拉氏变换存在,其表达式记作⑵拉氏变换基本定理线性定理位移定理延迟定理终值定理dtetfst0)(dtetftfLsFst0)()]([)()()()]()([22112211sFasFatfatfaL)()]([asFtfeLat)()]([sFetfLs)(lim)(lim0ssFtfst29数学工具-拉普拉斯变换与反变换续初值定理微分定理积分定理⑶拉氏反变换F(s)化成下列因式分解形式:a.F(s)中具有不同的极点时,可展开为)(lim)(lim0ssFtfst)0()(])([fssFdttdfL)0()0()(])(['222fsfsFsdttfdLsfssFdttfL)0()(])([1sfsfssFdttfL)0()0()(])([2212)())(()())(()()()(2121nmpspspszszszsksAsBsFnnpsapsapsasF2211)(kpskkpssAsBa)]()()([30b.F(s)含有共扼复数极点时,可展开为nnpsapsapspsasasF332121))(()(11)])(()()([][2121pspspspssAsBasac.F(s)含有多重极点时,可展开为)()()()()()(11111111nnrrrrrrpsapsapsbpsbpsbsF1])()()([1psrrpssAsBb111]})()()([{psrrpssAsBdsdb11]})()()([{!1psrjjjrpssAsBdsdjb1]})()()([{)!1(11111psrrrpssAsBdsdrb其余各极点的留数确定方法与上同。