常微分方程第三章

合集下载

第三章常微分方程的差分方法15

第三章常微分方程的差分方法15
第三章 常微分方程的差分方法
1.教学内容:
Euler方法:Euler公式,单步显式公式极其局部截断误 差;后退Euler公式,单步隐式公式极其局部截断误差;梯 形公式,预测校正公式与改进Euler公式。
2.重点难点:
Euler公式,预测校正公式与改进Euler公式
3.教学目标:
了解欧拉方法的几何意义、对给出的初值问题,能利 用Euler公式,改进Euler公式进行数值求解
科学技术当中常常需要求解常微分方程的定解问题。这类
问题的最简单的形式,是本章着重要考察的一阶方程的初值 问题:
y ' f x, y
y
x0
y0
(1) (2)
本章中我们假定右函数适当光滑以保证初值问题解的存
在唯一。虽然求解常微分方程有各种各样的解析方法,但求 解从实际问题中归结出来的微分方程要靠数值解法。
(其解析解为) y 2x 1
解:设步长 h=0.1,由改进的欧拉格式(10)有:
y
p
yn
h( yn
2xn ) yn

yc
yn
h( y p
2 xn1 ) yp
yn
1
1 2
(yp
yc )
n=0时
yp
y(xn ))
替代方程
y' (xn1) f (xn1, y(xn1))
中的导数项 y'xn1 再离散化,即可导出下列格式
yn1 yn hf xn1, yn1
(5)
该格式右端含有未知的 yn1 它实际上是个关于 yn1
的函数方程。故称该格式为隐式欧拉格式。
由于向前差商和向后差商具有同等精度,故隐式欧拉 格式也是一阶方法,精度与欧拉格式相当。但计算远 比显式格式困难得多。

常微分方程考研讲义第三章-一阶微分方程解的存在定理

常微分方程考研讲义第三章-一阶微分方程解的存在定理

第三章一阶微分方程解的存在定理[教学目标]1.理解解的存在唯一性定理的条件、结论及证明思路,掌握逐次逼近法,熟练近似解的误差估计式。

2.了解解的延拓定理及延拓条件。

3.理解解对初值的连续性、可微性定理的条件和结论。

[教学重难点] 解的存在唯一性定理的证明,解对初值的连续性、可微性定理的证明。

[教学方法] 讲授,实践。

[教学时间] 12学时[教学内容] 解的存在唯一性定理的条件、结论及证明思路,解的延拓概念及延拓条件,解对初值的连续性、可微性定理及其证明。

[考核目标]1.理解解的存在唯一性定理的条件、结论,能用逐次逼近法解简单的问题。

2.熟练近似解的误差估计式,解对初值的连续性及可微性公式。

3.利用解的存在唯一性定理、解的延拓定理及延拓条件能证明有关方程的某些性质。

§1 解的存在性唯一性定理和逐步逼近法微分方程来源于生产实践际,研究微分方程的目的就在于掌握它所反映的客观规律,能动解释所出现的各种现象并预测未来的可能情况。

在第二章介绍了一阶微分方程初等解法的几种类型,但是,大量的一阶方程一般是不能用初等解法求出其通解。

而实际问题中所需要的往往是要求满足某种初始条件的解。

因此初值问题的研究就显得十分重要,从前面我们也了解到初值问题的解不一定是唯一的。

他必须满足一定的条件才能保证初值问题解的存在性与唯一性,而讨论初值问题解的存在性与唯一性在常微分方程占有很重要的地位,是近代常微分方程定性理论,稳定性理论以及其他理论的基础。

例如方程dydx=过点(0,0)的解就是不唯一,易知0y =是方程过(0,0)的解,此外,容易验证,2y x =或更一般地,函数20 0() c<1x cy x c x ≤≤⎧=⎨-≤⎩ 都是方程过点(0,0)而且定义在区间01x ≤≤上的解,其中c 是满足01c <<的任一数。

解的存在唯一性定理能够很好地解释上述问题,它明确地肯定了方程的解在一定条件下的存在性和唯一性。

常微分方程第三章基本定理

常微分方程第三章基本定理

THANKS
感谢观看
线性化定理
总结词
线性化定理是将非线性常微分方程转化为线性常微分方程的方法,从而可以利用线性方程的解法来求解。
详细描述
线性化定理提供了一种将非线性常微分方程转化为线性常微分方程的方法。通过适当的变换,可以将非线性问题 转化为线性问题,从而可以利用线性方程的解法来求解。这个定理在解决复杂的非线性问题时非常有用,因为它 简化了问题的求解过程。
02
CATALOGUE
常微分方程的稳定性
稳定性定义
稳定性的定义
01
如果一个常微分方程的解在初始条件的小扰动下变化不大,那
么这个解就是稳定的。
稳定性的分类
02
根据稳定性的不同表现,可以分为渐近稳定、指数稳定、一致
稳定等。
稳定性判别方法
03
可以通过观察法、线性化法、比较法等方法来判断常微分方程
的解是否稳定。
龙格-库塔方法
总结词
龙格-库塔方法是常微分方程数值解法中一种更精确的 方法,它通过多步线性近似来逼近微分方程的解。
详细描述
龙格-库塔方法的基本思想是利用已知的初值和微分方 程,通过多步线性插值来逼近微分方程的解。具体来 说,龙格-库塔方法通过递推公式来计算微分方程的近 似解,公式如下:(y_{n+1} = y_n + h f(t_n, y_n) + frac{h^2}{2} f(t_{n-1}, y_{n-1}) - frac{h^2}{2} f(t_{n-2}, y_{n-2})) 其中 (h) 是步长,(t_n) 和 (y_n) 是已知的初值,(f) 是微分方程的右端函数。
存在唯一性定理表明,对于任意给定的初值问题,存在一个唯一的解,该解在某个区间内存在并连续 。这个定理是常微分方程理论的基础,为后续定理的证明提供了重要的依据。

