[LISP函数-矩阵]：求矩阵的逆矩阵

eachy

2. 
   
3. ;;|;
   
4. ;;获取wcs－》ucs的矩阵 3d转换
   
5. ;;求矩阵的逆矩阵;
   
6. (defun xd-matrix_inverse (matrix /        m0     m1     m2     m3
   
7.                           m00         m01        m02    m03    m10    m11
   
8.                           m12         m13        m20    m21    m22    m23
   
9.                           _a         ea:determinant              a11    a12
   
10.                           a21         a22        a33    a34    a41    a42
    
11.                           a43         a13        a23    a31    a32
    
12.                          )
    
13.   ;;行列式展开，计算代数余子式
    
14.   (defun ea:determinant        (a00 a01 a02 a10 a11 a12 a20 a21 a22 /)
    
15.     (- (+ (* a00 a11 a22)
    
16.           (* a01 a12 a20)
    
17.           (* a02 a21 a10)
    
18.        )
    
19.        (* a02 a11 a20)
    
20.        (* a01 a10 a22)
    
21.        (* a00 a21 a12)
    
22.     )
    
23.   )
    
24.   (setq        m0 (car matrix)
    
25.         m1 (cadr matrix)
    
26.         m2 (nth 2 matrix)
    
27.         m3 (last matrix)
    
28.   )
    
29.   (setq        m00 (car m0)
    
30.         m01 (cadr m0)
    
31.         m02 (nth 2 m0)
    
32.         m03 (last m0)
    
33.   )
    
34.   (setq        m10 (car m1)
    
35.         m11 (cadr m1)
    
36.         m12 (nth 2 m1)
    
37.         m13 (last m1)
    
38.   )
    
39.   (setq        m20 (car m2)
    
40.         m21 (cadr m2)
    
41.         m22 (nth 2 m2)
    
42.         m23 (last m2)
    
43.   )
    
44.   (setq        a11 (ea:determinant m11 m12 m13 m21 m22 m23 0.0 0.0 1.0)
    
45.         a12 (- (ea:determinant m10 m12 m13 m20 m22 m23 0.0 0.0 1.0))
    
46.         a13 (ea:determinant m10 m11 m13 m20 m21 m23 0.0 0.0 1.0)
    
47.         a21 (- (ea:determinant m01 m02 m03 m21 m22 m23 0.0 0.0 1.0))
    
48.         a22 (ea:determinant m00 m02 m03 m20 m22 m23 0.0 0.0 1.0)
    
49.         a23 (- (ea:determinant m00 m01 m03 m20 m21 m23 0.0 0.0 1.0))
    
50.         a31 (ea:determinant m01 m02 m03 m11 m12 m13 0.0 0.0 1.0)
    
51.         a32 (- (ea:determinant m00 m02 m03 m10 m12 m13 0.0 0.0 1.0))
    
52.         a33 (ea:determinant m00 m01 m03 m10 m11 m13 0.0 0.0 1.0)
    
53.         a41 (- (ea:determinant m01 m02 m03 m11 m12 m13 m21 m22 m23))
    
54.         a42 (ea:determinant m00 m02 m03 m10 m12 m13 m20 m22 m23)
    
55.         a43 (- (ea:determinant m00 m01 m03 m10 m11 m13 m20 m21 m23))
    
56.   )
    
57.   (list        (list a11 a21 a31 a41)
    
58.         (list a12 a22 a32 a42)
    
59.         (list a13 a23 a33 a43)
    
60.         '(0.0 0.0 0.0 1.0)
    
61.   )
    
62. )

