试验1拉格朗日插值法.doc

前言 一、上机实验目的 上机实验的目的是提高学生对算法的理解程度，并掌握用实用工具进行数值计算的方法，通过实践环节理解数值分析的应用和研究方法。 二、实验基本内容 本课程实验内容分为6个实验。学生可以在课内机时先完成指导书中给出的程序或验证算法，理解所学的知识，在此基础上再编写其他应用程序。指导书中的15个实验如下。 1．拉格朗日插值法。 2．最小二乘拟合。 3．数值积分。 4．范数计算和LU分解。 5．牛顿迭代法。 6．龙格库塔法。 三、实验任务与时间安排 本课程是一门实践性很强的课程，除了在课内安排的实验外，鼓励同学在课外用相关技术 进行编程练习。具体实践课时安排如下： 实验序号 周次 实验内容 拟实验时间 实验形式 所需教师数 所需实验仪器 备注 1 2，4 拉格郎日插值法 4 设计性 1 2 6 最小二乘拟合； 2 验证 1 3 8 数值积分； 2 验证 1 4 10 范数计算 2 验证 1 5 12 LU分解 2 验证 1 6 14 牛顿迭代法 2 验证 1 7 16 龙格库塔法 2 验证 1 实验指导书 实验1 拉格朗日插值法 一、实验目的[1] 了解的基本原理和方法； [2] 掌握用MATLAB的方法；[3] 编程实现lagrange插值 ?二、实验内容 x0=[0.4 0.5 0.6 0.7 0.8];%输入节点的X值 y0=[-0.916291 -0.693147 -0.510826 -0.356675 -0.223144];%输入节点的Y值 T=interp1(x0,y0,0.54,linear);　’cubic’,’ nearest’,’ spline’等。 x=0.2:0.05:1.0;%从0.2开始到1.0每隔0.05取一个待求点 y=interp1(x0,y0,x,linear);%用插值法求出待求点对应的函数值y plot(x0,y0,’r*’);%用红色星号标出插值点对 hold on;plot(x,y);%蓝色线为所得的插值函数图形 2． 三、实验步骤 1．开启软件平台——MATLAB； 2．根据各种数值解法步骤编写M文件 ．观察运行结果(数值或图形)； 四、实验要求与任务 一、实验目的[1] 了解最小二乘拟合的基本原理和方法； [2] 掌握用MATLAB作最小二乘多项式拟合和曲线拟合的方法； [3] 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。 [4] 了解各种参数辨识的原理和方法； [5] 通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程； ?二、实验内容 1．用MATLAB中的函数作一元函数的多项式拟合与曲线拟合； 2．用MATLAB中的函数作二元函数的最小二乘拟合，作出误差图； 3．针对预测和确定参数的实际问题，建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 ．观察运行结果(数值或图形)； 四、实验要求与任务 根据实验内容和步骤，完成以下具体实验 旧车价格预测某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？ xi 1 2 3 4 5 6 7 8 9 10 yi 2615 1943 1494 1087 765 538 484 290 226 204 a,s]=polyfit(xdata,ydata,n) 其中n表示多项式的最高阶数，xdata，ydata为将要拟合的数据，它是用数组的方式输入．输出参数a为拟合多项式的系数,s为误差。 另： 多项式在x处的值y可用下面程序计算．y=polyval(a,x)， 对上面给出的数据做多项式拟合，可取不同的n观测此时的误差，看取什么样的n较好。 并在matlab中画出拟合曲线(取拟合曲线上若干个点对，用plot命令实现) 实验3 数值积分 一、实验目的[1] 了解的基本原理和方法； [2] 掌握用MATLAB的方法； [3] 通过实例学习如何用。 [4] 了解； ?二、实验内容 1．用MATLAB中的函数函数的； ．用MATLAB中的函数。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 ．观察运行结果(数值或图形)； 四、实验要求与任务 根据实验内容和步骤，完成以下具体实验 实验提示： Quad函数的使用方法，参数意义 Q = QUAD(FUN,A,B) Fun为待求积的函数 A,B为积分区间 Quad函数是使用误差为10