研究生数值分析期末考试试题A答案

合集下载

上海海事大学1213数值分析试A卷答案

上海海事大学1213数值分析试A卷答案

上海海事大学2012---2013学年第 2 学期 研究生 数值分析 课程考试试卷A (答案)学生姓名: 学号: 专业:1. 利用Seidel 迭代法求解Ax=b 时,其迭代矩阵是))-1s U L D B -=(; 当系数矩阵A 满足 严格对角占优 时,Seidel 迭代法收敛 。

7. 反幂法是求可逆矩阵按模最小 特征值和特征向量的计算方法. 6. QR 法是计算 非奇异矩阵的 所有 特征值和特征向量的计算方法 1. 利用Jacobi 迭代法求解Ax=b 时,其迭代矩阵是)(1U L D B J +=-;当系数矩阵A 满足 严格对角占优 时,Jacobi 迭代法收敛 。

2. 对于求解Ax=b ,如果右端有b δ的扰动存在而引起解的误差为x δ,则相对误差≤xxδ bbA Cond δ)(3. 幂法是求矩阵 按模最大 特征值和特征向量的计算方法.Jacobi 法是计算 实对称矩阵的所有 特征值和特征向量的计算方法 六.设方程组Ax=b 有唯一解*x ,其等价变形构造的迭代格式为f Bx x k k +=+)()1(,如矩阵谱半径1)(>B ρ,但B 有一个特征值满足1<λ,求证:存在初始向量)0(x ,使得迭代产生的序列{})(x x 收敛于*x 。

(7分)证明: 由f Bx x k k +=+)()1(,f Bx x +=**()()*)0(1k *)(*)1(---x x B x x B x xk k ++== 对于B 的一个特征值满足1<λ,特征向量设为y ,,,11y y B y By k k ++==λλ故取初始向量y x x +=*)0(,有()y y B x x B x x k k 11k *)0(1k *)1(--++++===λ∞→→==+++k yy x xk k k ,0-11*)1(λλ,所以{})(x x 收敛于*x八.给定函数函数)(x f ,对于一切x ,存在)(x f ',且M x f m ≤'≤<)(0, 证明对于范围M20<<λ内的任意定数λ,迭代过程)(-1k k k x f x x λ=+均收敛于0)(=x f 的根。

数值分析期末考试和答案

数值分析期末考试和答案

数值分析期末考试和答案一、单项选择题(每题2分,共20分)1. 在数值分析中,下列哪个方法用于求解线性方程组?A. 插值法B. 迭代法C. 直接法D. 拟合法答案:C2. 以下哪个数值方法是用于求解非线性方程的?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 拉格朗日插值法答案:B3. 在数值积分中,梯形法则的误差与下列哪个因素无关?A. 被积函数的二阶导数B. 积分区间的长度C. 积分区间的划分数量D. 被积函数的一阶导数答案:D4. 以下哪个数值方法是用于求解常微分方程的?A. 欧拉方法B. 牛顿迭代法C. 拉格朗日插值法D. 高斯消元法答案:A5. 在数值分析中,下列哪个方法用于求解特征值问题?A. 高斯消元法B. 幂迭代法C. 牛顿迭代法D. 梯形法则答案:B6. 以下哪个数值方法是用于求解线性最小二乘问题的?A. 高斯消元法B. 梯形法则C. 正交分解法D. 牛顿迭代法答案:C7. 在数值分析中,下列哪个方法用于求解非线性方程组?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 欧拉方法答案:B8. 在数值分析中,下列哪个方法用于求解偏微分方程?A. 有限差分法B. 牛顿迭代法C. 线性插值法D. 梯形法则答案:A9. 在数值分析中,下列哪个方法用于求解优化问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 单纯形法答案:D10. 在数值分析中,下列哪个方法用于求解插值问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 拉格朗日插值法答案:D二、填空题(每题2分,共20分)1. 在数值分析中,求解线性方程组的直接法包括______消元法和______消元法。

答案:高斯;LU2. 牛顿迭代法的收敛速度是______阶的。

答案:二3. 梯形法则的误差与被积函数的______阶导数有关。

答案:二4. 欧拉方法是一种求解______阶常微分方程的数值方法。

答案:一5. 幂迭代法是求解______特征值问题的数值方法。

《数值分析》A卷期末考试试题及参考答案

《数值分析》A卷期末考试试题及参考答案

一、单项选择题(每小题3分,共15分) 1、用Simpson 公式求积分1401x dx +⎰的近似值为 ( ).A.2924 B.2429C.65D. 562、已知(1)0.401f =,且用梯形公式计算积分2()f x dx ⎰的近似值10.864T =,若将区间[0,2]二等分,则用递推公式计算近似值2T 等于( ). A.0.824 B.0.401 C.0.864 D. 0.8333、设3()32=+f x x ,则差商0123[,,,]f x x x x 等于( ).A.0B.9C.3D. 64的近似值的绝对误差小于0.01%,要取多少位有效数字( ). A.3 B.4 C.5 D. 25、用二分法求方程()0=f x 在区间[1,2]上的一个实根,若要求准确到小数 点后第四位,则至少二分区间多少次( ).A.12B.13C.14D. 15二、填空题(每小题4分,共40分)1、对于迭代函数2()=(3)ϕ+-x x a x ,要使迭代公式1=()ϕ+k k x x则a 的取值范围为 .2、假设按四舍五入的近似值为2.312,则该近似值的绝对误差限为 .3、迭代公式212(3)=,03++>+k k k k x x a x a x a收敛于α= (0)α>. 4、解方程4()530f x x x =+-=的牛顿迭代公式为 . 5、设()f x 在[1,1]-上具有2阶连续导数,[1,1]x ∀∈-,有1()2f x ''≤,则()f x 在[1,1]-上的线性插值函数1()L x 在点0处的误差限1(0)R ≤______.6、求解微分方程初值问题2(0)1'=-⎧⎨=⎩y xy yy ,0x 1≤≤的向前Euler 格式为 .7、设310131013A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭,则A ∞= .8、用梯形公式计算积分112-⎰dx x 的近似值为 . 9、设12A 21+⎡⎤=⎢⎥⎣⎦a 可作Cholesky 分解,则a 的取值范围为 . 10、设(0)1,(0.5) 1.5,(1)2,(1.5) 2.5,(2) 3.4f f f f f =====,若1=h ,则用三点公式计算(1)'≈f .三、解答题(共45分) 1、给定数据用复化Simpson 公式计算1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛. (5分)4、已知数据试对数据用最小二乘法求出形如=+y x b的拟合曲线. (8分) 5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (8分) 6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦一、单项选择题(每小题3分,合计15分) 1、A 2、D 3、C 4、C 5、D 二、填空题(每小题3分,合计30分) 1、0<<a ; 2、31102-⨯; 3;4、4135345++-=-+k k k k k x x x x x ; 5、14; 6、1(2)+=+-n n n n n y y h x y y ; 7、5;8、34-; 9、3>a ;10、1.2;三、计算题(合计55分) 1、给定数据用复化Simpson 公式计算 1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)解: 401024S [()4()()]6-=++x x f x f x f x ………… 1分 1.38 1.30(3.624 4.20 5.19)6-=+⨯+ 0.341= ………… 2分20422012234S [()4()()][()4()()]66--=+++++x x x xf x f x f x f x f x f x =0.342 ………… 6分2211[]15-≈-I S S S =-⨯40.6710 ………… 8分 2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 解:设111213212223313233u u u 123100135l 100u u 136l l 100u ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=*⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦………… 1分 111=u ,212=u ,313=u ,121=l ,131=l 122=u ,223=u ,132=l133=u ,133=l …………6分所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111011001L ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100210321U …………7分 由b Ly =得Ty )1,1,2(=;由y Ux =得Tx )1,1,1(-=. ………… 8分3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛.(6分)解:要使迭代序列具有平方收敛,则()0ϕ'*=x ………… 2分 而()()()ϕλ=+f x x x x ,即 ………… 3分 2()()()()10()λλλ''**-**+=*f x x x f x x …………4分 而()0*=f x 则有()1()λ'*=-*f x x ………… 5分所以()()23λ'=-=--x f x x ………… 6分4、已知数据试对数据用最小二乘法求出形如=+ay x b的拟合曲线. (8分) 解:因为11=+b x y a a ,令0111,,,====b a a y x x a a y……2分 则有法方程01461061410⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭a a ……5分解出014,1==-a a ,则1,4=-=-a b ……7分 所以1=4-y x……8分5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (7分)解:01()(2)8l x x x =- …………2分 211()(4)4l x x =-- …………4分21()(2)8l x x x =+ …………6分 2012()()(2)()(0)()(2)L x l x f l x f l x f =-++24=+x …………7分6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦解:100010001D ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,00010021002L ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦,10021002000U ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………3分1100211()0221002J B D L U -⎡⎤⎢⎥⎢⎥⎢⎥=+=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………5分 2102111()0222102J E B λλλλλλ⎡⎤-⎢⎥⎢⎥⎢⎥-=--=-=⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦…………6分()2J B ρ=…………7分 所以用Jacobi 迭代法求解方程组Ax b =收敛 …………8分。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案试题一:1. 简答题(共10分)a) 什么是数值分析?它的主要应用领域是什么?b) 请简要解释迭代法和直接法在数值计算中的区别。

2. 填空题(共10分)a) 欧拉方法是一种______型的数值解法。

