[原创]:[LSP函数]求逆矩阵。
;;求矩阵的逆矩阵;(defun ea:matrix_inverse (matrix /m0 m1 m2 m3 m00 m01m02 m03 m10 m11 m12 m13m20 m21 m22 m23 _a ea:determinant a11 a12 a21 a22a33 a34 a41 a42 a43 a13a23 a31 a32 ) ;;行列式展开,计算代数余子式 (defun ea:determinant(a00 a01 a02 a10 a11 a12 a20 a21 a22 /) (- (+ (* a00 a11 a22) (* a01 a12 a20) (* a02 a21 a10) ) (* a02 a11 a20) (* a01 a10 a22) (* a00 a21 a12) ) ) (setqm0 (car matrix)m1 (cadr matrix)m2 (nth 2 matrix)m3 (last matrix) ) (setqm00 (car m0)m01 (cadr m0)m02 (nth 2 m0)m03 (last m0) ) (setqm10 (car m1)m11 (cadr m1)m12 (nth 2 m1)m13 (last m1) ) (setqm20 (car m2)m21 (cadr m2)m22 (nth 2 m2)m23 (last m2) ) (setqa11 (ea:determinant m11 m12 m13 m21 m22 m23 0.0 0.0 1.0)a12 (- (ea:determinant m10 m12 m13 m20 m22 m23 0.0 0.0 1.0))a13 (ea:determinant m10 m11 m13 m20 m21 m23 0.0 0.0 1.0)a21 (- (ea:determinant m01 m02 m03 m21 m22 m23 0.0 0.0 1.0))a22 (ea:determinant m00 m02 m03 m20 m22 m23 0.0 0.0 1.0)a23 (- (ea:determinant m00 m01 m03 m20 m21 m23 0.0 0.0 1.0))a31 (ea:determinant m01 m02 m03 m11 m12 m13 0.0 0.0 1.0)a32 (- (ea:determinant m00 m02 m03 m10 m12 m13 0.0 0.0 1.0))a33 (ea:determinant m00 m01 m03 m10 m11 m13 0.0 0.0 1.0)a41 (- (ea:determinant m01 m02 m03 m11 m12 m13 m21 m22 m23))a42 (ea:determinant m00 m02 m03 m10 m12 m13 m20 m22 m23)a43 (- (ea:determinant m00 m01 m03 m10 m11 m13 m20 m21 m23)) ) (list(list a11 a21 a31 a41)(list a12 a22 a32 a42)(list a13 a23 a33 a43)'(0.0 0.0 0.0 1.0) ))