《矩阵的逆C语言代码及报告.docx》由会员分享,可在线阅读,更多相关《矩阵的逆C语言代码及报告.docx(3页珍藏版)》请在三一办公上搜索。
1、矩阵的逆C语言代码及报告矩阵的逆, 一,问题分析 设原矩阵是C,其逆矩阵是B 首先将B置为单位矩阵,然后用C和B构造一个新的矩阵(C,B), 用初等列变换, 将左边(即C)化为单位矩阵, 右半块(即B)就是原来C矩阵的逆。首先我们需要先检验矩阵是否可逆 二,源程序及注释。 #include #define M 3 #define N (2*M) int main int i,j,k; double aMM=1,2,3,2,2,1,3,4,3; double resultMM; double bMN; cout请输入矩阵的值:endl; for(i=0;iM;i+) for(j=0;jaij;
2、/输入矩阵 bij=aij; for(i=0;iM;i+) for(j=M;jN;j+) if(i=(j-M) bij=1; /在矩阵右侧加一个特征矩阵 else bij=0; for(i=0;iM;i+) if(bii=0) for(k=i;kM;k+) if(bkk!=0) for(int j=0;jN;j+) double temp; temp=bij; bij=bkj; bkj=temp; break; if(k=M) cout该矩阵不可逆!=i;j-) bij/=bii; for(k=0;kM;k+) if(k!=i) double temp=bki; for(j=0;jN;j+) bkj-=temp*bij; for(i=0;iM;i+) for(j=3;jN;j+) resultij-3=bij; for(i=0;iM;i+) for(j=0;jM;j+) coutresultij ; coutendl; return 0; 三,运行结果显示 四,调试和运行程序过程中产生的问题及采取的措施: