数值分析第7章试题
数值分析课后习题与解答
课后习题解答第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式()有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式()(1.2.3)则得有5位有效数字，其误差限，相对误差限有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？〔1〕〔2〕解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。
〔1〕〔2〕4.近似数x*=0.0310,是 3 位有数数字。
5.计算取，利用：式计算误差最小。
四个选项：第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1与n=2的Lagrange插值或Newton插值,并应用误差估计〔5.8〕。
线性插值时，用0.5与0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值误差限，故2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h 应取多少?解：用误差估计式〔5.8〕，令因得3. 若，求和.解：由均差与导数关系于是4. 若互异，求的值，这里p≤n+1.解：，由均差对称性可知当有而当P＝n＋1时于是得5. 求证.解：解：只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048与cos 0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式〔5.17〕得其中计算时用Newton后插公式〔 5.18)误差估计由公式〔5.19〕得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
李庆扬-数值分析第五版第7章习题答案(0824)汇编
第7章复习与思考题求f (X )= 0的零点就等价于求(x )的不动点，选择一个初始近似值X 0，将它代入X =「(X ) 的右端，可求得X 1 h%X °)，如此反复迭代有 X k 1 二(X k ), k =0,1,2,..., (X)称为迭代函数，如果对任何X 。
• [a,b]，由x k 卜h%x k ),k =0,1,2,...得到的序列〈X k 1有极限则称迭代方程收敛，且X* =®(x*)为®(X )的不动点 故称X k q 二(X k ), k =0,1,2,...为不动点迭代法。
5•什么是迭代法的收敛阶？如何衡量迭代法收敛的快慢？如何确定X k 1 二「(X k )(k =0,1,2,...)的收敛阶P219设迭代过程X k 1'h%X k )收敛于 (X)的根X*，如果当k > 时，迭代误差e k = x k - x *满足渐近关系式—t C,C =const 式 0 e/则称该迭代过程是 p 阶收敛的，特别点，当 p=1时称为线性收敛，P>1时称为超线性收敛， p=2时称为平方收敛。
以收敛阶的大小衡量收敛速度的快慢。
6•什么是求解f(x)=0的牛顿法？它是否总是收敛的？若 f(X*) =0，X*是单根，f 是光 滑，证明牛顿法是局部二阶收敛的。
牛顿法:当| f (X k )卜J 时收敛。
7•什么是弦截法？试从收敛阶及每步迭代计算量与牛顿法比较其差别。
在牛顿法的基础上使用 2点的的斜率代替一点的倒数求法。
就是弦截法。
收敛阶弦截法1.618小于牛顿法2 计算量弦截法 <牛顿法(减少了倒数的计算量)8•什么是解方程的抛物线法？在求多项式全部零点中是否优于牛顿法？ P229X-mX k 1 =X kf (X k ) f (X k )设已知方程f （x） = 0的三个近似根，X k，X k^，X k^2，以这三点为节点构造二次插值多项式p（x），并适当选取p2（x）的一个零点X k卅作为新近似根，这样确定的迭代过程称为抛物线法。
第7 章第 2 题作业答案 用乘子法求解下列问题
w(3) = max{0, w(2) − σ [(x1(2) − 2)2 − x2(2) ]} = 2.46659 v(3) = v(2) − σ (2x1(2) − x2(2) −1) = 1.52622 第 3 次迭代：令σ = 2, w(3) = 2.46659 v(3) = 1.52622。
+
1
2σ
{[max(0,
w
−
σ
(( x1
−
2)2
−
x2
)]2
−
w2}
−
v(2x1
−
x2
− 1)
+
σ
2
(2 x1
−
x2
− 1)2
第 7 章第 2 题作业答案 分两种情形：
φ(x, w,v,σ ) =
⎪⎧( x1
−
2)2
+
( x2
−
3)2
+
1
2σ
{[w
−
σ
(( x1
−
2)2
−
x2
)]2
−
w2}
⎪ ⎪ ⎪
第 7 章第 2 题作业答案
而且满足条件
( x1(1)
−
2)2
−
x2(1)
=
−0.31267
<
w(1)
σ
=1
对 w, v 进行修正：
w(2) = max{0, w(1) − σ [(x1(1) − 2)2 − x2(1) ]} = 2.62534 v(2) = v(1) − σ (2x1(1) − x2(1) −1) = 1.35641 第 2 次迭代：令σ = 2, w(2) = 2.62534, v(2) = 1.35641。
数值分析第二版(丁丽娟)答案
1.32
1.68
2.08
2.52
3.00
解答下列问题 （1）试列出相应的差分表； （2）写出牛顿向前插值公式； （3）用二次牛顿前插公式计算 f（0.225）；
例3已知当 x=-1，0，2，3时，对应的函数值为
，
，
，
，
，求 的四次 Newton 插值多项式。
例4 设 对 n=1,2,3时
，证明：
例5 设 （1）
6 2730.5000 5051.0000 5051.5000
7 10922.5000 23483.0000 23483.5000
8 43690.5000 80827.0000 80827.5000
21.000000000000000 17.000000000000000 16.238095238095237 16.058823529411764 16.014662756598241 16.003663003663004 16.000915583226515
第一章答案
第二章答案
第三章答案
0 0.5 0.5 1 1 2.5000
5.0000 5.5000
第四章答案
2 10.5000 19.0000 19.5000
3 42.5000 91.0000 91.5000
4 170.5000 315.0000 315.5000
5 682.5000 1467.0000 1467.5000
3、 用规范化幂法求
按模最大的特征值和对应的特征向量，取初值
。当特征值有3位小数稳定时停止。
4、 用反幂法求矩阵
练习五
，迭代7次。
的最接近于6 的特征值和对应的特征向量，取初值
完整版数值分析第7章答案
1数值分析第七章第七章非线性方程求根一、重点内容提要（一）问题简介求单变量函数方程f(x)?0(7.1)f(x*)?0x*x*x*为也称为方程的根是指求(7.1).（实数或复数），使得称的根,m f(x)?(x?x*)g(x)f(x)f(x)函数的零点.若可以分解为g(x)g(x)?0x*x*为单称m=1满足时,是方程(7.1)的根.,则当其中m为正整数,g(x)x*x*是方程(7.1)的m称,充分光滑,为m重根.若重根,则有根;当m>1时(m?1)(m)f(x*)?f'(x*)?...?f(x*)?0,f(x*)?0f(x)f(a)f(b)?0,则方程(7.1)在(a,b)[a,b]若上连续且内至少有一个实根,称在[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得.(二)方程求根的几种常用方法1.二分法f(x)f(a)f(b)?0f(x)?0f(x)?0*x在上连续,再设内有根,则设.在(a,b)在[a,b]1x?(a?b)a?a,b?bf(x)f(x)?0000计算和.,若则(a,b)内仅有一个根.令20000a?xb?b[a,b])f(a)f(x?0x*?x;,则令,结束计算;若若得新的有根区间,10,11001a?ab?x0)?(f(a)fx,得新,则令的有根区间0110,0011b?a?(b?a)x?(a?b)[a,b][a,b]?[a,b]f(x)0101111再令计算,.,.同上法得221110101[a,b],如此反复进行出新的有根区间,可得一有根区间套22...?[a,b]?[a,b]?...?[a,b]001?n1?nnn2数值分析第七章11a?x*?b,n?0,1,2,...,b?a?(b?a)?...?(b?a)0n0?1nnn?1nn且. 221lim(b?a)?0,lim x?lim(a?b)?x* nnnnn故2????n??nn1x?(a?b)f(x)?0nnn的近似根,可作为,且有误差估计因此21(b?a)|x?x*|?n1?n(7.2)22.迭代法?(x?)x等价变形为将方程式(7.1) (7.3)??(x*)?)(xf(x*)?0x**xx*的一个不动点为函数.;反之亦然则.若要求称满足?(x)的不动点由式(7.3)产生的不动点迭代关系式(也求方程(7.1)的根等价于求称简单迭代法)为?(x),k?0,1,2...x?(7.4)k1?k?(x),k??x0,1,2...?(x)称为迭代函数.函数如果对任意,由式(7.4)产生的序列??x有极限kk??k则称不动点迭代法(7.4)收敛.kk?1x?x*lim?(x)?C[a,b]满足以下两个条件: 定理7.1(不动点存在性定理)设?(x)??b;x?[a,b]a有1.对任意??(y)|?|x?y|?,y[a,b]|(x)?x 2.存在正常数使对任意, ,都有(7.5)1?L?(x)[a,b]x*.则在上存在惟一的不动点?(x)?C[a,b]满足定理7.2(定理不动点迭代法的全局收敛性定理)设7.1中的两个??x]b,?x[a?(x)并条件,由,(7.4),的不动点式得到的迭代序列则对任意到.收敛k0有误差估计式3数值分析第七章L|x?*|?x||x?x1kkk?(7.6)L1?k L|x?x*|?|x?x|1?kkk L1?(7.7)和??'(xx))(xx**的某,为设在的不动点定理7.3(不动点迭代法的局部收敛性定理)?'(x)|?|1,则迭代法(7.4)局部收敛个邻域连续,且.?(xx?)x*,的根如果迭代误差收敛阶的概念设迭代过程(7.4)收敛于方程e?x?x*k??时成产下列渐近关系式当kk e k?1?C(常数C?0)e(7.8) k则称该迭代过程是p阶收敛的.特别地,p=1时称线性收敛,p>1时称超线性收敛,p=2时称平方收敛.(K)?(x)x*的邻近连续,并定理7.4(收敛阶定理在所求根)对于迭代过程(7.4),如果且(p?1)???(x*)?...?*)?'(x*)?0''(x(p)?(x*)?0(7.9)*x的邻近是收敛的,则该迭代过程在点并有e1)(p?1k?*)x?lim(p!ep??k (7.10)k斯蒂芬森(Steffensen)迭代法当不动点迭代法(7.4)只有线性收敛阶,甚至于不收敛时,可用斯蒂芬森迭代法进行加速.具体公式为??(y?)(x),zy?kkkk2)?x(y kk x?x?kk?1z?2y?x kkk k?0,1,2,...(7.11)4数值分析第七章此法也可写成如下不动点迭代式?(x),kx??0,1,2,...kk?12?)?x(x)(?(x)?x????(x)?2?(x(x))(7.12)?(x)x**x是为式(7.12)中则的不动点7.5(定理斯蒂芬森迭代收敛定理)设,?(x)???1*)''(x)?'(x(x)*x的不动点,存在,的不动点;设则,则斯蒂芬森迭代法是(7.11)是2阶收敛的.3.牛顿迭代法牛顿迭代法是一种特殊的不动点迭代法,其计算公式为f(x)k,x?k?0,1,2,...?x k?k1)xf'(其迭代函数为(7.13)k f(x)??(x)?x f'(x)f(x*)?0,f'(x*)?0,f''(x*)?0时牛顿迭代法的收敛速度当,容易证f''(x*)??0*)?''(x 0'(x*)?ff'(x*),由定理,明,7.4知,牛顿迭代法是平方收敛的,且ef''(x*)1?k?lim2*)f'(ex2??k(7.14)k f(x)?0(m?2)*x时,迭代函数的m重顿重根情形的牛迭代法当根是f(x)1??x)?(x?'(x*)?1??0?'(x*)|?1|)xf'(*x.所以牛顿迭代法求处的导数在,且m x*的重数m知道,重根只是线性收敛.若则迭代式f(x)k,k?0,1,2,...??xx?m kk?1)'(xf(7.15)k f(x)??x()f'(x)*x此时迭代式,的单重零点一定是函数,未知时m当.求重根二阶收敛5数值分析第七章?(x)f(x)f'(x)kkk?xx??x?kk?1k?)f''(x)x)]?f(x'(x)[f'(kkkk k?0,1,2,...(7.16)也是二阶收敛的.f(x)k,?k?0,1,2,...x?x k1k?)xf'(如下迭代法简化牛顿法0称为简化牛顿法或平行弦法.牛顿下山法为防止迭代不收敛,可采用牛顿下山法.具体方法见教材.4.弦截法f'(x)xxf(x)在,处的一阶差商来代替,将牛顿迭代法(7.13)中的即可得弦用kkk?1截法f(x)k(xx?x??x)1kk?1k?k f(x)?f(x)(7.17)??x*|:|x??*x内具有二阶连续导数,的邻域在其零点定理7.6假设且对任1kk?)(xfx,x??10f'(x)?0?x?,又初值,,意则当邻域充分小时,有弦截法(7.17)将按阶?1?5?p?1.6182???1?0?*x2的正根收敛到是方程..这里p5.抛物线法(x,f(x)),(x?f(x))两点的直线方程的根近似替弦截法可以理解为用过kk?1kk?1xxx0x)?(fx)?0f(用,过三若的根.已知个近似根,的2kk?1k?(x,f(x)),(x,f(x)),(x,f(x))f(x)?0的根,的抛物线方程的根近似代替2??k?k121k?kkk所得的迭代法称为抛物线法,也称密勒(Muller)法.f(x)f'(x*)?0*x,则抛物线法局部收敛当,在,的邻近有三阶连续导数且收敛阶p?1.839?1.84. 为数值分析第七章二、知识结构图三、常考题型及典型题精解3上有一个实根x*,并用二分法2]在[1,?1?例7-1 证明方程x0?x-6-3,需二分区间[1,2]10.若要求|x-x*|?求这个根,要求|x-x*|?10kk多少次?3在[1,2],则f(1)=-1<0,f(2)=5>0,故方程f(x)=0x?解设f(x)=x1?2在[1,2]时,f'(x)>0,即f(x)=0-1,所以当x?上有根x*.又因f'(x)=3x上有惟一实根x*.用二分法计算结果如表7-1所示.[1,2]7-1表k abxf(x)的符号kkkk+ 2 0 1 1.5- 1.5 1 1 1.25+ 2 1.25 1.51.3751.3125 3 1.251.375 -1.375 1.3438 1.3125 4 +1.312551.3282+1.1341.3125-861.32041.32041.32827-1.32431.32431.32821.3263+87数值分析第七章9 1.3243 1.3282 1.3253 +1.32631-3-3,可以作为x*的近??10此时x=1.3253满足|x-x*|?10?0.97799102似值.1-6?6,只需|x10-x*|?-x*|即可,解得k+1?19.932, 若要求|x?10?kkk+12即只需把[1,2]二分20次就能满足精度要求.x=1,(1)确定有根区间[a,b];(2)构造不动e例7-2 已知函数方程(x-2)点迭代公式使之对任意初始近似x?[a,b],迭代方法均收敛;(3)用所构0?3.|?10造的公式计算根的近似值,要求|x?x1k k?xx因此区间[2,3]0,e解 (1)令f(x)=(x-2)-1>-1,由于f(2)=-1<0,f(3)=e x x)=-1,f(,lim,lim f(x)=+?是方程f(x)=0的一个有根区间.又因f'(x)=(x-1)e???xx???1-1<0,当x>1时f(x)单增,x<1时f(x)单减,故f(x)=0在(-?,+?)内f'(1)=-e有且仅有一根x*,即x*?[2,3].x?xx?.由于当?将(x-2)e[2,3].则=1等价变形为x=2+ee(x)=2+,x(2)2??x??<1'(x)|=|-e?e[2,3]x?时2?|(x)?3,|x?[2,3]均收敛.??故不动点迭代法x=2+e x,k=0,1,2,...,对k0k+1x?进行迭代计算,结果如表7-2所示.e(3)取x=2.5,利用x=2+k k+10表7-28数值分析第七章此时x已满足误差要求,即x*?x?2.120094976.44例7?3考虑求解方程2cos x?3x?12?0的迭代公式2 x=4+cos x,k=0,1,2,...k k+13(1)试证:对任意初始近似x?R,该方法收敛;0-3;10-x|?(2)取x=4,求根的近似值x,要求|x k0k+1k+1(3)所给方法的收敛阶是多少?2?(x)=4+cos x,解 (1)由迭代公式知,迭代函数322?(x)的值域介于(4-)与(4+由于)之间,且(??,??).x?3322?'(x)|=|-sin x|??1|33?(x)在(??,??)内存在惟一的故根据定理7.1,7.2知,??收敛于x*.x?x?R,迭代公式得到的序列不动点x*,且对k0(2) 取x=4,迭代计算结果如表7-3所示.0表7-3x*?xx?3.347529903已满足误差要求，即此时55?'(x*)?0.136323129?0，故根据定理7 .4）由于（3知方法是线性收敛的，并e?1k?'(x?*)lim e??k。
计算方法第7章／《数值分析》／清华大学／上海交通大学／西安交通大学
Euler-Maclaurin 公式
å ò Tn - I =
k
c2 j [ f (2 j-1) (b) - f (2 j-1) (a)]h2 j +
b a
P2k ( x) f (2k ) (x)dx
j =1
8
å c2 j
=
(-1)
j +1
1
(2p )2
j
¥ k =1
1 k2j
，|
P2k (x) |£ c2k h2k
同上
ò RS =
b f (4) (x ) (x - a)(x - c)2 (x - b)dx a 4!
=- b - a (b - a )4 f (4) (h) 180 2
5
复化公式及误差分析
由上述误差表达式可知，区间越小，绝对误差越小，复化梯形公式：
将积分区间
n
等分，节点是
xi
=
a
+ ih, h
=
值公式
pn (xk ) = f (xk ) = yk
利用 Lagrange 插
1
å Õ pn (x) =
nn
(
k=0 j=0 j¹k
x xk
- xj - xj
)yk
¬ 做代换 x = a + th,t
=
x-a
å Õ n n t - j
=
(
k =0
j=0
k
-
j )yk
h
j¹k
以 pn (x) 代 f (x) 得
k -1
<e
停止
输出 Tk(k )
»
I
。否则 h
Ü
h 2
数值分析-第7章 矩阵特征值问题的数值解法n
7
9 11 12
6.104716
6.026349 6.006637 6.003327
(-0.450275, -0.322058, 1.0)
(-0.445914, -0.318617, 1.0) (-0.444814, -0.31775, 1.0) (-0.444630, -0.317606, 1.0)
其中i为A的特征值，P的各列为相应于i的特征向量。
P -1 AP D
2
n
2
定理7.1.3 ARnn，1, …, n为A的特征值，则
（1）A的迹数等于特征值之和，即 tr ( A) aii i
i 1 i 1
n
n
（2）A的行列式值等于全体特征值之积，即
1 xi(k +1) / xi(k )
i 1,2,, n
可见，当k充分大时， ( k ) 近似于主特征值， ( k +1) 与x ( k )的对应非零分量的比值 x x 近似于主特征值。
在实际计算中需要对计算结果进行规 , 范化。因为当 1 1时，x (k ) 趋于零， 当1 1时, x ( k )的非零分量趋于无穷。 从而计算时会出现下溢 或上溢。
特征值的范围. 解 我们先分别求出各个圆盘区域。 D1 = {z：|z – 1|£0.6}；D2 = {z：|z – 3|£0.8} D3 = {z：|z + 1|£1.8}；D4 = {z：|z + 4|£0.6}. 易见D2和D4为 弧立圆盘分别 包含A的两个实 特征值.
数值分析习题集及答案[1](精)
](https://img.taocdn.com/s3/m/7eabbe90cc22bcd126ff0c65.png)
数值分析习题集（适合课程《数值方法A 》和《数值方法B 》）长沙理工大学第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式(3.3)求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少?6. 设028,Y =按递推公式1n n Y Y -= ( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字≈27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2?10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b c s a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n nn n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj jj x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少? 9. 若2nn y =,求4n y ∆及4n y δ. 10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆. 12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()b aS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =. 3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式. 4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式. 5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nTx 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hhf x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10x e dx-⎰并计算误差. 5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2baf f x dx b a f b b a 'η=---⎰; (3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n n nnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2.不动点迭代法总是线性收敛的 （）
3.牛顿法有可能不收敛 （）
4.二分法与牛顿法一样都可以推广到多维方程组求解 （）
5.非线性方程（或方程组）的解通常不唯一 （）
三．解答题
1.设Q(x)=x-p(x)f(x)-q(x)f2(x),试确定函数p(x)和q(x),使求解f(x)=0且以Q(x)为迭代函数的迭代法至少三阶收敛。
答案
一．填空题
1.- < <0
2.Xk+1= k=0,1,2…
3.6
二．ห้องสมุดไป่ตู้断题
TF T F T
三．解答题
1.p(x)= ,q(x)=
2.解: 显然,迭代函数为 ,且 ,即 是 的不动点.
又 ,所以 ,
由定理7.4知，迭代是二阶收敛的，且 ．
3. 解 对于 ,因此牛顿迭代法为
根据定理7.4知
对于 ,牛顿法公式为
X=- (Q(x)-3x)=g(x)
|g‘(x)|=| - (Q‘(x)-3x)|= |(Q‘(x)-3x)|< <1
故|g‘(x*)|<1,迭代公式
xk+1=g(xk)= - (Q‘(xk)-3xk)局部收敛
根据定理7.4知
4. 证法一 用数列的办法,因 由 知 ,且 .又由
故 ,即 单减有下界 .根据单调原理知, 有极限.易证起极限为 .
证法二 设 .易知 在 内有惟一实根 .对 应用牛顿迭代法,得
利用例7-9的结论知,当 时, 单减有下界 ,且 .当 时,
此时,从 起, 单减有下界 ,且极限为 .
5.解：由x=Q(x)得 x-3x=Q(x)-3x
数学分析第七章试题
出题人 徐嘉瞳
一．填空题
1.设 2-5),要使迭代过程xk+1= k)局部收敛到x*= ,则 的取值范围________
2.应用牛顿法于方程x3-a=0,导出求立方根 的迭代公式_______
3.用二分法求方程x3-x2-1=0的正根，要进行___次迭代(区间[1,1.5])
二．判断题
2.给定初值 以及迭代公式 ,常数 证明: 该迭代函数是二阶收敛的
3应用牛顿法于方程 和 ,分别导出求 的迭代公式,并求
4.研究求 的牛顿公式xk+1= (xk+ ),x0>0,证明对一切k=1,2,…xk 且序列x1,x2,…是递减的
5.已知方程x=Q(x)在[a,b]内有根x*,且在[a,b]上满足|Q‘(x)-3|<1,利用Q(x)构造一个迭代函数g(x)，使xk+1=g(xk) (k=0,1,2…) 局部收敛于x*