[编程申请]：求矩阵的广义逆阵VLISP完整源码

chlh_jd

令A是任意mxn矩阵，G是A的广义逆阵，那么G是个nxm的矩阵，则AG满足下列条件：
(1) AGA=A   ;  
(2)  存在H=GA为幂等矩阵，即H^2=H，,rank(H)=rank(A);
(3)  存在F=AG为幂等矩阵，即F^2=F，rank(F)=rank(A);

广义逆阵算法有二：
  算法一：满轶分解法
    步骤1   利用初等行变换将矩阵A化为阶梯型 [ G O]
    步骤2   对单位矩阵 I 执行逆行初等变换，得到 P^-1
        步骤3   利用逆矩阵P^-1的前r列构造矩阵F
        步骤4  书写满轶分解结果 A= F*G
   
     算法二  广义逆矩阵的计算
    步骤1  计算矩阵Amxn的满轶分解 A=F*G
        步骤2  求广义逆矩阵 A^-1=G^T(F^T*A*G^T)^-1*F^T


编码采用VLISP

狂刀

用lisp解矩阵效率是不高的，程序写起来也麻烦，不如用vb或c解后把结果返回cad