数值分析第五版习题答案清华大学出版社
数值分析第五版课后答案2篇

数值分析第五版课后答案2篇数值分析第五版课后答案(一)第一章1.1 机器精度的数值为2^-52 ≈2.22 × 10^-16。
1.2 Example 1.2设f(x) = (1 - cosx)/sinx,则f(0)的分母为0,无法进行数值计算。
1.3 Example 1.3设f(x) = (1 - cosx)/sinx,则f(0)的分子为0,因此有f(0) = 0。
1.4 Example 1.4(a) 将x的值从1.8改为1.799,则f(x)的值由-0.000000000000159为0.000000000000313,差值为0.000000000000472。
(b) 我们有f'(x) = sinx/(1 - cosx) - 1/sin^2x。
将x的值从1.8改为1.799,利用f(x)和f'(x)的值可以得到下面的近似式:f(x + Δx) ≈ f(x) + f'(x)Δx = -0.000000000000159 + 0.449787416887455×0.001 = -0.000000000000137。
与(a)中的结果相近。
1.5 Example 1.5(a) 当x很接近于0时,函数值的符号取决于cosx的符号,其中cosx接近于1。
因此,函数值为正。
(b) 当x很接近于π时,函数值的大小趋于无穷大,因为分母趋向于0,而分子不为0。
1.6 Example 1.6(a) 因为函数在x = 0处是奇函数,所以它的导数为偶函数。
(b) 首先,我们有f''(0) = -2,因此x = 0是最大值。
其次,我们有f''(x) = -2 - 8sin^2x。
由于-f''(x)在x = 0处是正的,我们有当x越接近0时,f''(x)越小,也就意味着函数在x = 0处是严格的最大值。
1.7 Example 1.7(a) 我们有f(x) = x^3 - 2x^2 - 5x + 6,f'(x) =3x^2 - 4x - 5和f''(x) = 6x - 4。
数值分析第五版习题答案全部清华大学出版社

6
12.计算 f ( 2 1) ,取 2 ,利用下列等式计算,哪一个得到的结果最好?
1 , ( 2 1) 6
(3 2 2) 3 ,
6
1 , 99 70 2 。 (3 2 2) 3
解:设 y ( x 1) , 若x 若通过
若通过 (3 2 2) 计算 y 值,则
3
y * (3 2x* )2 x* 6 y * x* * 3 2x y * x*
若通过
1 计算 y 值,则 (3 2 2) 3 1 x* * 4 (3 2 x )
1 1 1 10 4 10 3 10 3 2 2 2 3 1.05 10
* * * (2) ( x1 x2 x3 ) * * * * * * * * * x1 x2 ( x3 ) x2 x3 ( x1 ) x1 x3 ( x2 )
1 1 1 1.1021 0.031 10 1 0.031 385.6 10 4 1.1021 385.6 10 3 2 2 2 0.215
* * * *
*
*
*
* * *
*
*
1 ( x1* ) 10 4 2 1 * ( x2 ) 10 3 2 1 * ( x3 ) 10 1 2 1 * ( x4 ) 10 3 2 1 * ( x5 ) 10 1 2
* * * (1) ( x1 x2 x4 ) * * * ( x1 ) ( x2 ) ( x4 )
解: y0
2 1.41
1 ( y0 *) 10 2 2
又 yn 10 yn 1 1
数值分析第5版课后习题答案(清华大学出版社)-第一章