b) 二分法是一种______法则。

c) 梯形法则是一种______型的数值积分方法。

3. 计算题(共80分)将以下函数进行数值求解:a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。

b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。

c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。

试题二:1. 简答题(共10分)a) 请解释什么是舍入误差,并描述它在数值计算中的影响。

b) 请解释牛顿插值多项式的概念及其应用。

2. 填空题(共10分)a) 数值稳定性通过______号检查。

b) 龙格-库塔法是一种______计算方法。

c) 零点的迭代法在本质上是将方程______转化为______方程。

3. 计算题(共80分)使用牛顿插值多项式进行以下计算:a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4),求在 x = 0.5 处的插值多项式值。

b) 已知插值节点 (0, 1), (1, 2), (3, 7),求插值多项式,并计算在 x = 2 处的值。

c) 使用 4 阶龙格-库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。

答案:试题一:1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。

它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。

b) 迭代法和直接法是数值计算中常用的两种方法。

迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。

数值分析期末复习题答案

数值分析期末复习题答案

数值分析期末复习题答案一、选择题1. 以下哪个算法是用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 共轭梯度法D. 辛普森积分法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的主要区别是什么?A. 插值点的选取不同B. 插值多项式的构造方式不同C. 计算复杂度不同D. 适用的函数类型不同答案:B3. 在数值积分中,梯形法则和辛普森法则的主要区别是什么?A. 精度不同B. 适用的积分区间不同C. 计算方法不同D. 稳定性不同答案:A二、简答题1. 解释什么是数值稳定性,并举例说明。

