原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
插值法部习题
20.给定数据表如下：
0.25
0.30
0.39
0.45
0.53
0.5000
0.5477
0.6425
0.6708
0.7280
试求三次样条插值S（x），并满足条件：
解：（1）在编辑窗口输入：
x=[0.25,0.30,0.39,0.45,0.53];
y=[0.5000,0.5477,0.6245,0.6708,0.7280];
dx0=1.0000;dxn=0.6868;
s=csfit(x,y,dx0,dxn)
s =
1.8863 -1.0143 1.0000 0.5000
0.7952 -0.7314 0.9127 0.5477
0.6320 -0.5167 0.8004 0.6245
0.3151 -0.4029 0.7452 0.6708
1、已知函数在下列各点的值
0.2 0.4 0.6 0.8 1.0
0.98 0.92 0.81 0.64 0.38
试用4次牛顿插值多项式及三次样条函数（自然边界条件）对数据进行插值。用图给出{（，），}，及.
解：
（1）在编辑窗口输入：
x=[0.2,0.4,0.6,0.8,1.0];
y=[0.98,0.92,0.81,0.64,0.38];
N=newpoly(x,y)
N =
-0.5208 0.8333 -1.1042 0.1917 0.9800
由此可以得出牛顿插值多项式为：
=
S=ThrSample2(x,y,0,0,0)
S =
125*(46/25*t-69/125)*(t-3/5)^2+125*(567/500-81/50*t)*(t-2/5)^2+25*(-57/140*t+57/350)*(3/5-t)^2-25*(-81/200+27/40*t)*(t-2/5)^2/5)^2
（2）
i
0 1 11 10
0.2 0.28 1.08 1.0
0.98 0.9596 0.2403 0.38
在编辑窗口输入：
x=[0.2,0.28,0.4,0.6,0.8,1.0,1.08];
y=[0.98 0.9596 0.92 0.81 0.64 0.38 0.2403];
cs=spline(x,[0 y 0]);xx=linspace(0.2,1.08,101);
plot(x,y,o,xx,ppval(cs,xx),-);
x1=[0.2,0.28,0.4,0.6,0.8,1.0,1.08];
y1=[0.98 0.9596 0.92 0.81 0.64 0.38 0.2403];
polyfit(x1,y1,4)
ans =
-0.4126 0.5487 -0.8475 0.1017 0.9893
plot(x,y,o,xx,ppval(cs,xx),-,x1,y1,-.r),grid
3．下列数据点的插值
x
0
1
4
9
16
25
36
49
64
y
0
1
2
3
4
5
6
7
8
可以得到平方根函数的近似，在区间[0,64]上作图．
（1）用这9个点作8次多项式插值．
（2）用三次样条（自然边界条件）程序求．
解：
在编辑窗口输入：
x=[0 1 4 9 16 25 36 49 64];
y=0:1:8;y=sqrt(x);
x1=[0 1 4 9 16 25 36 49 64];y1=0:1:8;
m=polyfit(x,y,8)
Warning: Polynomial is badly conditioned. Remove repeated data points
or try centering and scaling as described in HELP POLYFIT.
In polyfit at 81
m =
Columns 1 through 4
-0.00000000032806 0.00000006712680 -0.00000542920942 0.00022297151886
Columns 5 through 8
-0.00498070758014 0.06042943489173 -0.38141003716579 1.32574370075005
Column 9
文档评论（0）