应用数值分析第四版第一章课后作业答案
应用数值分析第四版第一章课后作业答案
第一章1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈-=-=a ee x a e r π （2）0011.00143.0143.07/1≈=≈-=-=a ee x a e r （3）0127.000004.00031.01000/≈=≈-=-=aee x a e r π （4）001.00143.03.147/100≈=≈-=-=aee x a e r2. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.5019373、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
最新应用数值分析第四版第一章课后作业答案
第一章1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈-=-=a ee x a e r π （2）0011.00143.0143.07/1≈=≈-=-=a ee x a e r （3）0127.000004.00031.01000/≈=≈-=-=aee x a e r π （4）001.00143.03.147/100≈=≈-=-=aee x a e r2. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.5019373、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
数值分析课第一次作业答案answer1
计算机习题： 1. 作多项式 p，以 −1，0，1 为零点，首项系数为 2，并计算 p(3)。 4
答案：p = poly ([−1, 0, 1])，s = polyval(p, 3)。 2. 已知函数在下列各点的值为 xi 0.2 0.4 0.6 0.8 1.0
2
a 6 6e+154 0 1 1
b 10 10e+154 1 -1e+5 -4
c -4 -4e+154 1 1 3.999999
-1e+155 -7e+155 1e+155 答案：第二种方法更准确，因为第一种方法是一个累加的过程。 matlab 的 x = a : h : b 和 x = a + (0 : n) ∗ h 是第二种方法实现的。 代码： format long e a = 0; b = 8; n = 9; h = (b-a)/n; x(1) = a; y(1) = a; for j = 1:n, x(j+1) = x(j) + h; y(j+1) = y(1) + j*h; end [x',y',(a:h:b)',a+(0:n)’*h] 第二章 插值法 1. 当 x = 1, −1, 2 时，f (x) = 0, −3, 4，求 f (x) 的二次插值多项式。 （计算两遍，分别用拉格朗日插值和牛顿插值）
5
f (xi ) 0.98 0.92 0.81 0.64 0.38 求 4 次牛顿插值多项式 P4 (x) 并画图。 答案： 代码： x=0.2:0.2:1.0; y=[0.98,0.92,0.81,0.64,0.38]; n = length(y); if length(x)~=n, error('x and y are not compatible'); end D = zeros(n,n); D(:,1)=y(:); for j=2:n for i=j:n D(i,j) = (D(i,j-1)-D(i-1,j-1))/(x(i)-x(i-j+1)); end end p=D(1,1)*[zeros(1,n-1),1]; for k=2:n p=p+D(k,k)*[zeros(1,n-k),poly(x(1:k-1))]; end x=0.2:0.01:1.0; z=polyval(p,x); plot(x,z) 比较：p = polyf it(x, y, 4)。
数值分析_第四版_课后习题答案_李庆扬
第一章1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x xx x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而n x 的误差为nn x x nn xn x x n x x x **1***%2%2)()()()(*⋅=='=-=εε，相对误差为%2)()x ()x (*n *n*n x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++；[解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
应用数值分析(第四版)课后习题答案第2章
第二章习题解答1.（1） R n×n中的子集“上三角阵”和“正交矩阵”对矩阵乘法是封闭的。
（2）R n×n中的子集“正交矩阵”，“非奇异的对称阵”和“单位上（下）三角阵”对矩阵求逆是封闭的。
设A 是ｎ×ｎ的正交矩阵。
证明A -1也是ｎ×ｎ的正交矩阵。
证明：(1),n nA B A B R⨯∈证明：为上三角阵，为上三角阵，10(),0(),0(),,()(()),()()ij ij nij ik kj ij k n n T T T T T T T T T T a i j b i j C AB c a b c i j A B A B R AA A A E BB B B EAB AB ABB A E AB AB B A AB E AB =⨯∴=>=>==∴=>∴∈========∴∑则上三角阵对矩阵乘法封闭。
以下证明：为正交矩阵，为正交矩阵，为正交矩阵,故正交矩阵对矩阵乘法封闭。
（2）A 是ｎ×ｎ的正交矩阵∴A A -1 =A -1A=E 故（A -1）-1=A∴A -1（A -1）-1=（A -1）-1A -1 =E 故A -1也是ｎ×ｎ的正交矩阵。
设A 是非奇异的对称阵，证A -1也是非奇异的对称阵。
A 非奇异 ∴A 可逆且A -1非奇异 又A T =A ∴（A -1）T =（A T ）-1=A-1故A -1也是非奇异的对称阵设A 是单位上（下）三角阵。
证A -1也是单位上（下）三角阵。
证明：A 是单位上三角阵，故|A|=1，∴A 可逆，即A -1存在，记为（b ij ）n×n由A A -1=E ，则∑==nj ik jkij ba 1δ （其中0=ij a j ＞i 时，1=ii a ）故b nn =1, b ni =0 (n≠j)类似可得，b ii =1 (j=1…n) b jk =0 (k ＞j)即A -1是单位上三角阵综上所述可得。
数值分析第四版答案资料讲解
数值分析第四版答案第一章绪论1.设x 0,x的相对误差为，求In x的误差。
解：近似值x*的相对误差为* e* x* x =ex* x*而In x 的误差为el nx* Inx* In x e* x*进而有(In x*)2.设x的相对误差为2%，求 E x n的相对误差。
解：设f(x) x n，则函数的条件数为C p丨空^丨f(x)H n 1又 f '(x) nx n 1, C p | x nx | n1n—11又「((x*) n) C p r(x*)且e (x*)为2r((x*)n) 0.02 n3•下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指出它们是几位有效数字：x；1.1021,x2 0.031 , x3 385.6,沧56.430 ,x57 1.0.解：x1 1.1021是五位有效数字；x2 0.031是二位有效数字；x；385.6是四位有效数字；x4 56.430是五位有效数字；X；7 1.0.是二位有效数字。
4•利用公式(2.3)求下列各近似值的误差限：(1) x；x；x；,(2) x；x；x；,(3) x；/x；. 其中X1,X2,X3,x4均为第3题所给的数。
解:*1 (X 1)2 10(1) (X 1X 2 X 4)（X ；）（x 2） （x 4）11021.05 10(2) (x ；x ；x ；)(3) (X 2/X 4) * I **X 2I(X 4) X 4* 2 X40.031 1 3 13-10 56.430 — 102 2 10 5 56.430 56.4305计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 4 °解：球体体积为V - R 33*1（X 2） 2 10 * 1 （X 3） 2 10 * 1 （X 4） 2 10 * 1 （X 5— 123 131101103X 1X 2 (X 3) 1.1021 0.031 0.215X 2X 31 2101X 1X 3 (X 2)10.031 385.6 - 101.1021 385.6 1 103*(X 2)则何种函数的条件数为C R(4 R2 4 R3 3r(V*) Cp|「(R*)3 r (R*)又；r (V*)11故度量半径R 时允许的相对误差限为r (R*) - 1 0.33 36 •设 Y o 28，按递推公式 Y, Y n-1,783 (n=1,2,…) 100计算到丫100。
应用数值分析西安电子科技大学课后答案
应用数值分析西安电子科技大学课后答案1. 大数据中的小数据可能缺失、冗余、存在垃圾数据，但不影响大数据的可信数据，是大数据的（）的表现形式。
[单选题] *A. 价值涌现B.隐私涌现C. 质量涌现(正确答案)D. 安全涌现2. 数据科学基本原则中，基于数据的智能的主要特点是（）。
[单选题] *A. 数据简单，但算法简单B.数据复杂，但算法简单(正确答案)C. 数据简单，但算法复杂D. 数据复杂，但算法复杂3. （）是数据库管理系统运行的基本工作单位。
[单选题] *A. 事务(正确答案)B.数据仓库C. 数据单元D. 数据分析4. 目前，多数NoSQL 数据库是针对特定应用场景研发出来的，其设计遵循（）原则，更强调读写效率、数据容量以及系统可扩展性。
[单选题] *B. READC. BASE(正确答案)D. BASIC5. 数据可视化的本质是（）。
[单选题] *A. 将数据转换为知识(正确答案)B.将知识转换为数据C. 将数据转换为信息D.将信息转换为智慧6.下列不属于大数据在社会活动中的典型应用的是（）。
[单选题] *A. 美团实现了快速精准的送餐服务B. 共享单车、滴滴打车方便了人们的日常出行C. 快递实现了订单的实时跟踪D. 供电公司提供电费账单查询(正确答案)7.在空间维度上刻画数据连续性是数据的（）。
[单选题] *A. 可关联性(正确答案)B.可溯源性C. 可理解性D.可复制性8.将观测值分为相同数目的两部分，当统计结果为非对称分布时经常使用的是（）。
[单选题] *B.标准差C. 中位数(正确答案)D.均值9. （）的本质是将低层次数据转换为高层次数据的过程。
[单选题] *A. 数据处理B.数据计算C. 数据加工(正确答案)D.整齐数据10. 在抽样方法中，当合适的样本容量很难确定时，可以使用的抽样方法是（）。
[单选题] *A. 有放回的简单随机抽样B. 无放回的简单随机抽样C. 分层抽样D.渐进抽样(正确答案)11.下列关于基本元数据描述正确的是（）。
数值分析课后答案chap1
∴ ε r (( x*) ) ≈ 0.02 n
3．下列各数都是经过四舍五入得到的近似 数，即误差限不超过最后一位的半个单位， 试 指 出 它 们 是 几 位 有 效 数 字 ：
∫
N +1
N
9．正方形的边长大约为了 100cm，应怎样 测量才能使其面积误差不超过 1cm 2 ？
解：正方形的面积函数为 A( x ) = x 2
∴ y1 = 10 y0 − 1 ∴ ε ( y1*) = 10ε ( y0 *)
又∵ y2 = 10 y1 − 1
∴ ε ( A*) = 2 A *iε ( x*) .
* * * (1)ε ( x1 + x2 + x4 ) * * * = ε ( x1 ) + ε ( x2 ) + ε ( x4 )
又∵ f '( x) = nx n−1 , ∴ C p =| 又∵ ε r (( x*) n) ≈ C p ⋅ ε r ( x*) 且 er ( x*) 为 2
1 1 1 = ×10−4 + ×10−3 + ×10−3 2 2 2 −3 = 1.05 × 10
x = 56.430 , x = 7 ×1.0.
* 解： x1 = 1.1021 是五位有效数字； * x2 = 0.031 是二位有效数字； * x3 = 385.6 是四位有效数字； * x4 = 56.430 是五位有效数字；
≈
1 1 0.031× × 10−3 + 56.430 × ×10 −3 2 2 = 56.430 × 56.430 −5 = 10
数值分析课后习题答案
第一章习题解答1. 在下列各对数中，X 是精确值ａ的近似值（1） ａ=π，x=3.1 （2） ａ=1/7，x=0.143 （3） ａ=π/1000，x=0.0031 （4） ａ=100/7，x=14.3 试估计x 的绝对误差和相对误差。
解：（1） e=∣3.1-π∣≈0.0416， δr = e/∣x ∣≈0.0143 （2） e=∣0.143-1/7∣≈0.0143 δr = e/∣x ∣≈0.1 （3） e=∣0.0031-π/1000∣≈0.0279 δr = e/∣x ∣≈0.9 (4) e=∣14.3-100/7∣≈0.0143 δr = e/∣x ∣≈0.0012. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［-x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.497073. 设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
数值分析课后习题及答案
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
(完整版)数值分析第四版习题和答案解析
第四版数值分析习题第一章 绪 论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 -=…)计算到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 cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xxx ----==证明()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 j j x l x xk n =≡=∑ii)()()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 -≤≤上给出()x f x e =的等距节点函数表,若用二次插值求x e 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n 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 →∞时,()nx ϕ在[],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 =<<<=,则[][][]()()()()()()()()()baS 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. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f 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 使[]22sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义()(,)()();()(,)()()()();bbaaa 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 ()nn x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.2y a bx =+.用最小二乘拟合求.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)()h h f 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()()()()()2ba f f x dxb 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.1xedx-,要求误差不超过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. 设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 cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n 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 j j x l x xk n =≡=∑ii)()()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 b a x b f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()x f x e =的等距节点函数表,若用二次插值求x e 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n 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()()()bbbba a a a f 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 =<<<=,则[][][]()()()()()()()()()baS 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. 试证{}*()nT x 是在[]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 使[]22sin 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()sin 2f 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)()hh f 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()()()()()2ba f f x dxb 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 nn n ππππ=-+-试依据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相比较。
数值分析第一章绪论习题答案
第一章绪论e In X* =In X * -Inx ：丄e*X*进而有；(In X *):2. 设X 的相对误差为2% ,求X n 的相对误差。
解：设f(χZ ,则函数的条件数为Cp=l fX+n _1X nχ I Xn n又；r ((X*) n) C P 7(X *)且 e r (χ*)为 2.7((χ*)n) 0.02 n3. 下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指* * * * *出它们是几位有效数字： X 1 =1.1021, χ2 =0.031, χ3 =385.6, χ4 = 56.430,x 5 = 7".0.. *解：X I -1.1021是五位有效数字；X 2 = 0.031是二位有效数字；X 3 =385.6是四位有效数字；X 4 =56.430是五位有效数字；X 5 =7 1.0.是二位有效数字。
4. 利用公式(2.3)求下列各近似值的误差限： (1) X 1 X 2 X 4,(2) X 1 X 2X 3 ,(3) X 2 /X 4 .其中χl ,x 2,x 3,X 4均为第3题所给的数。
1设X 0, x 的相对误差为 解：近似值X*的相对误差为 、：，求InX 的误差。
e* X* -X而InX 的误差为 又 f '(χ) =nx n 」 C P解：* 1 4；(x 1) 102* 1 3 ;(x 2) 10 2* 1 1；(x 3) 10* 1 3；(x 4) 102* 1 1；(x 5) 102(1) ；(x ； x ； x *)* * *=；(%) ；(x 2) *x 4)1 A 12 1 j310 10 102 2 2 -1.05 10J 3* * *(2) S(X I X 2X 3)* * * * * * ** * =X1X 2 £(X 3)+ X 2X 3 ^(X J + X 1X 3 E (X 2):0.215 ⑶；(x 2/x ；)* Il * * I * X 2 E(X 4) + X 4 &(X 2)全 Γ"2X 41-3 1 30.031 10 56.430 10= ______________________ 256.430X56.430-10 54 3解：球体体积为V R3则何种函数的条件数为1.1021 0.031 11θ' 2 + 0.031X385.6 x 1><10* 2 +∣ 1.1021 X 385.6卜-×1^35计算球体积要使相对误差限为 1 ,问度量半径R 时允许的相对误差限是多少?C P 愕'： C P “(R*) 9(R*)又γ(V*) -11故度量半径R 时允许的相对误差限为 ；r (R*) 1 ： 0.3331 ____6.设 Y 0 =28,按递推公式 Yn =Ynd- ------- ： 783 (n=1,2,…)100计算到Y oo 。
应用数值分析习题答案
应用数值分析习题答案应用数值分析习题答案数值分析是一门研究利用数值方法解决数学问题的学科。
在实际应用中，我们常常会遇到一些复杂的数学问题，而利用传统的解析方法往往效果不佳。
因此，数值分析的方法成为了解决这些问题的重要手段。
在学习数值分析的过程中，习题是不可或缺的一部分，通过解答习题，我们可以更好地理解和掌握数值分析的方法和技巧。
习题一：求解非线性方程非线性方程的求解是数值分析中的重要问题之一。
考虑如下的非线性方程：f(x) = x^3 - 2x - 5 = 0我们可以利用数值方法，例如二分法、牛顿法等来求解这个方程。
二分法是一种简单而有效的方法，它的基本思想是通过不断缩小方程的解所在的区间来逼近方程的解。
牛顿法则是通过对方程进行线性化，利用切线的斜率来逼近方程的解。
习题二：插值与拟合插值与拟合是数值分析中的重要内容，它们在实际问题中的应用非常广泛。
插值是指通过已知的数据点，构造出一个函数，使得这个函数通过已知数据点，并且在这些数据点之间的取值也符合一定的规律。
拟合则是指通过已知的数据点，找到一个函数，使得这个函数与已知数据点之间的误差最小。
习题三：求解常微分方程常微分方程是数值分析中的重要内容之一。
常微分方程描述了许多实际问题的变化规律。
通过数值方法，我们可以求解常微分方程的近似解。
常用的数值方法有欧拉法、改进的欧拉法、龙格-库塔法等。
这些方法的基本思想是将微分方程转化为差分方程，通过迭代的方式逼近微分方程的解。
习题四：线性方程组的求解线性方程组是数值分析中的重要问题之一。
在实际问题中，我们经常会遇到大规模的线性方程组，通过传统的解析方法往往效率低下。
因此，利用数值方法来求解线性方程组成为了一种常用的手段。
常用的数值方法有高斯消元法、LU分解法、迭代法等。
习题五：数值积分数值积分是数值分析中的重要内容之一。
在实际问题中，我们经常需要计算曲线的面积、求解定积分等。
通过数值方法，我们可以将曲线的面积或者定积分的计算转化为离散数据点的求和。
(完整版)数值分析课后习题答案
第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式(1.2.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),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
应用回归分析第四版课后习题答案-全
实用回归分析第四版第一章回归分析概述1.3回归模型中随机误差项ε的意义是什么？答：ε为随机误差项，正是由于随机误差项的引入，才将变量间的关系描述为一个随机方程，使得我们可以借助随机数学方法研究y与x1,x2…..xp的关系，由于客观经济现象是错综复杂的，一种经济现象很难用有限个因素来准确说明，随机误差项可以概括表示由于人们的认识以及其他客观原因的局限而没有考虑的种种偶然因素。
1.4 线性回归模型的基本假设是什么？答：线性回归模型的基本假设有：1.解释变量x1.x2….xp是非随机的，观测值xi1.xi2…..xip是常数。
2.等方差及不相关的假定条件为{E(εi)=0 i=1,2…. Cov(εi,εj)=｛σ^23.正态分布的假定条件为相互独立。
4.样本容量的个数要多于解释变量的个数，即n>p.第二章一元线性回归分析思考与练习参考答案2.1一元线性回归有哪些基本假定?答：假设1、解释变量X是确定性变量，Y是随机变量；假设2、随机误差项ε具有零均值、同方差和不序列相关性：E(εi)=0 i=1,2, …,nVar (εi)=σ2i=1,2, …,nCov(εi,εj)=0 i≠j i,j= 1,2, …,n假设3、随机误差项ε与解释变量X之间不相关：Cov(X i, εi)=0 i=1,2, …,n假设4、ε服从零均值、同方差、零协方差的正态分布εi~N(0, σ2) i=1,2, …,n2.3 证明（2.27式），∑e i =0 ，∑e i X i=0 。
证明：∑∑+-=-=niiiniXYYYQ12121))ˆˆ(()ˆ(ββ其中：即： ∑e i =0 ，∑e i X i =02.5 证明0ˆβ是β0的无偏估计。
证明：)1[)ˆ()ˆ(1110∑∑==--=-=ni i xxi n i i Y L X X X Y n E X Y E E ββ )] )(1([])1([1011i i xx i n i i xx i ni X L X X X n E Y L X X X n E εββ++--=--=∑∑==1010)()1(])1([βεβεβ=--+=--+=∑∑==i xx i ni i xx i ni E L X X X nL X X X n E 2.6 证明 证明：)] ()1([])1([)ˆ(102110i i xxi ni ixx i ni X Var L X X X n Y L X X X n Var Var εβββ++--=--=∑∑== 222212]1[])(2)1[(σσxx xx i xx i ni L X n L X X X nL X X X n +=-+--=∑=2.7 证明平方和分解公式：SST=SSE+SSR证明：2.8 验证三种检验的关系，即验证： （1）21)2(r r n t --=；（2）2221ˆˆ)2/(1/t L n SSE SSR F xx ==-=σβ 01ˆˆˆˆi i i i iY X e Y Y ββ=+=-())1()1()ˆ(222122xx ni iL X n X XX nVar +=-+=∑=σσβ()()∑∑==-+-=-=n i ii i n i i Y Y Y Y Y Y SST 1212]ˆ()ˆ[()()()∑∑∑===-+--+-=ni ii ni i i i ni iY Y Y Y Y Y Y Y 12112)ˆˆ)(ˆ2ˆ()()SSESSR )Y ˆY Y Y ˆn1i 2ii n1i 2i +=-+-=∑∑==0100ˆˆQQββ∂∂==∂∂证明：（1）ˆt======（2）2222201111 1111ˆˆˆˆˆˆ()()(())(()) n n n ni i i i xxi i i iSSR y y x y y x x y x x Lβββββ=====-=+-=+--=-=∑∑∑∑2212ˆ/1ˆ/(2)xxLSSRF tSSE nβσ∴===-2.9 验证（2.63）式：2211σ)L)xx(n()e(Varxxii---=证明：0112222222ˆˆˆvar()var()var()var()2cov(,)ˆˆˆvar()var()2cov(,())()()11[]2[]()1[1]i i i i i i ii i i ii ixx xxixxe y y y y y yy x y y x xx x x xn L n Lx xn Lβββσσσσ=-=+-=++-+---=++-+-=--其中：222221111))(1()(1))(,()()1,())(ˆ,(),())(ˆ,(σσσββxxixxiniixxiiiniiiiiiiiLxxnLxxnyLxxyCovxxynyCovxxyCovyyCovxxyyCov-+=-+=--+=-+=-+∑∑==2.10 用第9题证明是σ2的无偏估计量证明：2221122112211ˆˆ()()()22()111var()[1]221(2)2n ni ii in niii i xxE E y y E en nx xen n n Lnnσσσσ=====-=---==----=-=-∑∑∑∑第三章2ˆ22-=∑neiσ1.一个回归方程的复相关系数R=0.99，样本决定系数R 2=0.9801，我们能判断这个回归方程就很理想吗？ 答：不能断定这个回归方程理想。
数学分析第四版上册答案
数学分析第四版上册答案第一章环境建立1.1 算术基础在数学分析中，我们需要对数学中的基本运算进行复习和巩固。
这包括四则运算、乘方和开方等。
在本节中，我们将回顾这些基本算术技巧，并解答一些相关问题。
1.1.1 四则运算四则运算是我们进行数学计算的基本方法。
它包括加法、减法、乘法和除法。
在本节中，我们将通过一些例题来练习四则运算，并解答相应的问题。
例题1.1.1计算下列算式的结果：a) 2 + 3 * 4b) (5 - 2) * 7c) 10 / 5 + 3d) 8 - 6 / 2解答：a) 2 + 3 * 4 = 2 + 12 = 14b) (5 - 2) * 7 = 3 * 7 = 21c) 10 / 5 + 3 = 2 + 3 = 5d) 8 - 6 / 2 = 8 - 3 = 5计算下列算式的结果：a) 5 + 6 * 2 - 3b) 8 / 2 * (4 + 3)c) 7 - 4 / 2 + 5 * 3解答：a) 5 + 6 * 2 - 3 = 5 + 12 - 3 = 14 - 3 = 11b) 8 / 2 * (4 + 3) = 4 * 7 = 28c) 7 - 4 / 2 + 5 * 3 = 7 - 2 + 15 = 201.1.2 乘方与开方乘方和开方是我们在数学中经常用到的运算符。
乘方表示多次相乘，开方则相反，表示求一个数的平方根。
在本节中，我们将练习一些乘方和开方的计算，并解答相关问题。
例题1.1.3计算下列算式的结果： a) 2^3b) 4^0.5c) (23)2d) (32)3解答：a) 2^3 = 2 * 2 * 2 = 8 b) 4^0.5 = √4 = 2 c) (23)2 = 8^2 = 64 d) (32)3 = 9^3 = 729计算下列算式的结果：a) √9b) √(4^2)c) √(3^2 + 4^2)d) (√2 + 1)^2解答：a) √9 = 3 b) √(4^2) = √16 = 4 c) √(3^2 + 4^2) = √(9 + 16) = √25 = 5 d) (√2 + 1)^2 = (1.414 + 1)^2 = 2.414^2 = 5.8291.2 方程与不等式在数学分析中，方程和不等式是我们经常遇到和解决的问题。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一章1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈-=-=a ee x a e r π （2）0011.00143.0143.07/1≈=≈-=-=a ee x a e r （3）0127.000004.00031.01000/≈=≈-=-=aee x a e r π （4）001.00143.03.147/100≈=≈-=-=aee x a e r2. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.5019373、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
解：设=()u f x ，()()()()()()||||||||||()||()||||()||()||||r r rx e u df x e x df x e x e u u dx u dx u x df x x df x x e x x dx u dx u δ=≈==≤()||10.2(())||()||ln ln ln r r r r df x x x x f x x x dx u x x x xδδδδ==⋅⋅==4、长方体的长宽高分别为50cm ，20cm 和10cm ，试求测量误差满足什么条件时其表面积的误差不超过1cm 2。
解：设2()S xy yz zx =++{}[]{}(,,)(,,)(,,)()||()||()||()(,,)(,,)(,,)||||||max (),(),()2()2()2()max (),(),()1S x y z S x y z S x y z e S e x e y e z x y zS x y z S x y z S x y z e x e y e z x y z y z z x x y e x e y e z ∂∂∂≤++∂∂∂⎛⎫∂∂∂≤++ ⎪∂∂∂⎝⎭=+++++<{}[]11max (),(),()2()2()2()4()110.0031254(502010)320e x e y e z y z z x x y x y z <=+++++++===++所以，测量误差小于0.00625时其表面积的误差不超过1cm 2。
5、设x 和y 的相对误差为0.001，则xy 的相对误差约为多少？解：由公式：i r ini n in ni i ir x x fx x f x x x f x x f u δδδ∂∂=∂∂=∑∑==1111),,(),,()(则有：002.0001.0001.0)()()()(=+=+=≤y x xy xy e r r r r δδδ xy 的相对误差约为0.002.6. 改变下列表达式，使计算结果更准确。
（1）1,||1x x x +-≥ （2）11,||1121xx x x--≤++ （3）(1cos ),0,||1x x x x -≠≤ （4）11,||1x x x x x+--≥解：（1）111x x x x+-=++ （2）2112121(12)(1)x x x x x x --=++++ （3）2(1cos )sin (1cos )x x x x x -=+（4）22112(11)x x x xx x x +--=++-7、计算6(21)-的近似值，取2 1.414≈。
利用以下四种计算格式，试问哪一种算法误差最小。
（1）61(21)+ （2）3(322)- （3）31(322)+ （4）99702-解：计算各项的条件数'()(())||()xf x cond f x f x = 11 1.41461(),(())| 3.5145(1)x f x cond f x x ===+ 322 1.414()(32),(())|49.3256x f x x cond f x ==-= 331.41431(),(())| 1.4558(32)xf x c o n d f x x ===+ 441.414()9970,(())|4949x f x x c o n d f x ==-= 由计算知，第一种算法误差最小。
n n 118 ∞=∑、考虑无穷级数，它是微积分中的发散级数。
在计算机上计算该级数的部分和，会得到怎样的结果？为什么？解：在计算机上计算该级数的是一个收敛的级数。
因为随着n 的增大，会出现大数吃小数的现象。
9、 通过分析浮点数集合F=（10，3，-2，2）在数轴上的分布讨论一般浮点数集的分布情况。
解：浮点数集合F=（10，3，-2，2）在数轴上离原点越近，分布越稠密；离原点越远，分布越稀疏。
一般浮点数集的分布也符合此规律。
10、试导出计算积分1(1,2,3,4)14n n x I dx n x ==+⎰的递推计算公式111()4n n I I n -=-，并分析此递推公式的数值稳定性。
解：111111110000141()14414414n n n n n n n x x x x x I dx dx x dx dx x x x ----+-===-+++⎰⎰⎰⎰ 111()4n n I I n -∴=- 分析计算的误差，设初值0I 的误差为000I I e -=，递推过程的舍入误差不计，并记n n n I I e -=，则有 。
显然，随着计算的递推，误差越来越小，因此递推公式是稳定的。
11、为减少乘除法运算次数，应将下面算式怎样改写？0114)1(...)(41e I I I I en nn n n n n -==--=-=--3217151318)()(-+-+-+=x x x y 解：令11-=x u ，则18357+++=u u u y ))(( 12、试推导求函数值)(x f y =的绝对误差和相对误差。
解：**(())()()'()()'()()e f x f x f x f x x f x e x ξ=-=-≈***()()'()()'()()'()(())()()()()()r r f x f x f x x f x x x f x e f x x x e x f x f x f x x f x ξ---==≈=13、试推导求函数值),(y x f 的条件数)),((y x f Cond 。
解：(,)(,)(,)(,)()()f x y f x y f x y f x y x x y y x y∂∂=+-+-+∂∂ ''(,)(,)(,)(,)()()(,)(,)(,)y x yf x y xf x y f x y f x y x x y y f x y f x y x f x y y---≈+''''(,)(,)(,)(,)()()(,)(,)(,)(,)(,)()()max ,(,)(,)y x yxyf x y xf x y f x y f x y x x y y f x y f x y x f x y y yf x y xf x y x x y y f x y f x y x y ---≤+⎛⎫⎧⎫--⎪⎪≤ +⎪⎨⎬ ⎪⎪⎪⎩⎭⎝⎭'(,)'(,)((,))(,)(,)x y xf x y yf x y Cond f x y f x y f x y ∴=+14、序列{}n x 满足递推公式⎩⎨⎧⋅⋅⋅=-=≈=+),,(.1057321310n n x x x n n 求计算到20x 的的误差，并讨论计算过程的稳定性。
解：000x x e -=000111555e x x x x e =-=-=021*********e x x x x e =---=-=)()(…9020201084545⨯≈=.e e误差逐渐增大，计算不稳定15、写出下面Matlab 程序所描述的数学表达式。
（1）for j=1: nfor i=1: my ( i )= A (i, j)*x(j)+y(i) end end (2)for j=1:ny=x(j)*A(:,j)+y end解：（1）m n n m R y R x R A y Ax y ∈∈∈+=⨯,,,（2）y Ax y += A 为n 列的矩阵，x 为n 行1列的列向量，y 为与A 有相同行数的列向量16、写出下面Matlab 程序所描述的数学表达式。
(1) for i=1: mfor j=1: nA(i,j)= A (i, j)+x(i)*y(j) end end(2) for j=1:mA(:,j)= A (:, j)+ y(j)* x(:) end解：(1) n m nm R y R x RA y x A A ⨯⨯⨯∈∈∈+=11,,,*(2) x y A A *+= A 为m 列的矩阵,y 为与A 有相同行数的列向量。