2010数值分析试题及答案
2010级数值分析试题
A ；
4.已知 f ( x ) 2 x 3 4 x 5 ，写出以 2, 0, 1 为插值节点的二次拉格朗日插值多项式； 5.写出用变步长梯形求积公式计算积分 f ( x )dx的计算公式.
a xk
3
r
产生的序列 x k 收敛到 4
a ，使其收敛阶尽可能高，
并说明该迭代公式的收敛阶.
10 a 0 五、（本题满分15分）设 A b 10 b （ A 的行列式 0 a 5
2 det( A) 0），给定方程组 AX 1 1
要求：（1）写出计算公式；（2）画出算法框图. 七、（本题满分10分）确定求积公式
3
1
f ( x )dx A0 f (0) A1 f (1) A2 f (2)
中的待定参数 A0 ， A1 ，与 A2，使其代数精度尽可能高， 并指明该求积公式所具有的代数精度.
八、（本题满分5分）
给定等距节点
a b
二、（本题满分10分） 用高斯消去法求解方程组
a11 x1 a12 x 2 a13 x 3 a1n x n b1 a x a x a x a x b 21 1 22 2 23 3 2n n 2 a n1 x1 a n 2 x 2 a n 3 x 3 a nn x n bn
1.根据迭代法收敛的充分必要条件确定a,b的取值范围， 使求解上述方程组的雅可比迭代法收敛.
. 2.写出求解上述方程组的雅可比迭代公式.
3.若用高斯—塞德尔迭代法求解上述方程组，画出高斯 —塞德尔迭代法的算法框图.
2010(A)数值分析试卷
一. 选择题(每空2分, 共20分)1．设 是真值 的近似值，则有 ________位有效数字. 2.2. 设，则差商（均差）1)(3−+=x x x f =]3,2,1,0[f ________，且_________. =]4,3,2,1,0[f 3. 求方程()x f x =根的牛顿迭代格式是__________________.)(,),(),(10x l x l x l n "4. 是以0，1，…，n 为插值节点的Lagrange 插值基函数，则 ________________. =∑=ni ix il 0)(5. 牛顿—柯特斯求积公式的系数和()0n n k k C==∑__________________.6. 设方程组b Ax =,其中,则Jacobi 代法的迭代矩阵是⎥⎦⎤⎢⎣⎡−=51112．A 迭__________________,Gauss-Seidel 法的迭代矩阵是__________________. 7. 设是区间[上权函数为的最高项系数为1的正交多项式族,其中{∞0)(x k ϕ}]1,02)(x x =ρ1)(0=x ϕ,则=∫dx x x )(3102ϕ__________________,=)(1x ϕ__________________. 二. 判断题(每小题2分, 共20分.正确的打√,错误的打×)1. 梯形求积公式和复化梯形公式都是插值型求积公式_____（对或错）. ( )2. 若x 为n 维向量,则0>x . ( )3. 幂法是求矩阵所有特征值及特征向量的一种向量迭代法． ( )4. 用x +1近似表示x e 产生舍入误差． ( )第 A1 页 共 3 页40194x =* 2.40315x =5. n 个求积节点的插值型求积公式的代数精度为n . ( )6. 321.750有5位有效数字,其误差限31021−×≤. ( )7. 求解微分方程初值问题的二阶龙格—库塔公式的局部截断误差为)(2h o .( )8. 高斯型求积公式的代数精度为12∑∫=≈nk k k b a x f A dx x f 0)()(+n . ( ) 9. 对0的充要条件是A 的某种算子范数lim ,=∈∀∞→×m m n n A R A 1<A .（ ）10. 方程组b Ax =得系数矩阵A 的条件数刻画了解对初始数据的灵敏程度，即A 的条件数越大，方程组的病态程度越严重． ( ) 三. 计算题(每小题10分, 共40分)1. 用矩阵的直接三角分解法（）解方程组LU A =b AX =：。
东北大学10数值分析B(研)答案
为什么？ 解 由于 f ( x, y ) xye y 关于 y 满足 Lipschitz 条件， 2分 5分 的差分公式：
所以，改进 Euler 法收敛。
h y n 1 y n 4 (3k1 k 2 ) k f (x , y ) n n 1 k 2 f ( x n 2h, y n 2hk1 ) y0
1 3 1 3 ， x2 2 6 2 6
4分
1 1 1 3 1 3 ) f( )] 积分公式为： f ( x)dx [ f ( 0 2 2 6 2 6
6分
解得： a 1 / 3, b 13 / 18 0.7222 ， 拟合曲线为: y
13 2 1 x 18 3
y ( x n 1 ) y ( x n ) y ( x n )h
y n hf n
h2 h3 y ( x n ) y ( x n ) O(h 4 ) 2 6
5分
h 2 f n f n ( fn ) 2 x y
2 fn 2 f n 2 f n f n f h3 2 f n [ 2 2 fn fn ( n ) 2 f h ] O(h 4 ) 2 6 x xy x y y y
。
解
3 2 4.(6 分)设 xk 1 xk axk bxk c, k 0,1,2,... 是求方程根 1 的迭代法，试确定
1/ 3 1/ 3 0 0 1 / 3 ,所以 B 由于 B 13 1/ 2 1/ 4 0
1 1 9 所以， H 3 ( x) ( x 2)( x 2 4 x 1) x 3 3x 2 x 1 2 2 2 9． 分）给定离散数据 （7
光机所2010数值分析试题
光机所2010数值分析试题一、(24分)(1) 设x j 为互异节点，试证明()∑=≡n j k j k j x x l x 0，k =0, …, n(2) 证明 ()()()∑==≡-n j j kj n k x l x x0,,1,0 (3) 求一个次数不高于4的多项式()x p 4，使它满足()()00044='=p p ，()()11144='=p p ，()124=p .二、(24分)(1) 什么叫求积公式的m 次代数精度？(2) 确定下述数值积分的参数a ，使其代数精度尽量高，并说明该公式有几次代数精度()()()()h f A f A h f A dx x f h h 101220++-=--⎰(3) 已知411=x ，212=x ，433=x ，推导以这三个点作为求积节点在[0, 1]上的插值型求积函数，并用所求公式计算⎰102dx x 三、(12分) 已知函数1, x, 312-x 在[-1, 1]上两两正交，并求一个三次多项式，使其与这三个函数两两正交。
四、解下列线性方程组（用某种方法解方程组，貌似主要考的是矩阵计算）五、(10分) 已知 ||x || = ||y ||, x ≠y ，证明存在初等反射阵变换H ，使Hx = y ，若x = (1, 1, 0)，y =(此处缺一系数) ||x|| e 1，从计算效果考虑应该取正号还是负号，并计算相应的矩阵H六、(8分) 证明1 – x - sin x = 0在[0, 1]内有一个根，使用二分法求误差不大于41021-⨯的解需迭代二分多少次？七、(12分) 设有迭代格式 ()()g Rx x k k +=-1, k = 1, 2, … 其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=05.02/15.005.02/15.00R ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=5.015.0g 试证该迭代格式收敛。
取初始向量为 (0, 0, 0) ，计算x。
福州大学2010-2011年数值分析考题及答案1
1、若向量 x (4, 2,3) ，则
T
x 2 =___ 29 _________
＝____ 6 ____，A 的
2、
1 1 A ， 则 A 的谱半径 -5 1
=____6____
3、 确定求积公式 尽量高,则 A0=_
1
1
f ( x)dx A0 f (1) A1 f (0) A2 f '(1) 中的待定参数，使其代数精度
0 2 0 5、设 B 2 1 2 ,试用平面旋转矩阵对矩阵 A 进行 QR 分解，其中 Q 为正交 0 2 1
矩阵，R 为上三角阵（8 分）
4
记A1 A, 先将A的第一列变得与e1平行 cos = 0 2 0,sin = 1 04 04 0 1 0 0 1 0 0 P A 2 P A1 1 12 12 0 0 0 1
3、
h 用二步法 yn1 yn [ f ( xn , yn ) f ( xn1 , yn1 )] 求解一阶常微分方程初值问题 2
y f ( x, y ) 问：如何选择参数 , 的值，才使该方法的阶数尽可能地高？写出 y ( x0 ) y0
此时的局部截断误差主项，并说明该方法是几阶的。 证明:局部截断误差为：
( x x )l ( x) 等于
i 0 i i
4
（ a ） 1 (c) 2 (d) 4
(a)
0
(b)
3、设 f ( x) 3x5 4 x 4 x 2 1 和节点 xk k / 2, k 0,1 则差商 f [ x0 , x1 x5 ] (a) 4 (b) 2 (c) 3 (d) 1 （ （ c ） c ）
《数值分析》参考答案
参考答案第1章一、选择题1. D2. C3. A4. B5. B二、填空题1. 函数题头 H1行 帮助信息 函数体 注释部分 函数题头2. nargin varargin3. A=rand(4)4. 单引号三、解答题1. for 语句和while 语句均可以实现循环执行的功能。
二者的区别在于，for 循环语句一般适用于已知道循环次数，而不知道循环运算的目标的问题，而while 循环语句则相反，一般适用于已知循环目标，而循环次数未知的问题。
2. 程序如下：function [highavg,weightavg]=avg_high_weight(varargin) n=length(varargin); highsum=0; weightsum=0; for i=1:n highsum=highsum+varargin{i}(1);weightsum=weightsum+varargin{i}(2);endhighavg=highsum/n; weightavg=weightsum/n;第2章一、选择题1. A2. B3. A4. C5. D二、填空题1. 1.7 1.73 1.7322. 3 13. 5％4. 3三、解答题1. 解：1*()()nn x nxx x ε-≈-1***()()n nr nxx x x x x nnxxε---≈=()0.02r ne x n ==2数值分析2. 解：*1 1.1021x =有五位有效数字；*20.031x =有两位有效数字；*3385.6x =有四位有效数字；*47 1.0x =⨯有一位有效数字。
3. 解：（1）*******124124()()()()x x x x x x εεεε++≤++433111101010222---=⨯+⨯+⨯3*1.0510ε-=⨯=（2）*********123231113()()x x x x x x x x x ε⋅⋅≈⋅-+⋅****221233()()x x x x x x -+⋅-*0.197ε≈=（3）******2242244**2441(/) |()()|()x x x x x x x xx ε≤---****2224**44|()()|r r x x x x xxεε=-***224*4||[|()||()|]r r x e x e x x≤+331110100.0312256.4800.03156.480--⎡⎤⨯⨯⎢⎥=+⎢⎥⎢⎥⎣⎦5*10ε-≤=4. 解：33**34433()43r R RV Rππεπ-=*2**2R R R R R RRR R-++=⋅*223R R R RR-≈⋅*3R R R-=⋅1%=故*1300R R R-=5. 解：设Y =*27.983Y =，*31102Y Y δ-=-≤⨯，028Y =，*028Y =，*0000Y Y δ=-=*111282827.983100Y Y ⎛⎛⎫-=---⨯ ⎪⎝⎝⎭1100δ≤，**22111127.983100100Y Y Y Y ⎛⎛⎫-=-⨯--⨯ ⎪⎝⎝⎭**111()()100Y Y Y Y =---112100100100δδδ≤+=仿此可得：*100n n n Y Y δ-≤则*31001001001101002Y Y δδ--≤==⨯即计算100Y 的误差界不超过31102-⨯参考答案 36. 解：解方程25610x x -+=得：28282x =±±由第5题知27.983具有五位有效数字，故可取：1282827.98355.983x =++=21280.0178655.983x =-≈=7. 解：设正方形的边长为x ，则其面积为2y x =。
数值分析试题及答案..(优选)
一、单项选择题（每小题3分，共15分）1. 3.142和3.141分别作为π的近似数具有（ ）和（ ）位有效数字.A ．4和3B ．3和2C ．3和4D ．4和42. 已知求积公式()()211211()(2)636f x dx f Af f ≈++⎰，则A ＝（ ）A ． 16B ．13C ．12D ．233. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（ ）A ．()00l x ＝0，()110l x = B ．()00l x ＝0，()111l x =C ．()00l x ＝1，()111l x = D ．()00l x ＝1，()111l x =4. 设求方程()0f x =的根的牛顿法收敛，则它具有（ ）敛速。
A ．超线性B ．平方C ．线性D ．三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到的第3个方程（ ）.A ．232x x -+= B ．232 1.5 3.5x x -+=C ．2323x x -+= D ．230.5 1.5x x -=-单项选择题答案1.A2.D3.D4.C5.B二、填空题（每小题3分，共15分）1. 设TX )4,3,2(-=, 则=1||||X ，2||||X = .2. 一阶均差()01,f x x =3. 已知3n =时，科茨系数()()()33301213,88C C C ===，那么()33C = 4. 因为方程()420x f x x =-+=在区间[]1,2上满足 ，所以()0f x =在区间内有根。
5. 取步长0.1h =，用欧拉法解初值问题()211yy yx y ⎧'=+⎪⎨⎪=⎩的计算公式 .填空题答案1. 9和292.()()0101f x f x x x --3. 18 4. ()()120f f < 5. ()1200.11.1,0,1,210.11k k y y k k y +⎧⎛⎫⎪ ⎪=+⎪ ⎪=+⎨⎝⎭⎪=⎪⎩得 分 评卷人三、计算题（每题15分，共60分）1. 已知函数211y x =+的一组数据：求分段线性插值函数，并计算()1.5f 的近似值.计算题1.答案1. 解[]0,1x ∈，()1010.510.50110x x L x x --=⨯+⨯=---()12x L x -=-所以分段线性插值函数为()10.50.80.3x x L x x x ⎧-∈⎪=⎨-⎪⎩()1.50.8L =2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩（1） 写出雅可比迭代公式、高斯－塞德尔迭代公式；（2） 对于初始值()()00,0,0X =，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算()1X（保留小数点后五位数字）.计算题2.答案1.解 原方程组同解变形为 1232133120.10.20.720.10.20.830.20.20.84x x x x x x x x x =++⎧⎪=-+⎨⎪=++⎩雅可比迭代公式为()()()()()()()()()1123121313120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x +++⎧=++⎪⎪=-+⎨⎪=++⎪⎩(0,1...)m = 高斯－塞德尔迭代法公式()()()()()()()()()1123112131113120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x ++++++⎧=++⎪⎪=-+⎨⎪=++⎪⎩ (0,1...)m =用雅可比迭代公式得()()10.72000,0.83000,0.84000X =用高斯－塞德尔迭代公式得()()10.72000,0.90200,1.16440X =3. 用牛顿法求方程3310x x --=在[]1,2之间的近似根 （1）请指出为什么初值应取2？（2）请用牛顿法求出近似根，精确到0.0001.计算题3.答案4. 写出梯形公式和辛卜生公式，并用来分别计算积分1011dx x +⎰.计算题4.答案确定下列求积公式中的待定系数，并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案一、 填空（共20分，每题2分）1. 设2.3149541...x *=，取5位有效数字，则所得的近似值x= .2.设一阶差商()()()21122114,321f x f x f x x x x --===---，()()()322332615,422f x f x f x x x x --===--则二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 则2||||X = ，=∞||||X 。
数值分析试卷及答案
数值分析试卷及答案数值分析试卷一、选择题（共10题，每题2分，共计20分）1. 数值分析的研究内容主要包括以下哪几个方面？A. 数值计算方法B. 数值误差C. 数值软件D. 数学分析答：A、B、C2. 下列哪种方法不属于数值积分的基本方法？A. 插值法B. 微积分基本公式C. 数值微积分D. 数值积分公式答：A3. 数值积分的目的是求解什么？A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答：D4. 数值微分的目的是求解什么？A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答：A5. 数值微分的基本方法有哪几种？A. 前向差分B. 后向差分C. 中心差分D. 插值法答：A、B、C6. 用数值方法求解方程的基本方法有哪几种？A. 迭代法B. 曲线拟合法C. 插值法D. 数值积分法答：A、B、C7. 用迭代法求方程的根时，当迭代结果满足何条件时可停止迭代？A. 当迭代结果开始发散B. 当迭代结果接近真实解C. 当迭代次数超过一定阈值D. 当迭代结果在一定范围内波动答：B8. 下列哪种插值方法能够确保经过所有给定数据点？A. 拉格朗日插值B. 牛顿插值C. 三次样条插值D. 二次插值答：A、B、C9. 数值解线性方程组的基本方法有哪几种？A. 直接法B. 迭代法C. 插值法D. 拟合法答：A、B10. 下列哪种方程求解方法适用于非线性方程？A. 直接法B. 迭代法C. 插值法D. 曲线拟合法答：B二、填空题（共5题，每题4分，共计20分）1. 数值积分的基本公式是_________。
答：牛顿-科特斯公式2. 数值微分的基本公式是_________。
答：中心差分公式3. 数值积分的误差分为_________误差和_________误差。
答：截断、舍入4. 用插值法求解函数值时，通常采用_________插值。
答：拉格朗日5. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析习题集及答案
数值分析习题集（适合课程《数值方法A 》和《数值方法B 》）长沙理工大学第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少6. 设028,Y =按递推公式1n n Y Y -=…)计算到100Y .(五位有效数字),试问计算100Y 将有多大误差7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字.8. 当N 充分大时,怎样求211Ndx x +∞+⎰9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝210. 设212S gt =假定g 是准确的,而对t 的测量有±秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大这个计算过程稳定吗12.计算61)f =,1.4≈,利用下列等式计算,哪一个得到的结果最好3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量 a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xxx ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3. 给出f (x )=ln x 的数值表用线性插值及二次插值计算ln 的近似值.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nkkj j j x l x xk n =≡=∑ii)()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少9. 若2nn y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i) 若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限.19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.24. 给定数据表如下:试求三次样条插值并满足条件 i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小r 是否唯一 9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式. 10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差. 15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积19. 用许瓦兹不等式估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26. 用最小二乘法求一个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均方误差.28. 在某化学反应里,根据实验所得分解物的浓度与时间关系如下:用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: (1)101()()(0)()hh f x dx A f h A f A fh --≈-++⎰;(2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰; (4),6n =.3. 直接验证柯特斯公式具有5次代数精度.4. 用辛普森公式求积分1xedx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2baf f x dx b a f a b a 'η=-+-⎰;(2)2()()()()()2baf f x dx b a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰.6. 证明梯形公式和辛普森公式当n →∞时收敛到积分()ba f x dx ⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)8.1xedx-⎰,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nnnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =,和处的导数值,并估计误差.()f x 的值由下表给出:第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
数值分析习题(含答案)
第一章 绪论姓名 学号 班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字？（有效数字的计算） 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。
3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字？（有效数字的计算）解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 解：已知δ=-**xx x ，则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。
数值分析习题集及答案[1](精)
](https://img.taocdn.com/s3/m/7eabbe90cc22bcd126ff0c65.png)
数值分析习题集（适合课程《数值方法A 》和《数值方法B 》）长沙理工大学第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式(3.3)求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少?6. 设028,Y =按递推公式1n n Y Y -= ( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字≈27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2?10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b c s a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n nn n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj jj x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少? 9. 若2nn y =,求4n y ∆及4n y δ. 10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆. 12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()b aS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =. 3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式. 4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式. 5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nTx 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hhf x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10x e dx-⎰并计算误差. 5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2baf f x dx b a f b b a 'η=---⎰; (3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n n nnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
数值分析试题(卷)和答案解析
试题__2009___年~__2010___年第 一学期课程名称： 数值分析 专业年级： 2009级（研究生） 考生学号： 考生： 试卷类型： A 卷 √ B 卷 □ 考试方式: 开卷 √ 闭卷 □………………………………………………………………………………………………………一. 填空题（本大题共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 =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。
（10分） 六．试用Doolittle 分解法求解方程组：12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦（10分） 七．请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式，并判断其是否收敛？（10分）八．就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
204-测试试卷-2010年数值分析考试题及解答
,
―――――12 分
余项为 R(x)=
f
(5) (x 5!
)
x2
(x
-1)2
(
x
-
2)
―――――14 分
解法 2: 用扩展的牛顿法―――视两个节点 0，1 为重合节点，并作差商表：
xi
yi
0
0
0
0
1
1
1
1
2
1
所以，得
一阶差 商
二阶差商
三阶差商
0
1
1
1
0
-1
0
-1
-0.5
0.25 ―――――8 分
N4 (x) = 0 + 0(x - 0) +1ﾴ (x - 0)2 + (-1) ﾴ x2 (x -1) + 0.25x2(x -1)2
=
16 15
S2
-
1 15
S1
=
0.494377
d0 = ( f ,j0 ) =
2 1
x ln xdx ﾻ 0.494377
再求 d1 = ( f ,j1) =
2 1
x ln x(x -1.5)dx （用龙贝格数值积分法求如下） a = 1, b = 2, f (x) = x
x ln x
T1
――――5 分
求解
1 2
0 1
0 0
y1 y2
=
1148
得
y
=
(
14,
-10,
-72) T
3 - 5 1 y3 20
――――8 分
求解
1 0
2 1
3 -4
x1 x2
(完整word版)三峡大学2010-2011数值分析试题
二、（10分）求下列超定线性方程组的最小二乘解. ⎪⎩⎪⎨⎧=+=+=+00212212121x x x x x x三、(10分) 在区间[]4,2上利用压缩映像原理判断迭代格式Λ,2,1,0,321=+=+k x x k k 的敛散性.四、（14分）设n n ij R a A ⨯∈=)(对称,顺序主子式),,2,1(0n i i Λ=≠∆则T LDL A =分解存在,其中L 为单位下三角形矩阵,D 为对角阵,试写出求方程组b Ax =解的计算步骤(用矩阵表示), 此法称为改进平方根法. 试用它求解方程组：⎩⎨⎧=+=+1710712752121x x x x五、(10分) 已知)(x f y =的函数值如下表: 28.44.54.79.77.888.74.52.10)(20181614121086420x f x利用所有数据,用复合辛普森（Simpson ）公式计算积分dx x f ⎰20)(的近似值.六、(10分) 取节点1,010==x x ,写出xe x y -=)(的一次插值多项式),(1x L 并估计插值误差.七、（10分）分别写出解线性方程组⎪⎩⎪⎨⎧=++-=-+=+-1874165321321321x x x x x x x x x 收敛的Jacobi 迭代格式和Seidel Gauss -迭代格式，并说明其收敛的理由.八、（10分）设初值问题:⎩⎨⎧=+='1)0(y yx y )10(≤≤x(1) 写出用欧拉（Euler ）方法、取步长1.0=h 解上述初值问题数值解的公式； (2) 写出用改进欧拉（Euler ）方法、取步长1.0=h 解上述初值问题数值解的公式.九、（10分）求线性代数方程组Ax b =的数值解法主要有矩阵的直接分解法(如LU 分解法、Crout 分解法、Cholesky 分解法等)和迭代法（如Jacobi 迭代法、Gauss-Seidel 迭代法）.请你简述求解线性代数方程组Ax b =的直接分解法和迭代法这两类方法的不同点和相同点.。
数值分析试题答案
数值分析试题答案一、选择题1. 以下哪个数值方法用于求解线性方程组？A. 牛顿法B. 高斯消元法C. 欧几里得算法D. 拉格朗日插值法答案：B2. 在数值分析中，舍入误差通常是由什么引起的？A. 人为计算错误B. 计算机表示数字的限制C. 测量误差D. 数据输入错误答案：B3. 插值和拟合的区别在于：A. 插值通过所有数据点，而拟合不通过B. 拟合通过所有数据点，而插值不通过C. 插值是线性的，拟合是非线性的D. 插值是精确的，拟合是近似的答案：A4. 以下哪种方法最适合求解非线性方程？A. 雅可比迭代法B. 牛顿-拉弗森方法C. 托马斯算法D. 布雷尔-史密斯算法答案：B5. 在数值分析中，条件数用于衡量什么？A. 方程组解的存在性B. 方程组解的唯一性C. 方程组解的稳定性D. 方程组解的精确性答案：C二、填空题1. 在数值分析中，__________误差指的是由于计算机舍入而产生的误差，而__________误差指的是由于数据不精确或截断而产生的误差。
答案：截断；舍入2. 线性方程组的矩阵表示为__________，其中A是系数矩阵，x是变量向量，b是常数向量。
答案：Ax = b3. 牛顿法求解非线性方程时，需要计算函数的__________。
答案：导数4. 拉格朗日插值法通过构建一个多项式来近似数据点，该多项式的每一段都与数据点的__________相匹配。
答案：切线5. 为了减少数值分析中的误差，通常采用__________方法来提高计算的精度。
答案：增量三、简答题1. 请简述高斯消元法的基本思想及其在求解线性方程组中的应用。
高斯消元法的基本思想是通过行变换将系数矩阵转化为阶梯形矩阵，进而简化方程组的求解过程。
在求解线性方程组时，首先将增广矩阵进行行变换，使得主元下方的元素为零，然后通过回代过程逐步求解出未知数。
2. 描述牛顿-拉弗森方法求解非线性方程的迭代过程。
牛顿-拉弗森方法是一种迭代求解非线性方程的方法。
数值分析习题集及答案
数值分析习题集（适合课程《数值方法A 》和《数值方法B 》）长沙理工大学第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式(3.3)求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少?6. 设028,Y =按递推公式1n n Y Y -= ( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字≈27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2?10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b c s a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n nn n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk k j j j x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2nn y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆. 12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()b aS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式. 10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()x f x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22.()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hhf x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10x e dx-⎰并计算误差. 5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰; (3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长. 10. 证明等式3524sin3!5!n nn n ππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准
2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准D第2页，共5页第3页，共5页第4页，共5页()()()1121212,,n n n n n n h y y k k k f t y k f t h y hk +⎧=++⎪⎪=⎨⎪=++⎪⎩，它是 2 阶方法。
二、解答下列各题(每小题12分，共24分) 1. 用LU分解法求解线性方程组1234230156715112133721182x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪=⎪ ⎪ ⎪- ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭；解：() 2 3 0 1 -53 -2 1 2 41 2 1 -2 41 2 -1 2 1 -11 -2 2 -3015671511|2 1 133 7201182LU A b ⎛⎫-⎛⎫ ⎪⎪ ⎪-⎪ ⎪=−−−→ ⎪- ⎪ ⎪ ⎪--⎝⎭ ⎪⎝⎭分解1 0 0 023 0 13 1 0 0 0 -2 1 21 2 1 0 0 0 1 -21 2 -1 1 0 0 0 154,,41L U y -⎛⎫⎛⎫⎛ ⎪ ⎪⎪ ⎪ === ⎪ ⎪ ⎪ ⎪ -⎝⎭⎝⎭⎝ 1-2 2-1,x ⎫⎛⎫⎪ ⎪⎪ ⎪=⎪ ⎪⎪ ⎪⎭⎝⎭.. LU 矩阵(或对应元素每算对两个给1分) 2. x19 25 31 38 44 y 19.0 32.3 49.0 73.3 97.8 求形如2y a bx =+的拟合函数及误差平方和，并估计变量y 在21x =处的值。
第5页，共5页解：令a cb ⎛⎫= ⎪⎝⎭，1361162519611144411936A ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭， ........................... 3分对应的正规方程组TTA Ac A y =为55327271.453277277699369321.5a b ⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭....................... 7分 解此方程得0.972578660.05003512a b ⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭...................................... 9分 即最终的拟合函数 (2)0.972578660.05003512y x =+ 10分误差平方和为 ................................ ()()4220.1502kkk y x y σ==-=∑ 11分2(21)0.972578660.050035122123.0380683y =+⨯= ....... 12分三(本题12分)（1）Gauss 求积公式(()11585()d 0.600.6999f x x f f f -≈-++⎰的代数精度是多少？（2）利用上述公式计算定积分111d 3x x-+⎰，并与真值进行比较；0.60.7745967≈）（3）将区间[1,1]-等分成两个小区间，对每个小区间分别利用（1）中的Gauss 求积公式进行计算，并将最后得到的111d 3x x-+⎰的近似值与（2）进行比较。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
4分
y( x n 1 ) y( x n ) y ( x n )h
5分
y n hfn
h 2 f n f n ( fn ) 2 x y
2 fn 2 f n 2 f n f n f n 2 h3 2 f n [ 2 2 fn f ( ) f h ] O(h 4 ) n 2 6 x xy x y y y
班
级
学
号
姓
名
… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
东 北 大 学 研 究 生 院 考 试 试 卷 2010 —2011 学年第 数值分析 一 学期
总分
1--3
4--6
7--9
10--12
13-15
3 15 3 15 ， x2 5 5 5 5
9． （7 分）给定离散数据 xi yi －2 -1 -1 1 0 2 1 3
Gauss 点为： x1
4分
1 1 15 15 积分公式为： x 2 f ( x)dx [ f ( ) f( )] 1 3 5 5
6分
试求形如 y a bx2 的拟合曲线。 解 由于基函数为： 0 ( x) 1,1 ( x) x 2 1分 3分 4分 12.(9 分)利用复化 Simpson 公式 S 2 计算定积分 I cos xdx 的近似值，并估计误
1
所以， xk 1 3 2 xk 3 , k 0,1,2,...对任意初值 x0 [1, 2] 都收敛。 3分
3 2 6.(6 分)设 xk 1 2xk axk bxk c, k 0,1,2,... 是求方程根 1 的迭代法，试
确定参数 a, b, c 使迭代法的收敛阶尽可能高，并指出阶是多少？ 5分 7分 解 由 1得1 2 a b c 令 (1) 6 2a b 0, (1) 12 2a 0 得 a 6, b 6, c 1 ， 此时，迭代法 3 阶收敛。 2分 5分 6分
2
n b 1 2 所以， Ak xk x 2 dx (b3 a 3 ) a 3 k 0
2分 4分 6分 7分
11.（6 分）对积分 x 2 f ( x)dx 建立两点 Gauss 公式。
1
1
则有： a b c 2,8a 4b 2c 1,3a 2b c 0 解得： a 1 / 2, c 7 / 2, b 1
1 1 7 所以， H 3 ( x) x( x 2 2 x 7) x 3 x 2 x 2 2 2
解
由于 P0 ( x) 1, P1 ( x) x
( P0 , x) P0 x ， ( P0 , P0 )
2分
P2 ( x) x 2
( P0 , x 2 ) (P , x 2 ) 3 P0 1 P1 x 2 ( P0 , P0 ) ( P1 , P1 ) 5
y n1 y n 2 f n f n h3 2 f ( f n ) ( 2n 2 fn f n ) O(h 4 ) 2 2 x y 2 x xy y
h2 h3 y ( x n ) y ( x n ) O(h 4 ) 2 6
10 ． （ 5 分）设求积公式 f ( x)dx Ak f ( xk ) , (n 2) 是插值型求积公式，求
b a k 0
n
A x
k 0 k
n
2 k
. 由于插值型求积公式代数精度至少是 n， 2分 5分
解 8． （7 分）求满足条件 f (0) 0, f (1) 2, f (2) 1, f (1) 0 的三次插值多项 式 H 3 ( x) 的表达式。 解 令 H 3 ( x) x(ax bx c)
… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
7． （6 分）设 f ( x) 3x 3 2x 5 ，求差商 f [0,1], f [1,2,3,4], f [1,2,3,4,5] 。 解 f[0,1]=(6-5)/1=1 f[1,2,3,4]=3 f[1,2,3,4,5]=0 2分 4分 6分
x1 x 2 2 x3 2 3.(6 分)解线性方程组 2 x1 3 x 2 3 x3 3 的 Gauss-Seidel 迭代法是否收敛，为 4 x 6 x 7 x 1 1 2 3
又由于 1 3 5 ( x) 3 2x 3 3 7 2 ，
(10 )
2/3
4分 8分
x
*
B
10
1 B
x (1) x ( 0)
211 0.034683059 310
5.(10 分)说明方程 x 3 2 x 3 0 在区间[1, 2]内有唯一根，并建立一个收敛的 迭代格式，使对任意初值 x0 [1, 2] 都收敛，说明收敛理由。 解 由 f (1) 4 0, f (2) 1 0, f ( x) 3x 2 2 0 ，知有唯一根。 3分 7分 9分
10 步的误差 x (10 ) x *
。
解
1/ 4 1/ 4 0 0 1 / 3 ,所以 B 2 / 3 , 由于 B 1 / 3 1/ 3 1/ 6 0
2分
又由于 x1 (1 / 4,2 / 3,1 / 2)T , 所以 x (1) x ( 0) 3分 所以， x 6分
y x sin xy 13．(5 分)求解初值问题 y(0) 2
为什么？ 解
0 x2
的改进 Euler 方法是否收敛？ 15.（5 分）证明矩阵谱半径 ( A) 不是矩阵范数。 2分 5分 证明
由于 f ( x, y) x sin xy 关于 y 满足 Lipschitz 条件，
2 2/3 1 3(2 x 3) 2 / 3
什么？
1 2 1 2 3 3 0 3 2 1 (3 2)(7 6) 解 令 0 2 4 6 7 0 0 7 6
得 (G) 6 / 7 1 , 所以，Gauss-Seidel 迭代法收敛。
7分 9分
3
于是， y( xn1 ) yn1 O(h3 ) ，此差分公式是 2 阶的。
解得： a 5 / 2, b 5 / 6 ， 拟合曲线为: y
5 5 2 x 2 6
＝0.909622804
M 4 max f
( 4)
( x) max cos x 1
7分
I S 2 | R( f ) |
25 0.000694444 2880 2 4
2
… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
所以，x 具有 4 位有效数字。
2 6 4 3 的 Crout 分解式 A TM . 2.(6 分)写出矩阵 A 3 6 2 1 5 2 2 6 4 2 3 3 3 3 3 解 由于 A 3 6 2 1 5 2 7 12 0 1 3 2 2 0 所以， A 3 3 0 0 1 3 2 7 12 0 0 1
课程名称：
(共 3 页）
4 x1 x 2 x3 1 0 (0) 4.(8 分)用 Jacobi 法解线性方程组 x1 3 x 2 x3 2 ，取 x 0 ，估计迭代 0 2 x x 6 x 3 2 3 1
1.(5 分)设近似值 x 12.25 近似 x * 的相对误差限为 0.00013， 问 x 具有几位有效数 字。 解 由于绝对误差限为:12.250.00013＝0.0015925<0.510-2 2分 5分
的差分公式：
h y n 1 y n 4 (3k1 k 2 ) k f (x , y ) n n 1 k 2 f ( x n 2h, y n 2hk1 ) y0
求此差分公式的阶。 解 由于 2分
2 f n f n 2 fn 2 fn 2 2 fn k 2 f n 2h( f n ) 2h ( 2 2 fn f n ) O(h 3 ) 2 x y xy x y
所以，改进 Euler 法收敛。
0 1 因为 ( A) 0 时，不一定有 A 0 , 例如 A 0 0 ,
2分 5分
所以， ( A) 不满足范数的非负性，不是范数。 14.（9 分）已知求解常微分方程初值问题：
y f ( x, y) , x [a, b] y ( a)
0 2
差。 解
1 1 3 I S 2 [cos 0 cos 2 2 cos 1 4 cos 4 cos ] 6 2 2
于是： 0 (1,1,1,1),1 (4,1,0,1), f (1,1,2,3)
3分 5分 6分 9分
4a 6b 5 正则方程组为： 6a 18b 0