荆楚理工数学建模答案解析.doc

荆楚理工学院 课程设计成果 学院:____________ 班 级: 学生姓名: 学 号: 设计地点（单位）____ ___________ _设计题目:______________________ 完成日期： 201年 月 日 指导教师评语: ______________ ____________________ _____ _________________________________________________________________________________________________________________________________________________________________________________________________________ __________ _ 成绩(五级记分制):______ __________ 教师签名:__________ _______________ 插值在面积估值中的应用 学 院专 业年级班别 学 号 学生姓名指导教师 2014年月日摘要据插值是指一组散乱（又称非均匀）分布的数据采样点在整个区域上构造一个基本的函数的过程拉格朗日、分段线性、三次样条三种插值的方法MATLAB中的算法程序，以及梯形公式计算数值积分 三次样条插值； x、y数据位于某零件端面的轮廓线上，Y1和Y2分别对应轮廓的上下边界线上的纵坐标。假设需要得到x坐标每改变0.1时的y坐标，试完成加工所需数据，画出曲线，并求加工端面的面积。 x 0 3 5 7 9 11 12 13 14 15 Y1 0 1.7 2.4 2.7 3.0 3.1 2.9 2.5 2.0 1.8 Y2 0 1.3 1.8 2.0 2.1 2.0 1.8 1.2 1.0 1.8 示意图: 二、问题的分析 题目所给出是一些离散的数据，并求出图形的近似面积，但难以找到它的解析式，直接求出其他一些点上的函数值是非常困难的，所以我们选择插值法。 三、模型假设 假设给定的点是完全精确的。 假设计算机模拟过程不出现差错。 四、符号及变量说明 拉格朗日插值法的相关变量： 1.记x对应的各点为(i=1,2...10)。 2.记Y1对应的各点为(j=1,2...10)。 3.记Y2对应的各点为(j=1,2...10)。 4.记由x和Y1两组数确定的Lagrange插值基函数为(j=1,2...10)。 5.记由x和Y2两组数确定的Lagrange插值基函数为(j=1,2...10)。 6.记由x和Y1两组数确定的Lagrange插值多项式为 7.记由x和Y2两组数确定的Lagrange插值多项式为 五、模型的建立 拉格朗日插值法： 一般地，若已知y=f(x)在互不相同 n+1 个点处的函数值( 即该函数过这n个点），则可以考虑构造一个过这n+1个点的、次数不超过n的多项式,使其满足： 解决这个问题时先从构造插值基函数入手,这里的插值基函数是n次多项式,且满足条件 这表明除xk以外的所有节点都是的零点,故 由确定其中的系数c,结果有 利用插值基函数容易得出拉格朗日插值多项式 分段线性插值： 所谓分段线性插值就是通过插值点用折线段连接起来逼近f（x）。当给定n +1个点x0x1xn上的函数值，先选取两个节点，使，后在区间上作线性插值得 三次样条插值法： 三次样条插值法，基于分段进行插值的思想。在给定的区间 [ a ,b]上有 n +1个节点满足，和这些点上的函数值，为样条节点。在每个区间上构造三次多项式 显然，解此多项式共需要 4 n 个条件，由节点提供了n +1 个条件；用每个内点的关系建立条件，又得到 3 n -3个条件；再附加两个边界条件，即可以唯一确定样条函数了。此样条函数的存在性和唯一性可以用待定系数法进行证明，这里略。 实际计算时还需要引入边界条件才能完成计算。边界通常有自然边界（边界点的二阶导为0），夹持边界（边界点导数给定），非扭结边界（使两端点的三阶导与这两端点的邻近点的三阶导相等）。一般的计算方法书上都没有说明非扭结边界的定义，但数值计算软件如Matlab都把非扭结边界条件作为默认的边界条件。 六、模型求解 先用插值法