<br /> ;;|;<br /> ;;获取wcs－》ucs的矩阵 3d转换<br /> ;;求矩阵的逆矩阵;<br /> (defun xd-matrix_inverse (matrix /        m0     m1     m2     m3<br />                           m00         m01        m02    m03    m10    m11<br />                           m12         m13        m20    m21    m22    m23<br />                           _a         ea:determinant              a11    a12<br />                           a21         a22        a33    a34    a41    a42<br />                           a43         a13        a23    a31    a32<br />                          )<br />   ;;行列式展开，计算代数余子式<br />   (defun ea:determinant        (a00 a01 a02 a10 a11 a12 a20 a21 a22 /)<br />     (- (+ (* a00 a11 a22)<br />           (* a01 a12 a20)<br />           (* a02 a21 a10)<br />        )<br />        (* a02 a11 a20)<br />        (* a01 a10 a22)<br />        (* a00 a21 a12)<br />     )<br />   )<br />   (setq        m0 (car matrix)<br />         m1 (cadr matrix)<br />         m2 (nth 2 matrix)<br />         m3 (last matrix)<br />   )<br />   (setq        m00 (car m0)<br />         m01 (cadr m0)<br />         m02 (nth 2 m0)<br />         m03 (last m0)<br />   )<br />   (setq        m10 (car m1)<br />         m11 (cadr m1)<br />         m12 (nth 2 m1)<br />         m13 (last m1)<br />   )<br />   (setq        m20 (car m2)<br />         m21 (cadr m2)<br />         m22 (nth 2 m2)<br />         m23 (last m2)<br />   )<br />   (setq        a11 (ea:determinant m11 m12 m13 m21 m22 m23 0.0 0.0 1.0)<br />         a12 (- (ea:determinant m10 m12 m13 m20 m22 m23 0.0 0.0 1.0))<br />         a13 (ea:determinant m10 m11 m13 m20 m21 m23 0.0 0.0 1.0)<br />         a21 (- (ea:determinant m01 m02 m03 m21 m22 m23 0.0 0.0 1.0))<br />         a22 (ea:determinant m00 m02 m03 m20 m22 m23 0.0 0.0 1.0)<br />         a23 (- (ea:determinant m00 m01 m03 m20 m21 m23 0.0 0.0 1.0))<br />         a31 (ea:determinant m01 m02 m03 m11 m12 m13 0.0 0.0 1.0)<br />         a32 (- (ea:determinant m00 m02 m03 m10 m12 m13 0.0 0.0 1.0))<br />         a33 (ea:determinant m00 m01 m03 m10 m11 m13 0.0 0.0 1.0)<br />         a41 (- (ea:determinant m01 m02 m03 m11 m12 m13 m21 m22 m23))<br />         a42 (ea:determinant m00 m02 m03 m10 m12 m13 m20 m22 m23)<br />         a43 (- (ea:determinant m00 m01 m03 m10 m11 m13 m20 m21 m23))<br />   )<br />   (list        (list a11 a21 a31 a41)<br />         (list a12 a22 a32 a42)<br />         (list a13 a23 a33 a43)<br />         '(0.0 0.0 0.0 1.0)<br />   )<br /> )

aeo

看不懂.
我是把移动，比例，旋转分别求出来再组合的.


(setq        m0 (car matrix)
        m1 (cadr matrix)
        m2 (nth 2 matrix)
        m3 (last matrix)
  )
  (setq        m00 (car m0)
        m01 (cadr m0)
        m02 (nth 2 m0)
        m03 (last m0)
  )
  (setq        m10 (car m1)
        m11 (cadr m1)
        m12 (nth 2 m1)
        m13 (last m1)
  )
  (setq        m20 (car m2)
        m21 (cadr m2)
        m22 (nth 2 m2)
        m23 (last m2)
  )

-->

(mapcar 'set '(m00 m01 m02 m03 m10 m11 m12 m13 m20 m21 m22 m23)(apply 'append matrix))
是这个意思吧

eachy

最初由 aeo 发布
[B]看不懂.
我是把移动，比例，旋转分别求出来再组合的.


(setq        m0 (car matrix)
        m1 (cadr matrix)
        m2 (nth 2 matrix)
        m3 (last matrix)
  )
  (setq        m00 (car m0)
        m01 (cadr m0)
        m02 (nth 2 m0)
        m0... [/B]

这是很早写的一个函数，以后就没有再改过。就是按照行列式展开写的，也许用Lisp有更简洁的写法。

cheungyh

任意阶的逆矩阵如何求?