求解线性方程组.docVIP

班级 学号：0811081037 姓名：袁志勇 数值分析实验报告(一) 一、实验名称：求解线性方程组 二、实验问题： 从2到9的每一个n值，解n阶方程组Ax=b。这里A和b如下定义： aij=(i+j-1)-1,bi=p(n+i-1)-p(i-1) 式中 p(x)=(x2/24)(2+x2(-7+n2(14+n(12+3n)))) 解释发生的现象。 三、计算公式（数学模型）： Gauss消去法是将方程组Ax=b，用逐次消去未知数的方法，化为等价的三角方程组，再回代求解的方法。 设Ax=b，A， 如果，则可通过Gauss消去法将Ax=b约化为等价的三角形方程组 且计算公式为： (a)消元计算(k=1,2,…,n) (b)回代计算： 四、结构程序设计： 1、数据流向图 入Gauss消元回代过程 入 Gauss消元 回代过程 出 2、模块关系和模块流程图 Gauss Gauss消元 回代过程 控制模块 回代模块 Gauss消元模块 五、计算程序 C++程序 //改进的平方根法 #includeiostream.h double p(double x,double n){ return x*x/24*(2+x*x*(-7+n*n*(14+n*(12+3*n)))); } void main(void){ cout*************************线性方程组求解的数值实验题*****************************\n; double A[9][9]={0}; double L[9][9]={0}; double B[9]={0}; double d[9]={0}; double x[9]={0}; double y[9]={0}; double dd=0; int n,i,j,k; for(n=2;n10;n++){ coutn阶方程组Ax=b的解为\n; for(i=0;in;i++){ for(j=0;jn;j++){ A[i][j]=1.0/(i+j+1); //赋Ax=b中系数A的初值(A是对称阵) } B[i]=p((double)n+i,(double)n)-p((double)i,(double)n); //赋Ax=b中系数b的初值(调用函数p) } d[0]=A[0][0]; //d[]用来记录A[k][k] for(i=0;in;i++) L[i][0]=A[i][0]/d[0]; //LDL分解中L[i][j]的第一列 for(i=1;in;i++){ for(k=0;ki;k++) dd+=L[i][k]*L[i][k]*d[k]; d[i]=A[i][i]-dd; //LDL分解中的d[i] dd=0; for(j=i;jn;j++){ for(k=0;kj;k++) dd+=L[i][k]*d[k]*L[j][k]; //LDL分解中除第一列外的L[i][j] L[i][j]=(A[i][j]-dd)/d[j]; dd=0; } } y[0]=B[0]; //对y的回代过程 for(i=1;in;i++){ for(k=0;ki;k++) dd+=L[i][k]*y[k]; y[i]=B[i]-dd; dd=0; } x[n-1]=y[n-1]/d[n-1]; //对x的回代过程 for(i=n-2;i=0;i--){ for(k=i+1;kn;k++) dd+=L[k][i]*x[k]; x[i]=y[i]/d[i]-dd; dd=0; } for(i=0;in;i++){ coutxi+1=x[i]\n; } cout\n; } } //在C++数组中是以 开始的，所以i，j，k取值相应变到从0开始