数值分析实验报告(同名8562).docx

《数值分析实验报告(同名8562).docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、河北工业大学《数值分析》课程实验实验报告专业：计算机科学与技术班级：计1305班姓名：岳帅鹏学号：133007完成日期：2015.1.11实验一一、实验要求用两种不同的顺序计算，分析其误差的变化。二、实验目的1、通过上机编程，复习巩固以前所学程序设计语言；2、通过上机计算，了解舍入误差所引起的数值不稳定性。3、通过上机计算，了解运算次序对计算结果的影响，从而尽量避免大数吃小数的现象。三、算法描述1）首先用float精度正序相加得到10000项的和；2）接着用float精度逆序相加得到10000项的和；3）然后用double 精度正

2、序相加得到10000项的和；4）再用double 精度逆序相加得到10000项的和；5）根据结果，进行比较分析问题。四、源程序#include#includevoidmain(){inti;floats1=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))

3、;printf("正向求和结果为%f，误差为%f",s1,d1);printf("反向求和结果为%f，误差为%f",s2,d2);if(d1

4、-2，x=0，x=2.75时的函数近似值2牛顿插值按下列数据x0.300.420.500.580.660.72y1.044031.084621.118031.156031.198171.23223作五次插值，并求x=0.46，x=0.55，x=0.60时的函数近似值.二、实验目的1、熟悉拉格朗日插值多项式和牛顿插值多项式，注意其不同特点；三、算法描述1、拉格朗日插值（1）输入x,y(i=0,1,2,¼,n),令L(x)=0;（2）对=0,1,2,¼,n计算l(x)=L¬L+l(x)y2、牛顿插值多项式（1）输入n,x,y(i=0,

5、1,2¼,n);（2）对k=1,2,3¼,n,i=1,2,¼,k计算各阶差商f(x,x¼,x);（3）计算函数值N(x)=f(x)+f[x,x](x-x)+¼+f[x,x,¼,x](x-x)(x-x)¼(x-x)一、源程序1.拉格朗日#include#include#includeusingnamespacestd;floatjisuanmokuai(int,float,float[],float[]);voidmain(){floata,x[50],y[50],l;intm,

6、n;cout<<"题目：按下列数据"<>n;cout<<"插入计算次数:"<>m;for(inti=0;i<=m;i++){for(intj=0;j<=n;j++){inth99=i+1,b99=j+1;cout<<"输入第"<<

7、h99<<"个数的第"<>x[j]>>y[j];}2.牛顿#include#defineM6floatw(floatX,intn,floatx[]);floatquotient(intk,inti,floatf[][M],floatx[],floaty[]);floatnewton(floatX,intn,floatf[][M],floatx[],floaty[]);voidmain(){floatx[M],y[M],f[M][M];x[0]=0.30f;y[

8、0]=1.04403f;x[1]=0.42f;y[1]=1.08462f;x[2]=0.50f;y[2]=1.11803f;x[3]=0.58f;y[3]=1.15603f;x[4]=0.66f;y[4]=1.19817f;x[5]=0.72f;y