第三章 常微分方程的差分方法

第三章 常微分方程的差分方法
P1 P1 P0 Pi+1 Pn y=y(x) Pi Pn Pi Pi+1
Euler法的求解过程是:从初始点 P0(即点(x0,y0))出发,作积分曲线 y=y(x)在P0点上切线 P0 P (其斜率 1 为 y( x0 ) f ( x0 , y0 ) ),与x=x1直线
x0
x1
xi
xi+1
自 动 化 工 程 学 院
School of Automation Engineering
第 三 章
P1 P1 P0
常微分方程的差分方法
Pi+1 Pn Pi Pi+1 Pi y=y(x) Pn
x0
x1
xi
xi+1
xn
由此获得了P2的坐标。重复以上过程,就可获得一系列的 点:P1,P1,…,Pn。对已求得点 Pn ( xn , y n ) 以 y ( xn ) = f ( xn , yn )为斜率作直线 当 x xn1 时,得 取 y( xn ) y n
第 三 章
常微分方程的差分方法
第三章 常微分方程的差分方法
引言
包含自变量、未知函数及未知函数的导数或微分的方
程称为微分方程。在微分方程中, 自变量的个数只有一个, 称为常微分方程。自变量的个数为两个或两个以上的微分 方程叫偏微分方程。微分方程中出现的未知函数最高阶导 数的阶数称为微分方程的阶数。如果未知函数y及其各阶导 数
对于初值问题
散化,建立求数值解的递推公式。递推公式通常有两类,一 类是计算yi+1时只用到xi+1, xi 和yi,即前一步的值,因此有了 初值以后就可以逐步往下计算,此类方法称为单步法;其代 表是龙格—库塔法。另一类是计算yi+1时,除用到xi+1,xi和yi以 外,还要用到 xi p , yi p ( p 1,2,, k ) ,即前面k步的值,此类 方法称为多步法;其代表是亚当斯法。

常微分方程(王高雄)第三版 3.3

