计算方法-解线性方程组的直接法实验报告方案.docVIP

计算方法实验二实验报告1．【实验题目2．【实验目的3．【实验内容实验#includeiostream #includecmath #define N 20 using namespace std; void load(); float a[N][N]; int m; int main() { int i,j; int c,k,n,p,r; float x[N],l[N][N],s,d; cout下面请输入未知数的个数m=; cinm; coutendl; cout请按顺序输入增广矩阵a:endl; load(); for(i=0;im;i++) { for(j=i;jm;j++) c=(fabs(a[j][i])fabs(a[i][i]))?j:i; /*找列最大元素*/ for(n=0;nm+1;n++) { s=a[i][n]; a[i][n]=a[c][n]; a[c][n]=s; } /*将列最大数防在对角线上*/ for(p=0;pm+1;p++) couta[i][p]\t; coutendl; for(k=i+1;km;k++) { l[k][i]=a[k][i]/a[i][i]; for(r=i;rm+1;r++) /*化成三角阵*/ a[k][r]=a[k][r]-l[k][i]*a[i][r]; } } x[m-1]=a[m-1][m]/a[m-1][m-1]; for(i=m-2;i=0;i--) { d=0; for(j=i+1;jm;j++) d=d+a[i][j]*x[j]; x[i]=(a[i][m]-d)/a[i][i]; /*求解*/ } cout该方程组的解为:endl; for(i=0;im;i++) coutx[i]=x[i]\t; // system(pause); return 0; } void load() { int i,j; for(i=0; im;i++) for(j=0; jm+1; j++) cina[i][j]; } 一般消去法 #includestdio.h void solve(float l[][100],float u[][100],float b[],float x[],int n) { int i,j; float t,s1,s2; float y[100]; for(i=1;i=n;i++) /* 第一次回代过程开始 */ { s1=0; for(j=1;ji;j++) { t=-l[i][j]; s1=s1+t*y[j]; } y[i]=(b[i]+s1)/l[i][i]; } for(i=n;i=1;i--) /* 第二次回代过程开始 */ { s2=0; for(j=n;ji;j--) { t=-u[i][j]; s2=s2+t*x[j]; } x[i]=(y[i]+s2)/u[i][i]; } } void main() { float a[100][100],l[100][100],u[100][100],x[100],b[100]; int i,j,n,r,k; float s1,s2; for(i=1;i=99;i++)/*将所有的数组置零，同时将L矩阵的对角值设为1*/ for(j=1;j=99;j++) { l[i][j]=0,u[i][j]=0; if(j==i) l[i][j]=1; } printf (input n:\n);/*输入方程组的个数*/ scanf(%d,n); printf (input array A:\n);/*读取原矩阵A*/ for(i=1;i=n;i++) for(j=1;j=n;j