数值分析题库1
数值分析复习题及答案
数值分析复习题一、选择题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. 设2.3149541...x *=，取5位有效数字，那么所得的近似值x=.2.设一阶差商()()()21122114,321f x f x f x x x x --===---，()()()322332615,422f x f x f x x x x --===-- 那么二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 那么2||||X = ，=∞||||X 。
4．求方程 21.250x x --= 的近似根，用迭代公式x =01x =， 那么 1______x =。
5．解初始值问题 00'(,)()y f x y y x y =⎧⎨=⎩近似解的梯形公式是 1______k y +≈。
数值分析试题集
..数值分析试题集（试卷一）一（ 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* 的近似值时，二阶局部收敛的牛顿迭代公式是--------------------------- 。
数值分析期末试题及答案
数值分析期末试题及答案试题一：1. 简答题（共10分）a) 什么是数值分析？它的主要应用领域是什么？b) 请简要解释迭代法和直接法在数值计算中的区别。
2. 填空题（共10分）a) 欧拉方法是一种______型的数值解法。
b) 二分法是一种______法则。
c) 梯形法则是一种______型的数值积分方法。
3. 计算题（共80分）将以下函数进行数值求解：a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。
b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。
c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。
试题二：1. 简答题（共10分）a) 请解释什么是舍入误差，并描述它在数值计算中的影响。
b) 请解释牛顿插值多项式的概念及其应用。
2. 填空题（共10分）a) 数值稳定性通过______号检查。
b) 龙格－库塔法是一种______计算方法。
c) 零点的迭代法在本质上是将方程______转化为______方程。
3. 计算题（共80分）使用牛顿插值多项式进行以下计算：a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4)，求在 x = 0.5 处的插值多项式值。
b) 已知插值节点 (0, 1), (1, 2), (3, 7)，求插值多项式，并计算在 x = 2 处的值。
c) 使用 4 阶龙格－库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。
答案：试题一：1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。
它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。
b) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解，直到满足所需精度为止；而直接法则通过一系列代数运算直接得到解。
数值分析期末试题及答案
数值分析期末试题及答案一、选择题（每题5分，共20分）1. 在数值分析中，下列哪个算法不是用于求解线性方程组的？A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案：B2. 插值法中，拉格朗日插值法属于：A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案：A3. 以下哪个选项不是数值分析中的误差来源？A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案：C4. 在数值积分中，梯形法则的误差项是：A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案：A二、填空题（每题5分，共20分）1. 牛顿插值法中，插值多项式的一般形式为：______。
答案：f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时，迭代公式为：x_{n+1} = x_n -f(x_n) / __________。
答案：f'(x_n)3. 在数值分析中，______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。
答案：误差4. 线性方程组的解法中，______ 法是利用矩阵的LU分解来求解。
答案：克兰特三、解答题（每题10分，共60分）1. 给定函数f(x) = e^(-x)，使用拉格朗日插值法，求x = 0.5时的插值值。
解答：首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1，对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。
拉格朗日插值多项式为：L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得：L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。
数值分析试卷及答案
数值分析试卷及答案数值分析试卷一、选择题（共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析练习题加答案(一)
数值分析期末考试一、 设80~=x ，若要确保其近似数的相对误差限为0.1%，则它的近似数x 至少取几位有效数字？（4分）解：设x 有n 位有效数字。
因为98180648=<<=，所以可得x 的第一位有效数字为8（1分） 又因为21101011000110821--⨯=<⨯⨯≤n ε，令321=⇒-=-n n ，可知x 至少具有3位有效数字（3分）。
二、求矩阵A 的条件数1)(A Cond （4分）。
其中⎥⎦⎤⎢⎣⎡=4231A 解：⎥⎦⎤⎢⎣⎡--=-5.05.1121A （1分） 1A =7（1分） 2711=-A （1分）249)(1=A Cond （1分）三、用列主元Gauss 消元法法求解以下方程组（6分）942822032321321321=++-=++--=+-x x x x x x x x x解：→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----5.245.2405.35.230914220321821191429142821120321 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---8175835005,245.24091425.33.2305.245.2409142(4分) 等价三角方程组为：⎪⎪⎩⎪⎪⎨⎧-=-=+-=++,8175835,5.245.24,942332321x x x x x x （1分）回代得1,3,5123==-=x x x （1分）四、设.0,2,3,1,103)(3210234=-===-+-=x x x x x x x x f 1）求以3210,,,x x x x 为节3次Lagrange 多项式；（6分） 2）求以3210,,,x x x x 为节3次Newton 多项式；（6分）3）给出以上插值多项式的插值余项的表达式（3分）解：由0,2,3,13210=-===x x x x 可得10)(,34)(,1)(,11)(3210-==-=-=x f x f x f x f即得： +------+------=))()(())()(()())()(())()(()()(312101320130201032103x x x x x x x x x x x x x f x x x x x x x x x x x x x f x L=------+------))()(())()(()())()(())()(()(23130321033212023102x x x x x x x x x x x x x f x x x x x x x x x x x x x f+-+--+-⨯-+-+--+-⨯-)03)(23)(13()0)(2)(1()1()01)(21)(31()0)(2)(3(11x x x x x x326610.)20)(30)(10()2)(3)(1()10()02)(32)(12()0)(3)(1(34x x x x x x x x x -+--=+--+--⨯-+---------⨯2）计算差商表如下：i x )(i x f 一阶差商 二阶差商 三阶差商1 -11 3 -1 5 -2 34 -7 4 0-10-225-1则=+-----+-+-=)2)(3)(1()3)(1(4)1(511)(3x x x x x x x N326610x x x -+--3）)2)(3)(1())()()((!4)()(3210)4(3+--=----=x x x x x x x x x x x x f x R ξ五、给定方程组b Ax =，其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100131w w w w A 。
《数值分析》期末复习题(1)
《数值分析》期末复习题一、单项选择题1. 数值x *的近似值x =0.32502×10-1，若x 有5位有效数字，则≤-*x x ( ).(A)21×10-3 (B) 21×10-4 (C) 21×10-5 (D) 21×10-6 2. 设矩阵A ＝10212104135⎡⎤⎢⎥⎢⎥⎣⎦，那么以A 为系数矩阵的线性方程组A X ＝b 的雅可比迭代矩阵为( )(A)00.20.10.200.40.20.60--⎡⎤⎢⎥--⎢⎥--⎣⎦(B)10.20.10.210.40.20.61⎡⎤⎢⎥⎢⎥⎣⎦(C) 00.20.10.200.40.20.60⎡⎤⎢⎥⎢⎥⎣⎦ (D)021204130⎡⎤⎢⎥⎢⎥⎣⎦3. 已知(1)1,(2)4,(3)9f f f ===，用拉格朗日2次插值，则(2.5)f =( )(A) 6.15 (B) 6.25 (C) 6.20 (D) 6.10 4. 抛物形求积公式的代数精度是( )A. 1，B. 2 ，C. 3，D. 45. 改进欧拉格式的 局部截断误差是( ). (),A O h 2. (),B O h 3. (),C O h 4. ().D O h二、填空题1、以722作为π的近似值，它有（ ）位有效数字； 2、经过)1,2( ),2,1( ),1,0(C B A 三个节点的插值多项式为（ ）； 3、用高斯-赛德尔迭代法解方程组⎩⎨⎧-=+-=+,10,232121x bx bx x 其中b 为实数，则方法收敛的充分条件是b 满足条件（ ）；4、取步长为1.0=h ，用欧拉法计算初值问题22',(0)0,y x y y ⎧=+⎨=⎩的解函数)(x y ，它在3.0=x 的近似值为（ ）；5、已知方程0sin 1=--x x 在)1,0(有一个根，使用二分法求误差不大于41021-⨯的近似解至少需要经过（ ）次迭代。
数值分析习题集及答案[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 第三章 函数逼近 6 第四章 数值积分与数值微分 10 第五章 解线性方程组的直接解法 13 第六章 解线性方程组的迭代解法 14 第七章 非线性方程求根 16 第九章 常微分方程初值问题的数值解法 19
第一章 绪论
1.1 要使的相对误差不超过0.1%，应取几位有效
解 对y=f(x)的反函数进行三次插值，插值多项式为
+ + + =, 于是有
。
第三章 函数逼近
3.1证明定义于内积空间H上的函数是一种范数。
证明： 正定性当且仅当时； 齐次性 设为数域K上任一数 三角不等式 ；
于是有 故是H上的一种范数。
3.2求，在空间上的最佳平方逼近多项式，并给出 误差。
解： 第一步：构造内积空间上的一组正交基，其中内积： 第二步：计算的二次最佳平方逼近多项式 从第一步已经知道，利用公式得： 误差为：
数字？
解：
的首位数字。 设有 n位有效数字，由定理知相对误差限 令， 解得，即需取四位有效数字.
1.2 序列满足关系式，若，计算到，误差有多
大？这个算法稳定吗？
解：，于是 ，一般地，因此计算到其误差限为，可见这个计算过程是不稳定的。
1.3 计算球的体积，要使相对误差限为1%，问测 量半径R时允许的相对误差限是多少？
4.1、计算积分，若用复化梯公式，问区间应分多 少等份才能使截断误差不超过？若改用复化辛普 森公式，要达到同样的精度，区间应分多少等 份？
解：由于，，,故对复化梯公式，要求 ，
即，.取，即将区间分为等份时，用复化梯公式计算，截断误差不超过. 用复化辛普森公式，要求 ，
即，.取，即将区间等分为8等份时，复化辛普森公式可达精度.
解：
建立改进欧拉的迭代公式 n=0,1,2,3,4
迭代计算如下：
.1 0.00550 .2 0.02193 .3 0.05015 .4 0.09094 .5 0.14500
0.00516 0.02127 0.04918 0.08968 0.14347
X
1
3/2
0
2
f(x)
3
13/4
解：由函数值表建立差商表如下：
x
f(x)
一阶差商
3 二阶差商
5/3 三阶差商
1
3
3/2 0 2
13/4 3 5/3
1/2 1/6 -2/3
1/3 -5/3
-2
进而得牛顿多项式为
2.2、已知试选用合适的插值节点利用Lagrange
二次插值多项式计算的近似值，使之精度尽可能
超过0.0001的近似根。
解：
(1)迭代格式为在[1，2]上满足,因此迭代发散。 （2），又当时，，因此迭代格式收敛，取，则有
。
7.4 应用牛顿法于方程，导出求的迭代公式，由
此求。
解：
建立迭代格式： ，，因此该迭代格式在附近局部收敛到，且是2阶收敛的。取，迭代计算如下 ，
7.5确定常数p,q,r使迭代公式产生的序列收敛
3.3 确定参数，使得积分 取得最小值。
解：
定义内积，其中
，切比雪夫多项式关于此内积正交，其前三项为：
并且有
利用正交基最佳平方逼近多项式公式得二次最佳平方逼近多项
式有：
由最佳平方逼近的定义知，当时，积分取最小值，此时参数
为：
3.4 已知函数值表
x -2
-1
0
1
2
y0
1
2
1
0
试用二次多项式拟合这组数据。
2.4、求f(x)=在[a,b]上的分段线形插值函数， 并估计误差。
解：设采用节点, 定义 在[]上的线性插值函数 分段线性插值函数
误差估计： 而
2.5、已知单调连续函数y=f(x)在如下采样点处 的函数值：
1.0 1.4 1.8 2.0 -2.0 -0.8 0.4 1.2 求方程在[1，2]内根的近似值，使误差尽可能 小。
解：（1），，根据定理7.3，当 =<1，亦即<C<0时收敛。
(2）由定理7.4知，当，即C=时迭代至少是二阶收敛的，收敛最快。 （3）分别取C=，，并取，迭代计算结果如下表所示。
k
0
1.2
k
0
1.2
1
1.48
1
6
1.413369586
2
12
1.414209303
3
13
1.414215327
4
1.397989899 1.414120505 1.414213559 1.414213562
4.4、用三点高斯公式计算。
解：积分区间[1，3]，要使用高斯公式，需先变换到[-1，1]上，为此，作 变换：
y=(a+b)+(b-a)t = t + 2 则当y[1，3]时， t[-1，1], 且 dy = dt, =. 所以，三点高斯公式 =0.555 555 6（+）+
0.888 888 9 = 1.098 039 3
高。
解：依据误差估计式，选为插值节点，拉格朗日插值基函数为
二次插值多项式为
L2(x)= f(x0) ＋f(x1) ＋f(x2) ＝＋2 ＋3 ＝
于是 f(-0.5) L2(－0.5)=1l0(－0.5)+ 2l1(－0.5)+ 3l2(－0.5) =
2.3、已知函数的如下函数值表，利用插值法计算
的近似值。
4.2、求高斯型求积公式 的系数，及节点。
解：上面的求积公式为高斯型求积公式，故对 故有
精确成立的。
写成矩阵形式
上两式右端写成等式，并化简有：
令
得
从而有
由前两个方程，得 所求高斯型求积公式为
4.3、定下列求积分公式的待定参数，使得求积分 公式的代数精度尽可能的高，并指明所确定的求 积分公式所具有的代数精度。 （1）（其中A、B、C为待定参数） （2）（其中、为待定参数）
解： ，即测量半径R时允许的相对误差限是。
1.4 试给出一种计算积分近似值的稳定递推算法
解： 对分部积分得到，进而解出， 由此公式可以看出，当有误差时，有误差，依次类推得到的误差为，随着递推过程进行，初始 误差对计算结果的影响越来越小，因而是一种稳定的运算公式。
第二章 函数插值
2.1、利用如下函数值表构造差商表，并写出牛顿 插值多项式。
解： （1） 求积公式中含有3个待定参数A、B、C，故令求积公式对精 确成立，即
解得： 又由 故具有三次代数精度。
（2） 求积公式中含有两个待定参数，，当时，有 故令求积公式对精确成立，即
解得 将代入求积公式，显然 故求积公式，当求积节点取或时，具有二次代数精度。 令 得n2.54。故取n=3，即需取2n+1=7个节点处的函数值。
此时都达到<。事实上，=1.414213562
7.2曲线 与 在点（1.6，1）附近相切，试用牛 顿迭代法求切点横坐标的近似值，使.
解：两曲线的导数分别为 和 ，两曲线相切，导数相等，故有
令，则,故区间[1，2]是的有根区间。又当时， ，因此在区间[1，2]上有唯一实根 ，对f(x) 应用牛顿迭代法，得计算公式
到，并使收敛的阶尽可能高。
解：
迭代函数，要使迭代序列收敛阶尽可高，应使由 即 （1），由 （2），同样由 （3），由（1），（2），（3）得到，而且由于知该迭代是三阶收敛的。
第九章 常微分方程初值问题的数值解法
9.1取步长h=0.2，用欧拉法求解初值问题
解：
建立迭代公式： 迭代计算得：
9.2取步长h=0.1, 用改进欧拉法求解初值问题
第五章 解线性方程组的直接解法
5.1、用分解法解方程组
解：设 由矩阵乘法得： 由 得 由 得
5.2设，，已知方程组的精确解为。 （1）计算条件数； （2）若近似解，计算剩余；
解：（1） （2）
第六章 解线性方程组的迭代解法
6.1对方程组，用雅可比迭代法求解是否收敛？若 收敛，取初始向量，迭代计算至。
解：雅可比迭代矩阵
故雅可比迭代法收敛。 雅可比迭代格式为 取初始值，计算如下表示：
k
0
0
0
0
1
1
3
5
2
5
-3
-3
3
1
1
1
4
1
1
1
由 ，故取
6.2对于方程组
（1）求证用jacobi迭代法和G-S迭代求解此方
程组均收敛.
（2）取,用jacobi迭代法和G-S迭代法求解要
求.
解：①，故雅可比迭代收敛。 ，所以G-S迭代收敛。 ②雅可比迭代格式为： 取初值迭代计算如下： 1 0.1 0.05 0.3333333 2 0.176667 0.1033333 0.4000000 3 0.200667 0.1253333 0.4611111 4 0.217289 0.1362450 0.4837777 … …. …. … 14 0.231081 0.1470411 0.5083260 15 0.231087 0.1470550 0.5083930 G-S迭代格式为： 取初值迭代计算如下：
由于，故取迭代计算一定收敛，计算结果如下表所示。
k
k
0
2.0 3
1 2.293 055 4 556
2 1.817 783 5 592
继续计算仍得 =1.7，故
1.706 815 287
1.700 025 611
1.7
7.3用简单迭代法求在[1，2]上的根，若迭代函
数分别取为和,问他们是否收敛？若收敛求误差不
x
0.4
0.5