- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Matlab程序三次样条插值函数
已知一组数据点,编写一程序求解三次样条插值函数满足
并针对下面一组具体实验数据
0.25
0.3
0.39
0.45
0.53
0.5000
0.5477
0.6245
0.6708
0.7280
求解,其中边界条件为.
解:Matlab计算程序为:
clear
clc
x=[0.25 0.3 0.39 0.45 0.53]
y=[0.5000 0.5477 0.6245 0.6708 0.7280]
n=length(x);
for i=1:n-1
h(i)=x(i+1)-x(i);
end
for i=1:n-2
k(i)=h(i+1)/(h(i)+h(i+1));
u(i)=h(i)/(h(i)+h(i+1));
end
for i=1:n-2
gl(i)=3*(u(i)*(y(i+2)-y(i+1))/h(i+1)+k(i)*(y(i+1)-y(i))/h(i));
end
g0=3*(y(2)-y(1))/h(1);
g00=3*(y(n)-y(n-1))/h(n-1);
g=[g0 gl g00];
g=transpose(g)
k1=[k 1];
u1=[1 u];
Q=2*eye(5)+diag(u1,1)+diag(k1,-1)
m=transpose(Q\g)
syms X;
for i=1:n-1
p1(i)=(1+2*(X-x(i))/h(i))*((X-x(i+1))/h(i))^2*y(i);
p2(i)=(1-2*(X-x(i+1))/h(i))*((X-x(i))/h(i))^2*y(i+1);
p3(i)=(X-x(i))*((X-x(i+1))/h(i))^2*m(i);
p4(i)=(X-x(i+1))*((X-x(i))/h(i))^2*m(i+1);
p(i)=p1(i)+p2(i)+p3(i)+p4(i);
p(i)=simple(p(i));
end
s1=p(1)
s2=p(2)
s3=p(3)
s4=p(4)
for k=1:4
for z=x(k):0.001:x(k+1)
q=eval(subs(p(k),X,z));
plot(z,q,b)
hold on
end
end
grid on
legend(Èý´ÎÑùÌõÇúÏß)
title(²åÖµ)
xlabel(x)
ylabel(p)
计算结果为:
s1 =
-(705394867539368680*X^3 - 529046150654530286*X^2 + 23087199381998953*X - 40023205577025431)/112589990684262400
s2 =
(257361898089296225*X^3)/136796838681378816 - (160081743506404901*X^2)/60798594969501696 + (404209705972252727*X)/202661983231672320 + 429142243010323951/3166593487994880000
s3 =
- (3495912536773825*X^3)/7599824371187712 + (1437374409830143*X^2)/13510798882111488 + (10427488839800859*X)/11258999068426240 + 12358231431982259943/45035996273704960000
s4 =
(38626753769033575*X^3)/18014398509481984 - (245666153971053021*X^2)/72057594037927936 + (3614707928905781673*X)/1441151880758558720 + 26732501105874704913/720575940379279360000
通过绘制曲线,发现为3次样条曲线,且数据拟合较好。
您可能关注的文档
- HND计算机大二·。.doc
- HSE脚手架 HSE SCAFFOLDING.doc
- HSE管理体系建设培训.ppt
- HTC One V 简体中文 word版 说明书.doc
- HSE管理体系建设模式.ppt
- htV甲醛检测仪使用说明.doc
- HTML资料大集合.doc
- hypermesh入门实例、网格划分注意事项、主要面板英文翻译.doc
- ICU综合征的影响因素及护理进展.doc
- ICU侵袭性真菌感染早期经验性治疗.ppt
- 2023年江苏省镇江市润州区中考生物二模试卷+答案解析.pdf
- 2023年江苏省徐州市邳州市运河中学中考生物二模试卷+答案解析.pdf
- 2023年江苏省苏州市吴中区中考冲刺数学模拟预测卷+答案解析.pdf
- 2023年江苏省南通市崇川区田家炳中学中考数学四模试卷+答案解析.pdf
- 2023年江西省吉安市中考物理模拟试卷(一)+答案解析.pdf
- 2023年江苏省泰州市海陵区九年级(下)中考三模数学试卷+答案解析.pdf
- 2023年江苏省苏州市高新二中中考数学二模试卷+答案解析.pdf
- 2023年江苏省南通市九年级数学中考复习模拟卷+答案解析.pdf
- 2023年江苏省南通市海安市九年级数学模拟卷+答案解析.pdf
- 2023年江苏省泰州市靖江外国语学校中考数学一调试卷+答案解析.pdf
文档评论(0)