【原创】《矩阵的史诗级玩法》连载十一：逆矩阵的代码实现及其在45度地图中的应用

本文链接：https://blog.csdn.net/iloveas2014/article/details/82944815

从搜索引擎过来这篇文章的朋友可能会有点失望，因为我没在标题上说明是多少阶矩阵的代码。不得不说，固定阶数，并且还只是3阶的求逆实在是太简单了，上篇说初中生都能看懂。而任意高阶数则需要借助诸如克拉默法则一类的定理进行实现，并且可能还得尝试用高斯消元法进行优化。然而这些我都没去做。

但不管如何，既然缘分让我们聚在一起，那新来的朋友不妨多留几分钟看下我的总目录，没准能找到你们想要的东西。

下面开始正题时间。

前面我把矩阵的缩放平移旋转变换写到了Matrix类以外的工具类MatrixUtil中，是考虑到这些变换是功能性的，属于业务层面，而非纯数学上的概念，至于乘法它是矩阵本身所具备的一种运算，它跟业务无关所以我就写到了Matrix类的自身。逆矩阵跟乘法一样，也是纯数学的东西，因此我会把逆矩阵的计算追加到Matrix中，代码如下。

/**
 * 进行求逆运算，会修改当前矩阵的数据
 *
 **/
function invert() {
	var a = this.a;
	var b = this.b;
	var c = this.c;
	var d = this.d;
	var tx = this.tx;
	var ty = this.ty;
	var D = a * d - b * c; //D是行列式的符号，可自行查阅线性代数教材
	this.a = d / D;
	this.b = - b / D;
	this.c = - c / D;
	this.d = a / D;
	this.tx = (c * ty - d * tx) / D;
	this.ty = (b * tx - a * ty) / D;
}
this.invert = invert;

然后我们回到45度地图上，把基于业务的矩阵逆变换改为调用纯数学的invert方法。

为便于对照，我先把之前的代码放到这里来。

var matrix = new Matrix();
MatrixUtil.scale(matrix, Math.sqrt(2) / 2, Math.sqrt(2) / 2);	
MatrixUtil.rotate(matrix, Math.PI / 4);
MatrixUtil.scale(matrix, 2, 1);
MatrixUtil.translate(matrix, 400, 0);
	
var matrixInvert = new Matrix();
MatrixUtil.translate(matrixInvert, -400, 0);
MatrixUtil.scale(matrixInvert, 0.5, 1);
MatrixUtil.rotate(matrixInvert, -Math.PI / 4);
MatrixUtil.scale(matrixInvert, Math.sqrt(2), Math.sqrt(2));

改动的部分只有matrixInert：

var matrix = new Matrix();
MatrixUtil.scale(matrix, Math.sqrt(2) / 2, Math.sqrt(2) / 2);	
MatrixUtil.rotate(matrix, Math.PI / 4);
MatrixUtil.scale(matrix, 2, 1);
MatrixUtil.translate(matrix, 400, 0);

var matrixInvert = matrix.clone();
matrixInvert.invert();

这样的话，以后不管原始的matrix做何种业务变换，步骤有多少步，求逆的时候一个invert方法就可以解决了。

当然了，写好还是要验证一下，截图似乎无法说明效果的正确与否，那我也形式性地给出一张，大家自行尝试吧。