计算数值方法实验报告 _太原理工大学(2)

计算数值方法实验报告

} { if(i==k) continue; l[i]=a0[i][k]/a0[k][k]; } for (j=k+1;j<=n;j++) a[k][j]=a0[k][j]/a0[k][k]; b[k]=b0[k]/a0[k][k]; for(i=1;i<=n;i++) { if(i==k) continue; for(j=k+1;j<=n;j++) a[i][j]=a0[i][j]-l[i]*a0[k][j]; b[i]=b0[i]-l[i]*b0[k]; } displayA(); for(i=1;i<=n;i++) { for(j=k+1;j<=n;j++) a0[i][j]=a[i][j]; b0[i]=b[i]; } if(k==n) break; k++; }while(1); for(i=1;i<=n;i++) printf(\getchar(); 6

计算数值方法实验报告

六、结果讨论和分析 第一题答案： 第二题答案： 第三题答案： 7

计算数值方法实验报告

在编写这个程序的时候我感觉到了在编写的过程中的不易，在编写的过程中我改变初值然后运行，发现在初值的改变对运行时间的影响不大，结果没有影响。 在调试的过程中发现自己还是很粗心，容易犯很简单的错误。在今后的学习中应该做编写程序。 实验地点

8

综合楼六层606室 指导教师 王峥 计算数值方法实验报告

学院名称 学生姓名 课程名称 计算机科学与技术 某某 数值计算方法 专业班级 实验日期 实验题目 计算机 2011-6-13 学号 成绩 1111111111 实验三 线性方程组的迭代解法 一、 课题名称 线性方程组的迭代解法 使用雅可比迭代法或高斯-赛德尔迭代法对下列方程组进行求解。 ?10x1?x2?2x3?7.2???x1?10x2?2x3?8.3??x?x?5x?4.2123? 二、目的和意义 学习使用雅可比迭代法或高斯-赛德尔迭代法 三、计算公式 雅克比迭代法： 设线性方程组 Ax=b 的系数矩阵A可逆且主对角元素a11,a22,…,ann均不为零，令 D=diag(a11,a22,…,ann) 并将A分解成 A=(A-D)+D 从而线性方程组可写成 Dx=(D-A)x+b 则有迭代公式 x其中，B1=I-D-1(k+1)=B1x+f1 (k)A,f1=Db。 -1 四、结构程序设计 /雅可比 #include #include #include using namespace std; #define M 100/*最大允许的迭代次数*/ double fanshu(double *p,int n)/*求向量无穷范数的函数*/ {

9

计算数值方法实验报告

double n1; int i; n1=fabs(p[0]); for(i=1;in1) n1=fabs(p[i]); } return(n1); } void jac(double **s,double *g,int size,double eps)/*雅可比迭代函数*/ { int i,j,k; double **b,*x0,*x1,*x1_x0,n2,t; /*b用来表示迭代矩阵,x1_x0表示迭代前后向量x0与x1的差向量*/ b=new double*[size]; /*动态分配二维数组*/ for(i=0;i

10