中南大学数值分析试题7
数值分析试题集
..数值分析试题集（试卷一）一（ 10 分）已知 x 1* 1.3409 ，x 2* 1.0125 都是由四舍五入产生的近似值， 判断 x 1*x 2* 及 x 1* x 2*有几位有效数字。
二（ 10 分）由下表求插值多项式x 01 2 y2 34 y1－ 1三（ 15 分）设 f ( x)C 4 [a,b] ， H （ x ）是满足下列条件的三次多项式H (a) f (a) , H (b) f (b) , H (c)f (c) , H (c) f (c)( a c b )求 f (x)H ( x) ，并证明之。
12四（ 15 分）计算13 dx ，10 2。
x五（ 15 分）在 [0，2]上取 x 0 0 , x 1 1 , x 22 ，用二种方法构造求积公式，并给出其公式的代数精度。
六（ 10 分）证明改进的尢拉法的精度是 2 阶的。
七（ 10 分）对模型 yy , 0 ，讨论改进的尢拉法的稳定性。
八（ 15分）求方程 x 34x 2 7x 1 0 在 -1.2 附近的近似值，10 3。
-----------------------------------------------------------------------------------------------------------------------------（试卷二）一填空（ 4*2 分）1 {k ( x) } k 0 是区间 [0， 1]上的权函数为( x) x 2 的最高项系数为 1 的正交多项式族，其中10 (x)1，则x0 ( x) dx ------------------- ， 1 ( x) ------------------。
2 12 A，则 A1 4----------- ，( A) ----------------- 。
a 1 2 时， A 可作 LU 分解。
3 设 A，当 a 满足条件 ---------------- 14..4 设非线性方程 f ( x) (x33x23x1)( x 3) 0 ，其根 x1* 3 ， x2*1，则求 x1* 的近似值时，二阶局部收敛的牛顿迭代公式是--------------------------- 。
数值分析第7章答案
第七章非线性方程求根一、重点内容提要 （一）问题简介 求单变量函数方程(7.1) 的根是指求（实数或复数），使得.称为方程(7.1)的根,也称为函数的零点.若可以分解为其中m 为正整数,满足,则是方程(7.1)的根.当m=1时,称为单根;当m>1时,称为m 重根.若充分光滑,是方程(7.1)的m 重根,则有(1)()(*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若在[a,b]上连续且,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法设在[a,b]上连续,,则在(a,b)内有根.再设在(a,b)内仅有一个根.令,计算和.若则,结束计算;若,则令,得新的有根区间;若,则令,得新的有根区间.,.再令计算,同上法得出新的有根区间,如此反复进行,可得一有根区间套且110011*,0,1,2,...,()...()22n n n n n n a x b n b a b a b a --<<=-=-==-.故因此,可作为的近似根,且有误差估计 (7.2) 2.迭代法将方程式(7.1)等价变形为 (7.3)若要求满足则;反之亦然.称为函数的一个不动点.求方程(7.1)的根等价于求的不动点由式(7.3)产生的不动点迭代关系式(也称简单迭代法)为 (7.4)函数称为迭代函数.如果对任意,由式(7.4)产生的序列有极限 则称不动点迭代法(7.4)收敛.定理7.1(不动点存在性定理)设满足以下两个条件: 1.对任意有2.存在正常数,使对任意,都有 (7.5) 则在上存在惟一的不动点.定理7.2(不动点迭代法的全局收敛性定理)设满足定理7.1中的两个条件,则对任意,由(7.4)式得到的迭代序列收敛.到的不动点,并有误差估计式 (7.6) 和 (7.7)定理7.3(不动点迭代法的局部收敛性定理)设为的不动点,在的某个邻域连续,且,则迭代法(7.4)局部收敛.收敛阶的概念 设迭代过程(7.4)收敛于方程的根,如果迭代误差当时成产下列渐近关系式(7.8)则称该迭代过程是p 阶收敛的.特别地,p=1时称线性收敛,p>1时称超线性收敛,p=2时称平方收敛.定理7.4(收敛阶定理)对于迭代过程(7.4),如果在所求根的邻近连续,并且 (7.9)则该迭代过程在点的邻近是收敛的,并有(7.10)斯蒂芬森(Steffensen)迭代法 当不动点迭代法(7.4)只有线性收敛阶,甚至于不收敛时,可用斯蒂芬森迭代法进行加速.具体公式为 (7.11) 此法也可写成如下不动点迭代式(7.12)定理7.5(斯蒂芬森迭代收敛定理) 设为式(7.12)中的不动点,则是的不动点;设存在,,则是的不动点,则斯蒂芬森迭代法(7.11)是2阶收敛的. 3.牛顿迭代法牛顿迭代法是一种特殊的不动点迭代法,其计算公式为 其迭代函数为 (7.13)牛顿迭代法的收敛速度 当时,容易证明,,,由定理7.4知,牛顿迭代法是平方收敛的,且(7.14)重根情形的牛顿迭代法 当是的m 重根时,迭代函数在处的导数,且.所以牛顿迭代法求重根只是线性收敛.若的重数m 知道,则迭代式 (7.15)求重根二阶收敛.当m 未知时,一定是函数的单重零点,此时迭代式1()()'()'()['()]()''()0,1,2,...k k k k k k k k k k x f x f x x x x x f x f x f x k μμ+=-=--= (7.16)也是二阶收敛的.简化牛顿法 如下迭代法 称为简化牛顿法或平行弦法.牛顿下山法 为防止迭代不收敛,可采用牛顿下山法.具体方法见教材. 4.弦截法将牛顿迭代法(7.13)中的用在,处的一阶差商来代替,即可得弦截法 (7.17)定理7.6假设在其零点的邻域内具有二阶连续导数,且对任意有,又初值,,则当邻域充分小时,弦截法(7.17)将按阶收敛到.这里p 是方程的正根. 5.抛物线法弦截法可以理解为用过两点的直线方程的根近似替的根.若已知的三个近似根,,用过的抛物线方程的根近似代替的根,所得的迭代法称为抛物线法,也称密勒(Muller)法.当在的邻近有三阶连续导数,,则抛物线法局部收敛,且收敛阶为.二、知识结构图10[1,2]1x x --=≤≤--∈3-3-6k k 32三、常考题型及典型题精解例7-1 证明方程x 在上有一个实根x*,并用二分法求这个根,要求|x -x*|10.若要求|x -x*|10,需二分区间[1,2]多少次?解 设f(x)=x ,则f(1)=-1<0,f(2)=5>0,故方程f(x)=0在[1,2]上有根x*.又因f'(x)=3x -1,所以当x [1,2]时,f'(x)>0,即f (x)=0在[1,2]上有惟一实根x*.用二分法计算结果如表7-1所示.k 0 1 2 3 4 5 6 7 8 9 1 1 1.25 1.25 1.3125 1.3125 1.3125 1.3204 1.3243 1.3243 2 1.5 1.5 1.375 1.375 1.13438 1.3282 1.3282 1.3282 1.32631.5 1.25 1.375 1.3125 1.3438 1.3282 1.3204 1.3243 1.3263 1.3253+ - + - + + - - + +610x e -≤≤⨯≤≤≤≤≥∈-3-39910-6k k k+101此时x =1.3253满足|x -x*|0.9771010,可以作为x*的近2似值.1若要求|x -x*|,只需|x -x*|10即可,解得k+119.932,2即只需把[1,2]二分20次就能满足精度要求.例7-2 已知函数方程(x-2)=1,(1)确定有根区间[a,b];(2)构造不动点迭代公式使之对任意初始近似x [a,b],31|10.k x ---<k 迭代方法均收敛;(3)用所构造的公式计算根的近似值,要求|x1lim lim x x x x x e e e e →+∞→-∞∞∞∞∈解 (1)令f(x)=(x-2)-1,由于f(2)=-1<0,f(3)=-1>0,因此区间[2,3]是方程f(x)=0的一个有根区间.又因f'(x)=(x-1),f(x)=+,f(x)=-1,f'(1)=--1<0,当x>1时f(x)单增,x<1时f(x)单减,故f(x)=0在(-,+)内有且仅有一根x*,即x*[2,3].2'k k x x x x x x e e e e e e e ϕϕϕ-----∈∈≤≤≤∀∈k+100k+1(2)将(x-2)=1等价变形为x=2+,x [2,3].则(x)=2+.由于当x [2,3]时2(x)3,|(x)|=|-|<1故不动点迭代法x =2+,k=0,1,2,...,对x [2,3]均收敛.(3)取x =2.5,利用x =2+进行迭代计算,结果如表7-2所示.473cos 3120cos c k x x x ϕ--+=∈≤4k+10-30k+1k+1k 例 考虑求解方程2的迭代公式2x =4+,k=0,1,2,...3(1)试证:对任意初始近似x R,该方法收敛;(2)取x =4,求根的近似值x ,要求|x -x |10;(3)所给方法的收敛阶是多少?2解 (1)由迭代公式知,迭代函数(x)=4+3{}os ,(,).|'sin |1(,)x x x ϕϕϕ∈-∞+∞≤<-∞+∞∀∈0k 022由于(x)的值域介于(4-)与(4+)之间,且3322(x)|=|-33故根据定理7.1,7.2知,(x)在内存在惟一的不动点x*,且对x R,迭代公式得到的序列x 收敛于x*.(2) 取x =4,迭代计算结果如表7-3所示.此时已满足误差要求，即（3）由于，故根据定理7 .4知方法是线性收敛的，并且有。
2009硕士生数值分析试题
2009年中南大学硕士研究生“数值分析”课程试题（闭卷，可自带计算器，120分钟）一、设分段多项式（8分）S(x)=⎩⎨⎧≤≤-++≤≤+21,1210,2323x cx bx x x x x是以x = 0，1，2为节点的三次样条函数，试确定系数b ，c 的值。
二、设)3()(2-+=Φx c x x ，问如何选取 c 能保证迭代法)(1k k x x Φ=+具有局部收敛性。
（8分）三、利用列主元素消去法解方程：（8分）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∙⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---453311294642321x x x四、求()xf x e =在[0，1]上的一次最佳平方逼近多项式。
（8分）五、已知函数表如下，试用抛物插值求125的近似值，估计截断误差。
（14分）利用Romberg 方法计算积分：dx xx⎰1sin （计算到第一个Romberg 值）。
七、利用改进的Euler 法求解如下初值问题：（14分）,1)0(6.00,2)('⎪⎩⎪⎨⎧=≤≤-=y x y x y x y 取步长 h =0.3 。
八、利用反幂法和矩阵的LU 分解技术求下列矩阵A 接近于 p = -6.4 的特征值及其特征向量（保留3位小数迭代计算2次，分解后取1(1,1,1)T Uv =），-1212-4111-6A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦。
（14分）九、方程组AX=b ， 其中10101a A a a a ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦。
试建立解线性方程组AX=b 的雅可比迭代法和高斯-赛德尔迭代法，并讨论a 取何值时迭代收敛。
（12分）。
数值分析复习题参考答案
x1 )
h
2
x 0 x x1 6
4
所以， R x
h 10
2
8
解得， h 0 . 000383
4. 习题(第二章) 7
5. 习题(第二章) 9
6. 习题(第二章) 11
7. 习题(第二章) 13
8. 习题(第二章) 14
9. 习题(第二章) 20
10. 习题(第四章) 1
2
, k 0 ,1, 2 2 3 2a 3x
3
此时， ( x )
2x a 3x
, '( x) 2a
所以， ' ( 3 a )
2 3
3(
3
a)
3
0 1, 所以该迭代公式收敛。
21. 习题(第七章) 13
本题没有给出精度要求， 但x3与x2之间的差为 已经很小了，足以满足 精度。
[ f ( x n , y n ) f ( x n 1 , y n 1 )]
( 3 ) 基于 Taylor 展开法：
y ( x n 1 ) y ( x n h ) y ( x n ) y ' ( x n ) h
h
2
2
y ''( xn )
取 y ( x n 1 ) y ( x n ) y ' ( x n ) h ，即 y n 1 y n hf ( x n , y n )
k 个点的值
求解隐式：先用欧拉公 求解多步法：单步法开
式求得一个初步的近似 表头，然后预报
修正 校正 修正。
( 其实只要给出公式会用
就行！！ )
数值分析考试题
山东科技大学2008-2009学年第一学期《数值分析》考试一、设x =9.1234, y =10.486均具有5位有效数字。
试分析x - y和x3 y啲绝对误差限和相对误差限。
二、求一条拟合3点A(0,1), B(1,3),C(2,2)的直线。
三、设n _ 2为正整数，c为正数，记x*二n.c1) 说明不能用下面的迭代格式1 _nx k 1 = cx k ,k =Q1,2：= = " =求x*的近似值。
2) 构造一个可以求x*的迭代格式，证明所构造迭代格式的收敛性，并指出收敛阶数四、给定线性方程组_4 -1 0卩1 一2〕-1 a 1 x2 = 64」］X3」：2J】0 1其中a为非零常数。
1) 写出Jacobi迭代格式与Gauss-Seidel迭代格式并分析其收敛性。
2) 分析a在什么范围取值时以上迭代格式收敛。
五、做一个5次多项式H (x)使得H(1) =3,H (2) = —1, H(4) =3,H'(1) =2, H'(2) =1, H *(2) =2,六、求f (x) =x2在区间0,1上的一次最佳一致逼近多项式。
七、给定积分公式:1f(x)d x ：Af (-1) Bf (0) f (1)■ -41) 试确定求积系数A,B,C，使其具有尽可能高的代数精度，并指出其代数精度。
2) 试判断该求积公式是否为高斯型求积公式，并说明理由。
3) ................................................................................................ 将区间-1,作n等分，并记h=2,X j =-1 ih,i =0,1,............................................................ ，n,利用该求积公式n 构造一个复化求积公式。
中南大学数值分析作业完整版~
数值分析数学实验报告姓名：XX学号：xx指导老师：***专业班级：xx目录1. 高斯消去法 (3)2. LU分解 (6)3. 用牛顿法求积分 (10)4. 用复化梯形法求积分 (12)5. 用复化辛普森法、复化辛普森变步长法求积分 (13)6. 节点加密复化梯形公式 (16)7. 龙贝格积分 (17)8. 欧拉方法、休恩方法、泰勒方法、龙格-库塔方法 (20)一.高斯消去法x(i)=(b(i)-sum)/A(i,i);end4：实验结果：（1）高斯消去法（2）高斯列主元消去法5：实验总结这两个程序让我对高斯消去法有了更深刻的理解，能更灵活的运用各种基础函数与矩阵的运算来进行求解，参考了书上的已有程序流程图，程思想需要进一步完善，需要对函数的进一步熟悉。
姓名：xx2012年3月25日二. LU分解m=j;endendif m~=ifor k=1:nc(k)=A(i,k);A(i,k)=A(m,k);A(m,k)=c(k);endt=b(i);b(i)=b(m);b(m)=t;endfor j=i:nfor k=1:i-1M(k)=L(i,k)*U(k,j);endU(i,j)=A(i,j)-sum(M);endfor j=i+1:nfor k=1:i-1M(k)=L(j,k)*U(k,i);endL(j,i)=(A(j,i)-sum(M))/U(i,i);endendx=U\(L\b);4：实验结果：（1）普通LU分解三. 用牛顿法求积分（2）列主元LU 分解5：实验总结L U 分解在上学期已经学习过，这次的实验让我对LU 分解有了更深的了解，又掌握了一种解线性方程组的好方法。
姓名：XX2012年3月29日 学号XX 班级 XX 姓名XX 指导教师 易昆南 实验题目 用牛顿法求积分评 分1、设计（实习）目的：1． 进一步了解牛顿法及其应用 2．进一步理解牛顿法求积分的思想2、实验内容：用牛顿法求函数x x x x f ++=23)(的积分 3.详细设计：function y=newton(a,b,n) x=a:(b-a)/n:b; %插值节 y=x.^3+x.^2+x;四.用复化梯形法求积分五. 用复化辛普森法、复化辛普森变步长法求积分六.节点加密复化梯形公式七．龙贝格积分八．欧拉方法、休恩方法、泰勒方法、龙格-库塔方法4.龙格-库塔方法>> [x1 y1] =lungkuta(1)[x2 y2] =lungkuta(1/2)[x3 y3] =lungkuta(1/4)[x4 y4] =lungkuta(1/8)plot(x1,y1,'-',x2,y2,'r',x3,y3,'g',x4,y4,'b')x1 =1 2 3y1 =0.769531250000000 1.043746948242188 1.615647614002228x2 =0.500000000000000 1.000000000000000 1.500000000000000 2.000000000000000 2.500000000000000 3.000000000000000y2 =0.935424804687500 0.466060072183609 0.450289419204637 0.558880619571974 0.701568099555500 0.853603753857624x3 =Columns 1 through 65：实验总结用数值分析中的方法编程求积分它能帮助我们简化繁琐又难以计算的数学问题。
数值分析习题(含标准答案)
]第一章 绪论姓名 学号 班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1若误差限为5105.0-⨯,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需！41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（ , ）之间的任意数，都具有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π=的绝对误差限与相对误差限。
（误差限的计算）解：*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ绝对误差限为πππ252.051.02052)5,20(),(2=⨯⋅+⨯⋅⋅⋅≤-v r h v相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6设x 的相对误差为%a ,求nx y =的相对误差。
数值分析作业答案(第7章part2)
7.2.为求方程0123=--x x 在5.10=x 附近的一个根，设将方程改写成下列等价形式，并建立相应的迭代公式。
(1).2/11x x +=，迭代公式21/11k k x x +=+；(2).123+=x x ，迭代公式3211+=+k k x x ； (3).112-=x x ，迭代公式1/11-=+k k x x 。
试分析每种迭代公式的收敛性，并选取一种公式求出具有四位有效数字的近似根。
解 考虑5.10=x 的领域]6.1,3.1[。
(1).当]6.1,3.1[∈x 时，]6.1,3.1[11)(2∈+=xx ϕ，1910.03.122)('33<=≈≤-=L x x ϕ，故迭代2111kk x x +=+在]6.1,3.1[上整体收敛。
(2).当]6.1,3.1[∈x 时，()]6.1,3.1[1)(3/12∈+=x x ϕ，1522.0)3.11(6.132)1(32)('3/223/22<=≈+<+=L x x x ϕ， 故迭代3211+=+k k x x 在]6.1,3.1[上整体收敛。
(3).当]6.1,3.1[∈x 时，11)(-=x x ϕ，1)16.1(21)1(21)('2/3>->--=x x ϕ，故迭代1/11-=+k k x x 发散。
7.4.给定函数)(x f ，设对一切x ，)('x f 存在且M x f m ≤<<)('0，证明对于范围M /20<<λ内的任意定数λ，迭代过程)(1x f x x k k λ-=+均收敛于0)(=x f 的根*x 。
证明 由于0)('>x f ，故)(x f 为单调函数因此方程0)(=x f 的根*x 是唯一的。
迭代函数)()(x f x x λϕ-=，)('1)('x f x λϕ-=。
由M x f m ≤<<)('0及M /20<<λ，得：2)('0<≤≤<M x f m λλλ11)('111<-≤-≤-<-m x f M λλλ故1}1,1max{)('<--=≤M m L x λλϕ因此可得0*0*1*→-≤≤-≤--x x L x x L x x k k k Λ )(∞→k即*lim x x k k =∞→。
数值分析试题及答案
数值分析试题及答案一、选择题（每题2分，共20分）1. 以下哪个算法是数值分析中用于求解线性方程组的直接方法？A. 牛顿法B. 高斯消元法C. 梯度下降法D. 蒙特卡洛方法答案：B2. 插值法中，拉格朗日插值法和牛顿插值法的共同点是：A. 都是多项式插值B. 都使用差商C. 都只适用于等距节点D. 都需要预先知道所有数据点答案：A3. 在数值积分中，辛普森(Simpson)公式比梯形公式的误差：A. 更大B. 更小C. 相同D. 无法比较答案：B4. 以下哪个是数值稳定性分析中常用的方法？A. 条件数B. 收敛性C. 收敛速度D. 误差分析答案：A5. 在求解常微分方程的数值解时，欧拉方法属于：A. 单步法B. 多步法C. 隐式方法D. 显式方法答案：A6. 以下哪个是数值分析中求解非线性方程的迭代方法？A. 高斯-约当消元法B. 牛顿-拉弗森方法C. 雅可比迭代法D. 高斯-赛德尔迭代法答案：B7. 线性插值公式中，如果给定两个点\( (x_0, y_0) \)和\( (x_1, y_1) \)，插值多项式是：A. \( y = y_0 + \frac{y_1 - y_0}{x_1 - x_0}(x - x_0) \)B. \( y = y_0 + \frac{y_1 - y_0}{x_0 - x_1}(x - x_0) \)C. \( y = y_0 + \frac{x - x_0}{x_1 - x_0}(y_1 - y_0) \)D. \( y = y_1 + \frac{x_1 - x}{x_1 - x_0}(y_0 - y_1) \)答案：C8. 以下哪个是数值分析中用于求解特征值问题的算法？A. 幂法B. 共轭梯度法C. 牛顿法D. 欧拉法答案：A9. 在数值微分中，使用有限差分法来近似导数时，中心差分法的误差：A. 与步长成正比B. 与步长的平方成正比C. 与步长的立方成正比D. 与步长的四次方成正比答案：B10. 以下哪个是数值分析中用于求解线性最小二乘问题的算法？A. 梯度下降法B. 牛顿法C. 奇异值分解法D. 共轭梯度法答案：C二、简答题（每题10分，共30分）1. 简述数值分析中病态问题的特点及其对算法的影响。
《数值分析》练习题及答案解析
《数值分析》练习题及答案解析第一章 绪论主要考查点：有效数字，相对误差、绝对误差定义及关系；误差分类；误差控制的基本原则；。
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[。
数值分析练习题附答案
目录一、绪论------------------------------------------------------------------------------------- 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 分量迭代格式以及相应的迭代矩阵，并求解。
数值分析试题库与答案解析
模 拟 试 卷（一）一、填空题（每小题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）内。
数值分析计算方法试题集及答案
数值分析复习试题第一章 绪论 一. 填空题 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,如果取01.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。
数值分析试题及答案
数值分析试题及答案一、选择题（每题3分，共30分）1. 下列关于数值分析的说法，错误的是（）。
A. 数值分析是研究数值方法的科学B. 数值分析是研究数值方法的数学理论C. 数值分析是研究数值方法的误差分析D. 数值分析是研究数值方法的数学理论、误差分析及数值方法的实现答案：B2. 在数值分析中，插值法主要用于（）。
A. 求解微分方程B. 求解积分方程C. 求解线性方程组D. 通过已知数据点构造一个多项式答案：D3. 线性方程组的解法中，高斯消元法属于（）。
A. 直接方法B. 迭代方法C. 矩阵分解方法D. 特征值方法答案：A4. 牛顿法（Newton's method）是一种（）。
A. 插值方法B. 拟合方法C. 迭代方法D. 优化方法答案：C5. 在数值分析中，下列哪种方法用于求解非线性方程的根？A. 高斯消元法B. 牛顿法C. 雅可比方法D. 斯托尔-温格尔方法答案：B6. 下列关于误差的说法，正确的是（）。
A. 绝对误差总是大于相对误差B. 相对误差总是小于绝对误差C. 误差是不可避免的D. 误差总是可以消除的答案：C7. 在数值分析中，下列哪个概念与数值稳定性无关？A. 条件数B. 截断误差C. 舍入误差D. 插值多项式的阶数答案：D8. 用泰勒级数展开函数f(x)=e^x，下列哪一项是正确的？A. f(x) = 1 + x + x^2/2! + x^3/3! + ...B. f(x) = 1 - x + x^2/2! - x^3/3! + ...C. f(x) = x + x^2/2 + x^3/6 + ...D. f(x) = x - x^2/2 + x^3/6 - ...答案：A9. 插值多项式的次数最多为（）。
A. n-1B. nC. n+1D. 2n答案：B10. 下列关于数值积分的说法，错误的是（）。
A. 梯形法则是一种数值积分方法B. 辛普森法则是一种数值积分方法C. 龙格法则是数值积分方法中的一种D. 数值积分方法总是精确的答案：D二、填空题（每题3分，共15分）1. 在数值分析中，条件数是衡量问题的______。
数值分析试题及答案
数值分析试题及答案一、单项选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。
数值分析试题及答案.
一、单项选择题（每小题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 +⎧⎛⎫⎪ ⎪=+⎪ ⎪=+⎨⎝⎭⎪=⎪⎩L得 分 评卷人三、计算题（每题15分，共60分）1. 已知函数211y x =+的一组数据：求分段线性插值函数，并计算()1.5f 的近似值.计算题1.答案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） 对于初始值()()0,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.答案3. 解()331f x x x=--，()130f=-<，()210f=>()233f x x'=-，()12f x x''=，()2240f=>，故取2x=作初始值迭代公式为()()3111112113133n n nn n nn nf x x xx x xf x x---------=-=-'-()312121()31nnxx--+-或，1,2,...n= 02x=，()3122311.88889321x⨯+==⨯-，()3222 1.8888911.879453 1.888891x⨯+==⨯-210.009440.0001x x-=>()3322 1.8794511.879393 1.879451x⨯+==⨯-，320.000060.0001x x-=<方程的根 1.87939x*≈4. 写出梯形公式和辛卜生公式，并用来分别计算积分111dxx+⎰.计算题4.答案4 解梯形公式()()()2bab af x dx f a f b-≈⎡+⎤⎣⎦⎰应用梯形公式得11111[]0.75 121011dxx≈+=+++⎰辛卜生公式为()()()[4()]62bab a a bf x dx f a f f b-+≈++⎰应用辛卜生公式得()() 111010[04()1] 162dx f f fx-+≈++ +⎰1111[4]16101112=+⨯++++25 36 =得 分 评卷人四、证明题（本题10分）确定下列求积公式中的待定系数，并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案证明：求积公式中含有三个待定系数，即101,,A A A -，将()21,,f x x x =分别代入求积公式，并令其左右相等，得1011123112()02()3A A A h h A A h A A h ---⎧⎪++=⎪--=⎨⎪⎪+=⎩得1113A A h -==，043hA =。
中南大学研究生入学考试数学分析试题
中南大学研究生入学考试数学分析试题中南大学 - 研究生考试数学分析试题一、求下列极限（1）lim ,(0)n nnnn x x x x x --→+∞->+； （2）1lim ()1xx x x →+∞+-；（3）01lim sin AA xdx A →∞⎰。
二、（共16分，每小题8分）设函数 ()sinf x xπ=，(0,1)x ∈（1）证明()f x 连续；（2）()f x 是否一致连续？（请说明理由）。
三、（共16分，每小题8分） （1）设ax by u e +=，求n 阶全微分n d u ；（2）设cos u x e θ=，sin u y e θ=，变换以下方程22220z zx y ∂∂+=∂∂。
四、（共20分，每小题10分）（1）求积分101ln 1dx x-⎰；（2）求曲面22az x y =+ (0)a >，和22z x y =+所围成的体积。
五、（共12分，每小题6分）设1cos 21p qn n n I nπ∞==+∑，(0)q > （1）求I 的条件收敛域； （2）求I 的绝对收敛域。
六、证明：积分2()0()x a F a e dx +∞--=⎰是参数a 的连续函数。
七、（8分）设定义于(,)-∞+∞上的函数()f x 存在三阶的导函数(3)()f x ，且(1)0f -=，(1)1f =，(1)(0)0f =证明：(3)(1,1)sup ()3x f x ∈-≥。
一、（共27分，每小题9分）求下列极限 （1）lim ()n n n n →+∞+-；（2）1220lim[3(cos )]xxxx t dt →+⎰；（3）设()f x 在[0,1]上可积，且1()1f x dx =⎰，求1121lim ()2n n k k f n n →+∞=-∑。
二、（共24分，每小题12分）设函数()f x 在[,)a +∞上连续， （1）证明：若lim ()x f x →+∞存在，则()f x 在[,)a +∞上一致连续；（2）上述逆命题是否成立？（请给出证明或举出反例）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2 2 −1
0
3 1 −2 , (2) A = 1 4 2 .
−2
1
2 1
10. 设矩阵A ∈ Rn×n 为Hessenberg形, 对QR变换 A = QR, B = QT AQ = RQ 证明矩阵Q 和B 都是Hessenberg形矩阵. 2
2=
1) 是A 的一个特征值及对应的特征
.
向量. 试证: 若有正交矩阵P 使得P x = e1 , 则有 P AP T = (b) 已知矩阵
λ
0
0 B
2 A= 10
10 5 −8
2
−8 .
2
11
的 一 个 特 征 值λ = 9和 对 应 的 特 征 向 量x = (2/3, 1/3, 2/3)T . 试 求 镜 面 反 射 矩 阵P 使得P x = e1 , 并计算P AP T . 7. 用正交相似变换将下列矩阵化为对称三对角矩阵:
1.0 1.0
0 1 3
0.5 0.25
0.25 .
0.5
2.0ห้องสมุดไป่ตู้
5. 设x = (1, 1, 1, 1)T , 用下列两种方法分别求正交矩阵P , 使得P x = ± (a) P 为平面旋转矩阵的乘积. (b) P 为镜面反射矩阵.
x
2
e1 .
1
6. (a) 设矩阵A ∈ Rn×n 为对称矩阵, λ 和x( x
习题七
1. 用幂法求下列矩阵的主特征值和主特征向量:
3 −2 −4 A= −2 6 −2 .
−4 −2
3
当特征值有3位小数稳定时迭代终止, 再对计算结果用Aitken外推加速. 2. 用反幂法求下列矩阵模最小的特征值和对应的特征向量:
3 −4 3 A= −4 6 3 .
3
3
1
3. 用反幂法求矩阵
6 2 1 A= 2 3 1 .
1 1 1
的最接近6的特征值和特征向量. 4. 用Jacobi方法计算下列矩阵的全部特征值和特征向量:
4 0 0 1.0 A1 = 0 3 1 ; A2 = 1.0
1 3 4 A= 3 1 2 .
4 2 1
8. 用镜面反射变换求下列矩阵的QR变换:
.
1 1 1 A= 2 −1 −1 2 −4
5
9. 用平面旋转变换对下列Hessenberg形矩阵作一步QR变换: