/**
* 求矩阵A的逆矩阵Ai
*@param A 源矩阵
*@param Ai 逆矩阵
*@param size 矩阵的大小
*@return 求解成功返回非零值,失败返回零
*/
int InverseMatrix(double **Ai, double **A, int size)
{
int i, j;
double *b, *x;
b = (double*)malloc(sizeof(double)*size);
x = (double*)malloc(sizeof(double)*size);
for (i = 0; i < size; i++)
{
memset(b, 0, sizeof(double)*size);
b[i] = 1;
if (!LinearEquation(A, x, b, size))
{
free(b);
free(x);
return 0;
}
for (j = 0; j < size; j++)
Ai[j][i] = x[j];
}
free(b);
free(x);
return 1;
}
① 凡本网注明稿件来源为"原创"的所有文字、图片和音视频稿件,版权均属本网所有。任何媒体、网站或个人转载、链接转贴或以其他方式复制发表时必须注明"稿件来源:我考网",违者本网将依法追究责任;
② 本网部分稿件来源于网络,任何单位或个人认为我考网发布的内容可能涉嫌侵犯其合法权益,应该及时向我考网书面反馈,并提供身份证明、权属证明及详细侵权情况证明,我考网在收到上述法律文件后,将会尽快移除被控侵权内容。