′ e * ( f1 ) = f1 e * (1.4) =
对于 f 2 = (3 − 2 2 ) 3 ,
1 1 ′ e * ( f 2 ) = f 2 e * (1.4) = 6(3 − 2 × 1.4) 2 × × 10 −1 = 0.12 × 10 −1 < × 10 −1 ,没有有效数 2 2 字;
*
1 1 1 = (0.031 × 385.6) × 10 − 4 + (1.1021 × 385.6) × 10 −3 + (1.1021 × 0.031) × 10 −3 ; 2 2 2 −3 −3 −3 = 0.59768 × 10 + 212.48488 × 10 + 0.01708255 × 10 = 213.09964255 × 10 −3 = 0.21309964255
′ PAP −1 Px Ax PAx ′ 6、证明: A max = = max = max = ′ Px Px x PAP −1 。
7、证明:由范数的等价性,存在常数 C1 和 C 2 ,使得 C1 x s ≤ x t ≤ C 2 x s ,则有
C1 Ax s ≤ Ax t ≤ C 2 Ax s ,并且
∂f e (x + x + x ) = ∑ k =1 ∂x k
* * 1 * 2 * 4
n
* * * * ε ( x k ) = ε ( x1 ) + ε ( x 2 ) + ε ( x 4 ) ;
*
=
1 1 1 × 10 − 4 + × 10 −3 + × 10 −3 = 1.05 × 10 −3 2 2 2
y 0 = 1.41 1 可知, ε * ( y 0 ) = × 10 − 2 , y n − y n = 10( y n −1 − y n −1 ) ,即 2 y n = 10 y n −1 − 1
数值分析第五版1-3章

p( x) a0 a1 x a2 x 2 an x n
2015-6-4
16
第2章 插值法
引言
拉格朗 日插值
牛顿插值
埃尔米 特插值
分段低 次插值
三次样 条插值
2015-6-4
x x
f ( x) f ( x* ) f ( x)
x x
xf ( x) f ( x)
C p 10 即认为是病态
f ( x) x n
9 第1章 数值分析与科学计算引论
研究对象 作用特点
数值计算 误差
误差分析 避免危害
数值计算 算法设计
数学软件
2. 算法的数值稳定性 定义3 一个算法如果输入数据有误差,而在计算过程中舍入 误差不增长,则称此算法是数值稳定的,否则称此算法为不 稳定 的。 例1.1:P.9 I n e
2015-6-4 7 第1章 数值分析与科学计算引论
1 10 ( n 1) 2(a1 1)
研究对象 作用特点
数值计算 误差
误差分析 避免危害
数值计算 算法设计
数学软件
3 数值运算的误差估计
* * * * 1. x1 与x2 为两近似数, 误差限为 ( x1 ), ( x2 ), 则 : * * * * ( x1 x2 ) ( x1 ) ( x2 ); * * * * * * ( x1 x2 ) x2 ( x1 ) x1 ( x2 ); * * ( x1 / x2 ) * * * * x2 ( x1 ) x1 ( x2 ) * x2 2
数值分析第五版课后答案

