![](https://csdnimg.cn/release/download_crawler_static/86385193/bg9.jpg)
使得 Hn 的第 i 列 (i+1) 行到第 n 行元素为 0。
因此可以求得
Hk ·Hk-1·。。。。H2·H1·Hilber(n) · x =Hk ·Hk-1·。。。。H2·H1·b
其中
Hk ·Hk-1·。。。。H2·H1·Hilber(n)
为上三角矩阵。特别地,对于
由于 Householder 矩阵具有性质
是 n+r 阶 Householder 矩阵,所以对于
第
i 列求得的(n-i)阶矩阵 Hi 可以利用这种方式扩充为 n 阶矩阵,从而进行矩阵
乘法。
Householder 变换后,x 的左侧为上三角矩阵。因此可以使用回带函数求解得到 x
2.算法实现
for i in range(n):
#求二范数
k=norm(a[i:n,i],type)