当前位置：文档库 › 李庆扬-数值分析第五版第7章习题答案(20130824)

李庆扬-数值分析第五版第7章习题答案(20130824)

第7章复习与思考题

习题

数值分析习题与答案

第一章绪论习题一?1．设x＞0，x＊的相对误差为δ,求f(ｘ)=ln x的误差限。解:求lnx的误差极限就是求f(x)=lnx的误差限，由公式(１.2.４）有已知x*的相对误差满足，而 ,故即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字，并给出其误差限与相对误差限。解:直接根据定义和式（１.2.２)(1.２.3）则得?有５位有效数字，其误差限,相对误差限有２位有效数字，有5位有效数字， 3．下列公式如何才比较准确？ (１)?(2）解:要使计算较准确，主要是避免两相近数相减，故应变换所给公式。

（1)?(2） 4.近似数x*=０.０３１0,是 3 位有数数字。 5.计算取,利用 :式计算误差最小。四个选项: 第二、三章插值与函数逼近习题二、三 1. 给定的数值表用线性插值与二次插值计算ln０．5４的近似值并估计误差限. 解：仍可使用n＝1及n＝２的Laｇrange插值或Neｗｔo ｎ插值，并应用误差估计（5.8)。线性插值时,用0．５及0．6两点，用Newton插值??误差限，因,

故? 二次插值时，用0.5，0．6,０.7三点,作二次Ｎewton插值 ?误差限，故? 2. 在－4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过,函数表的步长ｈ应取多少? 解：用误差估计式（５．8), ?令因?得 3. 若,求和.

解：由均差与导数关系 ?于是 4. 若互异,求的值,这里p≤n+1. 解:,由均差对称性可知当有?而当P＝n ＋1时 ?于是得 5. 求证. 解:解:只要按差分定义直接展开得 ? 6. 已知的函数表

数值分析习题

习题1 1．填空题 (1) 为便于算法在计算机上实现，必须将一个数学问题分解为的运算； (2) 在数值计算中为避免损失有效数字，尽量避免两个数作减法运算；为避免误差的扩大，也尽量避免分母的绝对值分子的绝对值； (3) 误差有四大来源,数值分析主要处理其中的和； (4) 有效数字越多,相对误差越； 2. 用例1.4的算法计算10,迭代3次,计算结果保留4位有效数字. 3. 推导开平方运算的误差限公式,并说明什么情况下结果误差不大于自变量误差. 4. 以下各数都是对准确值进行四舍五入得到的近似数,指出它们的有效数位、误差限和相对误差限. 95123450304051104000003346087510., ., , ., .x x x x x -==?===? 5. 证明1.2.3之定理1.1. 6. 若钢珠的的直径d 的相对误差为1.0%,则它的体积V 的相对误差将为多少。（假定钢珠为标准的球形） 7. 若跑道长的测量有0.1%的误差,对400m 成绩为60s 的运动员的成绩将会带来多大的误差和相对误差. 8. 为使20的近似数相对误差小于0.05%,试问该保留几位有效数字. 9. 一个园柱体的工件,直径d 为10.25±0.25mm,高h 为40.00±1.00mm,则它的体积V 的近似值、误差和相对误差为多少． 10 证明对一元函数运算有 r r xf x f x k x k f x εε'≈= () (())(),() 其中并求出157f x x x ==()tan ,.时的k 值，从而说明f x x =()tan 在2 x π ≈时是病态问题． 11. 定义多元函数运算 1 1 1,,(),n n i i i i i i S c x c x εε====≤∑∑其中求出S ε()的表达式，并说明i c 全为正数时，计算是稳定的，i c 有正有负时，误差难以控制． 12. 下列各式应如何改进,使计算更准确：

李庆扬数值分析第五版习题复习资料清华大学出版社

第一章绪论 1．设0x >,x 的相对误差为δ，求ln x 的误差。解：近似值* x 的相对误差为* **** r e x x e x x δ-= = = 而ln x 的误差为()1 ln *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 -=Q , 1 ||n p x nx C n n -?∴== 又((*))(*)r p r x n C x εε≈?Q 且(*)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．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中**** 1234,,,x x x x 均为第3题所给的数。解：

