householder实现Ax=b的x非首元素的消除（二）

最新推荐文章于 2024-08-30 11:28:00 发布
最新推荐文章于 2024-08-30 11:28:00 发布
阅读量31 收藏
点赞数
文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_74981381/article/details/133321899
版权
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
    评论
视觉slam14讲学习笔记（持续更新中）
DavidHan008
05-07 5万+
写在最前面： SLAM特指：特指搭载传感器的主体，在没有环境先验的信息情况下，在运动过程中建立环境模型，通过估计自己的运动。 SLAM的目的是解决两个问题：1、定位 2、地图构建 也就是说，要一边估计出传感器自身的位置，一边要建立周围环境的模型 最终的目标：实时地，在没有先验知识的情况下进行定位和地图重建。 当相机作为传感器的时候，我们要做的就是根据一张张连续运动的图像，从中估计出
期中考试题目1
08-04
2. Givens变换用于消除矩阵的非零元素，逐步得到QR分解。在 `B=QR` 中，`Q` 是正交矩阵，`R` 是上三角矩阵。 3. 使用QR分解进行一次迭代，得到 `A_2`，并证明 `A_2` 与 `B` 相似，即存在一个可逆矩阵 `P` 使得 `A_2 ...
数值分析常见基本算法及MATLAB代码总结
热门推荐
xiaoye_dlut的博客
12-30 2万+
刚考完研究生的数值分析课，打算整理一下平时的上机代码，做一个汇总。 数值分析大致八部分内容—— 解线性方程组的直接法：Gauss消去法与矩阵三角分解法（Doolittle分解法相比Crout分解法更常用）及其选择列主元的改进方法、Doolittle分解法的延伸（实对称正定矩阵利用Cholesky分解得到的平方根法、三对角矩阵作为线性方程组系数矩阵的追赶法） 解线性方程组的迭代法：Jacobi迭代法、Gauss-Seidel迭代法（利用前者每次迭代已得到的最新分量加速）、逐次超松弛（SOR，Succes
<Numerical Analysis>(by Timothy Sauer) Notes
weixin_34220834的博客
02-03 489
2ed, by Timothy Sauer DEFINITION 1.3A solution is correct within p decimal places if the error is less than 0.5 × 10$^{−p}$ .-P29Bisection Method的优点是计算次数（step）是确定的（interval<精度）。后面介绍的算法的interva...
Square Root SAM论文原理
weixin_41469272的博客
07-07 971
文章目录Square Root SAM论文原理核心原理SLAM问题的3种表示贝叶斯网络因子图（Factor graph）马尔科夫随机场(Markov Random Field, MRF)SLAM最小二乘问题&线性化因式分解 factorization矩阵与图(Matrices ⇔ Graphs)因式分解&变量消元(Factorization ⇔ Variable Elimination )Cholesky分解（或LDL分解）变量消除在Cholesky分解中的作用步骤QR分解变量消除在QR分解
视觉SLAM ch3—三维空间的刚体运动
Johaden的博客
08-10 719
配套高翔老师以及其他视频资源的视觉SLAM笔记。 里面的视频可以在指示下观看，四元数的两个视频很好，建议全部观看！！
《视觉 SLAM 十四讲》V2 ——第 3 讲 【编程环境+旋转表示的四种方式】
weixin_46034116的博客
09-28 231
待做： https://zhuanlan.zhihu.com/p/388386417 关于本笔记的说明： 最好跟着 原书 整理个人笔记，他人笔记仅适合参考部分内容。 ———————— B站链接 高翔博客链接 百度网盘链接：https://pan.baidu.com/s/1VsrueNrdqmzTvh-IlFBr9Q 提取码：vfhe github源码链接V2 《视觉SLAM十四讲》V2 2019年 动态定位和环境建模 SLAM 同时定位与地图构建 SLAM Simultaneous Localiz
SLAM面经（持续更新中，欢迎关注一起学习......）
qq_43200940的博客
07-03 2349
SLAM面经
slam笔记
danmeng8068的博客
11-21 3460
写在最前面：  SLAM特指：特指搭载传感器的主体，在没有环境先验的信息情况下，在运动过程中建立环境模型，通过估计自己的运动。  SLAM的目的是解决两个问题：1、定位 2、地图构建  也就是说，要一边估计出传感器自身的位置，一边要建立周围环境的模型  最终的目标：实时地，在没有先验知识的情况下进行定位和地图重建。  当相机作为传感器的时候，我们要做的就是根据一张张连续运动的图像
《数值分析》-3-特征值与特征矩阵
chococolate's coding space
03-02 2320
数值分析学习笔记
QR分解_qr分解_
09-28
Givens旋转通过一系列小角度的旋转矩阵来逐步消除非主对角线元素，而Householder反射则通过反射操作一次性消除一整列的大部分元素，相对更高效。在实际应用中，Householder反射更为常见。 在求解线性最小二乘问题时...
李小山-第2章作业1
08-03
其核心是通过行交换和行标量乘积来消除非主元列的非零元素，以达到简化系统的目的。 - **QR分解**：将矩阵A分解为一个正交矩阵Q和一个上三角矩阵R，使得A=QR。这种分解在解决线性方程组时非常有用，因为它允许我们...
Leetcode42接雨水（单调栈）
AGNING的博客
08-28 359
维护一个栈，满足从栈底到栈顶元素大小递减。如果当前元素高于栈顶，且栈内有至少2个元素，则形成了一个凹槽，及一个图中横着的红色矩形，将它的面积累加即可，第一次写难在维护单调栈中计算面积和。求出前缀最大和后缀最大，用两者较小值减去当前高度，累加即可，这个思路容易想到，这里不赘述。
趣味算法------试用 6 和 9 组成的最大数字
markingyi的博客
08-26 554
给出一个数字，反转其中一位使数字最大化，递归求解。
递归、搜索与回溯算法
无聊_人的博客
08-29 3117
所谓的递归其实就是函数自己调用自己的情况。
代码随想录算法训练营第62天|Floyd 算法精讲、A * 算法精讲 （A star算法）
Yinems的博客
08-28 1237
分两种情况，一种是节点 i 到节点 j 的最短路径经过节点 k，grid[i][j][k] = grid[i][k][k-1] + grid[k][j][k-1]；另一种情况，节点 i 到节点 j 的最短路径不经过节点 k，grid[i][j][k] = grip[i][j][k-1]；刚开始初始化 k 是不确定的，如果赋值为1，那么不能知道节点 i 到节点 j 经过1个节点的最短距离，所以只能把 k 赋值为0，本题中节点0是无意义的，节点是从1到n。本题求的是最小值，所以其余值初始化为最大值。
趣味算法------猴子吃桃(循环，递归双重解法)
markingyi的博客
08-26 454
猴子吃桃问题，已知结果求初始值，分别使用循环和递归两种方式解决。
精选算法编程题
最新发布
A18937260063的博客
08-30 296
数组编程
householder变换的实现步骤
05-13
下面是实现Householder变换的步骤： 1. 选择一个向量$v$，使得$v$与要进行变换的向量$x$在同一平面上，并且$v$的长度为1。可以选择$v$为$x$的一个单位正交化的垂直向量，即$v = \frac{x - ||x||e_1}{||x - ||x||e_1...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值