（）
（A）有关；（B）无关；（C）不一定。
3．幂法是用来求矩阵（）特征值及特征向量的迭代法。
（A）按模最大；（B）按模最小；
（C）任意一个；（D）所有的。
4．解代数线性方程组的松弛法收敛的必要条件是（）
（A）1
≤ω；
0<
0<
<ω；（B）1
（C）2
≤ω。
0≤
0<
<ω；（D）2 5．反幂法是用来求矩阵（）特征值及特征向量的迭代法。
（A）按模最大；（B）按模最小；
（C）任意一个；（D）所有的。
四、计算题
1．用简单迭代法（雅可比迭代法）
解线性方程组
???
??-=+--=+-=+8
41353321
32131
x x x x x x x x
取
(0)(0,0,0)T
=x ，列表计算三次，保
留三位小数。
2．用高斯—赛德尔迭代法解线性方程组
1
3123123
353148x x x x x x x x +=??
-+=-??-+=-?
取(0)
(0,0,0)T
=x
，列表计算三次，保留
三位小数。 3．用幂法求矩阵
??
??
?
?????---=210121004A 按模最大特
征值及相应特征向量，列表计算三次，取(0)
(1,1,1)T
=x ，保留两位小数。
4*．取46.1=ω，用松弛法解线性方程组
??????
?=+-=-+-=-+-=-0
412021
24343232121x x x x x x x x x x
取
(0)(0,0,0)T
=x ，列表计算三次，保留
三位小数。
5*．用雅可比方法求实对称矩阵
??
??
?
?????=110121014A 的特征值及相应特征向量（按四位小数计算，1.0=ε）。
6*．用QR 算法求矩阵??
??
??????=410131012A 的全部
特征值。
练 习 题 五
一、是非题
1.在求插值多项式时，插值多项式的次
数越高，误差越小。 ( ) 2.
120102()()()()
x x x x x x x x ----表示节点处的二次插值基函
数。 ( )
3.牛顿插值多项式的优点是：在计算时，高一级的插值多项式可利用前一次插值的结果。 ( )
4.在拉格朗日插值中，插值节点0
1
,,,n
x x x L 必须按顺序排列。 ( ) 5.利用等距节点的牛顿插值公式计算附近的，用后插公式。 ( )
二、填空题
1.已知，则三次插值基函数=_____________________。 +1个节点的拉格朗日插值基函数的和
∑==n
i i x l 0______
)(。
3.已知4
)(x x f =，取节点(0,1,2,
k
x
k k ==…），用线
性插值求)1.2(f 的近似值，其计算公式
1(2.1)(2.1)________________
f P ≈=。
插值不仅要求插值函数和被插值函数在节点取已知函数值而且取已知导数值。 5.已知(1)2,(0)1,(2)3,f f f -===则
=
-]0,1[f __________________，
=
]2,0[f ___________，[1,0,2]__________f -=，牛顿二
次插值多项式
2()N x =
_____________________________。
三、选择题
1．函数101
x x x x --表示线性插值( )点的基函数. (A) ； (B) ； (C) (D) 。
2．过点)4,2(),3,0(),1,1(-的二次插值多项式)(2x p 中的系数为( ).
(A) – (B) (C) 2 (D) -2
3．给定互异的节点01,,,,n x x x L 是以它们为插值节点的插值多项式，则是一个( ).
(A). n +1次多项式 (B). n 次多项式
(C). 次数小于n 的多项式 (D). 次
数不超过n 的多项式
4．差商,7503)(6
99
x x x x f -+-=(]2,,2,2,1[100
2
=Λf ) (A) 0 (B) -3 (C) 50 (D) -7 5．对于次数不超过n 的多项式
为次插值多项式它的)(),(x p n x f ( ).
(A) 任意n 次多项式 (B) 任意不超过n 次的多项式
(C) 本身 (D) 无法确定
四、计算题
1.已知,4)2(,3)1(,2)1(-===-f f f 求的牛顿插值多项式)(2x N ，及)5.1(f 的近似值，取三位小数。
2.证明：若f (x )二阶连续可微，则对于f (x )的以10,x x 为节点的一次插值多项式，插值误差
012
101()()()()max 8x x x x x f x P x f x ≤≤-''-≤
3.设12)(4
-+=x x x f ，利用拉格朗日插值余项求以-1，0，1，2为插值节点的
三次插值多项式。
4*
．已知函数)(x f y =的数据010)1(,)2(,)1(m f y f y f ='==，用基函数法求 f
(x )的二次插值多项式)(2x H 使20212
0(1),(2),(1)H y H y H m '===.
5*
．要给出()x
f x e =在区间[-2,2]上的等距节
点函数表，用分段三次Hermite 插值求的近似值x
e ，要使误差不超过，问函数
表的步长h应为多少
f(x)函数表
(1)求f (x)的二次插值多项式；
(2)用反插值求x，使f (x)=0。
练习题六
一、判断题
1．在等距节点的情况下，才能计算函数的差分。( ) 2．向前差分与向后差分不存在等量关系。( ) 3．已知观察值),(i i y x(,2,1,0 i…,n)，用最小二乘法求得的拟合多项式其次数为n次。() 4．利用最小二乘原理对一组数据找出合适的数学公式来拟合，首先应确定公式的类型。 ( ) 5．数据拟合的步骤首先是建立正规方程组。()
二、填空题
1．已知某函数的二阶向前差分为,则其二阶向后差分为_______。
2．利用牛顿前插公式计算某点的近似值，应首先确定公式中的t ，其计算公式为t =____________。
3．已知函数i i y x n b a x f y 处的函数值个节点上的在1],[)(+=，则其三次样条插值函数满足的条件为)(x s ________________________。
4．已知),(i i y x (,2,1=i …,30)，其线性拟合的正规方程组为_________。 5．用形如b ax x
y +=
的非线性拟合数据),(i i y x 做变换_____________后为线性拟
合=x b a +。 三．选择题
1. ( )是利用函数的值求自变量的值。 (A) 三次样条插值 (B) 反插值
(C) 分段插值 (D) 爱尔米特插值
2．记
*
,1,2,,i i i y y i n δ=-=L ,最小二乘法原理要求下列哪个为最小 ( ) (A) i
n
i δ≤≤1max (B)∑=n
i i
1
δ
(C) ∑=n
i i
1
2δ
(D)∑=n
i i
1
δ
3．当线性方程组满足 ( )时称为超定方程组。
(A) (A)未知数的个数等于方程的个数 (B) (B)未知数的个数大于方程的个数 (C) (C)未知数的个数小于方程的个数 (D) (D)未知数的个数与方程的个数大小任意
4．是超定方程组A =x b 的最小二乘解的充分必要条件是( ).
(A)
*T T A A A =x x b 是的解
(B)
*T T AA A =x x b 是的解
(C)
*T T A =x x b 是的解 (D) 三者都不对 5．勒让德多项式21d ()[(1)]2!d n n
n n n
P x x n x =-是 ( )
(A) 小于n 次的多项式 (B) 等于n 次的多项式
(C) 大于n 次的多项式 (D) 小于等于n 次的多项式 四、计算题
1．已知函数解答下列问题的函数表如下,)(x f y =
(1)(2)分别写出四次牛顿向前插值公式和牛顿向后插值公式； (3)用三次插值多项式求)32.0()04.0(f f 和的近似值。
2．已知0.20)1.3(,5.18)4.2(,4.17)6.1(,8.14)3.1(====f f f f ，按最小二乘原理求一次多项式拟合上述数据。
3．求超定方程组 ???
??=+=-=+11
2354223212121x x x x x x 的最小二乘解。
4．已知观察值
4
3
2
1
21012y y y y y y x i
i
--
利用的二次拟合多项式)(x f )0(),(2
f x p '求的近似值。
5．用形如b x a y +=ln 的函数拟合下列数据
一、填空题
1.已知(1) 1.1f =，(2) 1.2f =，(3) 1.5f =，则三点式高斯求积公式为31
()d f x x ≈?( )，用抛物线求积公式求得3
1
()d f x x ≈?（ ）。 2.已知()30=f ，()45.0=f ，()31=f ，则用三点式可求得
计算方法——第二章——课后习题答案刘师少
2.1 用二分法求方程013=--x x 在[1, 2]的近似根，要求误差不超过3102 1-?至少要二分多少？ 解：给定误差限ε＝0.5×10－3，使用二分法时，误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(2 11 a b k 即可，亦即 96678.912lg 10lg 35.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n ＝10. 2.3 证明方程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，1]单调减少，所以f (x ) 在区间 [0，1]内有唯一实根. 给定误差限ε＝0.5×10－4，使用二分法时，误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211 a b k 即可，亦即 7287.1312 lg 10lg 45.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n ＝14. 2.4 方程0123=--x x 在x =1.5附近有根，把方程写成四种不同的等价形式，并建立相应的迭代公式： （1）211x x +=，迭代公式2111k k x x +=+ （2）231x x +=，迭代公式3211k k x x +=+ （3）112-=x x ，迭代公式111-=+k k x x （4）13-=x x ，迭代公式131-=+k k x x 试分析每种迭代公式的收敛性，并选取一种收敛迭代公式求出具有四位有效数字的近似根。 解：（1）令211)(x x f + =，则3 2)(x x f -='，由于 159.05.112)(33<≈≤='x x f ，因而迭代收敛。 （2）令321)(x x f +=，则322)1(3 2)(-+='x x x f ，由于
简便运算的练习试题和答案
乘法交换律：a×b＝b×a 乘法结合律：（a×b）×c＝a×（b×c） 38×25×4 42×125×8 25×17×4 （25×125）×（8×4） 49×4×5 38×125×8×3 (125×25)×4 5 ×289×2 （125×12）×8 125×（12×4） 乘法交换律和结合律的变化练习 125×64 125×88 44×25 125×24 25×28 加法交换律：a＋b＝b＋a 加法结合律：（a＋b）＋c＝a＋（b＋c） 357＋288＋143 158＋395＋105 167＋289＋33 129＋235＋171＋165 378＋527＋73 169＋78＋22 58＋39＋42＋61 138＋293＋62＋107 乘法分配律：（a＋b）×c＝a×c＋b×c （80＋4）×25 （20＋4）×25 （125＋17）×8 25×（40＋4）15×（20＋3）
乘法分配律正用的变化练习： 36×3 25×41 39×101 125×88 201×24 乘法分配律反用的练习： 34×72＋34×28 35×37＋65×37 85×82＋85×18 25×97＋25×3 76×25＋25×24 乘法分配律反用的变化练习： 38×29＋38 75×299＋75 64×199＋64 35×68＋68＋68×64 其他的一些简便运算。 800÷25 6000÷125 3600÷8÷5 58×101－58 74×99
姓名： （1）125×15×8×4 （2）25×24 （3）125×16 （4）75×16 （5）125×25×32 （6）25×5×64×125 （7）125×64+125×36 （8）64×45+64×71-64×16 （9）21×73+26×21+21 姓名：（1）（720+96）÷24 （2）（4500-90）÷45 （3）6342÷21 （4）8811÷89 （5）73÷36+105÷36+146÷36 （6）（10000-1000-100-10）÷10 （7）238×36÷119×5 （8）138×27÷69×50 （9）624×48÷312÷8 （10）406×312÷104÷203
计算方法引论课后答案.
第一章 误差 1. 试举例,说明什么是模型误差,什么是方法误差. 解: 例如,把地球近似看为一个标准球体,利用公式2 4A r π=计算其表面积,这个近似看为球体的过程产生 的误差即为模型误差. 在计算过程中,要用到π,我们利用无穷乘积公式计算π的值: 12 222...q q π=? ?? 其中 11 2,3,... n q q n +?=?? ==?? 我们取前9项的乘积作为π的近似值,得 3.141587725...π≈ 这个去掉π的无穷乘积公式中第9项后的部分产生的误差就是方法误差,也成为截断误差. 2. 按照四舍五入的原则,将下列各数舍成五位有效数字: 816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236 3. 下列各数是按照四舍五入原则得到的近似数,它们各有几位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位 三位 六位 四位 4. 若1/4用0.25表示,问有多少位有效数字? 解: 两位 5. 若 1.1062,0.947a b ==,是经过舍入后得到的近似值,问:,a b a b +?各有几位有效数字? 解: 已知4311 d 10,d 1022 a b --
《数值计算方法》试题集及答案
《数值计算方法》复习试题 一、填空题： 1、????? ?????----=410141014A ，则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：， 3、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ， 拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； （ 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公式为
( )] ,(),([2111+++++=n n n n n n y x f y x f h y y )； 10、已知f (1)＝2，f (2)＝3，f (4)＝，则二次Newton 插值多项式中x 2系数为( )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精 度为( 5 )； 12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 13、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式 19992001-改写为 199920012 + 。 14、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间 为 ，1 ,进行两步后根的所在区间为 ， 。 15、 、 16、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 ，用辛卜 生公式计算求得的近似值为 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。 17、 求解方程组?? ?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为 ?????-=-=+++20/3/)51()1(1)1(2)(2)1(1 k k k k x x x x ，该迭 代格式的迭代矩阵的谱半径)(M ρ= 121 。 18、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿 插值多项式为 )1(716)(2-+=x x x x N 。 19、 求积公式 ?∑=≈b a k n k k x f A x x f )(d )(0 的代数精度以( 高斯型 )求积公式为最高，具 有( 12+n )次代数精度。
小学四年级简便运算的练习题和答案
运算定律练习题 （1）乘法交换律：a×b＝b×a 乘法结合律：（a×b）×c＝a×（b×c） 38×25×4 42×125×8 25×17×4 （25×125）×（8×4） 49×4×5 38×125×8×3 (125×25)×4 — 5 ×289×2 （125×12）×8 125×（12×4） (2) 乘法交换律和结合律的变化练习 | 125×64 125×88 44×25 125×24 25×28 （3）加法交换律：a＋b＝b＋a 加法结合律：（a＋b）＋c＝a＋（b＋c） 357＋288＋143 158＋395＋105 167＋289＋33 129＋235＋171＋165 ~ 378＋527＋73 169＋78＋22 58＋39＋42＋61 138＋293＋62＋107
（4）乘法分配律：（a＋b）×c＝a×c＋b×c 正用练习 （80＋4）×25 （20＋4）×25 （125＋17）×8 25×（40＋4）15×（20＋3） （5）乘法分配律正用的变化练习： 36×3 25×41 39×101 125×88 201×24 ( （6）乘法分配律反用的练习： 34×72＋34×28 35×37＋65×37 85×82＋85×18 25×97＋25×3 76×25＋25×24 ~ （7）乘法分配律反用的变化练习： 38×29＋38 75×299＋75 64×199＋64 35×68＋68＋68×64 ； ☆思考题：（8）其他的一些简便运算。 800÷25 6000÷125 3600÷8÷5 58×101－58 74×99
【思路导航】在除法里，被除数和除数同时乘或除以一个相同的数，商不变。 325÷25 =（325×4）÷（25×4） =1300÷100 =13 【练一练1】 （1）450÷25 （2）525÷25 （3）3500÷125 / （4）10000÷625 （5）49500÷900 （6）9000÷225 ! 【经典例题二】计算25×125×4×8 【思路导航】如果先把25与4相乘，可以得到100，同时把125与8相乘，可以得到1000；再把100和1000相乘就可以了。运用了乘法交换律和结合律。 25×125×4×8 =（25×4）×（125×8） =100×1000 =100000【练一练2】 （1）125×15×8×4 （2）25×24 （3）125×16 （4）75×16 （5）125×25×32 （6）25×5×64×125 (
计算方法上机作业
计算方法上机报告 姓名： 学号： 班级： 上课班级：
说明： 本次上机实验使用的编程语言是Matlab 语言，编译环境为MATLAB 7.11.0，运行平台为Windows 7。 1. 对以下和式计算： ∑ ∞ ? ?? ??+-+-+-+=0681581482184161n n n n S n ，要求： ① 若只需保留11个有效数字，该如何进行计算； ② 若要保留30个有效数字，则又将如何进行计算； （1） 算法思想 1、根据精度要求估计所加的项数，可以使用后验误差估计，通项为： 1421114 16818485861681 n n n a n n n n n ε??= ---<< ?+++++??； 2、为了保证计算结果的准确性，写程序时，从后向前计算； 3、使用Matlab 时，可以使用以下函数控制位数： digits(位数)或vpa(变量，精度为数) （2）算法结构 1. ;0=s ?? ? ??+-+-+-+= 681581482184161n n n n t n ； 2. for 0,1,2,,n i =??? if 10m t -≤ end; 3. for ,1,2,,0n i i i =--??? ;s s t =+
（3）Matlab源程序 clear; %清除工作空间变量 clc; %清除命令窗口命令 m=input('请输入有效数字的位数m='); %输入有效数字的位数 s=0; for n=0:50 t=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); if t<=10^(-m) %判断通项与精度的关系break; end end; fprintf('需要将n值加到n=%d\n',n-1); %需要将n值加到的数值 for i=n-1:-1:0 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s=s+t; %求和运算 end s=vpa(s,m) %控制s的精度 （4）结果与分析 当保留11位有效数字时，需要将n值加到n=7， s =3.1415926536； 当保留30位有效数字时，需要将n值加到n=22, s =3.14159265358979323846264338328。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。
计算方法课后题答案之习题二
习题二 1. 证明方程043 =-+x x 在区间[1，2]内有一个根。如果用二分法求它具有5位有效数字的根，需要 二分多少次。 证明: （1） 不妨令 4)(3-+=x x x f ，求得： 02)1(<-=f 06)2(>=f 又因为4)(3-+=x x x f 在区间[1，2]内是连续的，所以在区间[1，2]内有至少一个根。 又因为 13)(2'+=x x f 在区间[1，2]内013)(2'>+=x x f ，所以4)(3-+=x x x f 单调。 得证，043 =-+x x 在区间[1，2]内仅有一个根。 （2）具有5位有效数字的根，说明根可以表示成 5 4321.a a a a a ，所以绝对误差限应该是 5a 位上的 一半，即： 4105.0-?=ε。由公式： ε≤-+1 2 k a b 可得到, 14=k 迭代次数为151=+k 次。 ---------------------------------------------------------------------------------------------------------------------- 2. 用二分法求方程 0)2 (sin )(2=-=x x x f 在区间[1.5，2]内的近似根（精确到10-3）。 解：043499.05625.099749.0)25.1(5.1sin )5.1(2 >=-=-=f 009070.0190930.0)22(2sin )2(2 <-=-=-=f 所以0)2 (sin )(2 =-=x x x f 在区间[1.5，2]内有根，又 x cos )('-=x x f 在区间[1.5，2]内 0x cos )('<-=x x f 所以 0)2 (sin )(2=-=x x x f 在区间[1.5，2]内有根，且唯一。符合二分条件，可以用二分法，二分的 次数为：
六年级数学简便计算专项练习题(附答案+计算方法汇总)
六年级数学简便计算专项练习题（附答案＋计算方法汇总） 小学阶段（高年级）的简便运算，在一定程度上突破了算式原来的运算顺序，根据运算定律、性质重组运算顺序。如果学生没真正理解运算定律、性质，他只能照葫芦画瓢。在实际解题的过程当中，学生的思路不清晰，常出现这样或那样的错误。因此，培养学生思维的灵活性就显得尤为重要。 下面，为大家整理了8种简便运算的方法，希望同学们在理解的基础上灵活运用，不提倡死记硬背哟！ 1.提取公因式 这个方法实际上是运用了乘法分配律，将相同因数提取出来，考试中往往剩下的项相加减，会出现一个整数。 注意相同因数的提取。 例如： 0.92×1.41＋0.92×8.59 =0.92×（1.41+8.59） 2.借来借去法 看到名字，就知道这个方法的含义。用此方法时，需要注意观察，发现规律。还要注意还哦,有借有还，再借不难。 考试中，看到有类似998、999或者1.98等接近一个非常好计算的整数的时候，往往使用借来借去法。 例如： 9999+999+99+9 =9999+1+999+1+99+1+9+1－4 3.拆分法
顾名思义，拆分法就是为了方便计算把一个数拆成几个数。这需要掌握一些“好朋友”，如：2和5，4和5，2和2.5，4和2.5，8和1.25等。分拆还要注意不要改变数的大小哦。 例如： 3.2×12.5×25 =8×0.4×12.5×25 =8×12.5×0.4×25 4.加法结合律 注意对加法结合律 (a＋b)＋c=a＋(b＋c) 的运用，通过改变加数的位置来获得更简便的运算。 例如： 5.76＋13.67＋4.24＋ 6.33 =（5.76＋4.24）＋(13.67＋6.33) 5.拆分法和乘法分配律结合 这种方法要灵活掌握拆分法和乘法分配律，在考卷上看到99、101、9.8等接近一个整数的时候，要首先考虑拆分。 例如： 34×9.9 = 34×(10－0.1) 案例再现：57×101=？ 6.利用基准数 在一系列数种找出一个比较折中的数字来代表这一系列的数字，当然要记得这个数字的选取不能偏离这一系列数字太远。 例如： 2072+2052+2062+2042+2083
《徐翠微计算方法引论》
第二章 插值法 知识点：拉格朗日插值法，牛顿插值法，余项，分段插值。 实际问题中，时常不能给出f （x ）的解析表达式或f （x ）解析表达式过于复杂而难于计算，能采集的只是一些f （x ）的离散点值{xi,f(xi)}(i=0,1,2,…n )。因之，考虑近似方法成为自然之选。 定义：设f （x ）为定义在区间[a ，b]上的函数，x0,x1,…,xn 为[a ，b]上的互异点，yi=f （xi ）。若存在一个简单函数?（x ），满足 （插值条件）?（xi ）=f （xi ），i=0，1，…，n 。 则称 ?（x ）为f （x ）插值函数，f （x ）为被插函数,点x0,x1,…,xn 为插值节点,点{xi,f(xi)},i=0,1,2,…n 为插值点。 于是计算f （x ）的问题就转换为计算 ?（x ）。 构造插值函数需要解决：插值函数是否存在唯一；插值函数如何构造(L 插值)；插值函数与被插函数的误差估计和收敛性。 对插值函数 ?（x ）类型有多种不同的选择，代数多项式常被选作插值函数。 P23(2.18)和(2.19)指出，存在唯一的满足插值条件的n 次插值多项式p n (x )。但是需要计算范德蒙行列式，构造插值多项式工作量过大，简单表达式不易得到，实际中不采用这类方法。 插值法是一种古老的数学方法，拉格朗日（Lagrange ）、牛顿（Newton ）等分别给出了不同的解决方法。 拉格朗日插值 拉格朗日（Lagrange ）插值的基本思想：把插值多项式p n (x )的构造问题转化为n+1个插值基函数l i (x)(i=0,1,…,n)的构造。 (1)线性插值 ①构造插值函数 已知函数y =f (x )的两个插值点(x 0,y 0)，(x 1,y 1)，构造多项式y =p 1(x )，使p 1(x 0)=y 0,p 1(x 1)=y 1。 p n (x )≈f (x )
数值计算方法大作业
目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20)
4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28)
数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码： Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求f(x)=e2x-6x=0在x=0.5附近的根（ep=10-10）
1.2牛顿法 程序代码： Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求56的值。（ep=10-10）
用简便方法计算下面各题
用简便方法计算下面各题 4.8×0.25 2.4×12.5 1.25×1.6×2.5 4.76×99+4.76 58.5×101-58.5 18.7×99+18.7 2.85×99 4.23×101 5.8×102 5.4×10.1 6.8×9.9 2.5×10.2 12.5×（100+8）9.4×10.1 93.7×0.32+93.7×0.68 2.52×101 1.25×0.7+1.25×1.2+12.5 3.6×2.5 7.2×0.2＋2.4×1.4 12.7×9.9＋1.2710.7×16.1-151×1.07
1、学校图书室长9.7 m，宽5.3 m，用边长0.9 m的正方形瓷砖铺地，70块够吗？（不考虑损耗。） 2、某公司出租车的收费标准如下：收费标准4 km及以内10元，超出4 km （不足1 km按1 km计算）每千米1.2元，某乘客要乘出租车去30 km处的某地，应付车费多少元？ 3小强家的固定电话收费标准如下：前3分钟收费0.4元，超过3分钟每分钟收费0.12元（不足1分钟按1分钟计算）。小强给爷爷和奶奶打电话用时8分钟52秒，他这一次通话的费用是多少？ 4、某市自来水公司供水收费标准如下：每月用水在12吨及以内，每吨收费2.65元；超出12吨部分，每吨3.8元。王琼家八月份用水18吨，付给自来水公司收费人员100元，应找回多少钱？ 5、刘强从家骑车到学校要用0.4小时，刘强的家离学校有多远？如果他改为步行，每小时走4.8km，0.9小时能到学校吗？（骑车：12千米/时） 6、我市某出租车公司租车计费方法如下：乘车路程不超过4km，收费8.5元（起步价）；超过部分按每千米1.5元加收费（不足1km，按1km计算）。爸爸和小亮乘车回家的路程为14.1km，付给出租车司机100元，应找回多少元？
《数值计算方法》试题集及答案
《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ,则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案: ?? ????????--??????????--=1556141501 4115401411A 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 ,拉 格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式就是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 7、计算方法主要研究( 截断 )误差与( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 10、已知f (1)＝2,f (2)＝3,f (4)＝5、9,则二次Newton 插值多项式中x 2系数为( 0、15 ); 11、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式
数值计算方法习题答案(绪论,习题1,习题2)
引论试题（11页） 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 ）(2 2 11222k k k k k k k k x a x a x 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 ≥+??? ? ??-=???? ??+=+2 12121 6 证明： 若k x 有n 位有效数字，则n k x -?≤ -1102 1 8， 而() k k k k k x x x x x 28882182 1-=-???? ??+=-+ n n k k x x 21221102 1 5.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.0102 21 111=??≤--x x e 71.22=x ，有两位有效数字，相对误差限为 025.0102 21 122=??≤--x x e 3 2.718x =，有两位有效数字，其相对误差限为： 00025.0102 21 333=??≤--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 位有效数字的近似数都成立的正确结论，故他对误差限的估计偏大，但计算略简单些；而第二种方法给出较好的误差限估计，但计算稍复杂。 （2）采用第二种方法时，分子为绝对误差限，不是单纯的对真实值与近似值差值的四舍五入，绝对误差限大于或等于真实值与近似值的差。 11. 解: ......142857.3722≈，.......1415929.3113 255≈ 21021 722-?≤-∴ π，具有3位有效数字 6102 1 113255-?≤-π，具有7位有效数字
四年级下册简便方法计算练习题
四年级下册简便方法计算练习题126×6×8 600÷25÷4 55×36＋64×55 755－122－78 600÷25 （8＋80）×125 125×18 234×80×5 781－499 125×38＋125×30 25×32 4004×25 25×16－25×10 25×16×125 （125＋16）×8 79×99＋79 781×101－781 79×16＋79×78＋79×6 25×101
789×99 800÷125 1736＋403 2000÷125 65＋93×65＋6×65 9999＋999＋99＋9 158+262+138 375+219+381+225 5001－247－1021－232 （181+2564）+2719 378+44+114+242+222 276+228+353+219 (375+1034)+(966+125) (2130+783+270)+1017 99+999+9999+99999 7755－(2187+755) 2214+638+286 3065－738－1065 899+344
2370+1995 3999+498 1883－398 12×25 75×24 138×25×4 (13×125)×(3×8) (12+24+80)×50 704×25 25×32×125 32×(25+125) 88×125 102×76 58×98 178×101－178 84×36+64×84 75×99+2×75 83×102－83×2 98×199 123×18－123×3+85×123 50×(34×4)×3 25×（24+16） 178×99+178 79×42+79+79×57 7300÷25÷4 8100÷4÷75 16800÷120 30100÷2100 32000÷400 49700÷700
数值计算方法试题集和答案
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：， 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ， 拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝，则二次Newton 插值多项式中x 2系数为( )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精 度为( 5 )； 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式 19992001-改写为 199920012 + 。
计算方法习题答案
计算方法第3版习题答案 习题1解答 1.1 解：直接根据定义得 *411()102x δ-≤?*411()102r x δ-≤?*3*12211 ()10,()1026 r x x δδ--≤?≤?*2*5331()10,()102r x x δδ--≤?≤ 1.2 解：取4位有效数字 1.3解：433 5124124124 ()()() 101010() 1.810257.563 r a a a a a a a a a δδδδ----++++++≤≤=?++? 123()r a a a δ≤ 123132231123 ()()() a a a a a a a a a a a a δδδ++0.016= 1.4 解：由于'1(),()n n f x x f x nx -==，故***1*(())()()()n n n f x x x n x x x δ-=-≈- 故** * ***(()) (())()0.02()r r n f x x x f x n n x n x x δδδ-= ≈== 1.5 解： 设长、宽和高分别为 ***50,20,10l l h h εεωωεεεε=±=±=±=±=±=± 2()l lh h ωωA =++，*************()2[()()()()()()]l l l h h l h h εδωωδδδωδδωA =+++++ ***4[]320l h εωε=++= 令3201ε<，解得0.0031ε≤， 1.6 解：设边长为x 时，其面积为S ，则有2()S f x x ==，故 '()()()2()S f x x x x δδδ≈= 现100,()1x S δ=≤，从而得() 1 ()0.00522100 S x x δδ≈ ≤ =? 1.7 解：因S ld =，故 S d l ?=?，S l d ?=?，*****()()()()()S S S l d l d δδδ??≈+?? * 2 ()(3.12 4.32)0.010.0744S m δ=+?=, *** ** * () () 0.0744 ()0.55%13.4784 r S S S l d S δδδ= = = ≈ 1.8 解：（1）4.472 （2）4.47 1.9 解：(1) （B ）避免相近数相减 (2)（C ）避免小除数和相近数相减 (3)（A ）避免相近数相减 (3)（C ）避免小除数和相近数相减，且节省对数运算 1.10 解 （1）357sin ...3!5!7!x x x x x =-+-+ 故有357 sin ..3!5!7! x x x x x -=-+-， （2） 1 (1)(1)1lnxdx ln ln ln N+N =N N +-N N +N +-? 1 (1)1ln ln N +=N +N +-N 1.11 解：0.00548。 1.12解：21 16 27 3102 ()()() -? 1.13解：0.000021
计算方法引论课后答案
第一章 误差 1. 试举例,说明什么是模型误差,什么是方法误差. 解: 例如,把地球近似看为一个标准球体,利用公式2 4A r π=计算其表面积,这个近似看为球体的过程产生的误差即为模型误差. 在计算过程中,要用到π,我们利用无穷乘积公式计算π的值: 其中 我们取前9项的乘积作为π的近似值,得 这个去掉π的无穷乘积公式中第9项后的部分产生的误差就是方法误差,也成为截断误差. 2. 按照四舍五入的原则,将下列各数舍成五位有效数字: 816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236 3. 下列各数是按照四舍五入原则得到的近似数,它们各有几位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位 三位 六位 四位 4. 若1/4用0.25表示,问有多少位有效数字? 解: 两位 5. 若 1.1062,0.947a b ==,是经过舍入后得到的近似值,问:,a b a b +?各有几位有效数字? 解: 已知4311 d 10,d 1022 a b --< ?
数值分析计算方法试题集及答案
数值分析复习试题 第一章 绪论 一. 填空题 1.* x 为精确值 x 的近似值；() **x f y =为一元函数 ()x f y =1的近似值； ()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-： *** r x x e x -= ()()()*'1**y f x x εε≈? ()() () ()'***1**r r x f x y x f x εε≈ ? ()()()() ()* *,**,*2**f x y f x y y x y x y εεε??≈?+??? ()()()()() ** * *,***,**222r f x y e x f x y e y y x y y y ε??≈ ?+??? 2、 计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 舍入误 差 。 3、 分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有 6 位和 7 位；又取 1.73≈-21 1.73 10 2 ≤?。 4、 设121.216, 3.654x x ==均具有3位有效数字，则12x x 的相对误差限为 0.0055 。 5、 设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为 0.01 。 6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得 到,则相对误差限为 0.0000204 . 7、 递推公式,??? ? ?0n n-1y =y =10y -1,n =1,2, 如果取0 1.41y ≈作计算,则计算到10y 时,误 差为 81 10 2 ?;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3* =π，则近似值141.3*1=π和1415.3*2=π分别有 3
计算方法-刘师少版课后习题答案
1.1 设3.14, 3.1415, 3.1416分别作为π的近似值时所具有的有效数字位数 解 近似值x =3.14＝0.314×101,即m =1，它的绝对误差是 －0.001 592 6…，有 31105.06592001.0-*?≤=- x x . 即n =3，故x =3.14有3位有效数字. x =3.14准确到小数点后第2位. 又近似值x =3.1416，它的绝对误差是0.0000074…，有 5-1*10?50≤00000740=-.. x x 即m =1,n ＝5，x =3.1416有5位有效数字. 而近似值x =3.1415，它的绝对误差是0.0000926…，有 4-1*10?50≤00009260=-.. x x 即m =1,n ＝4，x =3.1415有4位有效数字. 这就是说某数有s 位数，若末位数字是四舍五入得到的，那么该数有s 位有效数字 1.2 指出下列各数具有几位有效数字，及其绝对误差限和相对误差限： 2.0004 －0.00200 9000 9000.00 解 （1）∵ 2.0004＝0.20004×101, m=1 绝对误差限：4105.0000049.020004.0-*?≤≤-=-x x x m -n =-4,m =1则n =5,故x =2.0004有5位有效数字 1x =2，相对误差限000025.010******** 1)1(1 =??=??=---n r x ε （2）∵ －0.00200= -0.2×10-2, m =-2 5105.00000049.0)00200.0(-*?≤≤--=-x x x m -n =-5, m =-2则n =3,故x =-0.00200有3位有效数字 1x =2，相对误差限3 110221 -??=r ε=0.0025 (3) ∵ 9000=0.9000×104, m =4， 0105.049.09000?<≤-=-*x x x m -n =0, m =4则n =4,故x =9000有4位有效数字 4 110921-??=r ε＝0.000056 (4) ∵9000.00=0.900000×104, m =4， 2105.00049.000.9000-*?<≤-=-x x x m -n =-2, m =4则n =6,故x =9000.00有6位有效数字 相对误差限为6 110921-??=r ε＝0.000 00056 由(3)与(4)可以看到小数点之后的0，不是可有可无的，它是有实际意义的. 1.3 ln2=0.69314718…，精确到310-的近似值是多少？ 解 精确到310-＝0.001，即绝对误差限是ε＝0.0005, 故至少要保留小数点后三位才可以.ln2≈0.693 2.1 用二分法求方程013=--x x 在[1, 2]的近似根，要求误差不超过 31021-?至少要二分多少？ 解：给定误差限ε＝0.5×10－3，使用二分法时，误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211a b k 即可，亦即 96678.912lg 10lg 35.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n ＝10. 2.3 证明方程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，1]单调减少，所以f (x ) 在区间 [0，1]内有唯一实根. 给定误差限ε＝0.5×10－4，使用二分法时，误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211a b k 即可，亦即 7287.1312lg 10lg 45.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n ＝14. 2.4 方程0123=--x x 在x =1.5附近有根，把方程写成四种不同的等价形式，并建立相应的迭代公式： （1）211x x +=，迭代公式2111k k x x +=+ （2）231x x +=，迭代公式3211k k x x +=+ （3）112-=x x ，迭代公式111-=+k k x x （4）13-=x x ，迭代公式131-=+k k x x 试分析每种迭代公式的收敛性，并选取一种收敛迭代公式求出具有四位有效数字的近似根。 解：（1）令211)(x x f +=，则32)(x x f -='，由于
