数值计算习题(1)
数值计算第一二章答案
第一章数值计算中的误差习题一1.1 下列各近似数的绝对误差限是最末位的半个单位，试指出它们各有几位有效数字。
1x =-3。
105 , 2x =0.001, 3x =0。
100, 4x =253。
40， 5x =5000， 6x =5⨯310.答案：4，1，3，6，4，1。
1。
2 设100〉*x >10,x 是*x 的有五位有效数字的的近似数,求x 的绝对误差限。
答案:当10<x 〈100时，因为有5位有效数字,所以绝对误差限为0。
005。
1。
3 求下列各近似数的相对误差限和有效数字位数： 1） 123x x x ++,2） 124x x x 3） 24x x 答案：()10.0005e x ≤()20.0005e x ≤()30.0005e x ≤ ()40.005e x ≤ ()50.5e x ≤ ()60.5e x ≤1）()()()()123123e x x x e x e x e x ++=++≤()()()123e x e x e x ++3221.5100.15100.510---≤⨯=⨯≤⨯2123()0.1510x x x ε-++=⨯123123123()()0.0004993...0.0004994r x x x e x x x x x x ε++++==≤++123x x x ++=-3。
004 精确到小数点后两位,所以有三位有效数字。
2）()()()()()()12424112424114224()e x x x x x e x x e x x x x e x x x e x x e x =+=++ =()()()241142124)x x e x x x e x x x e x ++()()()241142124x x e x x x e x x x e x ≤++ =660.5100.31050.0005 3.1050.510--⨯+⨯+⨯⨯ 所以43124() 1.71275100.510x x x ε--=⨯≤⨯124x x x =43.105100.0003105--⨯=-41241244124() 1.7127510()0.5515...3.10510r x x x e x x x x x x ε--⨯===⨯3）()()2222424244444()()1x x e x x e x e e x e x x x x x x ⎛⎫≈-≤+⎪⎝⎭325105420.5100.5100.197316100.77868100.1997100.510253.40253.40------⨯⨯=+=⨯+⨯≈⨯<⨯ 又由24x x 50.3946310-≈⨯知有0位有效数字 ∴522440.1997100.5r x e x x x -⎛⎫⨯≤≈ ⎪⎝⎭1。
完整word版,《数值计算方法》试题集及答案(1-6) 2..
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：2.367，0.252、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*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 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()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+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。
数值计算方法习题答案
《数值计算方法》课后题答案详解吉 林 大 学第一章 习 题 答 案1. 已知(1)2,(1)1,(2)1f f f −===，求()f x 的Lagrange 插值多项式。
解：由题意知：()01201212001020211012012202121,1,2;2,1,1()()(1)(2)()()6()()(1)(2)()()2()()(1)(1)()()3(1)(2)(1)(2)()2162nj j j x x x y y y x x x x x x l x x x x x x x x x x l x x x x x x x x x x l x x x x x x x x L x y l x ==−=====−−−−==−−−−+−==−−−−−+−==−−−−+−==×+×−∴∑()2(1)(1)131386x x x x +−+×=−+2. 取节点01210,1,,2x x x ===对x y e −=建立Lagrange 型二次插值函数，并估计差。
解11201201210,1,;1,,2x x x y y e y e −−======1)由题意知：则根据二次Lagrange插值公式得：02011201201021012202110.510.520.51()()()()()()()()()()()()()2(1)(0.5)2(0.5)4(1)(224)(43)1x x x x x x x x x x x x L x y y y x x x x x x x x x x x x x x x x e x x e e e x e e x −−−−−−−−−−−−=++−−−−−−=−−+−−−=+−+−−+22)Lagrange 根据余项定理，其误差为(3)2210122()1|()||()||(1)(0.5)|3!61max |(1)(0.5)|,(0,1)6()(1)(0.5),()330.5030.2113()61()0.2113(0.21131)(0.21130.5)0.008026x f R x x e x x x x x x t x x x x t x x x x t x R x ξξωξ−+≤≤==−−≤−−∈′=−−=−+=−==≤××−×−=∴取 并令 可知当时，有极大值3. 已知函数y =在4, 6.25,9x x x ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。
《数值计算方法》试题集及答案(1-6)#优选.
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：2.367，0.252、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*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 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()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+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。
数值计算方法复习题
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） 绝对误差 ， 相对误差 ， 有效数字 ， 截断误差 ， 舍入误差 。
数值习题
1.1 数值计算方法概述 选择题1. 易 (1分)1.41300的近似值，有( )位有效数字。
A 、3；B 、4；C 、5；D 、62. 易 (1分)0.026900x *=-作为x 的近似值，它的有效数字位数为( ) 。
A 、7；B 、3；C 、不能确定D 、5.3. 易 (1分)1.41300的近似值，有( )位有效数字。
A 、3；B 、4；C 、5；D 、6。
4. 中下 (1分)下列说法错误的是（ ）。
A 、如果一个近似数的每一位都是有效数字，则称该近似数为有效数B 、凡是经“四舍五入”得到的近似数都是有效数C 、数值方法的稳定性是指初始数据的扰动对计算结果的影响D 、病态问题是由数学问题本身的性质决定的，与数值方法有关5. 中下 (1分)下列说法中不属于数值方法设计中的可靠性分析的是（ ）。
A 、方法收敛性；B 、方法的稳定性；C 、方法的计算量；D 、方法的误差估计6. 难 (1分) 1.732≈计算41)x =，下列方法中哪种最好？（ ）A 、28-B 、2(4- CD 7. 中 (1分)近似数*0.231x =关于真值0.229x =有（ ）位有效数字。
A 、1B 、2C 、3D 、48. 易 (1分)-324.7500是舍入得到的近似值，它有( )位有效数字。
A 、5B 、6C 、7D 、89. 中下 (1分)设π的近似数*π有3位有效数字，则其相对误差限为（ ）。
A 、21102-⨯ B 、31102-⨯ C 、41102-⨯ D 、51102-⨯ 10. 中上 (1分)在近似计算中，要注意以下原则：（1）计算速度快 （2）避免大数“吃掉”小数， （3）防止溢出 （4）减少计算次数 列主元消元法解方程组Ax b =是（ ）.A 、(1)和(2)B 、(2)和(3)C 、(3)和(4)D 、(4)和(1)11. 难 （1分）已知94763.≈有三位有效数字，则方程01162=+-x x 的具有三位有效数字的较小根为（ ）。
数值计算方法答案
1数值计算方法习题一（2）习题二（6）习题三（15）习题四（29）习题五（37）习题六（62）习题七（70）2009．9，92习题一1．设x>0相对误差为2%，4x的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x xf x f xδδ∆=≈得（1）()f x=11()()*2%1%22x xδδδ≈===；（2）4()f x x=时444()()'()4()4*2%8%xx x x xxδδδ≈===2．设下面各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出他们各有几位有效数字。
（1）12.1x =；（2）12.10x =；（3）12.100x =。
解：由教材9P关于1212.m nx a a a bb b=±型数的有效数字的结论，易得上面三个数的有效数字位数分别为：3，4，53．用十进制四位浮点数计算（1）31.97+2.456+0.1352；（2）31.97+（2.456+0.1352）哪个较精确？解：（1）31.97+2.456+0.1352≈21((0.3197100.245610)0.1352)fl fl⨯+⨯+=2(0.3443100.1352)fl⨯+=0.3457210⨯（2）31.97+（2.456+0.1352）21(0.319710(0.245610))fl fl≈⨯+⨯= 21(0.3197100.259110)fl⨯+⨯=0.3456210⨯易见31.97+2.456+0.1352=0.345612210⨯，故（2）的计算结果较精确。
4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？3解：设该正方形的边长为x，面积为2()f x x=，由(())(())'()()()()f x xf x f x xf x f xδδ∆=≈解得(())()()'()f x f xxxf xδδ≈=2(())(())22f x x f xx xδδ==0.5%5．下面计算y的公式哪个算得准确些？为什么？（1）已知1x<<，（A）11121xyx x-=-++，（B）22(12)(1)xyx x=++；（2）已知1x>>，（A）y=，（B）y=；（3）已知1x<<，（A）22sin xyx=，（B）1cos2xyx-=；（4）（A）9y=（B）y=解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
数值计算方法复习题
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） 绝对误差 ， 相对误差 ， 有效数字 ， 截断误差 ， 舍入误差 。
数值计算方法习题.doc
第一章 绪论1．把下列各数按四舍五入规则舍入为有3位小数的近似数，并写出近似数的绝对误差和相对误差，指出近似数有几位有效数字： 93.1822 4.32250 15.9477 17.3675 2．按四舍五入原则，将下列数舍成五位有效数字：816.9567 6.000015 17.32250 1.235651 93.18213 0.015236233．设 **,671.3,6716.3x x x 则==有几位有效数字？ 4．若1/4用0.25来表示,问有多少位有效数字? 5．若 1.1062,0.947a b ==是经过舍入后得到的近似值,问:,a b a b +⨯各有几位有效数字?6．设120.9863,0.0062y y ==是经过舍入后作为12,x x 的近似值, 求11y 和21y 的计算值与真值的相对误差限及12y y 和得到真值的相对误差限. 7．设0,x x >的相对误差为δ,求ln x 的绝对误差.8．正方形的边长约为100cm ,应该怎样测量,才能使其面积的误差不超过12cm . 9．设x 的相对误差为a %,求x n 的相对误差.10．计算球的体积,为了使相对误差限为1%,问度量半径R 时允许的相对误差限如何?11．5631.2*=x 是经四舍五入得到的近似值，则其相对误差≤*r e __________ 12． 设 0000073.0,1416.3,1415926.3**=-==x x x x 则称_________误差13．设⎰+=1061dx xx I nn ，设计一个计算10I 的算法，并说明你的算法的合理性。
14．设028Y =，按递推公式1n n Y Y -= （1,2,n = ）， 计算到100Y27.982≈（5位有效数字），试问计算100Y 将有多大误差． 15．求方程25610x x -+=的两个根，使它至少具有4位有效数字27.982≈）．16．当N 充分大时，怎样求121d 1N N x x ++⎰？17．序列{}n y 满足递推关系101n n y y =- （1,2,n = ），若0 1.41y =≈（三位有效数字），计算到10y 时误差有多大？这个计算过程稳定吗？18．计算61)f =1.41≈，利用下列算式计算，哪一个得到的结果最好？，3(3-，99- 19．()ln(f x x =，求(30)f 的值，若开平方用6位函数表，问求对数时误差有多大？若改用另一等价公式ln(ln(x x =- 计算，求对数时误差有多大？第二章 解线性方程组的直接方法1．用高斯消去法解方程组123234011921261x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦ 2．用LU 分解，将上题系数矩阵分解为L 和U 的乘积，L 是对角线元素为1的下三角矩阵，U 是上三角矩阵。
最新数值计算课后答案1
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x -==≈⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯ 所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x --==≈-⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.15＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯ 所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈-L L 相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯ 有效数字：因为π＝3.14159265...=0.314159265 (10)22 3.1428571430.3142857143107==⨯，m=1。
数值计算方法习题答案(第二版)(绪论)
数值分析（p11页）4 试证：对任给初值x 0,0)a >的牛顿迭代公式112(),0,1,2,......k ak k x x x k +=+= 恒成立下列关系式：2112(1)(,0,1,2,....(2)1,2,......kk k x k x x k x k +-=-=≥=证明：（1）(21122k k k k k kx 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-=-⎪⎪⎭⎫⎝⎛+=-+ nnk 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 位有效数字的近似数都成立的正确结论，故他对误差限的估计偏大，但计算略简单些；而第二种方法给出较好的误差限估计，但计算稍复杂。
数值计算方法》试题集和答案(1_6)2
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：，2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*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 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝，则二次Newton 插值多项式中x 2系数为( )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()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+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 ，1 ,进行两步后根的所在区间为 ， 。
《数值计算方法》试题集和答案(1_6)2
《计算方法》期中复习试题一、填空题:1、已知/（1）= 1°，/⑵=12 /⑶= 1.3,则用辛普生（辛卜生）公式计算求2、/（1）= 7 /⑵=2, /（3） = 11则过这三点的二次插值多项式中疋的系数为 _____ ,拉格朗日插值多项式为 ________________________ o3、近似值^=0 231关于真值A = 0.229有（2 ）位有效数字;4、设/（X ）可微,求方程x = f^的牛顿迭代格式是（）;5、 对/U ）= X 3 + A + 1J 差商/[0,1,2,3] =（ 1/[0,1,2,3,4] =（）；6、 计算方法主要研究（截断）误差和（舍入）误差；7、用二分法求非线性方程f （x ）=0在区间（s,®内的根时，二分力次后的误差限b_a为（诃 ）;8、已知f （l ）=2, f （2）=3, f （4）=,则二次Newton 插值多项式中#系数为（ ）;f/Wdv打⑴血-+ 八2^）]11、两点式高斯型求积公式J 。
八^（Jo22血2厲 ），代数精度为（5 ）;“34 6 y = 10 H ------ 1 --------- ------------12、 为了使计算x-1 （x-1）- （x-1）-的乘除法次数尽量地少，应将V2ooi - 71999"改写为、/^55T+Vi^ 。
得丄"如 答案:，用三点式求得广⑴〜 ____________答案：心+1 =心答案1 一广（占）该表达式改写为 y = 10+（3 + （4 — 6小）人 t =x-l ,为了减少舍入误差，应将表达式13、用二分法求方程/（X）= Q + x_ 1 = 0在区间［0, 1］内的根，进行一步后根的所在区间为，1 ，进行两步后根的所在区间为，。
14、计算积分匚5低肚，取4位有效数字。
用梯形公式计算求得的近似值为—，用辛卜生公式计算求得的近似值为梯形公式的代数精度为1 ,辛卜生公式的代数精度为3 。
《数值计算方法》试题集及答案(1-6) 2
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：2.367，0.252、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*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 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()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+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。
数值计算(一)
[x,fval]=fminbnd('sin(x)',-pi,pi) %求sinx在[-π,π]的极小值点 x= -1.5708 fval = -1.0000
2014-6-28
天津农学院基础科学系 朱文新
22
[x,fval]=fminsearch(' (x(1)-1)^2+(x(2)x(1)^2)^2',[0,0])
级数求和、幂级数展开的Matlab命令
命令格式 symsum(f) symsum(f,x) 含义 求符号表达式f对默认变量的不定和 求符号表达式f对自变量x的不定和
symsum(f,a,b)
求符号表达式f对默认变量的a到b的有限和
taylor(f) taylor(f,n,x) taylor(f,n,x,a)
24
练习
1.
求下列极限
1 x 1 lim x 0 x
n
xn 1 lim m , ( n, m为常数) x 1 x 1
ln x lim ( ) x x 1
x 1 lim( ) x0 x 1 ln x
2014-6-28
天津农学院基础科学系 朱文新
25
2、计算下列导数 (1) (2) (3) (4)
MATLAB数值计算（一）
主讲： 朱文新
2014-6-28
天津农学院基础科学系 朱文新
1
实验内容：
1. 微积分基本运算 2.非线性方程组的解法
2014-6-28
天津农学院基础科学系 朱文新
2
一、微积分基本运算
实验目的：
学会利用matlab 求解函数的极限、导数、极值、积分与级数运算
《数值计算方法》试题集和答案(1_6)
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：2.367，0.252、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*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 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()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+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。
现代数值计算方法习题解答
E( x n ) x − x* ≈n = nE r ( x) = 0.01n . 所以 E r ( x ) = x xn
n
8 、解：
9 、证： E ( S ) = S − S * ≈ gt (t − t * ) = gtE (t ) Er (S ) = S − S * gt (t − t * ) 2 E (t ) ≈ = S t gt 2 / 2 由上述两式易知，结论.
* * 10 x 0 − 1 − 10 x0 + 1 |=1 x0 − x 0 | x1 − x1* |=| 0 | |< = 10δ
* * | x2 − x2 |=| 10 x1 − 1 − 10 x1 + 1 |=1 0 | x1 − x1* |< = 10 2 δ
类推有
* x10 − x10 | | < =1010 δ =
4
现代数值计算方法习题答案
是否大于零来判断.
3 2 3 3 2 所以系数矩阵是对称 = 2 > 0, 2 2 0 = 6 > 0, a11 ＝ 3 > 0, 2 2 3 0 12
正定的. 记系数矩阵为 A，则平方根法可按如下三步进行： 第一步 分解：A = L LT. 由公式计算出矩阵的各元素：
l11 = 3
3 2 1 3 2 所以系数矩阵是对称 = 2 > 0, 2 2 0 = 4 > 0, a11 ＝ 3 > 0, 2 2 1 0 3
正定的. 记系数矩阵为 A，则平方根法可按如下三步进行： 第一步 分解：A = L LT. 由公式计算出矩阵的各元素：
l11 = 3
l 21 =
2 3 3 6 3
l 22 =
所以数组 A 的形式为：
《数值计算》课后习题汇总
《数值计算》课后习题汇总习题⼀解答1．取3.14，3.15，227，355113作为π的近似值，求各⾃的绝对误差，相对误差和有效数字的位数。
2、⽤四舍五⼊原则写出下列各数的具有五位有效数字的近似数。
346．7854，7．000009，0．0001324580，0．6003003、下列各数都是对准确数进⾏四舍五⼊后得到的近似数，试分别指出他们的绝对误差限和相对误差限和有效数字的位数。
4.0.1％。
5、在计算机数系F(10,4,-77,77)中，对31120.14281100.31415910x x =?=-?与，试求它们的机器浮点数()(1,2)i fl x i =及其相对误差。
6、在机器数系F(10,8,L,U)中，取三个数4220.2337125810,0.3367842910,0.3367781110x y z -=?=?=-?，试按(),()x y z x y z ++++两种算法计算x y z ++的值，并将结果与精确结果⽐较。
7、某计算机的机器数系为F(10,2,L,U)，⽤浮点运算分别从左到右计算及从右到左计算10.40.30.20.040.030.020.01+++++++试⽐较所得结果。
8、对于有效数1233.105,0.001,0.100x x x =-==，估计下列算式的相对误差限21123212333,,x y x x x y x x x y x =++==9、试改变下列表达式，使其计算结果⽐较精确（其中1x 表⽰x 充分接近0，1x表⽰x 充分⼤）。
(1)1212ln ln ,x x x x -≈； (2)11,111x x x x---+；(3)1x ；(4)cot ,01x x xx-≠且。
10、⽤4位三⾓函数表，怎样算才能保证1cos 2-有较⾼的精度？11、利⽤27.982≈求⽅程25610x x -+=的两个根，使它们⾄少具有4位有效数字。
12、试给出⼀种计算积分11(0,1,2,3,...)n x n I ex e dx n -==?，近似值的稳定算法。
数值计算习题一
习题一1.表中各*x都是对准确值x进行四舍五入得到的近似值。
试分别指出解：（1）当x=±0.x1x2....x s10m作为数*x的近似值时，x1x2...x s不一定为x的有效数字。
但是用四舍五入取准确值*x的前n位作为近似值0.x1x2....x n10m，则x必有n个有效数字x1x2...x n。
因为各0.3012,30.12=0.3012×102都是对准确值x进行四舍五入得到的近似值，所以0.3012,30.12都有4位有效数字3012 而30.120=0.30120×102有5位有效数字30120。
（2）根据有效数字的定义：设数x*的近似值x=±0.x1x2....x s10m，其中x i （i=1,2,...,s）是0到9之间的任一个正整数，且x1≠0，n是正整数，m是整数，如果绝对误差e*的|e*|=|x*-x|≤½×10m-n则称x为x*的具有n位有效数字的近似值，x准确到第n位，x1x2...x n为x的有效数字。
所以，具有四位有效数字的数0.3012,30.12=0.3012×102的绝对误差限分别为|x*-0.3012|≤½×100-4=0.5×10-4，|x*-30.12|≤½×102-4=0.5×10-2。
具有五位有效数字的数30.120=0.30120×102的绝对误差限分别为|x*-30.120|≤½×102-5=0.5×10-3。
（3）根据定理：设数x*的近似值x=±0.x1x2....x n10m具有n位有效数字，则x的相对误差满足下列不等式|e r*|≤½x1×101-n所以，具有四位有效数字的数0.3012,30.12=0.3012×102的相对误差限都为|e r*|≤½x1×101-n=(1/2×3)×101-4=1/6×10-3。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第5章 MATLAB 数值计算
1.选择和填空。
（1）下列变量名中的________是合法变量。
A. char_1 , i , j
B. x*y , a.1
C. x\y , a1234
D. end , 1bcx
（2）已知x 为1个向量，计算ln(x)的运算为________。
A. ln(x)
B. log(x)
C. Ln(x)
D. log10(x)
（3）已知a=0:4，b=1:5，下面的运算表达式出错的为_______。
A. a+b
B. a./b
C. a ’*b
D. a*b
2.用“from:step:to ”方式和linspace 函数分别得到从0到4π，步长为0.4π的变量x1和从0到4π分成10点的变量x2。
3.输入矩阵A =123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
，使用全下标方式取出元素“3”
，使用单下标方式取出元素“8”，取出后2行子矩阵块，使用逻辑矩阵方式取出1379⎡⎤⎢⎥⎣⎦。
4.输入A 为3×3的魔方阵，Ｂ为3×3的单位阵，由小矩阵组成3×6的大矩阵C 和6×3的大矩阵D ，将D 矩阵的最后１行构成小矩阵E 。
5.输入字符串变量a 为“hello ”，将a 的每个字符向后移4个，例如“h ”变为“l ”，然后再逆序排放赋给变量b 。
6.求矩阵1234⎡⎤⎢⎥⎣⎦
的转置矩阵、逆矩阵、矩阵的秩、矩阵的行列式值、矩阵的三次幂、矩阵的特征值和特征向量。
7.求解方程组12341241
23412342x 3x x 2x 8x 3x x 6
x x x 8x 77x x 2x 2x 5
-++=⎧⎪++=⎪⎨-++=⎪⎪+-+=⎩ 8.计算数组A =123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
，B =111222333⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦的左除、右除以及点乘和点除。
9.
计算函数2()sin(4)-=t f t t 的值，其中t 的范围为0~2π，步长取0.1π；z 为0.707；1()f t 为()0>=f t 的部分，计算1()f t 的值。
作业题：3、5、9（写到作业纸上，待通知交时再交上来）
其余为练习题（大家上机练习一下，课堂上可能会提问）。