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

高斯消元法求N阶矩阵的逆矩阵题目：用高斯消元法求N阶矩阵的逆矩阵 C语言实现 #include #include long int const N=1000; //定义矩阵最大为1000阶 int n; //n表示矩阵的行数和列数。 double juzhen[N][N]; //定义一个1000阶矩阵 double danwei[N][N]; //定义一个单位矩阵 bool zhaozuidazhi(int s)//定义一个布尔型从s行到n行选择最大的元素作为主元的函数。 { int ...

题目：用高斯消元法求N阶矩阵的逆矩阵 C语言实现 #include #include long int const N=1000; //定义矩阵最大为1000阶 int n; //n 表示矩阵的行数和列数。 double juzhen[N][N]; //定义一个1000阶矩阵 double danwei[N][N]; //定义一个单位矩阵 bool zhaozuidazhi(int s)//定义一个布尔型从s行到n行选择最大的元素作为主元的函数。 { int i,j,A; double mas,temp; //temp为中间变量，实现行交换 mas=fabs(juzhen[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=0;i--) { mas=juzhen[i][s]; juzhen[i][s]=0; for(j=0;j=0;i--) { HH(i); } printf("\n原始矩阵的逆矩阵为:\n");//输出逆矩阵语句 for(i=0;i #include #include double *inv(double *a,int n) { int *is,*js,i,j,k,l,u,v; double d,p; is=(int*)malloc(n*sizeof(int)); //用于向内存申请空间，分配长度为n*sizeof(int)_bytes字节的内存块// //(int *)是用于强制类型转换为用户定义的指针变量的类型// js=(int*)malloc(n*sizeof(int)); for (k=0; k<=n-1; k++) { d=0.0; for (i=k; i<=n-1; i++) for (j=k; j<=n-1; j++) { l=i*n+j; p=fabs(a[l]); if (p>d) { d=p; is[k]=i; js[k]=j; } } if (d+1.0==1.0) { free(is); free(js); //内存不再使用时，应使用free()函数将内存块is和js释放。// printf("此矩阵没有逆矩阵\n"); return NULL; //返回空指针NULL// } if (is[k]!=k) for (j=0; j<=n-1; j++) { u=k*n+j; v=is[k]*n+j; p=a[u]; a[u]=a[v]; a[v]=p; } if (js[k]!=k) for (i=0; i<=n-1; i++) { u=i*n+k; v=i*n+js[k]; p=a[u]; a[u]=a[v]; a[v]=p; } l=k*n+k; a[l]=1.0/a[l]; for (j=0; j<=n-1; j++) if (j!=k) { u=k*n+j; a[u]=a[u]*a[l]; } for (i=0; i<=n-1; i++) if (i!=k) for (j=0; j<=n-1; j++) if (j!=k) { u=i*n+j; a[u]=a[u]-a[i*n+k]*a[k*n+j]; } for (i=0; i<=n-1; i++) if (i!=k) { u=i*n+k; a[u]=-a[u]*a[l]; } } for (k=n-1; k>=0; k--) { if (js[k]!=k) for (j=0; j<=n-1; j++) { u=k*n+j; v=js[k]*n+j; p=a[u]; a[u]=a[v]; a[v]=p; } if (is[k]!=k) for (i=0; i<=n-1; i++) { u=i*n+k; v=i*n+is[k]; p=a[u]; a[u]=a[v]; a[v]=p; } } free(is); free(js); return a; } void display(double M[],int n) //控制输出// { int i,j; double *result=NULL,*A=NULL; printf("\n\n输入的原始矩阵为:\n"); for(i=0;i /* Test */ void main() { double matrix[3][3]={{1,-5,-2},{-1,3,1},{3,-4,-1}}; double matrix1[3][3]={{0,0,0},{-1,3,1},{3,-4,-1}}; double matrix2[4][4]={{3,2,1,1},{2,3,5,9},{-1,2,5,-2},{1,0,-1,3}}; display(matrix[0],3); display(matrix1[0],3); display(matrix2[0],4); printf("push any key to be end ."); getch(); }

                    本文档为【高斯消元法求N阶矩阵的逆矩阵】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
                    图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。 
 该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。

                    [版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。

                    本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。

                    网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
                

下载需要：免费已有0 人下载

立即下载

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

你可能还喜欢