高斯消元法求N阶矩阵的逆矩阵

《高斯消元法求N阶矩阵的逆矩阵》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、题目：用高斯消元法求N阶矩阵的逆矩阵C语言实现#include#includelongintconstN=1000;//定义矩阵最大为1000阶intn;//n表示矩阵的行数和列数。doublejuzhen[N][N];//定义一个1000阶矩阵doubledanwei[N][N];//定义一个单位矩阵boolzhaozuidazhi(ints)//定义一个布尔型从s行到n行选择最大的元素作为主元的函数。{inti,j,A;doublemas,temp;//temp为中间变量，实现行交换mas=fabs(juzhen[

2、s][s]);//数学函数:fabs　　功能：求浮点数juzhen[s][s]的绝对值//计算juzhen[s][s],当juzhen[s][s]不为负时返回juzhen[s][s]，否则返回-juzhen[s][s]A=s;for(i=s+1;i

3、zhen[A][j]=temp;temp=danwei[s][j];danwei[s][j]=danwei[A][j];danwei[A][j]=temp;}returntrue;}voidjisuan(ints)//消元计算{inti,j;doublemas=juzhen[s][s],r;for(i=s+1;i

4、r;danwei[i][j]=danwei[i][j]-danwei[s][j]*r;}}}voidHH(ints){inti,j;doublemas;mas=juzhen[s][s];for(i=s;i=0;i--){mas=juzhen[i][s];juzhen[i][s]=0;for(j=0;j

5、mas*danwei[s][j];}}intmain(){inti,j;zl:printf("第一行输入矩阵的阶数，然后回车输入原始矩阵:");scanf("%d",&n);for(i=0;i

6、改指针指向的内容。所以，//把x取地址再传过去。&是取juzhen[i][j]的地址,%lf是说juzhen[i][j]是一个双精度浮点型.}}//将原始矩阵转成上三角矩阵for(i=0;i=0;i--

7、){HH(i);}printf("原始矩阵的逆矩阵为:");//输出逆矩阵语句for(i=0;i

8、225原始矩阵的逆矩阵为:5.00-2.00-2.001.00第一行输入矩阵的阶数，然后回车输