- 1、本文档共23页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
0.1算法
1、 (p.11,题1)用二分法求方程在[1,2]内的近似根,要求误差不超过10-3.
【解】 由二分法的误差估计式,得到.两端取自然对数得,因此取,即至少需二分9次.求解过程见下表。
符号
0
1
2
1.5
+
1
2
3
4
5
6
7
8
9
2、(p.11,题2) 证明方程在区间[0,1]内有唯一个实根;使用二分法求这一实根,要求误差不超过。
【解】 由于,则在区间[0,1]上连续,且,,即,由连续函数的介值定理知,在区间[0,1]上至少有一个零点.
又,即在区间[0,1]上是单调的,故在区间[0,1]内有唯一实根.
由二分法的误差估计式,得到.两端取自然对数得,因此取,即至少需二分7次.求解过程见下表。
符号
0
0
1
0.5
1
2
3
4
5
6
7
0.2误差
1.(p.12,题8)已知e=2.71828…,试问其近似值,,x2=2.71,各有几位有效数字?并给出它们的相对误差限。
【解】有效数字:
因为,所以有两位有效数字;
因为,所以亦有两位有效数字;
因为,所以有四位有效数字;
;
;
。
评 (1)经四舍五入得到的近似数,其所有数字均为有效数字;
(2)近似数的所有数字并非都是有效数字.
?
2.(p.12,题9)设,,均为经过四舍五入得出的近似值,试指明它们的绝对误差(限)与相对误差(限)。
【解】 ,;
,;
,;
评 经四舍五入得到的近似数,其绝对误差限为其末位数字所在位的半个单位.
3.(p.12,题10)已知,,的绝对误差限均为,问它们各有几位有效数字?
【解】 由绝对误差限均为知有效数字应从小数点后两位算起,故,有三位;有一位;而,也是有一位。
1.1泰勒插值和拉格朗日插值
1、(p.54,习题1)求作在节点的5次泰勒插值多项式,并计算和估计插值误差,最后将有效数值与精确解进行比较。
【解】由,求得;;;;;,所以
插值误差:,若,则
,而,精度到小数点后5位,
故取,与精确值相比较,在插值误差的精度内完全吻合!
2、(p.55,题12)给定节点,试分别对下列函数导出拉格朗日余项:
(1);
(2)
【解】依题意,,拉格朗日余项公式为
(1) → ;
(2)因为,所以
3、(p.55,题13)依据下列数据表,试用线性插值和抛物线插值分别计算的近似值并估计误差。
0
1
2
0.32
0.34
0.36
0.314567
0.333487
0.352274
【解】依题意,,拉格朗日余项公式为
线性插值
因为在节点和之间,先估计误差
;须保留到小数点后4为,计算过程多余两位。
抛物线插值
插值误差:
抛物线插值公式为:
经四舍五入后得:,与精确值相比较,在插值误差范围内完全吻合!
1.3分段插值与样条函数
1、(p.56,习题33)设分段多项式
是以0,1,2为节点的三次样条函数,试确定系数b,c的值.
【解】依题意,要求S(x)在x=1节点
函数值连续: ,
即:
一阶导数连续: ,
即:
解方程组(1)和(2),得,即
由于,所以S(x) 在x=1节点的二阶导数亦连续。
2、 已知函数 的一组数据,和,(1)求其分段线性插值函数;
(2)计算的近似值,并根据余项表达式估计误差。
【解】(1)依题意,将x分为[0,1]和[1,2]两段,对应的插值函数为,利用拉格朗日线性插值公式,求得
;
(2),而 ,实际误差为:。
由,可知,则余项表达式
1.4 曲线拟合
1、(p.57,习题35)用最小二乘法解下列超定方程组:
【解】 构造残差平方和函数如下:
,
分别就Q对x和y求偏导数,并令其为零:
: ,
: ,
解方程组(1)和(2),得
2、(p.57,习题37)用最小二乘法求形如 的多项式,使之与下列数据相拟合。
【解】令,则为线性拟合,根据公式(p.39,公式43),取m=2,a1=0,N=5,求得
;
依据上式中的求和项,列出下表
xi
yi
Xi (=xi2)
Xi2(=xi4)
Xi yi (=xi2yi)
19
19
361
130321
6859
25
32.3
625
390625
20187.5
31
49
961
923521
47089
38
73.3
1444
2085136
105845.2
44
97.8
1936
3748096
189340.8
∑
157
271.4
5327
7277699
369321.5
将所求得的系数代入方程组(1)和(2),得
;
;
即:。
2.1 机械求积和插值求积
1、(p.94,习题3)确定下列求积公式中的待定参数,使其代数精度尽量高,并指明求积公式所具
文档评论(0)