下载文档原格式
常微分方程初值问题的解法及应用常微分方程是数学中非常重要的一部分，它涉及了许多领域的模型建立和问题求解。
本文将介绍常微分方程初值问题的解法及其应用。
一、常微分方程初值问题的定义常微分方程初值问题是指给定一个常微分方程，以及它在某一点上的初始条件，求解该方程的解曲线。
通常，一个常微分方程初值问题可以表示为:y'(x) = f(x,y), y(x0) = y0,其中，y(x)是未知函数，f(x,y)是已知函数，y(x0) = y0是初始条件。
二、常微分方程初值问题的解法常微分方程初值问题的解法有多种，下面我们将介绍几种常用的方法。
1.欧拉法欧拉法是最简单的一种求解常微分方程初值问题的方法。
该方法基于初始条件，通过不断迭代计算得到近似解曲线。
具体步骤如下：步骤1：设定步长h，确定求解区间[x0, xn]，计算步数n。
步骤2：初始化，即确定初始点(x0, y0)。
步骤3：根据方程dy/dx = f(x,y)和初始点(x0, y0)，计算斜率k = f(x0, y0)。
步骤4：根据已知的斜率和步长h，计算下一个点的坐标(xi+1,yi+1)。
步骤5：重复步骤3和步骤4，直到达到步数n。
步骤6：得到近似解曲线。
2.改进的欧拉法（改进欧拉法）改进的欧拉法是对欧拉法的改进，其求解精度比欧拉法更高。
具体步骤如下：步骤1：设定步长h，确定求解区间[x0, xn]，计算步数n。
步骤2：初始化，即确定初始点(x0, y0)。
步骤3：根据方程dy/dx = f(x,y)和初始点(x0, y0)，计算斜率k1 =f(x0, y0)。
步骤4：根据已知的斜率k1和步长h/2，计算中间点的坐标(x0+h/2, y0+k1*h/2)。
步骤5：根据方程dy/dx = f(x,y)和中间点的坐标(x0+h/2, y0+k1*h/2)，计算斜率k2= f(x0+h/2, y0+k1*h/2)。
步骤6：根据已知的斜率k2和步长h，计算下一个点的坐标(xi+1,yi+1)。
浅谈常微分方程初值问题数值解法在自然科学、工程技术、甚至社会科学的一些领域中，常常会遇见一阶常微分方程的求解问题：()上述问题，寻求解的具体表达式十分困难，仅对一些特殊形式的才有可能找到解的解析表达式，在大多情况下，初值问题的解不能用初等函数表示出来即使可写出解的解析表达式，但因为这些表达式过于复杂，要计算它在某些点上的函数值也异常困难。
在实际问题中，经常需要的恰是解在某些点上的函数值，因此研究初值问题的数值解法十分必要。
1 常微分方程初值问题的数值解法常微分方程的近似解法大体可分成三大类：一类是图解法和器械法；第二类是解的近似法；第三类是数值解法，即通过离散化的方法直接求出函数在某些点上的近似值，此数值解仅为精确解的近似解。
其基本原理为：一阶常微分方程的初值问题的解是上变量的连续函数，因此求上述问题的数值解，就是在区间上的若干离散点上用离散化的方法将初值问题化成离散变量的相应问题，从而相应问题的解可作为初值问题理论解的近似值。
由常微分方程的理论可知，只要在区域内连续，且关于满足林普希兹条件，则方程的解存在且唯一。
初值问题的数值解法通常采取“步进法”，而“步进法”又可分为“单步法”和“多步法”两类。
（1）单步法。
所谓“单步法”是指在计算时，只用到前一步的有关信息。
其一般形式为：，主要包括下面三种方法：Euler方法，改进的Euler公式-梯形公式和Runge-Kutta法。
（2）线性多步法。
单步法没有用到前几步计算得到的信息，因此为了提高精度，需重新计算多个点处的函数数值，如RK方法，故计算量较大。
线性多步法的基本思想是充分利用前面的已知信息来构造精度高且计算量小的算法来计算。
多步法常用方法是线性多步法，求解公式为：构造的常用方法是Taylor展开和数值积分方法。
常用的线性多步公式有：四阶Adams显式公式：四阶Adams隐式公式：四阶Milne显式公式：三阶Hamming公式：（隐式公式）预测校正系统和预测校正修正法：一般地，同阶的隐式法比显式法精确，而且数值稳定性好，但隐式公式中的求解较难，需要用到迭代法，这就增加了计算量。
数值分析复习题及答案数值分析复习题及答案 Company number：【WTUT-WT88Y-W8BBGB-BWYTT-19998】数值分析复习题⼀、选择题1. 和分别作为π的近似数具有（）和（）位有效数字.A ．4和3B ．3和2C ．3和4D ．4和42. 已知求积公式()()211211()(2)636f x dx f Af f ≈++?，则A ＝（）A ． 16B ．13C ．12D ．233. 通过点()()0011,,,x y x y 的拉格朗⽇插值基函数()()01,l x l x 满⾜（）A ．()00l x ＝0，()110l x = B ．()00l x ＝0，()111l x =l x = D ．()00l x ＝1，()111l x =4. 设求⽅程()0f x =的根的⽜顿法收敛，则它具有（）敛速。
A ．超线性B ．平⽅C ．线性D ．三次5. ⽤列主元消元法解线性⽅程组1231231 220223332x x x x x x x x ++=??++=??--=?作第⼀次消元后得到的第3个⽅程（）.A ．232x x -+= B ．232 1.5 3.5x x -+= C ．2323x x -+=D ．230.5 1.5x x -=-⼆、填空1. 设2.3149541...x *=，取5位有效数字，则所得的近似值x= .2.设⼀阶差商()()()21122114,3=---，()()()322332615,422f x f x f x x x x --===-- 则⼆阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 则2||||X = ，=∞||||X 。
4．求⽅程2 1.250x x --= 的近似根，⽤迭代公式 1.25x x =+，取初始值 01x =，那么1______x =。
5．解初始值问题 00'(,)()y f x y y x y =??=?近似解的梯形公式是 1______k y +≈。
数值计算试题库----填空题（每小题3分）第一章1、数x *＝2.1972246···的六位有效数字的近似数的绝对误差限是。
2、取 3.142x =作为 3.141592654x =┅的近似值，则x 有位有效数字.3、已知96112168.≈有五位有效数字，则方程01262=+-x x 的具有五位有效数字的较小根为。
4、3*x 的相对误差约是*x 的相对误差的_____ 倍5、为了提高数值计算精度, 当正数x 充分大时, 应将)1ln(2--x x 改写为_______.6、. 已知数 e=2.718281828...,取近似值 x=2.7182,那麽x 具有的有效数字是＿＿＿位。
7、设 3149541.2*=x ，取5位有效数字，则所得的近似值=x _____.8、数值计算方法中需要考虑的误差为。
9、计算4.12,)12(6≈-=取f ，利用算式6)12(1+，3)223(-，3)223(1+，27099-计算，得到的结果最好的算式为。
10、sin1有2位有效数字的近似值840.的相对误差限是第二章11、已知函数()f x 的函数值()()()()()0,2,3,5,6f f f f f ，以及均差如下()()()()()00,0,24,0,2,35,0,2,3,51,0,2,3,5,60f f f f f =====那么由这些数据构造的牛顿插值多项式的最高次幂的系数是12、满足()a a f x x =，()b b f x x =，()c c f x x =的拉格朗日插值余项为。
13、二阶均差f (x 0,x 1, x 2) = _________________________________． 14、设1)(3-+=x x x f ，则差商[]3 ,2 ,1 ,0f =__________.15、通过四个互异节点的插值多项式p (x ),只要满足＿＿＿＿＿＿＿，则p (x )是不超过二次的多项式。
微分方程定解问题的基本概念微分方程是数学中的一个重要分支，它用来描述物理、经济、生物等学科中的现象和问题。
微分方程定解问题则是微分方程研究的重点，它对于解决实际问题具有非常重要的作用。
一、微分方程的基本概念微分方程是描述变量之间的变化关系的方程，其形式通常为：y′ = f(x, y)其中y′ 表示 y 对 x 的导数，f(x, y) 表示 x 和 y 的函数关系。
微分方程的解是一组函数，它满足微分方程和附加条件（称为初值条件或边界条件）。
二、定解问题的基本概念定解问题是指在微分方程中确定初始条件或边界条件，求得微分方程的解。
定解问题可以分为初值问题和边值问题。
初值问题是在一个点（通常为 x0）给出一个函数值（通常为y(x0)）和其导数值（通常为y′(x0)），求解函数在另一点的取值。
初值问题通常用初值问题解法求解。
边值问题是在一段区间内给出一个函数值和其导数值，求解函数在该区间的取值。
边值问题通常用曲线拟合法或数值法求解。
三、常见的定解问题常见的定解问题包括：1.一阶常微分方程的初值问题。
例如：y′ = f(x, y), y(x0) = y02.一阶常微分方程的边值问题。
例如：y′ = f(x, y), y(a) = ya, y(b) = yb3.二阶常微分方程的初值问题。
例如：y′′ = f(x, y, y′), y(x0) = y0, y′(x0) = y0′4.二阶常微分方程的边值问题。
例如：y′′ = f(x, y, y′), y(a) = ya, y(b) = yb四、定解问题的应用定解问题在物理、工程、金融等领域中有广泛的应用。
例如：1.物理学中的定解问题：在自然界中的各种物理现象中，微分方程定解问题经常被用于对各种现象和性质的研究和分析。
2.工程学中的定解问题：设计和分析各种工程系统时，微分方程定解问题经常被用于模型的建立和计算。
3.金融领域中的定解问题：在金融领域中，微分方程定解问题被用来分析各种金融产品的产生和变化，预测市场走势等。
二单项选择题1. 已知近似值1x ,2x ，则()12,x x ()=A. ()()2112x x x x + B 。
()()12x x +C. ()()1122x x x x + D 。
()()12x x2. 已知求积公式()()211211()(2)636f x dx f Af f ≈++⎰，则A ＝（ ） A ． 16 B 。
13 C 。
12 D. 233. 已知2112A ⎡⎤=⎢⎥⎣⎦，则化为A 为对角阵的平面旋转变换角θ＝（ ) A ．6π B 。
4π C 。
3π D. 2π 4. 设求方程()0f x =的根的切线法收敛，则它具有( ）敛速。
A ． 线性 B. 超越性 C 。
平方 D 。
三次5。
改进欧拉法的局部截断误差为( ）A ． ()5O h B. ()4O h C. ()3O h D 。
()2O h填空题1。
π的近似值3.1428是准确到 近似值。
2. 满足()a a f x x =，()b b f x x =,()c c f x x =的拉格朗日插值余项为 。
3。
用列主元法解方程组时，已知第2列主元为()142a 则()142a ＝ 。
4．乘幂法师求实方阵 的一种迭代方法。
5. 欧拉法的绝对稳定实区间为 。
计算题1. 用已知函数表求抛物插值多项式,并求1()2f 的近似值。
2. 用紧凑格式解方程组123410114130141x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦3. 已知方程组123210113110121x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦（1） 证明高斯－塞德尔法收敛；（2） 写出高斯－塞德尔法迭代公式; （3） 取初始值()()00,0,0TX=,求出()1X 。
4. 用4n =复化辛卜公式计算积分1011dx x +⎰，并估计误差。
5. 用一般迭代法求方程[]0,0.5内的根。
