昆明理工大学数值分析上机报告6.docxVIP

第 第 PAGE 1 页 昆 明 理 工 大 学 理 06 工科硕士 《数值分析》上机实验报告 专业： 材 料 学 姓 名 ： 学 号 ： 2006202044 任课教师： 作业完成实验室： 个人 PC 实验成绩： 理论描述 数值公式 程序流程图 数据和 讨论 源程序 总分 (10) (10) 和程序结构 结果 (20) (20) (100) (20) (20) 实验内容： 题目/要求： 三次样条插值法 一、问题提出 设已知数据如下： 0.2 0.4 0.6 0.8 1.0 0.9798652 0.9177710 0.8080348 0.6386092 0.3843735 1、满足自然边界条件S ??0.2?? S ??1.0 1、满足自然边界条件S ??0.2?? S ??1.0?? 0 ； 2、满足第一类边界条件S (0.2) ? 0.20271 ， S (1.0) ? 1.55741 。 二、要求 3 、 打 印 输 出 用 追 赶 法 解 出 的 弯 矩 向 量 (M , M ,..., M ) S (0.2 ? 0.1i) (i ? 0,1,...,8) 的值。﹡并画出 y ? S (x) 的图形﹡。 0 1 4 和 作业环境(包括选用的程序语言、运行环境) 程序语言 Turbo.C 2.0 运行环境 WINDOWS XP 3．数学（理论背景）描述 某些实际问题，如船体放样与机翼设计，要求插值曲线不仅连续而且处处平滑。甚至要求尽可能采用流线型，使气流沿机翼的表面能形成平滑的流线，以减少空气的阻力。换句话说，所谓光滑插值就是既要分段低次又要保证接头光滑。 为适应这类需求，借助曲线板来作图，以保证曲线在接头处的光滑。在工程技术如船体放样中，作图员常用“样条”在指定节点间做光滑曲 线。光滑插值就是这类作图方法的数学模拟。因此，这类插值称作样条 （Spline）插值。 样条插值实际上是一种改进的分段插值，它要求插值函数在各分段的衔接处能保持一定程度的光滑性。既要保持插值函数的导数的连续性。在实际应用中，最常用的是三次样条插值，其定义如下： 若函数S ?x?? C 2 ?a, b?，且在每个小区间?x j , x j ?1 ?上是三次多项式，其 中a ? x 0 ? x ? ? x ?1 n ? ? b 是给定节点，则称 S ?x?是节点 x 0 , x , , x 上的 ?1 n ? 三次样条函数。若在节点 x  上给定函数值 y j j ? f ?x j ??j ? 0,1,?, n?，并成 立S ?x j ?? y j ?j ? 0,1,?, n?，则称 S ?x?为三次样条插值函数。 数值计算公式 iii三次样条函数分段插值公式 i i i ? ? ?x ? x?3 ?x ? x ?3 ? h 2 ? x ? x ? h 2 ? x ? x ? ? S x ? i M ? 6h i?1 i i?1 6h i M ? ? y ? i ? i?1 6 M i?1 ? h ?i ? ? ? y ? ? i 6 M i ? h ?i ? i?1 1 ? ? ?x ? x?2 ?x ? x ? h ? ? y ? y ? ? S ? x ? ? i M ???i?1 M ? i M ? M ? i i?1 2 2h i?1 2h i i i 6 i i?1 h i ? ? x x x ? x ? ? S ? x ? i h i M i?1 ???i?1 M h i i 3 ，其中 ?2M ?? 0 ? ? M ? ? ? 0 1 0 1 ?? ? 1 M ? 2M 0 1 ? ? M ? ? 1 2 1 ? ? ??????? 4 ??1 ?? ?M ? 2M ? ? M ? ? ? n?1 n?2 n?1 n?1 n n?1 ??? M ? n  n?1 h ? 2M ? ? n n ? ? ? ???i?1 5 h ? x ? x i?i?1y i? i?1 y i?1 y ? ? i h i?1 ? y ? y i ? i?1 h i ? ? 6 i  h ? h i i?1  ?6? j j ?1 j 算法程序流程图 开始 开始 ? 读入数据 读入数据 ? 利用 利用（5）（6）计算系数 ? ,? , i ? 1,2,?, n ? 1 1 1 ? 用追赶法解方程组（ 用追赶法解方程组（4）得 M , M ,?, M 0 1 n ? 判定插值点 判定插值点 x 所在的分段 ?x , x ? i ?1 i ? 按公式 按公式（1）（2）和（3）计算所求函数值 S ?x?及其一阶二阶数 值导数