数学 旧 .com 迁移

常微分方程:相关资料 / 非线性微分方程组

从旧 .com 全量搬运的历史内容,来源路径:/math/课程/常微分方程/resources/非线性微分方程组/

迁移来源

非线性微分方程组

自治微分方程与非自治微分方程、动力系统 #

对于一般的 nn 阶非线性微分方程\begin{equation} y^{(n)} = G(t,y,y’,y”,\ldots,y^{(n-1)}) \end{equation} 可通过变换 xy,xy,,xn=y(n1)x_y,x_y',\ldots,x_n=y^{(n-1)} 化为如下一阶微分方程组

\wfenx1t=x2,,\wfenxn1t=xn,\wfenxnt=G(t,x1,x2,,xn).\wfen{x_1}{t}=x_2,\cdots,\wfen{x_{n-1}}{t}=x_n,\wfen{x_n} t=G(t,x_1,x_2,\ldots,x_n).

所以我们接下来研究更一般的一阶微分方程组

\begin{equation} \begin{cases} \wfen{x_1}t =f_1(t,x_1,x_2,\ldots,x_n),\

10pt] \wfen{x_2}t =f_2(t,x_1,x_2,\ldots,x_n),\

10pt] \vdots \\ \wfen{x_n}t =f_n(t,x_1,x_2,\ldots,x_n), \end{cases} \end{equation}

我们将上述方程组简记为向量形式 \begin{equation} \wfen{\bm x} t=\bm F(t,\bm x) \end{equation} 其中,

x=[x1x2xn],F(t,x)=[f1(t,x1,x2,,xn)f2(t,x1,x2,,xn)fn(t,x1,x2,,xn)].\bm x=\left[\begin{aligned} x_1\\x_2\\\vdots\\x_n \end{aligned}\right],\quad \bm F(t,\bm x)=\left[\begin{aligned} f_1(t,x_1,x_2,\ldots,x_n)\\f_2(t,x_1,x_2,\ldots,x_n)\\\vdots\\f_n(t,x_1,x_2,\ldots,x_n) \end{aligned}\right].

如果上述方程组有初值 \begin{equation} \bm x(t_0)=\bm x_0=(x_{01},x_{02},\ldots,x_{0n})^T. \end{equation} 则该初始值问题也存在类似定理 \ref{解的存在唯一性定理} 的解的存在唯一性定理.

definition

微分方程组 \eqref{一阶微分方程组} 在 n+1n+1 维空间 Rn+1={t,x1,x2,,xn}\mathbb{R}^{n+1}=\{t,x_1,x_2,\ldots,x_n\} 中确定了一个向量场, 而初始值问题 \eqref{一阶微分方程组向量},\eqref{微分方程组初始值问题} 的解 x(t,t0,x0)\bm x(t,t_0,\bm x_0) 就是向量场中的一条积分曲线. 当 \eqref{一阶微分方程组向量} 中函数 F\bm F 满足解的唯一存在性条件时, 向量场中任一点有且仅有一条积分曲线经过.

definition

如果把 tt 理解为时间参数, 只考虑 x1,x2,,xnx_1,x_2,\ldots,x_n 构成的空间 Rn\mathbb{R}^n, 我们将这个空间称为方程组 \eqref{一阶微分方程组向量} 的相空间, 积分曲线在相空间的投影曲线称为方程组的轨线.

definition

当方程组 \eqref{一阶微分方程组向量} 中的函数 F\bm F 显含 tt 时, 称该方程组为非自治微分方程组.

如果函数 F\bm F 中不显含 tt, 即 \begin{equation} \wfen{\bm x} t=\bm F(\bm x), \end{equation} 则称为自治微分方程组.

definition

系统 \eqref{一阶微分方程组向量} 的常数解 x=x\bm x=\bm x^* 称为系统的平衡点(\mydef[微分方程奇点]{奇点}或\mydef[微分方程驻点]{驻点})

definition

系统 \eqref{一阶微分方程组向量} 的解 x=x(t)\bm x=\bm x(t), 若存在常数 T>0T>0 满足 tR,s.t. x(t+T)=x(t)\forall t\in\mathbb{R},s.t.\ \bm x(t+T)=\bm x(t). 则称 x(t)\bm x(t) 是一个\mydef[微分方程周期解]{周期解}.

definition

设 \eqref{一阶微分方程组向量} 的右端函数 F(t,x)\bm F(t,\bm x) 对于 xGRn,tRx\in G\subset\mathbb{R}^n,t\in \mathbb{R} 连续, 关于 x\bm x 满足 \hyperref[lpxc条件]{\lpxc}且有一个解 x=Φ(t)\bm x=\bm\Phi(t).

现给定 t0Rt_0\in\mathbb{R} 并设 Φ0=Φ(t0)\bm\Phi_0=\bm\Phi(t_0). 如果对于任意的 ε>0\varepsilon>0, 存在至多依赖 ε,t0\varepsilon,t_0δ>0\delta>0, 使得对于 \eqref{一阶微分方程组向量} 的任意满足 x(t0)=x0x(t_0)=x_0 的解 x(t,t0,x0)x(t,t_0,\bm x_0), 只要 \begin{equation} \Vert \bm x_0-\bm \Phi_0\Vert<\delta \end{equation} 就有\begin{equation} \Vert \bm x(t,t_0,\bm x_0)-\bm\Phi(t)\Vert<\varepsilon,\quad \forall t\geqslant t_0 \end{equation} 就称解 x=Φ(t)x=\bm\Phi(t) 是 Lyapunov 意义下稳定的, 简称稳定的, 否则称不稳定的.

