数值分析实验报告(四).docxVIP

河北工业大学 《数值分析》课程实验 实 验 报 告 专业： 计算机科学与技术 班级： 计1305班 姓名： 岳帅鹏 学号： 133007 完成日期： 2015.1.11 实验一 实验要求 用两种不同的顺序计算，分析其误差的变化。 实验目的 1、通过上机编程，复习巩固以前所学程序设计语言； 2、通过上机计算，了解舍入误差所引起的数值不稳定性。 3、通过上机计算，了解运算次序对计算结果的影响，从而尽量避免大数吃小数的现象。 算法描述 1）首先用float精度正序相加得到10000项的和；2）接着用float精度逆序相加得到10000项的和；3）然后用double?精度正序相加得到10000项的和；4）再用double?精度逆序相加得到10000项的和；5）根据结果，进行比较分析问题。 源程序 #includestdio.h #includemath.h void main() { int i; float s1=0,s2=0,d1,d2; for(i=1;i=10000;i++) s1=s1+1.0f/(i*i); for(i=10000;i=1;i--) s2=s2+1.0f/(i*i); d1=(float)(fabs(1.644834-s1)); d2=(float)(fabs(1.644834-s2)); printf(正向求和结果为%f\n，误差为%f\n\n,s1,d1); printf(反向求和结果为%f\n，误差为%f\n\n,s2,d2); if(d1d2) printf(正向求和误差小于负向求和误差\n); else if(d1=d2) printf(正向求和误差等于负向求和误差\n); else printf(正向求和误差大于负向求和误差\n); } 运行结果 对算法的理解与感悟  实验二 实验要求 1、拉格朗日插值 按下列数据 x -3.0 -1.0 1.0 2.0 3.0 y 1.0 1.5 2.0 2.0 1.0 作二次插值，并求x=-2，x=0，x=2.75时的函数近似值 2牛顿插值 按下列数据 x 0.30 0.42 0.50 0.58 0.66 0.72 y 1.04403 1.08462 1.11803 1.15603 1.19817 1.23223 作五次插值，并求x=0.46，x=0.55，x=0.60时的函数近似值. 实验目的 1、熟悉拉格朗日插值多项式和牛顿插值多项式，注意其不同特点； 算法描述 拉格朗日插值 输入x,y(i=0,1,2,?,n),令L(x)=0; （2）对=0,1,2,?,n计算 l(x)= L? L+l(x)y 牛顿插值多项式 输入n,x,y (i=0,1,2?,n); 对k=1,2,3?,n, i=1,2, ?,k计算各阶差商f(x,x ?,x); 计算函数值 N(x)=f(x)+f[x, x](x- x)+?+f[x, x,?,x](x- x)(x- x)?(x-x) 源程序 拉格朗日 #includeiostream #includeiomanip #includemath.h using namespace std; float jisuanmokuai(int,float,float[],float[]); void main() { float a,x[50],y[50],l; int m,n; cout题目：按下列数据endl; coutx:-3.0 -1.0 1.0 2.0 3.0endl; couty:1.0 1.5 2.0 2.0 1.5endl; cout作二次插值，求x=-2,x=0,x=2.75时的函数近似值.endl; cout输入插值次数：endl; cinn; cout插入计算次数:endl; cinm; for(int i=0;i=m;i++) { for(int j=0;j=n;j++) { int h99=i+1,b99=j+1; cout输入第h99个数的第b99个节点及函数值:endl; cinx[j]y[j]; } 牛顿 #inclu