*4 1* 3 2* 13* 3 4* 1 51()1021()1021()1021()1021()102 x x x x x εεεεε-----=?=?=?=?=? *** 124***1244333 (1)()()()() 1111010102221.0510x x x x x x εεεε----++=++=?+?+?=? *** 123*********123231132143 (2)() ()()() 111 1.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈ ** 24**** 24422 *4 33 5 (3)(/) ()() 11 0.0311056.430102256.43056.430 10x x x x x x x εεε---+≈ ??+??= ?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？解：球体体积为343 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ===g g (*)(*)3(*)r p r r V C R R εεε∴≈=g 又(*)1r V ε=Q

数值分析第五版全答案chap1

第一章绪论 1．设0x >,x 的相对误差为δ，求ln x 的误差。解：近似值*x 的相对误差为* **** r e x x e 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 =,*5 7 1.0.x =? 解：*1 1.1021x =是五位有效数字； *20.031x =是二位有效数字； *3385.6x =是四位有效数字； *456.430x =是五位有效数字； *57 1.0.x =?是二位有效数字。 4．利用公式(2.3)求下列各近似值的误差限：(1) *** 124x x x ++,(2) ***123x x x ,(3) **24 /x x . 其中****1234,,,x x x x 均为第3题所给的数。解：

*4 1*3 2*13*3 4*1 51 ()102 1()102 1()102 1()102 1()102x x x x x εεεεε-----=?=?=?=?=? ***124***1244333 (1)() ()()() 111101010222 1.0510x x x x x x εεεε----++=++=?+?+?=? ***123*********123231132143 (2)() ()()() 1111.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈ **24****24422 *4 33 5 (3)(/)()() 110.0311056.430102256.43056.430 10x x x x x x x εεε---+≈??+??=?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？解：球体体积为343V R π= 则何种函数的条件数为 2 3 '4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε=

(完整版)数值分析第7章答案

第七章非线性方程求根一、重点内容提要（一）问题简介求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为函数()f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在 (a,b)内仅有一个根.令00,a a b b ==,计算0001 ()2x a b =+和0()f x .若0()0f x =则*x x =,结束计算;若00()()0f a f x >,则令10,1a x b b ==,得新的有根区间11[,]a b ;若 00()()0 f a f x <,则令 10,10 a a b x ==,得新的有根区间 11[,]a b .0011[,][,]a b a b ?,11001()2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得出新的有根区间22[,] a b ,如此反复进行,可得一有根区间套 1100...[,][,]...[,] n n n n a b a b a b --????

数值分析习题