数值分析第五版课后答案(ii )2/(x ) = Imr0.40.50.60.7 0.8 lar一 0.916 291 一 0.693 147 一 0・ 510826-0. 356 675-0.223 144用线性插值及二次插值计算InO. 54的近似值•解 依据插值误差估计式选距离0. 54较近的点为插值节点,并建立差商 表如下:一 0.693 147-0.510 826 - 0.916 291写出Newton 插值多项式M(H ) =- 0.693 147 + 1.823 210Q — 0.5)N2)= M (_r) + (—0.204 115〉(工一0. 5)匕一0・6)计算近似值Ni (0. 54) =一 0.693 147+ 1.823 210(0. 54 — 0. 5) =—0.620 218 6弘(0.54) = N 】(0.54) — 0.204 115(0. 54 - 0.5X0. 54-0.6) =-0.616 8394・设门为互异节点(j = 0.1 ■…山).求证:A(I )三卫(上=0, 1 ■…,Q;n(ii )心一工)铅(门三o 仏=1. 2. •••■" 证明 (i )令fS 』工X 若插值节点为X/7 - 0,1 则/<x )的n次播值多项武为["(工)=工球丿3插值余项为R”(王〉=/(X )— L n (X )= /—(/)(n + 1)!/X—Ti-CkXVZ又因为k < 所以严)(0 = 0,R 心)二 0x 0 = 0. 5 X] = 0. 6工2 = 0. 4二> -0.204 1151.823 2102. 027 325所以丿・0 1 -n L'? /xsr ("卜;(_"“(/〉 r —0 丿•()L ' r / SCOg ( . ) (一 x)k -'x' =(彳一 Qi 三 05.设 /(x) 6 C 2[a, 6]且 /(a) = fib) = 0.求证: max | f(x) £(b —a),max | /z (j) \ a^r^ib .O心疋6证明 令x = a 和工=人以此为插值节点•则插值多项式为Li (工)= /(a) -—; + f(b) Y —- 三 0<2—o b — a应用插值余项公式有y*7(^) (X — a)(.x — 6) Wmax | /(g) I max | (x — a)(x — b) | / Wb a<jfCZ> _(6 — a )2 max | fXx ) | O aM 临 b6.在一 4<x<4上给出r (T )= e 『的等距节点函数表,若用分段二次插值求e 「的近似值,要使截断误差不超过10一&,问使用函数表的步长h 应取多少?解 若插值节点为IT , r 和工沖则分段二次插值多项式的插值余项为式中Ml = Xi — h,工沖=$ +札\R :(r) l^ye 1max | (文—刀_)) (_r —兀)〈工—J7°j )丨 0插值点个数< W 6 得 A < 0.006 5&是奇数,故实际可采用的函数值表步长7•若必=2S 求及解 根据向前差分算子和中心差分算子的定义进行求解£(:)(-】〉1巧” =£(:)(-1)-皿=孑y” = (F? — F~T )°y” = (E"r )*(E — IYy n =「2$% = L (%) = g = 2—8.如果fl 工、是刃次多项式,记= f (j--T-h ) —/(T ).证明/(x )的 &阶差分Nfa )(0W 是rn-k 次多项式,并且A^7(T> = 0(/为正 蔓牧).证明 对加次多项式/(才)应用Taylor 公式有A/(x ) = /(z + A ) —/(j ) = /( J )A H- rr/^x ) + ••• 4- Jf"' (x )Z! 初!即△/(/)为m- 1次的多项式・= △(△/&)),对加一 1 > 0次多项式应用上述推理过程知 △(△/(工))=庄只工)是加一2次的多项式.依此过程递推,知A7<^X0<Xr<r«)为m-k 次多项武. 所以必工)为常数,故 s = 0(/为正整数).9. 证明 A (/*g* ) = /* Ag* 4-A/*.证明 A/igJ = /n-ign-i ~ Ag* = /n-igHi - fkgkn 十/*gi - fkgk = gtrl (人+1 — 人> + fk(g^l 一创)=g 屮+ 介厶®15.证明两点三次Hermite 摘值余项是尺3(刃='‘4 ;目(工—九)'(H —)?, E €(N ,才屮)并由此求出分段三次Henniw 猶值的课差限・证明 若工W [工―文屮]・且插值多项式满足条件円3 {竝)=/(竝几 H3(X H -1)=产(工屮)H ; ( Z* ) = f (x> ) * H3' (J T H -I )=(.r*41)1 4-(- 4) 十 0. 006 581 268 冬 1 217 旦 N4 —(—4) F T8T2162 0. 006 579知插值余项RQ) = /(文)一耳(工> 有二重零点g和文卄故设R(攵)=以文)0 —比)?(文一攵申〃确定函数恥才几当JC = X*或工屮时來工)取任何有限值均可I当才H忑,J•屮时“&(仏°文屮),构造关于变量t的函数g(r) == /([)—丹3(『)一总(才)(〔一=*)2((— X*+l )2 显然有g(文▲)= 0. g(i?) = 0. ) = og'(r*.〉= 0, £心屮)=0在S ,工]和Dr, z*+l J上对g(T)使用Rolle定理,存在® €(无,才〉及少W (w, x*-ti)使得&'(》)=0, g'(%)= 0在a ,巾),Cyl *罪),<72« x*+i)上对g'(=)使用Rolle定理,存在供| € 5,巾),巾?€(6,%)和阻屮6(%,XHI)使得g"(知)=g"(?!2)= g"(少.屮)=0再依次对g(0和g"(“使用Rolle定理,知至少存在(比,工屮)使得gW(E)=而g⑷⑺=一虹小4!,将"弋入•得到£€5 •工屮)推导过程表明W依赖于工点,及=•综合以上过程可知R(T)= “(&(a■一忑)2(工一卫^)2下面建立分段三次Hermite插值的谋差限.记h (小为/Cr)在[a,刃上的基于等距节点的分段三次Hermite插值函数.x k = a+kh 4 = 0, !••• ♦ n), h = b — a■n在小区间[去,/小]上有I /(x)— /A<x)| —右 | 严(£)|(X— X*)-(X— XH-1)2 <7f max \尸4)(力))max (_r —业)?(工一z屮尸而最值0 才=十妙 ]max (工一及)■(工一 z>+! 「L l 「• , maxs"(5 ― l )2h 4 = r k n<<<! 16进而得误差估计1 /(文)-越空简|八(如】6・求一个次数不离于4次的多项式PCr 〉•使它满足P(0) = P(0) = 0, P(l) =P71) = HP(2) = 1.解法一 利用Hermite 插值可得到次数不高于4的多项式几== 1;为==打 W f > = 0 •加I = 1H 3(X )=(才)+ /(文)◎(才)=(1 一 2 三「卫■)(才二空)2 =(1 + 2刃1)2氐—XI 竝一 4G&) = (1-2 J ~-r| )( - )2 = (3 — 2&)疋Jj —竝 XI — To仇(工)=兀(工一 1)?向=(工一 1)JT 2所以Hj (2) = (3 — 2x )x 2 + (1* — 1 )J -? =— T 3 + 2z~设 = H 3(X )4-A (T -^)2(J —T ))2,其中・A 为待定常数,令 F (2)=1得于是P3十一尸这样可写岀Newton 插值公式P (x ) = 0 + 0(乂一0)十 1(工一0)? — 1(広一0)?(工一 1) +— 0)'($ — l )? =— 1) + 4-工?(&一 1)?=解法二(带重节点的Newton 插值法)建立如下差商表:-124 4J-x 2 (r ~ 3): 417 •设f (.C 二厂丄g 在一 5€工€5上取"=10•按等矩节点求分施线1 f JT性插值函数ha )・计算各节点间中点处Z A (J -)与/(x >的值,并佑计课差.解 若 = 5,r lc = 5,则步长 A = ------------- ---- -- -- = I =— 5+ ih = — 5 +n2(ow?w 10).在区间Cx-上•分段线性插值瓯数为/1°(X )= /(X,)工汁】一広+工一 rTT7T F+不分段线性插值函数定义如下:各节点间中点处函数值及插值函数值如下所示:估计谋差:在区间[乙,刀+门上lf (jr )—击厂(。
李庆扬数值分析第五版习题答案清华大学出版社

