数值分析报告作业.docx

数值计算方法作业 实验名称 实验4. 3三次样条插值函数 (P126) 4. 5三次样条插值函数的收敛性 (P127) 实验时间 姓名 班级 学号 成绩 实验4?3三次样条差值函数 实验目的： 掌握三次样条插值函数的三弯矩方法 实验函数： -dt X 0. 0 0. 1 0. 2 0. 3 0.4 F(x) 0. 5000 0. 5398 0. 5793 「0.6179 0. 7554 求f (0. 13)和f (0.36)的近似值 实验内容： 编程实现求三次样条插值函数的算法，分别考虑不同的边界条件； 计算各插值节点的弯矩值； 在同一坐标系中绘制函数f (x),插值多项式，三次样条插值多项式的曲线 比较 插值结果。 实验4. 5三次样条差值函数的收敛性 实验目的： 多项式插值不一定是收敛的，即插值的节点多，效果不一定好。对三次样条插值函数如何 呢？理论上证明三次样条插值函数的收敛性是比较困难的，通过本实验 可以证明这一理论 结果。 实验内容： 按照一定的规则分别选择等距或非等距的插值节点，并不断增加插值节点的个数。 实验要求： 随着节点个数的增加，比较被逼近函数和三样条插值函数的误差变化情 况，分析 所得结果并与拉格朗日插值多项式比较； 三次样条插值函数的思想最早产生于工业部门。作为工业应用的例子，考 虑如下例子：某汽车制造商根据三次样条插值函数设计车门曲线，其中一段数 据如下： Xk 0 1 2 3 4 5 6 7 8 9 10 yk 0. 0 0. 79 1. 53 2. 19 2. 71 3. 03 3. 27 2. 89 3. 06 3. 19 3.29 yk 0.8 0. 2 算法描述: 拉格朗日插值： 错误！未找到引用源。 n (x _ X )其中错误！未找到引用源。是拉格朗日基函数，其表达式为：h(x) 丄 牛顿插值: Nn(x) =f (Xg) f[Xo, Xi] (X - xO) f[Xo, Xi, X2) (X - xO) (x - Xi) .... f[Xg, Xi, X2,??? Xn] (X -x ) (X -Xi)... (X -Xn J ) f[Xi,Xj]二竺士 x】-Xj f[Xi,Xj)Xk]^Xj，Xk]——= Xk -Xf[Xi,xj] 其中 f [Xg, Xi. . . Xn] =(f[Xi,X2, ...Xn]-f[Xo, Xi, ...Xn J ) / (Xn -Xg) 三样条插值: 所谓三次样条插值多项式Sn(x)是一种分段函数，它在节点 Xi (aXOXl……Xnb)分成的每个小区间上是三次多项式，其在此区间 上的表 达式如下： S(x)二M, (Xi x)3 ? Mi (x Xy)s. [ y yi 4 hi (Mi My) 6 hi 6hx ( hi 6 (7i hi x i 4 ) ? (hi x i hi 6 ( 6 ,X [x一n] 式中 Mi二 S (Xx). 因此，只要确定了 Mi的值，就确定了整个表达式， Mi的计算方法如下: — hi 」hi十hy令“二 丄匚 hi +h“ d 纟 (yi \ yi 二 一 h】)〃叽必宀 % +?屮 hi屮 则Mi满足如下n-1个方程： 丄 i M 訂 2M 珀 M i 4 二山，i -1,2,?..n -1 常用的边界条件有如下几类： 给定区间两端点的斜率 mo, mn, BP S(xo) = y。二 m , S(Xn)二 yn 二 m* 给定区间两端点的二阶导数MO ,Mn,即S”(x。)二y。二MS”(xj二y■二M . 假设y二f(x)是以bp为周期的周期函数，则要求三次样条插值函数 S(x)也 为周期函数，对S (x)加上周期条件S(p, (x0 - 0)二S(p) (x-0), p二0,1,26 八—y° 2Mo +Mi m hi 对于第一类边界条件有 Mn J 2Mn 早(mn-仏丛) hn hn 对于第二类边界条件有2M0 +a 对于第二类边界条件有 2M0 +aM AnMn「2Mn=dn 0 do 其中 6 ■ 〒(认丿Iilll刀、y 6巴 dn =— (m - f[Xn4, XnD 2(1 - Un) M n hn y y。二 y n, M 0二 M n, S(x。0)二 S(Xn - 0),由此推得 *2. 乂 0 rM° 一 dO 1 円221.???? Mi di M2 d2 4 2 …? 11 4■■厶? n 4 Mn_l s n 4 ■丿 Mn 一 - dn 一 那么解就可以为 对于第三类边界条件, 2Mo oMl %Mn T — do，北中 ?Fv ?〃代』VWX°，Z] 一 3,那么解就可以为: 2 .. ,,, ?0 ? X 0 ? -do i 己?2?人????? M, d