特别的, 如果 δ\delta 至多依赖 ε\varepsilon 而与 t0t_0 的取值无关, 那么称该解是 Lyapunov 一致稳定的.

definition

如果 \eqref{一阶微分方程组向量} 的解 x=Φ(t)\bm x=\bm\Phi(t) 是稳定的, 且存在一个常数 δ0>0\delta_0>0, 使得对一切满足 \begin{equation} \Vert \bm x_0 -\bm\Phi_0\Vert<\delta_0 \end{equation} 的解 x(t,t0,x0)\bm x(t,t_0,\bm x_0) 都有 \begin{equation} \lim\limits_{t\to+\infty}\Vert \bm x(t,t_0,\bm x_0)-\bm\Phi(t)\Vert=0. \end{equation} 则称该解是渐进稳定的.

definition

如果 \eqref{一阶微分方程组向量} 的解 x=Φ(t)\bm x=\bm\Phi(t) 是渐进稳定的且存在区域 D0D_0, 只要 x0D0\bm x_0\in D_0 就有

limt+x(t,t0,x0)Φ(t)=0.\lim\limits_{t\to+\infty}\Vert \bm x(t,t_0,\bm x_0)-\bm\Phi(t)\Vert=0.

则称 D0D_0 为该解的吸引域.

特别的, 如果某个解的吸引域是全空间, 则称此解是全局渐进稳定的.

info

在研究某个解的稳定性时, 总可以用变换 \begin{equation} \bm y(t)=\bm x(t)-\bm\Phi(t) \end{equation} 从而将 \eqref{一阶微分方程组向量} 化为 \begin{equation} \wfen{\bm y}{t}=\bm G(t,\bm y), \end{equation} 其中 G(t,y)=F(t,y+Φ)F(t,Φ)\bm G(t,\bm y)=\bm F(t,\bm y+\bm\Phi)-\bm F(t,\bm\Phi). 且显然有 G(t,0)=0\bm G(t,\bm 0)=0. 即该特解对应着新方程的零解, 所以我们接下来主要研究零解.

question

题目 #

试给出一阶微分方程

\wfenxt=a(t)x\wfen x t=a(t)x

的零解稳定或渐进稳定的充要条件. {{< admonition note “答案” false >}} 该方程的解为 x(t)=x(0)e0ta(s)dsx(t)=x(0)e^{\int_0^t a(s)\text{d} s}.

根据稳定性定义, 取 t0=0t_0=0, 则要求 x0<δ|x_0|<\deltax(t)<ε|x(t)|<\varepsilon

那么则需要 e0ta(s)dse^{\int_0^t a(s)\text{d} s} 有界.

渐近稳定, 又需满足 limt+x(t)=0\lim\limits_{t\to+\infty}\Vert \bm x(t)\Vert=0.

那么还需要条件 limt+e0ta(s)ds=0\lim\limits_{t\to+\infty}e^{\int_0^t a(s)\text{d} s}=0.

题目 #

给定极坐标系下的微分方程

\wfenθt=1,\wfenrt={r2sin1r,r>0,0,r=0.\wfen{\theta} t=1,\quad \wfen{r}{t}=\begin{cases} r^2\sin\dfrac 1 r, & r>0,\\ 0, & r=0. \end{cases}
  • (1) 证明平衡点 (0,0)(0,0) 是稳定的, 但不是渐近稳定的.
  • (2) 试作出 (0,0)(0,0) 邻域的相图.
  • (1) {{< admonition note “证明” false >}} 当 r(12kπ+π,12kπ)r\in(\dfrac{1}{2k\pi+\pi},\dfrac{1}{2k\pi}) 时, \wfenrt>0\wfen{r}{t}>0, 那么当 r0r_0 在这个区间内时, 根据 rr 的连续性且 r=12kπr=\dfrac 1{2k\pi}\wfenrt=0\wfen r t=0, 可推出 r(t)12kπr(t)\leqslant \dfrac{1}{2k\pi}.

类似的可以证明 r(t)12kπ+πr(t)\geqslant\dfrac{1}{2k\pi+\pi}.

那么只需取最大的 kk 满足 12kπ<ε\dfrac 1{2k\pi}<\sqrt{\varepsilon}, 那么当 r02<δ=12kπr_0^2<\delta=\dfrac 1{2k\pi} 时就有 r(t)2<εr(t)^2<\varepsilon. 进而说明 (0,0)(0,0) 是稳定的.

同时在上述过程中我们也说明了 r(t)r(t)t+t\to+\infty 时不是 00.

自治微分方程组解的性质 #

question

题目 #

note

解空间: x(t)=x0cost,y(t)=x0sintx(t)=x_0\cos t,y(t)=x_0\sin t.

轨线: x2+y2=x02x^2+y^2=x_0^2.

题目 #

note

{{< /admonition >}}

讨论

评论

正在加载评论...