一阶非线性常微分方程解的存在性定理—Picard-Lindelof定理

上一节简单介绍了可求解的一阶常微分方程的解法,因为大部分非线性方程是不可解的,所以需要给出解的存在性的证明。本节主要介绍一阶非线性常微分方程Cauchy问题
$$
(E)\,\,\,\,\,\frac{dy}{dx}=f(x,y),\,\,\,\,\,y(x_{0})=y_{0}.
$$
解的存在性定理Picard-Lindelof定理(有的书上称它为Cauchy-Lipschitz定理). 对一阶常微分方程解的存在性理论作出重要贡献的数学家有Cauchy、Lipschitz、Picard、Lindelof、Peano等,其中Picard提出的Picard迭代法尤其值得关注。据传Picard证明Picard—Lindelof定理的原始论文足足有三四百页,后来数学家Banach把Picard的方法抽象出来证明了著名的Banach不动点定理。Banach不动点定理是分析学中最重要的定理之一,也是用的最多的定理之一,它在线性方程组求解迭代方法的收敛性、常微分方程的两点边值问题、隐函数定理、Lax-Milgram定理甚至代数方程解的存在性等问题中均有重要应用。许多微分方程(组)通过转化为等价的积分方程再利用不动点理论来证明解的存在性。本节也采用这一框架来探索方程(E)解的存在性。为此,首先利用Picard迭代给出Banach不动点定理的证明。

定理1 (Banach)  设$X$为Banach空间(即完备的赋范空间,完备的意思指所有的Cauchy列均收敛),$f:X\to X$为压缩映射,即存在常数$k, 0<k<1$,对任意$x,y\in X$有

$$
\|f(x)-f(y)\|\leq k\|x-y\|,
$$
则映射$f: X\to X$有且只有一个不动点$x\in X.$

证明: 任取$x_{0}\in X$,构造Picard迭代
$$
x_{n+1}=f(x_{n}),\,\,\,\,n\geq 0.
$$

$$
\|x_{n+1}-x_{n}\|=\|f(x_{n})-f_{x_{n-1}}\|\leq k\|x_{n}-x_{n-1}\|\leq\cdots\leq k^n\|x_{1}-x_{0}\|.
$$
设$m>n\geq 0$,由三角不等式和上式得
$$
\|x_{m}-x_{n}\|\leq \sum_{p=n}^{m-1}\|x_{p+1}-x_{p}\|\leq \frac{k^n}{1-k}\|x_{1}-x_{0}\|,
$$
当$m,n\to \infty$时,$\|x_{m}-x_{n}\|\to 0$, 故序列$\{x_{n}\}$为Cauchy列,由$X$的完备性知存在$x_{\infty}\in X$使得$\lim_{n\to\infty}x_{n}=x_{\infty}.$ $f:X\to X$满足Lipschitz条件,显然连续.故
$$
x_{\infty}=\lim_{n\to\infty}x_{n+1}=\lim_{n\to\infty}f(x_{n})=f(\lim_{n\to\infty}x_{n})=f(x_{\infty}).
$$
存在性得证。

误差估计:
$$
\|x_{n}-x_{\infty}\|=\lim_{m\to\infty}\|x_{n}-x_{m}\|\leq \frac{k^n}{1-k}\|x_{1}-x_{0}\|.
$$
若$\lim_{n\to\infty}x_{n}=x_{c}$,由上式知
$$
\|x_{c}-x_{\infty}\|=0.
$$
唯一性得证。证毕。

------

