计算方法的课后答案解析

合集下载

计算方法习题二答案

计算方法习题二答案

计算方法习题二答案习题二1、利用二分法求方程f(x)=x3-2x-5=0,在2,3内根的近似值,并指出误差。

解:f(2)=-1<0 f(3)=19>0 f(2).f(3)<0f’(x)=3x2-2 在x∈2,3f’(x) >0所以在1,2上必仅有一根x=2 f(2)=-1 -x=3 f(3)=16 +x=2.5 f(2.5)=5.625 +x=2.25 f(2.25)=1.890625 +x=2.125 f(2.125) +x=2.0625 f(2.0625) -x=2.09375 f(2.09375) -x=2.109375 f(2.109375) +x=2.1015625 f(2.1015625) +所以x=2.109375+2010156252=2.097656252、证明方程1-x-sinx=0在0,1内有一个根,使用二分法求误差不大于12×10?4的根。

解:令f(x)=1-x-sinxf(0)=1f(1)=-sin1f(0).f(1)<0f’(x)=-1-cosx<0在0,1恒成立所以1-x-sinx=0在0,1内恒有一个根n≥ln1?0?ln?(12×10?4)ln2-1≈13.289所以n=14n a n b n x n+1f(x n+1)符号0 0 1 0.5 +1 0.5 1 0.75 +2 0.875 1 0.9375 +..143、能不能用迭代法求解下列方程,若不能时,将方程改写成能用迭代法的形式。

(1、)x=(cosx+sinx)/4 (2)x=4-2x解:(1、)f(x)=x=(cosx+sinx)/4f’(x)=?sinx+cosx4<1对x任何数恒成立所以可用迭代法设x0=0,则x1=0.25x2=0.2511x 3=0.2511所以x=0.251(2、)f(x)=4-2xf’(x)=x.2x ?1<0在x 为任意数不恒成立所以不能用迭代法令x=log 2(4?x)x 0=0x 1=2x 2=1x 3= |φ‘(x)|=|-14?x 1ln 2|对x ∈(1,2)<124、为求方程x 3-x 2-1=0在x 0=1.5附近的一个根,设将方程改写成下列等价形式,并建立相应的迭代公式。

《数值计算方法》 课后题 答案(曾金平)湖南大学

《数值计算方法》  课后题   答案(曾金平)湖南大学
1 1 1 1
= 0.105 × 10 2 − 0.144 × 10 2 + 0.657 × 10 1 − 1
= 0.167 × 101
g ( 2.19) = ((−0.81) × 0.219 × 101 + 3) × 0.219 × 101 − 1
= 0.123 × 10 1 × 0.219 × 10 1 − 1 = 0.169 × 101 1 1 即 f ( x) = 0.167 × 10 , g ( x) = 0.169 × 10
而当 x = 2.19 时 x 3 − 3 x 2 + 3 x − 1 的精确值为 1.6852,故 g ( x ) 的算法较正确。 8.按照公式计算下面的和值(取十进制三位浮点数计算) : (1)
∑3
i =1
6
1
i 6
;(2)
∑3
i =6
1
1
i

解: (1)
∑3
i =1 1
1
i
1 1 1 1 1 1 = + 2 + 3 + 4 + 5 + 6 = 0.333 + 0.111 + 0.037 + 0.012 + 0.004 + 0.001 3 3 3 3 3 3
xk
1 1.5 1.25 1.125 1.0625 1.09375 1.109375 1.1171875 1.11328125 1.115234375 1.1142578125 1.11376953125
f ( xk )
-0.1585 0.4962 0.1862 0.015051 -0.0718 -0.02835 -0.00664 0.004208 -0.001216 0.001496 0.001398 -0.000538

计算方法习题第一、二章答案

计算方法习题第一、二章答案

第一章 误差1 问,,722分别作为π的近似值各具有几位有效数字?分析 利用有效数字的概念可直接得出。

解 π= 592 65… 记x 1=,x 2=,x 3=722.由π- x 1= 59…= 40…知3411110||1022x π--⨯<-≤⨯ 因而x 1具有4位有效数字。

由π- x 2= 59…= 59…知2231021||1021--⨯≤-<⨯x π因而x 2具有3位有效数字。

由π-722= 59 … 85…= 26…知231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。

2 已知近似数x*有两位有效数字,试求其相对误差限。

分析 本题显然应利用有效数字与相对误差的关系。

解 利用有效数字与相对误差的关系。

这里n=2,a 1是1到9之间的数字。

