c语言求矩阵逆函数 c语言求矩阵的逆矩阵( 三 ) _矩阵

j1++;j2++;
if(j1==j)j1++;k1=-1;k2=-1;
while(k2n-1){
k1++;
k2++;
if(k1==i)k1++;
*(*(po+j2)+k2)=*(*(p+j1)+k1);
}
}
tramform[i][j]=sgnx(2+i+j) * danx(po,n-1);
}
}
for(i=0;in;i++){
for(j=0;jn;j++){
printf("%4d",tramform[i][j]);
}
printf("\n");
}
}
}
}
//运行结果//希望对你有帮助
C语言编程：编写一个函数求逆矩阵#include stdio.h
#include stdlib.h
#include malloc.h
void MatrixOpp(double *A, int m, int n, double* invmat);
void MatrixInver(double *A, int m, int n, double* invmat);
double Surplus(double A[], int m, int n);
int matrix_inv(double* p, int num, double* invmat);
void MatrixOpp(double A[], int m, int n, double* invmat)
{
int i, j, x, y, k;
double *SP = NULL, *AB = NULL, *B = NULL, X;
SP = (double *) malloc(m * n * sizeof(double));
AB = (double *) malloc(m * n * sizeof(double));
B = (double *) malloc(m * n * sizeof(double));
X = Surplus(A, m, n);
X = 1 / X;
for (i = 0; im; i++)
for (j = 0; jn; j++)
{
for (k = 0; km * n; k++)
B[k] = A[k];
{
for (x = 0; xn; x++)
B[i * n + x] = 0;
for (y = 0; ym; y++)
B[m * y + j] = 0;
B[i * n + j] = 1;
SP[i * n + j] = Surplus(B, m, n);
AB[i * n + j] = X * SP[i * n + j];
}
}
MatrixInver(AB, m, n, invmat);
free(SP);
free(AB);
free(B);
}
void MatrixInver(double A[], int m, int n, double* invmat)
{
int i, j;
double *B = invmat;
for (i = 0; in; i++)
for (j = 0; jm; j++)
B[i * m + j] = A[j * n + i];
}
double Surplus(double A[], int m, int n)
{
int i, j, k, p, r;
double X, temp = 1, temp1 = 1, s = 0, s1 = 0;
if (n == 2)
{
for (i = 0; im; i++)
for (j = 0; jn; j++)
if ((i + j) % 2)
temp1 *= A[i * n + j];
else
temp *= A[i * n + j];
X = temp - temp1;
}
else
{
for (k = 0; kn; k++)
{
for (i = 0, j = k; im, jn; i++, j++)
temp *= A[i * n + j];
if (m - i)
{
for (p = m - i, r = m - 1; p0; p--, r--)
temp *= A[r * n + p - 1];
}
s += temp;
temp = 1;
}
for (k = n - 1; k = 0; k--)
{
for (i = 0, j = k; im, j = 0; i++, j--)
temp1 *= A[i * n + j];
if (m - i)
{
for (p = m - 1, r = i; rm; p--, r++)
temp1 *= A[r * n + p];
}
s1 += temp1;
temp1 = 1;
}
X = s - s1;
}
return X;
}
int matrix_inv(double* p, int num, double* invmat)
{
if (p == NULL || invmat == NULL)
{
return 1;
}
if (num10)
{
return 2;
}
MatrixOpp(p, num, num, invmat);
return 0;
}
int main()
{
int i, j;
int num;
double *arr=NULL;
double *result=NULL;
int flag;
printf("请输入矩阵维数：\n");
scanf("%d",num);
arr=(double *)malloc(sizeof(double)*num*num);
result=(double *)malloc(sizeof(double)*num*num);
printf("请输入%d*%d矩阵：\n", num, num);
for (i = 0; inum; i++)
{
for (j = 0; jnum; j++)
{
scanf("%lf", arr[i * num + j]);
}
}
flag = matrix_inv(arr, num, result);
if(flag==0)
{
printf("逆矩阵为：\n");
for (i = 0; inum * num; i++)
{
printf("%lf\t ", *(result + i));
if (i % num == (num - 1))
printf("\n");
}
}
else if(flag==1)
{
printf("p/q为空\n");
}
else
{
printf("超过最大维数\n");
}
system("PAUSE");
free(arr);
free(result);
return 0;

c语言求矩阵逆函数 c语言求矩阵的逆矩阵( 三 )

推荐阅读

双硬盘设置第一启动盘双硬盘设置，双硬盘有什么好处

阳光玫瑰葡萄要剥皮吗

stata多元回归分析结果

文竹的繁殖方法

乌贼的功效与作用及食疗偏方

藕粉可以杀死幽门螺旋杆菌吗？

干茶树菇泡发窍门

我读毛姆的四大名篇|我读毛姆的四大名篇第三本（讲讲《面纱》的故事）

如何做软件代理商，软件代理销售

黄桃太硬怎么变软

西藏哪里好玩

决明子茶和苦荞茶能一起喝吗

聚财招财旺财山水画

草莓干跟草莓哪个更有营养

气色不好怎么调理

苹果14配置参数表图片，iPhone5iPhone4SiPhone4的参数是什么

【|【中原焦点团队】贾素凤网络初14 坚持分享第 368天

ultr三星旗舰阵容有变数？传 Galaxy Note 20 成为 “系列末代机”

精灵盛典魔法师怎么加点

生腌为什么不能配啤酒