数值分析、计算方法试题库及答案

一、单项选择题（每小题3分,共15分）1. 和分别作为π(de)近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和42. 已知求积公式()()211211()(2)636f x dx f Af f ≈++⎰,则A ＝（ ）A ． 16B ．13C ．12D ．233. 通过点()()0011,,,x y x y (de)拉格朗日插值基函数()()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 =(de)根(de)牛顿法收敛,则它具有（ ）敛速.A ．超线性B ．平方C ．线性D ．三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到(de)第3个方程（ ）.A ．232x x -+=B ．232 1.5 3.5x x -+=C ．2323x x -+=D ．230.5 1.5x x -=-二、填空题（每小题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 =,用欧拉法解初值问题()211y y yx y ⎧'=+⎪⎨⎪=⎩(de)计算公式 .0,1,2分 人三、计算题（每题15分,共60分）1. 已知函数211y x =+(de)一组数据：求分段线性插值函数,并计算()1.5f (de)近似值.1. 解 []0,1x ∈,()1010.510.50110x x L x x --=⨯+⨯=---[]1,2x ∈,()210.50.20.30.81221x x L x x --=⨯+⨯=-+--所以分段线性插值函数为()[][]10.50,10.80.31,2x x L x x x ⎧-∈⎪=⎨-∈⎪⎩ ()1.50.80.3 1.50.35L =-⨯=2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩（1） 写出雅可比迭代公式、高斯－塞德尔迭代公式；（2） 对于初始值()()00,0,0X =,应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算()1X （保留小数点后五位数字）.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之间(de)近似根（1）请指出为什么初值应取2 （2）请用牛顿法求出近似根,精确到. 3. 解()331f x x x =--,()130f =-<,()210f =>()233f x x '=-,()12f x x ''=,()2240f =>,故取2x =作初始值4. 写出梯形公式和辛卜生公式,并用来分别计算积分111dxx+⎰.四、证明题（本题10分）确定下列求积公式中(de)待定系数,并证明确定后(de)求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明：求积公式中含有三个待定系数,即101,,A A A -,将()21,,f x x x =分别代入求一、 填空（共20分,每题2分）1. 设2.3149541...x *=,取5位有效数字,则所得(de)近似值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 .4．求方程 21.250x x --= (de)近似根,用迭代公式 1.25x x =+,取初始值 01x =, 那么 1______x =。

一. 填空题（本大题共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 λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。

（10分）《数值分析》（A ）卷标准答案（2009－2010－1）一． 填空题（每小题3分，共12分） 1. ()1200102()()()()x x x x l x x x x x --=--； 2.7；3. 3，8；4. 2n+1。

数值计算方法试题一一、 填空题（每空1分，共17分）1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。

2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。

