豪斯霍德(Householder Reduction)转换成三对角形式(python,数值积分)

第三十一篇 豪斯霍德转换成三对角形式

三对角形式

三对角矩阵就是对角线、邻近对角线的上下次对角线上有元素，其他位置均为0的矩阵。

雅可比主对角化求特征值要求转化矩阵应该具有下面的性质

豪斯霍尔德技术要求选择p为

其中{w}是一个归一化的列向量，使其欧几里德范数等于单位，因此

例如，让

按照上面的乘积

然后可以得到

它具有下面的性质

为消去[A]在对角线外第一行的项，向量{w}应被取为

因此，假设[A]为3 × 3，对第一行的变换矩阵为

当乘积[P][A][P]被计算时，结果矩阵的第一行包含以下三项

让

上面方程的项可以写成

通过向量的点乘为点位1，得到

因此通过方程(1)后面两项的平方，再代换上面方程得到

然后可以写到

再取方程（2）的平方，然后把上面两个方程替换进去得到，

然后把向量写成用向量表示的形式

其中

于是，转化矩阵为

对于r的确定，根据之前的等式，应该选择r的符号与a12是相反的。
所以对于通常的第i行，{v}采取下面形式的向量