自动控制原理控制系统的数据模型第二章控制系统的.ppt

第二章 控制系统的数学模型数学模型：描述控制系统输入变量、输出变量和内部变量之间关系的数学表达式，称为系统的数学模型描述控制系统动态特性的数学模型，称为动态模型在静态条件下(即变量的各阶导数为零)，描述变量之间关系的代数方程称为静态模型常用数学模型：常用解析形式的动态模型有微分方程、差分方程、传递函数；常用图形形式的动态模型有动态结构图、信号流图、频率特性建立数学模型的目的：用于分析控制系统的性能和设计满足要求性能的控制系统不同的分析设计方法常采用不同的数学模型，同一系统可用不同的数学模型表示1第二章 控制系统的数学模型¡2.1 传递函数¡2.2 闭环控制系统的动态结构图¡2.3 动态结构图的等效变换¡2.4 反馈控制系统的传递函数¡2.5 典型环节的传递函数¡2.6 信号流图与梅逊公式22.1 传递函数¡微分方程是在时域中描述系统动态性能的数学模型，在给定外作用和初始条件下，解微分方程可以得到系统的输出响应系统结构和参数变化时分析较麻烦¡传递函数是在用拉氏变换求解线性常微分方程的过程中引申出来的概念¡用拉氏变化法求解微分方程时，可以得到控制系统在复数域的数学模型－传递函数3一、传递函数　　传递函数:线性定常系统的传递函数，定义为在零初始条件下，系统输出量的拉氏变换与输入量的拉氏变换之比。

设线性定常系统由下述n阶线性常微分方程描述： 4式中c(t)是系统输出量，r(t)是系统输入量， 和是与系统结构和参数有关的常系数于是，由定义得系统传递函数为：　　设r(t)和c(t)及其各阶系数在t=0时的值均为零，即零初始条件，则对上式中各项分别求拉氏变换，并令R(s)＝L[c(t)]，R(s)=L[r(t)]，可得s的代数方程为：5G(s)R(s)C(s)(动态)方框图：M(s) – 分子多项式N(s) – 分母多项式，又称特征多项式，它决定着系统响应的基本特点和动态本质零极点形式、首一多项式形式、伊万思形式)(时间常数形式、尾一多项式形式、伯德形式)6zi – 传递函数的零点，即M(s)=0的根pj – 传递函数的极点，即特征方程N(s)=0的根，又称特征根传递函数的阶：特征多项式的阶次n即为传递函数的阶次，对应的系统为n阶系统静态放大系数(静态增益)：输出量与输入量静态值之比根轨迹增益：7　　　　传递函数由系统的结构和参数确定，与输入信号的形式与大小无关　　　　传递函数是复变量s的有理真分式函数，所有的系数均为实常数，且m≤n。

性质1：性质2：二、传递函数的特点　　　　如果传递函数已知，那么可以研究系统在各种输入信号作用下的输出响应性质3：传递函数与微分方程之间可以互相转换性质4：传递函数的拉氏变换是系统的单位脉冲响应性质5：传递函数的局限性：传递函数只适用于描述线性定常单输入、单输出系统，只直接反映系统在零状态下的动态特性8三、传递函数的求法 1、根据系统的微分方程求传递函数求系统传递函数的几个步骤：• 选定系统或元件的输入量、输出量• 列写原始方程利用适当的物理定律、化学定律或其它学科的公式和定律—如牛顿定律、基尔霍夫电流和电压定律、能量守恒定律等）• 在零初始条件下进行拉氏变换，消去中间变量，得到传递函数9　　解： 根据克希霍夫定律，列写出方程　　例2-1 试求图2-3所示RLC网络的传递函数输入量为ur(t)，输出量为uc (t) 在零初始条件下将以上两式进行拉氏变换后，得消去中间变量I(s)后得则求得传递函数uruc＋－＋－图2-3 RLC网络RLCi10F(t)y(t)mkf 图2-4 机械平移系统m – 质量块的质量f – 阻尼器的阻尼系数k – 弹簧的弹性系数F(t) – 外力y(t) –质量块的位移　　例2-2 设有一质量-弹簧-阻尼器的机械平移系统，如图2-4所示。

外力F(t) 为输入量， 质量块的位移y(t) 为输出量，试求系统的传递函数G(s)解：弹簧的恢复力阻尼器的阻力根据牛顿第二定律，得到传递函数11RLC网络机械平移系统令，则令，，则T – 时间常数， ζ – 阻尼比， Ｋ – (静态)放大系数或增益12(3) 模拟技术：有了相似系统的概念，可以利用对一种系统的研究来代替对另一种系统的研究，这就是所谓的模拟技术特别是用电子模拟装置模拟机械系统及其它物理系统讨论：(1) 两个完全不同的系统可能具有相同的传递函数2) 相似系统：物理量不同的两个系统具有相同形式的微分方程(数学模型)，这种系统称为相似系统而在微分方程中占据相同位置的物理量称为相似量机械系统Fmfkyv电路系统urLR1/Cqi表:相似量(uc=q/C)RLC网络机械平移系统132、用复阻抗概念求电路的传递函数例2-3无源网络 试求图2-5所示电路的传递函数 ur和uc分别是电路的输入量和输出量uruc＋－＋－图2-5 RLC网络R2LCR1u1＋－Ur(s)Uc(s)＋－＋－R2Ls1/CsR1U1(s)＋－传递函数：解14例2-4 有源网络 试求图2-7所示PI(比例-积分)调节器的传递函数。