常微分方程(王高雄)第三版 3.3
初值问题
dy f ( x, y ) , dx y ( x0 ) y0 (3.1) '
的解y ( x, x0 , y0 )都在区间 [a, b]上存在, 并且 ( x, x0 , y0 ) ( x, x0 , y0 ) , x [a, b] 则称初值问题(3.1) '的解y ( x, x0 , y0 )在点( x0 , y0 )
前提 解存在唯一
y0 ( x0 , x, y )
证明 在(3.1)满足y ( x0 ) y0的解存在区间内任取一值x1 ,
y1 ( x1 , x0 , y0 ), 则由解的唯一性知, (3.1)过点( x1 , y1 )与过点( x0 , y0 )的解是同一条积分曲线 , 即此解也可写成: y ( x, x1 , y1 ), 且显然有: y0 ( x0 , x1 , y1 ),
2 定理1 (解对初值的连续依赖性定理)
方程 条件: I. f ( x , y ) 在G内连续且关于 y满足局部Lips.条件;
dy f ( x, y) , dx ( x, y) G R2 (1)
II. y ( x , x0 , y0 ) 是(1)满足( x0 , y0 ) G 的解,定义
C 时,有 S G G 覆盖定理,存在N,当G i i 1 对 0 ,记 y , S ), min , / 2 d (G
N
Ci
G
L max L1,, LN 则以 为半径的圆,当其圆心从S的
G
左端点沿S 运动到右端点时,扫过 的区域即为符合条件的要找区域D
0
义, 其中 a x0 b, 则对 0, ( , a, b) 0, 使当

常微分方程

常微分方程

第三章 全纯函数的积分表示§3.1习题1.计算积分23z dz z γ-⎰,其中,γ为(iii)沿圆周{:||2}z z ∈= 的正向.23z dz z γ-⎰=23dz dz zγγ-⎰⎰=0-6i 6i ππ=-2.计算|z|=1dz ,2z +⎰并证明,012cos 54cos d πθθθ++⎰=0.奇点在区域外积分为0|z|=1dz 2z +⎰=20d 2i i e e πθθ+⎰=20i 2i i e d e πθθθ+⎰=20(2)i 0(2)(2)i i i i e e d e e πθθθθθ--+=++⎰ 20202022020120122412cos 2sin 052cos 2sin 2cos()2sin()12cos 2sin 054cos 12cos 2sin 054cos 54cos 12cos ()54cos i i i e d e e i d i i i d d d i d πθθθπππππππθθθθθθθθθθθθθθθθθθθθθ-+⇒=+++++⇒=+++-+-++⇒=++⇒+=++++=+⎰⎰⎰⎰⎰⎰⎰由对称性2012cos 12cos 54cos 54cos 12cos 054cos d d d ππππθθθθθθθθθ++=+++=+⎰⎰⎰所以3.计算积分|z|321dz (1z z z =--⎰).|z|3|z|3|z|=3||321111dz=4i (1(1)z 1z z z z dz dz dz z z z z z π===--+=+=---⎰⎰⎰⎰)4.如果多项式Q(z)比多项式P(z)高两次,试证: ||(z)lim0()R Z RP dz Q z →∞==⎰证明:2(z)lim ()R P z A Q z →∞= (A 为某个常数)从而存在M>0,R 0>0使得当R ≥R 0时,有2()()P z z M Q z ≤ 因此2||||||()()2|||||||0()(z)()|z |z R z R z RP z P z M Mdz z dz R Q Q z R π===≤≤=→→∞⎰⎰⎰ 所以()lim 0()R P z dz Q z →∞=⎰6.设f 1(),C D ∈γ是域D 中分别以a 和b 为起点和终点的可求长曲线.证明:()()()()z f z f z dz d z f b f a z γ∂∂+=-∂∂⎰ z)1f ),2f fi dz dx idy z x y ∂∂∂=-=+∂∂∂(( )1),2f z f fi d z dx idy z x y∂∂∂=+=-∂∂∂(( ()()11()()()()22f z f z f f f fdz d z i dx idy i dx dy z x y x y zf fdx dy df x y∂∂∂∂∂∂+=-+++-∂∂∂∂∂∂∂∂=+=∂∂故()(){}()()f z f z dz d z df f b f a z z γγ∂∂+==-∂∂⎰⎰8.设γ是域D 中以a 为起点,以b 为终点的可求长曲线,f,g 1()().H D C D ∈⋂证明分部积分公式:'f ()()()()|z).ba z g z dz f z g z f dz γγ=-⎰⎰’(()'()'()()['()()()'()][()()]'()()|b af zg z dz f z g z dz f z g z f z g z dz f z g z dz f z g z γγγγ+=+==⎰⎰⎰⎰故()'()()()|'()()b a f z g z dz f z g z f z g z dz γγ=-⎰⎰9.设γ是正向可求简单曲线,证明:γ内部的面积为12i rzdz ⎰. 证明:由公式得rzdz ⎰=z(dz Dd z z ∂-∧∂⎰)=-Ddz d z ∧⎰ =-()(i Ddx idy dx idy +∧-∧⎰⎰⎰DD)=2dx dy=2idA=2iA 所以A=12i rzdz ⎰11.设f 在z 0处连续,证明:(i )200r 001lim()z 2i f z re d f πθθπ→+=⎰();(ii) 0000||r1()lim()2r z z f z dz f z i z z π→-==-⎰. 证明:(i)220000011()(|()(|d 22i i f z re d f z f z re f z ππθθθθππ+-≤+-⎰⎰)) 00||sup |()()|0(0)z z rf z f z r +-=≤-→→所以200r 001lim()()2i f z re d f z πθθπ→+=⎰.(ii) 0200||r 01()1()22i z zf z dz f z re d i z z πθθππ-==+-⎰⎰故00r 00||1()lim()2z z r f z dz f z i z z π→-==-⎰12.设D={00z :arg()}(02),z a a θθαπ∈<-<+<≤ f 在\{}D a 上连续,证明: (i )如果a z D\{a}lim ()(),z z a f z A →∈-=那么r 0|z a|lim();r z Df z dz i A α→-=∈=⎰(ii)如果z D\{a}lim ()(),z z a f z B →∞∈-=那么z ||limf ().z a R z Dz dz i B α→∞-=∈=⎰证明:(i )||||||1()i ()|()()|||z a r z a r z a rz z z A f z dz A f z dz z a f z A dz z a r α-=-=-=∈-=-≤---⎰⎰⎰ ||sup |(z )()|0r 0)z a rz a f z A +-=∈≤--→→((ii)略§3.2习题1.计算积分: (i )2|z|||,||||rdz a r z a =≠-⎰ 222|z|0||||||()()()()i i r z r dz rd r dzz a i re a re a z a r az πθθθ-====-----⎰⎰⎰ 22222||112||||||/z r ir r dz r a z a r a z r a π=-⎛⎫=-= ⎪----⎝⎭⎰ (iii )4|z|5;1zdzz =-⎰ 2242222|z|5|z|5|z|51111012(1)(1)411zdz dz dz z z z z z ===⎛⎫==-= ⎪--+-+⎝⎭⎰⎰⎰(iv )22|z|2,0.zae dz a z a =>+⎰ 22|z|2||2||2112ai 2z z z a z a z ae e e dz dz dz z a z ai ai z ai ====-+-+⎰⎰⎰ =112i 22sin 2ai 2ai ai i e i e a ai aπππ-⋅⋅-⋅=2.设f 在{:||}z r z ∈<<∞ 中全纯,且z lim (z).zf A →∞=证明:||()2,z Rf z dz iA π==⎰这里r R >证明:设',r R R <<<∞利用Cauchy 积分公式得'||||'||dz 2()z Rz R z R A f z iA f z dz dz z π===-=-⎰⎰⎰()|z|'1|()|||R zf z A dz R =≤-'⎰ ''||2sup |()|0()z R zf z A R π=≤⋅-→→∞4.设0,r R f <<在(0,)B R 中全纯.证明:201i (0)();2i f f re d πθθπ=⎰()||ii ().z rf f z dxdy π<⎰21()(0)=r 证明:(i )0,R δ∀>> 由Cauchy 积分公式得||||1()1()2i 2z r z f z f z dz dz z i z δππ===⎰⎰ 即220011()()22i i f re d f e d ππθθθδθππ=⎰⎰,令0δ→,则有201()(0)2i f re d f πθθπ=⎰(ii)2|z|1()r r f z dxdy π<⎰=r2201()r i d f e d πθρρρθπ⎰⎰ (1) 利用(i )得20()2(0)i f e d f πθρθπ=⎰(2)由(1)和(2)式得2||1(0)()z rf f z dxdy r π<=⎰5.设u 是B(0,R)中的调和函数,0.r R <<证明:201u(0)=().2i u red πθθπ⎰证明:单连通域上的调和函数是某个全纯函数f 的实部:Re u f =22011(0)()0Re (0)()22i i f f re d u f u re d ππθθθθππ=⇒==⎰⎰()6.设0 1.r <<证明:20log(12rcos r )0d πθθ-+=⎰证明:令U(z)=2ln |1|z -,则U 在B(0,1)上调和,由题5,0=U(0)=2222111()log(12cos )log(12cos )22i u re d r rd r rd πππθθθθθθπππ=-+=-+⎰⎰⎰.1.设D 是域,,()f g H D ∈.如果'fg 在D 上有原函数ϕ.证明:'fg 在上D 有原函数fg ϕ-证明:,()f g H D ∈.由'fg =ϕ,得''()fg fg ϕ-=.3.设()()n f C H ∈ ,并且()()0n f z ≡.证明:f 必为次数不大于n -1的多项式. 证明:归纳法 k=1时显然,设k=n-1时成立,取定C ξ∈,由()(1)()0()n n f z f z C -≡⇒=故(1)1()0(1)!n n Cz f z n --⎡⎤-≡⎢⎥-⎣⎦,从而1()(1)!n Cz f z n ---是次数不超过n-2的多项式.5.设f 是凸域D 上的全纯函数,如果对每点z D ∈,有{}'Re ()0f z >,那么f 是D 上的单叶函数.证明:12,z z D ∀∈,12z z ≠ 则 211''21121210()()()(())()z z f z f z f d f z t z z z z dt ξξ-==+--⎰⎰则1'21121021()()(())f z f z f z t z z dt z z -=+--⎰ 从而{}1'21121021()()Re (())0f z f z f z t z z dt z z -≥+->-⎰ 故12()()f z f z ≠,这表明f 是D 上的单叶函数.1. 计算下列积分:(1)2||1sin 1z zdz z =-⎰ 解:2||1|1|11sin 1sin sin 2sin11111z z z z z z dz dz i i z z z z ππ=-==⎛⎫=⋅== ⎪--++⎝⎭⎰⎰ (2)2||2;1z dz z =+⎰ 解:2||2||2111012z z dz dz z i z i z i ==⎛⎫=-= ⎪+-+⎝⎭⎰⎰ (4)223||2;(1)(4)z dzz z =++⎰解:与第二题类似,答案为0.(6)||()()nz Rdzz a z b =--⎰,n 为正整数,a ≠b 不在圆周|z|=R 上. 解:原式0,2()2()0,.nn a b i a b a i a b a a b ππ⎧⎪⎪⎪-=⎨⎪-⎪-⎪⎩均在圆外.在圆外,b 在圆内.在圆内,b 在圆外.均在圆内3.设D 是由有限条可求长简单闭曲线围成的域,1,n z z …,是D 中n 个彼此不同的点.如果()()f H D C D ∈ ,证明:()()1()()2()n n n D z f P z d i zωζωζζπωζζ∂-=-⎰是次数不超过n-1的多项式,并且1()(),1,2,.()()(-)k k n n P z f z k n z z z z z ω==⋯=-⋯,其中,. (提示:证明()()n n z zωζωζ--是z 的次数不超过n-1的多项式.)证明:由于1()()n z z z ω=-n …(z-z )从而()()n n z z ωζωζ--是z 的n-1次多项式,记()()h(,)()nn z z f zωζωξξζ-=-取0ε>充分小,由Cauchy 积分公式11k 11||11(,)()(,)()2()k nn n k k j n j k z j kjj h z P z d h z z z z i z ξεξξπξ-===-=≠===∏-∏-∑∑⎰ 因为()h(,)()n z z f zωξξζ-=-是z 的次数不超过n-1的多项式,故P(z)是关于z 的次数不超过n-1的多项式. 又1()()(),nn k j j j kz z z z z ω=≠-=-∏-故()()k k P z f z =是关于z 的次数不超过n-1的多项式.5.设((0,1))((0,1))f H B C B ∈ .证明:2202202(1)()cos 2(0)(0);22(2)()sin 2(0)(0)2i i f e d f f f e d f f πθπθθθπθθπ⎛⎫'=+ ⎪⎝⎭⎛⎫'=- ⎪⎝⎭⎰⎰(提示:分别计算积分1111112()2()22d d f f i i ζζζζζζζζπζζπζζ==⎛⎫⎛⎫++-- ⎪ ⎪⎝⎭⎝⎭⎰⎰和 即可.)证明:由Cauchy 公式,得220011()1()1(0)()222i i i i f f e f d ie d f e d i i e i θππθθθζζζθθπζππ====⎰⎰⎰,22011()1(0)()22i i f f d f e e d i πθθζζζθπζπ-='==⎰⎰②又由Cauchy 定理,1()0f d ζζζ==⎰即201()02i i f e e d πθθθπ=⎰③②+③得22211(0)()cos ()(2cos 1)2i i f f e d f e d ππθθθθθθππ'==-⎰⎰即22221()cos(0)()(0)2(0)2i i f e d f f e d f f ππθθθθθππ''=+=+⎰⎰6.利用上题结果证明:设((0,1))((0,1))f H B C B ∈ ,且(0)1,Re ()0f f z =≥,那么 2Re (0)2f '-≤≤. 证明:2202()cos (0)2(0)2i f e d f f πθθθπ'=+⎰两边取实部,即2202Re()cos 2Re (0)02i f e d f πθθθπ'=+≥⎰同理2202()sin 2(0)(0)2i f e d f f πθθθπ'=-⎰ 2202Re (0)2()sin 2022i f f e d πθθθπ'=-≤-=⎰所以2Re (0)2f '-≤≤.§3.5习题1.设f 是有界整函数,12,z z 是(0,)B r 中任意两点.证明:12()0()()z rf z dz z z z z ==--⎰并由 得出Liouyille 定理.证明:利用Cauchy 积分公式得 1212121212()()()1()()2()()z r z r f z f z f z f z f z dz dz i z z z z z z z z z z z z π==⎛⎫-=-= ⎪------⎝⎭⎰⎰ 另一方面, 由于f 有界,0,..(),M s t f z M z ∃>≤∀∈C由Cauchy 积分定理121212()()20()()()()()()()z rz R rf z f z M dz dz R R z z z z z z z z R z R z π==>=≤⋅→→∞------⎰⎰从而 121212()()20()()()f z f z i f z f z f z C z z π-=⇒=⇒=-2.设f 是整函数, 如果当z →∞时,()(||),0,f z O z αα=≥证明f 是次数不超过[]α的多项式.解:令[]1n α=+ ()()()111!()!()!||()2()2()2!0n n n n z R z Rz R nn f n f n M f z d d d i z z R R z n M R R αζζζαζζζζζζπζπζπ+++-=-=-==≤≤--+≤→→∞⎰⎰⎰ 故()()0n fz ≡4.设f 是整函数,如果{}();Im 0f z z ⊂∈> ,证明f 是一个常值函数. 证明:令()()()f z i g z f z i-=+,则|g(z)|<1,g 是整函数. 从而()()()()f z i g z c f z c f z i-==⇒=+5.设f 是整函数,如果()\[0,1]f ⊂ 证明f 是一个常值函数.证明:令()()()1f z g z f z =-,则 ()[)\0,g ⊂∞ .再令()h z =则(){}:Im 0h z z ⊂∈>由上题知h 常值,故f 常值.6.设f 在域D 上全纯,0z D ∈定义00000()(),\{};()(),.f z f z z D z z z F z f z z z -⎧∈⎪-=⎨⎪'=⎩证明:()F H D ∈证明:0000000()()()()lim lim ()z z z z f z f z F z f z F z z z →→-'===- 故F 在0z 点也连续.将F 限制在0(,)B z ε上,则F M ≤,对D 内任一简单闭曲线γ,可取一含于0(,)B z ε的简单闭曲线使得()()f z dz f z dz εγγ=⎰⎰,对0,ε∀>由此易得()0f z dz γ=⎰,从而F 在D 上全纯.7.设γ是可求长曲线,f 在域D 上连续,在\D γ上全纯.证明:f 在D 上全纯. 证明:任取D 中简单闭曲线0γ(1) 当γ含于0γ内部时,延长γ,交0γ于A,B 两点. 012()()()0AB BA f z dz f z dz f z dz γγγ++=+=⎰⎰⎰ (2) 同理,当γ,0γ相交时,0()0f z dz γ=⎰ 故由Morera 定理知f 在D 上全纯.1、设D 是由有限条可求长简单闭曲线围成的域,1()f C D ∈.证明: (i )()()D Df d d f d ξξξξξξ∂∂∧=∂⎰⎰⎰ (ii )()()D Df d d f d ξξξξξξ∂∂∧=∂⎰⎰⎰ 证明:()f ξ是域上的一个0次微分形式1()f C D ∈,根据定理3.6.1D Df df ∂=⎰⎰ 可得()()()()()()D D D D f f f f d df d d d d d d ξξξξξξξξξξξξξξξ∂∂∂∂=∧=+∧=∧∂∂∂⎰⎰⎰⎰ 同理()()()()()D D Df f f f d d d d d d ξξξξξξξξξξξξξ∂∂∂∂=+∧=∧∂∂∂⎰⎰⎰ 4、设D 是由有限条可求长简单闭曲线围成的域,f 在D 上全纯,z D ∈.证明:()()()()D Df f f d d d z z z ξξξξξξξξξ∂'∧=+---⎰⎰⎰ 证:因()f H D ∈,从而1()f C D ∈,则对f 用Pompeiu 公式1()1()1()22D D f f f z d d d i z i zξξξξξπξπξξ∂∂=+∧-∂-⎰⎰⎰ 从而两边取共轭得 1()1()()22D D f f f z d d d i i z zξξξξξππξξ∂'=--∧--⎰⎰⎰ 1()1()()22D D f f f z d d d i z i z ξξξπξπξ∂'=-+∧--⎰⎰⎰ (1) 1()()2D f f z d i zξξπξ∂=-⎰ (2) 由(1)(2)可得()()()()D D f f f d d d z z z ξξξξξξξξξ∂'∧=+---⎰⎰⎰.2.设D 是域,()f C D ∞∈, 证明:若0()u C D ∞∈是非齐次∂方程的解,即0u f z ∂=∂, 则该方程的解的全体为0()u H D +.提示:显然()0u H D +中的元都是解.另一方面,设u 是方程的一个解,则()00u u f f z∂-=-=∂,即0u u -满足C-R 方程, 又()u C D ∞∈,故()0u u H D -∈.。

3. 一阶常微分方程解的存在唯一性

3. 一阶常微分方程解的存在唯一性

由于ϕn(x) = ϕ0(x) + ϕ1(x) − ϕ0(x) + ϕ2(x) − ϕ1(x) + · · · + ϕn(x) − ϕn−1(x) ,

故只需证明无穷级数ϕ0(x) + [ϕn+1(x) − ϕn(x)]在I上一致收敛即可。采用数学 n=0
归纳法来证明:
特别,取ϕ0(x) = y0,则 |ϕ1(x) − ϕ0(x)| =
设ϕ(x)和ψ(x)都 是 微 分 方 程(3.1)在I上 的 解。 记M = max |ϕ(x) − ψ(x)|, 根 x∈I
据Lipschitz条件,当x ∈ I时,有
|ϕ(x) − ψ(x)| ≤
x
|f (t, ϕ(t)) − f (t, ψ(t))|dt
x0 x
≤ L |ϕ(t) − ψ(t)|dt
第三章 一阶常微分方程解的存在唯一性
本章主要介绍和证明一阶微分方程解的Picard存在和唯一性定理,解的延拓,解对 初值的连续性和可微性等概念。
3.1 Picard存在唯一性定理
3.1.1 一阶显式微分方程
考虑一阶显式常微分方程的初值问题

dy dx
=
f (x, y)
y|x=x0 = y0
(3.1)

LnM n!
x
|t − x0|ndt
x0
=
LnM (n + 1)!
|x

x0|n+1
特别,当|x − x0| ≤ h时,
|ϕn+1(x)

ϕn(x)|

LnM (n + 1)!
hn+1

由于正项级数

常微分第三章第4节(奇解)

常微分第三章第4节(奇解)

y cx f (c)
' x f ( p) 0 ' 如果 x f ( p) 0,则 y xp f ( p) 消去 p 得到方程的另一个解。
这里
c 是任意常数。
25
注意,求得此解的过程正好与从通解
y cx f (c)
可以验证,此解的确是通解 中求包络的过程一样。 的包络。
注1:包络一定包含在 c-判别曲线中。 注2:c-判别曲线不一定为包络。
2 充分但不必要。 2 x y 0
15
C 判别曲线法求方程奇解的一般步骤: (1)求出方程的通解(积分曲线族); (2)求积分曲线族的 c 判别曲线;
(3)检验 c 判别曲线是否为包络,若是,则 为方程的奇解。
16
y 1
其中
容易求得原方程的通解为 y sin( x c)
c 为任意常数。而 y 1 是通解的包络。
所以此两直线都是方程的奇解。
22
例4
求方程
dy dy 2 y 2x ( ) 的奇解。 dx dx
y 2 xp p 2 解 从 2 x 2 p 0
消去
这是克莱罗方程,因而它的通解是 1 y cx c 1
27
y2 4x
O
图(3.5)
28
例6
求一曲线,使在其上每一点的切线截割坐
标轴而成的直角三角形的面积都等于2 。 y A
O 图(3.6)
B
x
29
依题意有ab 4,而
dy 2 dy 得 ( y x ) 4 dx dx
现在求曲线族的包络,亦即微分方程的奇解。
30
y 2c c 2 x 从 中消去 c 得微分方程的奇解 1 cx 0

第三章 常微分方程的边值和本征值问题

第三章 常微分方程的边值和本征值问题

因此比 较明智的做法是,在每一个试验本征值上,由 xmax
出发向后直接积分产生另一个数值解 Ѱ>。 为了判断 这个试验本征值是不是一个能量本征值,可以在一
个接合点 xm上比较 Ѱ<和 Ѱ>,其中接合点 xm要这样选择, 使得两个积分都是准确的。这里接合点 xm 的一个方便的选 择是左转折点或右转折点。
问题转化为求下面方程的根
Φk (1)= 0
3.3 一维薛定谔方程的定态解
一维位势 V(x) 中一个质量为 m 的粒子的 量子力学定态
在 x = xmin 和 x = xmax 处两点位势变为无穷大,也就是说在这 两点上有刚壁,在 这两点之间则是一个势阱。
定解问题
其中
求使这个问题有非零解的能量本征值 E 及其相应的波函数
Ѱ<和 Ѱ>的归一化总是可以这样选择,使得两个函数值在
xm 上相等。这时如果 它们的微商在 xm上也相等,那么就可 以断言这个试验本征值就是能量本征值.
数学表达式为
这里的
提供了一个方便的标尺
打靶法的基本思想是将边值问题当作一个含可调参数 δ 的
初始问们就可以通过积分这个初始问
题得到 yδ (b) .
一般来说,由于可调参数 δ 的随意选择, yδ(b) 和 yb 很难相等。
打靶法就是通过使用一个搜索算法去调整参数 δ ,使得 yδ (b) 和 yb 在误差容忍范围内相等,从而达到数值求解边 值问题的目的. 问题转化为求下面方程的根
3.2 打靶法求解本征值问题
考虑一根密度均匀的绷紧的弦的振动,分离变量后,空间
部分满足的方程和边界条件可以写成
φ 是弦的横向位移, k 是波数 解析解为
相比边值问题,本征值问题多了一个待定参数 策略:我们先猜测一个试验本征值 k,同时任取一个非零数 δ , 把微分方程变化为一个初始值问题

第三章 常微分方程

第三章 常微分方程

(3) 日治愈率为 , 病人治愈后仍可被感染
推知:
di N Ni (1 i ) Ni dt
di i (1 i ) i dt s (t ) i (t ) 1
目录 上页 下页 返回 结束
SIS模型的解
当 时,有
1 i (t ) [ ( )e ( )t ]1 i0
目录
上页
下页
返回
结束
人口分类(在此用绝对人数)
1、易感类S,该类成员没有感染SARS,也没有免疫
能力,处于自由环境,可以被传染上SARS
2、潜伏期类E,该类成员已经感染了SARS病毒,但 处于潜伏期,还不是SARS患者,不能把病毒传 染给S 类成员
3、患病未被发现类Iu,该类成员已经成员真正的 SARS患者,可以把病毒传染给S 类成员
目录 上页 下页 返回 结束
二、模型的建立—机理分析法
1. SI模型
假设条件
() 1 人群分为易感者s(t )和已感者i(t )两类
(2) 每个病人每天有效接触的平均人数为
推知:
每个病人每天可使 s(t )个健康者变为病人,病人数
为Ni(t ),所以每天有 Ns(t )i(t )个健康者变为病人.
目录 上页 下页 返回 结束
关于SARS的假设
1、所研究地区的人口总量一定,不考虑人口迁入 迁出、自然出生和自然死亡 2、结合SARS的病理特性,假设其传播方式为接 触性传播,不与患病者接触就不会被感染 3、假设病人被感染后需先进入潜伏期,在潜伏期 内不具备传染性
目录 上页 下页 返回 结束
4、假设易感者与患病者在人群中混合充分均匀, 易感者感染的机会与接触患病者的机会成正比, 并且传染率为常数 5、假设SARS患者被发现后就立即被隔离,被隔离 者不具备传染性,SARS患者只在被发现前可以 传染他人 6、假设SARS康复者不会被再次感染,且无传染性

常微分方程--第三章 一阶微分方程的解的存在定理(3.1-3.2)_OK

常微分方程--第三章 一阶微分方程的解的存在定理(3.1-3.2)_OK

x
L x0 1( ) 0 ( )d
L
x x0
M (
x0 )d
ML 2
(x
x0 )2
其中第二个不等式是由Lipschitz条件得到的,
由Lipschitz条件
17
设对于正整数n, 有不等式
n (x) n1(x)
MLn1 n!
(x
x0
)n
,
则当x0 x x0 h时,由Lipschitz条件有
dy dx
f
(x, y), (3.1)
y(x0 ) y0
证明: 若y (x)为(3.1)的连续解,则
d ( x)
dx
f
( x, ( x)),
(x0 ) y0
对第一式从x0到x取定积分得
x

x (x) (x0 ) x0 f (x,(x))dx (x) y0 x0 f (x,(x))dx
x
f ( , ( )) f ( ,( )) d x0
x
x
L ( ) ( ) d L g( )d
x0
x0
令u(x) L
x
g( )d ,
x0
则u(x)是定义于[x0, x0 h]上连续可微函数,
且u(x0 ) 0,0 g(x) u(x), u'(x) Lg(x),于是
u(x) Lu(x), (u(x) Lu(x))eLx 0,
(4) (x)是积分方程(3.5)定义于[x0 h, x0 h]上连续解
且唯一.
9
下面分五个命题来证明定理,为此先给出
积分方程
如果一个数学关系式中含有定积分符号且在定积分符 号下含有未知函数, 则称这样的关系式为积分方程.

数学物理方法_第3章 二阶线性常微分方程的幂级数解法本征值问题

数学物理方法_第3章 二阶线性常微分方程的幂级数解法本征值问题

y ( x) 2 1a2 3 2a3 x (k 2)( k 1)ak 2 x k
把以上结果代入方程,比较系 数得 2 2
2 1a2 a0 0, 3 2a3 a1 0, 4 3a4 2 a2 0, 5 4a5 2 a3 0,
(2k 1)!
a1.
于是方程的级数解为
1 1 1 y( x) a0 1 ( x)2 ( x) 4 (1) k ( x) 2 k 4! (2k )! 2! 2 k 1 a1 1 1 3 5 k ( x) x ( x) ( x) (1) 3! 5! (2k 1)! a a0 cos x 1 sin x.
n 1


n 1
cn1 (n 1)( x x0 )n ,
n 0

可将式(3.1.4)写成
c
n 0 n n
n ( n 2)( n 1)( x x ) [ ( k 1) a c ]( x x ) n2 0 n k k 1 0 n n 0 k 0
y( x) an x n
n 0
(3.3.2)
于是
y( x) nan x
n 1 n 1
(k 1)ak 1 x k ,
k 0

y( x) n(n 1)an x
n2
2 k 0

n2
(k 2)(k 1)ak 2 x k ,
(1 l )(l 2) 3 (3 l )(1 l )(l 2)(l 4) 5 y1 ( x) x x x 3! 5! (2k 1 l )(2k 3 l ) (1 l )(l 2)(l 4) (l 2k ) 2k 1 x (2k 1)!

常微分方程的差分方法

常微分方程的差分方法

x, y x

h 2

f
xn, y xn f
xn1, y xn1
再离散化,即可得如下计算公式
yn1

yn

h 2

f
xn ,
yn

f
xn1,
yn1
与梯形求积公式相呼应的这一差分格式称为梯形格式。
第三章 常微分方程的差分方法
§ 1 欧拉方法 § 2 改进的欧拉方法 § 3 龙格-库塔方法 § 4 亚当姆斯方法 § 5 收敛性与稳定性 § 6 方程组与高阶方程的情形 § 7 边值问题
数值分析简明教程
4.1
安工大2004
引言
科学技术当中常常需要求解常微分方程的定解问题。这类问题的
最简单的形式,是本章着重要考察的一阶方程的初值问题:
数值分析简明教程
4.9
安工大2004
七桥问题
• 七桥问题Seven Bridges Problem
18世纪著名古典数学问题之一。在哥 尼斯堡的一个公园里,有七座桥将普雷 格尔河中两个岛及岛与河岸连接起来( 如图)。问是否可能从这四块陆地中任 一块出发,恰好通过每座桥一次,再回 到起点?欧拉于1736年研究并解决了 此问题,他把问题归结为如下右图的“ 一笔画”问题,证明上述走法是不可能 的。
,虽然他被别人从火海中救了出来,但他的书房和大量研究成果全部化为 灰烬了.

沉重的打击,仍然没有使欧拉倒下,他发誓要把损失夺回来.欧拉完
全失明以后,虽然生活在黑暗中,但仍然以惊人的毅力与黑暗搏斗,凭着
记忆和心算进行研究,直到逝世,竟达17年之久
数值分析简明教程
4.11
安工大2004
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
一个连续解 , 则 ( x ) ( x ), x [ x 0 , x 0 h ].
如果f ( x , y )在R上关于y的偏导数f y ( x , y )连续, 则f ( x , y )在R上关于y满足Lipschitz条件.
考虑一阶隐方程
F ( x, y , y ) 0,
x
x0 y0
f ( x0 , y0 ) exp( exp(
x
f ( x, ) y
dx)
x0
f ( x, ) y
dx)
x0

已知方程
dy dx sin xy
y0
exp(
x
f ( x, ) y
x0 0
0
dx)
x0
试求
[
'
(3 .5 )
如果在点
0 '
( x 0 , y 0 , y 0 )的某邻域中满足
( x , y , y ) 连续且存在连续偏导数
'
'
1 F ( x , y , y ) 对所有变元
,
2 F ( x0 , y0 , y0 ) 0,
0 '
3
0
F ( x0 , y0 , y0 )
'
y
'
0,
则方程(3.5)存在唯一解
y y ( x ), x x 0 h ( h 为足够小的正数 )
满足初始条件
y ( x0 ) y0 ; y0 ( x0 ) y0 ,
' '
(3 .8 )
近似计算和误差估计
求方程近似解的方法---Picard逐步逼近法,这里
0 0 n ( x) y0
(x) y
(即对( x, y, ) G , 以( x, y, )为中心球C G , 使 f ( x, y, )在C内对y满足Lipschitz 条件, L与无关)
解对初值的可微性
若函数 f ( x , y )以及 f y 都在区域 G 内连续 , 则方程
( 3 . 1)的解 y ( x , x 0 , y 0 ) 作为 x , x 0 , y 0的函数在它们存 在范围内是连续可微的 .
解对初值和参数的连续性定理
设 f ( x , y , ) 在区域 G 连续 , 且在 G 内一致地关于 局部 Lipschitz 作为 x , x 0 , y 0 , 的函数在它们存在范围 y 满足
条件 , 则方程 ( 3 . 1 ) 的解 y ( x , x 0 , y 0 , ) 内是连续的 .
方程
F ( x, y, dy dx ) 0, ( 3 . 34 )
的奇解包含在由方程组
F ( x, y, p ) 0 (3 .3 4 ) ' Fp ( x, y, p ) 0
消去参数 p 而得到的曲线 ( x , y ) 0 之中 ,
此曲线称为
( 3 . 34 )的
p 判别曲线.
推论1 对定义在平面区域
dy f ( x, y) , dx y ( x0 ) y0
G 上的初值问题
其中 ( x 0 , y 0 ) G .
若 f ( x , y ) 在 G 内连续且关于
y 满足局部
Lipschitz
条件 ,
则它的任一非饱和解均可延拓为饱和解. 推论2 设 y ( x )为初值问题
向 x 增大 ( 减少 ) 一方的延拓来说
(1 ) 解 y ( x ) 可以延拓到区间
[ x 0 , ) (( , x 0 ],
(2)
解 y ( x ) 可以延拓到区间
[ x 0 , m ) (( m , x 0 ],
其中 m 为有限数 , 当 x m 时 , 或者 y ( m ) 无界 , 或者 ( x , ( x )) G
( x , y ) R
, 计算
b 1 b
2
f ( x , y ) 1 b ; h min{ a ,
}
故可取 a 1, b 1
此时由定理得到初值问 1 2 x 1 2 . 题的解的存在唯一区间 是
局部李普希茨(Lipschitz)条件
若函数 f ( x , y ) 在区域 G 内连续 , 且对 G 内的每 一点 P , 有以 P 为中心完全含于 G 内的闭矩形 R P 存在 , ,
(x) y
(3 .7 )

x x0
f ( , n 1 ( ))d
x0 x x0 h
( n 1, 2 , )
命题2 对于所有 n 和 x [ x 0 , x 0 h ], n ( x ) 连续且满足
n ( x) y0 b,
(3 .8 )
b M ), M max
( x , y ) R
,
这里 h min( a ,
f (x, y)
命题1 初值问题(3.1)等价于积分方程
y y0

x
f ( t , y )dt
x0
(3 .5 )
构造Picard逐步逼近函数列{ n ( x )}
0 0 n ( x) y0
.
这里 F ( x , y , p ) 是 x , y , p 的连续可微函数
注: p 判别曲线是否为方程的奇解,尚需进一步讨论.
克莱罗(Clairaut)方程
形如
dy y x f dx dx dy
的方程,称为克莱罗(Clairaut)方程.
这里 f ( p ) 是 p 的连续可微函数 .
对于一阶微分方程 F ( x , y , y ) 0,设 ( x , y , c ) 0 是 线族 ( x , y , c ) 0 的包络存在, 条奇积分曲线,即奇解 存在的话)也是通解的 。 包络。
它的通解。如果积分曲 则其包络就是方程的一 反之,方程的奇解(若
奇解的求法
, 其中 c 是参数 , ( x , y , c ) 是 x , y , c 的连续可微函数
曲线族(3.23)的包络是指这样的曲线, 它本身不包含在 曲线(3.23)中,但过这曲线的每一点有(3.23)中的一条 曲线和它在这点相切.
包络的求法
( x, y, c) 0,
(3.23)
曲线族(3.23)的包络包含在下列两方程
y ( x, x0 , y0 ) x0
]y
x0 0
0
0
,
[
y ( x, x0 , y0 ) y0
]y
0
.
.
解 f x ( x , y ) y cos xy , f y ( x , y ) x cos xy 在 xy 平面上连续
方程 dy dx sin xy 的解 y ( x , x 0 , y 0 ) 作为 x , x 0 , y 0 .
解对初值的连续性定理
方程 条件: 结论:
dy dx
f ( x , y ) 在G内连续且关于 y 满足局部Lips.条件;
y ( x , x0 , y0 ), ( x0 , y0 ) G ,作为x , x0 , y0的函数
f ( x, y) ,
( x, y) G R
2
(1)
在它的存在范围内是连续的.
1 x
2
0
(易见 y 0 是原方程的解
exp(

x
e2 xdx )
x0 0
0
0

[
y ( x, x0 , y0 ) x0
]y
0
x0ຫໍສະໝຸດ f ( x0 , y0 ) exp(
x x0
x
f ( x, ) y
x 0
0
dx)
x0
[ f ( x 0 , y 0 ) exp(
h
n 1
,
( 3 . 19 )
例 求初值问题
dy dx 1 y ,
2
y (0) 0
解的存在唯一区间.
解 对任意 a , b , 函数 f ( x , y ) 均在矩形区域
R {( x , y ) | x a , y b }
内连续 , 且对 y 有连续的偏导数
M max
2

f ( x, ) y
dx )] y 0 0
f ( 0 , 0 ) exp( f ( 0 , 0 ) exp(

x
x cos( x ( x , 0 , 0 )) dx ) xdx ) 0
0 x

0
包络和奇解
定义1:对于给定的一个单参数曲线族:
( x, y , c ) 0, ( 3 . 23 )

x x0
f ( , n 1 ( )) d
x0 x x0 h
( n 1, 2 , )
对方程的第
n 次近似解
n ( x ) 和真正解
n
y ( x )在
[ x 0 h , x 0 h ]内误差估计为
n (x) (x)
ML
( n 1 )!
R : x x0 a, y y0 b,
( 3.2)
并且对y满足Lipschitz条件 :
即存在L 0, 使对所有 x, y1 ), ( x, y2 ) R常成立 (
f ( x , y1 ) f ( x , y2 ) L y1 y2
则初值问题
( 3 . 1 ) 在区间 x x 0 h 上的解存在且唯一
相关文档
最新文档