常微分方程的边值问题
- 格式:docx
- 大小:134.51 KB
- 文档页数:9
常微分方程基本概念常微分方程(Ordinary Differential Equations,简称ODE)是数学分析中的一个重要分支,研究的是一元函数的导数与自变量之间的关系。
它在物理学、工程学、生物学等领域具有广泛的应用。
本文将介绍常微分方程的基本概念和相关知识。
一、常微分方程的定义常微分方程是描述未知函数的导数与自变量之间关系的方程。
一般形式可以表示为:dy/dx = f(x, y)其中,y是未知函数,x是自变量,f(x, y)是已知函数。
二、常微分方程的阶数常微分方程根据未知函数的最高阶导数的阶数不同,可以分为一阶、二阶、高阶等不同阶数的微分方程。
1. 一阶微分方程一阶微分方程是指含有一阶导数的方程。
一般形式可以表示为:dy/dx = f(x, y)例如,y' = 2x + 1就是一个一阶微分方程,其中y'表示y对x的一阶导数。
2. 二阶微分方程二阶微分方程是指含有二阶导数的方程。
一般形式可以表示为:d²y/dx² = f(x, y, dy/dx)例如,y'' + y = 0就是一个二阶微分方程,其中y''表示y对x的二阶导数。
三、常微分方程的初值问题和边值问题常微分方程除了描述函数的导数与自变量之间的关系外,还可以给出一些初始条件或边界条件,从而确定唯一的解。
1. 初值问题初值问题是指在微分方程中给出了函数在某一点的初值条件,要求求解出满足该条件的解。
一般形式可以表示为:dy/dx = f(x, y),y(x₀) = y₀其中,y(x₀) = y₀表示在点(x₀, y₀)处给定了函数的初始值条件。
2. 边值问题边值问题是指在微分方程中给出了函数在多个点的边界条件,要求求解出满足这些条件的解。
一般形式可以表示为:dy/dx = f(x, y),y(a) = y_a,y(b) = y_b其中,y(a) = y_a和y(b) = y_b表示在点(a, y_a)和(b, y_b)处给定了函数的边界条件。
18—1 常微分方程数值解法2§1 引言§2 Euler 方法§3 Runge -Kutta 方法§4 单步法的收敛性与稳定性§5 线性多步法§6 方程组与高阶方程的情况§7 边值问题的数值解法3§1 引言微分方程:关于一个未知函数的方程,方程中含有未知函数的(偏)导数,以及自变量等,其中关于未知函数导数的最高次数称为微分方程的阶数.例如:0)()(')()(''=++−x c y x b y x a x y4实际中,很多问题的数学模型都是微分方程. 常微分方程作为微分方程的基本类型之一,在理论研究与工程实际上应用很广泛. 很多问题的数学模型都可以归结为常微分方程. 很多偏微分方程问题,也可以化为常微分方程问题来近似求解.微分方程的应用情况5对于一个常微分方程:'(,) ,[,]dy y f x y x a b dx==∈为了使解存在,一般要对函数f 施加限制条件,例如要求f 对y 满足Lipschitz 条件:1212(,)(,)f x y f x y L y y −≤−6同时,一个有解的微分方程通常会有无穷多个解例如cos() sin(),dyx y x a a R dx=⇒=+∀∈为了使解唯一,需要加入一个限定条件. 通常会在端点出给出,如下面的初值问题:(,),[,]()dyf x y x a b dx y a y ⎧=∈⎪⎨⎪=⎩7常微分方程的解是一个函数,但是,只有极少数特殊的方程才能求解出来,绝大多数是不可解的.并且计算机没有办法对函数进行运算. 一般考虑其近似解法,一种是近似解析法,如逼近法、级数解法等,另一种是本章介绍的数值解法.8§2 Euler 方法92-1 Euler 公式对常微分方程初值问题:⎩⎨⎧==00')(),(y x y y x f y 数值求解的关键在于消除其中的导数项——称为离散化. 利用差商近似逼近微分是离散化的一个基本途径.10现在假设求解节点为),,1,0(m i ih a x i "=+=,其中ma b h −=为步长,这些节点相应的函数值为)(,),(1m x y x y ". 在点n x 处,已知))(,()('n n n x y x f x y =用n x 的向前差商nn n n x x x y x y −−++11)()(近似代替)('n x y ,如§1,则得到所谓的Euler 公式1(,)n n n n y y hf x y +=+——单步、显式格式11Euler 公式的局部截断误差:假设)(n n x y y =情况下,11)(++−n n y x y 称为局部截断误差.'''2311''23()()()()()2()(,()(()))2n n n n n n n n n y x y x y y x hy x h O h y x h y x f x y x h O h ++−=+++−−=+故有)(2)(''211n n n x y h y x y ≈−++. 122-2 后退的Euler 公式同样对常微分方程初值问题,在1+n x 点,已知))(,()(111'+++=n n n x y x f x y ,如果用向后差商hx y x y n n )()(1−+代替)(1'+n x y ,则得到后退的Euler 公式:111(,)n n n n y y hf x y +++=+——单步、隐式格式13相对于以上可以直接计算1+n y 的Euler 公式(显式),上式是隐式公式. 一般来讲,显式容易计算,而隐式具有更好的稳定性.求解上述公式,通常使用迭代法:对于给定的初值)0(1+n y,计算(1)()111(,)(0,1,)k k n n n n y y f x y k ++++=+=", 如果)(1lim k n k y +∞→收敛,则其极限必满足上述后退Euler 公式.14局部截断误差:假设)(n n x y y =,则),()(111++++=n n n n y x hf x y y .由于)]()[,())(,(),(1111111+++++++−+=n n n y n n n n x y y x f x y x f y x f η且''''2111(,())()()()()n n n n n f x y x y x y x hy x O h +++==++15则有'2''31111(,)[()]()()()()n y n n n n n n y hf x y y x y x hy x h y x O h η++++=−++++将此式减去式2'''31()()()()()2n n n n h y x y x hy x y x O h +=+++ 可得,2''311111()(,)[()]()()2n n y n n n n h y x y hf x y x y y x O h η+++++−=−−+16考虑到21111(,)()1(,)y n y n hf x O h hf x ηη++=++−,则有22''3''11()()()()22n n n n h h y x y y x O h y x ++−=−+≈−172-3 梯形公式由于上述两个公式的局部截断误差绝对值相等,符号相反,故求其算术平均得到梯形公式:111[(,)(,)]2n n n n n n hy y f x y f x y +++=++——单步、隐式格式18梯形法同样是隐式公式,可用下列迭代公式求解:(0)1(1)()111(,)[(,)(,)]2n n n n k k n n n n n n y y hf x y h y y f x y f x y +++++⎧=+⎪⎨=++⎪⎩局部截断误差:类似于后退Euler ,可计算出)(12)('''311n n n x y h y x y −≈−++192-4 改进的Euler 公式上述用迭代法求解梯形公式虽然提高了精度,但计算量也很大. 实际上常采用的方法是,用Euler 公式求得初始值(预测),然后迭代法仅施行一次(校正)——改进的Euler 公式:1111(,)[(,)(,)]2n n n n n n n n n n y y f x y hy y f x y f x y ++++⎧=+⎪⎨=++⎪⎩20估计上式中第二式当1+n y 为准确值时的局部截断误差:''11113(3)()()(()[()()])2()12n n n n n n n hy x y y x y x y x y x hy x ++++−=−++≈−212-5 Euler 两步公式如果用中心差商hx y x y n n 2)()(11−+−代替)('n x y ,则得Euler 两步公式112(,)n n n n y y hf x y +−=+——两步、显式格式22假设1−n y 及n y 均为准确值,利用Taylor 展式容易计算Euler 两步公式的局部截断误差为:11113(3)()()(()2(,()))()3n n n n n n n y x y y x y x hf x y x h y x +++−−=−+≈23此式与梯形公式相结合,得到如下的预测-校正公式:111112(,)[(,)(,)]2n n n n n n n n n n y y hf x y hy y f x y f x y −++++⎧=+⎪⎨=++⎪⎩假设第一式中的1−n y 及n y ,以及第二式中的n y 及1+n y 均是准确值,则有,2441)()(1111−≈−−++++n n n n y x y y x y 从而可得以下的事后估计式,111111114()()51()()5n n n n n n n n y x y y y y x y y y ++++++++⎧−≈−−⎪⎪⎨⎪−≈−⎪⎩25可以期望,以上式估计的误差作为计算结果的补偿,可以提高计算精度.以n p 及n c 分别表示第n 步的预测值和校正值,则有以下的“预测-改进-校正-改进”方案(其中在1+n p 与1+n c 尚未计算出来的前提下,以n n c p −代替11++−n n c p :26预测:'112n n n hy y p +=−+预测的改进:)(5411n n n n c p p m −−=++计算:),(11'1+++=n n n m x f m校正:)(2'1'1++++=n n n n m y hy c校正的改进:)(511111++++−+=n n n n c p c y计算:),(11'1+++=n n n y x f y27例 用Euler 方法求解初值问题2'[0,0.6](0)1y y xy x y ⎧=−−∈⎨=⎩取0.2h =,要求保留六位小数. 解:Euler 迭代格式为2210.2()0.80.2k k k k k k k k y y y x y y x y +=+−−=−因此2821000(0.2)0.80.20.8y y y x y ≈=−= 22111(0.4)0.80.20.6144y y y x y ≈=−=23222(0.6)0.80.20.461321y y y x y ≈=−=29例 用改进的Euler 方法求解初值问题2'sin 0[0,0.6](0)1y y y x x y ⎧++=∈⎨=⎩取0.2h =,求(0.2),(0.4)y y 的近似值,要求保留六位小数.解:改进的Euler 格式为212211110.2(sin )0.2(sin sin )2k k k k k k k k k k k k k y y y y x y y y y x y y x +++++⎧=+−−⎪⎨=+−−−−⎪⎩30即,222110.820.08sin 0.1(0.80.2sin )sin k k k k k k k k y y y x y y x x ++=−−−则有1(0.2)0.807285y y ≈=,2(0.4)0.636650y y ≈=31§3 Runge -Kutta 方法Def.1如果一种方法的局部截断误差为)(1+p h O ,则称该方法具有p 阶精度. 323-2 Runge —Kutta 方法的基本思想上述的Taylor 级数法虽然可得到较高精度的近似公式,但计算导数比较麻烦. 这里介绍不用计算导数的方法.))(,()()()('1h x y h x f h x y hx y x y n n n n n θθθ++=+=−+——平均斜率.33如果粗略地以),(n n y x f 作为平均斜率,则得Euler 公式;如果以221K K +作为平均斜率,其中),(1n n y x f K =,),(112hK y x f K n n +=+,则得改进的Euler 公式.343-3 二阶的Runge -Kutta 方法对点n x 和)10(≤<+=+p ph x x n p n ,用这两点斜率的线性组合近似代替平均斜率,则得计算公式:11122121()(,)(,)n n n n n p n y y h K K K f x y K f x y phK λλ++⎧=++⎪=⎨⎪=+⎩35现确定系数p ,,21λλ,使得公式具有二阶精度. 因为,取n y 为()n y x ,则'1(,)(,())'()n n n n n nK f x y f x y x y x y === 再把2K 在),(n n y x 处展开,有36'21(,)(,)n p n n n n K f x y phK f x ph y phy +=+=++代入可得,'2''31122()()n n n n y y hy ph y O h λλλ+=++++'2(,)(,)(,)()n n x n n y n n n f x y f x y ph f x y phy O h =+⋅+⋅+'2(')(,)()n x y n n y ph f f y x y O h =+⋅+⋅+'''2()n n y ph y O h =+⋅+37相比较二阶Taylor 展开''2'12n n n n y h hy y y ++=+,有,⎪⎩⎪⎨⎧==+211221p λλλ满足此条件的公式称为二阶Runge -Kutta 公式.38可以验证改进的Euler 公式属于二阶Runge -Kutta 公式. 下列变形的Euler 公式也是二阶Runge -Kutta 公式:12121(,)(,)22n n n n n n y y hK K f x y h h K f x y K +⎧⎪=+⎪=⎨⎪⎪=++⎩393-4 三阶Runge -Kutta 公式同二阶Runge -Kutta 公式,考虑三点,,(01)n n p n q x x x p q ++≤≤≤试图用它们的斜率321,,K K K 的线性组合近似代替平均斜率,即有如下形式的公式:1112233121312()(,)(,)(,())n n n n n n n n y y h K K K K f x y K f x ph y phK K f x qh y qh rK sK λλλ+=+++⎧⎪=⎪⎨=++⎪⎪=+++⎩40把32,K K 在),(n n y x 处展开,通过与)(1+n x y 在n x 的直接Taylor 展式比较,可确定系数s r q p ,,,,,,321λλλ,满足下式,从而使得上述公式具有三阶精度,41特别地,2,1,1,21,32,61231=−======s r q p λλλ是其一特例.123232223311213161p q p q pqs r s λλλλλλλλ++=⎧⎪⎪+=⎪⎪⎪+=⎨⎪⎪=⎪⎪+=⎪⎩423-5 四阶Runge -Kutta 公式相同的方法,可以导出下列经典的四阶Runge -Kutta 公式:112341213243(22)6(,)(,)22(,)22(,)n n n n n n n n n n h y y K K K K K f x y h h K f x y K h h K f x y K K f x h y hK +⎧=++++⎪⎪=⎪⎪⎪=++⎨⎪⎪=++⎪⎪=++⎪⎩43例 用经典四阶Runge —Kutta 方法求解初值问题'83[0,0.4](0)1y y x y =−⎧∈⎨=⎩,取0.2h =,求(0.4)y 的近似值,要求保留六位小数.解:四阶Runge —Kutta 格式为44112341211123122241330.2(22)6(,)830.2(,)83(0.1) 5.6 2.120.2(,)83(0.1) 6.32 2.372(,0.2)83(0.2) 4.208 1.578k k k k k k k k k k k kk k k k ky y K K K K K f x y y K f x y K y K yK f x y K y K y K f x y K y K y ++++⎧=++++⎪⎪==−⎪⎪⎪=+=−+=−⎨⎪⎪=+=−+=−⎪⎪⎪=+=−+=−⎩则10.5494 1.2016k k y y +=+,45故12(0.2) 2.3004,(0.4) 2.4654y y y y ≈=≈=.注:由准确解382()33xy x e −=−可得(0.2) 2.300792,(0.4) 2.465871y y ==46§5 线性多步法基本思想:在计算1+i y 之前,已计算出一系列的近似值i y y ,,1",如果充分利用这些已知信息,可以期望会获得更高精度的)(1+i x y 的近似值1+i y .基本方法:基于数值积分与基于Taylor 展开的构造方法.475-1 基于数值积分的构造方法对方程),('y x f y =两边从i x 到1+i x 积分,则得∫++=+1),()()(1i ix x i i dxy x f x y x y 设)(x P r 是f (x , y )的插值多项式,由此可得以下的一般形式的计算公式:∫++=+1)(1i ix x r i i dxx P y y 48例 取线性插值))(,())(,()(11111+++++−−+−−=i i i i ii i i i i r x y x f x x x x x y x f x x x x x P ,则得到梯形法:)],(),([2111+++++=i i i i i i y x f y x f hy y495-2 Adams 显式公式在区间],[1+i i x x 上利用r +1个数据点),(,),,(),,(11r i r i i i i i f x f x f x −−−−"构造插值多项式)(x P r ,由牛顿后插公式(注意到:j i j i j f f −Δ=∇)j i jrj j i r f j t th x P −=Δ⎟⎟⎠⎞⎜⎜⎝⎛−−=+∑0)1()(其中!)1()1(j j s s s j s +−−=⎟⎟⎠⎞⎜⎜⎝⎛". 50可得10rj i i rj i jj y y h f αΔ+−==+∑——Adams 显式公式其中1(1)j j t dt j α−⎛⎞=−⎜⎟⎝⎠∫,它可写成:∑=−++=rj ji rj i i f h y y 01β515-3 Adams 隐式公式在区间],[1+i i x x 上利用r +1个数据点),(,),,(),,(1111+−+−++r i r i i i i i f x f x f x "构造插值多项式)(x P r ,由牛顿后插公式101)1()(+−=+Δ⎟⎟⎠⎞⎜⎜⎝⎛−−=+∑j i jrj ji r f j t th x P 可得*11rj i i rj i j j y y h f α+−+==+Δ∑——Adams 隐式公式52其中01(1)jj t dt j −−⎛⎞α=−⎜⎟⎝⎠∫,它又可写成: *11ri i rj i j j y y h f β+−+==+∑535-4 Adams 预测-校正公式以r =3时的Adams 显式与隐式公式为例. 此时,显式公式为)9375955(243211−−−+−+−+=i i i i i i f f f f hy y 利用Taylor 展式,容易计算局部截断误差为)(720251)5(5i x y h . 54)5199(242111−−+++−++=i i i i i i f f f f hy y 同样利用Taylor 展开可得,其局部截断误差为5(5)19()720i h y x −. 隐式公式为55⎪⎩⎪⎨⎧+−++=−+−+=−−+++−−−+)519),(9(24)9375955(24211113211i i i i i i i i i i i i i f f f y x f hy y f f f f h y y 注 利用2-5节的相同作法同样可以构造更精确的计算过程.可构造利用显式预测,隐式校正的计算公式:56§6 方程组与高阶方程的情形6-1 一阶方程组常微分方程初值问题为⎩⎨⎧==00)(),('y x y y x f y 此时T m y y y ),,(1"=,Tm f f f ),,(1"=. 此时上述的一切方法均可使用,只是注意y 与f 此时为向量.576-2 化高阶方程为一阶方程组解下列的m 阶方程()(1)'(1)(1)000000(,,',,)(),'(),,()m m m m y f x y y y y x y y x y yx y −−−⎧=⎨===⎩""令)1(21,,',−===m m y y y y y y ",则有58'12'23'1'12(,,,,)m m m m y y y y y yy f x y y y −⎧=⎪=⎪⎪⎨⎪=⎪⎪=⎩#"初始条件为:)1(00'002001)(,,)(,)(−===m m y x y y x y y x y "。
常微分方程的难点
常微分方程是数学分析中的一门重要课程,也是应用数学中的基础课程之一。
它是研究一阶或高阶导数与自变量关系的方程,涉及到函数的连续性、可微性、可积性等重要的数学概念。
然而,常微分方程的学习也是有难点的。
其中,常见的难点包括以下几个方面:
1. 初值问题和边值问题的区别和联系。
初值问题和边值问题是常微分方程的两种基本类型,它们的解法和理论基础都有所不同。
2. 高阶常微分方程的解法。
高阶常微分方程的解法需要掌握多种技巧和方法,如常数变易法、欧拉公式、拉普拉斯变换等。
3. 变量分离法和分步法的应用。
变量分离法和分步法是解常微分方程中常用的技巧,但其应用需要考虑到方程的特殊性质和形式。
4. 非线性常微分方程的解法。
非线性常微分方程的解法涉及到多种数学工具和方法,如相似变量、对称性、积分因子等,需要掌握较高的数学知识。
5. 常微分方程的应用。
常微分方程是应用数学中的重要工具,在物理、工程、生物等领域都有着广泛的应用。
但其应用需要考虑到实际问题的特殊性质和背景知识。
- 1 -。
微分方程是数学中非常重要的一类方程,广泛应用于物理、工程、经济等领域。
在求解微分方程时,我们通常需要找到其特解,以满足特定的初始条件或边界条件。
本文将介绍微分方程中的特解求法。
微分方程是包含未知函数及其导数的方程。
常见的微分方程包括常微分方程和偏微分方程。
常微分方程是只包含一元未知函数及其导数的方程,而偏微分方程则包含多元未知函数及其偏导数。
求解微分方程的一般方法是通过积分和代数操作来消去导数项,最终得到等式中只含有未知函数及其低阶导数的方程。
我们先从常微分方程中的特解求法入手。
常微分方程的特解分为两类:初值问题和边值问题。
初值问题是在微分方程中给定某一点的函数值以及其导数值,求解满足这些条件的特解。
边值问题是在微分方程中给定一些边界条件,求解满足这些条件的特解。
对于初值问题,我们可以使用分离变量、变量代换、常数变易法等方法来求解微分方程。
以一阶常微分方程为例,假设方程为dy/dx=f(x),同时给定y(x0)=y0的初始条件。
我们可以首先将方程两边同时乘以dx,然后再两边积分,得到∫dy=f(x)dx。
再对两边同时求定积分,即可得到特解y=f(x)+C,其中C为常数。
代入初始条件,得到具体的特解。
对于边值问题,我们可以使用分离变量、变量代换、极值原理等方法来求解微分方程。
以二阶常微分方程为例,假设方程为d²y/dx²=g(x),同时给定y(x0)=y0和y(x1)=y1的边界条件。
我们可以首先将方程两边同时乘以dx²,然后再两边积分两次,得到∫∫d²y=g(x)dx²。
然后,我们可以使用边界条件来确定积分常数,从而得到特解。
对于偏微分方程的特解求法稍有不同。
由于偏微分方程包含多元未知函数及其偏导数,所以特解形式较复杂。
在求解偏微分方程的特解时,我们常常需要使用提升技巧、变量分离法、特征方程法等方法。
这些方法的具体步骤较为复杂,在此不一一赘述。
总之,微分方程中的特解求法是数学中重要且复杂的问题。
常微分方程知识点常微分方程是微积分的一个重要分支,是描述物理、生物、经济等各类现象的一种数学模型。
常微分方程描述了未知函数与其导数之间的关系,在实际问题中具有广泛的应用。
下面将介绍常微分方程的基本概念、解的存在唯一性、一阶常微分方程和高阶常微分方程等知识点。
1.基本概念:常微分方程描述的是函数与其导数之间的关系。
常微分方程可以分为初值问题和边值问题。
初值问题是给定了函数在特定点的初始值和导数,要求求解函数在整个定义域上的表达式;边值问题是给定了函数在两个点的值,要求求解函数在这两个点之间的表达式。
2.解的存在唯一性:对于一阶常微分方程的初值问题,如果方程的右端函数在整个定义域上连续且满足利普希茨条件,那么方程存在唯一解。
其中利普希茨条件是指有一个正数L,使得对于任意t和s,满足,f(t)-f(s),≤L,t-s。
3.一阶常微分方程:一阶常微分方程描述的是未知函数y与其一阶导数y'之间的关系。
一阶常微分方程的一般形式为dy/dt = f(t, y),其中f(t, y)是已知函数。
一阶常微分方程的解可以通过分离变量、线性方程、齐次方程和恰当方程等方法求解。
4.高阶常微分方程:高阶常微分方程描述的是未知函数与其高阶导数之间的关系。
高阶常微分方程的一般形式为d^n y/dt^n = F(t, y, y', ..., y^n-1),其中F(t, y, y', ..., y^n-1)是已知函数。
高阶常微分方程的解可以通过代数法、特征方程和待定系数法等方法求解。
5.变量分离方法:当一阶常微分方程的右端可以写成g(y)·h(t)的形式时,可以使用变量分离方法求解。
将方程改写为1/g(y) dy = h(t) dt,然后对两边分别积分得到∫1/g(y) dy = ∫h(t) dt,从而求得y的表达式。
6.线性方程方法:当一阶常微分方程可以写成y'+p(t)y=q(t)的形式时,可以使用线性方程方法求解。
边界值问题的定义及求解方法边界值问题(Boundary Value Problem,简称BVP)是数学中经典问题之一,它被广泛应用于各种科学和工程领域的模型分析和数值计算中。
本文将为您介绍边界值问题的定义、求解方法以及应用实例。
一、边界值问题的定义边界值问题是一类微分方程求解问题,它要求在某个区域内已知微分方程的解,以及在区域边界上给出解的初值或者边界值条件,求解微分方程在整个区域内的解。
边界值问题一般分为两种:Dirichlet问题和Neumann问题。
Dirichlet问题即在区域边界上给出解的值,而Neumann问题则是在区域边界上给出解的导数值。
二、边界值问题的求解方法1. 差分法差分法是一种常见的数值解法,它利用微分方程的一阶或者高阶差分逼近微分算子,将微分方程转化为代数方程组。
然后采用迭代或者直接求解代数方程组的办法得到微分方程的解。
2. 有限元法有限元法是一种求解偏微分方程的数值计算方法,它使用有限维函数空间来逼近实际问题的解。
将区域分割成若干个单元,建立有限元函数空间,然后根据偏微分方程和边界条件构造代数方程组,最后采用数值计算方法求解。
3. 辛普森法辛普森法是一种求解积分的数值方法,利用区间端点、抛物线顶点和中点构成的近似抛物线来逼近被积函数,从而得到积分的近似值。
三、边界值问题的应用实例1. 电路问题电路问题是一种常见的边界值问题,求解电路问题可以将电路看作一个带有边界条件的微分方程模型。
通过差分法或者有限元法求解该微分方程,可以得到电路中电流、电压等物理量的数值解。
2. 热传导问题热传导问题是一种边界值问题,它描述了物体中的温度分布问题。
通过差分法或者有限元法求解该方程,可以得到物体中温度的分布以及热流分布,为物体的热力学分析提供了重要的数值计算方法。
3. 声波传播问题声波传播问题也是一种边界值问题,它描述了声波在介质中的传播。
通过有限元法求解该方程,可以得到声波的传播路径以及声压分布,为声学分析提供了重要的数值计算方法。
常微分方程第三版课本概述“常微分方程第三版课本”是一本由X编写的教材,主要介绍了常微分方程的基本概念、理论和解析方法。
本教材内容丰富、结构清晰,适用于高等院校的常微分方程课程教学,也可以作为自学的参考资料。
目录1.基本概念– 1.1 常微分方程的定义– 1.2 解的定义及存在唯一性定理– 1.3 初值问题和边值问题2.一阶常微分方程– 2.1 可分离变量方程– 2.2 齐次线性方程– 2.3 一阶线性常微分方程– 2.4 Bernoulli 方程和 Riccati 方程– 2.5 可降阶的高阶微分方程3.高阶线性常微分方程– 3.1 高阶常微分方程的一般理论– 3.2 同解、通解和特解– 3.3 常系数齐次线性方程– 3.4 常系数非齐次线性方程及其特解– 3.5 变系数线性方程4.线性常微分方程组– 4.1 二阶齐次线性方程组– 4.2 二阶非齐次线性方程组和线性常系数方程组– 4.3 三阶及三阶以上线性方程组内容简介基本概念本章介绍了常微分方程的基本概念,包括常微分方程的定义、解的定义及存在唯一性定理、初值问题和边值问题。
通过对这些概念的学习,读者可以对常微分方程有一个基本的认识。
一阶常微分方程本章主要介绍了一阶常微分方程的解析方法,包括可分离变量方程、齐次线性方程、一阶线性常微分方程、Bernoulli方程和 Riccati 方程、可降阶的高阶微分方程等。
通过对这些解析方法的学习,读者可以熟练地解决一阶常微分方程的问题。
高阶线性常微分方程本章主要介绍了高阶线性常微分方程的理论和方法。
包括高阶常微分方程的一般理论、同解、通解和特解、常系数齐次线性方程、常系数非齐次线性方程及其特解、变系数线性方程等。
通过对这些理论和方法的学习,读者可以掌握高阶线性常微分方程的解法。
线性常微分方程组本章主要介绍了线性常微分方程组的解法。
包括二阶齐次线性方程组、二阶非齐次线性方程组和线性常系数方程组、三阶及三阶以上线性方程组等。
-179-第十五章 常微分方程的解法建立微分方程只是解决问题的第一步,通常需要求出方程的解来说明实际现象,并加以检验。
如果能得到解析形式的解固然是便于分析和应用的,但是我们知道,只有线性常系数微分方程,并且自由项是某些特殊类型的函数时,才可以肯定得到这样的解,而绝大多数变系数方程、非线性方程都是所谓“解不出来”的,即使看起来非常简单的方程如22x y dxdy+=,于是对于用微分方程解决实际问题来说,数值解法就是一个十分重要的手段。
§1 常微分方程的离散化下面主要讨论一阶常微分方程的初值问题,其一般形式是⎪⎩⎪⎨⎧=≤≤=0)(),(y a y bx a y x f dxdy(1)在下面的讨论中,我们总假定函数),(y x f 连续,且关于y 满足李普希兹(Lipschitz)条件,即存在常数L ,使得|||),(),(|y y L y x f y x f -≤-这样,由常微分方程理论知,初值问题(1)的解必定存在唯一。
所谓数值解法,就是求问题(1)的解)(x y 在若干点b x x x x a N =<<<<= 210处的近似值),,2,1(N n y n =的方法,),,2,1(N n y n =称为问题(1)的数值解,n n n x x h -=+1称为由n x 到1+n x 的步长。
今后如无特别说明,我们总取步长为常量h 。
建立数值解法,首先要将微分方程离散化,一般采用以下几种方法: (i )用差商近似导数若用向前差商hx y x y n n )()(1-+代替)('n x y 代入(1)中的微分方程,则得),1,0())(,()()(1 =≈-+n x y x f hx y x y n n n n化简得))(,()()(1n n n n x y x hf x y x y +≈+如果用)(n x y 的近似值n y 代入上式右端,所得结果作为)(1+n x y 的近似值,记为1+n y ,则有),1,0(),(1 =+=+n y x hf y y n n n n (2)这样,问题(1)的近似解可通过求解下述问题 ⎩⎨⎧==+=+)(),1,0(),(01a y y n y x hf y y n n n n (3)得到,按式(3)由初值0y 可逐次算出 ,,21y y 。
本科科研训练论文常微分方程的边值问题学生姓名:郭骏学号:**********专业:数学及其应用数学年级:08级学院:理学院【摘要】边值问题是微分方程问题的一个类型。
在求解微分方程时,除了给出方程本身,往往还需给出一定的定解条件。
最常见的是初值问题,即给出的定解条件为初始条件;但也有一些情况,定解条件要考虑所讨论区域的边界,如在一个区间讨论时,定解条件在区间的两个端点给出,这种定解条件称为边界条件,相应的定解问题称为边值问题。
边值问题的提出和发展,与流体力学,材料力学,波动力学等密切相关;并且在现代控制理论等学科中有重要应用。
【关键词】常微分方程边值问题研究目录第一章引言1.1常微分方程的起源和发展1.2常微分方程的内容1.3常微分方程的应用1.4 常微分方程的实例第二章常微分方程边值问题的研究2.1 边值问题的提出2.2 二阶线性常微分方程边值问题的可解性2.3 特征值问题参考文献第一章引言1.1 常微分方程的起源微分方程研究的来源:它的研究来源极广,历史久远。
I.牛顿和G.W.莱布尼茨创造微分和积分运算时,指出了它们的互逆性,事实上这是解决了最简单的微分方程y┡=ƒ(x)的求解问题。
当人们用微积分学去研究几何学、力学、物理学所提出的问题时,微分方程就大量地涌现出来。
20世纪以来,随着大量的边缘科学诸如电磁流体学、化学流体力学、动力气象学、海洋动力学、等等的产生和发展,也出现不少新型的微分方程(特别是方程组)。
70年代随着数学向化学和生物学的渗透,又出现了大量的反应扩散方程。
常微分方程在我国的发展中华人民共和国建立后,微分方程得到了重视和发展。
培养了许多优秀的微分方程的工作者,在常微分方程稳定性、极限环、结构稳定性等方面做出了很多有水平的结果;在偏微分方程混合型刻画渗流问题的拟线性退缩抛物型、椭圆组和拟线性双曲组的间断解等方面做出了很多有水平的结果。
1.2 常微分方程的内容定义1 凡含有参数,未知函数和未知函数导数 (或微分) 的方程,称为微分方程,有时简称为方程,未知函数是一元函数的微分方程称作常微分方程,未知数是多元方程的微分方程称作偏微分方程.微分方程中出现的未知函数最高阶导数的阶数,称为微分方程的阶.定义式如下:F(x, y, y¢, ...., y(n)) = 0 。
定义2 任何代入微分方程后使其成为恒等式的函数,都叫做该方程的解.若微分方程的解中含有任意常数的个数与方程的阶数相同,且任意常数之间不能合并,则称此解为该方程的通解(或一般解).当通解中的各任意常数都取特定值时所得到的解,称为方程的特解。
1.3 常微分方程的应用现在,常微分方程在很多学科领域内有着重要的应用,自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等。
这些问题都可以化为求常微分方程的解,或者化为研究解的性质的问题。
应该说,应用常微分方程理论已经取得了很大的成就,但是,它的现有理论也还远远不能满足需要,还有待于进一步的发展,使这门学科的理论更加完善。
1.4 常微分方程的实例下列方程都是微分方程 (其中 y, v, q 均为未知函数):(1) y= kx, k 为常数;(2) ( y - 2xy) dx + x2 dy = 0;(3) mv (t) = mg - kv(t);第一章 常微分方程边值问题的研究2.1 边值问题的提出求常微分方程满足给定边界条件的解的问题。
亦即,设常微分方程为对区间I 上的点α1,α2,…,αk 及值y (αi ),y ┡(αi ),…,y (n-1)(αi )(i=1,2,…,k ,k >1),给定了一些条件,求此方程在 I 上的满足这些条件的解的问题。
这些条件称为边界条件,诸αi 及y (αi )、y ┡(αi )、…、y (n-1)(αi ) 称为边值或边界值。
当k =2,α1、α2是区间I 的端点时,称为两点边值问题。
边值问题的提出和发展,与流体力学、材料力学、波动力学以及核物理学等密切相关;并且在现代控制理论等学科中有重要应用。
因为常微分方程可以解析求解的类型甚少,所以求边值问题的解也是困难的。
为了适应实际问题的需求,不得不采用近似解法,这样,首先需要回答:边值问题的解是否存在?是否惟一?这就是边值问题的基本论题。
2.2 二阶线性微分方程边值问题的可解性我们将重点研究试射法。
二阶常微分方程一般可表示成如下的形式:),,()(y y x f x y '='', b x a ≤≤ (2.1)边值条件有如下三类:第一类边值条件α=)(a y , β=)(b y (2.2)第二类边值条件α=')(a y , β=')(b y (2.3)第三类边值条件[19]ααα='-)()(10a y a y , βββ='+)()(10b y b y (2.4)其中010≥αα, 010≥ββ, 010≠+αα, 010≠+ββ。
在对边值问题用数值方法求解之前,应该从理论上分析该边值问题的解是否存在,若问题的解不存在,用数值方法计算出来的数据没有任何意义。
下面的定理给出了边值问题存在唯一解的充分条件。
定理 设方程(2.1)中的函数f 及y f ∂∂,y f '∂∂在区域},,|),,{(∞<'<-∞≤≤'=Ωy y b x a y y x内连续,并且 (ⅰ),0),,(>∂'∂yy y x f Ω∈'∀),,(y y x ; (ⅱ)y y y x f '∂'∂),,(在Ω内有界,即存在常数M ,使得 M y y y x f ≤'∂'∂),,(, Ω∈'∀),,(y y x , 则边值问题(2.1)-(2.4)的解存在且唯一。
(1)试射法为了描述试射法的基本思想,先给出初值问题的概念。
由(2.5)中的二阶常微分方程以及初始条件α=)(a y , s a y =')( (2.6)构成的定解问题⎩⎨⎧='=≤≤=+'+''-s a y a y b x a x r y x q y x p y )(,)(,)()()(α (2.7)称之为初值问题。
对于边值问题(2.5)的求解,“试射法”的基本思想是将边值问题转化成初值问题来求解,即根据边界条件(2.2),寻求与之等价的初始条件(2.6),也就是说,反复调整初始时刻的斜率值s y =',使得初值问题(2.7)的积分曲线)(x y y =能“命中”β=)(b y .设能够提供s 的两个预测值1s ,2s ,我们按这两个斜率“试射”, 通过求解相应的初值问题(2.7)可以得到)(b y 的两个预测值分别为1β,2β。
若1β和2β都不满足预定的精度,则可用线性插值的方法校正1s ,2s 得到新的斜率值)(1121213ββββ---+=s s s s (2.8) 然后再按斜率3s 试射,求解相应的初值问题(2.7)又得到新的结果3)(β=b y .若ββ=3或εββ<-3,则可将3s 作为s 的近似值;否则,继续过程(2.8)直到找到计算结果)(b y 与β相当符合为止。
基于叠加原理的试射法设二阶线性常微分方程边值问题(2.5)的解存在并且唯一,并定义线性算子L :y x q y x p y Ly )()(:+'+''-=. (2.9)我们考虑如下的两个线性微分方程的初值问题:⎩⎨⎧==,)(),(αa u x r Lu 0)('=≤≤a u b x a (2.10) 和⎩⎨⎧==,0)(,0a v Lv 1)('=≤≤a v b x a (2.11) 设)(x u 和)(x v 分别为问题(2.10)和(2.11)的解,不难验证)()()()()(x v b v b u x u x y -+=β (2.12)是问题(2.5)的解,其中0)(≠b v . 通过上述描述,我们可以得到基于叠加原理的打靶法的基本步骤为:1. 根据边值问题(2.5)构造相应的初值问题(2.10)和(2.11);2. 分别求出两个初值问题(2.10)和(2.11)的解)(x u 和)(x v ;3. 将)(x u 和)(x v 按(2.12)式做组合,所得的函数)(x y 就是边值问题(2.5)的解.(2.10)和(2.11)均为二阶常微分方程初值问题,求解时可通过引入变量代换将其化成相应的一阶方程组初值问题。
如令:,,11v v u u == v v u u '='=22 (2.13) 则(2.10)式可以写成 ⎪⎩⎪⎨⎧==-+='=',0)(,)(),()()(,2112221a u a u x r u x q u x p u u u α (2.14) (2.11)式可以写成⎪⎩⎪⎨⎧==+='=',1)(,0)(,)()(,2112221a v a v v x q v x p v v v (2.15) 这样就可以利用Runge-Kutta 方法求解(2.14)和(2.15)。
对于更一般的线性边值问题:⎪⎩⎪⎨⎧≠+≥='+≥='-≤≤=+'+''-≡0,0,)()(,0,)()(,)()()(0010101010βαβββββαααααb y b y a y a y b x a x r y x q y x p y Ly (2.16) 用基于叠加原理的打靶法的步骤为:1. 根据(2.16)式,构造两个相应的初值问题:⎩⎨⎧-==,)(),(1αc a u x r Lu α0)(c a u b x a -='≤≤ (2.17) 和⎩⎨⎧==,)(,01αa v Lv 0)(α='≤≤a v b x a (2.18) 其中0c 和1c 是满足条件10110=-ααc c 的两个任意的常量.2. 求解初值问题(2.17)和(2.18)式,设其解分别为)(x u 和)(x v .3. 将)(x u 和)(x v 做线性组合)()()()]()([)()(1010x v b v b v b u b u x u x y '+'+-+=βββββ 由此计算得到的函数)(x y 就是(2.16)式的解。
(2)有限差分法将区间[a,b]进行等分:h=(b-a)/(n+1), x i =a+ih,i=0,1,…,n+1,设在x=x i ,i=0,1,…,n+1处得数值解为y i 。
用中心差分近似微分,即{y i ′≈y i+1−y i−12ℎy i "≈y i+1−2y i +y i−1h 2,i=0,1,…,n则离散化成差分方程 y i+1−2y i +y i−1=h 2f(x i ,y i ,y i+1−y i−12h ),i=0,1,…,n对应的边界条件也离散成ẏ(a)-α0y(a)=α1,ẏ(b)+β0y(b)=β1;第一类边界问题:y 0=α,y n+1=β第二类边界问题:y 1-y 0=h α,y n+1-y n−1=h β第三类边界问题:y 1-(1+α0h)y 0=α1h,(1+β0h )y n+1-y n =β1h若f(x,y,y ’)是y ,y ’的线性函数时,f 可以写成f(x,y,y ’)=p(x)y ’(x)+q(x)y(x)+r(x)其中p(x),q(x),r(x)为已知函数,则由常微分方程的理论知,通过变量替换总是可以消去方程中的y ’项,不妨假设变换后的方程为y ”(x)-q(x)y(x)=r(x);y(a)=α,y (b )=β则近似差分方程成离散差分方程为y i+1−2y i +y i−1h 2-q i y i =r i ;y 0=α,y n+1=β其中q i=q(x i) ,r i=r(x i),i=0,1,…,n 将以上方程合并同类型整理得方程组{y0=αy i−1y n+1=β−(2+q i h2)y i+y i+1=r i h2其中只要q i≫0,则方程组的系数矩阵为弱对角占优的三对角阵,方程组为三对角线方程组,可以用追赶法求解。