数值计算方法计算习题

合集下载

数值计算方法习题答案(第二版)(绪论)

数值计算方法习题答案(第二版)(绪论)

数值计算⽅法习题答案(第⼆版)(绪论)数值分析(p11页)4 试证:对任给初值x 0,0)a >的⽜顿迭代公式112(),0,1,2,......k a k k x x x k +=+= 恒成⽴下列关系式:2112(1)(,0,1,2,....(2)1,2,......k k k x k x x k x k +-=-=≥=证明:(1)(21122k k k k k k x a x x x x +-??-=+==? ??(2)取初值00>x ,显然有0>k x ,对任意0≥k ,a a x a x x a x x k k k k k ≥+-= +=+2121216 证明:若k x 有n 位有效数字,则n k x -?≤-110218,⽽()k k k k k x x x x x 288821821-=-???? ??+=-+ nn k k x x 2122110215.22104185.28--+?=??<-∴>≥ 1k x +∴必有2n 位有效数字。

8 解:此题的相对误差限通常有两种解法.①根据本章中所给出的定理:(设x 的近似数*x 可表⽰为m n a a a x 10......021*?±=,如果*x 具有l 位有效数字,则其相对误差限为()11 **1021--?≤-l a x x x ,其中1a 为*x 中第⼀个⾮零数)则7.21=x ,有两位有效数字,相对误差限为025.010221111=??≤--x x e 71.22=x ,有两位有效数字,相对误差限为025.010221122=??≤--x x e 3 2.718x =,有两位有效数字,其相对误差限为:00025.010221333=??≤--x e x ②第⼆种⽅法直接根据相对误差限的定义式求解对于7.21=x ,0183.01<-e x∴其相对误差限为00678.07.20183.011≈<-x e x 同理对于71.22=x ,有003063.071.20083.022≈<-x e x 对于718.23=x ,有00012.0718.20003.033≈<-x e x备注:(1)两种⽅法均可得出相对误差限,但第⼀种是对于所有具有n 位有效数字的近似数都成⽴的正确结论,故他对误差限的估计偏⼤,但计算略简单些;⽽第⼆种⽅法给出较好的误差限估计,但计算稍复杂。

数值计算方法复习题

数值计算方法复习题

fuxiti例1证明方程1-x-sin x=0在区间[0,1]内有一个根,使用二分法求误差不超过0.5×10-4的根要迭代多少次?证明令f(x)=1-x-sin x,∵f(0)=1>0,f(1)=-sin1<0∴f(x)=1-x-sin x=0在[0,1]有根.又f'(x)=1-c os x>0(x∈[0.1]),故f(x)=0在区间[0,1]内有唯一实根.给定误差限ε=0.5×10-4,有只要取n=14.例4选择填空题1. 设函数f(x)在区间[a,b]上连续,若满足,则方程f(x)=0在区间[a,b]一定有实根.答案:f(a)f(b)<0解答:因为f(x)在区间[a,b]上连续,在两端点函数值异号,由连续函数的介值定理,必存在c,使得f(c)=0,故f(x)=0一定有根.2. 用简单迭代法求方程f(x)=0的实根,把方程(x)=0表成x=ϕ(x),则f(x)=0的根是( )(A)y=x与y=ϕ(x)的交点(B) y=x与y=ϕ(x)交点的横坐标(C) y=x与x轴的交点的横坐标(D) y=ϕ(x)与x轴交点的横坐标答案:(B)解答:把f(x)=0表成x=ϕ(x), 满足x=ϕ(x)的x是方程的解,它正是y=x与y=ϕ(x)的交点的横坐标.3.为求方程x3―x2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( )(A)(B)(C)(D)答案:(A)解答:在(A)中故迭代发散.在(B)中,故迭代收敛.在(C)中,,故迭代收敛.在(D)中,类似证明,迭代收敛.例3填空选择题:1. 用高斯列主元消去法解线性方程组作第1次消元后的第2,3个方程分别为。

解答1. 选a21=2为主元,作行互换,第1个方程变为:2x1+2x2+3x3=3,消元得到是应填写的内容。

一、解答下列问题:1) 数值计算中,最基础的五个误差概念(术语)是 , , , , .2) 分别用 2.718281, 2.718282 作数e 的近似值 ,它们的有效位数分别有位, 位; 又取73.13≈ (三位有效数字),则≤-73.13 .3)为减少乘除法运算次数,应将算式32)1(7)1(51318---+-+=x x x y 改写成4)为减少舍入误差的影响,应将算式 9910- 改写成 5)递推公式 ⎪⎩⎪⎨⎧=-==-,2,1,110210n y y y n n如果取41.120≈=y 作计算,则计算到10y 时,误差有这个计算公式数值稳定不稳定 ?1) 绝对误差 , 相对误差 , 有效数字 , 截断误差 , 舍入误差 。

