工程数值方法2第二章.pptVIP

例题 解方程得 于是 例题 例题 §2.4 三弯矩插值法 三次样条函数S(x)在每个子区间[xj-1i,xj] (j=1,2,???,n)上是三次多项式，在[a,b]上二阶连续可微，故 三弯矩插值法 三弯矩插值法 由上式可解出： 代入,得 三弯矩插值法 三弯矩插值法 因此，只要能求出所有的Mj，就能求出样 条插值函数S(x).下面考虑Mj的求法 由S(x)在节点的一阶导数的连续性 三弯矩插值法 从而导出如下方程 其中 问题1 从而导出如下方程 问题2 直接有 问题3 从而导出如下方程 其中 三弯矩插值法 在周期插值情形，由于要求 ，因此出现在方程组中的yn应理解为 问题1，2，3得到的都是三对角带状非奇异方程组。由于Mi相应于材料力学中梁的弯矩，而每一个方程式中只出现相邻的三个Mi，故上述方程称为三弯矩法。 求解步骤 总结以上论述，可得求三次样条的步骤为： （1）确定边界条件，判定是问题类型； （2）根据所确定的条件计算各值，形成方程组； （3）解三对角方程组，求得M0, M1 , M2, ??? Mn ； （4）将求得的Mi值代回S(x)的表达式中，从而可求 得函数y=f(x)在任一点的近似值S(x)。 例题 例3．已知f(x)在若干点处的值为f(0)=0, f(1)=1, f(2)=1, f(3)=0,试