研究生数值分析试题

数值分析（研究生，2008-12-15）1.（10分）求函数⎩⎨⎧≤≤++<≤-+=10,101,1sin )(2x x x x x x f 在区间[-1，1]上的最佳平方逼近式x e a x a a x 210)(++=φ。

2．（15分）利用乘幂法计算下列矩阵的主特征值和相应的特征向量⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----110141012，初始向量为T x ]0,0,1[0=（要求结果有三位有效数字）。

同时计算该矩阵的1-条件数和谱条件数。

3.（15分）已知函数x x f sin )(=在36.0,34.0,32.0210===x x x 处的值分别为352274.0,333487.0,314567.0210===y y y 。

用Lagrange 插值多项式对3167.0=x 的函数值进行近似计算，并估计近似计算的误差界。

4.（15分）用Newton 迭代法求方程0ln 2=+x x 在区间（0，2π）内的解，选择你认为合适的初始点，计算方程的根，使得近似解具有四位有效数字。

请从理论上估计达到所需精度所需的迭代次数。

5.（15分）用Gauss-Seidel 迭代法解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---542834*********x x x 取初始近似向量0[0,0,0]Tx =，估计达到4位有效数字需要的迭代次数，并实际计算之。

就该具体问题分析计算过程中总的乘除法计算量。

6. （10分）应用拟牛顿法解非线性方程组⎪⎩⎪⎨⎧=-+=-+.12,2322112221x x x x x x 取T x ]1,0[)0(= ，终止容限210-=ε。

7.（10分） 求解矛盾方程组⎪⎪⎩⎪⎪⎨⎧=++=++=++=++232328.12221321321321321x x x x x x x x x x x x8. （10分）用复合Simpson 公式计算积分⎰=21sin )(xdx f I 讨论在误差要求不超过410-的条件下的步长。

------------------------------------------------ 装 ---------------------------------订 ---------------------------------线 ------------------------------------------------装 订 线 左 侧 不 要 书 写 内 容允许使用计算器一、 填空题 （本大题共10小题，每小题 2分，共 20分）1. 若2.71828x e == ，取近似值* 2.7180x =，则*x 具有 4 位有效数字。

2.为了提高数值计算精度，应将8格式进行计算。

3.已知n=3时牛顿—柯特斯系数(3)(3)(3)012133,,888C C C ===，那么(3)3C =18 。

4.设3()1f x x x =+-，则函数的四阶差商[0,1,2,3,4]f = 0 。

5. 用牛顿迭代法解方程0x x e --=在0.5x =附近的近似实根的牛顿迭代格式为)1,0(e 1e )()(1=+--='-=--+n x x x f x f x x nnx x n n n n n n6. 对给定的剖分01:n a x x x b ∆=<<<= ，当()s x 满足条件 ()s x 在[a,b]有2阶连续导数且在每个子区间上是个3次多项式 时是三次样条函数。

7．用最小二乘法拟合三点()()()0,1,1,3,2,2A B C 的直线是1322y x =+。

8.向量序列()211cos ,sin ,3Tk k x e k k k k -⎛⎫=+ ⎪⎝⎭ 的极限向量为()0,1,3T9.求积公式 10311()()(1)434f x dx f f ≈+⎰的代数精度为 2 。

10.若绝对误差限为31102-⨯，那么近似数0.03600有 2 位有效数字二、单项选择题（本大题共5小题，每小题 2 分，共 10分）1. 已知实验数据555521111(,)(1,2,3,4,5),15,31,55,105.5,k k k k kk k k k k k x y k x y x x y =========∑∑∑∑其中则用最小二乘法求近似公式01y a a x =+的法方程为（ C ）A 0101153155105.5a a a a +=⎧⎨+=⎩B 0101515551531105.5a a a a +=⎧⎨+=⎩C 0101515311555105.5a a a a +=⎧⎨+=⎩ D0101531153155105.5a a a a +=⎧⎨+=⎩ 2． 以下矩阵是严格对角占优矩阵的是（ B ）A 3210141011410012⎛⎫ ⎪ ⎪ ⎪⎪⎝⎭ B 2100131013610113-⎛⎫⎪--⎪ ⎪-- ⎪-⎝⎭C 5210113121410012-⎛⎫⎪--⎪ ⎪⎪⎝⎭D 4211141021411315⎛⎫⎪ ⎪⎪- ⎪⎝⎭3.已知两种递推公式11(1)35(1,2,,20)31(2)(20,,1)55n n n n I nI n I I n n n--=-==-= 则在数值计算过程中（ C ）。

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

二、填空题（4*5=20）1、4； 2、⎪⎪⎪⎭⎫ ⎝⎛323203*⎪⎪⎪⎪⎭⎫⎝⎛320323； 3、)]23()0()23([3f f f ++-∏；4、kk 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(41)(-=x x x P 。

.。

（2分） 四、（10分）设⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--=1001001*10010021321u 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 ,,11516)(>=J B ρ即Jacobi 迭代发散；。

允许使用计算器一、 填空题 （本大题共10小题，每小题 2分，共 20分） 1. 若 2.71828x e ==，取近似值* 2.7180x =，则*x 具有 4 位有效数字。

2.为了提高数值计算精度，应将8格式进行计算。

3.已知n=3时牛顿—柯特斯系数(3)(3)(3)012133,,888C C C ===，那么(3)3C =18 。

4.设3()1f x x x =+-，则函数的四阶差商[0,1,2,3,4]f = 0 。

5. 用牛顿迭代法解方程0xx e在0.5x 附近的近似实根的牛顿迭代格式为)1,0(e 1e )()(1=+--='-=--+n x x x f x f x x nnx x n n n n n n6. 对给定的剖分01:n a x x x b ∆=<<<=，当()s x 满足条件 ()s x 在[a,b]有2阶连续导数且在每个子区间上是个3次多项式 时是三次样条函数。

7．用最小二乘法拟合三点()()()0,1,1,3,2,2A B C 的直线是1322y x =+。

8.向量序列()211cos ,sin ,3Tk k xe k k k k -⎛⎫=+ ⎪⎝⎭的极限向量为()0,1,3T9.求积公式 10311()()(1)434f x dx f f ≈+⎰的代数精度为 2 。

10.若绝对误差限为31102-⨯，那么近似数有 2 位有效数字二、单项选择题（本大题共5小题，每小题 2 分，共 10分）1. 已知实验数据555521111(,)(1,2,3,4,5),15,31,55,105.5,k k k k kk k k k k k x y k x y x x y =========∑∑∑∑其中则用最小二乘法求近似公式01y a a x =+的法方程为（ C ）A 0101153155105.5a a a a +=⎧⎨+=⎩B 0101515551531105.5a a a a +=⎧⎨+=⎩C 0101515311555105.5a a a a +=⎧⎨+=⎩D 0101531153155105.5a a a a +=⎧⎨+=⎩2． 以下矩阵是严格对角占优矩阵的是（ B ）A 3210141011410012⎛⎫ ⎪ ⎪ ⎪⎪⎝⎭ B 2100131013610113-⎛⎫⎪--⎪ ⎪-- ⎪-⎝⎭C 5210113121410012-⎛⎫ ⎪--⎪ ⎪ ⎪⎝⎭ D 4211141021411315⎛⎫⎪ ⎪⎪- ⎪⎝⎭3.已知两种递推公式11(1)35(1,2,,20)31(2)(20,,1)55n n n n I nI n I I n n n--=-==-=则在数值计算过程中（ C ）。

研究⽣考试数值分析试题研究⽣2002级数值分析⼀（12分）、对于积分=+1,2,1,0,999n dx x x n。

（1）试推导递推公式 ,2,1,19991=+-=-n nI I n n ；（2）分析上述算法的数值稳定性；（3）若上⾯算法不稳定，请选择合适的算法，并分析其稳定性。

⼆（12分）、解⽅程组= 00001.8800001.626221x x 和?=00002.8800001.626221x x ，就所观察到的现象进⾏分析。

三（12分）、设⽅程组=--=+-=+-7989783212121x x x x x x x ；（1）适当调整⽅程的排列顺序，使得⽤Gauss-Seidel 迭代法求解时收敛？说明收敛原因。

（2）取初始向量()()Tx 0,0,00=，⽤Gauss-Seidel 迭代求近似解()2x，并求其()()k k x x-+1误差。

四（12分）、（1）已知函数()4xe xf =，在[0,1]内三点0，1/2，1的函数值，求其⼆次插值的余项；（2）三个节点如何安排能使其余项达最⼩，此时⼈余项为多少？五（12分）、对于⽅程()02ln =+-x x ，若求［-1.9，-1］内的根，分别选取迭代⽅程()2ln +=x x 和2-=x e x ，它们的收敛性如何？再写出⽜顿迭代公式。

六（10分）、设()?=>+-='100,5y x x y y ，解析解xe x y -+-=25262515，分别取45.0,4.0,2.0,1.0=h ，利⽤Euler ⽅法计算得y(10)的近似值分别为1.96，1.96，5.2851，142.8863，对此现象进⾏分析。

七（10分）、设()x e x f =，分别取步长0001.0,01.0,5.0=h ，⽤中⼼差商公式计算()0f '的近似值并求出误差，对结果作分析⽐较。

⼋（10分）、求不超过2次的多项式()x P 2，使其满⾜条件：()21=f ，()32=f ，()12='f ，并写出其误差估计。

数值分析（100分试题） 第 1 页 共 3 页一、、填空题（本大题共5小题，每小题4分，共20分）1、设*0.034x 为经过四舍五入后得到的近似数，则数*x 的有效数字位数是 。

2、设节点,0,1,2,3,,i x i n = ，(),0,1,2,,i l x i n = 是关于上述节点的Lagrange 插值基函数，则对于0,1,2,,k n = ，0()n k i i i x l x ==∑ 。

3、已知矩阵411141114A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦， 求||||A ∞= ；2()cond A = 。

4、给定方程22cos x x =-,求该方程根的Newton 迭代格式是 。

5、 步长为h 时，求常微分方程初值问题⎩⎨⎧=≤≤=-1)0(,10,0'3y x xy y 的改进的Euler 公式是 。

二、（10分）求一个3次多项式)(x p ,使其满足4)2('',3)2(,2)1(',1)1(====p p p p .三、（10分）给定线性方程组12310112013a x a a x a x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦其中a 为常数.写出求解上述线性方程组的Jacobi 迭代格式，并分析当a 取何值时Jacobi 迭代法收敛。

四、（10分）用列主元Gauss 消去法解线性方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡81213612002120321203214321x x x x 五、（10分）求3()f x x = 在区间[1,1]- 上关于()1x ρ= 的最佳平方逼近2次多项式。

六、（10分） 分析方程01224=---x x x 存在几个实根,并用迭代法求出其中一个实根,精确到3位有效数字。

七、（10分）已知求积公式 )53(95)0(98)53(95)(11f f f dx x f ++-≈⎰- 为Gauss 公式,试给出形如)()()()(221100x f A x f A x f A dx x f ba++≈⎰的求积公式,使其代数精度达到5.八、（10分）用初等反射矩阵将111211245A ⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦分解为QR 的形式，其中Q 为正交矩阵，R 为上三角矩阵。

武汉大学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学年第一学期硕士研究生期末考试试题科目名称：数值分析学生所在院：学号：姓名：注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。

数值分析试题及答案一、单项选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。

一.(6分) 已知描述某实际问题的数学模型为xy y x y x u 223),(+=，其中，y x ,由统计方法得到，分别为4,2==y x ,统计方法的误差限为0.01,试求出u 的误差限)(u ε和相对误差限)(u r ε.解：)(23)(6)(),()(),()(222y x y x x x y xy y y y x u x x y x u u εεεεε⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-=∂∂+∂∂≈6.016.044.001.0)412(01.0)448(=+=⨯++⨯-=0.010714566.03)()(22=≈+=xy y x u u r εε 二.(6分) 已知函数13)(3+=x x f 计算函数)(x f 的2阶均差]2,1,0[f ,和4阶均差]4,3,2,1,0[f .解：21142512)1()2(]2,1[,311401)0()1(]1,0[=-=--==-=--=f f f f f f9232102]1,0[]2,1[]2,1,0[=-=--=f f f0!4)(]4,3,2,1,0[)4(==ξf f三.(6分)试确定求积公式: )]1(')0('[121)]1()0([21)(10f f f f dx x f -++≈⎰的代数精度. 解：记⎰=10)(dx x f I )]1(')0('[121)]1()0([21f f f f I n -++=1)(=x f 时：1110==⎰dx I 1]00[121]2[21=-+=n Ix x f =)(时：2110==⎰xdx I 21]11[121]1[21=-+=n I2)(x x f =时：31102==⎰dx x I 31]20[121]1[21=-+=n I3)(x x f =时：41103==⎰dx x I 41]30[121]1[21=-+=n I4)(x x f =时：51104==⎰dx x I 61]40[121]1[21=-+=n I求积公式)]1(')0('[121)]1()0([21)(1f f f f dx x f -++≈⎰具有3次代数精度. 四.(12分) 已知函数122)(23-++=x x x x f 定义在区间[-1,1]上,在空间},,1{)(2x x Span x =Φ上求函数)(x f 的最佳平方逼近多项式.其中,权函数1)(=x ρ,154))(),((,1532))(),((,34))(),((210-==-=x x f x x f x x f ϕϕϕ. 解：0))(),(())(),((21))(),((1101101100=====⎰⎰--dx x x x x x dx x x ϕϕϕϕϕϕ32))(),(())(),(())(),((112110220====⎰-dx x x x x x x x ϕϕϕϕϕϕ0))(),(())(),((1131221===⎰-dx x x x x x ϕϕϕϕ52))(),((11422==⎰-dx x x x ϕϕ解方程组⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛1541532345203203203202210a a a 得⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛15161210a a a 则)(x f 的最佳平方逼近多项式为：1516)(2-+=x x x p 五.(16分) 设函数)(x f 满足表中条件:(1) 填写均差计算表((2) 分别求出满足条件22k k k k 的 2次 Lagrange 和 Newton 差值多项式.(3) 求出一个四次插值多项式)(4x H ,使其满足表中所有条件.并用多项式降幂形式表示.解：12)12)(02()1)(0()20)(10()2)(1()(22+-=----+----=x x x x x x x L12)1)(0(1)0)(1(1)(22+-=--+--+=x x x x x x N 令)2)(1()(12)(24--+++-=x x x b ax x x x H则)2()()2)(1)(()2)(1(22)('4-++--++--+-=x x b ax x x b ax x x ax x x H)1()(-++x x b ax由 ⎩⎨⎧-=+=+⇒⎩⎨⎧=-++-=-=-++-=1220)12(2)2(24)2('2)21)((22)1('44b a b a b a H b a H 解得 5,3=-=b a因此1820143)2)(1()53(12)(23424++-+-=--+-++-=x x x x x x x x x x x H 六.(16分)(1). 用Romberg 方法计算⎰31dx x ,将计算结果填入下表(*号处不填).(2). 试确定三点 Gauss-Legender 求积公式⎰∑-=≈11)()(k k k x f A dx x f 的Gauss 点k x 与系数k A ,并用三点 Gauss-Legender 求积公式计算积分: ⎰31dx x .解：过点（1，-1）和点（3，1）作直线得 y t x +=所以积分⎰⎰-+=11312dt t dx x由三次Legendre 多项式 )35(21)(33x x x p -= 得得Gauss 点：,515,0,515210==-=x x x再由代数精度得 ⎪⎪⎪⎩⎪⎪⎪⎨⎧==+==+-==++⎰⎰⎰---32535305155152111220112011210dt x A A dt x A A dt A A A即 ⎪⎩⎪⎨⎧=+=-=++9/10022020210A A A A A A A解得 ,95,98,95210===A A A所以三点Gauss-Legendre 求积公式为：()⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛-≈⎰-5159509851595)(11f f f dx x f 因此 79746.2515295298515295211=+++-≈+=⎰-dx t I 七.(14分)(1) 证明方程02ln =--x x 在区间(1,∞)有一个单根.并大致估计单根的取值范围. (2) 写出Newton 迭代公式,并计算此单根的近似值.(要求精度满足: 5110||-+<-k k x x ). 解：令 2ln )(--=x x x f),1(,011)('∞∈>-=x xx f > 即)(x f 在区间 ),1(∞ 单调增 又 04)(,02ln )2(22>-=<-=e e f f 所以 02ln =--x x 在区间 ),1(∞有一单根 ),1(20e x ∈Newton 迭代公式为1ln 112ln 1-+=----=+k kk k kk k k k x x x x x x x x x 令 20=x 计算得八. (12分) 用追赶法求解方程组:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛022112111131124321x x x x 的解.解： 由计算公式 ⎪⎩⎪⎨⎧-===+====-1,,2,,,2,,111111n i c n i b a c b i i ii i i i i i βααβγγβαα得 ,2,1,1,21,1,24321111======γγγββαα25211322212=⨯-=⇒=+ααβγb 52222222==⇒=αββαc c 53521133323=⨯-=⇒=+ααβγb 35333333==⇒=αββαc c 37352144434-=⨯-=⇒=+ααβγb因此 ⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛135152121137253125121211113112 即 LU A = 令 b Ly = 解 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-022137253125124321y y y y 得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛23753214321y y y y令 y Ux =解 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛237532113515212114321x x x x 得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫⎝⎛21104321x x x x九. (12分) 设求解初值问题⎩⎨⎧==00)(),('y x y y x f y 的计算格式为:)],(),([111--+++=n n n n n n y x bf y x af h y y ,假设11)(,)(--==n n n n y x y y x y ,试确定参数b a ,的值,使该计算格式的局部截断误差为二阶,即截断部分为: )(3h o .（注：原题中)(2h o 错误）解：)],(),([111--+++=n n n n n n y x bf y x af h y y )](')('[)(1-++=n n n x by x ay h x y])('''21)('')('[)(')(2++-++=n n n n n x y h x hy x y hb x hay x y ++-++=)('''21)('')(')()(32n n n n x by h x by h x y b a h x y 对比 ++++=+)('''61)(''21)(')()(321n n n n n x y h x y h x hy x y x y 得 ⎩⎨⎧==+2/11b b a ， 即 2/1==b a 时该计算格式具有二阶精度.。

2009级研究生《数值分析》试卷一.(6分) 已知描述某实际问题的数学模型为xyy x y x u 223),(+=，其中，y x ,由统计方法得到，分别为4,2==y x ,统计方法的误差限为0.01,试求出u 的误差限)(u ε和相对误差限)(u r ε.二.(6分) 已知函数13)(3+=x x f 计算函数)(x f 的2阶均差]2,1,0[f ,和4阶均差]4,3,2,1,0[f .三.(6分)试确定求积公式: )]1(')0('[121)]1()0([21)(10f f f f dx x f -++≈⎰的代数精度.四.(12分) 已知函数122)(23-++=x x x x f 定义在区间[-1,1]上,在空间},,1{)(2x x Span x =Φ上求函数)(x f 的最佳平方逼近多项式.其中,权函数1)(=x ρ,154))(),((,1532))(),((,34))(),((210-==-=x x f x x f x x f ϕϕϕ.五.(16分) 设函数)(x f 满足表中条件:(1) 填写均差计算表(标有*号处不填):(2) 分别求出满足条件)2,1,0(),()(),()(22===k x f x N x f x L k k k k 的 2次 Lagrange 和 Newton 差值多项式.(3) 求出一个四次插值多项式)(4x H ,使其满足表中所有条件.并用多项式降幂形式表示. 六.(16分)(1). 用Romberg 方法计算⎰31dx x ,将计算结果填入下表(*号处不填).(2). 试确定三点 Gauss-Legender 求积公式⎰∑-=≈112)()(k k k x f A dx x f 的Gauss 点k x 与系数k A ,并用三点 Gauss-Legender 求积公式计算积分: ⎰31dx x .七.(14分)(1) 证明方程02ln =--x x 在区间(1,∞)有一个单根.并大致估计单根的取值范围. (2) 写出Newton 迭代公式,并计算此单根的近似值.(要求精度满足: 5110||-+<-k k x x ).八. (12分) 用追赶法求解方程组:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛022112111131124321x x x x 的解.九. (12分) 设求解初值问题⎩⎨⎧==00)(),('y x y y x f y 的计算格式为:)],(),([111--+++=n n n n n n y x bf y x af h y y ,假设11)(,)(--==n n n n y x y y x y ,试确定参数b a ,的值,使该计算格式的局部截断误差为二阶,即截断部分为: )(3h o .2008年春季学期数值数学试题一．（10分）设给实数0a >，初值00x >：⑴试建立求1a的Newton 迭代公式，要求在迭代函数中不含除法运算；⑵证明给定初值0x ，迭代收敛的充分必要条件为020x a<<；⑶该迭代的收敛速度是多少？⑷取00.1x =，计算15的近似值，要求计算迭代三次的值（结果保留5位小数）。

研究生2002级数值分析一（12分）、对于积分⎰=+1,2,1,0,999n dx x x n。

（1）试推导递推公式 ,2,1,19991=+-=-n nI I n n ；（2）分析上述算法的数值稳定性；（3）若上面算法不稳定，请选择合适的算法，并分析其稳定性。

二（12分）、解方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡00001.8800001.626221x x 和⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡00002.8800001.626221x x ，就所观察到的现象进行分析。

三（12分）、设方程组⎪⎩⎪⎨⎧=--=+-=+-7989783212121x x x x x x x ；（1）适当调整方程的排列顺序，使得用Gauss-Seidel 迭代法求解时收敛？说明收敛原因。

（2）取初始向量()()Tx0,0,00=，用Gauss-Seidel 迭代求近似解()2x ，并求其()()k k x x -+1误差。

四（12分）、（1）已知函数()4xe xf =，在[0,1]内三点0，1/2，1的函数值，求其二次插值的余项；（2）三个节点如何安排能使其余项达最小，此时人余项为多少？五（12分）、对于方程()02ln =+-x x ，若求［-1.9，-1］内的根，分别选取迭代方程()2ln +=x x 和2-=x e x ，它们的收敛性如何？再写出牛顿迭代公式。

六（10分）、设()⎩⎨⎧=>+-='100,5y x x y y ，解析解x e x y -+-=25262515，分别取45.0,4.0,2.0,1.0=h ，利用Euler 方法计算得y(10)的近似值分别为1.96，1.96，5.2851，142.8863，对此现象进行分析。

七（10分）、设()xe xf =，分别取步长0001.0,01.0,5.0=h ，用中心差商公式计算()0f '的近似值并求出误差，对结果作分析比较。

数值分析试题一．（10分）设给实数0a >，初值00x >：⑴试建立求1a的Newton 迭代公式，要求在迭代函数中不含除法运算； ⑵证明给定初值0x ，迭代收敛的充分必要条件为020x a<<；⑶该迭代的收敛速度是多少？⑷取00.1x =，计算15的近似值，要求计算迭代三次的值（结果保留5位小数）。

二．（10分）试确定参数,,a b c ，使得下面分段多项式函数()s x 是三次样条函数。

332,01()1(1)(1)(1),132x x s x x a x b x c x ⎧≤≤⎪=⎨--+-+-+≤≤⎪⎩ ()s x 是否是自然样条函数？三．（10分）利用Dollite 三角分解方法求解方程组123121022331302x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦ 四．（10分）给定3阶线性方程组123122*********x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦讨论其Jacobi 迭代格式的收敛性五．（10分）推导出中矩形求积公式()()()2baa b f x dx b a f +≈-⎰ ，并求出其截断误差。

六．（10分用最小二乘法确定拟合公式bx y ae =中的参数,a b 。

七．（10分）根据已知函数表：建立不超过三次的Newton 插值项式。

八．（10分）试确定常数01,A A ,使求积公式1011()(f x dx A f A f -≈+⎰有尽可能高的代数精度,并指出代数精度是多少,该公式是否是Gauss 型?并用此公式计算积分311I dx x=⎰（结果保留5位小数）。

九．（10分）利用经典四阶Runge-Kutta 方法求初值问题：20,01(0)1y y x y '=-≤≤⎧⎨=⎩在0.2x =处的数值解（取步长0.1h =）。

10．（10分）讨论两步方法 11112(4)33n n n ny y y hy +-+'=-+ 的局部截断误差，求出它的局部阶段误差的首项（主部），它是多少阶的？ （在线性多步法的局部截断误差中10111[()()],2,3,!p prr r i i i i C i a r i b r r -==-⎧⎫=--+-=⎨⎬⎩⎭∑∑ ）。