﻿功能:采用初等行或列变换来求解逆

格式:Inv4(a)

说明:a是一个方阵,本函数的特点是经过初等变换后,尽量使结果接近矩阵的逆.如果方阵本身非奇异,则直接求得其逆

例子:

a =
[  1.00000000000000   2.00000000000000   3.00000000000000
   4.00000000000000   5.00000000000000   7.00000000000000
   2.00000000000000   4.00000000000000   6.00000000000000   ]

b=inv4(a)//回车后得到如下结果,注意,a是一个行列式为0的矩阵,看第一行与最后一行的关系
b =
[  0.00000000000000   0.66666666666666   -0.8333333333333
   0.00000000000000   -0.3333333333333   0.66666666666666
   1.00000000000000   0.00000000000000   -0.5000000000000   ]

mul(b,a)//我们检验下结果,你会看到,结果是尽量为单位矩阵
ans =
[  1.00000000000000   0.00000000000000   -0.3333333333333
   0.00000000000000   1.00000000000000   1.66666666666667
   0.00000000000000   0.00000000000000   0.00000000000000   ]