张飞 数值分析实验报告.docVIP

江西理工大学研究生院 《数值分析》实验报告 姓 名： 张飞 专 业： 机械工程 学 号： 6720150104 日 期： 2015年12月12日 目 录 实验一 函数插值方法 3 实验二 函数逼近与曲线拟合 7 实验四 线方程组的直接解法 17 实验五 解线性方程组的迭代法 24 实验六 非线性方程求根 26 实验七 矩阵特征值问题计算 28 实验八 常微分方程初值问题数值解法 32  实验一 函数插值方法 一、问题提出 对于给定的一元函数的n+1个节点值。试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。 数据如下： （1） 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式，和分段三次插值多项式，计算, 的值。（提示：结果为, ） （2） 1 2 3 4 5 6 7 0.368 0.135 0.050 0.018 0.007 0.002 0.001 试构造Lagrange多项式，计算的,值。（提示：结果为, ） 二、问题分析 1、 利用Lagrange插值公式 编写出插值多项式程序； 2、 给出插值多项式或分段三次插值多项式的表达式； 3、 根据节点选取原则，对问题（2）用三点插值或二点插值，其结果如何； 4、 对此插值问题用Newton插值多项式其结果如何。Newton插值多项式如下： 其中： 三、实验程序及注释 1.（1）程序一 function f=Lagrange(x,fx,inx) x=[0.4 0.55 0.65 0.8 0.95 1.05] fx=[0.41075 0.57815 0.69675 0.90 1.0 1.25382] inx=[0.596,0.99]; n=length(x); m=length(inx); for i=1:m; z=inx(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x(j))/(x(k)-x(j)); end end s=p*fx(k)+s; end f(i)=s; end plot(x,fx,O,inx,f) （2）运行结果： x = 0.4000 0.5500 0.6500 0.8000 0.9500 1.0500 fx =0.4108 0.5782 0.6967 0.9000 1.0000 1.2538 ans =0.6257 1.0542 2、（1）程序二 function f=Lagrange(x,fx,inx) x=[1 2 3 4 5 6 7] fx=[0.368 0.135 0.050 0.018 0.007 0.002 0.001] inx=[1.8 6.15]; n=length(x); m=length(inx); for i=1:m; z=inx(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x(j))/(x(k)-x(j)); （2）运行结果： x = 1 2 3 4 5 6 7 fx = 0.3680 0.1350 0.0500 0.0180 0.0070 0.0020 0.0010 ans= 0.1648 0.0013 四、实验数据结果及分析 1 . 五次Lagrange多项式的运行结果为 经过迭代达到了给定结果的精度实验图像如图像一 图像一 六次Lagrange多项式的运行结果为 经过迭代达到了给定结果的精度实验图像如图像二: 图像二 五、实验结论 一、问题提出 从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。 在某冶炼过程中，