常微分方程 $6 一阶微分方程解的存在唯一性
§6 解的存在唯一性
C1 存在唯一性定理
dydx=f(x,y)\frac{\mathrm{d}y}{\mathrm{d}x}=f(x,y)dxdy=f(x,y),fff是矩形域R:∣x−x0∣≤a,∣y−y0∣≤bR:|x-x_0|\le a,|y-y_0|\le bR:∣x−x0∣≤a,∣y−y0∣≤b上连续函数 式6.6.1
1)Lipschitz条件
:若∣f(x,y1)−f(x,y2)∣<L∣y1−y2∣|f(x,y_1)-f(x,y_2)|<L|y_1-y_2|∣f(x,y1)−f(x,y2)∣<L∣y1−y2∣,则称fff关于yyy满足Lipschitz条件,LLL称Lipschitz常数
- 存在连续偏导数∂f∂y⟹\frac{\partial f}{\partial y}\implies∂y∂f⟹fff关于yyy满足Lipschitz条件
2)存在唯一性定理:若f(x,y)f(x,y)f(x,y)在矩形域RRR上连续且关于yyy满足Lipschitz条件,则式6.1.1才能在唯一解,定义于∣x−x0∣≤min(a,bmaxf)|x-x_0|\le \min(a,\frac{b}{\max f})∣x−x0∣≤min(a,maxfb)上
在证明过程中,积分曲线被限制在如下区域中
若f(x,y)=P(x)y+Q(y)f(x,y) = P(x)y + Q(y)f(x,y)=P(x)y+Q(y),则满足上述条件,且此时积分曲线不再被限制
推广:若点(x0,y0,y0′)(x_0,y_0,y_0')(x0,y0,y0′)的某一邻域中:
- F(x,y,y′)F(x,y,y')F(x,y,y′)对所有变元(x,y,y′)(x,y,y')(x,y,y′)连续,且存在连续偏导数
- F(x0,y0,y0′)=0F(x_0,y_0,y_0')=0F(x0,y0,y0′)=0
- ∂F∂y′≠0\frac{\partial F}{\partial y'}\neq 0∂y′∂F=0
则隐式微分方程F(x,y,y′)=0F(x,y,y')=0F(x,y,y′)=0有唯一解y=y(x),∣x−x0∣<ϵy=y(x),|x-x_0|<\epsilony=y(x),∣x−x0∣<ϵ,满足初值条件y(x0)=y0,y′(x0)=y0′y(x_0)=y_0,y'(x_0)=y'_0y(x0)=y0,y′(x0)=y0′
3)逐步逼近法:求解式6.1.1,令ϕn+1(x)=y0+∫x0xf(x,ϕn(x))dx\phi_{n+1}(x)=y_0+\int_{x_0}^xf(x,\phi_n(x))\mathrm{d}xϕn+1(x)=y0+∫x0xf(x,ϕn(x))dx,则limn→∞ϕn(x)=ϕ(x)\lim\limits_{n\to\infin}\phi_n(x)=\phi(x)n→∞limϕn(x)=ϕ(x)为解
证明:分五个命题证明
- 命题1:y=yo+∫x0xf(x,y)dxy=y_o+\int_{x_0}^xf(x,y)\mathrm{d}xy=yo+∫x0xf(x,y)dx与dydx=f(x,y)\frac{\mathrm{d}y}{\mathrm{d}x}=f(x,y)dxdy=f(x,y)同解
- 命题2:ϕn(x)\phi_n(x)ϕn(x)在[x0,x0+h][x_0,x_0+h][x0,x0+h]上有定义,连续且∣ϕn(x)−b∣≤b|\phi_n(x)-b|\le b∣ϕn(x)−b∣≤b
- 命题3:{ϕn(x)}\{\phi_n(x)\}{ϕn(x)}在[x0,x0+h][x_0,x_0+h][x0,x0+h]上一致收敛,limn→∞ϕn(x)=ϕ(x)\lim\limits_{n\to\infin}\phi_n(x) = \phi(x)n→∞limϕn(x)=ϕ(x)
- 命题4:ϕ(x)\phi(x)ϕ(x)是积分方程y=yo+∫x0xf(x,y)dxy=y_o+\int_{x_0}^xf(x,y)\mathrm{d}xy=yo+∫x0xf(x,y)dx的连续解
- 命题5:解唯一
- 误差公式:ϕn(x)−ϕ(x)≤MLn(n+1)!hn+1\phi_n(x)-\phi(x)\le \frac{ML^n}{(n+1)!}h^{n+1}ϕn(x)−ϕ(x)≤(n+1)!MLnhn+1
C2 延拓
1)局部Lipschitz条件
:对于区域GGG内每一点都存在一个含于GGG的满足Lipschitz条件的矩形域
2)解的延拓定理:若式6.1.1中的fff在GGG上满足局部Lipschitz条件,则解的存在性可以一致延拓到点(x,ϕ(x))(x,\phi(x))(x,ϕ(x))接近GGG的边界为止。这里的边界也可以是开边界,或者无穷
3)解对初值的连续性定理:若式6.1.1中fff满足局部Lipschitz条件,则y=ϕ(x,x0,y0)y=\phi(x,x_0,y_0)y=ϕ(x,x0,y0)关于初值也是连续的
4)解对初值的可微性定理:若式6.1.1中∂f∂y\frac{\partial f}{\partial y}∂y∂f在GGG上连续则y=ϕ(x,x0,y0)y=\phi(x,x_0,y_0)y=ϕ(x,x0,y0)关于x0,y0x_0,y_0x0,y0连续可微,且
∂ϕ∂x0=−f(x0,y0)e∫x0x∂f(x,ϕ)∂ydx∂ϕ∂y0=e∫x0x∂f(x,ϕ)∂ydx\frac{\partial \phi}{\partial x_0}=-f(x_0,y_0)e^{\int_{x_0}^x\frac{\partial f(x,\phi)}{\partial y}\mathrm{dx}}\\\frac{\partial \phi}{\partial y_0}=e^{\int_{x_0}^x\frac{\partial f(x,\phi)}{\partial y}\mathrm{dx}} ∂x0∂ϕ=−f(x0,y0)e∫x0x∂y∂f(x,ϕ)dx∂y0∂ϕ=e∫x0x∂y∂f(x,ϕ)dx
C3 奇解
1)包络
:特殊积分曲线,不属于曲线族,但线上每一点都与曲线族中一条积分曲线相切,对应的解称奇解
- 单参数曲线族Φ(x,y,c)=0\Phi(x,y,c)=0Φ(x,y,c)=0的包络为{Φ(x,y,c)=0Φc′(x,y,c)=0\begin{cases} \Phi(x,y,c)=0\\\Phi'_c(x,y,c)=0\end{cases}{Φ(x,y,c)=0Φc′(x,y,c)=0消去ccc的某一解,称
c-判别曲线
- 方程F(x,y,y′)=0F(x,y,y')=0F(x,y,y′)=0的奇解为{F(x,y,p)=0Fp′(x,y,p)=0\begin{cases}F(x,y,p)=0\\F'_p(x,y,p)=0\end{cases}{F(x,y,p)=0Fp′(x,y,p)=0消去ppp的某一解,称
p-判别曲线
- 一般的曲线族不一定有包络,如同心圆族,平行线族
Clairaut微分方程:y=xp+f(p),p=dydxy=xp+f(p),p=\frac{\mathrm{d}y}{\mathrm{d}x}y=xp+f(p),p=dxdy,通解为直线族y=cx+f(c)y=cx+f(c)y=cx+f(c),包络为{x+f′(p)=0y=xp+f(p)\begin{cases}x+f'(p)=0\\y=xp+f(p)\end{cases}{x+f′(p)=0y=xp+f(p)消去p
常微分方程 $6 一阶微分方程解的存在唯一性相关推荐
- 向前欧拉公式 matlab_阻尼谐振子运动的欧拉数值解法
这篇文章讲一点数值求解阻尼振子运动的微分方程相关的东西,总体来说是一个比较trivial的东西(划水本性暴露),为啥要研究呢,是因为这算是我们计算物理课程的一个作业,因为我以前自己写过四阶龙格库塔法求 ...
- 建模软件matlab,方程建模与MATLAB软件
第1章MATLAB的基本使用方法 1.1MATLAB概述 1.1.1MATLAB发展史 1.1.2MATLAB帮助 1.2MATLAB基础知识 1.2.1MATLAB变量 1.2.2MATLAB数据类 ...
- 高等数学学习笔记——第一百讲——微分方程稳定性初步
一.问题引入--如何预测种群数量的变化规律?如何控制卫星运行的稳定性? 二.稳定性的概念 1. 一阶微分方程解的变化 2. 一阶微分方程组解的变化 3. 一阶方程的平衡点及稳定性(平衡点.奇点).稳定 ...
- 一阶非线性常微分方程解的存在性定理—Picard-Lindelof定理
上一节简单介绍了可求解的一阶常微分方程的解法,因为大部分非线性方程是不可解的,所以需要给出解的存在性的证明.本节主要介绍一阶非线性常微分方程Cauchy问题 $$ (E)\,\,\,\,\,\frac ...
- python求解一阶常微分方程
一.用python求解一阶常微分方程 1.求解微分方程需要用到scipy库,pycharm中安装即可,同时需要导入numpy库和matplotlib两个库 2.使用scipy.integrate.od ...
- 使用C++,用四阶Runge-Kutta的方法来求解一阶常微分方程
#include <iostream> using namespace std;/* dy/dx = y - 2x/y, 0< x <= 1 步长h = 0.2 */ cons ...
- 一阶常微分方程(一)| 存在与唯一性定理 + 变量分离 + 齐次方程
一.一阶常微分方程解的的存在与唯一性定理 导数已解出的一阶常微分方程可以表示为如下的一般形式: {dydx=f(x,y)y(x0)=y0(1)\begin{cases} \frac{dy}{dx}=f ...
- 【数理知识】神仙文章(回忆大学所学)常微分方程
常微分方程分类及求解 1. 基本概念 1.1. 微分方程 (Differential Equation, D.E.) 1.2. 常微分方程 1.3. 解 1.4. 通解 1.5. 通积分 1.5. 阶 ...
- 二阶齐次线性微分方程的通解公式_数学考研辅导(数二)063|高数常微分方程6
5 常微分方程 题型五:二阶常系数线性微分方程 2011年全国硕士研究生招生考试数学(二)试题 一.选择题:1~8小题,每小题4分,共32分,下列每小题给出的四个选项中,只有一项符合题目要求的,请将所 ...
- 高数_第5章常微分方程_二阶线性微分方程解的结构
二阶线性微分方程的一般形式是 y'' + p(x)y' + q(x)y = f(x) -------- (1) 这里所谓的线性是指未知函数 y 及其 一阶导y', ...
最新文章
- MyBatis参数传入集合之foreach动态sql
- docker安装influxdb、grafana及展示数据到grafana
- 关于java几种输出的区别
- java绘画imo,搞死帮忙看下这个:需要为class interface或enum imort java.ut
- maven打包失败:自定义项目工具类打包给其他微服务使用
- centos7 docker端口_docker相关知识
- Ubuntu下利用python实现自动打卡系统全过程
- google bert deeping watch
- 计算机控屏教程,电脑怎么控制手机操作
- 【国产RX解码】XS9922A 四通道多合一同轴高清解码芯片 功能对标TP9930
- 银行业务模拟系统(C/C++实现内含详细注释)
- Easy Mock全解及使用
- BuildAdmin后台管理系统配置QQ邮箱
- 银行测试的优势在哪里?
- 微信小程序合集源码I(机器人聊天+仿别踩白块儿小游戏+仿步步高电子词典+仿知乎+日记+汉字拼音+转盘抽奖)
- r语言中怎样查看函数源代码
- Win10 配置Java JDK 16 环境变量
- 为什么onenote一直在加载_超好用的笔记软件,Onenote是我的最爱
- 为什么批量注册Facebook账号需要使用指纹浏览器防关联
- LeetCode 606. 根据二叉树创建字符串题解