矩阵的逆C语言代码及报告.docVIP

矩阵的逆， 一，问题分析 设原矩阵是C，其逆矩阵是B首先将B置为单位矩阵，然后用C和B构造一个新的矩阵(C,B), 用初等变换, 将左边(即C)化为单位矩阵, 右半块(即B)就是原来C矩阵的逆。 二，源程序及注释。 #include iostream.h#define M 3#define N (2*M)int main(){ int i,j,k; double a[M][M]={1,2,3,2,2,1,3,4,3}; double result[M][M]; double b[M][N]; cout请输入矩阵的值（默认大小为3*3的矩阵）：endl; for(i=0;iM;i++) { for(j=0;jM;j++) { cina[i][j]; b[i][j]=a[i][j]; } } for(i=0;iM;i++) { for(j=M;jN;j++) { if(i==(j-M)) { b[i][j]=1; //在矩阵右侧加一个特征矩阵（单位矩阵） } else { b[i][j]=0;  } } } for(i=0;iM;i++) { if(b[i][i]==0) { for(k=i;kM;k++) { if(b[k][k]!=0) { for(int j=0;jN;j++) { double temp; temp=b[i][j]; b[i][j]=b[k][j]; b[k][j]=temp; } break; } } if(k==M) { cout该矩阵不可逆！endl; } } for(j=N-1;j=i;j--) { b[i][j]/=b[i][i]; } for(k=0;kM;k++) { if(k!=i) { double temp=b[k][i]; for(j=0;jN;j++) { b[k][j]-=temp*b[i][j]; } } } } for(i=0;iM;i++) { for(j=3;jN;j++) { result[i][j-3]=b[i][j]; } } for(i=0;iM;i++) { for(j=0;jM;j++) { coutresult[i][j] ; } coutendl; } return 0;} 三，运行结果显示 四，调试和运行程序过程中产生的问题及采取的措施：