数值分析报告.docVIP

数值分析实验报告 姓名：*** 学号：*** 任课教师：*** 实验一 .线性方程组的数值实验 第一题 一、问题 对从2到9的每一个n值，解n阶方程组 。在这里和如下定义： ， 其中 解释发生的现象。 实验原理 Gauss消去法： 设矩阵，向量，把方程组写成 （1）顺序消去过程 假设，令 做运算将[]变换为 [] =, 其中 这样就得到与等价的方程组.其中[]第一行与 第一行相同，而其第一列对角线以下元素已经消为0。 假设有≠0，k=1,2,…，n-1,则经n-1步消元后得到 [|]= 对应等价方程组=。 由于上述初等变换系数矩阵行列式的值不变，有，所以。 回代过程 解=从最后一个方程开始，用逐次向后代入的方法，这称为“回代过程”，计算公式是： ， 如上述过程，由顺序消去过程和回代过程合起来解方程组Ax=b的方法，就是Gauss消去法。 三、实验步骤 定义题目中的函数p(x)； 按照题目要求初始化系数矩阵； 实现Gauss消去法顺序消元过程； 实现Gauss消去法回代过程； 输出不同的n值的情况下得到的线性方程组的解向量； 6、在求出一个解后，设，令，输出向量 ，可以观察从2到9的每一个n值，用Gauss消去法解方程的误差。 四、数据流向图 五．模块关系图 六．符号引用表 标志号 数学符号 类型 入 出 作用 N n 整型 方程组阶数 A A 双精度型 系数矩阵 A C 双精度型 常数B B B 双精度型 解X的中间值 B D 双精度型 存放中间值 B X 双精度型 * 解X B Y 双精度型 * 解的误差向量 源程序代码 #includeiostream #includestdio.h #includeiomanip using namespace std; #includemath.h #define N 9 // 矩阵的最大阶数 double p(double x, int n){ //定义p(x) double z1 = 14+n*(12+3*n); double z2 = -7+(n*n)*z1; double z3 = 2+x*x*z2; double z4 = x*x/24*z3; return z4; } double A[N+1][N+1]; //系数矩阵A double B[N+1]; //常数B double C[N+1][N+1]; //解X double D[N+1]; double X[N+1]; double Y[N+1]; void matrix(int n) //对于不同的值，初始化系数矩阵 { int i,j; for(i=1;i=n;i++) { for(j=1;j=n;j++) { A[i][j]=double(1)/double((i+j-1)); B[i]=p(i+n-1,n)-p(i-1,n); C[i][j]=A[i][j]; D[i]=B[i]; } } } void Gauss(int n) //高斯消去法 { int i,j,k; double s; for(k=1;k=n;k++) { if(A[k][k]==0) break; for(i=k+1;i=n;i++) { s=A[i][k]/A[k][k]; A[i][k]-=A[k][k]*s; B[i]-=B[k]*s; for(j=k+1;j=n;j++) A[i][j]-=A[k][j]*s; } } } void reback(int n) //回代求解X { int i,j; X[n]=B[n]/A[n][n]; for(i=n-1;i=1;i--) { double x =0; for(j=i+1;j=n;j++) x+=A[i][j]*X[j]; X[i]=(B[i]-x)/A[i][i]; } cout当n=n时,方程组解为:endl; cou