原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析大作业-------第二次
校区：徐汇 班级：B1329191
学号：1132919012 姓名：杨小亮
第二次大作业
1．给定插值条件如下：
i 0 1 2 3 4 5 6 7
Xi 8.125 8.4 9.0 9.485 9.6 9.959 10.166 10.2
Yi 0.0774 0.099 0.280 0.60 0.708 1.200 1.800 2.177
作三次样条函数插值，取第一类边界条件
Y0’=0.01087 Y7’=100
画出插值曲线的图像。
2．逆时针旋转座标轴45o 保持（1．）中结点和边界条件的几何关系不变，再次作三次样条函数插值，画出插值曲线的图像。
3．比较（1）、（2）的结果，能得到什么结论？
由实验结果知，（1）中三次样条曲线出现了振荡，而（2）中曲线形状较平滑。比较1,2的结果，可知三次样条差值图像会随着坐标轴的变化而发生形状改变。
本实验用到的程序：
clc；
clear；
x=[8.125 8.4 9.0 9.485 9.6 9.959 10.166 10.2];
y=[0.0774 0.099 0.280 0.60 0.708 1.200 1.800 2.177];
Y=[0.01087 y 100];
xx = 8.125:0.01:10.2;
yy = spline(x,Y,xx);
plot(xx,yy,b-,x,y,bo,MarkerFaceColor,g)
xlabel(x);ylabel(y);
grid on
axis equal
axis([8 10.3 0 2.8])
%旋转后
a=-pi/4;
r=[cos(a) -sin(a);sin(a) cos(a)];
%插值点旋转
xyr=r*[x; y];
xr=xyr(1,:);
yr=xyr(2,:);
dd=r*[1 1;0.01087 100];
xl1=dd(2,1)/dd(1,1);
xl2=dd(2,2)/dd(1,2);
Yr=[xl1 yr xl2];
xxr=min(xr):0.01:max(xr);
yyr=spline(xr,Yr,xxr);
figure(2)
plot(xxr,yyr,r-,xr,yr,bo,MarkerFaceColor,g)
view([-45 90])
grid on
xlabel(x);ylabel(y);
axis equal
文档评论（0）