3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x 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 插值基函数，则∑==nk kx l0)(( )，∑==nk k jk x lx 0)(( )，当2≥n 时=++∑=)()3(204x l x xk k n k k( )。

5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f。

6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。

7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ϕ，则⎰=14)(dx x x ϕ 。

8、给定方程组⎩⎨⎧=+-=-221121b 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 满足（ ）条件时，这种分解是唯一的。

数值分析试卷及答案数值分析试卷一、选择题（共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。

数值计算方法试题一一、 填空题（每空1分，共17分）1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。

2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。

3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x 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 插值基函数，则∑==nk kx l0)(( )，∑==nk k jk x lx 0)(( )，当2≥n 时=++∑=)()3(204x l x xk k n k k( )。

5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f。

6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。

7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ϕ，则⎰=14)(dx x x ϕ 。

8、给定方程组⎩⎨⎧=+-=-221121b 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 满足（ ）条件时，这种分解是唯一的。

数值分析复习试题第一章 绪论 一. 填空题 1.*x为精确值x 的近似值；()**x f y =为一元函数()x f y =1的近似值；()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-：***r x xe 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≈-211.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 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 . 8、精确值 14159265.3*=π，则近似值141.3*1=π和1415.3*2=π分别有 3位和 4 位有效数字。

9、若*2.71828x e x =≈=,则x 有 6 位有效数字，其绝对误差限为1/2*10-5。

10、 设x*的相对误差为2％，求(x*)n的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；12、计算方法主要研究( 截断 )误差和( 舍入 )误差； 13、为了使计算 ()()2334610111y x x x =++---- 的乘除法次数尽量地少，应将该表达式改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+。

数值计算方法试题一一、 填空题（每空1分，共17分）1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。

2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。

3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x 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 插值基函数，则∑==nk kx l0)(( )，∑==nk k jk x lx 0)(()，当2≥n 时=++∑=)()3(204x l x xk k nk k ( )。

5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f。

6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。

7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ϕ，则⎰=14)(dx x x ϕ 。

8、给定方程组⎩⎨⎧=+-=-221121b 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 满足（ ）条件时，这种分解是唯一的。

数值分析复习试题第一章绪论一.填空题1.为精确值的近似值；为一元函数的近似值；*xx ()**x f y =()x f y =1为二元函数的近似值，请写出下面的公式：：()**,*y x f y =()y x f y ,2=**e x x =-***r x xe 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≈-211.73 10 2-≤⨯4、设均具有3位有效数字，则的相对误差限为 0.0055 。

121.216, 3.654x x ==12x x 5、设均具有3位有效数字，则的误差限为 0.01 。

121.216, 3.654x x ==12x x +6、已知近似值是由真值经四舍五入得到,则相对误差限为0.0000204 .2.4560A x =T x 7、递推公式如果取作计算,则计算到时,误差为,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,0 1.41y =≈10y ;这个计算公式数值稳定不稳定 不稳定 .8110 2⨯8、精确值，则近似值和分别有 3 位和14159265.3*=π141.3*1=π1415.3*2=π4 位有效数字。

9、若,则x 有 6 位有效数字，其绝对误差限为1/2*10-5 。

*2.71828x e x =≈=10、 设x*的相对误差为2％，求(x*)n 的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；12、计算方法主要研究( 截断 )误差和( 舍入 )误差；13、为了使计算 的乘除法次数尽量地少，应将该表达式()()2334610111y x x x =++----改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+。

目录一、绪论------------------------------------------------------------------------------------- 2-2二、线性方程组直接解法列主元高斯LU LDL T GG T-------------------- 3-6二、线性方程组迭代法----------------------------------------------------------------- 7-10 三、四、非线性方程组数值解法二分法不动点迭代---------------------- 11-13五、非线性方程组数值解法牛顿迭代下山弦截法----------------- 14-15六、插值线性插值抛物线插值------------------------------------------------ 16-18七、插值Hermite插值分段线性插值-----------------------------------------19-22八、拟合------------------------------------------------------------------------------------ 23-24九、数值积分----------------------------------------------------------------------------- 25-29十、常微分方程数值解法梯形欧拉改进----------------------------------- 30-32 十一、常微分方程数值解法龙格库塔------------------------------------------ 33-35绪论1-1 下列各数都是经过四舍五入得到的近似值 ,试分别指出它们的绝对误差限,相对误差限和有效数字的位数.X 1 =5.420, X 2 =0.5420, X 3 =0.00542, X 4 =6000, X 5 =0.6×105注：将近似值改写为标准形式X 1 =（5*10-1+4*10-2+2*10-3+0*10-4）*101 即n=4,m=1 绝对误差限｜△X 1｜=｜X *1-X 1｜≤ 12×10m-n ＝12×10-3 相对误差限｜△r X 1｜= ｜X∗1−X1｜｜X∗1｜≤｜X∗1−X1｜｜X1｜= 12×10-3/5.4201-2 为了使101/2 的相对误差小于0.01%, 试问应取几位有效数字?1-3 求方程x 2 -56x+1=0的两个根, 使它们至少具有4位有效数字( √783≈27.982)注：原方程可改写为(x-28)2=783线性方程组解法（直接法）2-1用列主元Gauss消元法解方程组解：回代得解：X1=0 X2=-1 X3=12-2对矩阵A进行LU分解，并求解方程组Ax=b,其中解：(注：详细分解请看课本P25)A=(211132122)→(211(1/2)5/23/2(1/2)3/23/2)→(2111/25/23/21/2(3/5)3/5)即A=L×U=(11/211/23/51)×(2115/23/23/5)先用前代法解L y＝P b 其中P为单位阵（原因是A矩阵未进行行变换）即L y＝P b 等价为(11/211/23/51)(y1y2y3)＝(111)(465)解得 y 1=4 y 2=4 y 3=35再用回代解Ux ＝y ，得到结果x即Ux ＝y 等价为(2115/23/23/5)(x 1x 2x 3)＝(y 1y 2y 3)＝(443/5) 解得 x 1=1 x 2=1 x 3=1即方程组Ax=b 的解为x ＝(111)2-3 对矩阵A 进行LDL T 分解和GG T 分解，求解方程组Ax=b,其中A=(164845−48−422) ， b ＝(123)解：（注：课本 P 26 P 27 根平方法）设L=(l i j )，D=diag(d i )，对k=1,2,…,n,其中d k =a kk －∑l kj 2k−1j=1d jl ik =(a ik −∑l ij l kj k−1j=1d j )/ d k 即d 1=a 11－∑l 1j 20j=1d j =16－0=16因为 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=a 21/ d 1=416=14 所以d 2=a 22－∑l 2j 21j=1d j =5－(14)2d 1=4同理可得d 3=9 即得 D=(1649)同理l 11=(a 11−∑l ij l 1j 0j=1d j )/ d 1=1616=1=l 22=l 33 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=416=14 l 31=(a 31−∑l 3j l 1j 0j=1d j )/ d 1=816=12 l 32=(a 32−∑l 3j l 2j 1j=1d j )/ d 2=−4−12×14×164=−64=－32即L=(114112−321) L T=(114121−321) 即LDL T分解为A=(114112−321)(1649)(114121−321)解解：A=(164845−48−422)→(41212−32−33)故得GG T分解：A=(4122−33)(4122−33) LDL T分解为A=(114112−321)(1649)(114121−321) 由(114112−321)(y 1y 2y 3)＝(123) ，得(y 1y 2y 3)＝(0.250.8751.7083)再由(4122−33)(x 1x 2x 3)＝(0.250.8751.7083) ，得(x 1x 2x 3)＝(−0.54511.29160.5694)2-4 用追赶法求解方程组：解：(4−1−14−1−14−1−14−1−14)→(4−14−1154−415−15615−1556−120956−56209−1780209)由(4−1154−15615−120956−1780209)(y1y2y3y4y5)＝(100200)，得(y1y2y3y4y5)＝(256.66671.785700.4784753.718)再由(1−141−4151−15561−562091)(x1x2x3x4x5)＝(256.66671.785700.4784753.718)，得(x1x2x3x4x5)＝(27.0518.20525.769314.87253.718)线性方程组解法（迭代法）2-1 设线性方程组{4x 1−x 2+2x 3=1−x 1−5x 2+x 3=22x 1+x 2+6x 3=3（1） 写出Jacobi 法和SOR 法的迭代格式（分量形式） （2） 讨论这两种迭代法的收敛性（3） 取初值x (0)=(0，0，0)T ，若用Jacobi 迭代法计算时，预估误差 ｜｜x*-x (10)｜｜∞ （取三位有效数字）解：（1）Jacobi 法和SOR 法的迭代格式分别为Jacobi 法迭代格式SOR（2）因为A 是严格对角占优矩阵,但不是正定矩阵,故Jacobi 法收敛,SOR 法当0<ω≤1时收敛.⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=-+-=+-=+++216131525151412141)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x xx x ⎪⎪⎪⎩⎪⎪⎪⎨⎧-++-=+-+-=+-+-+=++++++)216131()525151()412141()(3)1(2)1(1)(3)1(3)(3)(2)1(1)(2)1(2)(3)(2)(1)(1)1(1k k k k k k k k k k k k k k k x x x x x x x x x x x x x x x ωωω（3）由(1)可见||B ||∞=3/4,且取x (0)=(0,0,0)T ,经计算可得x (1)=(1/4,-2/5,1/2)T ,于是||x (1)-x (0)||∞=1/2,所以有2-2 设方程组为{5x 1+2x 2+x 3=−12−x 1+4x 2+2x 3=202x 1−3x 2+10x 3=3试写出其Jacobi 分量迭代格式以及相应的迭代矩阵，并求解。

一． 单项选择题（每小题2分，共10分）1． 在下列四个数中，有一个数具有4位有效数字，且其绝对误差限为 51021-⨯，则该数是（ ） A 0.001523 B 0.15230 C 0.01523 D 1.52300 2． 设方阵A 可逆，且其n 个特征值满足：n λλλ>≥> (21)，则1-A 的主特征值是（ ）A11λ B nλ1 C1λ或n λ D 11λ或nλ13． 设有迭代公式→→+→+=fxB x k k )()1(。

若||B|| > 1，则该迭代公式（ ）A 必收敛B 必发散C 可能收敛也可能发散4． 常微分方程的数值方法，求出的结果是（ ）A 解函数B 近似解函数C 解函数值D 近似解函数值 5． 反幂法中构造向量序列时，要用到解线性方程组的（ ） A 追赶法 B LU 分解法C 雅可比迭代法D 高斯—塞德尔迭代法二． 填空题（每小题4分，共20分）1． 设有方程组⎪⎩⎪⎨⎧=+-=+-=+02132432132132x x x x x x x x ，则可构造高斯—塞德尔迭代公式为⎪⎩⎪⎨⎧2． 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=111112101A ，则=∞A3． 设1)0(,2'2=+=y y x y ，则相应的显尤拉公式为=+1n y4． 设1)(+=ax x f ,2)(x x g =。

若要使)(x f 与)(x g 在[0，1]上正交，则a =5． 设T x )1,2,2(--=→，若有平面旋转阵P ，使P →x 的第3个分量为0，则P =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡ 三． 计算题（每小题10分，共50分）1． 求27的近似值。

若要求相对误差小于0.1%，问近似值应取几位有效数字？2． 设42)(x x x f -=,若在[-1，0]上构造其二次最佳均方逼近多项式,请写出相应的法方程。

3． 设有方程组⎪⎩⎪⎨⎧=++=++=-+1221122321321321x x x x x x x x x ，考察用雅可比迭代解此方程组的收敛性。

数值计算方法试题一一、 填空题（每空1分，共17分）1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。

2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。

3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x 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 插值基函数，则∑==nk kx l0)(( )，∑==nk k jk x lx 0)(( )，当2≥n 时=++∑=)()3(204x l x xk k n k k( )。

5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f 。

6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。

7、{}∞=0)(k k x ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ϕ，则⎰=14)(dx x x ϕ 。

8、给定方程组⎩⎨⎧=+-=-221121b 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 满足（ ）条件时，这种分解是唯一的。

模 拟 试 卷（一）一、填空题（每小题3分，共30分）1．有3个不同节点的高斯求积公式的代数精度是 次的.2．设152210142-⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦A ，342⎛⎫ ⎪=- ⎪ ⎪⎝⎭x ，则 ∞A = .， 1x = ______.3．已知y =f (x )的均差（差商）01214[,,]3f x x x =，12315[,,] 3f x x x =，23491[,,]15f x x x =，0238[,,] 3f x x x =, 那么均差423[,,]f x x x = .4．已知n =4时Newton －Cotes 求积公式的系数分别是：,152,4516,907)4(2)4(1)4(0===C C C 则)4(3C ＝ .5．解初始值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进的Euler 方法是 阶方法；6．求解线性代数方程组123123123530.13260.722 3.51x x x x x x x x x --=⎧⎪-++=⎨⎪++=⎩的高斯—塞德尔迭代公式为 ，若取(0)(1,1,1)=-x, 则(1)=x .7．求方程()x f x =根的牛顿迭代格式是 . 8．1(), (),, ()n x x x 是以整数点01, ,, ,n x x x 为节点的Lagrange 插值基函数，则()n kjk k xx =∑= .9．解方程组=Ax b 的简单迭代格式(1)()k k +=+xBx g 收敛的充要条件是 .10．设(-1)1,(0)0,(1)1,(2)5f f f f ====，则()f x 的三次牛顿插值多项式为 ，其误差估计式为 .二、综合题（每题10分，共60分）1．求一次数不超过4次的多项式()p x 满足：(1)15p =，(1)20p '=，(1)30p ''=(2)57p =，(2)72p '=.2．构造代数精度最高的形式为10101()()(1)2xf x dx A f A f ≈+⎰的求积公式，并求出 其代数精度.3．用Newton 法求方程2ln =-x x 在区间) ,2(∞内的根, 要求8110--<-kk k x x x .4．用最小二乘法求形如2y a bx =+的经验公式拟合以下数据：5．用矩阵的直接三角分解法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡71735 30103421101002014321x x x x .6 试用数值积分法建立求解初值问题0(,)(0)y f x y y y '=⎧⎨=⎩的如下数值求解公式1111(4)3n n n n n hy y f f f +-+-=+++，其中(,),1,,1i i i f f x y i n n n ==-+.三、证明题（10分）设对任意的x ，函数()f x 的导数()f x '都存在且0()m f x M '<≤≤,对于满足20Mλ<<的任意λ，迭代格式1()k k k x x f x λ+=-均收敛于()0f x =的根*x .参考答案一、填空题1．5； 2. 8, 9 ; 3.9115； 4. 1645； 5. 二； 6. (1)()()123(1)(1)()213(1)(1)(1)312(330.1)/5(220.7)/6(12)*2/7k k k k k k k k k x x x x x x x x x ++++++⎧=++⎪=+-⎨⎪=--⎩, ，，0.1543)7. 1()1()k k k k k x f x x x f x +-=-'-; 8. j x ; 9. ()1B ρ<; 10.32(4)11,()(1)(1)(2)/24(1,2)66x x x f x x x x ξξ+-+--∈-二、综合题 1．差商表：233234()1520(1)15(1)7(1)(1)(2)5432p x x x x x x x x x x =+-+-+-+--=++++其他方法：设233()1520(1)15(1)7(1)(1)()p x x x x x ax b =+-+-+-+-+ 令(2)57p =，(2)72p '=，求出a 和b. 2．取()1,f x x =，令公式准确成立，得：0112A A +=,011123A A +=, 013A =, 116A =. 2()f x x =时，公式左右14=；3()f x x =时，公式左15=, 公式右524=∴ 公式的代数精度2=.3．此方程在区间) ,2(∞内只有一个根s ，而且在区间（2，4）内。

数值分析试题一、 填空题（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 )的特点是：=∑=ni 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 。

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 ）。

1. 和分别作为π的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．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 -=-单项选择题答案二、填空题（每小题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.答案()101x L x -=-()12x L x -=-()10.8L 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）请用牛顿法求出近似根，精确到.计算题3.答案4. 写出梯形公式和辛卜生公式，并用来分别计算积分1011dx x +⎰.计算题4.答案四、证明题（本题10分）确定下列求积公式中的待定系数，并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案证明：求积公式中含有三个待定系数，即101,,A A A -，将()21,,f x x x =分别代入求积公式，并令其左右相等，得一、 填空（共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 。

数值分析试题及答案一、单项选择题（每题3分，共30分）1. 线性代数中，矩阵A的逆矩阵记作（）。

A. A^TB. A^-1C. A^+D. A*答案：B2. 插值法中，拉格朗日插值多项式的基函数是（）。

A. 多项式B. 指数函数C. 正弦函数D. 余弦函数答案：A3. 在数值积分中，梯形规则的误差是（）阶的。

A. O(h^2)B. O(h^3)C. O(h)D. O(1/h)答案：A4. 求解线性方程组时，高斯消元法的基本操作不包括（）。

A. 行交换B. 行乘以非零常数C. 行加行D. 行除以非零常数答案：D5. 非线性方程f(x)=0的根的迭代法中，收敛的必要条件是（）。

A. f'(x)≠0B. f'(x)=0C. |f'(x)|<1D. |f'(x)|>1答案：C6. 利用牛顿法求解非线性方程的根时，需要计算（）。

A. 函数值B. 函数值和导数值C. 函数值和二阶导数值D. 函数值、一阶导数值和二阶导数值答案：B7. 矩阵的特征值和特征向量是（）问题中的重要概念。

A. 线性方程组B. 特征值问题C. 线性规划D. 非线性方程组答案：B8. 在数值分析中，条件数是衡量矩阵（）的量。

A. 稳定性B. 可逆性C. 正交性D. 稀疏性答案：A9. 利用龙格现象说明，高阶插值多项式在区间端点附近可能产生（）。

A. 振荡B. 收敛C. 稳定D. 单调答案：A10. 雅可比迭代法和高斯-塞德尔迭代法都是求解线性方程组的（）方法。

A. 直接B. 迭代C. 精确D. 近似答案：B二、填空题（每题4分，共20分）11. 线性代数中，矩阵A的行列式记作________。

答案：det(A) 或 |A|12. 插值法中，牛顿插值多项式的基函数是________。

答案：差商13. 在数值积分中，辛普森规则的误差是________阶的。

答案：O(h^4)14. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。

数值分析所有常考例题及详细答案第二章线性方程组的直接解法 (2)第三章解线性方程组的迭代法 (4)第五章非线性方程和方程组的数值解法 (7)第六章插值法与数值微分 (11)第七章数据拟合与函数逼近 (16)第八章数值积分 (20)第九章常微分方程的数值解法 (25)第二章 线性方程组的直接解法1、用LU 分解法求如下方程组的解（1）3351359059171⎛⎫⎛⎫ ⎪ ⎪X = ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭，（2）3235220330127X ⎡⎤⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦解：（1）13351124522133A L U ⎛⎫⎛⎫⎪⎪ ⎪⎪== ⎪⎪ ⎪⎪⎪⎪⎝⎭⎝⎭4(101)(1,1,)339(,,2)22T TTL Y Y UX Y X =⇒=-=⇒=-（2）132332352222012333301271313b ⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥-⎣⎦⎣⎦15521133371311y y ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⇒=-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥-⎣⎦⎣⎦ 3235121123321313X X ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-=-⇒=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦2121311()21()44254213142541425421310212127127350624r r r r r r +-↔+-⎡⎤⎢⎥-⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥-⎣⎦→→ 32344254102127210084r r +⎡⎤⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥-⎣⎦→得同解方程组1232334254121272184x x x x x x ⎧⎪++=⎪⎪-+=-⎨⎪⎪-=⎪⎩回代求解得(9,1,6)TX =--②212131112312323111011323231110523523111011323032323122112215747012323r r r r r r +↔+⎡⎤⎢⎥----⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥----⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎣⎦⎣⎦⎢⎥-⎣⎦→→323252()57231110231110574757470101232323235235193223030023235757r r r r +-↔⎡⎤⎡⎤⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦→→得同解方程组12323323110574701232319322300()5757x x x x x x ⎧⎪-++=⎪⎪++=-⎨⎪⎪++-=⎪⎩回代得(0.212435,0.549222, 1.15544)T X =-4、用Jordan 消去法解矩阵方程，AX B =其中：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=112221111A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=011001B 解：容易验证0A ≠，故A 可逆，有1X A B -= .因此，写出方程组的增广矩阵，对其进行初等变换得111101111011110122010111101111211100313000263---⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-→--→--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦⎣⎦100211002110111101022330013001322⎡⎤⎡⎤⎢⎥-⎢⎥-⎢⎥⎢⎥⎢⎥→--→-⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥-⎣⎦⎣⎦121122332X A B -⎡⎤⎢⎥-⎢⎥⎢⎥∴==-⎢⎥⎢⎥⎢⎥-⎣⎦5、用LU 分解法求解如下方程组12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦解：100256210037341004A LU -⎡⎤⎡⎤⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦12312311021193413010,19201,34304(10,1,4)TLy by y y y y y y =⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦==-=-=-==-(1)解得即 123321(2)25610371441,2,3(3,2,1)T Ux yx x x x x x x =-⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦====解解得：所以方程组的解为。