原创力文档- 1、本文档共13页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
南京工业大学工程数学第5章作业解析
学号：642081701046 姓名：苗洁 班级：化工2班
第五章作业
1.解：x=[0.520,3.1,8.0,17.95,28.65,39.62,50.65,78,104.6,156.6,208.6,260.7,312.5,364.4,416.3,468,494,507,520];
y=[5.288,9.4,13.84,20.20,24.90,28.44,31.10,35,36.9,36.6,34.6,31.0,26.34,20.9,14.8,7.8,3.7,1.5,0.2];
pp2=csape(x,y,complete,[1.86548,-0.046115])
pp2 =
form: pp
breaks: [0.5200 3.1000 8 17.9500 28.6500 39.6200 50.6500 78 104.6000 156.6000 208.6000 260.7000 312.5000 364.4000 416.3000 468 494 507 520]
coefs: [18x4 double]
pieces: 18
order: 4
dim: 1
pp2.coefs
ans =
-0.0029 -0.0979 1.8655 5.2880
0.0080 -0.1201 1.3031 9.4000
-0.0003 -0.0029 0.7007 13.8400
0.0003 -0.0128 0.5447 20.2000
0.0000 -0.0039 0.3655 24.9000
0.0000 -0.0038 0.2803 28.4400
0.0000 -0.0030 0.2052 31.1000
-0.0000 -0.0009 0.0990 35.0000
0.0000 -0.0013 0.0403 36.9000
-0.0000 -0.0000 -0.0301 36.6000
0.0000 -0.0004 -0.0527 34.6000
0.0000 -0.0002 -0.0816 31.0000
0.0000 -0.0001 -0.0979 26.3400
-0.0000 -0.0001 -0.1109 20.9000
-0.0000 -0.0002 -0.1254 14.8000
-0.0000 -0.0002 -0.1462 7.8000
0.0001 -0.0009 -0.1748 3.7000
0.0000 0.0030 -0.1467 1.5000
fnplt(pp2)
再求x=2,30,130,350,515各点上的函数值，则输入
xi=[2,30,130,350,515];
yi=ppval(pp2,xi)
yi =
7.8252 25.3862 37.2138 22.4751 0.5427
三次样条插值函数为
=++[-]+[-]
由以上结果可得：
当[0.520，3.1]时：
当[28.65，39.62]时：
当[104.6，156.6]时：
当[312.5，364.4]时：
当[507，520]时：
综上，计算结果为：
2 30 130 350 515 7.8252 25.3862 37.2138 22.4751 0.5427 4.1121 0.3296 -0.1053 -0.0991 -0.0574
2.解：
（1）复化公式的余项估计,
==≤≤
35
采用复化的要求具有9位有效数字则至少要取个点
利用MATLAB程序计算
先求的四阶导数： syms x
f=4/(1+x^2) %定义函数f（x）
n=input(输入所求导数阶数:)
f2=diff(f,x,n)
f =
4/(x^2 + 1)
输入所求导数阶数:4
n =
4
f2 =
96/(x^2+1)^3-(1152*x^2)/(x^2+1)^4+(1536*x^4)/(x^2+1)^5
复化辛普森程序：
syms x
f=inline((4/(1+x^2)),x);
f2=inline((96/(x^2 + 1)^3 - (1152*x^2)/(x^2 + 1)^4 + (1536*x^4)/(x^2 + 1)^5),x);
文档评论（0）