欢迎光临散文网 会员登陆 & 注册

Householder矩阵计算公式的简单推导

2023-03-16 21:53 作者:散佚的魔导书  | 我要投稿

Householder矩阵计算公式的简单推导的记录,如果我忘记了就可以再看一下这个记录...

如果对任何人有帮助就再好不过了,如果其中有错误,请一定要告诉我,谢谢啦! 

Householder矩阵的几何意义是将向量按照一个平面镜像。

如果有一个Householder矩阵 H 和一个向量 A,那么 使用矩阵变换此向量,也就是

A_%7Bmirrored%7D%3DHA

这样可以得到镜像后的向量 A_%7Bmirrored%7D

如下图所示

通过镜像A得到A_mirrored,其中B是平面的法向量,未归一化


看书上提到Householder矩阵的计算公式是

H%3DI-2ww%5ET

其中w是归一化的平面法向量,若按照上图的说明,就有

w%3D%5Cfrac%7BB%7D%7B%7CB%7C%7D%20

下面是对此公式的简单推导:

如下图所示,可以看出向量 A 和它的镜像之间相差2倍的投影长度,此投影长度是向量 A 对向量 B 的投影的长度。

所以有

A_%7Bmirrored%7D%3DA-2A_%7Bprojected%7D

带入向量对向量的投影公式 A_%7Bprojected%7D%3D%5Cfrac%7BA%5Ccdot%20B%7D%7BB%5Ccdot%20B%7D%20B ,可得

A_%7Bmirrored%7D%3DA-2%5Cfrac%7BA%5Ccdot%20B%7D%7BB%5Ccdot%20B%7DB%20

向量也可以被看作是长或宽为1的矩阵,向量点积是也可化成矩阵乘积的形式,也就是

A%5Ccdot%20B%3DA%5ETB%3DB%5ETA

B%5Ccdot%20B%3DB%5ETB

带入上式可得

A_%7Bmirrored%7D%3DA-2%5Cfrac%7BB%5ETA%7D%7BB%5ETB%7DB%20

由于 B%5ETA 和 B%5ETB 都代表常数,可整理得

A_%7Bmirrored%7D%3DA-2%5Cfrac%7B1%7D%7BB%5ETB%7D%20B(B%5ETA)

虽然矩阵乘法不满足交换律,但是满足结合律,所以有

A_%7Bmirrored%7D%3DA-%5Cfrac%7B2%7D%7BB%5ETB%7D%20(BB%5ET)A

A_%7Bmirrored%7D%3D(I-%5Cfrac%7B2%7D%7BB%5ETB%7D%20BB%5ET)A

其中,B%5ETB%3DB%5Ccdot%20B%3D%7CB%7C%5E2

至此,若不要求平面法向量归一化,上述公式已经推导完毕了。

如果要得到关于归一化法向量的公式,也就是有  w%3D%5Cfrac%7BB%7D%7B%7CB%7C%7D%20,其中向量 w 是归一化的法向量,有 ww%5ET%3D%5Cfrac%7BB%7D%7B%7CB%7C%7D%5Cfrac%7BB%5ET%7D%7B%7CB%7C%7D%3D%5Cfrac%7BBB%5ET%7D%7BB%5ETB%7D%20%20%20

带入上式可以得到,

A_%7Bmirrored%7D%3D(I-%5Cfrac%7B2%7D%7BB%5ETB%7D%20BB%5ET)A

A_%7Bmirrored%7D%3D(I-2%5Cfrac%7BBB%5ET%7D%7BB%5ETB%7D%20)A

A_%7Bmirrored%7D%3D(I-2ww%5ET%20)A

其中 I-2ww%5ET 就是镜像变换矩阵Householder矩阵。

Householder矩阵计算公式的简单推导的评论 (共 条)

分享到微博请遵守国家法律