插值算法实现.pdf

计算方法——插值算法的MATLAB 实现 专业：电气工程及其自动化 计算方法实验报告 姓名：____齐琛森_ ______ 学号：____3120104723___ 日期：_2014-4-22 ___ 实验名称： 插值算法的MATLAB 实现 __ 1 引言 在自然科学和工程中有很多问题的解决归结为利用实验中得到的一些数据来作为根据求 解几个变量之间的函数关系。在这种情况下，利用这些数据点进行分段分析来得到一个分段 函数关系式就会至关重要。由分段线性与非线性的插值法就是其中最为主要的代表，也是计 算机进行拟合的重要选择。其中，线性插值是利用两点确定一条直线的原理，利用两点之间 一个直线来实现差值，得到拟合曲线。而二次差值和拉格朗日差值，则是通过几点之间高次 差值，在解一个高次方程组的原理下，得到一个基本的拟合曲线。。本实验所用的差值方式就 是分段线性差值与分段二次差值和全区间上的拉格朗日插值。 2 实验目的和要求 （一）实验目的 （1） 了解基本差值的过程和基本算法原理。 （2 ） 对现有数据进行一下简单的整理 （3 ） 根据自己对插值手法的理解设计插值计算的流程图 （4 ） 根据流程图设计插值程序 （5 ） 将课本所给的数据算例带入程序中进行运行 （6 ） 将结果与原本的数据点进行比较，分析结果和。 （7 ） 更换插值方式数据进行插值分析，观察不同的之间有什么样的区别和联系函数图 象。 （8 ） 感受差值方法对于解决数字分析的有什么样的优势 （二）实验要求 （1）了解插值算法的要求，对数据的要求和工作条件 （2 ）会用以及了解matlab 的使用方法和数据处理 （3 ）总结分析实验结果 3 算法原理与流程图 （1） 线性插值算法原理： 2 n 已知n 次插值多项式P (x)=a +a x+a x +……+a x n 0 1 2 n 但是对于一个实验中得到的结果来说，但并不是每一次都会使用高次插值进行分 析，因为， f (n1) () R (x) f (x)  p (x)  (x) n n n1 (n 1)! 计算方法——插值算法的MATLAB 实现 有时候截断误差并不一定会随着n 的增大而减小，这时候往往会使用低次插值， 比如一次差值和二次差值。如果我们取到Pn(x)=a0+a1x 就是一个线性插值（一次 差值）的结果。 我们线性插值就是使用待求变量附近的两个数据点，通过两点之间得到一条直线 的原理，得到待求变量所在的直线，再根据直线方程得到待求变量的函数值，最 终以函数值为返回值，绘制函数图象。 流程图 输入x ，y i i 输入x P=xj Yes