数值计算方法复习题08.6

数值计算方法复习题08.6

数值计算方法复习题 一、填空题 1、以n + 1个 整 数 点k ( k =0,1,2,…,n ) 为 节 点 的 Lagrange 插 值 基 函 数 为l k (x)( k =0,1,2,…,n ),则()_____.nkk klx ==∑()nkk klx x ==∑以n + 1个 整 数 点k ( k =1,2,…,n ,n+1) 为 节 点 的 Lagrange 插 值 基 函 数 为l k (x)( k =1,2,…,n,n+1),则111(0)__________.n n kk lk ++==∑(1)(1)!n n -+2、序列{}n=0n y ∞满足递推关系:1101,(1,2,...)n n y y n -=-=,若0y 有误差, 这个计算过程是否稳定?____________.不稳定3、42()23,[1,2,3,4,5,6]_____.f x x x f =+-=若则[1,2,3,4,5]_____.f =[1,2,3,4,5,6]0f = [1,2,3,4,5]f = 4、形如()()nbk k ak f x dx A f x =≈∑⎰的插值型求积公式,其代数精度至少可达______次,至多可达______次。

n , 2n+15、已知x =62.1341是由准确数a 经四舍五入得到的a 的近似值,试给出x 的绝对误差界_______________.41102-⨯ 6、设x 和y 的相对误差均为0.001,则xy 的相对误差约为____________.0.002二、计算题 1、(1)设3211A ⎡⎤=⎢⎥-⎣⎦,求1().cond A 解:1121513A--⎡⎤=⎢⎥⎣⎦,1111()414cond A A A -==⨯=。

(2) 设121111A ⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦,求2().cond A 解:1232,7,2,26T T A A A A σσ⎡⎤===⎢⎥⎣⎦的特征值为,122()cond A σσ===2、给出数据点: 013419156i i x y =⎧⎨=⎩(1)用012,,x x x 构造二次Lagrange 插值多项式2()L x ,并计算 1.5x =的近似值2(1.5)L 。

数值计算方法课后习题答案

数值计算方法课后习题答案

第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。