%5101211021|*||*||)(|1211*=⨯⨯≤⨯≤-=+-+-n ra x x x x ε3 已知近似数的相对误差限为%,问x*至少有几位有效数字?分析 本题利用有效数字与相对误差的关系。

解 a 1是1到9间的数字。

1112*10)1(2110)19(21102110003%3.0)(--⨯+≤⨯+⨯=⨯<=a x r ε 设x*具有n 位有效数字,令-n+1=-1,则n=2,从而x*至少具有2位有效数字。

4 计算,问要取几位有效数字才能保证相对误差限不大于%。

分析 本题应利用有效数字与相对误差的关系。

解 设取n 位有效数字,由=…,故a 1=9。

411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式411101021-+-≤⨯n a 知取n=4即可满足要求。

5 计算76017591-,视已知数为精确值,用4位浮点数计算。

解 =-76017591 8×10-2-0.131 6×10-2=×10-5结果只有一位有效数字,有效数字大量损失,造成相对误差的扩大,若通分后再计算:56101734.0105768.01760759176017591-⨯=⨯=⨯=- 就得到4位有效数字的结果。

数值计算方法课后习题答案

数值计算方法课后习题答案

第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。

[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%,求n x 的相对误差。

[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x 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 。

数值计算方法课后习题答案

数值计算方法课后习题答案

习题一1.设x >0相对误差为2%4x 的相对误差。

解:由自变量的误差对函数值引起误差的公式:(())(())'()()()()f x xf x f x x f x f x δδ∆=≈得 (1)()f x =11()()*2%1%22x x δδδ≈===;(2)4()f x x =时444()()'()4()4*2%8%x x x x x x δδδ≈===4.计算正方形面积时,若要求面积的允许相对误差为1%,测量边长所允许的相对误差限为多少?解:设该正方形的边长为x ,面积为2()f x x =,由(())(())'()()()()f x xf x f x x f x f x δδ∆=≈ 解得(())()()'()f x f x x xf x δδ≈=2(())(())22f x x f x x xδδ= =0.5%5.下面计算y 的公式哪个算得准确些?为什么?(1)已知1x <<,(A )11121x y x x -=-++,(B )22(12)(1)x y x x =++; (2)已知1x >>,(A)y =,(B)y =(3)已知1x <<,(A )22sin x y x=,(B )1cos 2xy x -=;(4)(A)9y =(B)y =解:当两个同(异)号相近数相减(加)时,相对误差可能很大,会严重丧失有效数字;当两个数相乘(除)时,大因子(小除数)可能使积(商)的绝对值误差增大许多。

故在设计算法时应尽量避免上述情况发生。

(1)(A )中两个相近数相减,而(B )中避免了这种情况。

故(B )算得准确些。

(2)(B )中两个相近数相减,而(A )中避免了这种情况。

故(A )算得准确些。

(3)(A )中2sin x 使得误差增大,而(B )中避免了这种情况发生。

故(B )算得准确些。

(4)(A )中两个相近数相减,而(B )中避免了这种情况。

计算方法 课后习题答案

计算方法 课后习题答案
解:因为第一列中10最大,因此把10作为列主元素
得到方程组
3。举例说明一个非奇异矩阵不一定存在LU分解。
例如:设
与题设相矛盾,所以一个非奇异矩阵不一定存在LU分解。
4。下列矩阵能否分解为LU(其中L为单位下三角矩阵,U为上三角矩阵)?若能分解,那么分解是否唯一?
解:
设 B可以进行LU分解,则B=
计算得
其中。 。
解:(1)由题意,可设 ,由Lagrange插值余项公式得
(2)由(1)式可知,
15.给定数据表:
1
0
2
3
构造出函数 的差商表,并写出它的三次 插值多项式.
解:利用Newton插值公式:
先作出差商表
一阶差商
二阶差商
三阶差商
0
1
3
1
3/2
13/4
1/2
2
0
3
1/6
1/3
3
2
5/3
-2/3
-5/3
证明:据题4可知,
令 ,则有 。注意到
(证明见王能超数值简明教程145页题6)
令 即有 。
9.已知 ,求差商 和 。
解:根据差商与微商的关系,有
10.已知 互异,求 。其中 。(此题有误。)(见王能超《教程》P149-题2)
解:因为 ,则
由差商性质 可知,
11.设首项系数为1的n次式 有n个互异的零点 ,证明
解:1)用梯形公式有:
事实上,
2)Simpson公式
事实上,
3)由Cotes公式有:
事实上,
2.证明Simpson公式 具有三次代数精度。
证明:
而当 时
左侧:
右侧:

计算方法第二版课后练习题含答案

计算方法第二版课后练习题含答案

计算方法第二版课后练习题含答案前言本文将为大家提供计算方法第二版课后练习题的答案,旨在帮助读者更好地学习和掌握计算方法的知识。

本文全部内容均为作者整理,尽可能保证每一题的答案正确性。

读者可以借助本文的答案,检验自己的练习成果,加强对计算方法知识的理解和掌握程度。

同时,读者也应该注意切勿直接复制答案,本文的答案仅供参考,希望读者能够通过自己的思考和探索,获得更深层次的学习感悟。

第一章引论1.1 计算方法的基本概念和思想练习题 1写出计算方法的三要素,并分别简要解释。

答案计算方法的三要素为:模型、算法、误差分析。

•模型:计算方法所涉及的实际问题所对应的数学模型,是解决问题的基础;•算法:根据模型,构造相应的计算程序,即算法;•误差分析:计算结果与实际应用中所需的精度之间的差异,称为误差。

误差分析是对计算结果质量的保障。

1.2 算法的误差练习题 2写出二分法算法,并解释其误差。

答案算法:function binarySearch(a, target) {let low = 0;let high = a.length - 1;while (low <= high) {let midIndex = Math.floor((low + high) / 2);let midValue = a[midIndex];if (midValue === target) {return midIndex;} else if (midValue < target) {low = midIndex + 1;} else {high = midIndex - 1;}}return -1;}误差:二分法算法的误差上界为O(2−k),其中k为迭代次数。

在二分法被成功应用时,k取决于与目标值x的距离,即 $k=\\log _{2}(\\frac{b-a}{\\epsilon})$,其中[a,b]是区间,$\\epsilon$ 是目标值的精度。

最优化计算方法课后习题答案----高等教育出社。施光燕

最优化计算方法课后习题答案----高等教育出社。施光燕

习题二包括题目: P36页 5(1)(4)5(4)习题三包括题目:P61页 1(1)(2); 3; 5; 6; 14;15(1) 1(1)(2)的解如下3题的解如下5,6题14题解如下14. 设22121212()(6)(233)f x x x x x x x =+++---, 求点在(4,6)T-处的牛顿方向。

