大连理工大学矩阵与数值分析2017年考题
大连理工大学《矩阵与数值分析》学习指导与课后参考答案第三章、逐次逼近法
第三章 逐次逼近法1.1内容提要1、一元迭代法x n+1=φ(x n )收敛条件为：1）映内性x ∈[a,b]，φ(x) ∈[a,b] 2）压缩性∣φ(x) -φ(y)∣≤L ∣x-y ∣其中L ＜1，此时φ为压缩算子，在不断的迭代中，就可以得到最终的不动点集。
由微分中值定理，如果∣φ’∣≤L ＜1，显然它一定满足压缩性条件。
2、多元迭代法x n+1=φ(x n )收敛条件为：1）映内性x n ∈Ω，φ(x n ) ∈Ω 2）压缩性ρ（▽φ）＜1，其中▽φ为x n 处的梯度矩阵，此时φ为压缩算子，在不断的迭代中，就可以得到最终的不动点集。
3、当φ(x )= Bx+f 时，收敛条件为，ρ（B ）＜1，此时x n+1= Bx n +f ，在不断的迭代中，就可以得到线性方程组的解。
4、线性方程组的迭代解法，先作矩阵变换 U L D A --= Jacobi 迭代公式的矩阵形式 f Bx b D x U L D x n n n +=++=--+111)(Gauss-Seidel 迭代公式的矩阵形式 f Bx b L D Ux L D x n n n +=-+-=--+111)()( 超松弛迭代法公式的矩阵形式f Bx b L D x U D L D x k k k +=-++--=--+ωωωωω111)(])1[()(三种迭代方法当1)(<B ρ时都收敛。
5、线性方程组的迭代解法，如果A 严格对角占优，则Jacob 法和Gauss-Seidel 法都收敛。
6、线性方程组的迭代解法，如果A 不可约对角占优，则Gauss-Seidel 法收敛。
7、Newton 迭代法，单根为二阶收敛 2211'''21lim)(2)(lim---∞→+∞→--=-==--k k k k k k k k x x x x f f c x x ξξαα8、Newton 法迭代时，遇到重根，迭代变成线性收敛，如果知道重数m ， )()('1k k k k x f x f m x x -=+仍为二阶收敛 9、弦割法)()())((111--+---=k k k k k k k x f x f x x x f x x 的收敛阶为1.618，分半法的收敛速度为（b-a ）/2n-110、Aitken 加速公式11211112)(),(),(+----+-+--+---+---===k k k k k k k k k k k x x x x x x x x x x x ϕϕ1.2 典型例题分析1、证明如果A 严格对角占优，则Jacob 法和Gauss-Seidel 法都收敛。
大连理工大学 矩阵与数值分析 第4章-4.2非线性方程的迭代解法
敛呢？不管非线性方程 f (x) = 0 的形式如何，总可以构造
ϕ(x) = x − k(x)
x
(k(x) ≠ 0)
（4-25）
作为方程（4-17）求解的迭代函数。 因为
ϕ′(x) = 1− k′(x) f (x) − k(x) f ′(x)
可知 | ϕ ′(x) | 在根 α 附近越小时 ,其局部收敛速度越快，
则迭代法 xk+1 = ϕ (xk ) 是 p 阶收敛。
练习1 取迭代函数
ϕ (x) = x + a(x2 −5)
要使迭代法收敛到 x* = 5, 则a应取何值？
且其收敛阶是多少？
解： ϕ′( x) = 1+ 2a x , 令
( ) ϕ′ 5 = 1+ 2a 5 < 1, 即有
−1 < 1+ 2a
x = x − f (x) = ϕ(x)
f ′(x)
（4-24）
建立的迭代格式至少是平方收敛。
证 根据定理4.6, 只需证明 ϕ′(α ) = 0 。 因为
ϕ′(α
)
=
⎡ ⎢
x
−
⎣
f f
(x) ′( x)
⎤ ⎥ ⎦
' x=α
=
⎡ ⎢1 ⎣
−
(
f
′(
x))2 − ( f ′(
f (x) x))2
f
′′( x)
x1 = 2× 0 −1 = −1 , x2 = 2(−1)3 −1 = −3 , x3 = 2(−3)3 −1 = −55 , L
显然, 当 k → ∞时, xk → −∞ ,故迭代法发散。 上述例子表明，迭代法的收敛与发散，依赖于迭代
大连理工大学2017年研究生矩阵与数值分析考试
大连理工大学2017年研究生矩阵与数值分析考试大连理工大学2017年研究生矩阵与数值分析考试考试日期：2017年6月5日一、填空题（50分，每空2分）1.a=0.3000经过四舍五入具有4位有效数字，则，2.已知X=(1,5,12)T，Y=(1,0,a)T，则由X映射到Y的Householder矩阵为：，计算||H||2=，cond2(H)=3.根据3次样条函数的性质（后面-前面=a（x-x0）3），一个求其中的参数b==4.，写出隐式Euler格式：梯形法格式：5.已知A=XXT，其中X为n维列向量，则||A||2=，||A||F=，矩阵序列的极限：=6.A=LU，其解为，写出一步迭代后的改善格式：7.，请问通过幂法与反幂法计算出的特征值分别是，8.，=，=，=，=，=9.是Newton-cotes公式，则=，具有代数精度=10.f(x)=7x7+6x6+…+x，f[20,21,22….,28]=11.，=12.f(0)=1,f(1)=-1,f(2)=1,f(3)=19,请问对该节点进行插值后最高次的系数=还有2空没有回忆出来，但是比上面题目还简单，因此不用担心。
二、，（1）计算LU分解（2）利用LU求逆矩阵（3）写出G-S格式（12分）三、给出，计算该迭代式收敛到某个值，收敛阶（8分）答案：收敛到，且收敛阶为3，因为，，而四、y=ae-bx，利用最小二乘法计算。
（8分）x-1012ye-11ee2数据可能有错，但是不影响计算思路。
五、计算权函数为1，区间[-1,1]的二次正交多项式，并且据此计算的具有三次代数精度求积公式（8分）六、已知线性2步3阶法（14分）（1）写出局部截断误差（必须含有主项）（2）判断收敛性（3）写出绝对稳定区间答：提示：上面公式的与书上的不是同一个，注意计算的时候区分。
大连理工大学(已有10试题)
大连理工大学应用数学系数学分析2001——2005，2009（2005有答案）高等代数2000——2005、2007（2005有答案）物理系数学物理方法2000——2005量子力学2000，2002——2005热力学与统计物理2000，2002——2005电动力学2000，2002——2005普通物理2000——2005光学（几何光学与波动光学）2000晶体管原理2000半导体材料2004——2005半导体器件2004——2005半导体物理2001——2002，2004——2005神经科学基础2004——2005生物统计学2004——2005生物物理学2004——2005工程光学2005微电子技术2003——2005离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005模拟电子技术2001——2005工程力学系材料力学1999——2001，2003——2005，2010（2010为回忆版）理论力学1995，1999——2001，2003——2005理论力学（土）2000土力学1999——2005自动控制原理(含现代20%) 1999——2005杆系结构静力学1998，2000弹性力学（不含板壳）1999——2004流体力学1999——2005流体力学（土）2004——2005流体力学基础2002——2005岩石力学1999——2000钢筋混凝土结构1999——2000工程流体力学2001，2004——2005水力学1999——2000，2002，2004——2005机械工程学院机械设计2001——2005（2001——2005有答案）机械原理1999——2000，2003——2005画法几何及机械制图2003——2005控制工程基础2001，2003——2005微机原理及应用（8086）1999——2000微机原理及应用（机）2004——2005微机接口与通讯及程序设计1999——2000模拟电子技术2001——2005离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005过程控制（含计算机控制）2000杆系结构静力学1998，2000微电子技术2003——2005系统工程概论1999——2002晶体管原理2000系统工程概论1999——2005管理基础知识1999——2001，2003——2005（2003——2005有答案）计算机组成原理（软）2005管理学基础2004——2005（2004——2005有答案）管理学2010（回忆版）材料力学1999——2001，2003——2005，2010（2010为回忆版）自动控制原理(含现代20%) 1999——2005材料科学与工程学院材料科学基础2003——2005，2010（2010为回忆版）机械设计2001——2005（2001——2005有答案）模拟电子技术2001——2005微电子技术2003——2005物理化学2004物理化学及物理化学实验1991——1993，2000，2002——2005（2002——2004有答案）胶凝材料学2001——2005硅酸盐物理化学2001——2002，2005杆系结构静力学1998，2000金属学2000金属热处理原理2000金属材料学2000钢筋混凝土结构1999——2000晶体管原理2000土木水利学院材料力学（土）2000，2003——2005材料力学1999——2001，2003——2005，2010（2010为回忆版）土力学1999——2005结构力学2000——2001，2003——2005水力学1999——2000，2002，2004——2005杆系结构静力学1998，2000理论力学（土）2000弹性力学（不含板壳）1999——2004流体力学1999——2005流体力学（土）2004——2005流体力学基础2002——2005岩石力学1999——2000钢筋混凝土结构1999——2000工程流体力学2001，2004——2005系统工程概论1999——2005工程经济学2004——2005无机化学2003——2005传热学2002，2004——2005工程力学2004——2005工程项目管理2004——2005建筑材料2005工程热力学2001——2002，2004——2005热工基础（含工程热力学和传热学）2003化工学院无机化学2003——2005物理化学2004物理化学及物理化学实验1991——1993，2000，2002——2005（2002——2004有答案）有机化学及实验2001，2003——2005高分子化学及物理2002——2005化工原理及化工原理实验2001——2005材料力学1999——2001，2003——2005，2010（2010为回忆版）工程流体力学2001，2004——2005硅酸盐物理化学2001——2002，2005热力学基础2005天然药物化学2005药剂学2005生物化学及生物化学实验1999——2005船舶工程学院船舶动力装置2002——2005船舶设计原理2001——2005水声学原理2002——2005船舶静力学2001——2005杆系结构静力学1998，2000电子与信息工程学院模拟电子技术2001——2005信号与系统（含随机信号20%）1999——2005 自动控制原理(含现代20%) 1999——2005工程光学2005通信原理2004——2005离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005离散数学与计算机组成原理2005离散数学与数据库原理2004——2005数据结构与计算机组成原理2004——2005计算机组成原理与计算机体系结构2004——2005 计算机组成原理与数字逻辑2000计算机组成原理（软）2005编译方法1999——2000操作系统1999——2001高等代数2000——2005过程控制（含计算机控制）2000微电子技术2003——2005微机接口与通讯及程序设计1999——2000系统工程概论1999——2005晶体管原理2000能源与动力学院汽车理论2000——2005机械原理1999——2000，2003——2005自动控制原理(含现代20%) 1999——2005化工原理及化工原理实验2001——2005普通物理2000高等代数2000——2005离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005运筹学基础及应用2004——2005计算机信息管理1999——2001，2004——2005 微电子技术2003——2005杆系结构静力学1998，2000系统工程概论1999——2005晶体管原理2000信息管理与信息系统2010（回忆版）管理学院计算机信息管理1999——2001，2004——2005 运筹学基础及应用2004——2005离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005公共经济学基础2004——2005，2010（2010为回忆版）过程控制（含计算机控制）2000微电子技术2003——2005系统工程概论1999——2002政治学原理2004——2005行政管理学2004——2005，2010（2010为回忆版）经济学基础2001——2005（2001——2005有答案）运筹学基础及应用2004——2005公共管理学2005社会保障学2004——2005管理学2010（回忆版）信息管理与信息系统2010（回忆版）人文社会科学学院经济学基础2001——2005（2001——2005有答案）管理基础知识1999——2001，2003——2005（2003——2005有答案）管理学基础2004——2005（2004——2005有答案）管理学2010（回忆版）系统工程概论1999——2002现代科学技术基础知识1999——2000，2004——2005思想政治教育学2004——2005马克思主义哲学原理2004——2005马克思主义哲学2001——2002西方哲学史2005哲学概论2004——2005科学技术史（含命题作文）2004——2005科学史、技术史、命题作文2001——2003政治学原理2004——2005行政管理学2004——2005，2010（2010为回忆版）传播学2004——2005新闻传播实务2004——2005民法学2004——2005法理学与商法总论2004——2005政治学2004——2005中外教育史2004——2005教育学2005中国近现代史2004——2005世界近现代史2004——2005电气工程及应用电子技术系电路理论2002——2005自动控制原理(含现代20%) 1999——2005过程控制（含计算机控制）2000微电子技术2003——2005系统工程概论1999——2005晶体管原理2000外国语学院二外德语2002，2004二外俄语2002——2004二外法语2004——2005二外日语2002——2004专业基础英语2003英汉翻译2003，2005英汉翻译与写作2004英语水平测试2004——2005二外英语2002——2005日语水平测试2004——2005翻译与写作（日）2004——2005专业基础日语2002——2003外国语言学与应用语言学（日语）专业综合能力测试2002——2003体育教学部运动生物力学2005人体测量与评价2004——2005生物学基础2005体质学2004——2005建筑艺术学院建筑设计（8小时）2000，2004——2005建筑设计原理1999——2000，2003建筑设计理论综合2004——2005城市建设史2002——2003中国与外国建筑史2000建筑构造与建筑结构1999——2000城市规划历史与理论2004——2005城市规划原理2003城市设计2002规划设计(8小时)2004-2005素描(8小时)2005泥塑(8小时)2005色彩(4小时)2005软件学院离散数学及应用2000——2001（2000有答案）离散数学与数据结构2002——2005离散数学与计算机组成原理2005离散数学与数据库原理2004——2005数据结构与计算机组成原理2004——2005计算机组成原理与计算机体系结构2004——2005计算机组成原理与数字逻辑2000计算机组成原理（软）2005编译方法1999——2000操作系统1999——2001环境与生命学院物理化学2004物理化学及物理化学实验1991——1993，2000，2002——2005（2002——2004有答案）化工原理及化工原理实验2001——2005硅酸盐物理化学2001——2002，2005基因工程原理2004——2005微生物学2004——2005细胞生物学2005环境化学2004——2005环境工程原理2004——2005，2010（2010为回忆版）分子遗传学2004——2005环境微生物2002经济系经济学基础2001——2005（2001——2005有答案）公共经济学基础2004——2005，2010（2010为回忆版）高科技研究院数学分析2001——2005，2009（2005有答案）高等代数2000——2005数学物理方法2000——2005量子力学2000，2002——2005热力学与统计物理2000，2002——2005电动力学2000，2002——2005物理化学2004物理化学及物理化学实验1991——1993，2000，2002——2005（2002——2004有答案）硅酸盐物理化学2001——2002，2005微电子技术2003——2005。
大连理工大学矩阵与数值分析上机作业
end
case2%2-范数
fori=1:n
s=s+x(i)^2;
end
s=sqrt(s);
caseinf%无穷-范数
s=max(abs(x));
end
计算向量x，y的范数
Test1.m
clearall;
clc;
n1=10;n2=100;n3=1000;
x1=1./[1:n1]';x2=1./[1:n2]';x3=1./[1:n3]';
xlabel('x');ylabel('p(x)');
运行结果：
x=2的邻域：
x =
1.6000 1.8000 2.0000 2.2000 2.4000
相应多项式p值：
p =
1.0e-003 *
-0.2621 -0.0005 0 0.0005 0.2621
p(x)在 [1.95,20.5]上的图像
程序：
[L,U]=LUDe.(A);%LU分解
xLU=U\(L\b)
disp('利用PLU分解方程组的解:');
[P,L,U] =PLUDe.(A);%PLU分解
xPLU=U\(L\(P\b))
%求解A的逆矩阵
disp('A的准确逆矩阵:');
InvA=inv(A)
InvAL=zeros(n);%利用LU分解求A的逆矩阵
0 0 0.5000 -0.2500 -0.1250 -0.0625 -0.0625
0 0 0 0.5000 -0.2500 -0.1250 -0.1250
0 0 0 0 0.5000 -0.2500 -0.2500
大连理工大学【精品】2016-2017学年第2学期高等数学A期末考试试卷
装订线大连理工大学2016-2017学年第2 学期高等数学A期末考试试卷2016~2017学年第2 学期考试科目：高等数学A考试类型：（闭卷）考试考试时间：120 分钟学号姓名年级专业题号一二三四总分得分评阅人一、填空题（本大题共5小题，每小题3分，共15分）1．二元函数2ln(21)z y x=-+的定义域为。
2. . 设向量设向量(2,1,2)a=，(4,1,10)b=-，c b al=-，且a c^，则l=。
3．经过(4,0,2)-和(5,1,7)且平行于x轴的平面方程为。
4．设yzu x=，则du=。
5．级数11(1)npnn¥=-å，当p满足条件时级数条件收敛。
二、单项选择题（本大题共5小题，每小题3分，共15分）1．微分方程2()'xy x y y+=的通解是（）A．2xy Ce=B．22xy Ce=C．22yy e Cx=D．2ye Cxy=2．求极限(,)(0,0)24limx yxyxy®-+=（）A．14B．12-C．14-D．12得分得分1 2 3．直线:327x y z L ==-和平面:32780x y z p -+-=的位置关系是的位置关系是 （ ）A ．直线L 平行于平面pB ．直线L 在平面p 上C ．直线L 垂直于平面pD ．直线L 与平面p 斜交斜交4．D 是闭区域2222{(,)|}x y a x y b £+£，则22Dx y d s +=òò （ ）A ．33()2b a p - B ．332()3b a p - C ．334()3b a p - D ．333()2b a p - 5．下列级数收敛的是下列级数收敛的是 （（ ））A ．11(1)(4)n n n ¥=++å B ．2111n nn¥=++å C ．1121n n ¥=-å D ．311(1)n n n ¥=+å三、计算题（本大题共7小题，每小题7分，共49分） 1. 求微分方程'xy y e +=满足初始条件0x =，2y =的特解。
大连理工大学数值分析历年真题与答案(研究生期末卷)
位
. ,
A 2＝
4 2 （3）设 A 2 4 , 则 A 1＝ 谱半径 ( A) ＝
,
A ＝
,
A F＝
, .
, 2-条件数 cond 2 ( A) ＝
, 奇异值为
线
（4）设 A C 44 ,特征值 1 2 2, 3 4 3 ,特征值 2 是半单的，而特征值 3 是 亏损的，则 A 的 Jordan 标准型 J
x 3 ( x [1,1]) 的二次最佳平方逼近多项式, 构造 Gauss 型求积公式 f ( x )dx A0 f ( x0 ) A1 f ( x1 ) , 并验证
1
1
其代数精度.
A－3
大
连
理 工
计算方法 数学系
大
学 2006 年试题
试卷： A 考试形式： 闭卷 试卷共 8 页
A－5
1 3 四、 （4 分）求 Householder 变换矩阵将向量 x 2 化为向量 y 0 . 2 0
五、 （12 分）写出解线性方程组的 Jacobi 法，G-S 法和超松弛（SOR）法的矩阵表示形式， 并根据迭代法 x ( k 1) Bx ( k ) f 对任意 x ( 0) 和 f 均收敛的充要条件为 ( B) 1 , 证明若线性方 程组 Ax b 中的 A 为严格对角占优矩阵, 则超松弛（SOR）法当松弛因子 (0,1] 时收敛.
师：张宏伟
一、填空（每一空 2 分，共 42 分） 1．为了减少运算次数，应将表达式.
装
16 x 5 17 x 4 18 x 3 14 x 2 13 x 1 x 4 16 x 2 8 x 1
大连理工大学2000-2017年数学分析真题
大连理工大学2000年数学分析真题 (2)大连理工大学2001年数学分析真题 (4)大连理工大学2002年数学分析真题 (6)大连理工大学2003年数学分析真题 (8)大连理工大学2004年数学分析真题 (10)大连理工大学2005年数学分析真题 (12)大连理工大学2006年数学分析真题 (14)大连理工大学2008年数学分析真题 (16)大连理工大学2009年数学分析真题 (18)大连理工大学2010年数学分析真题 (20)大连理工大学2011年数学分析真题 (22)大连理工大学2013年数学分析真题 (24)大连理工大学2014年数学分析真题 (25)大连理工大学2015年数学分析真题 (28)大连理工大学2016年数学分析真意 (30)大连理工大学2017年数学分析真题 (32)大连理工大学2000年数学分析真题一.从以下的第一到第八题中选取6题解答，每题10分 1.证明：()xx f 1=于区间()10，δ（其中0<0δ<1）一致连续，但是于（0,1）内不一致连续。
2.证明：若()x f 于[a ，b]单调，则()x f 于[a ，b]内Riemann 可积。
3.证明：Dirichlet 函数：()()⎪⎩⎪⎨⎧==有理数为无理数q px q x x f ,1,0在所有无理点连续，在有理点间断。
4.证明：若()()b a C x f ,∈，（指（a ，b ）上的连续函数，且任意()()b a ,,⊂βα，()⎰=βα0dx x f ，那么()()b a x x f ,0∈≡，。
5.证明：∑∞=-1n nx ne 于（0,+∞）不一致收敛，但是对于0>∀δ，于[)+∞,δ一致收敛。
6.证明：()⎪⎩⎪⎨⎧=≠=0,00,1sin 4x x xx x f ，在0=x 处有连续的二阶导数。
7.利用重积分计算三个半长轴分别为a,b,c 的椭球体的体积。
8.计算第二类曲面积分：⎰⎰∑++zdxdy ydzdx xdydz ，其中，∑是三角形()10,,=++>z y x z y x ，，法方向与z y x ,,轴成锐角为正。
大连理工大学矩阵与数值分析试卷-2013
1 0 0 0
3 ⎞ ⎛2 5 ⎟ T ⎟ ; LL 分解中 L= ⎜ ⎜3 4 ⎜ − ⎟ ⎟ ⎝2 5⎠
1 1 2 2
0 ⎞ ⎟ 7 ⎟。 ⎟ 2 ⎠
Gauss 求 积 公 式 ， 则
1 ∫ x + 1 f (x ) dx ≈ A f (x ) + A f (x ) + A f (x ) 为
2）为使二点数值求积公式 积节点和求积系数应为 (A) x0 = −
∫
1
f ( x) 1 − x2
.
−1
dx ≈ A0 f ( x0 ) + A1 f ( x1 ) 具有最高的代数精度，其求
B
2 2 π 1 1 1 , x1 = ; A0 = A1 = ； (B) x0 = − , x1 = ; A0 = A1 = ； 2 2 2 2 2 2
⎛ ⎜ 即 V = ( v1 v2 ) = ⎜ ⎜ ⎜ ⎝ ⎛ ⎜ V1 = V = ⎜ ⎜ ⎜ ⎝ 1 2 1 2 1 ⎞ ⎛ ⎟ ⎜ 2⎟ 或 V = ( v1 v2 ) = ⎜ −1 ⎟ ⎜ ⎟ ⎜ 2⎠ ⎝ 1 2 1 2 ⎛ ⎜ ⎜ ⎜ ⎜ ⎝ −1 ⎞ 2⎟ ⎟ ，因 rank(A)=1，故有 1 ⎟ ⎟ 2⎠ 1 ⎞ ⎛ 1 ⎞ ⎟ ⎜ ⎟ 2⎟ (1) = ⎜ 2 ⎟ , 由 U = (U1U 2 ) , 则 1 ⎟ ⎜ 1 ⎟ ⎟ ⎜ ⎟ 2⎠ ⎝ 2⎠
17）. 为了减少运算次数，应将表达式.
4 x3 − 3x 2 − 2 x − 1 改写为 x4 + x2 + x − 1
( ( 4 x − 3) x − 2 ) x − 1 ； ( ( ( x + 0 ) x + 1) x + 1) x − 1
大连理工大学 矩阵与数值分析 第5章5.2.4Hermite与分段低次2017(春)
=
sin
⎛ ⎜⎝
π
2
⎞ ⎟⎠
=
1
,
f
′ ⎛⎜⎝
π
2
⎞ ⎟⎠
=
cos
⎛ ⎜⎝
π
2
⎞ ⎟⎠
=
0
,
H3 (x) =
f
⎛ ⎜⎝
π
4
⎞ ⎟⎠
⎛ ⎜ ⎜1 ⎜
−
2
x
π
−π
4
−π
⎝ 42
⎞ ⎟ ⎟ ⎟
⋅
⎛ ⎜ ⎜ ⎜
x
π
−π
2
−π
⎠⎝4 2
⎞2 ⎟ ⎟ ⎟ ⎠
+
f
′
⎛ ⎜⎝
π
4
⎞ ⎟⎠
⎛ ⎜⎝
x
−
π
4
⎞ ⎟⎠
⎛ ⎜ ⋅⎜ ⎜
二点三次Hermite插值多项
设已知函数 f(x) 在 2 个互异点 x1, x2 处的函数值和导数值：
f ( x1 ) , f ′( x1 ) ; f ( x2 ) , f ′( x2 ) ;
构造一个三次多项式 H3 ( x) = ax3 + bx2 + cx + d
使其满足插值条件
H3 ( x1 ) = f ( x1 ) , H3′ ( x1 ) = f ′( x1 ) ;
答案是否定的。原因在于插值节点增多导致插值多项式 的次数增高，而高次多项式的振荡次数增多有可能使插值多 项式在非节点处的误差变得很大。
例如，对于函数
f
(x
)
=
1
5 +x
2
在[-5,5]上构造等距节
xk
数值分析试题及答案
数值分析试题及答案一、单项选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。
大连理工大学矩阵与数值分析上机作业13388
大连理工大学矩阵与数值分析上机作业课程名称：矩阵与数值分析研究生姓名：交作业日时间：2016 年12 月20日1.1程序：Clear all;n=input('请输入向量的长度n：')for i=1:n;v(i)=1/i;endY1=norm(v,1)Y2=norm(v,2)Y3=norm(v,inf)1.2结果n=10 Y1 =2.9290Y2 =1.2449Y3 =1n=100 Y1 =5.1874Y2 =1.2787Y3 =1n=1000 Y1 =7.4855Y2 =1.2822Y3 =1N=10000 Y1 =9.7876Y2 =1.2825Y3 =11.3 分析一范数逐渐递增，随着n的增加，范数的增加速度减小；二范数随着n的增加，逐渐趋于定值，无群范数都是1.2.1程序clear all;x(1)=-10^-15;dx=10^-18;L=2*10^3;for i=1:Ly1(i)=log(1+x(i))/x(i); d=1+x(i);if d == 1y2(i)=1;elsey2(i)=log(d)/(d-1);endx(i+1)=x(i)+dx;endx=x(1:length(x)-1);plot(x,y1,'r');hold onplot(x,y2);2.2 结果2.3 分析红色的曲线代表未考虑题中算法时的情况，如果考虑题中的算法则数值大小始终为1，这主要是由于大数加小数的原因。
第3题3.1 程序clear all;A=[1 -18 144 -672 2016 -4032 5376 -4608 2304 -512];x=1.95:0.005:2.05;for i=1:length(x);y1(i)=f(A,x(i));y2(i)=(x(i)-2)^9;endfigure(3);plot(x,y1);hold on;plot(x,y2,'r');F.m文件function y=f(A,x) y=A(1);for i=2:length(A); y=x*y+A(i); end;3.2 结果第4题4.1 程序clear all;n=input('请输入向量的长度n：')A=2*eye(n)-tril(ones(n,n),0);for i=1:nA(i,n)=1;endn=length(A);U=A;e=eye(n);for i=1:n-1[max_data,max_index]=max(abs(U(i:n,i))); e0=eye(n);max_index=max_index+i-1;U=e0*U;e1=eye(n);for j=i+1:ne1(j,i)=-U(j,i)/U(i,i);endU=e1*U;P{i}=e0;%把变换矩阵存到P中L{i}=e1;e=e1*e0*e;endfor k=1:n-2Ldot{k}=L{k};for i=k+1:n-1Ldot{k}=P{i}*Ldot{k}*P{i};endendLdot{n-1}=L{n-1};LL=eye(n);PP=eye(n);for i=1:n-1PP=P{i}*PP;LL=Ldot{i}*LL;endb=ones(n,2);b=e*b; %解方程x=zeros(n,1);x(n)=b(n)/U(n,n);for i=n-1:-1:1x(i)=(b(i)-U(i,:)*x)/U(i,i);endX=U^-1*e^-1*eye(n);%计算逆矩阵AN=X';result2{n-4,1}=AN;result1{n-4,1}=x;fprintf('%d:\n',n)fprintf('%d ',AN);4.2 结果n=51.0625 -0.875 -0.75 -0.5 -0.06250.0625 1.125 -0.75 -0.5 -0.06250.0625 0.125 1.25 -0.5 -0.06250.0625 0.125 0.25 1.5 -0.0625-0.0625 -0.125 -0.25 -0.5 0.0625n=101.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 1.125 -0.75 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 1.25 -0.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 0.0625 0.125 0.25 1.5 -0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625 -0.0625 -0.125 -0.25 -0.5 0.0625同样的方法可以算出n=20,n=30时的结果，这里就不罗列了。
大连理工数学分析试题及解答
大连理工数学分析试题及解答大连理工大学硕士生入学考试数学分析试题一. 从以下的1到8题中选答6题1. 证明:2()f x x =在区间[0,]M 内一致连续(M 为任意正数),但是在[0,)+∞不一致连续2. 证明:若()f x 在[,]a b 内连续,那么()f x 在[,]a b 内Riemann 可积.3. 证明:若1α>,那么广义积分1sin x dx α+∞收敛4. 证明:若()f x ,()g x 为区间(,)a b 上的连续函数,对任意的(,)(,)a b αβ?有:()()f x dx g x dx ββαα=??,那么, ()()f x g x ≡于(,)a b5. 证明:若1nn a∞=∑收敛,那么1nxn n a e∞-=∑在[0,)+∞一致收敛6. 已知：2,0()0,0x e x f x x -?≠?=?=??，求"(0)f7. 已知：()()1(,)()22x atx at x at x at u x t d aφφψαα+-++-=+. 其中, ψ和φ分别是可以求导一次和求导两次的已知函数,计算22222(,)(,)u x t u x t a t x ??-??8. 计算,半径为R 的球的表面积二. 从9到14题中选取6题9.已知: lim '()0x f x →∞=,求证: ()lim0x f x x→∞=10.证明: ()af x dx +∞收敛,且lim ()x f x λ→+∞=,那么0λ=11.计算曲面积分: 333SI x dydz y dzdx z dxdy =++??, 其中S 为旋转椭球面2222221x y z a b c++=的外侧12.设()[0,1]f x C ∈,(0)0f =,(1)1f =,0()1f x ≤<. 求证: ()()n n S x f x =对于任意小于1的正数δ,在区间(0,1]δ-一致收敛,但是不在(0,1)一致收敛13.设()[0,1]f x C ∈,(0)0f =,(1)1f =,0()1f x ≤<. 求证: 1lim ()0n n f x dx →∞=?14.证明:若()[,]n u x C a b ∈，1,2,...,...n =且1()n n u b ∞=∑发散,那么1()n n u x ∞=∑不在[,)a b 一致收敛大连理工大学2001年硕士生入学考试数学分析试题解答一.1. 证利用定义证明(1) 对于0ε?>,21M εδ?=+,12||x x δ?-<,那么12121212|()()||()()|2||2f x f x x x x x M x x M δε-=-+<-<<(2) 任取1ε=,0δ?>,1211,22x x δδδ==+, 1212121|()()||()()|1f x f x x x x x δδ-=-+>=,推出矛盾,从而命题得证■2. 证利用一致连续的定义和Riemann 可积的定义来做因为函数在闭区间内连续,所以一致连续. 根据一致连续的定义对0ε?>,δ?,12||x x δ?-<,12|()()|f x f x ε-<考虑可积的定义,对于一个[,]a b 分割112:...n a a a a b ?=<<<=,11max ||i i i na a λ+≤<=-下面证明:振幅函数121110,[,]1()limmax {()}()i i n i i x x a a i w x f x a a λ+-+→∈==-∑=0当λδ<时,12111110,[,]110()limmax {()}()()i i n n i i i i x x a a i i w x f x a a a a b λεε+--++→∈==≤=-≤-=∑∑.根据夹逼定理,不难得到()0w x =. 从而,命题得证■3. 证利用莱布尼兹交错级数:假设;n a n π=,1sin nn a n a s x dx α-=?考虑:111|||||sin ||sin |n nnn a a n n a a s s x dx x dx αα+-+-=-?1111[|sin ||sin |]n n n n x x dx xx dx ππππααππα--+-=+??1111[|sin |(2)|sin |]n n n n xx dx n x x dx ππππααππαπ--++=--??1111[(2)]|sin |0n n x n x x dx ππααπαπ--+=--<?11lim |||sin |||lim ||0nnn n a a n n a a n n s x dx dx n n s αααπππ--→∞→∞=≤=--?=??如此,不难看出1sin x dx α+∞是一个莱布尼兹交错级数,从而命题得证■4. 证不妨设:2a bc +=()()x c F x f t dt =?,那么()()F x G x =于(,)x a b ∈因为()f x ()g x 都是(,)x a b ∈上的连续函数,所以()'()()()f x F x G xg x ===■5. 证利用A-D 判别法做，也可以通过Abel 求和公式出发推导1nxn n a e∞-=∑中nxn b e-=,现在,根据原题:1n n a∞=∑收敛,1nxnb e -=≤一致有界所以,根据Abel 判别法,知该函数项级数在定义域一致收敛. ■6. 解题目有问题,在零点不连续■7. 解不断利用链式求导法则()()1(,)()22x atx at x at x at u x t d aφφψαα+-++-=+(,)()()()()()()()()()()22'()'()'()'()22u x t xx at x at x at x at x at x at x at x at x at x x at x x x ax at x at x at x at aφφψψφφψψ+?+?-?-?+?-++--?+??-=+++-+--=+22'()'()()()(,)()()()()22"()"()'()'()22x at x at x at x at u x t x at x at x at x at x ax at x at x at x at aφφψψφφψψ?+?-?+?-+-+?-?+?-=+++-+--=+同理:(,)()()()()()()()()()()22'()'()'()'()22u x t tx at x at x at x at x at x at x at x at x at t x at t t t aa x at a x at x at x at φφψψφφψψ+?+?-?-?+?-++--?+??-=++--++-=+222'()'()()()(,)()()()()"()"()'()'()22x at x at x at x at aa u x t x at x at x at x at x x at x at x at x at a aφφψψφφψψ?+?-?+?--++?-?+?-=+++-+--=+22222(,)(,)0u x t u x t a t x ??-=??■8. 解方法很多,此处介绍一种比较简单的假设:()V R 为半径R 为的球的体积2234()()3R R V R R x dx R ππ-=-=?假设: ()S R 为半径R 为的球的表面积20()()()'()4RV R S x dx S R V R R π=?==?■二9. 证L ’Hosptial 法则因为x →+∞，()'()lim lim lim '()0'x x x f x f x f x x x →∞→∞→∞===■10. 证反证法如果命题不成立,即0λ≠,那么,根据极限的定义,G ?,当x G >的时候, |()|||2f x λ>()Gf x dx +∞→∞?和收敛矛盾,从而命题得证■11. 解利用Gauss 定理加换元3332223()VSI x dydz y dzdx z dxdy x y z dxdydz =++=++换元sin cos sin sin ,[0,1],[0,2),[0,]cos x ar y br r z cr ?θ?θθπ?π?=??=∈∈∈??=?4222222223sin (sin sin sin cos cos )VI abc r a b c drd d ??θ?θ?θ?=++22322322200033(sin sin sin cos )cos sin 55abc abc a b d d c d πππ?θ?θ?θ=++332232203646()sin ()5555abc abc abc abc a b d a b πππ??=++=++?■12. 证首先由于在闭区间内连续,所以函数在闭区间内一致连续(1)(0,1]x δ?∈-,根据确界存在定理,存在上确界,且上确界不等于1,否则和题意矛盾不妨设:(0,1]sup ()1x f x m δ∈-=<根据定义,对于0ε?>,ln ln N mε=,当n N >,|()||()|n n n S x f x m ε=≤< 从而知一致收敛于0(2)首先,根据前半题,显然()n S x 于(0,1)x ?∈收敛于0由于(1)1f =,且函数一致收敛，存在一组数列：12...a a <<,1()1i fa n=- 如此,考虑11lim ()lim ()lim(1)0nnn n n n n n S a f a ne→∞→∞→∞==-=≠,从而不是一致收敛的. ■13. 证利用前一小题的结论因为()nf x 内闭一致收敛,对于0ε?>,2εδ?=,当n 足够大的时候:10()2n f x dx δε-<又1111|()|||2n f x dx dx δδε--<=所以,1111()()()n n n f x dx f x dx f x dx δδε--=+<?从而命题得证. ■14. 证反证法:假设命题不成立,那么1()n n u x ∞=∑在[,)a b 一致收敛.即0ε?>,N ?,,m n N ?>,(,)x a b ?∈,|()|m n nu x ε<∑因为|()|lim |()|m mn n x bnnu b u b ε→=≤∑∑,否则与()[,]n u x C a b ∈矛盾而1|()|n n u b ∞=∑发散,所以|()|n n Nu b ∞=∑发散,与|()|lim |()|m m n n x bnnu b u b ε→=≤∑∑矛盾从而命题得证. ■。
大连理工大学矩阵与数值分析上机作业13388
共享知识分享快乐大连理工大学矩阵与数值分析上机作业课程名称：矩阵与数值分析研究生姓名：12 交作业日时间：日20 月年2016卑微如蝼蚁、坚强似大象．共享知识分享快乐第1题1.1程序：Clear ;all n=input('请输入向量的长度n：') for i=1:n;v(i)=1/i;endY1=norm(v,1)Y2=norm(v,2)Y3=norm(v,inf)1.2结果n=10 Y1 =2.9290Y2 =1.2449Y3 =1n=100 Y1 =5.1874Y2 =1.2787Y3 =1n=1000 Y1 =7.4855Y2 =1.2822Y3 =1N=10000 Y1 =9.7876Y2 =1.2825Y3 =11.3 分析一范数逐渐递增，随着n的增加，范数的增加速度减小；二范数随着n的增加，逐渐趋于定值，无群范数都是1.卑微如蝼蚁、坚强似大象．共享知识分享快乐第2题2.1程序;clear all x(1)=-10^-15;dx=10^-18;L=2*10^3; i=1:L fory1(i)=log(1+x(i))/x(i); d=1+x(i); d == 1ify2(i)=1;elsey2(i)=log(d)/(d-1);endx(i+1)=x(i)+dx;end x=x(1:length(x)-1););'r'plot(x,y1,on holdplot(x,y2);卑微如蝼蚁、坚强似大象．共享知识分享快乐2.2 结果2.3 分析红色的曲线代表未考虑题中算法时的情况，如果考虑题中的算法则数值大小始终为1，这主要是由于大数加小数的原因。
第3题3.1 程序;clear all A=[1 -18 144 -672 2016 -4032 5376 -4608 2304 -512];x=1.95:0.005:2.05; i=1:length(x);for y1(i)=f(A,x(i)); y2(i)=(x(i)-2)^9;end figure(3);plot(x,y1);;on hold卑微如蝼蚁、坚强似大象．共享知识分享快乐);'r'plot(x,y2,F.m文件y=f(A,x)function y=A(1); i=2:length(A);for y=x*y+A(i);;end3.2 结果第4题卑微如蝼蚁、坚强似大象．共享知识分享快乐4.1 程序;clear all n=input('请输入向量的长度n：')A=2*eye(n)-tril(ones(n,n),0); i=1:n for A(i,n)=1;end n=length(A);U=A; e=eye(n);for i=1:n-1[max_data,max_index]=max(abs(U(i:n,i))); e0=eye(n);max_index=max_index+i-1; U=e0*U; e1=eye(n); j=i+1:n fore1(j,i)=-U(j,i)/U(i,i);endU=e1*U;中把变换矩阵存到P P{i}=e0;% L{i}=e1; e=e1*e0*e;endk=1:n-2for Ldot{k}=L{k}; i=k+1:n-1forLdot{k}=P{i}*Ldot{k}*P{i};endend Ldot{n-1}=L{n-1};LL=eye(n);PP=eye(n); i=1:n-1for PP=P{i}*PP;LL=Ldot{i}*LL;endb=ones(n,2);解方程 %b=e*b;x=zeros(n,1);x(n)=b(n)/U(n,n); i=n-1:-1:1for卑微如蝼蚁、坚强似大象．共享知识分享快乐x(i)=(b(i)-U(i,:)*x)/U(i,i);end计算逆矩阵%X=U^-1*e^-1*eye(n);AN=X'; result2{n-4,1}=AN;result1{n-4,1}=x;,n)'%d:\n'fprintf(fprintf('%d ',AN);4.2 结果n=51.0625 -0.875 -0.75 -0.5 -0.0625-0.0625 0.0625 -0.75 1.125 -0.5-0.0625 0.125 0.0625 1.25 -0.5-0.0625 0.1250.25 0.06251.50.0625-0.5-0.25-0.0625 -0.125n=101.0625 -0.875 -0.75 -0.5 -0.0625 1.0625 -0.875 -0.75 -0.5 -0.0625 -0.0625 1.125 0.0625 -0.75 -0.5 -0.5 0.0625 1.125 -0.75 -0.0625 -0.0625 0.0625 0.125 1.25 1.25 -0.0625 -0.5 0.0625 0.125 -0.5-0.0625 0.250.250.0625 0.1251.5 1.5 -0.0625 0.1250.06250.0625 -0.0625 -0.125 -0.25 0.0625 -0.5 -0.0625 -0.125 -0.25 -0.5 -0.0625 -0.75 1.0625 -0.5 -0.0625 -0.875 -0.5 -0.75 1.0625 -0.875 -0.0625 -0.5 0.0625 1.125 -0.5 0.0625 1.125 -0.75 -0.0625 -0.75 1.25 0.125 0.0625 -0.0625 -0.0625 -0.5 -0.5 0.0625 0.125 1.250.25-0.0625 -0.0625 1.50.1250.0625 0.0625 0.250.1251.5-0.0625 -0.125 -0.25 0.0625-0.5 0.0625 -0.0625 -0.125 -0.25-0.5同样的方法可以算出n=20,n=30时的结果，这里就不罗列了。
大连理工大学程名松矩阵与数值分析计算方法-第3章
1 (1 − ρ ( A ) ) > 0 一定存在 2 一种相容的矩阵范数 ⋅ ，使得 A ≤ ρ ( A) + ε 。
充分性 根据定理2.8，对于 ε =
又根据相容矩阵范数的性质， 再注意到上述关系式有 1 ρ ( A ) + ε = (1 + ρ ( A ) ) < 1 2 那么
A k ≤ A ≤ ( ρ ( A) + ε ) ≤ q k < 1
k→∞
k→∞
1⎞ lim ⎛ 1 + ⎟ k →∞ ⎜ ⎝ k⎠
k
lim sin k k →∞ k lim
lim
k →∞
k →∞
e− k
k
k
⎞ ⎛e 0 ⎞ ⎟ ⎜1 0 ⎟= A ⎟ =⎜ ⎟ ⎟ ⎜1 1 ⎟ ⎠ ⎠ ⎝
由矩阵序列极限的定义可以看出，矩阵序列收敛的性质和数 列收敛性质相似。 由定义可见，C m × n 中的矩阵序列的收敛相当于mn个数列同时 收敛。 因此可以用初等分析的方法来研究它。
∞ ∞
并且 则
lim Ak = A ， k →∞
lim Bk = B
k →∞
lim A k B k = AB
k →∞
证 由
Ak Bk − AB = Ak Bk − Ak B + Ak B − AB
≤ B ⋅ Ak − A + Ak ⋅ Bk − B
由定理1和推论可知，结论成立。
性质3
n n lim Ak = A 并且 设 {Ak }k=1∈C × 中的矩阵序列，
∞
⎛ 0.1 0.7⎞ 。 由于 A ∞ = 0.9 < 1 ，故 计算 ∑ A ，其中 A = ⎜ ⎟ k =0 ⎝ 0.3 0.6⎠ ∞ k ρ ( A ) < 1，从而 ∑ A 收敛，且有
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
大连理工大学2017年研究生矩阵与数值分析考试
考试日期：2017年6月5日
一、填空题（50分，每空2分）
1.a=0.3000经过四舍五入具有4位有效数字，则x a a
-≤ ，ln ln x a -≤ 2.已知X=(1,5,12)T ，Y=(1,0,a)T ，则由X 映射到Y 的Householder 矩阵为： ，计算||H||2= ，cond 2(H)=
3.根据3次样条函数的性质（后面-前面=a （x -x0）3），一个求其中的参数b==
4. 2'3u u t =，写出隐式Euler 格式： 梯形法格式：
5.已知A=XX T ，其中X 为n 维列向量，则||A||2= ，||A||F = ，矩阵序列的极限：2lim k
k A A →∞⎛⎫ ⎪ ⎪⎝⎭
= 6.A=LU ，其解为x ，写出一步迭代后的改善格式： 7. 531A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭
，请问通过幂法与反幂法计算出的特征值分别是 ，
8. 1111A ⎛⎫ ⎪= ⎪ ⎪⎝⎭
，sin A = ，823A A A +-= ，At e = ，d d At e t = ，2
1
At e dt ⎰=
9. ()()()()2
0120
012f x dx A f A f A f =++⎰是Newton -cotes 公式，则1A = ，具有代数精度=
10. f(x)=7x 7+6x 6+…+x ，f[20, 21,22…., 28]=
11. 0.40.200.5A ⎛⎫= ⎪⎝⎭，1
k k A ∞=∑= 12.f(0)=1, f(1)=-1, f(2)=1, f(3)=19,请问对该节点进行插值后最高次的系数= 还有2空没有回忆出来，但是比上面题目还简单，因此不用担心。
二、121232352A -⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭，121b ⎛⎫ ⎪= ⎪ ⎪-⎝⎭
（1）计算LU 分解
（2）利用LU 求逆矩阵
（3）写出G -S 格式（12分）
三、给出()
21233k k k k x x a x x a ++=+，计算该迭代式收敛到某个值，收敛阶（8分）
3，
因为()1k x x φ+=
，
0φφ'''==
，而302a φ'''=-≠
-bx 五、计算权函数为1，区间[-1,1]的二次正交多项式，并且据此计算20
(x)dx f ⎰的具有三次代数精度求积公式（8分）
六、已知线性2步3阶法（14分） ()21102210812
n n n n n n h u u u f f f ααββ++++++=++ （1）写出局部截断误差（必须含有主项）
（2）判断收敛性
（3）写出绝对稳定区间
答：提示：上面公式的2β与书上的2β不是同一个，注意计算的时候区分。