第六章_常微分方程初值问题的数值解法_习题课

数值分析习题参考解答 江世宏编

1 第一章 绪论

姓名 学号 班级

习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1 若误差限为5105.0,那么近似数0.003400有几位有效数字？（有效数字的计算）

解：2*103400.0x，325*10211021xx

故具有3位有效数字。

2 14159.3具有4位有效数字的近似值是多少？（有效数字的计算）

解：10314159.0，欲使其近似值*具有4位有效数字，必需

41*1021，3*310211021，即14209.314109.3*

即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。

3 已知2031.1a，978.0b是经过四舍五入后得到的近似值，问ba，ba有几位有效数字？（有效数字的计算）

解：3*1021aa，2*1021bb，而1811.2ba，1766.1ba

2123****102110211021)()(bbaababa

故ba至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(bbaaabbaab故ba至少具有2位有效数字。

4 设0x，x的相对误差为，求xln的误差和相对误差？（误差的计算）

解：已知**xxx，则误差为 ***lnlnxxxxx

则相对误差为 ******lnln1lnlnlnxxxxxxxx

5测得某圆柱体高度h的值为cmh20*，底面半径r的值为cmr5*，已知cmhh2.0||*，cmrr1.0||*，求圆柱体体积hrv2的绝对误差限与相对误差限。（误差限的计算）

微分方程初值问题练习题求解微分方程的初值问题

微分方程是数学中的一个重要分支，它描述了变量之间的关系及其变化率。初值问题是指在给定一个微分方程及初始条件的情况下，求解出一个特定的解。本文将通过练习题的形式，来介绍如何求解微分方程的初值问题。

1. 练习一：一阶线性常微分方程

考虑以下一阶线性常微分方程：

\[ \frac{dy}{dx} + P(x)y = Q(x) \]

其中，\(P(x)\) 和 \(Q(x)\) 是给定的函数。已知初值条件 \(y(x_0) =

y_0\)，求解出该微分方程的解。

解答：

首先将原方程变形为标准形式：

\[ \frac{dy}{dx} = -P(x)y + Q(x) \]

接下来使用积分因子法来求解该微分方程，积分因子定义为：

\[ \mu(x) = e^{\int -P(x) dx} \]

对原方程两边同时乘以积分因子，得到：

\[ \mu(x) \frac{dy}{dx} + \mu(x)P(x)y = \mu(x)Q(x) \]

由于左边是积分的导数，可以写成： \[ \frac{d}{dx}(\mu(x)y) = \mu(x)Q(x) \]

对上式两边同时积分，得到：

\[ \int \frac{d}{dx}(\mu(x)y) dx = \int \mu(x)Q(x) dx \]

应用积分的基本性质，化简上式得到：

\[ \mu(x)y = \int \mu(x)Q(x) dx + C \]

其中，\(C\) 是常数。

最后将 \(y\) 解出来，得到：

\[ y(x) = e^{-\int P(x) dx}(\int e^{\int P(x) dx}Q(x) dx + C) \]

将初值条件 \(y(x_0) = y_0\) 代入上式，可以求解出常数 \(C\) 的值，从而得到特定的解。

2. 练习二：二阶线性常微分方程

考虑以下二阶线性常微分方程：

常微分方程中的初值问题

一、介绍

初值问题是在微积分学中一个非常基础的概念，在常微分方程（ODEs）中也有很重要的应用。我们从初值问题开始，逐步深入探讨ODEs的相关知识。

二、什么是初值问题？

在ODEs的求解中，我们通常需要给出一个初值条件，也就是某个时刻的初始条件。通常我们把这个条件称之为初值问题（Initial Value Problem, IVP）。

例如，我们可以假设现在有一个物体在运动。如果我们想要得到它在任意时间点上的位置和速度，就需要知道它在某个时刻的位置和速度，这个时刻就称为初值。

三、ODEs的解与初值问题

ODEs的求解通常与初值问题密切相关。在求解ODEs时，我们通常需要设定初值条件，从而得到方程的一组解。

举个例子来说，如果一个物体在力的作用下做匀加速运动，那么我们可以得到ODEs如下：

$\frac{d^2x}{dt^2}=a$

这里，x表示物体的位移，t代表时间，a代表加速度。我们可以通过对此方程积分，得到如下解：

$x(t)=\frac{1}{2}at^2+C_1t+C_2$

其中，C1和C2都是常数，需要通过初值条件来确定。

假设我们知道在t=0时，这个物体的位移为 $x_0$ ，速度为

$v_0$ 。那么我们就可以得到初始条件：

$x(0)=x_0,C_2=x_0$

$\frac{dx}{dt}(0)=v_0,C_1=v_0$

通过这两个初始条件，我们就可以得到这个物体在任意时刻的位移和速度。

四、初值问题的数值求解

除了解析求解以外，初值问题在实际工程中还有很多数值求解的方法。

在给出数值解之前，首先需要对微分方程进行离散化。一种简单的离散化方式是欧拉法。

对于ODEs：$\frac{dy}{dt}=f(t,y)$

我们可以将它离散化为：

$\frac{y_{i+1}-y_i}{h}=f(t_i,y_i)$

其中，h是离散化的步长，i表示当前离散点的下标。这个式子可以帮助我们递推地求出 $y_{i+1}$ 的值。

第六章常微分方程的数值解法

第六章 常微分方程的数值解法

在自然科学研究和工程技术领域中，常常会遇到常微分方程的求解问题。传统的数学分析方法仅能给出一些简单的、常系数的、经典的线性方程的解析表达式，不能处理复杂的、变系数的、非线性方程，对于这些方面的问题，只能求诸于近似解法和数值解法。而且在许多实际问题中，确确实实并不总是需要精确的解析解，往往只需获得近似的解或者解在若干个点上的数值即可。在高等数学课程中介绍过的级数解法和逐步逼近法，能够给出解的近似表达式，这一类方法称为近似解法。还有一类方法是通过计算机来求解微分方程的数值解，给出解在一些离散点上的近似值，这一类方法称作为数值方法。

本章主要介绍常微分方程初值问题的数值解法，包括Euler 方法、Runge-Kutta 方法、线性多步法以及微分方程组与高阶微分方程的数值解法。同时，对于求解常微分方程的边值问题中比较常用的打靶法与有限差分法作了一个简单的介绍。

§1 基本概念

1.1 常微分方程初值问题的一般提法

常微分方程初值问题的一般提法是求解满足如下条件的函数，，b

x a x y ≤≤)(

=<<=α

)(),(a y b

x a y x f dx

dy

， (1.1) 其中),(y x f 是已知函数，α是给定的数值。

通常假定上面所给出的函数),(y x f 在给定的区域},),{(+∞<≤≤=y

b x a y x D 上面满足如下条件：

(1) 函数),(y x f 在区域D 上面连续；

(2) 函数),(y x f 在区域D 上关于变量y 满足Lipschitz(李普希茨)条件： 212121,),(),(y y b x a y y L y x f y x f ?≤≤?≤?，， (1.2)

其中常数L 称为Lipschitz(李普希茨)常数。

由常微分方程的基本理论可以知道，假如(1.1)中的),(y x f 满足上面两个条件，则常微分方程初值问题(1.1)对于任意给定的初始值α都存在着唯一的解，，b x a x y ≤≤)(并且该唯一解在区间[a,b]上是连续可微的。