一阶微分方程解的存在定理

一阶微分方程解的存在定理

1 解的存在唯一性定理

1.1 Lipschitz条件

函数$f(x,y)$称为在$R$上关于$y$满足利普希茨条件，如果存在常数$L>0$，使得不等式
$$\left|f\left(x,y_1\right)-f\left(x,y_2\right)\right|⩽L\left|y_1-y_2\right|$$
对于所有的$(x,y_1),(x,y_2)\in R$都成立，$L$称为利普希茨常数.

1.2 解的存在唯一性定理

定理： 如果$f(x,y)$在矩形域$R$上连续且关于$y$满足利普希茨条件，则方程$\frac{dy}{dx}=f(x,y)$存在唯一的解$y=\phi (x)$，定义于区间$\left|x-x_0\right|\le h$上，且满足初值条件
$$\phi (x_0)=y_0$$
这里$h=\min (a,\frac{b}{M}),M={\max }_{(x,y)\in R}|f(x,y)|.$

2 逐步逼近法

考虑初值问题
$$\{\begin{array}{l}\frac{dy}{dx}=f(x,y)\\ y(x_0)=y_0\end{array}$$
我们可以构造一个逐步逼近函数列
$$\begin{array}{rl}{\phi }_0(x)& =y_0\\ {\phi }_1(x)& =y_0+{\int }_{x_0}^{x}f(x,{\phi }_0(x))dx\\ {\phi }_2(x)& =y_0+{\int }_{x_0}^{x}f(x,{\phi }_1(x))dx\\ ...& \\ {\phi }_n(x)& =y_0+{\int }_{x_0}^{x}f(x,{\phi }_{n-1}(x))dx\end{array}$$
如果${\phi }_n(x)={\phi }_{n+1}(x)$，则${\phi }_n(x)$就是方程的解。${\phi }_n(x)$称为初值问题的$n$次近似解.

${\phi }_n(x)$是定义在$|x-x_0|\le h$上，$h=min\left(a,\frac{b}{M}\right)，M={\max }_{(x,y)\in R}|f(x,y)|$

3 误差分析

第$n$次近似解${\phi }_n(x)$和真实解$\phi (x)$之间的误差为
$$\left|{\phi }_n(x)-\phi (x)\right|\le \frac{M{L}^n}{(n+1)!}{h}^{n+1}$$

在这里$h=min\left(a,\frac{b}{M}\right)，M={\max }_{(x,y)\in R}|f(x,y)|,L$为利普希茨常数

例题
方程$\frac{dy}{dx}=x^2+y^2$定义在矩形域$R:-1\le x\le 1,-1\le y\le 1$上，试利用存在唯一性定理确定经过点$(0,0)$的解的存在区间，并求在此区间上与真正解的误差不超过0.05的近似解的表达式.

解 这里$M={\max }_{(x,y)\in R}|f(x,y)|=2,h=\min (a,\frac{b}{M})=\frac{1}{2}$，在$R$上函数$f(x,y)=x^2+y^2$的利普希茨常数可取为$L=2$.因为对于$\forall x\in [-1,1],y_1,y_2\in [-1,1]$
$$|f(x,y_1)-f(x,y_2)|=|y_1^2-y_2^2|\le (|y_1|+|y_2|)|y_1-y_2|\le 2|y_1-y_2|$$
因为误差限定在不能超过0.05
$$|{\phi }_n(x)-\phi (x)|\le \frac{M{L}^n}{(n+1)!}{h}^{n+1}=\frac{M}{L}\frac{1}{(n+1)!}(Lh{)}^{n+1}=\frac{1}{(n+1)!}\le 0.05$$
所以可以取$n=3$.
我们构造逐步逼近序列
$$\begin{array}{rl}{\phi }_0(x)& =0,\\ {\phi }_1(x)& ={\int }_0^x[{\xi }^2+{\phi }_0^2(\xi )]d\xi =\frac{x^3}{3},\\ {\phi }_2(x)& ={\int }_0^x[{\xi }^2+{\phi }_1^2(\xi )]d\xi =\frac{x^3}{3}+\frac{x^7}{63},\\ {\phi }_3(x)& ={\int }_0^x[{\xi }^2+{\phi }_2^2(\xi )]d\xi \\ & ={\int }_0^x({\xi }^2+\frac{{\xi }^6}{9}+\frac{2{\xi }^{10}}{189}+\frac{{\xi }^{14}}{3969})d\xi \\ & =\frac{x^3}{3}+\frac{x^7}{63}+\frac{2x^{11}}{2079}+\frac{x^{15}}{59535}\end{array}$$
则${\phi }_3(x)$即为所求.