李庆扬数值分析第五版习题答案清华大学出版社Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()n f x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅ 且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字;*20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少解:球体体积为343V R π= 则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=故度量半径R 时允许的相对误差限为1(*)10.333r R ε=⨯≈ 6.设028Y =,按递推公式1n n Y Y -=(n=1,2,…) 计算到100Y27.982≈(5位有效数字),试问计算100Y 将有多大误差解:1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =,若取27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
数值分析第五版习题答案清华大学出版社

1.设 x 0 , x 的相对误差为 ,求 ln x 的误差。
解:近似值
*
x 的相对误差为
* e* x * x
= er
x*
x*
而 ln x 的误差为 e ln x *
ln x * ln x
1 e*
x*
进而有 (ln x*)
2.设 x 的相对误差为 2%,求 xn 的相对误差。
解:设 f ( x )
cos x 的近似
值时,采用的线性插值法插值余项不为 0,也会有一定的误差。因此,总误差界的计算应综
合以上两方面的因素。
当 0 x 90 时,
令 f ( x ) cos x
1 取 x0 0, h ( )
60
1 60 180
令 xi x 0 ih , i 0,1,..., 5400
10800
则 x5400
3
2.给出 f ( x ) ln x 的数值表
X
0.4
lnx
-0.916291
用线性插值及二次插值计算
解:由表格知,
0.5
0.6
-0.693147
-0.510826
ln 0.54 的近似值。
x0 0.4, x1 0.5, x 2 0.6, x3 0.7, x4 0.8; f ( x 0 ) 0.916291, f ( x1 ) 0.693147 f ( x 2 ) 0.510826, f ( x3 ) 0.356675 f ( x 4 ) 0.223144
2 1.41
( y 0 *)
1
2
10
2
又 y n 10 yn 1 1
y1 10 y0 1
( y1*) 10 ( y 0 *) 又 y 2 10 y1 1
清华大学第五版《数值分析》课后答案

