算法报告实例.docVIP

《数值分析》实验报告 成绩:班级：信息12?2 成绩: 小组成员： 方 MM (121246402XX) 卜 MM (121246402XX) 谭 GG(121246402XX) 2014年12月25日 广东石油化工学院理学院信息系 目录 一.课题：线性方程组的矩阵分解法 TOC \o 1-5 \h \z HYPERLINK \l bookmark18 \o Current Document \h 运行代码 3 HYPERLINK \l bookmark8 \o Current Document \h 实验题目 4 HYPERLINK \l bookmark10 \o Current Document \h 运行结果 5 二、课题：线性方程组的迭代法 7 HYPERLINK \l bookmark6 \o Current Document \h 1.运行代码 7 HYPERLINK \l bookmark20 \o Current Document \h 2 ?实验题目 10 103?运行结果. 10 三、课题:拉格朗日插值法 三、课题: 拉格朗日插值法 15 运行代码 15 实验题目 16 17运行结果. 17 心得 20 、课题：线性方程组的矩阵分解法 1、运行代码: #include stdio.h #include math.h #define MAX_N 20 int main(int argc, char* argv[]) { int n; //未知数个数 int i,j,k; static double a[MAX_N][MAX_N], b[MAX_N], x[MAX_N], y[MAX_N]; static double l[MAX_N][MAX_N], u[MAX_N][MAX_N]; printf(\nlnput n value(dim of Ax=b):); scanf(%dH, n); if(n MAX_N) { printf(The input n is larger than MAX_N, please redefine the MAX_N.\rT); return 1; } if(n = 0) { printf(Please input a number between 1 and %d.\n, MAX_N); return 1; } // {{程序输入 printf(Now input the matrix a(i, j), i, j = 0,%d:\n, n-1); for (i=0; in; i++) for(j=0;jn;j++) scanf(H%lf\ a[i]|j]); printf(Now input the matrix b(i), i = 0,%d:\n, n-1); for(i=0; ivn; i++) scanf(n%ir, b[i]); // }}程序输入 for(i=0; ivn; i++) u[i][i] = l；// for(k=0; kn; k++) { for(i=k; in; i++) //计算L的第k列元素 l[i][k] = a[i][k]; for(j=0; j=k-l; j++) l[i][k] -= (l[i]U]*u[j][k]);for(j=k+l; jvn; j++) //计算U的第k彳亍元素 { u[k]U] = a[k][j]; for(i=0; i=k-l; i++) u[k]U] -= (l[k][i]*u[i]Ul)； u[k]U] /=l[k][k]; } } for(i=0; in; i++) // 解 Ly = b y[i]=b[i]； for(j=0; j=i-1; j++) y[i] -= (i[i]U]*yUl)； y[il /=l[i][i]； } for(i=n-l; i=0; i-) // 解 UX = Y { x[i]=y[i]； for(j=i+l; jn; j++) x[i] -= (u[i]U]*x[j]); } printf(MSolve...x_i = \nn); // 输出结果 for(i=0; in; i++) printf(n%f\n\ x[i]); return 0; } 2、实验题目 1) 4)、 4 1 0 1 4 2 0 2 4 0 0-1 (4,321) 6 0 1 2 _ 4 _ 2 -8 3 0 -2 兀2 11 -1 0 13 0 5 兀3 = 17 1 0 -4 9 2 ￡ -10 2 2 -1 1 7 _ 玉. _ 5 _ 3)、 (1丄丄1) 9 -4 (1 丄 0) 3 2 11 1 3 1 4 2 3 2 -7 -6 0 1 2). (100,1) _7