定理2 (Picard—Lindelof)  设初值问题
$$
(E)\,\,\,\,\,\frac{dy}{dx}=f(x,y),\,\,\,\,\,y(x_{0})=y_{0}.
$$
$f: Q_{a,b}\to \mathbb{R}$为连续函数,并且对第二个变量满足Lipschitz条件,即
$$
|f(x,y_{1})-f(x,y_{2})|\leq L |y_{1}-y|
$$
其中矩形区域$Q_{a,b}$ 为
$$
Q_{a,b}=\{(x,y):|x-x_{0}|\leq a, |y-y_{0}|\leq b\}.
$$
$\,a,\, b,\,L$均为常数。则(E)在局部区间$I=[x_{0}-h,x_{0}+h]$上有解,其中常数
$$
h=\min\left\{a,\frac{b}{M},\frac{1}{L}\right\},\,\,\,\,M=\max_{(x,y)\in \mathbb{R}}f(x,y)
$$
证明: 由微积分基本定理知,方程(E)等价于积分方程
$$
y(x)=y_{0}+\int_{x_{0}}^{x}f(s,y(s))ds.
$$
取区间
$$
I_{\varepsilon}=[x_{0}-\varepsilon,x_{0}+\varepsilon] \sub [x_{0}-a,x_{0}+a].
$$

$$
J_{\varepsilon}=[y_{0}-M\varepsilon,y_{0}+M\varepsilon]\sub [y_{0}-b,y_{0}+b]
$$

其中$\varepsilon$ 为待定常数,

定义映射
$$
F(y)=y_{0}+\int_{x_{0}}^{x}f(s,y(s))ds,
$$

