线性时滞动力学问题的精细积分解法.pdfVIP

http :// 线性时滞动力学问题的精细积分解法 1 1 田象滔 任传波 1 山东理工大学交通与车辆工程学院 255049 E-mail: happytxt@ 摘要：具有时滞的动力系统广泛存在于各工程领域。本文给出了一种新的解决时滞系统动力 学的方法-精细积分法。把时滞项看作激励项进行处理，通过对时滞项积分上下限的讨论， 采用了线性插值和Romberg积分法。是同其它一些方法相比，用精细积分法解决时滞问题具 有过程简单，结果精确的优点。 关键词：精细积分，时滞，动力学，Romberg 积分法 1.引言 随着科学技术的发展，最近几十年来，自然科学和社会科学的领域提出了大量的时滞动 力学问题，如机械系统，电路系统，神经网络，光学系统，流行病学，生态系统等，社会科 学方面主要是各种经济现象的时滞描述，如市场供求关系问题，财富分布理论，运输调度问 题等，因此，时滞系统越来越受到广泛的重视。时滞现象在控制系统中尤其常见，甚至严格 说来是不可避免的。因此，时滞系统也就成为了优化控制和预测的一个非常重要的典型动力 学系统。现在已经有很多求解时滞动力学算法[1-8] ，但这些算法有些求解问题不太成功，有 些求解过程非常繁琐。本文采用了一种新的算法来解决时滞问题，这就是精细积分算法。 2. 时滞方程的精细积分算法 对于一般的线性时滞问题可以用下列式子来描述 y (t ) Ay (t) =+By (t −τ) +f (t) t ≥0 (1) y (t ) g (t) =−τ ≤t ≤0 (2) 对于(1)式可以将等式右边的后两部分看作是激励，则其解可以写为 t y (t) exp(At )y (0) =+ exp[A(t −s)][By (s −τ) +f (s)]ds ∫0 (3) t t exp(At)y (0) =+ exp[A(t −s)]By (s −τ)ds + exp[A(t −s)]f (s)ds ∫0 ∫0 取步长∆t tk +1 −tk ，对式(3)进行数值离散，有 t t k +1 k +1 y k +1 exp(A=∆t)y k + exp[A(tk +1 −s)]By (s −τ)ds + exp[A(tk +1 −s)]f (s)ds (4) ∫t ∫t k k 上式可进一步写成 (k +1)∆t y k +1 exp(A=∆t)y k +∫k∆t