ur和uc分别是有源网络的输入量和输出量传递函数解uruc＋－＋－图2-7 PI 调节器R2CR1＋－比例-积分调节器Kp=R2/R1 – 调节器的增益Ti=R2C – 调节器的时间常数153、非线性微分方程的线性化•现实系统中的元部件几乎都具有程度不同的非线性•解析法求解非线性微分方程非常困难在工程上，将非线性微分方程简化为线性微分方程具有很大的实际意义非线性微分方程的线性化：在工程上，将非线性微分方程在一定条件下转化为线性微分方程这一过程，称为非线性微分方程的线性化小偏差法：将非线性特性在(平衡)工作点附近(微小邻域)展开成泰勒级数，去掉高阶项而取其线性部分16例 非线性方程 y=f(x)，工作点(x0，y0)，求其小偏差线性化模型　　在工作点(x0，y0)处, 台劳级数展开式为0xyf(x)x0x0+Δxy0y0+Δy　　增量较小时略去其高次幂项而取一次近似式，则有 令 式中，Δy=KΔx则线性化方程为K为工作点处f(x)的一阶导数值，即该点的切线斜率　　忽略增量符号，可写成y=Kx (但应明确它是一个增量方程，y、x均为对平衡工作点的增量)几何意义：以工作点处的切线代替工作点邻域的曲线。

17例 流体运动系统hQ1Q2Q1 – 输入流量Q2 – 输出流量h – 液位高度质量守恒定律(AF – 贮槽的横截面积)节流阀的流量公式(Cv –节流阀的流量系数)则在工作点h0、 Q10、 Q20处(非线性项)那么18在工作点处的平衡关系线性化方程式可以省略Δ，简写成那么，传递函数19四、常见元部件的传递函数EKsθ1Ksθ2θ1θ2us+-图2-10 角度误差检测器θmax – 电位器的最大转角E – 施加在电位器两端的电压Ks=E/ θmax – 比例系数传递函数输入量：输出量：θe (t)us(t)1、角度误差检测器20u(a) 永磁式直流测速发电机＋－磁铁转子u＋－~(b) 交流测速发电机输出绕组激磁绕组ω(t) --被测物体的转速，输入量；u(t) – 发电机的输出电压，输出量　　直流测速发电机与交流测速发电机的输出电压均正比于电机转子的转速，即2、测速发电机：用来测量旋转物体的转速或传递函数图2-13 测速发电机213、机械转动系统：JMf图 2-14 机械转动系统M – 转矩，输入量；ω – 转速，输出量J – 旋转物体的转动惯量；f – 阻尼器的粘性摩擦系数；力矩平衡方程或传递函数224、齿轮系：用来减速或增大力矩。

电动机轴的力矩平衡方程J1、f1Mm图 2-15 齿轮系J2、f2M2M1负载电动机Z1Z2Mm – 电磁力矩，输入量；θ1 – 电动机轴的角位移，输出量Z1、Z2 – 主、从动轮的齿数；i – 齿轮传动的传动比；机械原理负载轴的力矩平衡方程23或传递函数消去M1、M2、θ2得(折算到电动机轴上的总等效转动惯量)其中(折算到电动机轴上的总等效粘性摩擦系数)245、电枢控制的他激直流电动机：广泛用作执行机构Eb、Kb – 电动机的反电势与反电势系数；θm – 电动机轴的角位移，输出量＋－＋If=const.－uaRaLa－＋Eaθm、ω、MLJ、f负载图 2-16 电枢控制的直流电动机Mm、Cm – 电动机的电磁力矩与力矩系数；ua – 电枢电压，输入量；25＋－＋If=const.－uaRaLa－＋Eaθm、ω、MLJ、f负载图 2-16 电枢控制的直流电动机电动机的反电势电动机的电磁力矩电枢回路电压平衡方程电动机轴上的力距平衡方程26消去Mm、Eb、ia得令La=0、ML=0得式中或传递函数(在位置控制系统中取转角)(在速度控制系统中取转速)(电动机的机电时间常数)(电动机的传递系数)27ω(t) – 转速，输出量。

ua(t) – 控制电压，输入量；6、两相伺服电动机：小功率交流执行机构图2-17 两相伺服电动机及其特性曲线控制绕组激磁绕组uaub转子0Mωua=75V50V25VM –输出转矩；Cω=dM/dω – 阻尼系数；Ms – 堵转转矩(在一定控制电压下)；线性化后的转矩—速度曲线方程堵转转矩与控制电压的关系CM – 比例系数，额定电压下的堵转转矩与额定控制电压的比值28转矩平衡方程消去Ms 、 M得J、f – 折算到电动机轴上的转动惯量和粘性摩擦系数式中(时间常数)(比例系数)或传递函数297、直流发电机Rf、Lf – 激磁绕组的电阻和电感ifω0图 2-18 直流发电机－＋ua－＋ufRfLfω0 – 原动机转速uf – 激磁电压，输入量ua – 电枢电压，输出量激磁回路电压平衡方程ω0恒定时，电枢电压(K为比例系数)传递函数(原动机转速恒定)(激磁恒定，与直流发电机相同)Tf= Lf / Rf – 电磁时间常数Kf= K/ Rf – 电磁时间常数30。