《数值计算方法》复习试题 一、填空题： 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 )次代数精度。
重要值的计算方法 以综合数值表示植物物种在群落中的相对重要值。 重要值=相对多度+相对频度+相对显著度 或，重要值=（相对多度+相对频度+相对显著度)/3 补充： 针对乔木而言：重要值=（相对密度【即相对多度】+相对频度+相对显著度【即相对优势度】）/3 针对灌草而言：重要值=（相对密度【即相对多度】+相对频度+相对盖度【即相对优势度】）/3 注： 频度：是指一个种在所作的全部样方中出现的频率.相对频度指某种在全部样方中的频度与所有种频度和之比。 相对频度=(该种的频度/所有种的频度总和)×100% 显著度【优势度】：指样方内某种植物的胸高断面积除以样地面积。 相对显著度【相对优势度】=(样方中该种个体胸面积和/样方中全部个体胸面积总和)×100% 密度（D）=某样方内某种植物的个体数/样方面积 相对密度（RD）=（某种植物的密度/全部植物的总密度）×100 =（某种植物的个体数/全部植物的个体数）×100 盖度（cover degree，或coverage）指的是植物地上部分垂直投影面积占样地面积的百分比，即投影盖度。后来又出现了“基盖度”的概念，即植物基部的覆盖面积。对于草原群落，常以离地面1英寸（2.54cm）高度的断面计算；对森林群落，则以树木胸高（1.3m处）断面积计算。基盖度也称真盖度。乔木的基盖度特称为显著度（dominant）。盖度可分为种盖度（分盖度）、层盖度（种组盖度）、总盖度（群落盖度）。林业上常用郁闭度来表示林木层的盖度。通常，分盖度或层盖度之和大于总盖度。群落中某一物种的分盖度占所有分盖度之和的百分比，即相对盖度。某一物种的盖度占盖度最大物种的盖度的百分比称为盖度比（cover ratio）。
数值计算方法试题一 一、 填空题（每空1分，共17分） 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。 2、迭代格式)2(2 1-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。 3、已知?????≤≤+-+-+-≤≤=31)1()1()1(211 0)(2 33x c x b x a x x x x S 是三次样条函数，则 a =( )， b =（ ）， c =（ ）。 4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数，则 ∑== n k k x l )(( )，∑== n k k j k x l x 0 )(( )，当2≥n 时= ++∑=)()3(20 4x l x x k k n k k ( )。 5、设1326)(2 47+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=?07 f 。 6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。 7、{}∞ =0)(k k x ?是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族， 其中1)(0=x ?，则?=1 04)(dx x x ? 。 8、给定方程组?? ?=+-=-2211 21b x ax b ax x ，a 为实数，当a 满足 ，且20<<ω时， SOR 迭代法收敛。 9、解初值问题00(,)()y f x y y x y '=??=?的改进欧拉法?????++=+=++++)],(),([2),(] 0[111] 0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是 阶方法。 10、设 ?? ????????=11001a a a a A ，当∈a （ ）时，必有分解式T LL A =，其中L 为下三角阵，当其对角线元素)3,2,1(=i l ii 满足（ ）条件时，这种分解是唯 一的。 二、 二、选择题（每题2分） 1、解方程组b Ax =的简单迭代格式g Bx x k k +=+) () 1(收敛的充要条件是（ ）。 （1）1)(A ρ, (4) 1)(>B ρ 2、在牛顿-柯特斯求积公式： ? ∑=-≈b a n i i n i x f C a b dx x f 0 )() ()()(中，当系数) (n i C 是负值时，
《计算方法》期末考试试题 一 选 择（每题3分，合计42分） 1. x* ＝ 1.732050808，取x ＝1.7320，则x 具有 位有效数字。 A 、3 B 、4 C 、5 D 、6 2. 取7 3.13≈（三位有效数字），则 ≤-73.13 。 A 、30.510-? B 、20.510-? C 、10.510-? D 、0.5 3. 下面_ _不是数值计算应注意的问题。 A 、注意简化计算步骤，减少运算次数 B 、要避免相近两数相减 C 、要防止大数吃掉小数 D 、要尽量消灭误差 4. 对任意初始向量)0(x 及常向量g ，迭代过程g x B x k k +=+)() 1(收敛的充分必要条件是_ _。 A 、11< B B 、1<∞ B C 、1)(数值计算方法试题及答案
数值计算方法试题一 一、填空题（每空1分，共17分） 1、如果用二分法求方程在区间内的根精确到三位小数，需对分（）次。 2、迭代格式局部收敛的充分条件是取值在（）。 3、已知是三次样条函数，则 =( )，=（），=（）。 4、是以整数点为节点的Lagrange插值基函数，则 ( )，( )，当时( )。 5、设和节点则 和。 6、5个节点的牛顿-柯特斯求积公式的代数精度为，5个节点的求积公式最高代数精度为。 7、是区间上权函数的最高项系数为1的正交多项式族，其中，则。 8、给定方程组，为实数，当满足，且时，SOR迭代法收敛。 9、解初值问题的改进欧拉法是 阶方法。 10、设，当（）时，必有分解式，其中为下三角阵，当其对角线元素满足（）条件时，这种分解是唯一的。 二、二、选择题（每题2分） 1、解方程组的简单迭代格式收敛的充要条件是（）。（1）, (2) , (3) , (4) 2、在牛顿-柯特斯求积公式：中，当系数是负值时，公式的稳定性不能保证，所以实际应用中，当（）时的牛顿-柯特斯求积公式不使用。 （1），（2），（3），（4）， （1）二次；（2）三次；（3）四次；（4）五次 4、若用二阶中点公式求解初值问题，试问为保证该公式绝对稳定，步长的取值范围为（）。 (1), (2), (3), (4)
三、1、 2、（15 (1)(1) 试用余项估计其误差。 （2）用的复化梯形公式（或复化 Simpson公式）计算出该积分的近似值。 四、1、（15分）方程在附近有根，把方程写成三种不同的等价形式（1）对应迭代格式；(2)对应迭代格式；（3）对应迭代格式。判断迭代格式在的收敛性，选一种收敛格式计算附近的根，精确到小数点后第三位。选一种迭代格式建立Steffensen迭代法，并进行计算与前一种结果比较，说明是否有加速效果。 2、（8分）已知方程组，其中 ， （1）（1）列出Jacobi迭代法和Gauss-Seidel迭代法的分量形式。 （2）（2）求出Jacobi迭代矩阵的谱半径，写出SOR 迭代法。 五、1、（15分）取步长，求解初值问题用改进的欧拉法求的值；用经典的四阶龙格—库塔法求的值。 2、（8分）求一次数不高于4次的多项式使它满足 ,,,, 六、（下列2题任选一题，4分） 1、1、数值积分公式形如 （1）（1）试确定参数使公式代数精度尽量高；（2）设，推导余项公式，并估计误差。 2、2、用二步法 求解常微分方程的初值问题时，如何选择参数使方法阶数尽可能高，并求局部截断误差主项，此时该方法是几阶的。 数值计算方法试题二 一、判断题：（共16分，每小题２分） １、若是阶非奇异阵，则必存在单位下三角阵和上三角阵，使唯一成立。（）
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 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)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 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 ，为了减少舍入误差，应将表达式1999 2001-
《数值计算方法》复习试题 一、填空题: 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．找两地的经度差： (1)如果已知地和要求地同在东经或同在西经，则： 经度差=经度大的度数—经度小的度数 (2)如果已知地和要求地不同是东经或西经，则： 经度差=两经度和（和小于180°时） 或经度差=（180°—两经度和）。（在两经度和大于180°时） 2．把经度差转化为地方时差，即： 地方时差=经度差÷15°/H 3．根据要求地在已知地的东西位置关系，加减地方时差，即：要求点在已知点的东方，加地方时差；如要求点在已知点西方，则减地方时差。 二．东西位置关系的判断： （1）同是东经，度数越大越靠东。即：度数大的在东。 （2）是西经，度数越大越靠西。即：度数大的在西。 （3）一个东经一个西经，如果和小180°，东经在东西经在西；如果和大于180°，则经度差=（360°—和），东经在西，西经在东；如果和等于180，则亦东亦西。 三．应用举例： 1、固定点计算 【例1】两地同在东经或西经 已知：A点120°E,地方时为10：00，求B点60°E的地方时。 分析：因为A、B两点同是东经，所以，A、B两点的经度差=120°－60°=60° 地方时差=60°÷15°/H=4小时 因为A、B两点同是东经，度数越大越靠东，要求B点60°E比A点120°E小，所以，B点在A点的西方，应减地方时差。 所以，B点地方时为10：00—4小时=6：00 【例2】两地分属东西经 A、已知:A点110°E的地方时为10:00,求B点30°W的地方时. 分析:Ａ在东经，Ｂ在西经，110°+30°=140°＜180°，所以经度差=140°，且A点东经在东，B点西经在西，A、B两点的地方时差=140°÷15°/H=9小时20分，B点在西方， 所以，B点的地方时为10：00—9小时20分=00：40。 B、已知A点100°E的地方时为8：00，求B点90°W的地方时。 分析：A点为东经，B点为西经，100°+90°=190°＞180°， 则A、，B两点的经度差=360°—190°=170°，且A点东经在西，B点西经在东。 所以，A、B两点的地方时差=170°÷15°/H=11小时20分，B点在A点的东方， 所以B点的地方时为8：00+11小时20分=19：20。 C、已知A点100°E的地方 8：00，求B点80°W的地方时。 分析：A点为100°E，B点为80°W，则100°+80°=180°，亦东亦西，即：可以说B点在A点的东方，也可以说B点在A点的西方，A，B两点的地方时差为180÷15/H=12小时。所以B点的地方时为8：00+12小时=20：00或8：00—12小时，不够减，在日期中借一天24小时来，即24小时 +8：00—12小时=20：00。 2、变化点计算 【例1】一架飞机于10月1日17时从我国上海（东八区）飞往美国旧金山（西八区），需飞行14小时。到达目的地时，当地时间是（） A. 10月2日15时 B. 10月2日3时 C. 10月1日15时 D. 10月1日3时
数值分析复习要点 引论 1 数值计算研究的对象与特点 计算方法研究的对象是专门研究各种数学问题的计算机解法（数值解法）,包括方法的构造和求解过程的理论分析及软件实现,包括方法的收敛性、稳定性以及误差分析等. 计算方法即具有纯数学的抽象性与严密性的特点,又具有应用的广泛性与实验的技术性特点. 2 误差的概念 2.1 误差的来源 模型误差：数学模型的解与实际问题的解之间出现的误差,称为模型误差. 测量误差：在测量具体数据时产生的误差称为测量误差. 截断误差：数学模型的准确解与数值方法的准确解之间的误差称为截断误差. 舍入误差：由于计算机字长的限制而产生的误差,称为舍入误差. 2.2误差的度量 (1).绝对误差与绝对误差限 (2).相对误差与相对误差限 (3). 有效数字 2.3 误差的传播 和、差的误差限 不超过各误差限的和. 积、商的相对误差限不超过各相对误差限的和. 3 数值计算的若干原则 避免两相近数相减和绝对值太小的除数、简化计算步骤、使用数值稳定的算法 方程求根 1 二分法 用二分法求方程0)(=x f 的实根*x 的近似值,其主要思想是:将含有根*x 的隔离区间二分,通过判断二分点与边界点函数值的符号,逐步对半缩小隔离区间,直到缩小到满足精度要求为止,然后取 最后二分区间的中点为根*x 的近似值. 2 迭代法 一般地,为了求一元非线性方程 0)(=x f 的根,可以先将其转换为如下的等价形式 ()x x ?=然后构造迭代公式.()k k x x ?=+1 2,1,0=k 3 收敛性和收敛速度 （收敛性基本定理）的条件和结论 收敛速度的快慢可用收敛阶来衡量.（收敛阶）设序列{}∞=0k k x 收敛到 *x ,并记误差
《计算方法》练习题一 练习题第1套参考答案 一、填空题 1． 14159.3=π的近似值3.1428，准确数位是（ 2 10- ）。 2．满足d b f c a f ==)(,)(的插值余项=)(x R （ ))((!2) (b x a x f --''ξ ） 。 3．设)}({x P k 为勒让德多项式，则=))(),((22x P x P （5 2 ）。 4．乘幂法是求实方阵（按模最大 ）特征值与特征向量的迭代法。 5．欧拉法的绝对稳定实区间是（ ]0,2[-）。 二、单选题 1．已知近似数,,b a 的误差限)(),(b a εε，则=)(ab ε（Ｃ ）。 A ．)()(b a εε Ｂ．)()(b a εε+ Ｃ．)()(b b a a εε+ Ｄ．)()(a b b a εε+ 2．设x x x f +=2 )(，则=]3,2,1[f （ Ａ ）。 Ａ．１ Ｂ．２ Ｃ．３ Ｄ．４ 3．设Ａ＝?? ? ? ??3113，则化Ａ为对角阵的平面旋转=θ（ Ｃ ） ． Ａ． 2π Ｂ．3π Ｃ．4π Ｄ．6 π 4．若双点弦法收敛，则双点弦法具有（Ｂ ）敛速． Ａ．线性 Ｂ．超线性 Ｃ．平方 Ｄ．三次 5．改进欧拉法的局部截断误差阶是（ Ｃ ）. A ．)(h o Ｂ．)(2 h o Ｃ．)(3 h o Ｄ．)(4 h o 三、计算题 1．求矛盾方程组：??? ??=-=+=+2 42321 2121x x x x x x 的最小二乘解。 2 212 212 2121)2()42()3(),(--+-++-+=x x x x x x x x ?， 由 0,021=??=??x x ? ?得：???=+=+9 629232121x x x x ， 解得14 9 ,71821== x x 。
山西大学计算机与信息技术学院实验报告 姓名XXX 学号XXX 专业班级（2012）计算机科学与技术 课程名称计算方法实验日期10.28 成绩指导老师李桂成批改日期 实验名称实验一方程求根 一、实验目的 用各种方法求任意实函数方程f(x)=0在自变量区间[a,b]上，或某一点附近的实根。并比较方法的优劣。 二、实验方法 (1)二分法 对方程f(x)=0在[a,b]内求根，将所给区间二分，在分点x= 2a b 判断是否f(x)=0，若是，则有根x=(b+a)/2。否则，继续判断是否f(a)* f(x)<0，若是，则令b=x，否则令a=x、重复此过程直至求出方程f(x)=0在[a,b]中的近似根为止。 (2)迭代法 将方程f(x)=0等价变换为x=φ（x）形式，并建立相应的迭代公式x k+1=φ（x k）。 (3)牛顿法 若已知方程f(x)=0的一个近似根x0，是函数f(x)在点x0附近可用一阶泰勒多项式 p1(x)=f(x )+f`(x0)(x-x0)来近似，因此方程f(x)=0可近似表示为f(x0)+f`(x0)(x-x0)=0。设f`(x0)≠0，则x=x0-f(x0)/f`(x0)。则x作为原方程新的近似根x1，然后x1将作为x0带入上式。迭代公式为x k+1=x k -f(x k)/f`(x k)。 三、实验内容 （1）在区间[0，1]上用二分法求方程ex+10x-2=0的近似根，要求误差不超过0.5x10-3。 （2）取初值x0=0，用迭代公式，x k+1= 10 k x e- 2，（k=0,1,2···）求方程e x+10x-2=0的 近似根。要求误差不超过0.5x10-3。 （3）取初值x0=0，用牛顿迭代法求方程e x+10x-2=0（k=0,1,2···）的近似根。要求误差不超过0.5x10-3。
计算方法模拟试题 一、 单项选择题（每小题3分，共15分） 1．近似值210450.0?的误差限为（ ）。 A ． 0.5 B. 0.05 C ． 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为（ ）。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足（ ）时，则存在唯一单位下三角阵L 和上三角阵R ，使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ，则=∞A （ ）。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ，则乘幂法计算公式1e =（ ）。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题（每小题3分，共15分） 1. 14159.3=π，具有4位有效数字的近似值为 。 2. 已知近似值21,x x ，则=-?)(21x x 。 3．已知1)(2-=x x f ，则差商=]3,2,1[f 。 4．雅可比法是求实对称阵 的一种变换方法。