在数学 领域中, 时滞微分方程 , 或延时微分方程 (DDE) 是一类微分方程 , 其中未知函数的在确定时刻的导数由先前时刻函数所决定.
对于 x ( t ) ∈ R n {\displaystyle x(t)\in R^{n}} , 时滞微分方程方程的一般形式是:
d d t x ( t ) = f ( t , x ( t ) , x t ) , {\displaystyle {\frac {d}{dt}}x(t)=f(t,x(t),x_{t}),} 其中 x t = { x ( τ ) : τ ≤ t } {\displaystyle x_{t}=\{x(\tau ):\tau \leq t\}} 表示过去时间的解轨道. 在这个方程中, f {\displaystyle f} 是一个从 R × R n × C 1 {\displaystyle R\times R^{n}\times C^{1}} 到 R n {\displaystyle R^{n}\,} 的泛函算子.
例子 编辑
d d t x ( t ) = f ( t , x ( t ) , ∫ − ∞ 0 x ( t + τ ) d μ ( τ ) ) {\displaystyle {\frac {d}{dt}}x(t)=f\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )d\mu (\tau )\right)} d d t x ( t ) = f ( t , x ( t ) , x ( t − τ 1 ) , … , x ( t − τ n ) ) {\displaystyle {\frac {d}{dt}}x(t)=f(t,x(t),x(t-\tau _{1}),\ldots ,x(t-\tau _{n}))} for τ 1 > … > τ n ≥ 0 {\displaystyle \tau _{1}>\ldots >\tau _{n}\geq 0} . d d t x ( t ) = A 0 x ( t ) + A 1 x ( t − τ 1 ) + … + A m x ( t − τ m ) {\displaystyle {\frac {d}{dt}}x(t)=A_{0}x(t)+A_{1}x(t-\tau _{1})+\ldots +A_{m}x(t-\tau _{m})} 其中 A 0 , … , A m ∈ R n × n {\displaystyle A_{0},\ldots ,A_{m}\in R^{n\times n}} . d d t x ( t ) = a x ( t ) + b x ( λ t ) , {\displaystyle {\frac {d}{dt}}x(t)=ax(t)+bx(\lambda t),} 其中 a , b 且 λ 为常数 0 < λ < 1. 这一方程及其广义形式以电车上的集电弓 命名.
时滞微分方程求解 编辑
时滞微分方程通常用分步的方法求解. 例如考虑如下具有单一时滞的时滞微分方程
d d t x t = f ( x ( t ) , x ( t − τ ) ) {\displaystyle {\frac {d}{dt}}x_{t}=f(x(t),x(t-\tau ))} 及初始条件 ϕ : [ − τ , 0 ] → R n {\displaystyle \phi :[-\tau ,0]\rightarrow R^{n}} . 那么在区间 [ 0 , τ ] {\displaystyle [0,\tau ]} 上的解 ψ ( t ) {\displaystyle \psi (t)} 就是以下非齐次初值问题 的解
d d t ψ ( t ) = f ( ψ ( t ) , ϕ ( t − τ ) ) {\displaystyle {\frac {d}{dt}}\psi (t)=f(\psi (t),\phi (t-\tau ))} ,且 ψ ( 0 ) = ϕ ( 0 ) {\displaystyle \psi (0)=\phi (0)} . 这样就可以利用前面区间的解作为非齐次项一步步求得整个区间上的解. 在实际的计算中, 初值问题通常采用数值计算.
例子 编辑 假设 f ( x ( t ) , x ( t − τ ) ) = a x ( t − τ ) {\displaystyle f(x(t),x(t-\tau ))=ax(t-\tau )} 且 ϕ ( t ) = 1 {\displaystyle \phi (t)=1} . 那么初值问题可由积分求得,
x ( t ) = a ∫ s = 0 t ϕ ( t − τ ) d t + C {\displaystyle x(t)=a\int _{s=0}^{t}\phi (t-\tau )\,dt+C} ,即, x ( t ) = a t + 1 {\displaystyle x(t)=at+1} , 其中我们取 C = 1 {\displaystyle C=1} 以满足初值条件 x ( 0 ) = ϕ ( 0 ) {\displaystyle x(0)=\phi (0)} . 类似的对于区间 t ∈ [ τ , 2 τ ] {\displaystyle t\in [\tau ,2\tau ]} 我们积分并且使其满足初始条件可以求得 x ( t ) = a t 2 / 2 + t + D {\displaystyle x(t)=at^{2}/2+t+D} 其中 D = ( a − 1 ) τ + 1 − a τ 2 / 2 {\displaystyle D=(a-1)\tau +1-a\tau ^{2}/2} .
简化为常微分方程(ODE) 编辑
在某些情况下, 时滞微分方程等价于一个常微分方程组 (由常微分方程 组成的系统).
d d t x ( t ) = f ( t , x ( t ) , ∫ − ∞ 0 x ( t + τ ) e λ τ d τ ) . {\displaystyle {\frac {d}{dt}}x(t)=f\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )e^{\lambda \tau }d\tau \right).} 引入函数 y ( t ) = ∫ − ∞ 0 x ( t + τ ) e λ τ d τ {\displaystyle y(t)=\int _{-\infty }^{0}x(t+\tau )e^{\lambda \tau }d\tau } , 可得到一个常微分方程组 d d t x ( t ) = f ( t , x , y ) , d d t y ( t ) = x − λ y . {\displaystyle {\frac {d}{dt}}x(t)=f(t,x,y),\quad {\frac {d}{dt}}y(t)=x-\lambda y.} d d t x ( t ) = f ( t , x ( t ) , ∫ − ∞ 0 x ( t + τ ) cos ( α τ + β ) d τ ) {\displaystyle {\frac {d}{dt}}x(t)=f\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )\cos(\alpha \tau +\beta )d\tau \right)} 等价于 d d t x ( t ) = f ( t , x , y ) , d d t y ( t ) = cos ( β ) x + α z , d d t z ( t ) = sin ( β ) x − α y , {\displaystyle {\frac {d}{dt}}x(t)=f(t,x,y),\quad {\frac {d}{dt}}y(t)=\cos(\beta )x+\alpha z,\quad {\frac {d}{dt}}z(t)=\sin(\beta )x-\alpha y,} 其中 y = ∫ − ∞ 0 x ( t + τ ) cos ( α τ + β ) d τ , z = ∫ − ∞ 0 x ( t + τ ) sin ( α τ + β ) d τ . {\displaystyle y=\int _{-\infty }^{0}x(t+\tau )\cos(\alpha \tau +\beta )d\tau ,\quad z=\int _{-\infty }^{0}x(t+\tau )\sin(\alpha \tau +\beta )d\tau .} 特征方程 编辑
同常微分方程 (ODE)类似, 可以通过分析线性时滞微分方程的特征方程 [1] 来分析和研究解的性质.具有离散时滞的线性时滞微分方程
d d t x ( t ) = A 0 x ( t ) + A 1 x ( t − τ 1 ) + … + A m x ( t − τ m ) {\displaystyle {\frac {d}{dt}}x(t)=A_{0}x(t)+A_{1}x(t-\tau _{1})+\ldots +A_{m}x(t-\tau _{m})} 的特征方程是
d e t ( − λ I + A 0 + A 1 e − τ 1 λ + … + A m e − τ m λ ) = 0 {\displaystyle det(-\lambda I+A_{0}+A_{1}e^{-\tau _{1}\lambda }+\ldots +A_{m}e^{-\tau _{m}\lambda })=0} .特征方程的根 λ 被称为特征根或特征值, 解集通常被称为谱. 与常微分方程不同, 时滞微分方程的特征方程含有指数, 具有无限个特征值, 使得谱分析 变得很困难, 但是谱对于 DDE 的分析仍然具有一些很好的性质. 例如, 虽然具有无限个特征值, 但是只有有限个特征值位于复平面的右侧.
特征方程是一个非线性特征问题 , 有许多计算谱的数值方法[2] . 少数的特殊情况可以显式地求解特征方程. 例如, 时滞微分方程
d d t x ( t ) = − x ( t − 1 ) . {\displaystyle {\frac {d}{dt}}x(t)=-x(t-1).} 的特征方程是
− λ − e − λ = 0. {\displaystyle -\lambda -e^{-\lambda }=0.\,} 这个方程对于变量 λ 有无穷多个复数解. 复解可表示为
λ = W K ( − 1 ) {\displaystyle \lambda =W_{K}(-1)} ,其中 W K {\displaystyle W_{K}} 是朗伯W函数 的第 K 个分支.
注释 编辑 参考文献 编辑 外部链接 编辑