数值分析中的误差.
数值分析中的名词解释
数值分析中的名词解释数值分析是一门研究如何利用计算机进行数值计算和模拟的学科,它在科学计算、工程领域以及许多其他领域中都有广泛的应用。
本文将通过解释数值分析中的一些重要名词,来介绍这个领域的基本概念和方法。
一、误差与精度在数值分析中,误差是指数值计算和实际结果之间的差异。
由于计算过程中存在舍入误差、截断误差等,数值计算很难得到完全准确的结果。
为了度量误差的大小,我们需要引入精度的概念。
精度表示了计算结果的准确程度,通常使用绝对误差或相对误差来衡量。
绝对误差是计算结果与实际结果的差值,而相对误差则是绝对误差与实际结果的比值。
二、插值与外推插值是指根据已知数据点的数值,通过某种方法去估算出未知点的数值。
常用的插值方法有拉格朗日插值、牛顿插值等。
而外推则是利用已知数据点的数值,通过推算来估计未知点的数值。
插值和外推在数值分析中常常用于构建函数的近似表达式或预测未来数据的趋势。
三、数值积分与数值微分数值积分是指通过数值方法来近似求解定积分。
由于很多函数的原函数无法用解析算式表示,或者求解困难,因此数值积分成为了一种常用的求解方法。
常见的数值积分方法有梯形法则、辛普森法则等。
而数值微分则是通过数值方法来近似求解微分。
数值微分的目的是通过逼近导数的定义来估算导数值,通常使用数值差商或有限差分来实现。
四、线性方程组的解法在科学计算中,线性方程组的求解是一个核心问题。
数值分析中有各种不同的算法和方法可以用来解决线性方程组,如高斯消元法、追赶法、迭代法等。
这些方法的基本思想是通过对系数矩阵进行操作或迭代运算来求解未知数的值。
线性方程组的求解在很多科学和工程问题中都非常重要,比如力学模拟、电路分析等。
五、常微分方程的数值解法常微分方程是描述自然界中许多现象的数学模型。
然而,绝大部分的常微分方程都无法用解析算式求解,因此需要使用数值方法来近似求解。
数值分析中有许多不同的方法可以用于求解常微分方程,如欧拉法、龙格-库塔法、四阶龙格-库塔法等。
数值分析中的误差分析与收敛性
数值分析中的误差分析与收敛性数值分析是一门研究利用计算机进行数学计算和问题求解的学科,它在科学计算、工程设计、金融分析等领域中具有广泛的应用。
然而,在数值计算过程中,由于计算机的有限精度和数值算法的近似性质,误差问题成为了一个不可避免的挑战。
因此,了解误差的来源和性质,以及数值计算方法的收敛性,对于保证计算结果的准确性和可靠性非常重要。
本文将探讨数值分析中的误差分析与收敛性问题。
1. 误差的来源及分类在数值计算中,误差可以分为四类:舍入误差、截断误差、模型误差和舍入误差。
舍入误差是由于计算机内部使用有限位数表示实数导致的误差,它来源于将实数近似为计算机可表示的数值。
截断误差是在计算过程中采取舍入法或截断法将无限级数或无限小量等进行有限近似所引入的误差。
模型误差是将实际问题用数学模型进行近似所引入的误差,它包括了模型的简化和不完全描述等因素。
舍入误差是由于使用有限位数存储和运算导致的误差。
2. 误差的度量方法误差的度量方法包括绝对误差和相对误差。
绝对误差是指数值近似解与真实解之间的差值,它可以用来度量数值计算的准确度。
相对误差是绝对误差除以真实解的绝对值后得到的比值,它可以用来度量数值计算的相对准确度。
通过对误差进行度量和分析,可以评估数值计算方法的准确性,并选择合适的数值方法来解决实际问题。
3. 收敛性在数值计算中,所谓的收敛性是指数值方法的逼近解序列以某种方式趋近于真实解。
一个数值方法是收敛的,意味着当步长趋于0时,逼近解趋近于真实解。
收敛性的评估是数值分析中一个重要的问题,它关系到数值方法的稳定性和可靠性。
常见的收敛性分析方法包括局部截断误差、阶、收敛速度等。
局部截断误差是用来评估数值方法在每个步长上的近似误差,阶是用来度量数值方法逼近真实解的速度。
4. 提高数值计算的准确性与可靠性为了提高数值计算的准确性与可靠性,我们可以采取多种方法。
首先,选择合适的数值方法和算法,确保其满足问题的数学性质和准确性要求。
数值分析误差
I k −1
11 ( k = n, n − 1,…,2,1) = − Ik 5k
(1 − 3)
依式( 依式(1-3)计算
* 0
的近似值。 I n −1 , I n − 2 ,…, I 1 , I o 的近似值。
* 14
1 1 1 分别取 I = 0.18232155, I = + ≈ 0.01222222 2 6 × 15 5 × 15 按算法1、算法 2的计算结果见下屏表 1 − 1:
逆向递推公式在数学上完全等价,却导致两种完全不同的 逆向递推公式在数学上完全等价, 算法。对于实数序列的递推由于初始误差的存在,可以一 算法。对于实数序列的递推由于初始误差的存在,
种方向的递推会使误差扩大, 种方向的递推会使误差扩大,而另一方向的递推会使得误 差逐步减小。在设计(选用) 差逐步减小。在设计(选用)算法时要用使初始误差不增 长的算法。 长的算法。
1 3 1 5 作近似计算, 取 S = x − x + x ,作近似计算,则 3! 5! 为其截断误差。 为其截断误差。
R = sin x − S
条 件 问 题
计算方法中有一类问题称为条件问题, 公式) 条件问题是一个算法 (公式)由于初始 数据或者中间某些数据微小摄动对计算结 果产生影响的敏感性的问题。舍入误差、 果产生影响的敏感性的问题。舍入误差、 观测误差都属初始数据的摄动。研究坏条 观测误差都属初始数据的摄动。 件问题的计算方法是十分重要的课题, 件问题的计算方法是十分重要的课题,有 的时候,一些问题的条件并不坏, 的时候,一些问题的条件并不坏,但由于 算法不恰当, 算法不恰当,初始数据的微小摄动或舍入 误差在计算过程中不断被放大,而可能导 误差在计算过程中不断被放大, 致计算结果的精度大大降低, 致计算结果的精度大大降低,甚至使计算 失去意义。
数值分析误差限的计算公式
数值分析误差限的计算公式1、误差x∗为 x 一个近似值绝对误差:e∗=x∗−x相对误差:e∗r=e∗x=x∗−xx,由于真值 x 总是不知道的,通常取e∗r=e∗x∗=x∗−xx∗误差限:|x∗−x|≤ε∗相对误差限:ε∗r=ε∗|x∗|ε(f(x∗))≈|f′(x∗)|ε(x∗)2、插值法记ωn+1(x)=(x−x0)(x−x1)⋯(x−xn)Lagrange 插值多项式系数:lk(xk)=(x−x0)⋯(x−xk−1)(x−xk+1)⋯(x−xn)(xk−x0)⋯(xk−xk−1)(x −xk+1)⋯(x−xn)Lagrange 插值多项式:Ln(x)=∑k=0nlk(x)yk=∑k=0nykωn+1(x)ω′n+1(xk)(x−xk) 余项:记 Mn+1=maxa≤x≤b|fn+1(x)|R(x)=fn+1(ξ)ωn+1(x)(n+1)!≤Mn+1(n+1)!|ωn+1(x)|均差与 NewTon 插值多项式一阶均差:f[x0,xk]=f(xk)−f(x0)xk−x0k 阶均差:f[x0,x1,⋯,xk]=f[x0,⋯,xk−2,xk]−f[x0,⋯,xk−2,xk−1]xk−xk−1f[x0,x1,⋯,xn]=f(n)(ξ)n!(x0,x1,⋯,xn,ξ∈[a,b])f[x0,x1,⋯,xk]=∑j=0kf(xj)ω′k+1(xj)NewTon 插值多项式:Pn(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+⋯+f[x0,x1,⋯,xn](x−x0)(x−x1)⋯(x−xn−1)余项:R(x)=f[x0,x1,⋯,xn]ωn+1(x)Hermite 插值Taylor 多项式:Pn(x)=f(x0)+f′(x0)(x−x0)+⋯+f(n)(x0)n!(x−x0)n余项:R(x)=f(n+1)(ξ)(n+1)!(x−x0)n+1若已知 f(x0),f′(x1),f(x1),f(x2):P(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+A(x−x0)(x−x1)(x−x2)其中 A 由 P′(x1)=f′(x1) 可得余项:R(x)=14!f(4)(ξ)(x−x0)(x−x1)2(x−x2)两点三次 Hermite 插值多项式:H3(x)=αk(x)yk+αk+1(x)yk+1+βk(x)mk+βk+1(x)mk+1其中 mk=f′(xk),mk+1=f′(xk+1)⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧αk(x)=(1+2x−xkxk+1−xk)(x−xk+1xk−xk+1)2αk+1(x)=(1+2x−xk+1xk−xk+1)(x−xkxk+1−xk)2⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧⎧βk(x)=(x−xk)(x−xk+1xk−xk+1)2βk+1(x)=(x−xk+1)(x−xkxk+1−xk)2余项:R(x)=f(4)(ξ)4!(x−xk)2(x−xk+1)2分段低次插值h=b−an对每个小区间使用对应插值公式求 Ih(x)余项对分段线性插值函数:maxa≤x≤b|f(x)−Ih(x)|≤M28h2对分段三次埃尔米特插值:maxa≤x≤b|f(x)−Ih(x)|≤M4384h43、数值积分代数精度定义:如果某个求积公式对于次数不超过 m 的多项式均能够准确成立,但对于 m+1 次多项式就不准确成立,则称该公式具有 m 次代数精度梯形公式公式与中矩形公式梯形公式:∫baf(x)dx≈b−a2f(a)+b−a2f(b)余项:R[f]=−(b−a)312f′′(η)(η∈(a,b))矩形公式:∫baf(x)dx≈(b−a)f(a+b2)余项:R[f]=(b−a)324f′′(η)(η∈(a,b))Newton-Cotes 公式将积分区间 [a,b] 分成 n 等分Simpson 公式(n=2):∫baf(x)dx≈b−a6f(a)+b−a6f(b)+2(b−a)3f(a+b2)余项:R[f]=−(b−a)5180∗24f(4)(η)(η∈(a,b))Cotes 公式(n=4):C=b−a90[7f(x0)+32f(x1)+12f(x2)+32f(x3)+7f(x4)]余项:R[f]=−2(b−a)7945∗46f(6)(η)(η∈(a,b))复合求积公式积分区间 [a,b] 分成 n 等分,步长 h=b−an复合梯形公式:Tn=h2[f(a)+2∑k=0n−1f(xk)+f(b)]余项:Rn(f)=−b−a12h2f′′(η)复合 Simpson 求积公式:Sn=h6[f(a)+2∑k=0n−1f(xk)+4∑k=1n−2f(x(k+1)/2)+f(b)] 其中 x(k+1)/2=xk+h2Rn(f)=−b−a180(h2)4f(4)(η)龙贝格求积算法T(0)0=h2[f(a)+f(b)]求梯形值 T0(b−a2k),利用递推公式求 T(k)0,递推公式:T2n=12Tn+h2∑k=0n−1f(xk+12)求加速值:T(k)m=4m4m−1Tk+1m−1−14m−1T(k)m−1k=1,2,⋯高斯-勒让德求积公式积分区间为 [−1,1]∫1−1f(x)dx≈∑k=0nAkf(xk)余项:n=1 时,R1[f]=1135f(4)(η)4、解线性方程组的直接方法列主元高斯消去法在每次消元时,选取列主元在最前面,列主元为该列最大值矩阵三角分解法如果 n 阶矩阵 A 的各阶顺序主子式 Dk(k=1,2,⋯,n−1) 均不为零,则必有单位下三角矩阵 L 和上三角矩阵 U,使得 A=LU,并且 L 和 U 是唯一的。
数值分析中的误差估计理论
数值分析中的误差估计理论数值分析是研究通过数值计算方法来解决数学问题的学科。
在数值计算过程中,由于计算机本身的限制以及数值计算方法的局限性,必然会引入一定的误差。
误差估计理论是数值分析中的重要内容,它的主要任务是评估数值计算结果的准确性,并为我们提供合理的结果判断依据。
一、误差类型在进行误差估计之前,我们首先需要了解误差的分类。
在数值计算中,误差可以分为截断误差和舍入误差两种类型。
1. 截断误差:截断误差是由于数值计算方法的有限步骤导致的近似解与准确解之间的差距。
通常情况下,我们使用有限级数或多项式来近似某个函数,但是由于级数或多项式只能截取有限的项数,从而无法精确地表示原函数,所以会引入截断误差。
2. 舍入误差:舍入误差是由于计算机在存储和表示数值时的有限精度所引起的误差。
计算机只能存储有限位数的数字,而且在计算过程中会进行舍入操作,从而导致精确数字的丢失和近似数字的产生。
二、误差估计的方法误差估计的方法主要包括局部误差估计和全局误差估计两种。
1. 局部误差估计:局部误差估计方法是通过分析数值计算方法的近似性质,对每一步计算过程的误差进行估计。
通常情况下,我们会使用泰勒级数展开来近似求解函数值,然后通过对级数剩余项的估计来获得局部误差的上界。
2. 全局误差估计:全局误差估计方法是通过分析数值计算方法的整体性质,对整个计算过程的误差进行估计。
该方法通常使用数值稳定性定理或者收敛速度分析来评估数值计算的精度,从而给出全局误差的上界。
三、误差控制策略在数值计算中,确保误差控制是非常重要的。
误差控制策略通过采用合适的数值计算方法和调整计算过程的步骤,减小误差并控制误差的传播,从而提高结果的准确性。
1. 精确算法选择:在进行数值计算之前,我们需要评估不同数值计算方法的精确性和稳定性,并选择适合的方法。
合适的数值计算方法可以最大程度地减小误差的产生。
2. 步长控制:对于迭代算法或差分方法,我们可以通过控制步长的大小来控制误差。
数值分析 误差知识与算法知识
一、误差的来源与分类 二、 绝对误差、相对误差与有效数字
三、误差估计的基本方法
四、算法的计算复杂性 五、数值运算中的一些原则
1.2误差知识与算法知识
一、误差的来源与分类 模型误差 (描述误差 ) ( 测量误差) (方法误差 ) ( 计算误差 )
观测误差
截断误差 舍入误差
建模过程中 产生的误差
三、误差估计的基本方法 (一)误差估计的一般运算 一元函数:
e( f (a)) f (a) e(a)
二元函数:
( f (a)) f (a) (a)
f (a, b) f (a, b) e( f (a, b)) e(a) e(b) x y
f (a, b) f (a, b) ( f (a, b)) ( a) (b) x y
Tn an 秦九韶算法 Tk xTk 1 ak , k n 1, n 2,,1,0 p ( x) T 0 n
加法次数: n
n(n 1) 乘法次数: 2
pn ( x) a0 x(a1 x(a2 x(an1 xan ) )
有效数字=可靠数字+存疑数字
(3)有效数字 有效数字的定义: 设a是x的近似值,如果a的误差绝对值不超过x 的第k位小数的半个单位,即
则称近似值a准确到小数点后第k位。 从这个小数点后第k位数字直到最左边非零数 字之间的所有数字都叫有效数字。
1 k x a 10 2
1 1 2 (2.18) 10 (2.1200) 10 4 2 2
例8 设有三个近似数
a=2.31, b=1.93, c=2.24 它们都有三位有效数字,试计算 p a bc, ( p), r ( p), 并问:p的计算结果能有几位有效数字? 教材例4
数值分析中的误差
第9章 数值分析中的误差 典型问题解析考试知识点:误差、有效数字。
(6%)学习要点:误差、有效数字。
典型问题解析:一、误差绝对误差e :e =x -x *(设精确值x *的近似值x , 差e =x -x *称为近似值x 的绝对误差(误差))。
绝对误差限ε:ε≤-=*x x e(绝对误差限ε是绝对误差e 绝对值的一个上界。
)相对误差e r :***-==x x x x e e r (绝对误差e 与精确值x *的比值,常用x e e r =计算) 相对误差限r ε:r r e ε≤(相对误差e r 绝对值的一个上界),r r x x x x e εε=≤-=||||||***,*xr εε=,常用x ε计算. 绝对误差限的估计式:(四则运算中))()()(2121x x x x εεε+=± )()()(122121x x x x x x εεε+≈22122121+=x x x x x x x )()()(εεε 二、有效数字有效数字:如果近似值x 的误差限ε 是它某一个数位的半个单位,我们就说x 准确到该位. 从这一位起到前面第一个非0数字为止的所有数字称为x 的有效数字.(1)设精确值x *的近似值x ,若m n a a a x 10.021⨯±=a 1,a 2,…,a n 是0~9之中的自然数,且a 1≠0,n l x x l m ≤≤110⨯50=≤--,.*ε 则x 有l 位有效数字.例1 设x *= π=3.1415926…,若x *的近似值x 为3.14,3.1415,3.143,求x 的有效数字位数.解:若x =3.14=0.314×101,(m =1)31105.06592001.0-*⨯≤=- x x (l =3)故x =3.14有3位有效数字。
若x =3.1415=0.31415×101,(m =1)41105.00000926.0-*⨯≤=- x x (l =4)故x =3.1415有4位有效数字。
数值分析中的复化梯形法误差分析
数值分析中的复化梯形法误差分析数值分析中的复化梯形法误差分析在数值分析中,复化梯形法是一种常用的数值积分方法。
它使用梯形规则进行近似求解定积分,通过将定积分区间分割成若干个小区间,并在每个小区间上使用梯形规则进行求解,最后将各个小区间上的积分结果相加得到整个定积分的近似值。
本文将对复化梯形法进行误差分析。
1. 复化梯形法原理复化梯形法的原理是将定积分区间[a, b]等分为n个小区间,令h=(b-a)/n为小区间长度,梯形法的近似结果T可以表示为:T = h/2 * (f(a) + 2*f(x1) + 2*f(x2) + ... + 2*f(x(n-1)) + f(b))其中,f(x)为被积函数在x点处的取值。
2. 复化梯形法误差分析复化梯形法的误差主要包括局部误差和全局误差。
2.1 局部误差在每个小区间上,我们使用梯形规则进行积分计算,其误差可以通过泰勒展开进行推导。
设f(x)在[a, b]区间上具有充分高阶连续导数,则对于每个小区间[xk, x(k+1)],我们有如下局部误差公式:E_local = - (h^3/12) * f''(ξ)其中,ξ为[xk, x(k+1)]上的某点,f''(ξ)为f(x)的二阶导数在ξ点的取值。
2.2 全局误差全局误差是指整个区间[a, b]上的积分近似与真实积分之差。
复化梯形法的全局误差可以通过对各个小区间上的局部误差进行累加得到。
假设积分的真实值为I,则全局误差E_global可以表示为:E_global = (b-a) * (h^2/12) * f''(ξ)其中,ξ为[a, b]区间上的某点,f''(ξ)为f(x)的二阶导数在ξ点的取值。
3. 误差分析实例为了更好地理解复化梯形法的误差特点,我们以一个具体的例子进行分析。
考虑定积分∫(0, 1)sin(x)dx的近似求解,将积分区间等分为4个小区间进行计算。
数值分析中的误差传播理论
数值分析中的误差传播理论
误差传播理论是数值分析领域中的重要理论之一,它用于研究数值
计算中的误差如何随着计算过程的进行而逐步累积或减小。
在实际的
数值计算中,由于浮点数表示精度有限、截断误差、舍入误差等原因,误差是无法完全避免的。
因此,了解误差传播理论对于提高数值计算
的准确性和稳定性至关重要。
误差传播理论主要包括截断误差和舍入误差两个方面。
截断误差是
由于数值计算中采用截断近似导致的误差,舍入误差是由于计算机在
表示实数时采用有限精度浮点数表示法引入的误差。
在数值计算过程中,这两种误差会相互影响,相互传播,从而影响最终计算结果的准
确性。
当进行多步数值计算时,每一步的计算结果都会带有误差。
这些误
差在后续计算中会被传播,并随着计算的进行而累积。
因此,为了减
小误差的传播,需要采取一些措施,比如适当选择计算方法、合理设
计算法、增加计算精度等。
误差传播理论对于数值计算中的算法设计、收敛性分析和稳定性评
价起着重要作用。
通过对误差传播过程的深入研究,可以帮助我们更
好地理解数值计算中的误差来源和传播规律,有助于提高数值计算方
法的精度和效率。
总之,误差传播理论是数值计算中不可或缺的理论基础,只有充分
理解和掌握误差传播规律,才能有效地提高数值计算的准确性和稳定
性,确保计算结果的可靠性。
希望通过对误差传播理论的深入研究和应用,能够为数值计算领域的发展和应用带来更多的启发和帮助。
数值分析中的误差分析
E ( x) = x − X
*
*
x*
| E ( x) |=| x − x* |<= η
此时,称为近似值的绝对误差限,简称误差限或精度
• 相对误差与相对误差限 E ( x) x − x* Er( x) = = 绝对误差与精度值之比,即称 x X * X 的相对误差.在实际中,由于精确值x一般无 为近似值 x − x* * 法知道,因此往往取 Er ( x) = 作为近似值的相对误差.
x*
类似于绝对误差的情况,若存在 δ >0 ,使得 x − x* * | Er ( x) |=| * |<= δ 则称 δ 为近似值 X 的相对误差限, x 相对误差是无量刚的数,通常用百分比表示,称为百分误 差.
• 有效大小,又能表示其精确程度,于是需要引 进有效数字的概念.再实际计算中,当准 确值x有很多位时,我们常按四舍五入得到 的近似值. |若近似值的绝对误差限
数值分析中的误差分析
误差与数值计算的误差估计
误差可以分为以下四种 • • • • 模型误差 观测误差 截断误差 舍如误差
误差与有效数字
• 绝对误差与绝对误差限 设某一量的精确值为x,其近似值为 X * ,则称 为近似值 X 的绝对误差,简称误差 当E(x)>0时,称为弱近似值或亏近似值,当E(x)<0时,称 X *为强近似值或盈近似值. 一般的,某一量的精确值x是不知道的,因而E(x)也无法求 出,但往往可以估计出E(x)的上界,即存在,使得
数值分析中的误差分析与收敛性
数值分析中的误差分析与收敛性数值分析是一门研究使用计算机进行数值计算的学科,它广泛应用于工程、科学和金融等领域。
在数值计算中,误差分析和收敛性是两个重要的概念。
本文将深入探讨数值分析中的误差分析和收敛性,并介绍它们的应用和意义。
一、误差分析在数值计算中,由于使用的是有限的计算机资源和近似的计算方法,无法得到完全准确的结果。
因此,误差分析成为一项必不可少的工作。
误差可以分为绝对误差和相对误差两种。
绝对误差是指数值计算的结果与真实值之间的差别,常用符号表示为Δx。
相对误差是指绝对误差与真实值之比,常用符号表示为εx。
绝对误差和相对误差可以通过以下公式计算:绝对误差:Δx = |x - x*|相对误差:εx = |(x - x*)/x*|其中,x表示近似值,x*表示真实值。
误差分析的目的是评估数值计算的精度和稳定性。
当误差较小且符合预期范围时,可以认为数值计算结果是可靠的。
二、收敛性在数值分析中,收敛性是指使用逼近方法得到的数值序列逐渐接近于准确值的性质。
收敛性分析是评估逼近方法有效性的重要手段。
常见的收敛性准则包括绝对收敛和相对收敛。
绝对收敛是指逼近序列的差值趋近于零,即对于任意给定的正数ε,存在正整数N,对于所有n>N,有|xn+1 - xn| < ε。
相对收敛是指逼近序列的比值趋近于一,即对于任意给定的正数ε,存在正整数N,对于所有n>N,有|(xn+1 -xn)/xn| < ε。
收敛性分析可以帮助我们评估数值计算方法的有效性和稳定性。
当逼近序列满足收敛准则时,可以认为该方法是可靠且收敛的。
否则,需要重新评估和改进计算方法。
三、误差分析与收敛性的应用误差分析和收敛性是数值分析中不可或缺的工具,其应用广泛且重要。
1. 误差分析在数值模拟中的应用数值模拟是利用数值方法来模拟和求解物理问题的过程。
在数值模拟中,误差分析可以帮助我们判断计算结果的可靠性,评估模拟的精度和稳定性。
通过分析误差来源和大小,可以优化计算方法,提高模拟结果的准确性。
《数值分析》第一章 数值计算中的误差
值,其绝对误差限等于该近似值末位的半个单位。
14
§2 舍入方法与有效数字
2.2 舍入方法
2.2.2四舍五入法
• 例:设a=-2.18和b=2.1200是分别由准确值x和y 经过四舍五入而得到的近似值,问: a、b的绝 对误差限、相对误差限各是多少?
解: (a) 0.005 0.5 102
(b) 0.00005 0.5104
n位
≤ 0 . 0 … 0 999... < 0 . 0 … 0 1=1×10-n
n位
n-1位
• 截断法产生的绝对误差限不超过近似数a最末位 的1个单位。
11
§2 舍入方法与有效数字
2.2 舍入方法
2.2.2四舍五入法
• 四舍情况,
A=a0 a1 … am . am+1 … am+n
• 当am+n+1 =0,1,2,3,4时,
4
§2 舍入方法与有效数字
5
§2 舍入方法与有效数字
2.1 绝对误差与相对误差
• 近似数a的绝对误差 , 简称误差 设a是精确值A的近似值,
=a-A
• 绝对误差限 ||=|a-A|<(上界)
• 由上式可推知 a- <A<a+,也可表示为A=aAFra biblioteka-a
a+
6
§2 舍入方法与有效数字
2.1 绝对误差与相对误差
• 相对误差 : 绝对误差与精确值之比 =/A。 • 实际计算/a。
代替后误差
a A 1 2
A a Aa
Aa
• 相对误差限 ||=|/a |< /|a|= (上界)
• 绝对误差是有量纲的量,相对误差没有量纲,有时 亦用百分比、千分比表示。
数值分析中的梯形法误差估计技巧
数值分析中的梯形法误差估计技巧梯形法(Trapezoidal Rule)是数值分析中用于近似计算定积分的一种方法,其基本原理是将被积函数的图像划分成若干个梯形,通过对这些梯形面积的求和来逼近定积分的值。
然而,梯形法并不是完全准确的,其近似值与真实值之间会存在一定的误差。
因此,精确地估计梯形法误差对于数值分析的准确性和可靠性具有重要意义。
本文将介绍数值分析中的梯形法误差估计技巧。
一、梯形法的基本原理梯形法是通过将被积函数的图像划分成一系列梯形来近似计算定积分的值。
给定被积函数 f(x),并将其在区间 [a, b] 上等距离地划分成 n 个小区间,记作 h = (b - a) / n,其中 n 为划分的小区间数。
梯形法的基本思想是,将每个小区间上的函数值视为一个梯形的高度并计算其面积,然后将所有梯形的面积求和,即可得到定积分的近似值。
具体而言,梯形法可以表示为如下的公式:∫[a,b] f(x)dx ≈ h/2 * [f(a) + 2f(a+h) + 2f(a+2h) + ... + 2f(a+(n-1)h) +f(a+nh)]二、梯形法误差估计的原理由于梯形法是一种近似计算方法,其得到的近似值与真实值之间会存在误差。
为了准确评估这种误差,我们需要引入梯形法误差估计公式。
根据数值分析中的原理,使用梯形法进行近似计算的误差估计公式为:其中 |E| 是梯形法近似值与真实值之间的误差绝对值,f''(ξ) 是被积函数 f(x) 在区间 [a, b] 上的二阶导数,ξ 是 [a, b] 中的某个点。
由此可见,梯形法的误差与步长 h 的平方成正比,也与被积函数的二阶导数有关。
因此,在实际应用中,我们可以通过减小步长 h 或者选择合适的划分区间来减小误差。
三、梯形法误差估计的应用示例为了更好地理解梯形法误差估计的应用,我们给出以下实际问题的例子:问题:使用梯形法计算函数 f(x) = x^3 在区间 [0, 2] 上的定积分,要求保留三位小数,并估计其误差。
数值分析中的误差分析方法
数值分析中的误差分析方法数值分析是一门研究离散数据逼近和连续函数求解的学科,广泛应用于科学、工程和金融等领域。
在数值计算过程中,误差是不可避免的,因此准确评估和分析误差是至关重要的。
本文将介绍数值分析中常用的误差分析方法,以帮助读者更好地理解误差来源和影响,从而提高数值计算的准确性和可靠性。
一、绝对误差和相对误差绝对误差是指数值计算结果与真实值之间的差异。
在数值分析中,我们往往无法得知真实值,因此无法直接计算绝对误差。
相对误差则是相对于近似值的误差,它可以更好地反映计算结果的准确性。
二、截断误差截断误差是由于采用有限的计算步骤或取舍了一些无限级数的项而引入的误差。
在数值计算中,我们通常使用近似方法,如级数展开和数值积分等。
由于截断误差的存在,我们得到的结果与真实值之间会有一定的差距。
截断误差的大小取决于所采用的数值方法和步长,可以通过逐步减小步长来减小截断误差。
三、舍入误差舍入误差是由于对无限精度数进行有限舍入导致的误差。
计算机中的数值表示是有限的,而真实数值通常是无限的。
因此,在计算机中进行数值计算时,会存在一定程度的舍入误差。
舍入误差可以通过采用更高精度的数据类型或者使用舍入误差分析技术来减小。
四、传播误差传播误差是由于输入数据的不确定性或测量误差在数值计算过程中扩散而引入的误差。
在实际问题中,输入数据通常带有不确定性,例如测量误差或近似值。
这些不确定性会随着计算的进行而传播,影响到计算结果的准确性。
传播误差需要通过敏感性分析等方法来进行评估和控制。
五、误差估计误差估计是通过数值分析方法来评估近似解与真实解之间的误差。
常用的误差估计方法包括残差估计、收敛性分析和算例分析等。
残差估计法通过计算数值解与原方程的残差来估计误差的大小。
收敛性分析则通过逐步减小步长和比较不同精度下的数值解来判断数值方法是否收敛。
算例分析是通过计算实际问题的已知解或近似解来评估数值方法的误差。
六、误差限制和误差控制误差限制和误差控制是保证数值计算结果准确性和可靠性的重要手段。
解释数值分析第五版中的截断误差和舍入误差。
解释数值分析第五版中的截断误差和舍入
误差。
截断误差和舍入误差是数值计算中很常见的两种误差。
截断误差
截断误差是指由于在计算中使用了有限数量的项或数字,而导致结果与精确结果之间的差异。
例如,在泰勒级数中,当我们只考虑其中有限的几项时,就会产生截断误差。
截断误差通常通过增加计算中使用的项数来减小。
然而,当项数变得非常大时,计算机的时间和内存可能无法处理。
舍入误差
舍入误差是指由于计算机在进行浮点数运算时必须对结果进行舍入而导致的误差。
计算机使用二进制来表示浮点数,而二进制无法准确表示某些
十进制小数。
例如,0.1 在二进制中是无限循环的。
这种情况下,
计算机需要对结果进行舍入,从而导致了舍入误差。
舍入误差可以通过使用高精度算法或减小计算的次数来减少。
在某些情况下,可以使用符号计算或数值计算软件来避免舍入误差。
总的来说,截断误差和舍入误差是数值计算中常见的误差类型,了解和掌握它们的特点和规律可以帮助我们更好地进行数值计算。
数值分析要点
− −
y2 )( y − y2 )( y0
y3 ) − y3
)
+
f
−1
(
y1
)
(
(y y1
− −
y0 )( y y0 )( y1
− −
y2 )( y − y3 ) y2 )( y1 − y3
)
+
f
−1
(
y
2
)
(
(y − y2 −
y0 )( y y0 )( y2
− −
y1 )( y − y1 )( y2
于求例解如的:问函题数,这f (种x) 简用化泰带勒入(T误ay差lo称r) 多为项方式法误差或截断误差。
Pn (x) =
f (0) +
f ′(0) x + 1!
f ′′(0) x 2 2!
+⋯+
f (n) (0) x n n!
近似代替,则数值方法的截断误差是泰勒余项。
4、舍入误差 计算机只能处理有限数位的小数运算,初始参 数或中例间如结:果用都必3 .须141进59行近四似舍代五替入运π ,算产,生这的种误误差差称为舍入误差。
⎢ ⎢
y1
⎥ ⎥
⎢y2 ⎥ ⎢⎢⋯ ⎥⎥
⎢⎣ yn ⎥⎦
其系数行列式为范德蒙 ( Vandemonde )矩阵行列式
1 x0 x02 ⋯ x0n
∏ 1
A= ⋯
x1 ⋯
x12 ⋯
⋯ ⋯
x1n ⋯
= (xi
0≤ j<i≤n
− xj) ≠ 0
1
xn
x
2 n
⋯ xnn
因此,线性方程组的解 a0,a1,⋯,an 存在且唯一。
数值分析误差及分析
数值分析误差及分析数值分析是一种通过数学方法和计算机模拟来处理和解决实际问题的方法。
然而,由于计算机的运算能力和存储能力有限,以及问题本身的复杂性,数值分析往往会引入一定的误差。
误差是指数值计算结果与真实值之间的差异,它分为截断误差和舍入误差两种类型。
截断误差是由于在数值分析过程中对无限小量和无限级数的截取而产生的误差。
无限小量是指小到可以忽略不计的量,无限级数是指由无限多个项相加的数列。
在实际计算过程中,为了获得可计算的结果,人们往往只考虑有限项的计算,这就导致了截断误差的出现。
截断误差的大小与问题本身的性质以及截止条件的选择有关。
舍入误差是由于计算机内部的浮点数表示方式而引入的误差。
计算机内部使用有限的位数来表示实数,这就不可避免地导致了浮点数的精度问题。
当计算结果需要表示的位数超过了计算机所能表示的范围时,就会发生舍入误差。
舍入误差的大小与计算机的表示精度以及计算过程中的计算次数有关。
为了减小误差,提高数值分析的精度,可以采取以下方法:1.增加计算机的位数:增加计算机的位数可以扩大浮点数的表示范围,从而减小舍入误差的发生概率。
2.使用更高精度的数据类型:在一些特殊情况下,为了提高计算结果的精度,可以使用更高精度的数据类型,如使用双精度浮点数代替单精度浮点数。
3.改进算法:优化算法可以减小截断误差的影响,例如使用数值积分的自适应算法、迭代法等。
4.选择合适的截止条件:在数值分析过程中,需要选择适当的截止条件。
截止条件的选择既不应过于严格,以免造成大的截断误差,也不应过于宽松,以免在计算机内部引入较大的舍入误差。
5.进行误差分析:在数值分析过程中,应该对误差进行分析和估计。
可以通过理论方法、数值试验和统计方法等途径来估计误差的上界或下界,从而评估计算结果的可靠性。
总而言之,数值分析误差是不可避免的,但可以通过增加计算机位数、改进算法、选择合适的截止条件、使用高精度数据类型和进行误差分析等方法来减小误差,提高数值分析的精度和可靠性。
《数值分析误差》PPT课件
绝对误差往往是未知的,而只知道它的一个上
限,此上限|e*|=|x-x*|记为ℇ*, 称为绝对误差
限(accuracy)。
工程上常记为x=x*±ℇ* ,例如
1ex2dx 0.740 3.006 0
相对误差 (relative error)
e r * e x * * x x * x * 或 e r * e x * x x x *
1.3 方法误差 (截断误差 Truncation Error〕
在数学模型〔包括参数值〕确定以后,就常要考虑 选用某种数值方法具体进展计算,许多数值方法都 是近似方法,故求出的结果与准确值之间是有误差 的,该误差称为截断误差或方法误差。例如,函数 f(x)用Taylor多项式
f ( x ) p n ( x ) f ( 0 ) f 1 '( ! 0 ) x f ' 2 '( ! 0 ) x 2 f ( n n ) ! ( 0 ) x n
1. 来源与分类 ( Source & Classification )
• 模型误差 • 参数误差(观测误差) • 方法误差(截断误差) • 舍入误差
1.1 模型误差 (Modeling Error)
用计算机解决实际问题时,首先要建立数学 模型,各种实际问题是十分复杂的,而数学 模型是对被描述的实际问题进展抽象、简化 而得到的,往往忽略了一些次要因素,因而 是近似的,我们把数学模型与实际问题之间 出现的这种误差称为模型误差。如自由落体 公式
|e r *| (x y x )* (x y * * y * ) 0 .0 0 .0 0 0 .0 1 3 1 7% 0
可以看到相对误差比较大.
数值分析中的误差传播理论
数值分析中的误差传播理论数值分析是一门应用数学领域的重要学科,旨在通过利用数值方法解决实际问题。
然而,在数值计算中,由于各种因素的影响,我们无法完全避免误差的产生。
误差传播理论是数值分析中一种重要的理论工具,旨在帮助我们理解和控制误差的产生和传播过程。
一、误差及其分类在数值计算中,由于测量和计算过程中的不确定性,我们得到的结果往往与真实值存在一定的偏差,我们称之为误差。
根据误差产生的原因和性质,误差可以分为以下几类:1. 舍入误差:由于计算机存储空间的有限性,无法精确表示某些实数,从而导致舍入误差的产生。
舍入误差是最常见的一类误差,通常通过保留足够的有效数字或采用更高的精度计算来减小。
2. 截断误差:当我们使用近似方法计算某个函数或数值时,由于截断计算过程中的无穷级数或无穷小量,截断误差会被引入。
减小截断误差的常用方法是增加计算的步骤或者使用更高阶的近似方法。
3. 模型误差:数值计算往往需要建立数学模型来描述实际问题,而模型本身的不准确性会引入模型误差。
模型误差可以通过改进数学模型或采用更适合实际情况的模型来减小。
二、误差传播的基本原理误差传播理论是基于线性近似的思想,它假设误差在传播过程中是线性累积的。
根据误差传播的基本原理,我们可以通过对误差的传播规律进行研究,从而评估计算结果的可靠性。
误差传播理论的基本公式为:δf = |∂f/∂x₁|δx₁ + |∂f/∂x₂|δx₂ + ... + |∂f/∂xₙ|δxₙ其中,δf表示函数f的误差,δx₁、δx₂、...、δxₙ表示自变量x₁、x₂、...、xₙ的误差,|∂f/∂x₁|、|∂f/∂x₂|、...、|∂f/∂xₙ|表示函数f对自变量的偏导数。
该公式表明,函数f的误差δf由自变量的误差δx₁、δx₂、...、δxₙ以及函数f对各自变量的偏导数共同决定。
三、利用误差传播理论进行数值计算在实际的数值计算中,我们可以利用误差传播理论来评估计算结果的误差范围,并采取相应的措施来减小误差。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第9章 数值分析中的误差练习9.1 (B) 1.B 2. A 3.量纲 4. 半个 5.D 6. 3位 7.C 8. B 9. ≤ 10. 0.5mm 练习9.2(B) 1. )()(21x x εε+ )()(21x x εε+2. )()(1221x x x x εε+,221221)()(x x x x x εε-3. 舍入误差不增加4. )()(x e x f '5. 使用数值稳定的算法;防止两个相近数相减;简化计算步骤,减少运算次数;避免除数的绝对值远小于被除数的绝对值;防止大数“吃掉”小数 习题9 1. 0.00005 0.017% 四位有效数字 0.005 0.017% 四位有效数字 0.0005 0.0017% 五位有效数字 0.5 0.0017% 五位有效数字 0.00005×105 0.017% 四位有效数字 2. 5.5(Ω) 0.2375(Ω) 4.32% 3. (1) B (2)A (3)A4.)()()(2121+=±x x x x εεε )()()(2212121121±±±=±x x x x x x x x x x r r r εεε)()()(122121+≈x x x x x x εεε )()()(1221+≈x x x x r r r εεε22122121+=x x x x x x x )()()(εεε )()()(2121+=x x x x r r r εεε 5. 0.0056. 0.00333…7. 取四位. 利用定理2.第10章 线性方程组的数值解法练习10.1(A)1. (2,1,-1)T ,2. (1,2,3)T ,3. (1,-1,2)T, 4. T)3678.0,05113.0,4900.0(--≈*X(B)1. C2. 系数矩阵的各阶主子式均不为0.3. B4. 矩阵A 是严格对角占优矩阵5. 见教材第10章公式(1.6)6. 见教材第10章公式(1.8) 练习10.2 (A)1.(1,2,-1,3)T,2. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=200140111112010001U L3. Y =(14,-10,-72)T,X =(1,2,3)T. 4. Y =(2.4493,11.247,85.254)T,X =(1,0,23)T.5. X =(5,4,3,2)T .6. X =(2,2,1)T.7. Y =(3,0.5,-1)T ,X =(2,1,-1)T. (B)1. 初等矩阵 上(下)2.单位下三角形 上三角形3. A *4.D5. 三对角线矩阵6. ⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=10001000001000010000132nl l lL ,⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=--n n n u r u u r u r u 000000000000001132211U 7.C 8. 对称正定矩阵练习10.3 (A)1. T T )1,2,3(,)9998813.0,999838.1,00032.3()10(==*X X2. 精确解(1,2,3)雅可比迭代法:X (9)=(0.9998, 1.9998, 2.9998)T高斯−−赛德尔迭代法:X (5)=(0.9997, 1.9999, 2.9999)T3. 精确解(-3,3,1) (B)1.X =B 0X +f X (k +1)=B 0X (k )+f k =0,1,2,… 2.≠0 3.B 4.B D 5.A 6.D 习题101. 顺序消去法X =(1.000015,1.00001,1.000006,0.9999858)T列主元消去法X =(1.000000,1.00000,1.00000,1,0000000)T2. X =(1.9272,-0.69841,0.90038)T3. Y =(3.77,-3.06,4. 66,4.22)TX =(-0.329,0.322,2.37,1.04)T4. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=113012001L ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=700720412U5.⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡--=137910161015161013100010L ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡7419100109-10370031323100-26=011U Y =(6,-3,23/5,-191/74)T X =(1,-1,1,-1)T6. Y =(5,0,3)T X =(2,2,1)T7.略 8. 雅可比迭代法T )9999990.2,0000040.1,9999950.1,9999941.0()24(--=X高斯-赛德尔迭代法T )9999990.2,0000040.1,9999970.1,9999966.0()14(--=X9. 雅可比迭代法收敛,高斯-赛德尔迭代法发散. 10.略 11.∣a ∣>2第11章 函数插值与最小二乘拟合练习11.1(A)1. 6671≈35=33331≈34=22+31=111.)(.)()()(P P x x P 2. 6671=2≈27+-61=222.)()()(P x x x P3. P 2(x)=232371x x -+ 4. 4.7943(0.6-x)+5.6464(x -0.5) sin0.57891≈0.546675. 1+2x(B)1. 节点; 插值多项式; 被插值函数2. C.3. B4.))(())(())(())(())(())((120210210120201021------------x x x x x x x x x x x x x x x x x x x x x x x x 5. B练习11.2 (A)1. f(x 0,x 1)=-5, f(x 1,x 2)=-1,f(x 2,x 3)=9;f(x 0,x 1,x 2)=2,f(x 1,x 2,x 3)=5,f(x 0,x 1,x 2,x 3)=12. f(x)=x 3+x 2+x+1 3. 39.0625(用牛顿插值多项式,))(()()(6-5-+5-11+25=3x x x x P )4. ∆y 0=0.02119,∆y 1=0.02020,∆y 2=0.01931,∆y 3=0.01848,∆y 4=0.01773,∆2y 0=-0.00099,∆2y 1=-0.00089,∆2y 2=-0.00083,∆2y 3=-0.00075,∆3y 0=0.00010,∆3y 1=0.00006,∆3y 2=0.00008, ∆4y 0=-0.00004,∆4y 1=0.00002,∆5y 0=0.00006 5. N 4(x)=(x 3-4x 2+3)-)4)(3)(2)(1(247----x x x x (B)1.),()()(344343--x x f x x x f x f 2. B 3.A 4.D5.))(())(())((120222101120100--+--+--x x x x y x x x x y x x x x y 6. C 练习11.3 (A)1. 2.666 67 2. 略3. ⎪⎩⎪⎨⎧3≤≤239+52-666723+333332≤≤166673-12+33338-21≤≤0+66672+66671-=232323x x x x x x x x x x x x x S ......)(4.⎪⎩⎪⎨⎧5≤≤433-6525+6255+37504≤≤27+254-8751+250-2≤≤11-751+3750+1250-232323x x x x x x x x x x x x .........(B)1.C2.见教材第11章公式(3.1) 3.A4. S(x)在[a,b]上具有2阶连续导数 S(x j )=y j (j=0,1,2,…,n) 在每个子区间[x k ,x k+1](k=0,1,2,…,n-1)上,S(x)是3次的多项式 .5.B 练习11.4 1. y=-1.43+6.43p 2. S=5.34+0.30t 3. y=5.045-4.043x +1.009x 2 4.y=11.6789e (-1.1109/x)(B)1. B2.⎪⎪⎩⎪⎪⎨⎧=+=+∑∑∑∑∑1=1=21=1=1=nk k k n k k n k k nk k n k k y x x b x a y x b na 3. (lnx k ,y k )(k=1,2,…,n) 习题111. P 3(x)=)924152(3123-+--x x x 2. P 3(x)=0.2(x 3-13x 2+69x-92) 3. P 1(x)= 2x-1 4. N 3(x)= x 3+x 2+x+17. 0=2221=2228171),...,,(,),...,,(f f8. 0.603 144 9.0.64310.S (x )=⎪⎪⎩⎪⎪⎨⎧530≤≤450450-10889+530-71257-450-37021-530-66581450≤≤390390-188311+450-418310-390-22112-450-86972390≤≤300300-95446+390-10786-300-91321-390-70942300≤≤250250-96610+300-0110-250-8774-300-7617633333333..).(.).(.).(.).(...).(.).(.).(.).(...).(.).(.).(.).(...).(.).(.).(.).(.x x x x x x x x x x x x x x x x x x x x 11.y=0.0213e 0.272x实习作业:1. P(0.5635)=826116 2. N n (0.596)=0.631918; N n (0.895)=1.01937第12章 数值积分与微分练习12.1(A)1. (1)二次代数精度,1563,561±==βα ; (2)二次代数精度,a 0=1/3, a 1=4/3, a 2=1/32. 五次代数精度(B) 1. A 2.B 练习12.2 (A)1.略2.略3.1.71828 用梯形公式分476个子区间,用抛物线公式分6个子区间(B )1.C k n kn()=∑=01 (归一公式),)(n k C 与a ,b 无关,且)()(n k n n k C C -=(称为对称性)2. B 练习12.3 (A)1. 1.15472. 2.0013889 (B)1. 2n+12. B3. 1==3110A A ,练习12.4 (A)1. 取平均为)..(.47216+4861321=97914 2. -0.4 -0.2(B)1. B2.)()()()()()(210220121003+4-21≈'+-21≈'-4+3-21≈'y y y hx f y y hx f y y y hx f 习题121. 0.311 680 24; 0.310 248 53; 0.310 268 842. 1.148 714 467; 1.147 792 8573.4. 略5. 用复化梯形公式,4等分区间,近似值为T 4=5.058 337 用复化抛物线公式,2等区间,近似值为S 2=5.033 002 用科茨公式 近似值为C =5.03292 6.3 7. A 0=2 A 1=-3 A 2=3第13章 方程求根练习13.1(A )1. 0.9212.3.1455(B ) 1.B 2.C 3.ln()ln ln b a --ε24.[1.5,2]练习13.2 (A )1.(1) 0.7391 (2) 0.20391 (3)2.09455 2.0.090525(B )1.A2.x x =ϕ()的形式3.B4.Cϕϕ()(~)(~)~x x x x x x x xn n n n n n n+++++---+1112112 练习13.3 (A )1.x 02= x 1=1.5829 x 2=1.5009 x 3=x 4=1.49732. x 1=1.5970149 x 2=1.5945637 x 3=x 4=1.59456213. 3.317(B)1.线性化2.x f x f x n n n n ----'=11112()()(,,) ,0≠'1-)(n x f 3.C 4.B5.,...),,()(210=5+231=21-1-n x x x n n n6.D练习13.4(A )1. 2.09455 2. 3.14619 (B ) 1. ,...),()()()()(32=---=2-1-2-1-1-1-n x x x f x f x f x x n n n n n n n 2.A 3.D 4.C练习13.5(A)1. x ***..=⎡⎣⎢⎤⎦⎥=-⎡⎣⎢⎤⎦⎥x x 120254098360497267762.x x 1121125225()()..⎡⎣⎢⎤⎦⎥=⎡⎣⎢⎤⎦⎥,x x 1222120277778()().⎡⎣⎢⎤⎦⎥=⎡⎣⎢⎤⎦⎥(B)1. 线性化2.B3.C习题131.(1) -0.56714 (2)1.267172.(1)1.9405 (2) 0.4502 (3)-1.154173.(1) 3.63197 (2) 0.011754.(1)0.33767 1.30749 (2)1.386165. (1) (0.8260,0.5636)T(2) (0.9991,-0.0888)T第14章 常微分方程的数值解法 练习14.1(A) 1. 略2.将解列入表14-1中。