[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%,求n x 的相对误差。

[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。

3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。

[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。

4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。

(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。

数值计算方法倪勤习题答案

数值计算方法倪勤习题答案

数值计算方法倪勤习题答案数值计算方法倪勤习题答案数值计算方法是一门研究如何利用计算机进行数值计算的学科。

它在科学计算、工程计算、金融计算等领域中有着广泛的应用。

倪勤的《数值计算方法》是该领域的经典教材之一,其中的习题是帮助学生巩固所学知识的重要资源。

下面是一些数值计算方法倪勤习题的答案,供大家参考。

一、插值与拟合1. 设有下列数据点:(0, 0),(1, 1),(2, 4),(3, 9)。

试用拉格朗日插值多项式求x=2.5处的函数值。

解答:拉格朗日插值多项式的表达式为:P(x) = ∑[f(xi) * l(x)] / ∑[l(xi)]其中,l(x) = ∏[(x - xj) / (xi - xj)],i ≠ j根据给定的数据点,可以得到:l0(x) = (x - 1)(x - 2)(x - 3) / (0 - 1)(0 - 2)(0 - 3) = -x(x - 1)(x - 2) / 6l1(x) = (x - 0)(x - 2)(x - 3) / (1 - 0)(1 - 2)(1 - 3) = x(x - 2)(x - 3) / 2l2(x) = (x - 0)(x - 1)(x - 3) / (2 - 0)(2 - 1)(2 - 3) = -x(x - 1)(x - 3) / 2l3(x) = (x - 0)(x - 1)(x - 2) / (3 - 0)(3 - 1)(3 - 2) = x(x - 1)(x - 2) / 6代入公式,得到:P(x) = 0 * l0(x) + 1 * l1(x) + 4 * l2(x) + 9 * l3(x)= -x(x - 1)(x - 2) / 6 + 4x(x - 1)(x - 3) / 2 + 9x(x - 1)(x - 2) / 6= -x(x - 1)(x - 2) / 6 + 2x(x - 1)(x - 3) + 3x(x - 1)(x - 2) / 2= x^3 - 3x^2 + 3x将x=2.5代入上式,得到:P(2.5) = 2.5^3 - 3 * 2.5^2 + 3 * 2.5 = 2.375因此,当x=2.5时,函数值为2.375。

数值计算方法习题答案

数值计算方法习题答案

数值分析(p11页)4 试证:对任给初值x 0,0)a >的牛顿迭代公式 112(),0,1,2,......ka k kx x x k +=+=恒成立下列关系式:2112(1)(,0,1,2,....(2)1,2,......kk kx k x x k x k +-=-=≥=证明: (1)(2211222k k k k k k k kx a x ax x x x x +⎫⎛-+=+-==⎪ ⎝⎭(2) 取初值0>x ,显然有0>kx,对任意0≥k ,a a x a x x a x x k k k k k ≥+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛+=+2121216 证明:若kx 有n 位有效数字,则n kx-⨯≤-110218,而()k k k kk x x x x x 288821821-=-⎪⎪⎭⎫⎝⎛+=-+nnk k x x 2122110215.22104185.28--+⨯=⨯⨯<-∴>≥Θ1k x +∴必有2n 位有效数字。

8 解:此题的相对误差限通常有两种解法. ①根据本章中所给出的定理:(设x 的近似数*x 可表示为mna a a x 10 (02)1*⨯±=,如果*x 具有l 位有效数字,则其相对误差限为()11**1021--⨯≤-l a xx x ,其中1a 为*x 中第一个非零数) 则7.21=x ,有两位有效数字,相对误差限为 025.010221111=⨯⨯≤--x x e 71.22=x ,有两位有效数字,相对误差限为 025.010221122=⨯⨯≤--x x e 3 2.718x =,有两位有效数字,其相对误差限为: 00025.010221333=⨯⨯≤--x e x ②第二种方法直接根据相对误差限的定义式求解对于7.21=x ,0183.01<-e x∴其相对误差限为00678.07.20183.011≈<-x e x 同理对于71.22=x,有 003063.071.20083.022≈<-x e x对于718.23=x,有00012.0718.20003.033≈<-x e x备注:(1)两种方法均可得出相对误差限,但第一种是对于所有具有n 位有效数字的近似数都成立的正确结论,故他对误差限的估计偏大,但计算略简单些;而第二种方法给出较好的误差限估计,但计算稍复杂。

数值计算方法课后习题答案(李庆扬等) (修复的)

数值计算方法课后习题答案(李庆扬等) (修复的)

,。

第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。

[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x xx x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%,求n x 的相对误差。

[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。

3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。

[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。

4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。

(1)*4*2*1x x x ++;[解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。

数值计算方法习题@中科大

数值计算方法习题@中科大

上机作业题程序1 编写程序计算下列向量范数∑==+++=ni i n x x x x X1211||||||||∑==+++=ni inxx x x X12222212或),(2X X X=|}{|max |}|,|,||,max{|121i n i n x x x x X≤≤∞==输入:向量的阶数n ,向量X 的元素 输出:向量范数程序2 编写程序计算下列矩阵范数||||max ||A a j nij i n111=≤≤=∑Aa i ni j j n∞≤≤==∑max ||11||||||A a E i n ij j n=⎛⎝⎫⎭⎪⎪⎪==∑∑12112输入:方程组的阶数n ,矩阵A 的元素 输出:矩阵范数程序3 编写程序计算如下级数,要求误差小于1.06e -11()()k x k k x ∞=ψ=+∑并计算0.1,0.2,,1.0,10.0,20.0,,300.0x = 的值 输入:无 输出:,()x x ψ的值程序4 下面给出美国从1920年到1970年的人口表:的人口。

在1910年的实际人口约为91772000,请判断插值计算得到的1965年和2002年的人口数据准确性是多少?程序5 数据同上表,用牛顿插值估计: (1)1965年的人口数; (2)2002年的人口数。

程序6 数据同上表,用自然样条函数预测在1910,1965和2002年的人口数。

请比较以上三种方法所求值的效果。

那一种方法最优?程序7 给定1+n 个插值节点,构造n 次拉格朗日插值多项式,并计算)(x f 。

输入:插值节点数n ,插值点{}n i x f x i i ,,2,1,0)(, =,;要计算的函数点x 。

输出:)(x L n 的值。

程序8 对函数21() , [5,5](1)f x x x =∈-+以如下两组节点为插值节点构造插值函数,(1)105,0,1,i x i i N N =-+= (2)215cos ,0,1,22i i x i NN π+=-=+并用式子55max ()()max ()(),5,0,10010i i i x iif x p x f y p y y i -≤≤-≈-=-=对5,10,20,40N =估计两组节点的误差输入:无。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1.已知ln(
2.0)=0.6931;ln(2.2)=0.7885,ln(2.3)=0.8329, 试用线性插值和抛物插值计算.ln2.1的值并估计误差(牛顿插值和拉格朗日插值)
2.已知函数y=sinx 的数表如下,分别用前插和后插公式计算sin0.57891的值,并估算误差。

i x
0.4 0.5 0.6 0.7 )(i x f
0.38942
0.47943
0.56464
0.64422
3. 已知
i x
-2 -1 0 1 2 )(i x f
4
2
1
3
5
求)(x f 的二次拟合曲线)(2x p ,并求)0(f '的近似值。

4. 数值积分公式形如
⎰'+'++=≈1
)1()0()1()0()()(f D f C Bf Af x S dx x xf 试确定参数D C B A ,,,使公式代数
精度尽量高;(2)设]1,0[)(4
C x f ∈,推导余项公式⎰-=1
)
()()(x S dx x xf x R ,并估计
误差。

5. 已知数值积分公式为:
)]
()0([)]()0([2)(''20
h f f h h f f h
dx x f h
-++≈⎰
λ,试确定积分公式中的参数
λ,使其代数精确度尽量高,并指出其代数精确度的次数。

6. 用复化Simpson 公式计算积分
()⎰
=1
0sin dx x x I 的近似值,要求误差限为
5105.0-⨯。

7. 已知012113,,4
2
4
x x x ===,给出以这3个点为求积节点在[]0.1上的插值型求积公
式。

8. 给出 900,cos ≤≤x x 的函数表,步长
)60/1(1='=h ,若函数具有5位有
效数字,研究用线性插值求x cos 近似值时的总误差界。

9. 求一个次数不高于4次的多项式)(x P ,使它满足0)0()0(='=P P ,
1)1()1(='=P P ,1)2(=P 。

10. 单原子波函数的形式为bx
ae y -=,试按照最小二乘法决定参数a 和b ,已
知数据如下:
X 0 1 2 4 y
2.010 1.210 0.740 0.450
11. 分别用梯形公式和辛普森公式计算下列积分:⎰
+1
02
4dx x x。

并估算误差。

12. 用矩阵的克劳特和克利特尔三角分解法求解方程组:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛7173530103421101002014321x x x x。

相关文档
最新文档