解:已知 (1)(4,6)T x=-,由题意得121212212121212(6)2(233)(3)()2(6)2(233)(3)x x x x x x x f x x x x x x x x +++-----⎛⎫∇= ⎪+++-----⎝⎭∴ (1)1344()56g f x -⎛⎫=∇=⎪⎝⎭21212122211212122(3)22(3)(3)2(233)()22(3)(3)2(233)22(3)x x x x x x x f x x x x x x x x +--+--------⎛⎫∇= ⎪+--------+--⎝⎭∴ (1)2(1)1656()()564G x f x --⎛⎫=∇=⎪-⎝⎭(1)11/8007/400()7/4001/200G x --⎛⎫= ⎪--⎝⎭∴ (1)(1)11141/100()574/100d G x g -⎛⎫=-=⎪-⎝⎭15(1)解如下15. 用DFP 方法求下列问题的极小点(1)22121212min 353x x x x x x ++++解:取 (0)(1,1)T x=,0H I =时,DFP 法的第一步与最速下降法相同2112352()156x x f x x x ++⎛⎫∇= ⎪++⎝⎭, (0)(1,1)T x =,(0)10()12f x ⎛⎫∇= ⎪⎝⎭(1)0.07800.2936x -⎛⎫= ⎪-⎝⎭, (1)1.3760() 1.1516f x ⎛⎫∇= ⎪-⎝⎭以下作第二次迭代(1)(0)1 1.07801.2936x x δ-⎛⎫=-= ⎪-⎝⎭, (1)(0)18.6240()()13.1516f x f x γ-⎛⎫=∇-∇= ⎪-⎝⎭0110111011101T T T TH H H H H γγδδδγγγ=+- 其中,111011126.3096,247.3380T T TH δγγγγγ===111.1621 1.39451.3945 1.6734Tδδ⎛⎫= ⎪⎝⎭ , 01101174.3734113.4194113.4194172.9646T TH H γγγγ⎛⎫== ⎪⎝⎭所以10.74350.40560.40560.3643H -⎛⎫= ⎪-⎝⎭(1)(1)1 1.4901()0.9776dH f x -⎛⎫=-∇= ⎪⎝⎭令 (2)(1)(1)1xx d α=+ , 利用 (1)(1)()0df x d d αα+=,求得 10.5727α=-所以 (2)(1)(1)0.77540.57270.8535xx d⎛⎫=-= ⎪-⎝⎭ , (2)0.2833()0.244f x ⎛⎫∇= ⎪-⎝⎭以下作第三次迭代(2)(1)20.85340.5599x x δ⎛⎫=-= ⎪-⎝⎭ , (2)(1)2 1.0927()()0.9076f x f x γ-⎛⎫=∇-∇= ⎪⎝⎭22 1.4407T δγ=- , 212 1.9922T H γγ=220.72830.47780.47780.3135T δδ-⎛⎫=⎪-⎝⎭1221 1.39360.91350.91350.5988T H H γγ-⎛⎫= ⎪-⎝⎭所以22122121222120.46150.38460.38460.1539T T T TH H H H H δδγγδγγγ-⎛⎫=+-= ⎪-⎝⎭(2)(2)20.2246()0.1465d H f x ⎛⎫=-∇= ⎪-⎝⎭令 (3)(2)(2)2xxdα=+ , 利用(2)(2)()0df x d d αα+=,求得 21α=所以 (3)(2)(2)11x x d ⎛⎫=+=⎪-⎝⎭, 因为 (3)()0f x ∇=,于是停止 (3)(1,1)T x =-即为最优解。

计算方法第二章作业答案参考

