西南交大 数值分析题库
西南交通大学2018-2019数值分析Matlab上机实习题

西南交通⼤学2018-2019数值分析Matlab上机实习题数值分析2018-2019第1学期上机实习题f x,隔根第1题.给出⽜顿法求函数零点的程序。
调⽤条件:输⼊函数表达式()a b,输出结果:零点的值x和精度e,试取函数区间[,],⽤⽜顿法计算附近的根,判断相应的收敛速度,并给出数学解释。
1.1程序代码:f=input('输⼊函数表达式:y=','s');a=input('输⼊迭代初始值:a=');delta=input('输⼊截⽌误差:delta=');f=sym(f);f_=diff(f); %求导f=inline(f);f_=inline(f_);c0=a;c=c0-f(c0)/f_(c0);n=1;while abs(c-c0)>deltac0=c;c=c0-f(c0)/f_(c0);n=n+1;enderr=abs(c-c0);yc=f(c);disp(strcat('⽤⽜顿法求得零点为',num2str(c)));disp(strcat('迭代次数为',num2str(n)));disp(strcat('精度为',num2str(err)));1.2运⾏结果:run('H:\Adocument\matlab\1⽜顿迭代法求零点\newtondiedai.m')输⼊函数表达式:y=x^4-1.4*x^3-0.48*x^2+1.408*x-0.512输⼊迭代初始值:a=1输⼊截⽌误差:delta=0.0005⽤⽜顿法求得零点为0.80072迭代次数为14精度为0.00036062⽜顿迭代法通过⼀系列的迭代操作使得到的结果不断逼近⽅程的实根,给定⼀个初值,每经过⼀次⽜顿迭代,曲线上⼀点的切线与x轴交点就会在区间[a,b]上逐步逼近于根。
上述例⼦中,通过给定初值x=1,经过14次迭代后,得到根为0.80072,精度为0.00036062。
西南交通大学研究生数值分析总复习

记x*表示x的近似值,若x* 0.a1a2 an 10m , (ai 是0,1,,9中的一个数字,a1 0),
*
1 mn 如果 x x 10 , 则称x *近似x时具有n位有效数字。 2
返回
前进
3. 记近似值x*=0.a1a2…an×10m,若要保留五位有效数 字(这是 以后常会用到的),即要求误差限ε<0.5×10m-n, 则n=5;
1 这即要求出满足: 10( n 1) 0.01%的n 2a1
例3(续)
1 由a1 5 10( n 1) 0.01% 0.0001 25 10( n 1) 0.001 n 1 lg 0.001 3 n 4 1 因此,只要对 0.052631578 的近似值取四位 19 1 有效数字为 0.05263 ,则其相对误差限就不 超过0.01% 19
返回
前进
§2 绝对误差、相对误差和有效数字
2.1 绝对误差与相对误差 设 x *为准确值的近似值,记
e xx
*
e x x* er x x
分别称e为近似值x *的绝对误差或误差, er为x*的相对误差。
一般情况下,准确值是不知道的,从而也不能算出绝 对误差e的准确值,但往往可以根据测量工具或计算的情 况估计出e 的取值范围,即估计出绝对误差的一个上界ε :
返回
前进
迭代法是一种重要的逐次逼近法,其基本思想是: 设方程f (x) = 0在区间[a, b]内有一根x*,将方程化为等价 方程x = (x),并在[a, b]内任取一点x0作为初始近似值, 然后按迭代公式计第二章 非线性方程求解算: x ( x ), (k 0,1,2,) (2 - 3)
返回
数值分析试题与答案

