控制)系统仿真6.ppt

(控制)系统仿真 (CONTOL) SYSTEM SIMULATION 第六章 simulink 五、S函数 System Function 扩展simulink的功能 状态方程的语句描述 可用m语言,C等编写 主要2个阶段:初始化(参数设置)\运行阶段(循环直到结束) Toolbox\simulink\blocks\ sfuntmpl1.m 一般连续系统,csfunc.m状态空间的连续系统,vsfunc变步长模块的模板. mdlinitializeSizes(flag=0) mdlderivatives(flag=1) mdlupdate(flag=2) mdloutputs(flag=3) mdlgettimeofnextvarhit(flag=4) [例6.10] 利用matlab仿真时变系统，A=[-0.09 -0.01；t2 0];B=[ 1 –7；0 -2];C=[ 0 0； 1 -5];D=[-3 0；1 0]; function [sys,x0,str,ts] = csfunc(t,x,u,flag) A=[0 0 0 0]; B=[ 1 -7 0 -2]; C=[ 0 0 1 -5]; D=[-3 0 1 0]; function sys=mdlDerivatives(t,x,u,A,B,C,D) a=t^2; A=[-0.09 -0.01 a 0]; sys = A*x + B*u; %%%over 发动机电子喷油控制 汽车发动机工作时，可简单分为启动、怠速（不踩油门）、正常行驶、急加速、减速、大负荷等工况，不同工况下，电喷系统喷油量不同（控制喷油器开启时间的脉冲信号宽度不同） 电子节气门 ECT为电子节气门模型 CONTROL为控制器 六 基于SIMULINK和NI的在环仿真 例，以ABS仿真为例： 2、硬件在环仿真 第七章 数值积分法 一、 连续系统仿真特点和方法 仿真是要将连续系统模型转化为离散时间模型，如差分方程、Z传递函数、离散状态方程等，最终用差分方程的形式来实现 由于转换方法、采样周期等的选择不同，同一系统所转换后的离散模型可能并不一样。如何保证离散模型和原连续模型等价或者相似，是数字仿真要解决的首要问题 三个基本要求。 ① 稳定性 原来稳定的系统，离散化后得到的仿真模型也应该是稳定的。 ② 准确性 准确性即仿真精度，其最基本的评价准则是绝对误差和相对误差要小于某个给定值。 ③ 快速性 快速性是指数字仿真的推进速度或计算所花的时间，包括时间的大小和确定性。 数值积分法 连续系统最基本的时域模型： 微分方程，高阶微分方程，微分方程组，状态方程----》微分方程（组）--》解： 解析法，低阶、特殊的微分方程 数值求解法，通用情况 数值积分法有很多具体方法，是一种适用范围很广的方法。由于其采用时域模型，容易理解，根据初值采用递推方式运算，可以方便地应用于非线性、时变系统，这是其它很多方法所不具备的优点； 缺点是计算量比较大，特别是为了提高精度时增加算法的复杂性以及减少步长对运算量影响较大，以及在循环中计算量大，有时不适合实时仿真。 二 欧拉法 根据泰勒级数展开，可得： 欧拉法 矩形欧拉法的几何意义 可以看出其误差 梯形欧拉法 问题：dy=f(y,t) 而y(n+1)未知 采用预估—校正方法，矩形法估计 y=y0;t=0;h=a; For i=0:100 t=t+h;dy1=f(y,t,u);y=y+h*dy1; %预估 dy2=f(y,t,u); y=y+h*(dy1+dy2)/2; %校正 end 通过例7。1，分析数值积分法的误差： （1）.原理误差或截断误差。根据泰勒展开式，欧拉法去掉了式中的高阶，从算法基本原理上来说就存在截断误差，且步长越大，误差也越大。 （2）.累积误差。递推的每一步都存在误差，如初值y(0)是精确的，但y(1)并不精确，因此在此基础上计算的导数也不精确，经过多次的循环，误差越来越大，严重情况下甚至出现不稳定的现象。 （3）.舍入误差。数字计算机能表达的精度或软件系统采用的数值类型使计算中的数据和结果不能足够精确地表达实际数据。如采用一般的浮点数类型，其精度较低。当循环次数多时，这些误差累计将达到一定程度。 例7 抛射物体的运动分析 　问题：子弹出膛速度为600m/s，计算各角度的水平射程（也可为任意目标角度的射角，设忽略g的变化，空气密度的变化、风的因素，子弹轴向总保持速度方向而不发生偏转翻滚，直径6.2mm，质量4克） 背景：在空中运动的物体，受到空气的阻力，在空气中如果速度低于2