原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
样条插值法绘制公路
东南大学《数学实验》报告
学号 姓名 袁骏杰 成绩
实验内容： 样条值法绘制公路利用样条插值法，根据已知坐标，绘制整条曲线。 0 30 50 70 80 90 120 148 170 180 Y/米 80 64 47 42 48 66 80 120 121 138 X/米 202 212 230 248 268 271 280 290 300 312 Y/米 160 182 200 208 212 210 200 196 188 186 X/米 320 340 360 372 382 390 416 430 478 440 Y/米 200 184 188 200 202 240 246 280 296 308 X/米 420 380 360 340 320 314 280 240 200 Y/米 334 328 334 346 356 360 392 390 400
对于上表给出的数据，估计公路的长度。
三 实验原理
估测公路函数满足三次样条插值条件。
公路在(478，296)处折返，因而整条曲线不是函数曲线，故将公路在折点处分段。
由于h(x)未知，根据表中数值估测h’(x)的值，其中
h’(x0+)= h’(0+) =-16/30 h’(x28-)= h’(478-) =16/48
h’(x28+)= h’(478+) =-12/30 h’(x38-)= h’(200-) =-10/40
估测公路长度时，以x轴的一米为间隔微分公路函数，求其总和。
四 主要实验代码及结果
for i=2:n1-1
lmd1(i)=h1(i)/(h1(i-1)+h1(i));
mu1(i)=1- lmd1(i);
d1(i)=6*((Y1(i+1)-Y1(i))/h1(i)-(Y1(i)-Y1(i-1))/h1(i-1))/(h1(i-1)+h1(i));
end % 计算hj,μj,λj,dj
for k=1:m1
for i=1:n1-1
if Z1(k)=X1(i)Z1(k)=X1(i+1)
S1(k)=M1(i)*(X1(i+1)-Z1(k))^3/(6*h1(i))+M1(i+1)*(Z1(k)-X1(i))^3/(6*h1(i))+(Y1(i)-M1(i)*h1(i)^2/6)*(X1(i+1)-Z1(k))/h1(i)+(Y1(i+1)-M1(i+1)*h1(i)^2/6)*(Z1(k)-X1(i))/h1(i);
break
end
end
end % 获得S1(x)各点值
for k=1:m2
for i=1:n2-1
if Z2(k)=X2(i)Z2(k)=X2(i+1)
S2(k)=M2(i)*(X2(i+1)-Z2(k))^3/(6*h2(i))+M2(i+1)*(Z2(k)-X2(i))^3/(6*h2(i))+(Y2(i)-M2(i)*h2(i)^2/6)*(X2(i+1)-Z2(k))/h2(i)+(Y2(i+1)-M2(i+1)*h2(i)^2/6)*(Z2(k)-X2(i))/h2(i);
break
end
end
end % 获得S2(x)各点值
% 估算公路长度
L=0;
for t=1:477
L=L+((Z1(t)-Z1(t+1))^2+(S1(t)-S1(t+1))^2)^0.5;
end
for t=1:277
L=L+((Z2(t)-Z2(t+1))^2+(S2(t)-S2(t+1))^2)^0.5;
end
L
运行结果
公路长度约为1016.3米
L =
1.0163e+003
本次实验采用样条插值的方法来拟合公路，从实验中熟悉了曲线拟合的方法，并对一些特殊情况作出了相应的处理，不局限于实验所提出的要求，着力于使求解出来的模型更加符合实际情况，而非仅仅完成一项实验。
3
文档评论（0）