第一章绪论习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5 105.0-?,那么近似数0.003400有几位有效数字？（有效数字的计算） 2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算） 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5* =，已知 cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差限与相对误差限。（误差限的计算） 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大？（函数误差的计算） 8 设?-=1 1 dx e x e I x n n ，求证：（1）)2,1,0(11 =-=-n nI I n n （2）利用（1）中的公式正向递推计算时误差逐步增大；反向递推计算时误差逐步减小。（计算方法的比较选择）

第二章插值法习题主要考察点：拉格朗日插值法的构造，均差的计算，牛顿插值和埃尔米特插值构造，插值余项的计算和应用。 1 已知1)2(,1)1(,2)1(===-f f f ，求)(x f 的拉氏插值多项式。（拉格朗日插值） 2 已知9,4,10=== x x x y ，用线性插值求7的近似值。（拉格朗日线性插值） 3 若),...1,0(n j x j =为互异节点，且有 ) ())(())(()())(())(()(11101110n j j j j j j j n j j j x x x x x x x x x x x x x x x x x x x x x l ----------= +-+- 试证明 ),...1,0()(0 n k x x l x n j k j k j =≡∑=。（拉格朗日插值基函数的性质） 4 已知352274.036.0sin ,333487.034.0sin ,314567.032.0sin ===，用抛物线插值计算3367.0sin 的值并估计截断误差。（拉格朗日二次插值） 5 用余弦函数x cos 在00=x ，4 1π =x ，2 2π = x 三个节点处的值，写出二次拉格朗日插值多项式, 并近似计算6 cos π 及其绝对误差与相对误差，且与误差余项估计值比较。（拉格朗日二次插值） 6 已知函数值212)6(,82)4(,46)3(,10)1(,6)0(=====f f f f f ，求函数的四阶均差 ]6,4,3,1,0[f 和二阶均差]3,1,4[f 。（均差的计算） 7 设)())(()(10n x x x x x x x f ---= 求][1,0p x x x f 之值，其中1+≤n p ，而节点 )1,1,0(+=n i x i 互异。（均差的计算） 8 如下函数值表建立不超过三次的牛顿插值多项式。（牛顿插值多项式的构造） 9求一个次数小于等于三次多项式)(x p ，满足如下插值条件：2)1(=p ，4)2(=p ， 3)2(='p ，12)3(=p 。（插值多项式的构造）

东南大学_数值分析_第七章_偏微分方程数值解法

第七章偏微分方程数值解法 ——Crank-Nicolson 格式 ****（学号） *****（姓名）上机题目要求见教材P346,10题。一、算法原理本文研究下列定解问题(抛物型方程) 22(,) (0,0)(,0)() (0) (0,)(), (1,)() (0)u u a f x t x l t T t x u x x x l u t t u t t t T ?αβ???-=<<≤≤???? =≤≤??==<≤?? (1) 的有限差分法，其中a 为正常数，,,,f ?αβ为已知函数，且满足边界条件和初始条件。关于式(1)的求解，采用离散化方法，剖分网格，构造差分格式。其中，网格剖分是将区域{}0,0D x l t T =≤≤≤≤用两簇平行直线 (0) (0)i k x x ih i M t t k k N τ==≤≤?? ==≤≤? 分割成矩形网格，其中,l T h M N τ==分别为空间步长和时间步长。将式(1)中的偏导数使用不同的差商代替，将得到不同的差分格式，如古典显格式、古典隐格式、Crank-Nicolson 格式等。其中，Crank-Nicolson 格式具有更高的收敛阶数，应用更广泛，故本文采用Crank-Nicolson 格式求解抛物型方程。 Crank-Nicolson 格式推导：在节点(,)2 i k x t τ +处考虑式(1)，有 22(,)(,)(,)222 i k i k i k u u x t a x t f x t t x τττ??+-+=+?? (2) 对偏导数 (,)2 i k u x t t τ ?+?用中心差分展开 []2311+13 1(,)(,)(,)(,) ()224k k i k i k i k i i k i k u u x t u x t u x t x t t t t ττηητ++??+=--<

数值分析第七章非线性方程求根习题答案

第七章非线性方程求根（一）问题简介求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为函数() f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在(a,b)内仅有一个根.令00,a a b b ==,计算0001()2x a b =+和 0()f x .若0()0f x =则*x x =,结束计算;若 00()()0 f a f x >,则令 10,1a x b b ==,得新的有根区间 11[,] a b ;若 00()()0 f a f x <,则令 10,10a a b x ==,得新的有根区间11[,]a b .0011[,][,]a b a b ?,11001 () 2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得出新的有根区间22[,]a b ,如此反复进行,可得一有根区间套 1100...[,][,]...[,] n n n n a b a b a b --???? 且110011 *,0,1,2,...,()...() 22n n n n n n a x b n b a b a b a --<<=-=-==-. 故 1 lim()0,lim lim ()* 2n n n n n n n n b a x a b x →∞→∞→∞-==+=

数值计算方法第七章习题 2013

计算方法第七章习题复习与思考题 1．设f ∈C [a , b ]，写出三种常用范数2 1 f f 及∞ f 。 2．f , g ∈C [a , b ]，它们的内积是什么？如何判断函数族{? 0, ? 1, …, ? n }∈C [a , b ]在[a ,b ]上线性无关？ 3．什么是函数f ∈C [a , b ]在区[a , b ]上的n 次最佳一致逼近多项式？ 4．什么是f 在[a , b ] 上的n 次最佳平方逼近多项式？什么是数据{}m i f 0的最小二乘曲线拟合？ 5．什么是[ a , b ]上带权ρ (x )的正交多项式？什么是[ -1, 1 ]上的勒让德多项式？它有什么重要性质？ 6．什么是切比雪夫多项式？它有什么重要性质？ 7．用切比雪夫多项式零点做插值得到的插值多项式与拉格朗日插值有何不同？ 8．什么是最小二乘拟合的法方程？用多项式做拟合曲线时，当次数n 较大时为什么不直接求解法方程？ 9．哪种类型函数用三角插值比用多项式插值或分段多项式插值更合适？ 10．判断下列命题是否正确？（1）任何f (x ) ∈C [a , b ]都能找到n 次多项式P n (x ) ∈ H n ，使| f (x ) - P n (x ) | ≤ ε ( ε 为任给的误差限)。（2）n n H x P ∈)(* 是f (x )在[ a , b ]上的最佳一致逼近多项式，则)()(lim * x f x P n n =∞ →对 ],[b a x ∈?成立。（3）f (x ) ∈C [a , b ]在[a , b ]上的最佳平方逼近多项式P n (x ) ∈ H n 则)()(lim x f x P n n =∞ →。（4）)(P ~ x n 是首项系数为1的勒让德多项式，Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 ? ? --1 1 21 1 2d )(d )](P ~ [x x Q x x n n 。（5）)(T ~ x n 是[-1 , 1]上首项系数为1的切比雪夫多项式。Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 .)(max )(~ max 1 11 1x Q x T n x n x ≤≤-≤≤-≤ （6）当数据量很大时用最小二乘拟合比用插值好。

数值分析第七章上机题

数值分析第七章计算机实习题写一程序实现下面问题的牛顿算法——求解方程组： ?? ???=--=-+.0)1sin(,18)7)(3(12321x e x x x 源程序如下： function [x,it,hist] = newton2(x0,f,g,maxit,tol) % Newton method for eqation systerm % INPUTS: % x0 initial point % f function % g gradient % maxit maximum iteration % tol tolerance for convergence % OUTPUTS: % x solution % it iteration % hist history of iteration format long ; if nargin<5, tol = 1e-7; if nargin<4, maxit = 100; if nargin<3, error('too few input!!'); end end end flag = 1; x0 = [0;0]; x = x0; hist = x; it = 0; for k = 1:maxit, x = x0 - feval(g,x0(1),x0(2))\feval(f,x0(1),x0(2)); if norm(x0-x)>=tol, x0 = x; else fprintf('\nNewton Iteration successes!!\n') return end it = it + 1;

hist = [hist x]; end flag = 0; fprintf('\nNewton Iteration fails!!\n'); 在命令窗口输入： >>f = inline('[(x1+3)*(x2^3-7)+18;sin(x2*exp(x1)-1)]','x1','x2'); >>g = inline ('[x2^3-7,3*x2^2*(x1+3);x2*exp(x1)*cos(x2*exp(x1)-1),exp(x1)*cos(x2*exp(x1)-1)]','x1','x2'); >> [x,it,hist] = newton2([0;0],f,g) 得到如下运行结果： >> [x,it,hist] = newton2([0;0],f,g) Newton Iteration successes!! x = -0.000000000000000 1.000000000000000 it = 5 hist = 0 -0.428571428571429 -0.141348392468100 -0.002875590925150 0.000000056935424 -0.000000000000101 0 1.557407724654902 1.087738055836075 1.001269946612821 1.000000431005363 1.000000000000127 由以上运行结果可知：该方程组采用牛顿迭代法迭代5步可到足够精度，解为??? ? ??=10x .

数值分析第五版答案

第一章绪论 p19 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 ε∴≈ 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？解：球体体积为343 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又 (*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε= ?≈ 7．求方程2 5610x x -+=的两个根，使它至少具有427.982 =）。解：2 5610x x -+= ，故方程的根应为1,228x =故 128 2827.98255.982x = ≈+= 1x ∴具有5位有效数字 211 280.0178632827.98255.982 x =-= ≈ =≈+ 2x 具有5位有效数字

9．正方形的边长大约为了100cm ，应怎样测量才能使其面积误差不超过2 1cm ？解：正方形的面积函数为2 ()A x x = p7 当*100x =时，若(*)1A ε≤, 则21 (*)102 x ε-≤ ? 故测量中边长误差限不超过0.005cm 时，才能使其面积误差不超过2 1cm 第二章插值法p48 1．当1,1,2 x =-时，()0,3,4f x =-, 分别用单项式基底、拉格朗日基底、牛顿基底求() f x 的二次插值多项式。解： 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23537623 l x l x x x x x x x =-+=---+-+=+- 2．给出()ln f x x =的数值表用线性插值及二次插值计算的近似值。解：由表格知，

数值分析习题集及答案

数值分析习题集（适合课程《数值方法A》和《数值方法B》）长沙理工大学第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2％,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少? 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差? 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求? 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程稳定吗? 12.计算,取,利用下列等式计算,哪一个得到的结果最好? 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足第二章插值法 1.根据定义的范德蒙行列式,令证明是n次多项式,它的根是,且 .

2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式. 3. 4.给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界. 5.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函数表的步长应取多少? 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误

数值分析第五版全答案

第四章数值积分与数值微分 1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度： 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解：求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。（1）若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =，则 1012h A A A -=++ 令()f x x =，则 110A h A h -=-+ 令2 ()f x x =，则 3 221123 h h A h A -=+ 从而解得 01 1431313A h A h A h -?=?? ?=?? ?=?? 令3 ()f x x =，则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 故 101()()(0)()h h f x dx A f h A f A f h --=-++? 成立。令4 ()f x x =，则

数值分析第五版答案

第一章绪论 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．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234 ,,,x x x x 均为第3题所给的数。解： *4 1* 3 2* 13* 3 4* 1 51 ()1021()1021()1021()1021()102 x x x x x εεεεε-----=?=?=?=?=? *** 124***1244333 (1)()()()() 1111010102221.0510x x x x x x εεεε----++=++=?+?+?=? *** 123*********123231132143 (2)() ()()() 111 1.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈

** 24**** 24422 * 4 33 5 (3)(/) ()() 11 0.0311056.430102256.43056.430 10x x x x x x x εεε---+≈ ??+??= ?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？解：球体体积为34 3 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε=?≈ 6．设028Y = ，按递推公式1n n Y Y -= （n=1,2,…）计算到100Y 27.982≈（5位有效数字），试问计算100Y 将有多大误差？解：1n n Y Y -=- 10099Y Y ∴=- 9998Y Y = 9897Y Y =-…… 10Y Y =- 依次代入后，有1000100Y Y =- 即1000Y Y = 27.982, 100027.982Y Y ∴=-

第七章习题答案数值分析.doc

第七章习题解答 2、试确定系数a ，b 的值使 220[()cos ]ax b x dx π +-?达到最小解：设22 0(,)[()cos ]I a b ax b x dx π=+-? 确定a ，b 使(,)I a b 达到最小，必须满足 0,0I I a b ??==?? 即 3222222000022222000012[cos ]0cos 248212[cos ]0cos 8 2a b ax b x xdx a x dx b xdx x xdx a b ax b x dx a xdx b dx xdx πππππππππππππ???+=-+-=+=??????????????+=+-=+=??????????????解得：0.6644389, 1.1584689a b ≈-≈ 5、试用Legendre 多项式构造()f x x =在[-1, 3]上的二次最佳平方逼近多项式解：作变量代换，将区间[-1, 3]变为[-1, 1]，令21x t =+，即12x t -= 则()()(21)21(11)F t f x f t t t ==+=+-≤≤ 对()F t 利用Legendre 多项式求其在}{21,,span t t 上的最佳平方逼近多项式 20()()j j j S t C P t ==∑，其中11(,)21()()(0,1,2)(,)2j j j j j P f j C F t P t dt j P P -+===? 则有： 1121012 1121112122121212 15[(21)(21)]24311[(21)(21)]285(31)(31)45[(21)(21)]22264C t dt t dt C t dt t dt t t C t dt t dt ---------=--++==--++=--=--++=?????? 01151145()()()()4864 S t P t P t P t ∴=++则()f x 在[-1, 3]上的最佳二次逼近多项式 2* 0121511114514051061081()()()()()()2421226421536x x x x x x S t S t S P P P -----+===++= 7、确定一条经过原点的二次曲线，使之拟合下列数据

数值分析复习题要答案

第一章 1、ln2＝0.69314718…，精确到 10－ 3 的近似值是多少？解精确到 10－ 3＝0.001，即绝对误差限是 ε＝0.05％，故至少要保留小数点后三位才可以。 ln2≈0.693。 2、设115.80,1025.621≈≈x x 均具有5位有效数字，试估计由这些数据计算21x x ， 21x x +的绝对误差限解：记126.1025, 80.115x x ==%% 则有11232 4 1110, | 102|||2 x x x x --≤?-≤?-%% 所以 121212121212211122||||||||||||x x x x x x x x x x x x x x x x x x -=-+-+≤--%%%%%%%%% 3411 80.11610 6.10102522 0.007057-==??+≤?? 12121122 4 3|()|||11|10100.0005522 |x x x x x x x x --≤≤?+?=+-+-+-%%%% 3、一个园柱体的工件,直径d 为10.25±0.25mm,高h 为40.00±1.00mm,则它的体积V 的近似值、误差和相对误差为多少。解： ()() 22222222 4314210254000000330064 221025400002510251002436444 3300624362436 0073873833006 ,.....; ()()()......, ..().()..% .r d h V d h V mm d h V dh d d h V mm V V V πππππεεεεε= ≈=??===+=???+?==±====第二章： 1、分别利用下面四个点的Lagrange 插值多项式和Newton 插值多项式N 3(x )，计算L 3(0.5)及N 3(-0.5) 解：(1)先求332211003)()()()()(y x l y x l y x l y x l x L +++= （1分） =----+---+=------=)12)(02)(12() 1)(0)(1())()(())()(()(3020103210x x x x x x x x x x x x x x x x l x x x )1)(1(6 1-+-，（2分）

数值分析第五版_李庆扬

数值分析第五版_李庆扬一、课程基本信息课程中文名称：数值分析课程英文名称： Numerical Analysis 课程类别：专业基础课开课学期：秋适用专业：信息与计算科学；应用数学总学时： 86学时（其中理论课56学时，上机实习30学时）总学分： 5（理论课3学分；上机实习2学分）预修课程（编号）：数学分析，高等代数，常微分方程课程简介：本课程是大学本科信息与计算科学和应用数学专业的一门基础课，也是工科研究生的必修课。本课程的主要内容是研究各种数学问题的数值计算方法的设计、计算误差分析以及有关理论和具体实现的一门数学课程。是应用数学的重要分支之一。建议教材：《计算方法》（二版）（邓建中、刘之行），西安，西安交通大学出版社，2001 参考书： [1]数值分析学习指导，关治编，出版社：清华大学出版社，出版时间：2008年； [2]数值分析，何汉林，梅家斌，科学出版社，2007年； [3]《数值计算引论》白峰杉高等教育出版社 2005年 [4]《数值分析》（第五版）李庆扬易大义等清华大学出版社 2008年 [5]Numerical Analysis，R.Kress，世界图书出版公司2003 6、数值分析学习辅导习题解析，李宏、徐长发编，华中科技大学出版社，2001年。二、理论课程教育目标通过本课程的教学使学生能了解现代科学计算中常用的数值计算方法及其基本理论，系统掌握数值分析的基本概念和分析问题、解决问题的基本方法，为运用数值分析的理论知识并为掌握更复杂的现代计算方法打好。三、理论教学内容与要求（含学时）第一章：计算方法的一般概念（4学时）本章教学内容：理解计算方法的意义、研究内容与方法，理解并掌握误差的概念（包括误差的来源、绝对误差、相对误差），掌握有效数字及舍入误差对计算的影响。第二章：解线性方程组的直接法（8学时）

李庆扬-数值分析第五版第7章习题答案(20130824)

数值分析习题与答案

数值分析习题

李庆扬数值分析第五版习题复习资料清华大学出版社

数值分析第五版全答案chap1

(完整版)数值分析第7章答案

最新数值分析课程第五版课后习题答案(李庆扬等)1

数值分析习题

东南大学_数值分析_第七章_偏微分方程数值解法

数值分析第七章非线性方程求根习题答案

数值计算方法第七章习题 2013

数值分析第七章上机题

数值分析第五版答案

数值分析习题集及答案

数值分析第五版全答案

数值分析第五版答案

第七章习题答案数值分析.doc

数值分析复习题要答案

数值分析第五版_李庆扬