答案:数值稳定性指的是数值方法在计算过程中对于舍入误差的敏感程度。

例如,在求解线性方程组时,如果系数矩阵的条件数很大,则该方程组的数值解对舍入误差非常敏感,即数值稳定性差。

2. 说明数值微分与数值积分的区别。

答案:数值微分是估计函数在某一点的导数,而数值积分是估计函数在某个区间上的积分。

数值微分通常用于求解函数的局部变化率,而数值积分用于求解函数在一定区间内的累积效果。

三、计算题1. 给定一组数据点:(1, 2), (2, 3), (3, 5), (4, 6),请使用拉格朗日插值法构造一个三次插值多项式。

答案:首先写出拉格朗日插值基函数,然后根据数据点构造插值多项式。

具体计算过程略。

2. 给定函数 f(x) = x^2,使用牛顿-科特斯公式中的辛普森积分法在区间 [0, 1] 上估计积分值。

答案:首先确定区间划分,然后应用辛普森积分公式进行计算。

具体计算过程略。

四、论述题1. 论述数值分析中误差的来源及其控制方法。

答案:误差主要来源于舍入误差和截断误差。

舍入误差是由于计算机在进行浮点数运算时的精度限制造成的,而截断误差是由于数值方法的近似性质导致的。

控制误差的方法包括使用高精度的数据类型、选择合适的数值方法、增加计算步骤等。

五、综合应用题1. 给定一个线性方程组 Ax = b,其中 A 是一个 3x3 的矩阵,b 是一个列向量。

数值分析试题(卷)和答案解析

数值分析试题(卷)和答案解析

【试题__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 3i y2 4 12 <3i y '并估计误差。