计算方法第二章作业答案参考
main()
{
double mid=0,x1=0.3,x2=0.4;
while(1)
{mid=(x1+x2)/2;
if(f(mid)<0) {x2=mid;printf("x1=%10f,x2=%10f\n",x1,x2);}
else if(f(mid)>0) {x1=mid;printf("x1=%10f,x2=%10f\n",x1,x2);}
main()
{
float x1=1.5,a=0,b=0;
int i;
for(i=0;i<20;i++)
{
b=x1;
a=f(x1);
x1=f(a);
x1=(x1*b-a*a)/(x1-2*a+b);
printf("%10f",x1);
}
printf("\n因此,x=%.4f\n",x1);
}
9.用牛顿法求方程 在 附近的根,要求准确到小数点后第3位
main()
{
float x1=1.3;
int i;
for(i=0;i<20;i++)
{
x1=f(x1);
printf("%10f",x1);
}
printf("\n因此,x=%f\n",x1);
}
(3) ,
#include<stdio.h>
#include<math.h>
#define f(x) pow(x-1,-1*1.0/2)
main()
{
float x1=1.3;

计算方法教程(第2版)习题答案

计算方法教程(第2版)习题答案

《计算方法教程(第二版)》习题答案第一章 习题答案1、浮点数系),,,(U L t F β共有 1)1()1(21++---L U t ββ 个数。

3、a .4097b .62211101110.0,211101000.0⨯⨯c .6211111101.0⨯4、设实数R x ∈,则按β进制可表达为:,1,,,3,2,011)11221(+=<≤<≤⨯++++++±=t t j jd d l t t d t t d dd x βββββββ按四舍五入的原则,当它进入浮点数系),,,(U L t F β时,若β211<+t d ,则 l tt d dd x fl ββββ⨯++±=)221()(若 β211≥+t d ,则 l tt d d d x fl ββββ⨯+++±=)1221()(对第一种情况:t l lt l t t d x fl x -++=⨯≤⨯+=-βββββ21)21(1)()(11对第二种情况:t l lt l t t d x fl x -++=⨯≤⨯--=-ββββββ21)21(1)(11就是说总有: tl x fl x -≤-β21)( 另一方面,浮点数要求 β<≤11d , 故有l x ββ1≥,将此两者相除,便得t x x fl x -≤-121)(β 5、a . 5960.1 b . 5962.1 后一种准确6、最后一个计算式:00025509.0原因:避免相近数相减,避免大数相乘,减少运算次数7、a .]!3)2(!2)2(2[2132 +++=x x x yb .)21)(1(22x x x y ++=c .)11(222-++=x x x yd . +-+-=!2)2(!6)2(!4)2(!2)2(2642x x x x y e .222qp p q y ++=8、01786.098.5521==x x9、 m )10(m f - 1 233406.0- 3 20757.0- 5 8.07 710计算宜采用:])!42151()!32141()!22131[()(2432 +⨯-+⨯-+⨯--=x x x f第二章 习题答案1、a .Tx )2,1,3(= b .Tx )1,2,1,2(--= c .无法解 2、a .与 b .同上, c .T T x )2188.1,3125.0,2188.1,5312.0()39,10,39,17(321---≈---=7、a .⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛---14112111473123247212122123211231321213122 b . ⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫⎝⎛----333211212110211221213231532223522121⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=111211212130213221219、T x )3415.46,3659.85,1220.95,1220.95,3659.85,3415.46(1= T x )8293.26,3171.7,4390.2,4390.2,3171.7,8293.26(2= 10、T LDL 分解:)015.0,579.3,9.1,10(diag D =⎪⎪⎪⎪⎪⎭⎫⎝⎛=16030.07895.05.018947.07.019.01L Cholesky 分解⎪⎪⎪⎪⎪⎭⎫⎝⎛=1225.01408.10833.15811.18918.12333.12136.23784.18460.21623.3G 解:)1,1,2,2(--=x 12、16,12,1612111===∞A A A611,4083.1,61122212===∞A A A2)(940)()(12111===∞A Cond A Cond A Cond 524)(748)()(22221===∞A Cond A Cond A Cond⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=--180.0000 180.0000- 30.0000 180.0000- 192.0000 36.0000- 30.0000 36.0000- 9.0000 ,0.0139 0.1111- 0.0694- 0.1111- 0.0556 0.1111- 0.0694- 0.1111- 0.0139 1211A A1151.372,1666.0212211==--A A15、 1A :对应 Seidel Gauss - 迭代收敛,Jacobi 迭代不收敛; 2A :对应 Seidel Gauss - 迭代收敛,Jacobi 迭代不收敛; 3A :对应 Seidel Gauss - 迭代收敛,Jacobi 迭代收敛;第三章 习题答案1、Lagrange 插值多项式:)80.466.5)(20.366.5)(70.266.5)(00.166.5()80.4)(20.3)(70.2)(00.1(7.51)66.580.4)(20.380.4)(70.280.4)(00.180.4()66.5)(20.3)(70.2)(00.1(3.38)66.520.3)(80.420.3)(70.220.3)(00.120.3()66.5)(80.4)(70.2)(00.1(0.22)66.570.2)(80.470.2)(20.370.2)(00.170.2()66.5)(80.4)(20.3)(00.1(8.17)66.500.1)(80.400.1)(20.300.1)(70.200.1()66.5)(80.4)(20.3)(70.2(2.14)(4--------⨯+--------⨯+--------⨯+--------⨯+--------⨯=x x x x x x x x x x x x x x x x x x x x x LNewton 插值多项式:)80.4)(20.3)(70.2)(00.1(21444779.0)20.3)(70.2)(00.1(527480131.0)70.2)(00.1(855614973.2)00.1(117647059.22.14)(4----+------+-+=x x x x x x x x x x x N2、设)(x y y =,其反函数是以y 为自变量的函数)(y x x =,对)(y x 作插值多项式:)1744.0)(1081.0)(4016.0)(7001.0(01253.0)1081.0)(4016.0)(7001.0(01531.0)4016.0)(7001.0(009640.0)7001.0(3350.01000.0)(----+---+--+--=y y y y y y y y y y y N3376.0)0(=N 是0)(=x y 在]4.0,3.0[中的近似根。

计算方法课后习题答案第四章作业

计算方法课后习题答案第四章作业

(五)课后习题4.1 对于积分⎰-aadx x f )(,以a x x a x ==-=210,0,为节点,构造形如⎰-++≈aax f A x f A x f A dx x f )()()()(221100的插值型求积公式,并讨论所得公式的代数精度。

解答:⎰⎰--=------=----=aa a a a dx a a a a x x dx x x x x x x x x A 31))(0())(0())(())((2010210⎰⎰--=-+-+=----=aa a a a dx a a a x a x dx x x x x x x x x A 34)0)(0())(())(())((2101201⎰⎰--=-+-+=----=aa a a a dx a a a x a x dx x x x x x x x x A 31)0)(()0)(())(())((1202102易知为Simpson 公式,因此代数精度为34.2 确定 下列求积公式中的待定参数,使其代数精度尽量高,并指出所得公式的代数精度。

(1)⎰++≈2210)2()1()0()(f A f A f A dx x f(2)⎰-⋅++≈hh f f h h f f hdx x f 0''2)]()0([)]()0([2)(α解答:(1)令2,,1)(x x x f =,假定求积公式均准确成立,从而有: ⎰++==202102A A A dx 21022102⋅+⋅+⋅==⎰A A A xdx22212022210038⋅+⋅+⋅⋅==⎰A A A dx x 解以上三元线性方程组从得:34,31120===A A A ,显然仍为Simpson 公式,因此代数精度为3(2)求积公式中只含一个待定参数α,当x x f ,1)(=时,有 ⎰++=hh dx 00]11[2,⎰-++=h h h hxdx 02)11(]0[2α故令2)(x x f =时求积公式准确成立,即⎰-⨯++=hh h h h dx x 0222]202[]0[2α,解得121=α将3)(x x f =代入上述确定的求积公式,有:⎰-++=hh h h h dx x 02233]30[12]0[2,这说明求积公式至少有3次代数精度,再令 4)(x x f =,代入求积公式时有:⎰-++≠hh h h h dx x 03244]40[12]0[2故所建求积公式为⎰-++≈hh f f h h f f h dx x f 0''2)]()0([2)]()0([2)(4.3 对于xxx f sin )(=,利用下表数据,计算8,4=n 时的复合梯形公式84,T T ,以及4=n 复合Simpson 公式4S 的值。

计算方法练习题集和答案解析

计算方法练习题集和答案解析

练习题与答案练习题一练习题二练习题三练习题四练习题五练习题六练习题七练习题八练习题答案练习题一一、是非题1.*x=–12.0326作为x的近似值一定具有6位有效数字,且其误差限£41021-⨯。

( )2.对两个不同数的近似数,误差越小,有效数位越多。

()3.一个近似数的有效数位愈多,其相对误差限愈小。

()4.用212x-近似表示c o s x产生舍入误差。

() 5. 3.14和 3.142作为π的近似值有效数字位数相同。

() 二、填空题1.为了使计算()()2334912111yx x x=+-+---的乘除法次数尽量少,应将该表达式改写为;2.*x=–0.003457是x舍入得到的近似值,它有位有效数字,误差限为,相对误差限为;3.误差的来源是;4.截断误差为;5.设计算法应遵循的原则是。

三、选择题1.*x=–0.026900作为x的近似值,它的有效数字位数为( ) 。

(A) 7; (B) 3;(C) 不能确定 (D) 5.2.舍入误差是( )产生的误差。

(A) 只取有限位数 (B) 模型准确值与用数值方法求得的准确值(C) 观察与测量 (D) 数学模型准确值与实际值3.用 1+x近似表示e x所产生的误差是( )误差。

(A). 模型 (B). 观测 (C). 截断 (D). 舍入4.用s *=21g t 2表示自由落体运动距离与时间的关系式 (g 为重力加速度),s t 是在时间t 的实际距离,则s t s *是( )误差。

(A). 舍入 (B). 观测 (C). 模型 (D). 截断 5.1.41300作为2的近似值,有( )位有效数字。

(A) 3; (B) 4; (C) 5; (D) 6。

四、计算题1. 3.142,3.141,227分别作为π的近似值,各有几位有效数字?2. 设计算球体积允许的相对误差限为1%,问测量球直径的相对误差限最大为多少?3. 利用等价变换使下列表达式的计算结果比较精确:(1)1||,11211<<+-++x x x x , (2) 1||1112<<+⎰+x dt t x x (3) 1||,1<<-x e x , (4) 1)1ln(2>>-+x x x4.真空中自由落体运动距离s 与时间t 的关系式是s =21g t 2,g 为重力加速度。

计算方法答案王能超

计算方法答案王能超

计算方法答案王能超【篇一:计算方法习题集及实验指导书】s=txt>计算机科学与技术系檀明2008-02-10课程性质及目的要求(一)课程性质自计算机问世以来,科学计算一直是计算机应用的一个重要领域,数值计算方法是解决各种复杂的科学计算问题的理论与技术的基础。

《计算方法》课程讨论用于科学计算中的一些最基本、最常用的算法,不但具有数学的抽象性与严密的科学性的特点,而且具有应用的高度技术性的特点。

它对于培养从事计算机应用的科技人才有着重要的作用,是计算机应用专业(本科段)的一门重要的技术基础课程。

(二)目的要求通过本课程的学习和上机实验,了解用计算机解决科学计算问题的方法特点,掌握计算方法中的一些基本概念、基本公式和相应的算法流程,提高根据算法描述设计高级语言程序并进行验证的技能。

在学习过程中,应注重理解和应用,在搞清基本原理和基本概念的基础上,通过习题、编程和上机等环节,巩固和加深已学的内容,掌握重要的算法及其应用。

注重理论与算法的学习和应用相结合,强调编程及上机计算的技能培养,是本课程不同于一般数学课程的重要特点。

(三)学习方法指导1.循序渐进逐章学习本课程从第二章开始,每章都讨论一个大类的算法。

虽然各算法是相对独立的,但是也存在相互联系与前后继承的关系。

前面的概念和算法学好了,后面的内容也就容易学,越学越感到容易。

前面的内容没有学好,后面就会感到难学,甚至会出现越来越感到困难、失去学习信心的情况。

2.稳扎稳打融会贯通学习要扎实、要讲求实效。

每一个重要的概念和公式,都会搞清楚,做到融会贯通。

只有这样,才能取得学习的学习效果。

3.多学练勤做习题教材及本习题集中的每一章都附有适量的习题,可以帮助考生巩固和加深理解所学的知识,提高解题能力。

因此,在学习过程中,应当适合习题进行思考,应当尽可能多做习题,遇到某些不会做的题,应三思之后再请老师给予提示。

4.抓住特点前后联系本课程只讲了五大类算法。

每类算法都是针对一类特定的计算问题,都有其自身的特点。

数值计算方法课后习题答案(李庆扬等)

数值计算方法课后习题答案(李庆扬等)

第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。

[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%,求n x 的相对误差。

[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x 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 。

计算方法 课后习题答案

计算方法 课后习题答案

计算方法课后习题答案计算方法课后习题答案计算方法是一门重要的学科,它涉及到数值计算、算法设计和数据处理等方面的内容。

在学习计算方法的过程中,课后习题是不可或缺的一部分。

通过解答习题,我们可以巩固所学的知识,提高自己的计算能力。

下面是一些计算方法课后习题的答案,希望对大家的学习有所帮助。

1. 矩阵的转置矩阵的转置是将矩阵的行和列互换得到的新矩阵。

对于一个m×n的矩阵A,它的转置记作A^T。

转置后的矩阵A^T的行数和列数分别为原矩阵A的列数和行数。

例如,对于一个3×2的矩阵A,它的转置A^T是一个2×3的矩阵。

2. 矩阵的加法和减法矩阵的加法和减法是对应位置上的元素进行相加或相减得到的新矩阵。

对于两个相同大小的矩阵A和B,它们的和记作A+B,差记作A-B。

加法和减法的运算规则是相同位置上的元素进行相应的运算。

3. 矩阵的乘法矩阵的乘法是指将两个矩阵相乘得到一个新矩阵的运算。

对于两个矩阵A和B,它们的乘积记作AB。

矩阵乘法的运算规则是矩阵A的行与矩阵B的列进行相乘,并将结果相加得到新矩阵的对应位置上的元素。

4. 矩阵的逆矩阵的逆是指对于一个可逆矩阵A,存在一个矩阵B,使得AB=BA=I,其中I是单位矩阵。

如果一个矩阵A存在逆矩阵,则称其为可逆矩阵或非奇异矩阵。

求解矩阵的逆可以使用伴随矩阵和行列式的方法。

5. 线性方程组的求解线性方程组是指由一组线性方程组成的方程组。

求解线性方程组的方法有很多,包括高斯消元法、LU分解法、迭代法等。

其中,高斯消元法是一种常用的求解线性方程组的方法,它通过消元和回代的过程,将线性方程组转化为上三角形矩阵或对角矩阵,从而求解出方程组的解。

6. 数值积分的方法数值积分是指通过数值计算的方法来求解定积分的近似值。

常用的数值积分方法包括梯形法则、辛普森法则和龙贝格法则等。

这些方法都是基于将定积分转化为离散求和的形式,通过计算离散点上的函数值来估计定积分的近似值。

计算方法-刘师少版第二章课后习题完整答案

计算方法-刘师少版第二章课后习题完整答案

0 < λf ′(x) < 2
− 2 < −λf ′(x) < 0
−1 < 1 − λf ′(x) < 1
1 − λf ′(x) < 1
即 ϕ ′(x) < 1 ,所以 xk+1 = ϕ (xk ) = xk − λf (xk ) 收敛于 f (x) =0 的根。
2.7 试用牛顿迭代法导出下列各式的迭代格式:
应的迭代公式:
(1)x
=1+
1 x2
,迭代公式
xk
+1
=1+
1
x
2 k
(2)x3 = 1 + x 2 ,迭代公式 xk+1 = 3 1 + xk2
(3) x 2
=
x
1 −
1
,迭代公式
xk
+1
=
1 xk −1
(4) x = x3 − 1 ,迭代公式
xk+1 = xk3 − 1
试分析每种迭代公式的收敛性,并选取一种收敛迭代公式求出具有四位有效数字的近似
x8 = 1.4656344
x9 = 1.4656000
x9
− x8
≤ 1 ×10−4 , 2
x9 = 1.4656000
2.5 对于迭代函数ϕ (x) = x + C(x 2 − 2) ,试讨论:
(1) 当 C 取何值时, xk+1 = ϕ (xk ), (k = 0,1,2,L) 产生的序列 {xk }收敛于 2 ;
6 6x2
63
ϕ ′(3 a ) == 5 − a (3 a )−3 = 5 − 1 = 1 ≠ 0
63

计算方法引论课后答案.

计算方法引论课后答案.

计算⽅法引论课后答案.第⼀章误差1. 试举例,说明什么是模型误差,什么是⽅法误差.解: 例如,把地球近似看为⼀个标准球体,利⽤公式24A r π=计算其表⾯积,这个近似看为球体的过程产⽣的误差即为模型误差.在计算过程中,要⽤到π,我们利⽤⽆穷乘积公式计算π的值:12222...q q π=?其中112,3,...n q q n +?=??==?? 我们取前9项的乘积作为π的近似值,得3.141587725...π≈这个去掉π的⽆穷乘积公式中第9项后的部分产⽣的误差就是⽅法误差,也成为截断误差.2. 按照四舍五⼊的原则,将下列各数舍成五位有效数字:816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236 3. 下列各数是按照四舍五⼊原则得到的近似数,它们各有⼏位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位三位六位四位4. 若1/4⽤0.25表⽰,问有多少位有效数字? 解: 两位5. 若 1.1062,0.947a b ==,是经过舍⼊后得到的近似值,问:,a b a b +?各有⼏位有效数字?解: 已知4311d 10,d 1022a b --()433211110100.551010222d a b da db da db ----+=+≤+=?+?=?所以a b +有三位有效数字;因为0.1047571410a b ?=?,()43321110.94710 1.1062100.600451010222所以a b ?有三位有效数字.6. 设120.9863,0.0062y y ==,是经过舍⼊后作为12,x x 的近似值.求1211,y y 的计算值与真值的相对误差限及12y y ?与真值的相对误差限. 解: 已知-4-41112221211d ,d ,d =10,d 1022x y x x y x x x =+=+?=?, ()44111111110d d 12dr dr 0.50100.9863x xx x x y --==≈=≈? ???;()42222222110d d 12dr dr 0.81100.0062x xx x x y --==≈=≈? ???;()()()4221212dr dr dr 0.50100.81100.8210x x x x ---?=+≈?+?≈?.7. 正⽅形的边长约为100cm,应该怎样测量,才能使其⾯积的误差不超过1cm 2.解: 设正⽅形⾯积为S,边长为a,则S=a 2.所以要使:2d d 2d 1s a a a ==≤,则要求211d 0.5102200a a -≤==?.所以边长的误差不能超过20.510-?cm.8. ⽤观测恒星的⽅法求得某地维度为4502'''o(读到秒),试问:计算sin ?将有多⼤误差?解: ()()1d sin cos d cos 45022*''?'''==o.9 . 真空中⾃由落体运动距离s 与时间的关系由公式212s gt =确定,g 是重⼒加速度.现在假设g 是准确的,⽽对t 的测量有0.1s ±的误差,证明t 增加时,距离的绝对误差增加⽽相对误差却减⼩.证明: 因为:221d d d d d d d ;2.122s gt t gt t t s gt gt t s s t gt ??d s 与t 成正⽐,d s s与t 成反⽐,所以当d t 固定的时候, t 增加时,距离的绝对误差增加⽽相对误差却减⼩.10. 设0x >,x 的相对误差为δ,求ln x 的绝对误差. 解: 已知d x x δ=,所以ln x 的绝对误差()d d ln x x x δ==.11. 设x 的相对误差为%α,求nx 的相对误差.解: 1d d d %n n n n x nx x n xn x x xα-===.12. 计算球的体积,为了使相对误差限为1%,问度量半径R 时允许的相对误差限如何? 解: 已知34 3V R π=,设()d dr R R a R ==,则要使得 ()()3d dr dln d ln 3d ln 3d ln 3dr 31%V V V R R R R a V ========,则11%3a =?.第⼆章插值法与数值微分1.设y =在100,121,144x =三处的值是很容易求得的,试以这三个点建⽴y =的⼆次插值多项式,,且给出误差估计.⽤其中的任意两点,构造线性插值函数,⽤得到的三个线性插值函数,,并分析其结果不同的原因.解: 已知012012100,121,144;10,11,12x x x y y y ======,建⽴⼆次Lagrange 插值函数可得:()()()()21211441001441011100121100144121100121144121100 12144121144100x x x x L x x x ----= +------+--()211510.7228L ≈=.误差()()()()()()2012012,,,,3!f R x x x x x x x x x x ξξξ'''=---∈,所以20.00065550.001631R <<利⽤前两个节点建⽴线性插值函数可得:()()()()()11211001011100121121100x x L x --=+--()111510.7143L ≈=.利⽤后两个节点建⽴线性插值可得:()()()()()11441211112121144144121x x L x --=+--()111510.7391L ≈=.利⽤前后两个节点建⽴线性插值可得:()()()()()21441001012100144144100x x L x --=+()111510.6818L ≈=.,⼆次插值⽐线性插值效果好,利⽤前两个节点的线性插值⽐其他两个线性插值效果好.此说明,⼆次插值⽐线性插值效果好,插⽐外插效果好.2. 利⽤(2.9)式证明()()()0121001max ,8x x x x x R x f x x x x ≤≤-''≤≤≤证明: 由(2.9)式()()()()0101,2!f R x x x x x x x ξξ''=--<<当01x x x <<时,()()01max x x x f f x ξ≤≤''''≤,()()()01201101max 4x x x x x x x x x ≤≤--≤- 所以()()()0121001max ,8x x x x x R x f x x x x ≤≤-''≤≤≤3. 若()0,1,...,j x n 为互异节点,且有()()()()()()()()()011011............j j n j jj j j j j n x x x x x x x x l x xx x x x x x x -+-+----=证明()0,0,1,...,nk kj j j x l x xk n =≡=∑证明: 由于() 1 ;0 .j i ij i j l x i j δ=?==?≠? 且()0nk j j j x l x =∑和kx都为k 次多项式,⽽且在k+1个不同的节点处的函数值都相同0,1,...,k n =, 所以马上有()0,0,1,...,nk kj j j x l x xk n =≡=∑.4. 设给出sin x 在[],ππ-上的数值表,⽤⼆次插值进⾏计算,若希望截断误差⼩于5 10-,问函数表的步长最⼤能取多少? 解: 记插值函数为p(x),则()()()()()11sin sin 3!i i i x p x x x x x x x ξ-+'''-=--- 所以()()()()11cos max sin 3!i i i x x p x x x x x x ππξ-+-≤≤--=---()()()[]3112,0,2i g x th h t t t t -+=--∈⼜()()()[]12,0,2t t t t t ?=--∈的最⼤值为10.3849??= ?,所以有 350.3849max sin 106x x p h ππ--≤≤-≤< 所以 0.0538h ≤.5. ⽤拉格朗⽇插值和⽜顿插值找经过点()()()()3,1,0,2,3,2,6,10---的三次插值公式. 解: Lagrange 插值函数:()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()12302330101020310121301301223202123303132 31033101622731033 .2781/5x x x x x x x x x x x x L x y y x x x x x x x x x x x x x x x x x x x x x x x x y y x x x x x x x x x x x x x x x x x x x x x x x x ------= +------------++--------+--=++-+-++⽜顿插值: ⾸先计算差商3 10 2 13 2 1.333 0.38896 104 0.8889 0.1420-----()()()()()3130.38893 1.142033.N x x x x x x x =-++-+++-也可以利⽤等距节点构造,⾸先计算差分。

最优化计算方法课后习题集答案解析

最优化计算方法课后习题集答案解析
(1)
解:取 , 时,DFP法的第一步与最速下降法相同
, ,

以下作第二次迭代

其中,

所以
令 , 利用 ,求得
所以 ,
以下作第三次迭代


所以
令 , 利用 ,求得
所以 , 因为 ,于是停止
即为最优解。
习题四
包括题目: P95页 3;4;8;9(1);12选做;13选做
3题解如下
3.考虑问题 ,其中
X1,x2,x3≥0 (3)
求出点(1,1,0)处的一个下降可行方向.
解:首先检查在点(1,1,0)处哪些约束为有效约束。检查易知(1),X3≥0为有效约束。设所求可行方向d=(d1,d2,d3)T。根据可行方向d的定义,应存在a>0,使对∀t∈(0,a)能有
X+td=(1+td1,1+td2,0+td3)T
(1)
s.t.
(2)
s.t.
(1)解:非线性规划的K-T条件如下:
(1)
(2)
(3)
再加上约束条件 (4)
为求出满足(1)~(4)式的解,分情况考虑:
①若(4)式等号不成立,即 ,那么由(2)式得 ,将 代入(1)式解得 , ,所得值不满足 的条件,故舍去。
②若(4)式等号成立,由(1)式可以解得 , ,代入(4)式有:
JBi
1
2
3
4
5
6
7
8
9
di0
1
1
0
-5/6
-1/6
1
10/6
4
0
0
38/6
2
0
1
-9/6

沪教版三年级数学下册整十数除两、三位数的计算方法含答案及解析

沪教版三年级数学下册整十数除两、三位数的计算方法含答案及解析

数学测试卷小学基础知识达标重点难点过关沪教版第5课时《整十数除两、三位数的计算方法》一、单选题。

1.下面除法算式中,正确的是()A.570÷80=60 (90)B.440÷70=6 (20)C.140÷30=4 (2)D.33÷20=1 (3)2.三个车间组装电脑,第一车间40天组装了360台,第二车间30天组装了300台,第三车间50天组装了450台,()车间组装得快。

A.第一车间B.第二车间C.第三车间D.一样快3.用竖式计算268÷40,商应写在()位上。

A.个B.十C.百D.不确定4.439÷50是一道有余数的除法,余数应该()A.大于50B.等于50C.小于50D.小于4395.有水果店运来248千克的苹果,每40千克装一袋,可以装()袋,还剩下()千克。

A.8,6B.8,8C.6,6D.6,86.在下面的算式中,商是一位数的算式有()个。

148÷20, 372÷30, 186÷20,479÷40A.1B.2C.3D.47.48个小朋友排队坐船,每条船坐10人,至少要()条船。

A.2B.3C.4D.58.750÷70=10…()A.5000B.500C.50D.5二、判断题。

1.判断对错:230÷70=3……20()2.□67÷50,要使商是一位数,方框里最大填5。

()三、填空题。

1.口算540÷90时,想90×_____=540,所以540÷90=_____。

2.163里面最多有______个40;560是8的______倍。

3.490连续减去70,减______次得0;961连续减去_______个60后还余1。

4.在下面的□里填上合适的数。

□÷30=7...17□÷20=6...2□÷40=3 (36)5.小巧在计算一道除法算式时,错把被除数428看成了420,结果得出的商是7,正确的商是_____,余数是_____。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《计算方法》习题答案第一章 数值计算中的误差1.什么是计算方法?(狭义解释)答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。

2.一个实际问题利用计算机解决所采取的五个步骤是什么?答:一个实际问题当利用计算机来解决时,应采取以下五个步骤: 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果4.利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值,并编程获得解。

解:400)(2345-+⋅+-⋅+=x x x x x x P ,从而所以,多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。

5.叙述误差的种类及来源。

答:误差的种类及来源有如下四个方面:(1)模型误差:数学模型是对实际问题进行抽象,忽略一些次要因素简化得到的,它是原始问题的近似,即使数学模型能求出准确解,也与实际问题的真解不同,我们把数学模型与实际问题之间存在的误差称为模型误差。

(2)观测误差:在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的,由于仪器的精密性,实验手段的局限性,周围环境的变化以及人们的工作态度和能力等因素,而使数据必然带有误差,这种误差称为观测误差。

(3)截断误差:理论上的精确值往往要求用无限次的运算才能得到,而实际运算时只能用有限次运算的结果来近似,这样引起的误差称为截断误差(或方法误差)。

(4)舍入误差:在数值计算过程中还会用到一些无穷小数,而计算机受机器字长的限制,它所能表示的数据只能是一定的有限数位,需要把数据按四舍五入成一定位数的近似的有理数来代替。

这样引起的误差称为舍入误差。

6.掌握绝对误差(限)和相对误差(限)的定义公式。

答:设*x 是某个量的精确值,x 是其近似值,则称差x x e -=*为近似值x 的绝对误差(简称误差)。

若存在一个正数ε使ε≤-=x x e *,称这个数ε为近似值x 的绝对误差限(简称误差限或精度)。

把绝对误差e 与精确值*x 之比***xx x x e e r -==称为近似值x 的相对误差,称*x εη=为近似值x 的相对误差限η≤r e ,由于真值*x 是未知的,所以常常用xe x x x e r =-=*来表示相对误差,于是相对误差可以从绝对误差求出。

7.近似值的规格化表示形式如何?答:一般地,对于一个精确值*x ,其近似值x 的规格化形式为mp x x x x 10.021⨯±= ,其中{}),2,1(9,2,1,0,01p i x x i =∈≠,p 为正整数,m 为整数。

8.有效数字的概念是什么?掌握有效数字与误差的关系。

答:若近似值x 的(绝对)误差限是它的某一位的半个单位,也就是说该近似值准确到这一位,且从该位起直到前面第一个非零数字为止的所有数字都称为有效数字。

若近似值x 的(绝对)误差限为n m x x e -⨯≤-=1021*,则称x 为具有n 位有效数字的有效数,或称它精确到nm -10位,其中的每一位数字n x x x ,,21都是x 的有效数字。

设精确值*x 的近似值x 的规格化形式为m p x x x x 10.021⨯±= ,若x 具有n 位有效数字,则其相对误差限为n r x e -⨯≤111021;反之,若x 的相对误差限为n r x e -⨯+≤1110)1(21,则x 至少有n 位有效数字。

9.下列各数都是对真值进行四舍五入后获得的近似值,试分别写出它们的绝对误差限,相对误差限和有效数字的位数。

(1)024.01=x (2)4135.02=x (3)50.573=x (4)600004=x (5)55108⨯=x ;解:(1)0005.01*1≤-=x x e ;0021.0*≤=-=xex x x e r ;有三位有效数字。

(2)00005.02*2≤-=x x e ;000121.0*≤=-=xex x x e r ;有四位有效数字。

(3)005.03*3≤-=x x e ;000087.0*≤=-=xex x x e r ;有四位有效数字。

(4)5.04*4≤-=x x e ;0000084.0*≤=-=xex x x e r ;有五位有效数字。

(5)5.05*5≤-=x x e ;000000625.0*≤=-=xex x x e r ;有六位有效数字。

10.为了使19的相对误差≤0.1%,问至少应取几位有效数字?解:由19的首位数是4.设近似数*x 有n 位有效数字,由定理4.1可知,相对误差001.010421)(1*≤⨯⨯≤-n r x e ,解得097.3≥n ,即取4位有效数字,近似数的相对误差不超过0.1%。

11.已知33,3100,1150)(*2==-+==x x x x x P y ,计算)3100(*p y =及)33(P y =,并求x 和y 的相对误差。

解: 55555.51150)3100()3100()3100(2*-≈-+==p y 281150)33()33()33(2-=-+==P y 333.0)(*≈-=x x x e 0101.0)()(≈=xx e x e r 44444.22)(*≈-=y y y e 801587.0)()(≈=yy e y e r 12.写出误差估计的一般公式(以二元函数),(y x f z =为例)。

解:二元函数),(y x f z =的绝对误差: )(|)(|)(),(),(y e yfx e x f z e y x y x ⋅∂∂+⋅∂∂≈二元函数的相对误差: z y e y f z x e x f z z e z e y x y x r )(|)(|)()(),(),(⋅∂∂+⋅∂∂≈=)(|)(|),(),(y e yfz y x e x f z x r y x r y x ⋅∂∂⋅+⋅∂∂⋅=13.用电表测得一个电阻两端的电压和流过的电流范围分别为V V 2220±=,A I 1.010±=,求这个电阻的阻值R ,并估算其绝对误差和相对误差。

解:2)(≤V e ,1.0)(≤I e ,又2,1,IV I R I V R I V R -=∂∂=∂∂=。

所以: 42.01.01002202101)(|)(|)(|)(|)(),(),(),(),(=⨯+⨯=⋅∂∂+⋅∂∂≤⋅∂∂+⋅∂∂≈I e I R V e V R I e IRV e V R R e I V I V I V I V21099.1)()(-⨯≈=RR e R e r 。

14.若01.045.0,01.003.1*2*1±=±=x x ,计算22121x e x y +=的近似值,并估计)(y e 及其上界。

解:45.0221)03.1(e y +≈ )(21))(()21()21()(2*22*21*11*11*1*x x x x e e x x x x e x e x y y y e -++-=+-+=-= ),(,01.0211006.2)(21))((*2221*11*12*2x x e e e x x x x x x ∈⨯⨯+⨯=-++-≤-ξξ15.已测得某场地长为m l 110=,宽d 的值为m d 80=,已知m l l l e 2.0)(*≤-=,m d d d e 1.0)(*≤-=,试求面积ld s =的绝对误差限和相对误差限。

解:由ld s =,l ds d l s =∂∂=∂∂,,m l l l e 2.0)(*≤-=,m d d d e 1.0)(*≤-=。

可得:301.0802.0110)(|)(|)(|)(|)(),(),(),(),(=⨯+⨯=⋅∂∂+⋅∂∂≤⋅∂∂+⋅∂∂≈d e ds l e l s d e d s l e l s s e d l d l d l d l 3104.3)()(-⨯≈=ss e s e r 。

16.掌握二元函数的加、减、乘、除和开方运算的绝对误差和相对误差估计公式。

解:(1)加、减运算:由于()1/=∂+∂x y x ()()(),1/,1/,1/-=∂-∂=∂-∂=∂+∂y y x x y x y y x ,所以()()()()()()()()()()()()()()()()()()()()()|||/|||/|||,//,,//,y e y x y x e y x x y x e y e y x y x e y x x y x e y e x e y x e y e y x y x e y x x y x e y e x e y x e r r r r r r r r r ⋅-+⋅-≤-⨯--⨯-≈--≈-⨯++⨯+≈++≈+从而有(2)乘法运算: 由于()(),x yxy y x xy =∂∂=∂∂,所以()()()()()()y e x e xy e y xe x ye r r r +≈+≈,x y e ,从而()()()|||||||||y e x x e y xy e ⋅+⋅≤(3)除法运算: 由于2)(,1)(yx y y x y x yx-=∂∂=∂∂,所以)()(1)(2y e yxx e y y xe -≈,)()()(y e x e yxe r r r -≈(4)乘方及开方运算:由于()1-=∂∂n nnx xx ,所以()()()()x ne x e x e nx x e r n r n n ≈≈-,1 17.求方程01562=+-x x 的两个根,使它至少具有4位有效数字(982.27783≈)。

解:782.55982.272812114)56(5621=+≈⨯⨯⨯--+=x017863.0782.55112≈==x c x 19.求方程01162=+-x x 的较小正根,要求有3位有效数字。

解:937.15937.7812114)16(1621=+≈⨯⨯⨯--+=x062747.0937.15112≈==x c x 所以较小正根为062747.02≈x 。

20.设4110,,2,1,0, ==⎰n dx e xI x nn 。

(1)证明:4110,,2,1,0, =-=-n nI e I n n ;(2)给出一个数值稳定的算法,并证明算法的稳定性。

(1)证明:11111---=-===⎰⎰⎰n x n xn xnn nI e x d e nx e e d x dx e xI(2))(11n n I e nI -=- 设n n n I I e -=*,则n nnn n n n n n n e n I I e e nI I e e n I I e 1110*0022*221*11=-==-==-=------当n 无限大时,n e 越小,所以该算法稳定。

相关文档
最新文档