数值分析试题及答案汇总

数值分析试题及答案汇

总

https://www.docsj.com/doc/a010621664.html,work Information Technology Company.2020YEAR

数值分析试题

一、填空题（2 0×2′） 1.

??

????-=?

?????-=32,1223X A 设x =0.231是精确值x *=0.229的近似值，则x 有 2 位有效数字。

2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]=

0 。

3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____，

‖AX ‖∞≤_15_ __。

4. 非线性方程f (x )=0的迭代函数x =

(x )在有解区间满足 |

’(x )| <1 ，则使用该迭代

函数的迭代解法一定是局部收敛的。

5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。

6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商

公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。

7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n

i i x a 0)( 1 ；所以当系数

a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使

20的近似值的相对误差小于0.1%，至少要取 4 位有效数字。

9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收

敛于方程组的精确解x *的充分必要条件是 (B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。

x 0 0.5 1 1.5 2 2.5 y =f (x )

-2

-1.75

-1

0.25

2

4.25

11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。

12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差

r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。

13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且

f (x )的二阶导数不变号，则初始点x 0的选取依据为 f(x0)f ”(x0)>0 。 14. 使用迭代计算的步骤为建立迭代函数、 选取初值 、迭代计算。 二、判断题（10×1′） 1、

若A 是n 阶非奇异矩阵，则线性方程组AX ＝b 一定可以使用高斯消元法求解。

( × ) 2、 解非线性方程f (x )=0的牛顿迭代法在单根x *附近是平方收敛的。 ( √ ) 3、 若A 为n 阶方阵，且其元素满足不等式

)

,...,2,1( 1n i a a n

i

j j ij ii =≥∑≠=

则解线性方程组AX ＝b 的高斯——塞德尔迭代法一定收敛。 ( × ) 4、 样条插值一种分段插值。 ( √ )

5、 如果插值结点相同，在满足相同插值条件下所有的插值多项式是等价的。 ( √ )

6、

从实际问题的精确解到实际的计算结果间的误差有模型误差、观测误差、截断

误差及舍入误差。 ( √ ) 7、 解线性方程组的的平方根直接解法适用于任何线性方程组AX ＝b 。 ( × ) 8、

迭代解法的舍入误差估计要从第一步迭代计算的舍入误差开始估计,直到最后一

步迭代计算的舍入误差。 ( × ) 9、

数值计算中的总误差如果只考虑截断误差和舍入误差，则误差的最佳分配原则

是截断误差＝舍入误差。 ( √ ) 10、插值计算中避免外插是为了减少舍入误差。 ( × ) 三、计算题（5×10′）

1、用列主元高斯消元法解线性方程组。

???

??=++-=+--=+-112123454 3

21321321x x x x x x x x x 解答：

（1，5，2）最大元5在第二行，交换第一与第二行：

???

??=++-=+--=+-1124 123453

21321321x x x x x x x x x L 21=1/5=0.2,l 31=2/5=0.4 方程化为：

???

??=--=+--=+-8.152.06.26.1 0.4 2.0123453

232321x x x x x x x （-0.2,2.6）最大元在第三行，交换第二与第三行：

???

??-=+-=--=+-6.1 0.4 2.08.152.06.2123453

232321x x x x x x x L32=-0.2/2.6=-0.076923,方程化为：

???

??-==--=+-38466.00.38462 8.152.06.2123453

32321x x x x x x

回代得：

???

??-===00010.1 99999.500005.33

21x x x 2、用牛顿——埃尔米特插值法求满足下列表中插值条件的四次插值多项式P 4(x )，并写出其截断误差的表达式(设f (x )在插值区间上具有直到五阶连续导数)。

解答： 做差商表

P4(x)=1-2x-3x(x-1)-x(x-1)(x-1)(x-2) R4(x)=f(5)(ξ)/5!x(x-1)(x-1)(x-2)(x-2)

3、对下面的线性方程组变化为等价的线性方程组，使之应用雅克比迭代法和高斯——赛德尔迭代法均收敛，写出变化后的线性方程组及雅克比迭代法和高斯——赛德尔迭代法的迭代公式，并简单说明收敛的理由。

???????=-+-=-+=+-=+-3 38 4 65 1 2321432431421x x x x x x x x x x x x

解答：

交换第二和第四个方程，使系数矩阵为严格对角占优：

雅克比迭代公式：

《计算机数学基础(2)》数值分析试题

一、单项选择题(每小题3分，共15分)

1. 已知准确值x *与其有t 位有效数字的近似值x ＝0.0a 1a 2…a n ×10s (a 1≠0)的绝对

误差∣x *－x ∣≤( )．

(A) 0.5×10 s －1－t (B) 0.5×10 s －t (C) 0.5×10s ＋1－t (D) 0.5×10 s ＋t

2. 以下矩阵是严格对角占优矩阵的为( )．

(A) ????

?????

???------21001210012100

12，

(B)?

?

???

????

???2100141101410125 (C) ??

???????

???--2100

14121241

0125 (D) ??

???

??

??

???-513

114120141112

4 ??????

?=+-=-+=-+-=+-6

5 8 4 3 3 1 2431432321421x x x x x x x x x x x x ??????

?=+-=-+=-+-=+-6

5 8 4 3 3 1

2431432321421x x x x x x x x x x x x

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

第一章 绪论 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

北师大网络教育 数值分析 期末试卷含答案

注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考北师大网络教育——数值分析——期末考试卷与答案 一．填空题（本大题共4小题，每小题4分，共16分） 1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为 。 2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。 3.设110111011A -????=--????-??，233x ?? ??=?? ???? ，则1A ＝ ，1x = 。 4. 1n +个节点的高斯求积公式的代数精确度为 。 二．简答题（本大题共3小题，每小题8分，共24分） 1. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？ 2. 什么是不动点迭代法？()x ?满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ?的不动点？ 3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥ ，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。 三．求一个次数不高于3的多项式()3P x ，满足下列插值条件： i x 1 2 3 i y 2 4 12 i y ' 3 并估计误差。（10分） 四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1 01 1I dx x =+? 。（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。（10分） 六．试用Doolittle 分解法求解方程组：

注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考 12325610413191963630 x x x -?????? ??????-=?????? ??????----?????? （10分） 七．请写出雅可比迭代法求解线性方程组1231231 23202324 812231530 x x x x x x x x x ++=?? ++=??-+=? 的迭代格式，并 判断其是否收敛？（10分） 八．就初值问题0(0)y y y y λ'=??=?考察欧拉显式格式的收敛性。（10分）

数值分析学期期末考试试题与答案(A)

期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩

清华大学高等数值计算(李津)实践题目一(共轭梯度CG法,Lanczos算法与MINRES算法)

高等数值计算实践题目一 1. 实践目的 本次计算实践主要是在掌握共轭梯度法，Lanczos 算法与MINRES 算法的基础上，进一步探讨这3种算法的数值性质，主要研究特征值特征向量对算法收敛性的影响。 2. 实践过程 （一）生成矩阵 （1）作5个100阶对角阵i D 如下： 1D 对角元：1,1,...,20,1+0.1(-20),21,...,100j j d j d j j ==== 2D 对角元：1,1,...,20,1+(-20),21,...,100j j d j d j j ==== 3D 对角元：,1,...,80,81,81,...,100j j d j j d j ==== 4D 对角元：,1,...,40,41,41,...,60,41+(60),61,...,100j j j d j j d j d j j =====-= 5D 对角元：,1,...,100j d j j == 记i D 的最大模特征值和最小模特征值分别为1i λ和i n λ，则i D 特征值分布有如下特点： 1D 的特征值有较多接近于i n λ，并且1/i i n λλ较小， 2D 的特征值有较多接近于i n λ，并且1/i i n λλ较大， 3D 的特征值有较多接近于1i λ，并且1/i i n λλ较大， 4D 的特征值有较多接近于中间模特征值，并且1/i i n λλ较大， 5D 的特征值均匀分布，并且1/i i n λλ较大 （2）随机生成10个100阶矩阵j M ： (100(100))j M fix rand = 并作它们的QR 分解，得j Q 和j R ，这样可得50个对称的矩阵T ij j i j A Q DQ =，其中i D 的对角元就是ij A 的特征值，若它们都大于0，则ij A 正定，j Q 的列就是相应的特征向量。结合（1）可知，ij A 都是对称正定阵。

数值分析期末考试复习题及其答案.doc

数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知,n=6 5.01021 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620* 21021,6,0,10325413.0-?=-=-=?=ε绝对误差限n k k X 2分 2. 已知?????=001A 220 - ???? ?440求21,,A A A ∞ （6分） 解： {},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=001A A T 420 ?? ?? ? -420?????001 220 - ?????440=?????001 080 ???? ?3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A 3. 设3 2 )()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2 解： ①Newton 迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3分

②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ，其中：? ??=1 3A ??? 22，??????-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收 敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --? ??--=-=ααααα21231A I B 2分 其特征方程为 0) 21(2)31(=----= -αλα ααλλB I 2分 即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(

数值分析试卷及其答案2

1、（本题5分）试确定7 22作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22=3.142857…=1103142857.0-? π=3.141592… 所以 3 12 10 2 110 21005.0001264.07 22--?= ?= <=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3 10 2 10005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：???? ? ??=????? ??????? ??--654131321 112321x x x ； 解 设???? ? ? ?????? ? ?????? ??===????? ? ?--11 1 11113 1321 11232312132 1 32 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,21527,25,2323121321- == - == -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23,97,910( ,)5 63, 7,4(== （3分） 3、（本题6分）给定线性方程组??? ? ? ??=++-=+-+=-+-=-+17722238231138751043214321 321431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为

清华大学贾仲孝老师高等数值分析报告第二次实验

高等数值分析第二次实验作业

T1.构造例子特征值全部在右半平面时, 观察基本的Arnoldi 方法和GMRES 方法的数值性态, 和相应重新启动算法的收敛性. Answer: （1） 构造特征值均在右半平面的矩阵A ： 根据实Schur 分解，构造对角矩阵D 由n 个块形成，每个对角块具有如下形式，对应一对特 征值i i i αβ± i i i i i S αββα-?? = ??? 这样D=diag(S 1,S 2,S 3……S n )矩阵的特征值均分布在右半平面。生成矩阵A=U T AU ，其中U 为 正交阵，则A 矩阵的特征值也均在右半平面。不妨构造A 如下所示： 2211112222 /2/2/2/2N N A n n n n ?-?? ? ? ?- ? = ? ? ? - ? ?? ? 由于选择初值与右端项：x0=zeros(2*N,1);b=ones(2*N,1); 则生成矩阵A 的过程代码如下所示： N=500 %生成A 为2N 阶 A=zeros(2*N); for a=1:N A(2*a-1,2*a-1)=a; A(2*a-1,2*a)=-a; A(2*a,2*a-1)=a; A(2*a,2*a)=a; end U = orth(rand(2*N,2*N)); A1 = U'*A*U; （2） 观察基本的Arnoldi 和GMRES 方法 编写基本的Arnoldi 函数与基本GMRES 函数，具体代码见附录。 function [x,rm,flag]=Arnoldi(A,b,x0,tol,m) function [x,rm,flag]=GMRES(A,b,x0,tol,m) 输入:A 为方程组系数矩阵，b 为右端项，x0为初值，tol 为停机准则，m 为人为限制的最大步数。 输出:x 为方程的解，rm 为残差向量，flag 为解是否收敛的标志。 外程序如下所示： e=1e-6; m=700;

北航2010-2011年研究生数值分析期末模拟试卷1-3

数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数=________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ，则],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(2 323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________. 4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2 x q ________. 5 设???? ??????=11001a a a a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当 其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的. 二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有? ∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？

数值分析试题A卷10.1

中国石油大学（北京）2009--2010学年第一学期 研究生期末考试试题A (闭卷考试) 课程名称：数值分析 注：计算题取小数点后四位 一、填空题（共30分，每空3分） 1、 已知x =是由准确数a 经四舍五入得到的近似值，则x 的绝对误差 界为_______________。 2、数值微分公式()() '()i i i f x h f x f x h +-≈ 的截断误差为 。 3、已知向量T x =，求Householder 变换阵H ，使(2,0)T Hx =-。 H = 。 4、利用三点高斯求积公式 1 1 ()0.5556(0.7746)0.8889(0)0.5556(0.7746)f x dx f f f -≈-++? 导出求积分 4 0()f x dx ?的三点高斯求积公式 。 5、4 2 ()523,[0.1,0.2,0.3,0.4,0.5]_____.f x x x f =+-= 若则 6、以n + 1个互异节点x k ( k =0,1,…,n ),（n >1）为插值节点的 Lagrange 插值基函数为l k (x)( k =0,1,…,n )，则 (0)(1)__________.n k k k l x =+=∑ 7、已知3()P x 是用极小化插值法得到的cos x 在[0,4]上的三次插值多项式，则3()P x 的 截断误差上界为3()cos ()R x x P x =-≤_________.

8、已知向量(3,2,5)T x =-，求Gauss 变换阵L ，使(3,0,0)T Lx =。L =_________. 9、设3 2 ()(7)f x x =-, 给出求方程()0f x =根的二阶收敛的迭代格式_________。 10、下面M 文件是用来求解什么数学问题的________________________. function [x,k]=dd （x0） for k=1:1000 x=cos (x0); if abs(x-x0)<, break end x0=x; end 二、（15分）已知矛盾方程组Ax=b ，其中11120,1211A b ???? ????==???????????? ， （1）用施密特正交化方法求矩阵A 的正交分解，即A=QR 。 （2）用此正交分解求矛盾方程组Ax=b 的最小二乘解。 三、（10分）已知求解线性方程组Ax=b 的分量迭代格式 1 (1) (1) ()1 +1 /, 121,,i n k k k i i ij j ij j ii j j i x b a x a x a i n n -++===-- =-∑∑（），， （1）试导出其矩阵迭代格式及迭代矩阵； （2）若11a A a ?? = ??? ，推导上述迭代格式收敛的充分必要条件。 四、（15分）（1）证明对任何初值0x R ∈，由迭代公式11 1sin ,0,1,2, (2) k k x x k +=+ = 所产生的序列{}0k k x ∞ =都收敛于方程1 1sin 2 x x =+ 的根。 （2）迭代公式11 21sin ,0,1,2, (2) k k k x x x k +=-- =是否收敛。 五、（15分）用最小二乘法确定一条经过原点(0,0)的二次曲线，使之拟合下列数据

数值分析期末试卷

数值分析2006 — 2007学年第学期考试 课程名称：计算方法 A 卷 考试方式：开卷[] 闭卷[V ] 半开卷[] IV 类 充要条件是a 满足 二、（18分）已知函数表如下 1?设 f(0) = 0, f (1) =16 , f( 2) =46，则 f [0,1]= ,f[0,1,2]二 2 ?设 AJ <2 -3 -1 ，则X ，A := A 1 1 j — 3 ?计算积分 xdx ，取4位有效数字。用梯形公式求得的近似值为 "0.5 （辛普森）公式求得的近似值为 ，用 Spsn 4?设f （x ）二xe x -3，求方程f （x ） =0近似根的牛顿迭代公式是 ，它的收 敛阶是 5 ?要使求积公式 1 1 [f (x)dx 拓一(0) + A , f (x 1)具有2次代数精度，则 捲= _________________ , 0 4 6 ?求解线性方程组 x 1 ax 2 = 4 , 12_3 （其中a 为实数）的高斯一赛德尔迭代格式收敛的 10 11 12 13 In x 2.3026 2.3979 2.4849 2.5649

三、(20分)构造如下插值型求积公式，确定其中的待定系数，使其代数精度尽可能高, 并指出所得公式的代数精度。 2 f (x)dx ： A o f (0) A f (1) A2f(2) o

X 2 4 6 8 y 2 11 28 40 五、（14分）为求方程X ’ -X 2 -1 =0在X o =1.5附近的一个根，将方程改写为下列等价 形式，并建立相应的迭代公式: 试问上述两种迭代公式在 x 0 =1.5附近都收敛吗？为什么？说明理由。 （1）X =1 ?丄，迭代公式 X 1 X k 1 = 1 - X k （2） X 2二1 ，迭代公式 X —1 2 （X k ）； X k 1

数值分析期末试题

数值分析期末试题 一、填空题（20102=?分） （1)设??? ? ? ??? ??---=28 3 012 251A ，则=∞ A ______13_______。 （2)对于方程组?? ?=-=-3 4101522121x x x x ，Jacobi 迭代法的迭代矩阵是=J B ?? ? ? ??05.25.20。 （3)3*x 的相对误差约是*x 的相对误差的 3 1倍。 （4）求方程)(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 。 （6）设n n ?矩阵G 的特征值是n λλλ,,,21 ，则矩阵G 的谱半径=)(G ρi n i λ≤≤1max 。 （7）已知?? ? ? ??=1021 A ，则条件数=∞ )(A Cond 9 （8）为了提高数值计算精度，当正数x 充分大时,应将)1ln(2 -- x x 改写为 )1ln(2 ++ -x x 。 （9）n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。 （10）拟合三点))(,(11x f x ，))(,(22x f x ，))(,(33x f x 的水平直线是)(3 1 3 1 ∑== i i x f y 。 二、（10分）证明：方程组? ?? ??=-+=++=+-1 211 2321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明：Jacobi 迭代法的迭代矩阵为 ???? ? ?????---=05 .05 .01015.05.00J B J B 的特征多项式为

数值分析期末试题

一、（8分）用列主元素消去法解下列方程组： ??? ??=++-=+--=+-11 2123454 321321321x x x x x x x x x 二、（10分）依据下列数据构造插值多项式：y(0)=1,y(1)= —2,y '(0)=1, y '(1)=—4 三、（12分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式并利用复化的梯形公式、复化的辛普生公式计算下列积分： ? 9 1dx x n=4 四、（10分）证明对任意参数t ，下列龙格－库塔方法是二阶的。 五、（14分）用牛顿法构造求c 公式，并利用牛顿法求115。保留有效数字五位。 六、（10分）方程组AX=B 其中A=????????? ?10101a a a a 试就AX=B 建立雅可比迭代法和高斯-赛德尔迭代法，并讨论a 取何值时 迭代收斂。 七、（10分）试确定常数A,B,C,a,使得数值积分公式?-++-≈2 2 ) (}0{)()(a Cf Bf a Af dx x f 有尽可能多的 代数精确度。并求该公式的代数精确度。 八、{6分} 证明： A ≤ 其中A 为矩阵，V 为向量. 第二套 一、（8分）用列主元素消去法解下列方程组： ??? ??=++=+-=+3 2221 43321 32132x x x x x x x x 二、（12分）依据下列数据构造插值多项式：y(0)=y '(0)=0, y(1)=y '(1)= 1,y(2)=1 三、（14分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式,并利用复化的梯形公式、 复化的辛普生公式及其下表计算下列积分： ?2 /0 sin πxdx ????? ? ? -+-+=++==++=+1 3121231)1(,)1(() ,(),()(2 hk t y h t x f k thk y th x f k y x f k k k h y y n n n n n n n n

数值分析整理版试题及复习资料

例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1） 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为

[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 ()()()()()()()()1 1 200110 1 1 2011000 1 210 1 ,11, ,3 1 23 ,,, ,3226 9,324 dx x dx xdx f x x dx f x x x dx ??????????==== ====++= =++= ????? 所以，法方程为 011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6a =

数值分析试题及答案

数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差 商公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以 当系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。

清华大学高等数值计算(李津)实践题目二(SVD计算及图像压缩)(包含matlab代码)

第1部分 方法介绍 奇异值分解（SVD ）定理： 设m n A R ?∈，则存在正交矩阵m m V R ?∈和n n U R ?∈，使得 T O A V U O O ∑??=?? ?? 其中12(,, ,)r diag σσσ∑=，而且120r σσσ≥≥≥>，(1,2, ,)i i r σ=称为A 的 奇异值，V 的第i 列称为A 的左奇异向量，U 的第i 列称为A 的右奇异向量。 注：不失一般性，可以假设m n ≥，（对于m n <的情况，可以先对A 转置，然后进行SVD 分解，最后对所得的SVD 分解式进行转置，就可以得到原来的SVD 分解式） 方法1：传统的SVD 算法 主要思想： 设()m n A R m n ?∈≥，先将A 二对角化，即构造正交矩阵1U 和1V 使得 110T B n U AV m n ?? =?? -?? 其中1200n n B δγγδ??? ???=?????? 然后，对三角矩阵T T B B =进行带Wilkinson 位移的对称QR 迭代得到:T B P BQ =。 当某个0i γ=时，B 具有形状12B O B O B ?? =? ??? ，此时可以将B 的奇异值问题分解为两个低阶二对角阵的奇异值分解问题；而当某个0i δ=时，可以适当选取'Given s 变换，使得第i 行元素全为零的二对角阵，因此，此时也可以将B 约化为两个低 阶二对角阵的奇异值分解问题。 在实际计算时，当i B δε∞≤或者() 1j j j γεδδ-≤+（这里ε是一个略大于机器精度的正数）时，就将i δ或者i γ视作零，就可以将B 分解为两个低阶二对角阵的奇异值分解问题。

2019年数值分析第二学期期末考试试题与答案A

卷）期末考试试卷（A2007学年第二学期考试科目：数值分析分钟考试时间：120 年级专业学号姓 名 题号一2二三0四总分 分）分，共10一、判断题（每小题210001?n）（ 1. 用计算机求时，应按照从小到大的顺序相加。1000n1n?219992001?为了减少误差2. ,应将表达式进行计算。（改写为）19992001?） （ 3. 用数值微分公式中求导数值时，步长越小计算就越精确。） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（系数矩阵及其演变方式有用迭代法解线性方程组时，5. 迭代能否收敛与初始向量的选择、） （关，与常数项无关。 分）二、填空题（每空2分，共36_________. ________，相对误差限为已知数a的有效数为0.01，则它的绝对误差限为1. 0?110??????????xA?Ax,0?21,x??5A?_____. 则设______，_____，2. ????21?????1?130????53f(x)?2x?4x?5x,f[?1,1,0]?f[?3,?2,?1,1,2,3]? 3. 已知则, . 331?)?Af(0)?Af(f(x)dx?Af(?)的代数精度尽量高，应使4. 为使求积公式321331?A?A?A?，此时公式具有，，次的代数精度。312 ?nA)(A的关系是 5. A阶方阵的谱半径与它的任意一种范数. (k?1)(k)BAX??N(k?XMX?0,1,2,)产时，使迭代公式用迭代法解线性方程组6. ??)k(X . 生的向量序列收敛的充分必要条件是

AX?BAL和上三角矩7. 使用消元法解线性方程组系数矩阵时，可以分解为下三角矩阵1 4?2??BAX?.A?LUU?A，则阵若采用高斯消元法解的乘积，即，其中??21??L?U?AX?B，则，______________；若使用克劳特消元法解_______________u?lu BAX?的大小关系为_____（选填：则____；若使用平方根方法解>与，，111111<，=，不一定）。 ??x?yy?8. 以步长为1的二阶泰勒级数法求解初值问题的数值解，其迭代公式为 ?y(0)?1?___________________________. 三、计算题（第1～3、6小题每题8分，第4、5小题每题7分，共46分） 32?x01??3x?xf(x)?2)(1, 1.在区间为初值用牛顿迭代法求方程内的根，要求以0证明用牛顿法解此方程是收敛的；（1）,xx,计算结果（2）给出用牛顿法解此方程的迭代公式，并求出这个根（只需计算21位）。取到小数点后4 2 2.给定线性方程组 x?0.4x?0.4x?1?312?0.4x?x?0.8x?2?321?0.4x?0.8x?x?3?312（1）分别写出用Jacobi和 Gauss-Seidel迭代法求解上述方程组的迭代公式； （2）试分析以上两种迭代方法的敛散性。

数值分析试题及答案汇总

数值分析试题及答案汇 总 TTA standardization office【TTA 5AB- TTAK 08- TTA 2C】

数值分析试题 一、填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =(x )在有解区间满足 |’(x )| <1 ，则使用该迭代函数 的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差 商公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当系数 a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 (B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。

清华大学高等数值分析作业李津1——矩阵基础

20130917题目 求证：在矩阵的LU 分解中，1 11n n T n ij i j j i j L I e e α-==+??=- ??? ∑∑ 证明： 在高斯消去过程中，假设0jj a ≠ ，若a=0，可以通过列变换使得前面的条件成立，这里不考虑这种情况。 对矩阵A 进行LU 分解，()() () ()()1 11 1111L M n M M M n ---=-=??-………… ， 其中()1n T n ij i j i j M j I e e α=+??=+ ??? ∑ ，i e 、j e 为n 维线性空间的自然基。 ()M j 是通过对单位阵进行初等变换得到， 通过逆向的变换则可以得到单位阵，由此很容易得到()M j 的逆矩阵为1n T n ij i j i j I e e α=+??- ???∑。故111n n T n ij i j n j i j L I e e I α-==+?? ??=- ? ? ????? ∏∑ 上式中的每一项均是初等变换，从右向左乘，则每乘一次相当于对右边的矩阵进行一次 向下乘法叠加的初等变换。由于最初的矩阵为单位阵，变换从右向左展开，因而每一次变换不改变已经更新的数据，既该变换是从右向左一列一列更新数据，故 11n n T n ij i j j i j L I e e α==+??=- ??? ∑∑。 数学证明：1n T ij i j i j e e α=+?? ???∑具有 ,0 00n j j A -?? ??? 和1,1000n j n j B -+-+?? ?? ? 的形式，且有 +1,-11,10000=000n j j n j n j A B --+-+???? ?????? ? 而1 1n n T ij i j j k i j e e α-==+?? ??? ∑∑具有1,1000n k n k B -+-+?? ???的形式，因此： 1 311111211121==n n n n n n T T T n ij i j n ij i j n ik i k j i j j i j k n i k n n T n i i n ik i i i k L I e e I e e I e e I e e I e ααααα---==+==+=-=+==+??????????????=---?? ? ? ? ? ? ? ? ???????????????????????=-- ? ? ?????∏∑∏∑∑∑∑∑……11211n n n T T k n ik i k k k i k e I e e α--===+????=- ?? ?????? ∑∑∑#

数值分析试题及答案解析

数值分析试题 一、 填空题（2 0×2′） 1. ??????-=? ?????-=32,1223X A 设x =0.231是精确值x *=0.229的近似值，则x 有 2 位有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用 该迭代函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所 以当系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于0.1%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 ρ(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。