数值逼近课程设计（2）——D1样条插值

最新推荐文章于 2023-07-01 08:00:00 发布

胜天半子·浩

最新推荐文章于 2023-07-01 08:00:00 发布

阅读量899

点赞数 1

文章标签： matlab MATLAB 三次样条插值法三次样条插值，三弯矩

本文链接：https://blog.csdn.net/m0_46330603/article/details/107070156

版权

本文介绍了数值分析中的D1样条插值方法，通过MATLAB编程求解，有效避免了高阶多项式插值的龙格现象。在4到6之间的插值结果为1.609770287689208，插值函数保持了一阶和二阶导数的连续性。

摘要由CSDN通过智能技术生成

问题背景

在这里插入图片描述
在数值分析这个数学分支中，样条插值是使用一种名为样条的特殊分段多项式进行插值的形式。由于样条插值可以使用低阶多项式样条实现较小的插值误差，这样就避免了使用高阶多项式所出现的龙格现象，所以样条插值得到了流行。
该问题的方法很好地解决了第一问的龙格现象，得到的插值函数有着很好的性质。
样条的概念来自船体放样中中所用的木头长条。一条船的外曲线要求很光顺，又要通过一些点，这些点称为型值点，工人技术人员就用这种样条来画曲线。两个相邻型值点之间的曲线，近似地为一个三次多项式。整个曲线通过各型值点处，其一阶导数、二阶导数是连续的。

数值演算

由matlab编写程序，利用D1样条插值函数在４到６之间的多项式表达式求值，得到的结果为：1.609770287689208
用D1样条插值得到的曲线为：
在这里插入图片描述

源代码

clear;clc;
format long
x=[1 2 3 4 6];
f=[log(1),log(2),log(3),log(4),log(6)];
f11=1;
f55=1/6;
f12=(f(2)-f(1))/(x(2)-x(1));
f23=(f(3)-f(2))/(x(3)-x(2));
f34=(f(4)-f(3))/(x(4)-x(3));
f45=(f(5)-f(4))/(x(5)-x(4));
f123=(f23-f12)/