一. 填空题(本大题共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()f x dx ⎰,要把区间[0,1]一般要等分 41 份才能保证满足误差小于0.00005的要求(这里(2)()1f x ∞≤);如果知道(2)()0f x >,则 用复化梯形公式计算积分1()f x dx ⎰此实际值 大 (大,小)。
在以10((),())()(),(),()[0,1]g x f x xf x g x dx f x g x C =∈⎰为内积的空间C[0,1]中,与非零常数正交的最高项系数为1的一次多项式是 23x -3. (15分)导出用Euler 法求解 (0)1y yy λ'=⎧⎨=⎩的公式, 并证明它收敛于初值问题的精确解解 Euler 公式 11,1,,,k k k xy y h y k n h nλ--=+==L -----------(5分) ()()1011kk k y h y h y λλ-=+==+L ------------------- (10分)若用复化梯形求积公式计算积分1x I e dx =⎰区间[0,1]应分 2129 等分,即要计算个 2130 点的函数值才能使截断误差不超过71102-⨯;若改用复化Simpson 公式,要达到同样精度区间[0,1]应分12 等分,即要计算个 25 点的函数值1.用Romberg 法计算积分 232x e dx -⎰解 []02()()2b aT f a f b -=+= 9.6410430E-003 10221()222b a a bT T f -+=+= 5.1319070E-00310022243T T S -== 4.6288616E-00322T = 4.4998E-003 21122243T T S -== 4.E-0031002221615S S C -== 4.6588636E-00332T = 4.7817699E-00332222243T T S -== 4.1067038E-0032112221615S S C -== 4.5783515E-0031002226463C C R -== 4.7358037E-0032.用复合Simpson 公式计算积分232x e dx -⎰(n=5)解 44501()4()2()(),625k k h h b aS f a f a kh f a kh f b h ==⎡⎤-=++++++=⎢⎥⎣⎦∑∑5S =4.3630653 E-0033、 对于n+1个节点的插值求积公式()()bnk k k af x dx A f x =≈∑⎰ 至少具有 n 次代数精度. 4、 插值型求积公式()()bnk k k af x dx A f x =≈∑⎰的求积系数之和0nk k A =∑=b-a 5、 证明定积分近似计算的抛物线公式()()4()()22bab a a b f x dx f a f f b -+⎡⎤≈++⎢⎥⎣⎦⎰具有三次代数精度 证明 如果具有4阶导数,则()()4()()22bab a a b f x dx f a f f b -+⎡⎤-++⎢⎥⎣⎦⎰=)(f 2880)a b ()4(5η--(η∈[a,b])因此对不超过3次的多项式f(x)有()()4()()022bab a a b f x dx f a f f b -+⎡⎤-++=⎢⎥⎣⎦⎰即()()4()()22bab a a b f x dx f a f f b -+⎡⎤=++⎢⎥⎣⎦⎰精确成立,对任一4次的多项式f(x)有 因此定积分近似计算的抛物线公式具有三次代数精度 或直接用定义证.6、 试确定常数A ,B ,C 和a ,使得数值积分公式有尽可能高的代数精度。
数值分析试卷及答案

数值分析试卷及答案数值分析试卷一、选择题(共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)题目一: (1)1.1计算结果 (1)1.2结果分析 (1)题目二: (2)2.1计算结果 (2)2.2结果分析 (3)题目三: (4)3.1计算结果 (4)3.2结果分析 (5)总结 (5)附录 (6)Matlab程序: (6)题目一: (6)第一问Newton法: (6)第二问Newton法: (6)第一问Steffensen加速法: (7)第二问Steffensen加速法: (7)题目二 (8)1、Jacobi迭代法 (8)2、Causs-Seidel迭代法 (8)题目三: (9)题目一:分别用牛顿法,及基于牛顿算法下的Steffensen 加速法(1)求ln(x +sin x )=0的根。
初值x0分别取0.1, 1,1.5, 2, 4进行计算。
(2)求sin x =0的根。
初值x0分别取1,1.4,1.6, 1.8,3进行计算。
分析其中遇到的现象与问题。
1.1计算结果求ln(x +sin x )=0的根,可变行为求解x-sinx-1=0的根。
1.2结果分析从结果对比我们可发现牛顿—Steffensen 加速法比牛顿法要收敛的快,牛顿法对于初值的选取特别重要,比如第(1)问中的初值为4的情况,100次内没有迭代出来收敛解,而用Steffensen 加速法,7次迭代可得;在第(2)问中的初值为1.6的情况,收敛解得31.4159,分析其原因应该是x x f cos )('=,x0=1.62π≈,0)('≈x f ;迭代式在迭代过程中会出现分母趋近于0,程序自动停止迭代的情况,此时得到的x 往往非常大,而在第一问中我们如果转化为用x+sinx=1,则可以收敛到结果。
用雅格比法与高斯-赛德尔迭代法解下列方程组Ax=b,研究其收敛性,上机验证理论分析是否正确,比较它们的收敛速度,观察右端项对迭代收敛有无影响。
(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4];b1=[-3,2,4]T,b2=[100,-200,345]T,(2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1];b1=[3,2,1]T,b2=[5,0,-10]T,(3)A行分别为A1=[1,3],A2=[-7,1];b=[4,6]T2.1计算结果初值均为0矩阵带入(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4];b1=[-3,2,4]T,b2=[100,-200,345]T2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1];b1=[3,2,1]T,b2=[5,0,-10]TT2.2结果分析ρ小于1,故方程组雅可比迭代收第一小题的经计算谱半径为5427B(=).0敛。
西南交大数值分析题库插值逼近题库

