矩阵求逆算法(Chapter28)

1. /**
2. *   求矩阵A的逆矩阵Ai
3. *@param     A       源矩阵
4. *@param     Ai      逆矩阵
5. *@param     size    矩阵的大小
6. *@return     求解成功返回非零值，失败返回零
7. */
8. int InverseMatrix(double **Ai, double **A, int size)
9. {
10.     int i, j;
11.     double *b, *x;
12.     b = (double*)malloc(sizeof(double)*size);
13.     x = (double*)malloc(sizeof(double)*size);
15.     for (i = 0; i < size; i++)
16.     {
17.         memset(b, 0, sizeof(double)*size);
18.         b[i] = 1;
19.         if (!LinearEquation(A, x, b, size))
20.         {
21.             free(b);
22.             free(x);
23.             return 0;
24.         }
25.         for (j = 0; j < size; j++)
26.             Ai[j][i] = x[j];
27.     }
28.     free(b);
29.     free(x);
30.     return 1;
31. }