() 原连续模型y=f(y,u,1) (tn)≈0 仿真模型y=f(y,,tn) u(t vtn) 图21相似原理
u(t) h y(t) - + 图2.1 相似原理 原连续模型 y f y u t = ( , , ) 仿真模型 ˆ ( , , ) ˆ ˆ n y f y u t = ˆ( ) n y t ey (t n ) 0 ˆ( ) n u t
对仿真建模方法三个基本要求 冷(1)稳定性:若原连续系统是稳定的,则离散化 后得到的仿真模型也应是稳定的。 (2)准确性:有不同的准确性评价准则,最基本 的准则是: 绝对误差准则:、()=(n)-y()≤6 相对误差准则:k,(,)=5n)156 j(tn)-y(tn )I 令其中δ规定精度的误差量
对仿真建模方法三个基本要求 ❖ (1)稳定性:若原连续系统是稳定的,则离散化 后得到的仿真模型也应是稳定的。 ❖ (2)准确性:有不同的准确性评价准则,最基本 的准则是: ❖ 绝对误差准则: ❖ 相对误差准则: ❖ 其中 规定精度的误差量。 ey (t n ) = y ˆ(t n ) − y(t n ) − = ˆ( ) ˆ( ) ( ) ( ) n n n y n y t y t y t e t
对仿真建模方法三个基本要求(续) 令3)快速性:若第n步计算对应的系统时间间 隔为h=m-1计算机由计算需要的时间为Tm, 若Tm=hn称为实时仿真,Tn<hn称为超实 时仿真7m>hn称为亚实时仿真,对应于离 线仿真
对仿真建模方法三个基本要求(续) ❖ 3)快速性:若第n步计算对应的系统时间间 隔为 计算机由计算需要的时间为Tn, 若 Tn=hn 称为实时仿真,Tnhn称为超实 时仿真 Tnhn 称为亚实时仿真,对应于离 线仿真 , n n 1 n h = t − t +
数值积分算法 对y=f(y,,1),已知系统变量y的初始条 件y(x)=y。,要求y随时间变化的过程一一初值问 题 今计算过程:由初始点y)=y的f(t0,y) y() 0+O f(t, y)dt 欧拉法 yi=yt,@yo+ Dt?f(to, yo) 冷对任意时刻tn+1 ym1=y(uuD@y,+(tI-t,)?f(t,, ym,) 截断误差正比于h f(t, y f(t
数值积分算法 ❖ 对 ,已知系统变量y的初始条 件 ,要求y随时间变化的过程――初值问 题 ❖ 计算过程:由初始点 的 ❖ 欧拉法 ❖ 对任意时刻tn+1 ❖ 截断误差正比于 y f y u t = ( , , ) 0 0 y t y ( ) = 0 0 y t y ( ) = ( ) 0 0 f t ,y 0 0 ( ) ( , ) t t y t y f t y dt = + ò 1 1 0 0 0 y y t y t f t y = @ + D ? ( ) ( ) , 1 1 1 ( ) ( ) ( ) n n n n n n n y y t y t t f t y + + + = @ + - ? , h 2 f(t,y) f(t0 ,yo ) t t t0 t1
数值积分算法续) 令梯形法:m=y(n)+21(nyn)+f(m,ym) 是隐函数形式。预报一一欧拉法估计初值, 校正一一用梯形法 校正 校正公式m”yn+2(n,y,)+f(m,ym 今预报公式m≡十h,f(n,yn) 反复迭代,直到满足pm-ymls 经典的数值积分法可分为两类: 单步法与多步法
数值积分算法(续) ❖ 梯形法: ❖ 是隐函数形式。预报-—欧拉法估计初值, 校正-—用梯形法校正: ❖ 校正公式 ❖ 预报公式 ❖ 反复迭代,直到满足 ❖ 经典的数值积分法可分为两类: ❖ 单步法与多步法 [ ( , ) ( , )] 2 1 ( ) n+1 = n+1 n + n n + n+1 n+1 y y t y h f t y f t y [ ( , ) ( , )] 2 1 ( ) 1 1 ( 1) 1 i n n n n n i n y y h f t y f t y + + + + + + ( , ) ( ) 1 n n n i n y y + h f t y + − + + + i n i n y y 1 1 1