$$
F: C(I_{\varepsilon};J_{\varepsilon})\to C(I_{\varepsilon};J_{\varepsilon}).
$$
事实上,
$$
|F(y)-y_{0}|=\left|\int_{x_{0}}^{x}f(s,y(s)ds\right|\leq M\varepsilon.
$$
取$C(I_{\varepsilon};J_{\varepsilon})$的上确界范数,压缩条件
$$
\|F(y_{1})-F(y_{2})\|=\sup_{x\in I_{\varepsilon}}\left|\int_{x_{0}}^{x}f(s,y_{1}(s))-f(s,y_{2}(s))ds\right|\leq L\varepsilon \|y_{1}-y_{2}\|
$$
故当
$$
\varepsilon <\frac{1}{L},\,\,\varepsilon<a,\,\,and,\,\,\varepsilon<\frac{b}{M}.
$$
时,由Banach不动点定理知存在唯一的$y\in C(I_{\varepsilon};J_{\varepsilon})$使得$F(y)=y$,即为原微分方程等价的积分方程的唯一解。

------

定理3 (改进的Picard-Lindelof)  设初值问题
$$
(E)\,\,\,\,\,\frac{dy}{dx}=f(x,y),\,\,\,\,\,y(x_{0})=y_{0}.
$$
$f: Q_{a,b}\to \mathbb{R}$为连续函数,并且对第二个变量满足Lipschitz条件,即
$$
|f(x,y_{1})-f(x,y_{2})|\leq L |y_{1}-y|
$$
其中矩形区域$Q_{a,b}$ 为
$$
Q_{a,b}=\{(x,y):|x-x_{0}|\leq a, |y-y_{0}|\leq b\}.
$$
$\,a,\, b,\,L$均为常数。则(E)在局部区间$I=[x_{0}-h,x_{0}+h]$上有解,其中常数
$$
h=\min\left\{a,\frac{b}{M}\right\},\,\,\,\,M=\max_{(x,y)\in \mathbb{R}}f(x,y)
$$
注:这个定理与上个定理的不同在于$h$的范围变大了一些。证明它的工具为以下推广的Banach不动点定理。

定理4 (推广的Banach不动点定理)  设$X$为Banach空间(即完备的赋范空间,完备的意思指所有的Cauchy列均收敛),$F^{n}:X\to X, n\geq 1$为压缩映射,则映射$F: X\to X$有且只有一个不动点$x\in X.$

定理4的证明 :

不妨设$n\geq 2$,由Banach不动点定理知存在唯一的$x_{\infty}\in X,F^{n}x_{\infty}=x_{\infty}$,又
$$
F^{n}(F(x))=F^{n+1}(x)=F(F^{n}x)=(F(x)).
$$
上式表明$F(x)$也是$F^{n}:X\to X$的一个不动点,由唯一性知$F(x)=x.$ 证毕.

定理3的证明:符号设定均与定理2的证明相同。


$$
\forall y,z\in C(I_{\varepsilon},J_{\varepsilon}),y_{1}=y,z_{1}=z.
$$

$$
y_{k+1}=y_{0}+\int_{x_{0}}^{x}f(s,y_{k}(s))ds,\,\,k\geq 1.
$$

$$
z_{k+1}=y_{0}+\int_{x_{0}}^{x}f(s,z_{k}(s))ds,\,\,k\geq 1.
$$

有估计式
$$
|y_{2}(x)-z_{2}(x)|=\left|\int_{x_{0}}^{x}f(s,y(s))-f(s,z(s))\right|\leq L\|y-z\|\cdot |x-x_{0}|,
$$
依次递推,
$$
|y_{3}-z_{3}|=\left|\int_{x_{0}}^{x}f(s,y_{2}(s))-f(s,z_{2}(s))\right|\leq\left|\int_{x_{0}}^{x}L|y_{2}-z_{2}|\right|
$$

$$
\leq L^{2}\|y-z\|\left|\int_{x_{0}}^{x}|s-x_{0}|ds\right|=\frac{L^2}{2!}\|y-z\|\cdot(x-x_{0})^{2}
$$

$$
\cdots
$$

$$
\|F^{n+1}y-F^{n+1}z\|\leq \frac{L^{n}\varepsilon^{n}}{n!}\|y-z\|
$$


$$
\lim_{n\to\infty}\frac{L^{n}\varepsilon^n}{n!}=0.
$$
也就是说存在$p\in \mathbb{N}^{+},s.t.\,\,F^{p}$为压缩映射,从而根据推广的Banach定理知映射$F: X\to X$有且只有一个不动点$x\in X.$ 这里对$\varepsilon$的限制为
$$
\varepsilon<a\,\,\,\,and,\,\,\,\,\varepsilon<\frac{b}{M}.
$$
证毕。

原文地址:https://www.cnblogs.com/zhangwenbiao/p/9563207.html

时间: 2024-08-02 09:19:44

一阶非线性常微分方程解的存在性定理—Picard-Lindelof定理的相关文章

poj 1659 Frogs&#39; Neighborhood Havel-Hakimi定理 可简单图定理

作者:jostree 转载请注明出处 http://www.cnblogs.com/jostree/p/4098136.html 给定一个非负整数序列$D=\{d_1,d_2,...d_n\}$,若存在一个无向图使得图中各点的度与此序列一一对应,则称此序列可图化.进一步,若图为简单图,则称此序列可简单图化. 可图化的判定为:$d_1+d_2+ \cdots +d_n=0(mod2)$.即把奇数度的点配对,剩下的变为自环.可简单图化的判定,即Havel-Hakimi定理: 我们把序列$D$变换为非

【POJ1548】Robots Dilworth定理(偏序集定理2)

题意: 有一些位置有垃圾,让机器人从左上角开始走,只能往右或者往下,问最少走多少次可以清理完所有垃圾. 题解: 一看就是网络流经典题,或者说是二分图-最小路径覆盖:但是现在毕竟是在做一些贪心,这道题用的是一种贪心相关定理,Dilworth定理. 这道题可以理解为部分两点之间有偏序(可走的关系),呃,可以视为当xa<=xb&&ya<=yb时有偏序,那么姑且认为反之则为反偏序,那么定义一条链为由n-1个偏序连接起来的n个点,那么答案就是"最长反链的大小. 比如题中的数据1

poj1061-青蛙的约会-(贝祖定理+扩展欧几里得定理+同余定理)

青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions:132162   Accepted: 29199 Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能

IT定理:摩尔定理,安迪-比尔定理,反摩尔定理

前两天在网上不经意间搜到了一本吴军的<浪潮之巅>,讲的是现代国际上计算机界的各大公司的兴衰沉浮,包括AT&T公司与IBM等等,把它当作IT历史书看,到现在已经看了一部分了.其中,我对计算机工业的生态链这一章很有兴趣,之中主要你讲了三个IT定理,分别是摩尔定理,安迪-比尔定理和反摩尔定理.还有一个很有意思的现象是,这些大公司基本上是在美国,而且美国的反垄断法对这些公司也是很有威胁的,目前我看到的好几家大公司的发展史上,都或多或少的被反垄断法所困扰. 摩尔定理 最早发现这个定理的是英特尔的

【bzoj3782】上学路线 dp+容斥原理+Lucas定理+中国剩余定理

题目描述 小C所在的城市的道路构成了一个方形网格,它的西南角为(0,0),东北角为(N,M).小C家住在西南角,学校在东北角.现在有T个路口进行施工,小C不能通过这些路口.小C喜欢走最短的路径到达目的地,因此他每天上学时都只会向东或北行走:而小C又喜欢走不同的路径,因此他问你按照他走最短路径的规则,他可以选择的不同的上学路线有多少条.由于答案可能很大,所以小C只需要让你求出路径数mod P的值. 输入 第一行,四个整数N.M.T.P. 接下来的T行,每行两个整数,表示施工的路口的坐标. 输出 一

赫尔布朗特定理开启自动定理证明的大门

1930年,年仅22岁的法国"小毛头"赫尔布朗特(Jacques Herbrand,1908-1931)在登山时遇难的前一年给我们留下了一条数理逻辑的基本定理:赫尔布朗特定理.这条定理有什么意义呢? 大家知道,在数理逻辑里面,引入量词符号"?"与"?",是很重要的,但是,对于数学自动鼎定理证明而言,量词符号"?"与"?"就是累赘.怎么办呢? 在所谓"一阶逻辑"里面,符号"?&q

cf451E Devu and Flowers 卢卡斯定理+容斥定理

题目:http://codeforces.com/problemset/problem/451/E 题意:有n个盒子(n<=20),每个盒子中有10^12个小球,现从每个盒子中取出若干球(可为0),求共取出s个小球(s<=10^14)的方案数. 组合数学问题,求C(n,m).但n,m过大时,可用卢卡斯定理. 卢卡斯定理:C(n,m) %p = C(n/p,m/p) * C(n%p,m%p) 从n个盒子中取出s个球的方案数,相当于插板,即 C(s+n-1,n-1).注意这是没有限制条件的情况.

[数分笔记]用Dedekind切割定理证明确界定理

1.定理内容 Dedekind切割定理:设是实数集的一个切割,则或者有最大数,或者有最小数. 确界定理:非空有上界的数集必有上确界,非空有下界的数集必有下确界. 2.证明过程 设非空数集有上界 记,即是上界的集合 令的补集为,即 从而形成实数集的一个切割 由Dedekind定理知,要么有最大数,要么有最小数 若有最大数,设是的最大数 由于,所以不是的上界 从而,s.t 那么,从而也不是的上界,故 与是的最大数矛盾,从而没有最大数 所以有最小数 即有最小上界,即上确界 #

BZOJ-1951-古代猪文-SDOI2010-费马小定理+欧拉函数+lucas定理+中国剩余定理

描述 =>G∑(ni),i|nmodP 分析 k=∑Cin,i|n(modP) G?(P)≡1(modP),?(p)=p?1 P′=P?1 =>GP′≡1(modP) Gk≡GkmodP′(modP) 如何求k? lucas定理 (nm)=(nmodP′mmodP′)?(n/P′m/P′) P'不是素数, lucas定理不适用. 所以把P'-1拆成2*3*4679*35617再用中国剩余定理来解. 一下子用这么多不熟悉的定理和方法感觉这个题好厉害. 终于知道当被模的数不是质数该怎么用中国剩余定