河北工业大学数值分析直接交.docVIP

《数值分析》 实验报告 学院：计算机科学与软件学院 姓名：Mr 班级：计算机999班 学号：77777 实验一：舍入误差与数值稳定性 实验目的： 通过上机编程，复习巩固以前所学程序设计语言； 通过上机计算，了解舍入误差所引起的数值不稳定性。 通过上机计算，了解运算次序对计算结果的影响，从而尽量避免大数吃小数的现象。 实验内容：用两种不同的顺序计算，分析其误差的变化。 实验流程图： 实验源程序: #include iostream using namespace std; double abs(double x) { if (x=0) return x; else return -x; } int main() { double X=1.644834; float sum1=0.0,sum2=0.0; for(float i=1.0;i10001;i++) sum1=sum1+1.0/(i*i); for(float j=10000.0;j0;j--) sum2=sum2+1.0/(j*j); coutsum1=sum1 E=((abs(sum1-X))/X)*100%endl; coutsum2=sum2 E=((abs(sum2-X))/X)*100%endl; return 0; } 实验结果： 实验分析：在运算中，如果参加运算的数的数量级相差很大，而计算机位数有限，如不注意运算次序就可能出现大数“吃掉”小数的现象，进而影响计算结果的可靠性，所以计算过程中要注意运算次序，避免出现这种现象。 实验二：拉格朗日插值法和牛顿插值法 实验目的：分别用拉格朗日差值和牛顿插值解决数学问题，并比较各方法的优劣。 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时的函数近似值。 实验流程图： 实验源程序: #include iostream using namespace std; int main() { int i; do{ float X,Y,X1,X2,X3,Y1,Y2,Y3,a,b,c; cout请输入(X1,Y1)：endl; cinX1Y1; cout请输入(X2,Y2)：endl; cinX2Y2; cout请输入(X3,Y3)：endl; cinX3Y3; cout请输入插值数X：endl; cinX; a=(X-X2)*(X-X3)/(X1-X2)/(X1-X3); b=(X-X1)*(X-X3)/(X2-X1)/(X2-X3); c=(X-X1)*(X-X2)/(X3-X1)/(X3-X2); Y=Y1*a+Y2*b+Y3*c; coutY=Yendl; coutcontinue?1 or 0endl; cini; }while(i==1); return 0; } 实验结果： 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时的函数近似值。 实验源程序: #includeiostream using namespace std; void main() { float x[11],y[11][11],xx,temp,a; int i,j,n,J; cout请输入插值的次数(n11) :n=; cinn; cout请输入n+1组值:endl; for(i=0;in+1;i++) { coutxi=; cinx[i]; coutyi=; ciny[0][i]; } for(i=1;in+1;i++) for(j=i;jn+1;j++) { if(i1) y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-i]); else y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-1]); } do{ temp=1;a=y[0][0]; printf(请输入要运算的值:x=); cinxx; for(i=1;in+1;i++) {