图像处理之基础---矩阵求逆实现

最近做一个加密算法遇到需要计算矩阵的逆，闲着无聊，记录一下，以后免得再麻烦。 [cpp] view plaincopyprint? #include #include #include #define MAX 20 #define E 0.000000001 /** * 计算矩阵src的模 */ double calculate_A( double src[][MAX], int n ) { int i,j,k,x,y; double tmp[MAX][MAX], t; double result = 0.0; if( n == 1 ) { return src[0][0]; } for( i = 0; i < n; ++i ) { for( j = 0; j < n - 1; ++j ) { for( k = 0; k < n - 1; ++k ) { x = j + 1; y = k >= i ? k + 1 : k; tmp[j][k] = src[x][y]; } } t = calculate_A( tmp, n - 1 ); if( i % 2 == 0 ) { result += src[0][i] * t; } else { result -= src[0][i] * t; } } return result; } /** * 计算伴随矩阵 */ void calculate_A_adjoint( double src[MAX][MAX], double dst[MAX][MAX], int n ) { int i, j, k, t, x, y; double tmp[MAX][MAX]; if( n == 1 ) { dst[0][0] = 1; return; } for( i = 0; i < n; ++i ) { for( j = 0; j < n; ++j ) { for( k = 0; k < n - 1; ++k ) { for( t = 0; t < n - 1; ++t ) { x = k >= i ? k + 1 : k ; y = t >= j ? t + 1 : t; tmp[k][t] = src[x][y]; } } dst[j][i] = calculate_A( tmp, n - 1 ); if( ( i + j ) % 2 == 1 ) { dst[j][i] = -1*dst[j][i]; } } } } /** * 得到逆矩阵 */ int calculate_A_inverse( double src[MAX][MAX], double dst[MAX][MAX], int n ) { double A = calculate_A( src, n ); double tmp[MAX][MAX]; int i, j; if ( fabs( A - 0 ) <= E ) { printf("不可能有逆矩阵！\n"); return 0; } calculate_A_adjoint( src, tmp, n ); for( i = 0; i < n; ++i ) { for( j = 0; j < n; ++j ) { dst[i][j] = (double)( tmp[i][j] / A ); } } return 1; } /** * 输出矩形查看 */ void print_M( double M[][MAX], int n ) { int i, j; for ( i = 0; i < n; ++i ) { for ( j = 0; j < n; ++j ) { printf("%lf ", M[i][j]); } printf("\n"); } } /** * main */ int main() { /** * 测试矩阵 */ double test[MAX][MAX], dst[MAX][MAX]; int n = 3; int is_exist; /** * 构造最简单的： * 1, 0, 0, * 0, 2, 0, * 0, 0, 5 */ memset( test, 0, sizeof( test ) ); test[0][0] = 1.0; test[1][1] = 2.0; test[2][2] = 5.0; is_exist = calculate_A_inverse( test, dst, n ); if ( is_exist ) { print_M(dst, n); } else { printf("不存在!\n"); } return 0; } http://blog.csdn.net/shanshanpt/article/details/16820325