(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 λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。

复旦大学《数值分析》2020-2021学年第一学期期末试卷A卷

复旦大学《数值分析》2020-2021学年第一学期期末试卷A卷

课程编号:A071001复旦大学2020-2021学年第一学期数值分析期末试题A一.解下列各题(每小题6分)1.求极限n n nn )111(lim 2++∞→.2..已知f 是可导函数,且x x f dx d 11(arctan =,求4(πf '.微分法,可以补用考虑微分次数,不断向下推。

导数法,比需两边对同一变量求导。

3.求出23||ln )(2+-=x x x x f 的间断点,并指出是第几类间断点.4.已知2)13(lim 2=++-+∞→bx ax x x ,试确定其中常数b a ,.二.解下列各题(每小题7分)1.设⎩⎨⎧+=+-=23)1ln(tt y t t x ,求22dx y d .2.试确定常数b a ,的值,使点)3,1(是曲线34bx ax y +=的拐点,并求出曲线的凹凸区间.3.求由方程0sin 21=+-y y x 所确定的隐函数)(x y y =的二阶导数.4.已知2112sin )(1lim30=--+→x x e x x f ,求)(lim 0x f x →.复合函数与函数求导公式可以一起用。

三.(9分)设数列}{n x 满足010<<-x ,),2,1,0(221 =+=+n x x x n nn ,证明}{n x 收敛,并求n n x ∞→lim .四.(9分)设)(x f 有二阶连续导数,0)0(=f ,⎪⎩⎪⎨⎧='≠=0),0(0,)()(x f x x x f x g ,求)(x g '并讨论)(x g '的连续性.五.(9分)一个体积给定的观察站底部是一个直圆柱,顶部是一个半球形,如果顶部单位面积的造价是侧面单位面积造价的二倍,问圆柱的底半径r 与高h 分别为多少时可使总造价最低?六.(8分)证明,当1>x 时,11ln +-≥x x x .七.(9分)(1)已知当0→x 时,2cos x e x -与k cx 是等价无穷小,求c 与k 的值;(2)求极限222sin )(cos 112lim 2xe x x x x x -+-+→.八.(4分)设)(xf 在],[b a 上连续,在),(b a 内可导,0)(≠'x f ,证明存在),(,b a ∈ηξ,使ηηξ---=''e ab e e f f a b )()(.最后一道题一定要会拼与凑。

060708研究生数值分析试卷(A).doc

060708研究生数值分析试卷(A).doc

武汉大学2006〜2007学年第一学期硕士研究生期末考试试题(A 卷)科H 名称:数值分析 学生所在院: 学号: 姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。

一、(12分)设方程组Ax = 0为■1、 (1\J 1>(1)用Doolittle 分解法求解方程组;(2) 求矩阵A 的条件数Cwd(A)g 二、(12分)设A 为n 阶对称正定矩阵,A的n 个特征值为山 < 心< .•. V 九,为 求解方程组Ax = b,建立迭代格式求出常数s 的取 值范围,使迭代格式收敛。

三、(12分)已知数据试用二次多项式p ⑴=ax 1 2+hx + c 拟合这些数据。

四、(14分)已知y = /(x)的数据如下:取得最小值。

六、 (12)确定常数片,使求积公式1求f (x)的Hermite 插值多项式W 3(x);2 为求\\f{x)dx 的值,采用算法:•⑴必:=「久3)击+ R 试导出截断误差R五、(12分)确定常数。

,b 的值,使积分r I.2I(a,b) = J 0(czx + /?-/) dxc 2^f{x)dx a A/(0) + A2/(l) + A3/(2)的代数精度尽可能高,并问是否是Gauss型公式。

七、(12分)设伊⑴导数连续,迭代格式x M =(p{x k)—阶局部收敛到点x*。

对于常数人,构造新的迭代格式:A 1 ,、队=一从+ 一心)1 +2 1 + 人问如何选取人,使新迭代格式有更高的收敛阶,并问是儿阶收敛。

八、(14分)对于下面求解常微分方程初值问题」方= 的单步法:Mo) = JoA)'〃+】=儿 + hk2< k、=(1)验证它是二阶方法;(2)确定此单步法的绝对稳定区域。

武汉大学2007~2008学年第一学期硕士研究生期末考试试题科目名称:数值分析学生所在院:学号:姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2010年秋研究生数值分析期末考试试题答案
一、单选题(4*5=20分)1、D; 2、B ; 3、D ; 4、B ; 5、D 。

二、填空题(4*5=20)1、4; 2、
⎪⎪⎪⎭⎫ ⎝⎛323
203*⎪⎪⎪⎪


⎝⎛320
323; 3、)]23
()0()23([3f f f ++-∏;4、k
k k k x x x x 2221--=+;5、9.605。

三、(10分)由两点三次Hermite 插值多项式公式秋得:
)2()(23x x x H -=,设所求多项式223)1()()(-+=x Ax x H x P ,。

(4分) 由P(2)=1,得A=1/4,。

(4分) 故22
)3(4
1)(-=
x x x P 。

.。

(2分) 四、(10分)设⎪⎪⎪


⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--=10
01
001*10010021321
u u l l l A ,由追赶法公式求得, 15/56,15/4,4/15,4/1,432211=-==-==l u l u l ,。

(4分) 由Ly=d,求得T y )77.0,87.0,25.0(=,(3分) 由Ux=y,求得,T x )5179.0,0714.1,7679.0(=(3分)
五、(10分)Jacobi 迭代计算格式:⎪⎩⎪⎨⎧++-=--=--=+++3/)221(5/)327(24)(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 x x x 。

(2分) G-S 迭代计算格式: ⎪⎩⎪⎨⎧++-=--=--=++++++3/)221(5/)327(24)
1(2)1(1)1(3
)
(3)1(1)
1(2)
(3)(2)1(1k k k k k k k k k x x x x x x x x x 。

(2分) 由于016415)(3=-+=-λλλJ B I del ,,115
16
)(>=
J B ρ即Jacobi 迭代发散;。

(2分) 由于0)1215()(2
=-=-λλλG B I del ,,15
2)(<=
G B ρ即G-S 迭代收敛;。

(2分)
SOR 迭代格式:⎪⎩⎪⎨⎧++-+-=--+-=--+-=++++++3/)221()1(5/)327()1()24()1()
1(2)1(1)(3)1(3
)
(3)1(1)(2)1(2)
(3)(2)(1)1(1k k k k k k k k k k k k x x w x w x x x w x w x x x w x w x 由于2
)(112J opt B w ρ-+=,而
1)(>J B ρ,所以最佳松弛因子w 不存在。

(2分)
六、(10分)设2
1
)(,1)(x x f x a x f =
'-=,于是Newton 迭代为:.。

(3分) )2(11
1k k k
k
k k ax x x x a x x -=-
-
=+,。

(4分) 由于a
x 1
*
=
是f(x)=0的单根,故为平方收敛。

(3分) 七、(10分)取02/)35()(33=-=x x x P 的零点,5
3
,0,53210=
=-
=x x x 。

(3分) 构造公式:

-++-
=1
1
210)5
3
()0()53()(f A f A f A dx x f ,具有5次代数精度。

(2分) 令2,,1)(x x x f =,可得:⎪⎪⎪

⎪⎪
⎪⎨⎧
=
++=++-=++3253*0*53*053*0*)53(*22102
10210A A A A A A A A A ,。

(3分)
求得:9
5
,98,95210===
A A A 。

(2分) 八、(10
分)改进的Euler 公式计算格式为:⎪⎪⎪

⎪⎪⎪⎨⎧+=++-+=++-+=+++)(21]100[*]100[*113
13
c p n n n p n c n n
n n p y y y x x y h y y x x y h y y 。

(3分) 4阶RK 公式的计算格式为:⎪⎪⎪⎪

⎩⎪
⎪⎪⎪⎪⎨⎧
+++++-=+++++-=+++++-=++-=++++=+)
()()(100)2
()2()2(100)2()2()2(100100)22(6334323
3123
143211h x h x hK y K h x h x K h y K h x h x K h y K x x y K K K K K h y y n n n n n n n
n n n
n n n n 。

(3分)
依题,100-=λ由于Euler 方法稳定域为:02<<=-h λ; 因此当步长满足:02.00<=<h 时计算稳定。

(2分) 而4阶RK 方法稳定域为:0785.2<<=-h λ, 因此当步长满足:02785.00<=<h 时计算稳定。

(2分)。

相关文档
最新文档