xkj lk (0)
1, 0, ( 1) n x0 x1...xn
n k 0 n 1 xk lk ( x) n k 0
j 0 j 1,2,...,n j n 1
n f ( n 1) () wn 1 ( x) 其中,wn+1(x)= ( x (n 1)! j 0
n k 0
证明: f ( x)
f (3) () 2 (x 3! k 0 xk )
(3). 三次样条插值与一般分段 3 次多项式插值的区别是_____ (三次样条连续且光 滑,一般分段 3 次连续不一定光滑。) §2. 计算题 (1). (a10 分)依据下列函数值表,建立不超过 3 次的 lagrange 插值多项式 L3(x). x 0 1 2 3 f(x) 1 9 23 3 解:基函数分别为
xn
1
(x
j 0
xj )
n 1 xk lk ( x)
wn 1 ( x) 可见其为 n 次多项式,并且可得其最高次系数为
(x0+…+ xn) (5). 设函数 f(x)是 k 次多项式,对于互异节点 x1,…, xn,, 证明当 n>k 时,差商 f [x, x1,…,xn]0,当 nk 时,该差商是 k-n 次多项式。 证明:因 f [ x0 , x1 , , xn ]
多项式 P1(x)在子区间[a,b]上的余项估计式,再估计最值即可。
f ( x) P 1 ( x) f () ( x a)( x b) 2!
x3 , c( x 1)
2 3
hi2 max f // ( x) a 8 x b
0 x x 1 2
x [a, b]
(12). s(x)=
2
已知 s(x)是[0,2]上的已知自然边界条件的三次样条函数,试确定
《数值分析》练习题及答案解析

《数值分析》练习题及答案解析第一章 绪论主要考查点:有效数字,相对误差、绝对误差定义及关系;误差分类;误差控制的基本原则;。
1. 3.142和3.141分别作为π的近似数具有( )和( )位有效数字.A .4和3B .3和2C .3和4D .4和4 答案:A2. 设 2.3149541...x *=,取5位有效数字,则所得的近似值x=___________ .答案:2.31503.若近似数2*103400.0-⨯=x 的绝对误差限为5105.0-⨯,那么近似数有几位有效数字 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
4 . 14159.3=π具有4位有效数字的近似值是多少?解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。
第二章 非线性方程求根 主要考查点:二分法N 步后根所在的区间,及给定精度下二分的次数计算;非线性方程一般迭代格式的构造,(局部)收敛性的判断,迭代次数计算; 牛顿迭代格式构造;求收敛阶;1.用二分法求方程012=--x x 的正根,要求误差小于0.05。
(二分法)解:1)(2--=x x x f ,01)0(<-=f ,01)2(>=f ,)(x f 在[0,2]连续,故[0,2]为函数的有根区间。
"(1)计算01)1(<-=f ,故有根区间为[1,2]。
(2)计算041123)23()23(2<-=--=f ,故有根区间为]2,23[。
(3)计算0165147)47()47(2>=--=f ,故有根区间为]47,23[。
(4)计算06411813)813()813(2>=--=f ,故有根区间为]813,23[。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
考试目标及考试大纲本题库的编纂目的旨在给出多套试题,每套试题的考查范围及难度配置均基于“水平测试”原则,按照教学大纲和教学内容的要求,通过对每套试题的解答,可以客观公正的评定出学生对本课程理论体系和应用方法等主要内容的掌握水平。
通过它可以有效鉴别和分离不同层次的学习水平,从而可以对学生的学习成绩给出客观的综合评定结果。
本题库力求作到能够较为全面的覆盖教学内容,同时突显对重点概念、重点内容和重要方法的考查。
考试内容包括以下部分:绪论与误差:绝对误差与相对误差、有效数字、误差传播分析的全微分法、相对误差估计的条件数方法、数值运算的若干原则、数值稳定的算法、常用数值稳定技术。
非线性方程求解:方程的近似解之二分法、迭代法全局收敛性和局部收敛定理、迭代法误差的事前估计法和事后估计法、迭代过程的收敛速度、r 阶收敛定理、Aitken加速法、Ne w to n法与弦截法、牛顿局部收敛性、Ne w to n收敛的充分条件、单双点割线法(弦截法)、重根加速收敛法。
解线性方程组的直接法:高斯消元法极其充分条件、全主元消去法、列主元消去法、高斯-若当消元法、求逆阵、各种消元运算的数量级估计与比较、矩阵三角分解法、Doolittle 和Crout三角分解的充分条件、分解法的手工操作、平方根法、Cholesky分解、改进的平方根法(免去开方)、可追赶的充分条件及适用范围、计算复杂性比较、严格对角占优阵。
解线性方程组迭代法:向量和矩阵的范数、常用向量范数的计算、范数的等价性、矩阵的相容范数、诱导范数、常用范数的计算;方程组的性态和条件数、基于条件数误差估计与迭代精度改善方法;雅可比(Jacobi)迭代法、Gauss-Seidel迭代法、迭代收敛与谱半径的关系、谱判别法、基于范数的迭代判敛法和误差估计、迭代法误差的事前估计法和事后估计法;严格对角占优阵迭代收敛的有关结论;松弛法及其迭代判敛法。
插值法:插值问题和插值法概念、插值多项式的存在性和唯一性、插值余项定理;Lagrange插值多项式;差商的概念和性质、差商与导数之间的关系、差商表的计算、牛顿(Newton)插值多项式;差分、差分表、等距节点插值公式;Hermite插值及其插值基函数、误差估计、插值龙格(Runge)现象;分段线性插值、分段抛物插值、分段插值的余项及收敛性和稳定性;样条曲线与样条函数、三次样条插值函数的三转角法和三弯矩法。
曲线拟合和函数逼近:最小二乘法原理和多项式拟合、函数线性无关概念、法方程有唯一解的条件、一般最小二乘法问题、最小二乘拟合函数定理、可化为线性拟合问题的常见函数类;正交多项式曲线拟合、离散正交多项式的三项递推法。
最佳一致逼近问题、最佳一致逼近多项式、切比雪夫多项式、切比雪夫最小偏差定理、切比雪夫多项式的应用(插值余项近似极小化、多项式降幂)。
本段加黑斜体内容理论推导可以淡化,但概念需要理解。
数值积分与微分:求积公式代数精度、代数精度的简单判法、插值型求积公式、插值型求积公式的代数精度;牛顿一柯特斯(Newton-Cotes)公式、辛卜生(Simpson)公式、几种低价牛顿一柯特斯求积公式的余项;牛顿一柯特斯公式的和收敛性、复化梯形公式及其截断误差、复化Simpson公式及其截断误差、龙贝格(Romberg)求积法、外推加速法、高斯型求积公式、插值型求积公式的最高代数精度、高斯点的充分必要条件。
正交多项式的构造方法、高斯公式权系数的建立、Gauss-Legendre公式的节点和系数。
本段加黑斜体内容理论推导可以淡化,但概念需要理解。
常微分方程数值解:常微分方程初值问题数值解法之欧拉及其改进法、龙格—库塔法、阿当姆斯方法。
本套题库均采用闭卷考试,卷面总分为100分。
试题形式分为判别正误、多项选择、填空、解答和证明等多种题型。
其中判断题、多项选择题和填空题覆盖整个内容范围,题量多而广,重点集中在基本概念、公式和方法的构建与处理思想等方面,此类题型主要用于考查学生对整体内容的理解与掌握情况;解答题重点放在主要的计算技术和方法的具体实现过程,主要考查学生对主要计算技术、技巧和方法理解与掌握情况;证明题主要集中在主要的计算技术和方法的分析过程,主要考查学生的理论分析能力和知识的综合运用能力。
本课程的考试方法与要求:期末闭卷考试,按时完成上机习题。
学习合格条件:考试卷面成绩 60且上机习题符合要求,二者缺一不可。
综合成绩:原则上=卷面成绩,但可参考上机习题完成情况作微调。
填空题1 绪论(1). 要使20的近似值的相对误差限≤0.1%, 应至少取___4____位有效数字。
20=0.4…⨯10, a 1=4, εr ≤121a ⨯10-(n-1)< 0.1% ,故可取n ≥4, 即4位有效数字。
(2). 要使20的近似值的相对误差限≤0.1%, 应至少取___4___位有效数字,此时的绝对误差限为31102(3). 设y =f (x 1,x 2) 若x 1,x 2,的近似值分别为x 1*, x 2*,令y *=f (x 1*,x 2*)作为y 的近似值,其绝对误差限的估计式为: ε ≤| |f (x 1*,x 2*)|x 1-x*1|+ |f (x 1*,x 2*)|x 2-x*2| (4). 计算 f=(2-1)6 , 取2=1.4 , 利用下列算式,那个得到的结果最好?答:__C_____.(A)6121)(-, (B) (3-22)2, (C)32231)(+, (D) 99-702(5). 要使17的近似值的相对误差限≤0.1%, 应至少取_________位有效数字?17=0.4…⨯10, a 1=4, εr ≤121a ⨯10-(n-1)< 0.1% 故可取n ≥3.097, 即4位有效数字。
(6). 设x =3.214, y =3.213,欲计算u =y x -, 请给出一个精度较高的算式u =.u=yx y x +-(7).设x =3.214, y =3.213,欲计算u =y x -, 请给出一个精度较高的算式u = .u=yx y x +-(8).设y =f (x 1,x 2) 若x 1,x 2,的近似值分别为x 1*, x 2*,令y *=f (x 1*,x 2*)作为y 的近似值,其绝对误差限的估计式为: ε ≤| |f (x 1*,x 2*)|x 1-x*1|+ |f (x 1*,x 2*)|x 2-x*2|;2 方程根(9). 设迭代函数ϕ(x )在x *邻近有r (≥1)阶连续导数,且x * = ϕ(x *),并且有ϕ(k )(x *)=0(k =1,…,r -1),但ϕ(r ) (x *)≠0,则x n +1=ϕ(x n )产生的序列{ x n }的收敛阶数为___r___(10). 称序列{x n }是p 阶收敛的如果c x x x x pn n n =--+∞→**lim1(11). 用牛顿法求 f (x)=0 的n 重根,为了提高收敛速度,通常转化为求另一函数u(x)=0的单根,u(x)=()()f x f x '(12). 用Newton 法求方程f (x )=x 3+10x -20=0 的根,取初值x 0= 1.5, 则x 1= ________ 解x 1=1.5970149 (13). 用牛顿法解方程0123=--x x 的迭代格式为_______________ 解 kk k k k k x x x x x x 2312231----=+ (14). 迭代过程)(1k k x x ϕ=+收敛的充分条件是)(x ϕ' ≤ 1.___(15). 用Newton 法求方程f(x)=x 3+10x-20=0 的根,取初值x 0= 1.5, 则x 1= 1.5970149 (16). 用牛顿法解方程0123=--x x 的迭代格式为(17). 用Newton 法求方程f (x )=x 3+10x -20=0 的根,取初值x 0= 1.5, 则x 1= ________ 解x 1=1.5970149(18). 迭代公式x k +1=x k (x k 2+3a )/(3x k 2+a )是求a 1/2的 (12) 阶方法3方程组(19). 矩阵的 LU 分解中L 是一个 _为单位下三角阵,而U 是一个上三角阵____。
(20). 设线性方程组的系数矩阵为A =⎪⎪⎪⎪⎪⎭⎫⎝⎛-6847153131483412,全主元消元法的第一次可选的主元素为 -8,或8___,第二次可选的主元素为 8+7/8或-8-7/8 ____. 列主元消元法的第一次主元素为 _-8_________;第二次主元素为(用小数表示) 7.5_____;(21). 在方阵A 的LU 分解中, 方阵A 的所有顺序主子不为零,是方阵A 能进行LU 分解的充 分 (充分,必要)条件;严格行对角占优阵 能__(能,不能)进行LU 分解;非奇异矩阵___不一定___(一定,不一定)能进行LU 分解。
(22). 设A 是正定矩阵,则A 的cholesky 的分解 唯一 (唯一,不唯一).(23). 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=2021012a a A ,为使A 可分解为A=LL T ,其中L 是对角线元素为正的下三角形矩阵,则a 的取值范围是 ,取a=1,则L= 。
(24). 解 )3,3(-∈a ,⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡32320023210024迭代(1). ⎥⎦⎤⎢⎣⎡-=3211A ,则=1||||A ,=2||||A ,=∞||||A ; 答:4,3.6180340,5;(2). 已知方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡2121132.021b b x x ,则解此方程组的Jacobi 迭代法___是___收敛(填“是”或“不”)。
(3). 给定方程组 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--111 211111112321x x x 记此方程组的Jacobi 迭代矩阵为B J =(a ij )3⨯3,则a 23= -1; , 且 相应的Jacobi 迭代序列是__发散_____的。
(4). 设3()1f x x ,则()f x 关于[0,1]C的f 1, 2f(5). ⎥⎦⎤⎢⎣⎡-=1301A ,则)1,)1(|(|1)(,4||||2,121=-=-==λλλρA I A A (6). R n 上的两个范数||x||p , ||x||q 等价指的是_∃C,D ∈R,_C_||x||q _≤||x||p ≤D ||x||q _; R n 上的两个范数_一定____是等价的。
(选填“一定”或“不一定”)。
(7). Tx )12,4,0,3(-=,则=1||||x 19 ,=2||||x 13____,=∞||||x ____12 ; (8). 已知方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡2121132.021b b x x ,则解此方程组的Jacobi 迭代法___收敛(填“收敛”或“发散”),12∞解(10). 已知方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡2121132.021b b x x ,则解此方程组的Jacobi 迭代法_____________收敛(填“是”或“不”),解 (3)因⎥⎦⎤⎢⎣⎡=132.021A 的Jacobi 迭代矩阵⎥⎦⎤⎢⎣⎡=032.020B ,8.0)(=B ρ,故Jacobi 迭代是收敛的,(11). 已知方程组⎩⎨⎧=-=+26203825y x y x ,其雅可比法的迭代矩阵是______________,高斯-塞德尔法的迭代格式是________________;解 ⎪⎪⎩⎪⎪⎨⎧+=+-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+++10132035852,0203520)1()1()()1(k k k k x y y x (12). 已知方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡2121132.021b b x x ,则解此方程组的Jacobi 迭代法_____________收敛(填“是”或“不”),解 因⎥⎦⎤⎢⎣⎡=132.021A 的Jacobi 迭代矩阵⎥⎦⎤⎢⎣⎡=032.020B ,8.0)(=B ρ,故Jacobi 迭代是收敛的,(13). 已知方程组⎩⎨⎧=-=+26203825y x y x ,其雅可比法的迭代矩阵是______________,高斯-塞德尔法的迭代格式是________________;解(14). ⎥⎥⎦⎤⎢⎢⎣⎡=21010a A ,要使0lim =∞→k k A ,a 应满足___________; 解 1<a12∞⎥⎦⎤⎢⎣⎡-=1301A ,则=1||||A ,=)(A ρ 。