伺服驱动与控制建模与Matlab仿真分析

1会计学伺服驱动与控制建模与伺服驱动与控制建模与Matlab仿真分析仿真分析1.1 模型分类1.2 建模基本方法l 机理模型l 统计模型 （1）频率响应 （2）系统辨识 (1) 物理模型：采用实物作为模型，可以按比例搭建； (2) 数学模型：以数学公式作为仿真对象； (3) 混合模型：既有物理模型也有数学模型。1.1 模型分类 1.2.1 机理模型法 （1）定义：采用由一般到特殊的推理演绎方法，对已知结构、参数的物理系统，运用相应的定律或定理，经过合理分析简化而建立起来的描述系统各物理量动、静态变化性能的数学模型。 l主要通过理论分析推导建立数学模型，常用到的理论知识包括：物质不灭定律、能量守恒定律、牛顿第二定律、基尔霍夫定律等。l提取主要因素、忽略次要因素。抓住对系统模型具有决定性影响的物理量及相互关系，舍弃次要。l注意系统的线性化。通过合理简化将非线性因素近似为线性系统。1.2 建模基本方法（2）机理建模实例 一阶倒立摆一阶倒立摆结构原理图1)摆杆绕其重心的转动方程为2)摆杆重心的水平运动可能描述为3)摆杆中心在垂直方向上的运动可描述为4)小车水平方向运动可描述为sincosyxJF lF l22(sin )xdFmxldt22( cos )ydFmgmldt202xd xFFmdt运动学与动力学分析建模：2222 222 2202 2202 2220()()sinsincos()()coscossincos()lgsincos()()JmlFlm Jmlm l gxJmlmmm lmlFm lmm mm lJmlmm 精确模型：22 22000200()()()lg()JmlFm l gxJ mmm mlmm mmlFJ mmm ml若只考虑在工作点附近 附近000010101cossin02双37高炮（2）机理建模实例 高炮炮塔随动系统微分方程（建模）拉氏变换（代数）稳定性分析系统性能分析PID控制器22+mmmddJfC idtdtdEKdtaadiULiREdt2+mmmJsfsC isEKaaUL siR iE0( )1( )pIDU sG sKKK se ss（ ）12( )mamammCG sUs L J sL f sKC（0101( )( )( )1( )( )( )fiVG s G s H sVG s G s H s闭环传函aR忽略( )ife sV V（2）机理建模实例 高炮炮塔随动系统系统框图开环传函1.2.2 统计模型法定义：采用由特殊到一般的逻辑、归纳方法，根据一定数量在系统运行过程中实测、观察的物理量数据，运用统计规律、系统辨识等理论合理估计出反映系统各物理量相互制约关系的数学模型。 由于其主要依据来自实测数据，又称为实验测定法。 常用于黑箱或灰箱问题，根据测得的系统输入、输出数据来建立实际系统的数学描述。 1.2 建模基本方法 激励信号的选择：l白噪声白噪声是指在较宽的频率范围内，各相同带宽频带所含的噪声能量相等的噪声。 l伪随机信号近似的白噪声l正弦扫频信号l多频正弦信号组合lchirp信号1( )sin2Nkkkkx taf t统计模型法 频率特性法 (1) 由已知数据绘制该系统开环频率响应bode图(2) 用20dB/dec及其倍数的折线逼近幅频特性，得到两个转折频率121/ ,2.85/rad srad s相应的惯性环节时间常数为12121110.35Ts Ts(3) 由低频幅频特性可知0( )0,1LK频率特性法建模实例 (4)由高频段相频特性知，该系统存在纯滞后环节，系统的开环传递函数应为以下形式122( )(1)(1)(1)(0.351)ssKeeG sT sT sss(5) 确定纯滞后时间111/, ()86rads 时11180()arctan1arctan 0.3586 112.85/, ()169rads 时再查图中22180()arctan 2.85arctan(0.352.85)2.85169 120.352s(6) 最终求得该系统的开环传递函数模型G(s)为0.35122( )(1)(1)(1)(0.351)ssKeeG sTsT sss频率特性法建模实例 系统辨识方法是现代控制理论中常用的方法，可根据系统的输入输出响应估计系统的动态模型。响应信号包括：频率响应、阶跃响应、伪随机响应、白噪声响应等。下图为系统辨识原理框图。统计模型法 系统辨识法 系统辨识的方法有许多种，这里主要讲述两种：Levy法和ARX法。（1）Levy法对连续系统的模型进行辨识Levy法源于Levy提出的对复数曲线进行拟合的一种方法 Complex-curve fitting J,IRE transactions on AC,1959.011( )1rrmmssG sss假设对象的传递函数为：, mr为待定系数通过实验可以获取对象的频率响应特性()iiG jPjQ其中i为采样点，, QP为采样点处的幅值与相位0112112()()1()rrmmjjBjBGjjjAjA问题：如何确定待定系数？从幅频特性的角度考虑所假定的对象传递函数，则有：/ 2 1/ 2 2211222100(1),(1)rriiiiiiiiBB/ 2 / 2 2211222100(1),(1)mmiiiiiiiiAA（1）Levy法对连续系统的模型进行辨识定义如下优化目标函数，以求最优的, mr21min()()()NJD jG jG j0,0,;0,1,;iiJJirimXB通过求解上述多元一次方程组，就可以得到最优的待定系数, mr（1）Levy法对连续系统的模型进行辨识Matlab中，提供了相应的工具，可以辨识连续系统的数学模型。num,den=invfreqs(H,w,r,m)H为通过实验数据获取的系统的频率响应数据，其格式为P-jQw为实验中所对应的频率点， 为角频率r为期望的对象模型分子的阶次，需要通过先验知识确定m为期望的对象模型分母的阶次，需要通过先验知识确定num为辨识出的对象模型的分子多项式系数den为辨识出的对象模型的分母多项式系数（2）Levy法对连续系统的模型进行辨识例：假设实际对象模型：3243272424( )10355024sssGsssss根据对象模型计算出频域响应数据（1）Levy法对连续系统的模型进行辨识调用B,A=invfreqs(H,w,3,4); G1=tf(B,A),可得：324324.0757.98131.411( )7.07410.4425.7331.25sssGsssss通过对比可知，其精度不高！（1）Levy法对连续系统的模型进行辨识（2）ARX法对离散系统的模型进行辨识一般离散系统的传递函数可写为：它对应的差分方程为：ARX模型辨识法就是通过辨识上述差分方程的系数而获取对象模型的。111211( )( )( )1mdmnzbb zb zy tG zzu ta za z112( )(1)()()(1)(1)nmy ta y ta y tnb u tdb u tdb y tdm上述差分方程的形式又被称作ARX（自回归遍历）模型。假设已知一组对象的输入输出数据根据ARX模型可得： (1),(2),(),(1),(2),()uuuu myyyy m111111(1)(0)(1)(1)()(1)(2)(1)(2)(2)(1)(2)()()()()()()nmnmnmya ya ynbudb y m dya ya ynbudb ym dy Ma y Ma y M nbu m db y M m dM(,)yy u 11,nmaabb目标：找出一组，使得最小（2）ARX法对离散系统的模型进行辨识定义如下的优化指标2min( )MiJi1TTy 0,0,;0,1,;iiJJinimab 上述系数的获取是通过使残差平方和最小而获得的， 又被称作最小二乘法。（2）ARX法对离散系统的模型进行辨识Matlab中，提供了相应的工具，可以辨识离散系统的数学模型。T=arx(y,u,m,n,d)y为对象的输出向量u为对象的输入向量m为期望的对象模型分子的阶次，需要通过先验知识确定n为期望的对象模型分母的阶次，需要通过先验知识确定d为期望的对象的纯滞后时间，需要通过先验知识确定T.A为辨识出的对象模型的分子多项式系数T.B为辨识出的对象模型的分母多项式系数（2）ARX法对离散系统的模型进行辨识例：实际对象模型324320.31240.57430.38790.0889( )3.2333.98692.22090.4723zzsG zzzzz根据对象模型得到的响应数据（2）ARX法对离散系统的模型进行辨识调用T=arx(y,u,4,4,1); G=tf(T),可得：1234112340.31260.58120.39580.091791()13.2554.0422.2680.486zzzzGzzzzz通过对比可知，其精度较高（2）ARX法对离散系统的模型进行辨识p m函数p 图形绘制功能p Simulink简介p S函数p Real-Time-Workshop & dSPACEp m文件与Simulink程序的混合编程p Matlab工具箱p 面向对象编程与GUI 利用上述系统，可以完成从系统建模、分析、离线仿真直到在线实验的全过程，可使开发人员全神贯注于控制方案的构思而不必再在琐碎的杂务上花费相当多的时间，从而可以大大缩短开发周期。目前已经成为进行快速控制原型验证和半实物仿真的首选实时平台。矩阵A所有特征根具有负的实部Ate有界系统的极点全部在左半平面判断准则 间接方法 Routh判据 Hurwitz判据 Lyapunov判据 matlab工具 eig（G） pzmap（G）o 稳定性分析三、基于Matlab的控制系统分析21 ,ncTB AB A BAB构造矩阵Matlab工具：Tc=ctrb(A,B) r=rank(Tc)如果矩阵Tc满秩，则系统完全可控，其秩为系统可控状态的个数。o 可控可观性分析21 ,nToTC CA CACA构造矩阵Matlab工具：To=obsv(A,C) r=rank(To)如果矩阵To满秩，则系统完全可观，其秩为系统可观状态的个数。三、基于Matlab的控制系统分析y,t=step(G) 自动选择时间向量对阶跃响应进行分析y,t=step(G，tf) 设置系统的终止时间，对阶跃响应进行分析y=step(G,t) 用户自己设置时间向量，对阶跃响应进行分析step(G1,-g,G2,-.b,G3,:r) 在同一图像窗口绘制多个系统的响应曲线，可设置它们的线型与颜色 o 线性系统阶跃响应分析三、基于Matlab的控制系统分析例1：已知系统传函28(1)(2)(3)( )(3.5)(5)(4)(1)(1)ssssG sessssjsjG=zpk(-1, -2, -3,-3.5, -5, -4, -1-j, -1+j,8,ioDelay,2)step(G)step(G,10)step(G, :r,10)三、基于Matlab的控制系统分析通过窗口图形对曲线进行编辑、数据分析三、基于Matlab的控制系统分析y,x=impulse (sys) 自动选择时间向量对脉冲响应进行分析y,x=impulse (sys,tf) 设置系统的终止时间，对脉冲响应进行分析y=impulse (sys,t) 用户自己设置时间向量，对脉冲响应进行分析step(sys1,-g,sys2,-.b,sys3,:r) 在同一图像窗口绘制多个系统的响应曲线，可设置它们的线型与颜色o 线性系统脉冲响应分析三、基于Matlab的控制系统分析424)(2sssG例：已知系统传函三、基于Matlab的控制系统分析 lism(sys,iu,t,x0)y,x= lism(sys,iu,t,x0)sys: 系统数学描述iu：输入向量t: 时间向量X0：初始状态o 线性系统任意输入响应分析三、基于Matlab的控制系统分析11)(ssG三、基于Matlab的控制系统分析 (1) 零极点图绘制 pzmap（sys） 极点用表示，零点用o表示 H=tf(2,5,1,1,2,3); pzmap(H)o 根轨迹分析三、基于Matlab的控制系统分析(2) r=rlocus(sys，k) r,k=rlocus(sys)1235 . 1)(23ssssGo 根轨迹分析三、基于Matlab的控制系统分析 (1) bode(num,den) mag,phase,w= bode(num,den) mag,phase= bode(num,den,w)o 频率特性响应分析1235 . 1)(23ssssGbode(1.5,1 3 2 1)三、基于Matlab的控制系统分析 margin（）函数，根据频率响应数据计算幅值和相角裕度margin（sys）margin（mag,phase,w）gm,pm,wcg,wcp=margin（mag,phase,w）mag: 幅值； phase：相角； w：角频率gm：幅值裕度；pm：相角裕度；wcg：相角交界频率；wcp：截止频率； o 频率特性响应分析三、基于Matlab的控制系统分析) 110)(15(10)(ssssGG=tf(10,conv(1,0,conv(5,1,10,1);bode(G,0.01,100);margin(G);gm,pm,wcg,wcp=margin(G)三、基于Matlab的控制系统分析双37高炮（2）机理建模实例 高炮炮塔随动系统 (1) 由已知数据绘制该系统开环频率响应bode图(2) 用20dB/dec及其倍数的折线逼近幅频特性，得到两个转折频率121/ ,2.85/rad srad s相应的惯性环节时间常数为12121110.35Ts Ts(3) 由低频幅频特性可知0( )0,1LK频率特性法建模实例 (4)由高频段相频特性知，该系统存在纯滞后环节，系统的开环传递函数应为以下形式122( )(1)(1)(1)(0.351)ssKeeG sT sT sss(5) 确定纯滞后时间111/, ()86rads 时11180()arctan1arctan 0.3586 112.85/, ()169rads 时再查图中22180()arctan 2.85arctan(0.352.85)2.85169 120.352s(6) 最终求得该系统的开环传递函数模型G(s)为0.35122( )(1)(1)(1)(0.351)ssKeeG sTsT sss频率特性法建模实例例1：已知系统传函28(1)(2)(3)( )(3.5)(5)(4)(1)(1)ssssG sessssjsjG=zpk(-1, -2, -3,-3.5, -5, -4, -1-j, -1+j,8,ioDelay,2)step(G)step(G,10)step(G, :r,10)三、基于Matlab的控制系统分析y,x=impulse (sys) 自动选择时间向量对脉冲响应进行分析y,x=impulse (sys,tf) 设置系统的终止时间，对脉冲响应进行分析y=impulse (sys,t) 用户自己设置时间向量，对脉冲响应进行分析step(sys1,-g,sys2,-.b,sys3,:r) 在同一图像窗口绘制多个系统的响应曲线，可设置它们的线型与颜色o 线性系统脉冲响应分析三、基于Matlab的控制系统分析11)(ssG三、基于Matlab的控制系统分析 (1) 零极点图绘制 pzmap（sys） 极点用表示，零点用o表示 H=tf(2,5,1,1,2,3); pzmap(H)o 根轨迹分析三、基于Matlab的控制系统分析