The invention relates to a matrix inverse operation method. The method comprises the steps of 1, conducting column pivoting LU decomposition, wherein a source matrix A is decomposed into a unit lower triangular matrix L, an upper triangular matrix U and a permutation matrix P according to the formula PA=LU; 2, conducting triangular matrix inversion, wherein the inverse matrix L-1 of the matrix L is obtained through matrix inversion, and matrix inversion is conducted on the transposed matrix of the matrix U and then transposition is conducted to obtain U-1; 3, finally conducting matrix multiplication, wherein the matrix U-1 and the matrix L-1 are multiplied, and column transformation is conducted on the matrix multiplication result according to the permutation matrix P to obtain a source matrix A-1. The method has the advantages that by using the column pivoting LU decomposition algorithm, the time complexity of the matrix inversion algorithm is effectively reduced, parallelizability of matrix inversion operation is improved, time for matrix inversion operation is shortened, matrix inversion operation of any order can be conducted, and the number of hardware resources can be increased or reduced according to count requirements of operation so that practical application requirements can be better met.