#include<stdio.h> #include<math.h> #include<stdlib.h> #define dx 100 struct bb { int m; int n; int hl[dx][dx]; int jk[dx][dx]; }; double det(struct bb *A, int n); double algebraic_cofactor(struct bb *A, struct bb *B, int row, int col); void adjoint(struct bb *A, struct bb *B); void inverse(struct bb *A,double inv[][dx],int N); int main() { struct bb A; int m,n; printf("输入几行几列：\n"); scanf("%d %d",&m,&n); A.m = m; A.n = n; printf("请输入矩阵：\n"); for(int i = 0; i < A.m; i++) { for(int j = 0; j < A.n; j++) { scanf("%d", &A.hl[i][j]); } } double inv[dx][dx]; int N = A.m; // Assuming square matrix // 计算逆矩阵 inverse(&A, inv, N); // 输出结果 printf("逆矩阵：\n"); for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { printf("%.2f ", A.hl[i][j]); } printf("\n"); } return 0; } double det(struct bb *A, int n) { double sum=0; if(n==1) { return A->hl[0][0]; } else if(n==2) { return A->hl[0][0]*A->hl[1][1]-A->hl[0][1]*A->hl[1][0]; } int i,j,k; struct bb *mybb = (struct bb *)malloc(sizeof(struct bb)); for(k=0;k<n;k++) { double b[dx][dx]; for(i=1;i<n;i++) { for(j=0;j<n;j++) { if(j<k) { b[i-1][j]=A->hl[i][j]; } else if(j>k) { b[i-1][j-1]=A->hl[i][j]; } } } mybb->m = n - 1; mybb->n = n - 1; for(i=0;i<mybb->m;i++) { for(j=0;j<mybb->n;j++) { mybb->hl[i][j] = b[i][j]; } } double detb=det(mybb,n-1); sum+=A->hl[0][k]*pow(-1,k)*detb; } free(mybb); return sum; } double algebraic_cofactor(struct bb *A, struct bb B, int row, int col) { int i,j,m=0,n=0,M=A->m; double sign; if((row+col)%2==0) { sign=1; } else { sign=-1; } for(i=0;i<M;i++) { if(i!=row) { for(j=0;j<M;j++) { if(j!=col) { B->jk[m][n]=A->hl[i][j]; n++; } } m++; n=0; } } double detb=det(B,M-1); return signdetb; } void adjoint(struct bb *A, struct bb *B) { int i,j,M=A->m; for(i=0;i<M;i++) { for(j=0;j<M;j++) { B->hl[j][i]=algebraic_cofactor(A,B,i,j); //注意这里是 hl[j][i] 而不是 hl[i][j] } } } void inverse(struct bb *A,double inv[][dx],int N) { // 构造伴随矩阵 struct bb B; B.m = N; B.n = N; adjoint(A, &B); // 计算行列式的值 double dets=det(A,N); // 判断行列式是否为零 if(dets == 0) { printf("该矩阵不可逆！\n"); return; } // 计算逆矩阵 for(int i=0;i<N;i++) { for(int j=0;j<N;j++) { inv[i][j] = B.hl[i][j] / dets; } } }修改这个代码找出错误并改正
3301
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
