河北工业大学数分析Romo.docxVIP

河北工业大学 《数值分析》课程实验 实 验 报 告 题目： 数值分析实验报告 专业： 计算机 班级： 计666666 姓名： Ray 学号： 666666 完成日期： 2014.01.10 实验一 用两种不同的顺序计算，分析其误差的变化 开始算法流程 开始 输入sum的初始值以及n=1 输入sum的初始值以及n=1 n=10000 n=10000? 否是 否 是 利用公式求sum 利用公式求sum n=n+1 n=n+1 结束 结束 #includestdio.h int main() { float sum1=0.0,sum2=0.0,n=1.0; for(n=1.0;n=10000;n++) { sum1=sum1+1.0/(n*n); } printf(sum1=%f\n,sum1); for(n=10000;n=1;n--) { sum2=sum2+1.0/(n*n); } printf(sum2=%f\n,sum2); } 实验二 插值 拉格朗日插值 按下列数据 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时的函数近似值 开始 开始 分别输入各点坐标 分别输入各点坐标 定义变量L1，L2，L3,X 定义变量L1，L2，L3,X 带入数值计算L 带入数值计算L 结束 结束 源程序 #includestdio.h #includemath.h void main(void) { float x1,y1,x2,y2,x3,y3,x,y; float L,L1,L2,L3; printf(请输入x1,y1\n); scanf(%f%f,x1,y1); printf(请输入x2,y2\n); scanf(%f%f,x2,y2); printf(请输入x3，y3\n); scanf(%f%f,x3,y3); printf(请输入x\n); scanf(%f,x); { L1=(x-x2)*(x-x3)/(x1-x2)/(x1-x3); L2=(x-x1)*(x-x3)/(x2-x1)/(x2-x3); L3=(x-x1)*(x-x2)/(x3-x1)/(x3-x2); L=y1*L1+y2*L2+y3*L3; } printf(%f\n,L); } 牛顿插值 按下列数据 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时的函数近似值. 开始 开始 定义数组x[11],y[11][11] 定义数组x[11],y[11][11] 输入插入次数n和要运算的数 输入插入次数n和要运算的数 i=1 i=1 否in+1 否 in+1 是 是 i=i+1 i=i+1,j=1 否jn+1 否 jn+1 是是 是 是 计算y[i][j] 计算y[i][j] j=j+1 j=j+1 i=1 i=1 否in+1 否 in+1 是 是 a=a+y[ a=a+y[i][j] i=i+1 i=i+1 输出a 输出a 结束 结束 源程序 #includestdio.h void main() { float x[11],y[11][11],xx,temp,a; int i,j,n; printf(牛顿插值:\n请输入要运算的值:x=); scanf(%f,xx); printf(请输入插值的次数(n11):n=); scanf(%d,n); printf(请输入%d组值:\n,n+1); for(i=0;in+1;i++) { printf(x%d=,i); scanf(%f,x[i]); printf(y%d=,i); scanf(%f,y[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]); } temp=1;a=y[0][0]; for(i=1;in+1;i++) { temp=temp*(xx-x[i-1]); a=