东南大学数值分析第四章三次样条插值.docVIP

第四章 多项式插值与函数最佳逼近 ——曲线拟合之3次样条插值 *****（学号） *****（姓名） 上机题目要求见教材P195,37题。 一、算法原理 题目要求编写第一边界条件的3次样条插值函数的通用程序，同时根据汽车门曲线值点构造三次紧压样条曲线函数。其基本原理如下 定义 设有N+1个点，其中。如果存在N个三次多项式，系数为满足如下性质： 则成为三次样条函数。 现证明其存在： 由于是分段三次多项式，其二阶导数是在区间内是分段线性的。根据线性拉格朗日插值可以表示为： 用代入上式，得 将上式积分两次，会引入两个积分常数，可得到如下形式： 将代入上式，并利用可得两个方程： 求解，并将所得的结果带入方程得 求式（7）的导数，并化简得 由上述方程可得如下方程 重组上述方程，得三角线性方程组，表示为 该式具有严格对角优势。算出系数后可由如下公式计算的样条系数。 二、流程图 用已知的N+1个点构造三次紧压样条曲线的问题。其通用程序流程图1所示。 图1关于求解三阶样条曲线算法流程图 具体步骤如下： 计算，， 联系端点约束条件求解样条函数的系数。 三、计算代码 核心代码 for k=2:N-1 temp=A(k-1)/B(k-1); B(k)=B(k)-temp*C(k-1); U(k)=U(k