《C语言求矩阵的逆》由会员分享，可在线阅读，更多相关《C语言求矩阵的逆（8页珍藏版）》请在金锄头文库上搜索。
1、#include#includevoidjiafa（）intm,n;floata2020,b2020,c2020;inti,j;printf（”请输入矩阵行数：”）；scanf（%d,&m）;printf（请输入矩阵列数：”）；scanf（%d,&n）；printf（请输入第一个矩阵：”）；for（i=0；im；i+）for（j=0；jn；j+）scanf（%f,&aij）；printf（请输入第二个矩阵：”）；for（i=0；im；i+）for（j=0；jn；j+）scanf（%f,&bij）；printf（矩阵相加的结果为：”）；for（i=0；im；i+）for（j=0；jn；j+）ci
2、j=aij+bij；printf（%4f,cij）；printf（n）；voidjianfa（）intm,n；floata2020,b2020,c2020；inti,j；printf（请输入矩阵行数：”）；scanf（%d,&m）;printf（请输入矩阵列数：”）；scanf（%d,&n）;printf（请输入第一个矩阵：”）；for（i=0；im；i+）for（j=0；jn；j+）scanf（%f,&aij）；printf（请输入第二个矩阵：”）；for（i=0；im；i+）for（j=0；jn；j+）scanf（%f,&bij）；printf（矩阵相减的结果为：）；for（i=0；im；
3、i+）for（j=0；jn；j+）cij=aij-bij；printf（%4f,cij）；printf（n）；voidchengfa（）intm,n；floats；floata2020,b2020,c2020；inti,j,k；printf（请输入矩阵行数：”）；scanf（%d,&m）；printf（请输入矩阵列数：”）；scanf（%d,&n）；printf（请输入第一个矩阵：”）；for（i=0；im；i+）for（j=0；jn；j+）scanf（%f,&aij）；printf(请输入第二个矩阵：”);for(i=0;in;i+)for(j=0;jm;j+)scanf(%4f,&bij)
4、;for(i=0;im;i+)for(j=0;jm;j+)s=0;for(k=0;kn;k+)s=s+aik*bkj;cij=s;for(i=0;im;i+)for(j=0;jm;j+)printf(%4f,cij);printf(n);voidzhuanzhi()intm,n;floata2020,b2020;inti,j;printf(请输入矩阵行数：”)；scanf(%d,&m);printf(请输入矩阵列数：”)；scanf(%d,&n);printf(请输入一个矩阵：”)；for(i=0;im;i+)for(j=0；jn；j+)scanf(%f,&aij)；for(i=0;im;i+
5、)for(j=0;jn;j+)bij=aji;printf(%4f,bij);printf(n);voidqiuni()intN;printf(输入矩阵的阶数N:n);scanf(%d,&N);floata1010,b1020,c1010,t;inti,j,m;printf(请输入行列式不为0的矩阵A(%d阶):n,N);矩阵A的各元素存入二维数组a中。for(i=0;iN;i+)for(j=0;jN;j+)scanf(%f,&aij);增广矩阵(A|E)存入二维数组b中for(i=0;iN;i+)for(j=0;jN;j+)bij=aij;for(i=0;iN;i+)for(j=N;j2*N
6、;j+)bij=0;for(i=0;iN;i+)biN+i=1;for(m=0;mm)/实现交换。bim=t;/交换其它各列相应位置的元素for(j=0;jm;j+)t=bmj;bmj=bij;bij=t;for(j=m+1;j2*N;j+)t=bmj;bmj=bij;bij=t;for(i=m+1;i=m;j-)bij-=bim*bmj/bmm;m=0时，将第一行的-biO/bOO倍加到以下各行。这样以下每行第一个元素bi0就为0。for(j=2*N-1;j=m;j-)bmj/=bmm;对第m行作行变换，同除以bmm，使bmm为1。printf(第一步变换后得到的增广矩阵为：n);for(i
7、=0;iN;i+)for(j=0;j0)for(i=0;i=m;j-)/千万注意，此处j必须递减，否则bim先变为0，后面的计算就无效！bij-=bim*bmj;m-;printf(最后得到的增广矩阵为：n);for(i=0;iN;i+)for(j=0;j2*N;j+)printf(%3.5f,bij);printf(n);实现了：每个i对应一个换行。for(i=0;ivN;i+)将逆矩阵存入二维数组c中。for(j=0;jN;j+)cij=biN+j;printf(故逆矩阵为：n);for(i=0;ivN;i+)for(j=0;jvN;j+)printf(%3.5f,cij);printf(n);实现了：每个i对应一个换行。main()intw;printf(1矩阵加法n);printf(2矩阵减法n);printf(3矩阵乘法n”)；printf(4矩阵转置n);printf(5矩阵求逆n”)；printf(n);printfC请选择要进行的运算：)；scanf(%d,&w);switch(w)case1:jiafa();break;case2:jianfa();break;case3:chengfa();break;case4:zhuanzhi();break;case5:qiuni();break;return0;