第1章绪论内容提要#〜误差度量1数值分析研究两类误差:舍入误差和截断误差,由于计算机字长的有限性,对相关数据进行存储表示时便产生舍入误差,计算机必须在有限的时间内得到运行结果,于是无穷的运算过程必须截断为有限过程,由此产生截断误差,2,误差的度量分式有:绝对误差(限)、相对误差(限〗和有效数字,设?是真值工的一个近似,绝对误差为一:!相对误差为& ,绝对误差限〉和相对误X X差限6^ 〉分别是〉 |和^(:^ ^|的上限,3^对于非零近似值^的如下规格化标准形式X^ ^ 10^ X0#!1X2'&X&,&!' ?X I ^0 〈1. 1〉如果存在尽可能大的&,使得〉| & ^乂10"-",则称?有"位有效数字.进而当&^》时,称X,是有效数.4,有效数字和相对误差的关系定理1. 1 如果形如式〈1. V的有&位有效数字,则定理1.2如果形如式〈1. 0的:^的相对误差满足^|《"二" X化1-"则纟^至少有&位有效数字,二、浮点数系统对于5+ ^ + 2位的浮点数系0表示二进制阶码数值的二进制位数〃表示尾数的二进制位数,其他两位表示阶码和尾数的符号〉,机器数绝对值的范围是2-21〜22'-、实数表示的相对舍入误差限是2-'.当数据的绝对值大于22'-1时,计算机非正常停机,称之为上溢,当非零数据的绝对值小于2-2',用机器零表示,精度损失,称之为下溢,、误差传播如果在运算过程中舍入误差能够得到控制,或者舍入误差的增长不影响产生可靠的结果则纟称该算法是数值稳定的,函数值绝对误差传播公式如下^/(^" 丫) ## /(;:)〉 1 2〉^(/(^" ^-^:》#亡"";二…、^ 〉(丄门)!.^^")〉#| /'(?) |〈1.4〉、数值稳定性不同的教材对数值方法稳定性的定义有所不同,有的要求随计算过程的深入误差不增长,有的则要求误差增长速度不能太快^只要不影响产生具有有效数字的近似值即认为是稳定的,读者应注意教材中的定义.随着学习的深入,会针对各种具体算法给出稳定性的确切定义,^ 2 ^典型例题与解题技巧【例1】求!&的近似值,使其绝对误差限精确到1乂1。
数值分析作业答案——第五版

数值分析第一次作业及参考答案1. 已测得函数()y f x =的三对数据:(0,1),(-1,5),(2,-1),(1)用Lagrange 插值求二次插值多项式。
(2)构造差商表。
(3)用Newton 插值求二次插值多项式。
解:(1)Lagrange 插值基函数为0(1)(2)1()(1)(2)(01)(02)2x x l x x x +-==-+-+-同理 1211()(2),()(1)36l x x x l x x x =-=+故2202151()()(1)(2)(2)(1)23631i i i p x y l x x x x x x x x x =-==-+-+-++=-+∑(2)令0120,1,2x x x ==-=,则一阶差商、二阶差商为0112155(1)[,]4,[,]20(1)12f x x f x x ---==-==-----0124(2)[,,]102f x x x ---==-22()1(4)(0)1*(0)(1)31P x x x x x x =+--+-+=-+2. 在44x -≤≤上给出()x f x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?解:()40000(),(),[4,4],,,, 1.xk x f x e fx e e x x h x x h x x th t ==≤∈--+=+≤考察点及(3)200044343()()[(()]()[()]3!(1)(1)(1)(1)3!3!2.(4,4).6fR x x x h x x x x ht t tet h th t h e heξξ=----+-+≤+⋅⋅-=≤∈-则436((1)(1)100.006.t t th h--+±<<在点取到极大值令 得3.求2()f x x=在[a,b]上的分段线性插值函数()hI x,并估计误差。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
若,则
则
则法方程组为
从而解得
故
均方误差为
21。在某佛堂反应中,由实验得分解物浓度与时间关系如下:
时间
0 5 10 15 20 25 30 35 40 45 50 55
浓度
0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.62 4.64
用最小二乘法求。
计算到。若取(5位有效数字),试问计算将有多大误差?
解:
……
依次代入后,有
即,
若取,
的误差限为。
7.求方程的两个根,使它至少具有4位有效数字()。
解:,
故方程的根应为
故
具有5位有效数字
具有5位有效数字
8.当N充分大时,怎样求?
解
设。
则
9.正方形的边长大约为了,应怎样测量才能使其面积误差不超过?
解:正方形的面积函数为
证明:
若
令,可得
当时,
当时,
又,故
得证。
10。证明切xx多项式满足微分方程
证明:
切xx多项式为
从而有
得证。
11。假设在上连续,求的零次最佳一致逼近多项式?
解:
在闭区间上连续
存在,使
取
则和是上的2个轮流为“正”、“负”的偏差点。
由切xx定理知
P为的零次最佳一致逼近多项式。
12。选取常数,使达到极小,又问这个解是否唯一?
解:
令
则在上为奇函数
又的最高次项系数为1,且为3次多项式。
与0的偏差最小。
从而有
13。求在上的最佳一次逼近多项式,并估计误差。
解:
于是得的最佳一次逼近多项式为
即
误差限为
14。求在上的最佳一次逼近多项式。
解:
于是得的最佳一次逼近多项式为
15。求在区间上的三次最佳一致逼近多项式。
解:
令,则
且
令,则
若为区间上的最佳三次逼近多项式应满足
第一章 绪论
1.设,的相对误差为,求的误差。
解:近似值的相对误差为
而的误差为
进而有
2.设的相对误差为2%,求的相对误差。
解:设,则函数的条件数为
又,
又
且为2
3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:,, , ,
解:是五位有效数字;
是二位有效数字;
是四位有效数字;
是五位有效数字;
是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ,(2) ,(3) .
其中均为第3题所给的数。
解:
5计算球体积要使相对误差限为1,问度量半径R时允许的相对误差限是多少?
解:球体体积为
则何种函数的条件数为
又
故度量半径R时允许的相对误差限为
6.设,按递推公式 (n=1,2,…)
解:
若
则步长
在小区间上,分段线性插值函数为
各节点间中点处的与的值为
当时,
当时,
当时,
当时,
当时,
误差
又
令
得的驻点为和
18.求在上分段线性插值函数,并估计误差。
解:
在区间上,
函数在小区间上分段线性插值函数为
误差为
19.求在上分段埃尔米特插值,并估计误差。
解:
在区间上,
令
函数在区间上的分段埃尔米特插值函数为
.
当时,若,
则
故测量中边长误差限不超过时,才能使其面积误差不超过
10.设,假定g是准确的,而对t的测量有秒的误差,证明当t增加时S的绝对误差增加,而相对误差却减少。
解:
当增加时,的绝对误差增加
当增加时,保持不变,则的相对误差减少。
11.序列满足递推关系 (n=1,2,…),
若(三位有效数字),计算到时误差有多大?这个计算过程稳定吗?
依此类推,在内至少有一个零点。
记为使
又
其中依赖于
分段三次埃尔米特插值时,若节点为,设步长为,即
在小区间上
16.求一个次数不高于4次的多项式P(x),使它满足
解:利用xx特插值可得到次数不高于4的多项式
设
其中,A为待定常数
从而
17.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与值,并估计误差。
解
,
设
则
故
若改用等价公式
则
此时,
第二章 插值法
1.当时,,求的二次插值多项式。
解:
则二次xx日插值多项式为
2.给出的数值表
X
0.4
0.5
0.6
0.7
0.8
lnx
-0.916291
-0.693147
-0.510826
-0.356675
-0.223144
用线性插值及二次插值计算的近似值。
解:由表格知,
若采用线性插值法计算即,
解:
观察所给数据的特点,采用方程
两边同时取对数,则
取
则
则法方程组为
从而解得
因此
22。给出一张记录用FFT算法求的离散谱。
解:
则
0 1 2 3 4 5 6 7
4 3 2 1 0 1 2 3
4 4 4 4 0 4
8 4 0 4 8 0
16 0 0 0
23,用辗转相除法将化为xx分式。
解
24。求在处的阶xx逼近。
又且
又
即计算值比准确值大。
其几何意义为,为下凸函数,梯形面积大于曲边梯形面积。
8。用xx求积方法计算下列积分,使误差不超过.
解:
0
0.7717433
1
0.7280699
0.7135121
2
0.7169828
0.7132870
0.7132720
3
0.7142002
0.7132726
0.7132717
0.7132717
当
时,多项式与零偏差最小,故
进而,的三次最佳一致逼近多项式为,则的三次最佳一致逼近多项式为
16。,在上求关于的最佳平方逼近多项式。
解:
若
且,则
则法方程组为
解得
故关于的最佳平方逼近多项式为
17。求函数在指定区间上对于的最佳逼近多项式:
解:
若
且,则有
则法方程组为
从而解得
故关于的最佳平方逼近多项式为
若
且,则有
解:
采用复化梯形公式时,余项为
又
故
若,则
当对区间进行等分时,
故有
因此,将区间213等分时可以满足误差要求
采用复化xx公式时,余项为
又
若,则
当对区间进行等分时
故有
因此,将区间8等分时可以满足误差要求。
7。如果,证明用梯形公式计算积分所得结果比准确值大,并说明其几何意义。
解:采用梯形公式计算积分时,余项为
时间t(s)
0
0.9
1.9
3.0
3.9
5.0
距离s(m)
0
10
30
50
80
110
求运动方程。
解:
被观测物体的运动距离与运动时间大体为线性函数关系,从而选择线性方程
令
则
则法方程组为
从而解得
故物体运动方程为
20。已知实验数据如下:
19
25
31
38
44
19.0
32.3
49.0
73.3
97.8
用最小二乘法求形如的经验公式,并计算均方误差。
当时,
令
取
令
则
当时,线性插值多项式为
插值余项为
又在建立函数表时,表中数据具有5位有效数字,且,故计算中有误差传播过程。
总误差界为
4.设为互异节点,求证:
(1)
(2)
证明
(1)令
若插值节点为,则函数的次插值多项式为。
插值余项为
又
由上题结论可知
得证。
5设且求证:
解:令,以此为插值节点,则线性插值多项式为
(1)
得证。
+
得证。
14.求及。
解:
若
则
15.证明两点三次埃尔米特插值余项是
解:
若,且插值多项式满足条件
插值余项为
由插值条件可知
且
可写成
其中是关于的待定函数,
现把看成上的一个固定点,作函数
根据余项性质,有
由xx定理可知,存在和,使
即在上有四个互异零点。
根据xx定理,在的两个零点间至少有一个零点,
故在内至少有三个互异零点,
(1)若
令,则
令,则
令,则
从而解得
令,则
故成立。
令,则
故此时,
故
具有3次代数精度。
(2)若
令,则
令,则
令,则
从而解得
令,则
故成立。
令,则
故此时,
因此,
具有3次代数精度。
(3)若
令,则
令,则
令,则
从而解得
或
令,则
故不成立。
因此,原求积公式具有2次代数精度。
(4)若
令,则
令,则
令,则
故有
令,则
令,则
故此时,
=
插值余项为
6.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函数表的步长h应取多少?
解:若插值节点为和,则分段二次插值多项式的插值余项为