逆矩阵。 « CAD家园

[原创]:[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)  ))

标签: , ,

郑重声明:资讯 【逆矩阵。 « CAD家园】由 发布,版权归原作者及其所在单位,其原创性以及文中陈述文字和内容未经(企业库qiyeku.com)证实,请读者仅作参考,并请自行核实相关内容。若本文有侵犯到您的版权, 请你提供相关证明及申请并与我们联系(qiyeku # qq.com)或【在线投诉】,我们审核后将会尽快处理。
—— 相关资讯 ——