householder实现Ax=b的x非首元素的消除(二)

1.推算

v=x\pm ||x||^2e1(这里假设只取负号,如果想更精确,可以编写sgn符号函数),定义\omega =\nu /||\nu ||_2=(x-\alpha e1)/||x-\alpha e1||_2||\omega ||_2=1,

H=I-2\omega \omega \tau=I-2\nu \nu ^\tau /(||\nu ||_2||\nu ||_2)=I-2\nu \nu ^\tau /(\nu ^\tau \nu )=I-\beta \nu \nu ^\tau ,\beta =2/(\nu ^\tau \nu )

这里为了避免溢出现象,用x/||x||_\infty代替x,

第一步:v=x-||x||_2e1=(v1,v2,...,vn)=(x(1)-||x||_2,x(2),x(3),...,x(n))

第二步:\beta =2/(\nu^\tau \nu )

第三步:H=I-\beta \nu \nu ^\tau

第四步:y=H*x,转变x的非首元素为零

2.代码

实现结果如下:

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值