补充推导步骤,重写 Matrix Computations 5.1.2 节

本来的内容有点小小的跳跃,补一下跳跃的部分,下次推导时省点时间,备忘

1. 补充后的内容

2. 代码

LaTeX code:

\documentclass{article}
\title{Matrix Computations 5.1.2 time saving revision}
\date{}
\begin{document}
\maketitle
Let $v \in \mathbf{R}^m$ and $v \ne 0$
\begin{equation}\mathbf{P} = \mathbf{I}-\beta vv^T, \beta = \frac{2}{v^Tv}
\end{equation}Then $\mathbf{P}$ is $Householder\,\, reflection$.$y = \mathbf{P}x$ is the reflection of $x$ with supper plane: $\mathbf{span}\{v\}^\bot$\\\\
\begin{equation} \label{eq:eps2}
\mathbf{P}x =
\left(\mathbf{I} - \frac{2vv^T}{v^Tv}
\right)x
=
x - \frac{2v^Tx}{v^Tv}v
\end{equation}
To reflect x onto ${e}_1 = I_m(:,1)$, that means
\begin{equation} \label{eq:eps3}\mathbf{P}x = w*e_1, w \in \mathbf{R}
\end{equation}we can know that $v \in \mathbf{span}\{x, e_1 \}$.
Let\begin{equation} \label{eq:eps4}
v = x + \alpha e_1
\end{equation}We know that $v$ is under the restriction of $x$.Do $x^T * (\ref{eq:eps3})$:
$$x^Tv = x^Tx + \alpha x^Te_1$$
then
\begin{equation}
v^Tx = x^Tx + \alpha x_1
\end{equation}Do $(v^Tv) = [(\ref{eq:eps4})^T*(\ref{eq:eps4})]$:$$v^Tv = (x+\alpha e_1)^T(x+ \alpha e_1) = (x^T + \alpha {e_1}^T)(x+\alpha e_1)$$
then
\begin{equation}
v^Tv = x^Tx + 2\alpha x_1 + \alpha^2
\end{equation}From (\ref{eq:eps2}):
$$\mathbf{P}x = x - \frac{2v^Tx}{v^Tv}v$$
\begin{equation}\label{eq:eps7}
\mathbf{P}x = x - \frac{2v^Tx}{v^Tv}(x + v - x)
=x - \frac{2v^Tx}{v^Tv}x + \frac{2v^Tx}{v^Tv}(v - x)
\end{equation}and from (\ref{eq:eps4})
$$v = x + \alpha e_1$$
thus
\begin{equation} \label{eq:eps8}(v - x) = \alpha e_1
\end{equation}
Take (\ref{eq:eps8}) into (\ref{eq:eps7}):
\begin{equation} \label{eq:eps9}
\mathbf{P}x = (1-\frac{2v^Tx}{v^Tv})x + 2\alpha \frac{v^Tx}{v^Tv}e_1= \left(\frac{\alpha^2 - \parallel x \parallel _2^2}{x^Tx + 2\alpha x_1 + \alpha^2} \right)x - 2\alpha \frac{v^Tx}{v^Tv}e_1
\end{equation}Compare (\ref{eq:eps9}) with (\ref{eq:eps3}), we want:
\begin{equation}\mathbf{P}x = 0*x + we_1; w \in \mathbf{R}
\end{equation}thus
$$\frac{\alpha^2 - \parallel x \parallel _2^2}{x^Tx + 2\alpha x_1 + \alpha^2} = 0$$that means
$$\alpha^2 - \parallel x \parallel _2^2 = 0$$
\begin{equation}
\alpha = \pm \parallel x \parallel_2
\end{equation}Take the $\alpha$ into (\ref{eq:eps4}), we have\begin{equation}
v = x \pm \parallel x \parallel_2e_1
\end{equation}st.$$\mathbf{P}x = \left(\mathbf{I} - 2\frac{vv^T}{v^Tv}\right)x = \mp\parallel x \parallel_2e_1; \,\,\mathbf{where}\,\,v=x\pm \parallel x \parallel_2e_1$$When $v = x \pm \parallel x \parallel_2e_1$, $\mathbf{P}x$ reflects $x$ onto $e_1$.\end{document}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/249300.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

CSRF靶场练习

简述:CSRF漏洞实际很少;条件限制很多;局限性很大;实验仅供参考,熟悉csrf概念和攻击原理即可 Pikachu靶场 CSRF GET 登录用户vince的账户可以看到用户的相关信息; 点击修改个人信息,发现数据包…

[office] excel2010双向条形图制作 #经验分享#微信

excel2010双向条形图制作 本教程为大家介绍一下excel2010中excel2010双向条形图制作方法。 1.选中工作区域 2.点击插入-->图表,选择条形图 3.为美观可将中间竖线可去掉 4.方法是选中竖线,右击-->删除 5.接下来将图例靠上,选中图例,右击-->设置图例格式-->图例选项…

STM32——感应开关盖垃圾桶

STM32——感应开关盖垃圾桶 1.定时器介绍 软件定时 缺点:不精确、占用CPU资源 void Delay500ms() //11.0592MHz {unsigned char i, j, k;_nop_();i 4;j 129;k 119;do{do{while (--k);} while (--j);} while (--i); }定时器工作原理 使用精准的时基&#xff…

【Tomcat与网络8】从源码看Tomcat的层次结构

在前面我们介绍了如何通过源码来启动Tomcat,本文我们就来看一下Tomcat是如何一步步启动的,以及在启动过程中,不同的组件是如何加载的。 一般,我们可以通过 Tomcat 的 /bin 目录下的脚本 startup.sh 来启动 Tomcat,如果…

如何用MapTalks IDE来发布网站?

简介 MapTalks IDE 全称 MapTalks集成设计环境(Integrated Design Environment),是由MapTalks技术团队开发的新一代web地图设计软件。 通过MapTalks IDE,您可以自由的创建二维和三维地图,在其中载入或创建地理数据&a…

计算机语言的发展历史

计算机编程语言的发展,是随着计算机本身硬件发展而发展的。硬件速度越快、体积越小、成本越低,应用到人类社会的场景就会越多,那么所需要的算法就会越复杂,也就要求计算机编程语言越高级。最初重达几十吨但一秒只能运算5000次的EN…

【JavaSE篇】——内部类

目录 🎓内部类 🎈内部类的分类 🚩实例内部类 一.如何实例内部类对象 二.实例内部类中为什么不能有静态成员变量 (用final解决) 三.在实例内部类对象时,如何访问外部类当中相同的成员变量?…

linux中常用的命令

一:tree命令 (码字不易,关注一下吧,w~~w) 以树状形式查看指定目录内容。 tree --树状显示当前目录下的文件信息。 tree 目录 --树状显示指定目录下的文件信息。 注意: tree只能查看目录内容,不能…

基于MongoDB实现聊天记录的存储

一、mongodb简介 1.1 mongodb简介 MongoDB是一个基于分布式文件存储的数据库,使用C语言编写。它旨在为WEB应用提供可扩展的高性能数据存储解决方案。MongoDB介于关系数据库和非关系数据库之间,是非关系数据库当中功能最丰富、最像关系数据库的。 Mong…

C#用正则表达式判断字符串是否纯数字vs用Char.IsDigit 方法遍历字符数组是否纯数字

目录 一、使用的方法 1.正则表达式 2.Char.IsDigit 方法 二、源码 1.源代码 2.生成效果 一、使用的方法 1.正则表达式 在程序运行过程中,经常需要用户输入数字信息,如输入员工年龄、工资等。使用正则表达式Regex类的IsMatch方法,可以有…

【ASP.NET Core 基础知识】--Web API--创建和配置Web API(一)

一、简介 Web API(Web Application Programming Interface)的重要性在于其在现代软件开发中扮演着关键的角色。以下是一些关于Web API重要性的方面: 跨平台交互: Web API允许不同平台、不同技术栈的应用程序进行通信。无论是Web…

如何本地搭建Emby影音管理服务并结合内网穿透实现远程访问本地影音库

文章目录 1.前言2. Emby网站搭建2.1. Emby下载和安装2.2 Emby网页测试 3. 本地网页发布3.1 注册并安装cpolar内网穿透3.2 Cpolar云端设置3.3 Cpolar内网穿透本地设置 4.公网访问测试5.结语 1.前言 在现代五花八门的网络应用场景中,观看视频绝对是主力应用场景之一&…

负载均衡下的webshell连接

一、环境配置 1.在Ubuntu上配置docker环境 我们选择用Xshell来将环境资源上传到Ubuntu虚拟机上(比较简单) 我们选择在root模式下进行环境配置,先将资源文件复制到root下(如果你一开始就传输到root下就不用理会这个) …

微分几何——梅向明第四版学习笔记(一) 向量函数和曲线论

目录 引出向量函数曲线论简单曲线定义曲线的向量参数表示 曲线的切线【重要】曲线的法面【重要】曲线的自然参数表示 空间曲线曲线的密切平面空间曲线的基本三棱形【重要】单位切向量主法向量副法向量Frenet标架螺旋线的案例 曲线的曲率和曲率半径曲率的几何意义 曲线的挠率挠率…

顺序表与链表,栈与队列

名词辨析:指针 1.什么是指针,想必大家都不陌生,但是,在这部分的知识中,包含着一类特殊的指针,表面上它只是单个的数字,但它其实代表了作为栈或者队列载体的数组的下标,在实际题目中…

Golang语言异常机制解析:错误策略与优雅处理

前些天发现了一个巨牛的人工智能学习网站,通俗易懂,风趣幽默,忍不住分享一下给大家。点击跳转到网站https://www.captainbed.cn/kitie。 前言 作为开发者来说,我们没办法保证程序在运行过程中永远不会出现异常,对于异常…

Java多线程编程中的异常处理策略

第1章:引言 大家好,我是小黑,咱们今天聊聊异常处理。想必大家在写代码的时候都遇到过各种各样的异常吧?有时候,一个小小的异常如果处理不当,就可能导致整个程序崩溃。特别是在多线程环境下,异常…

科技云报道:云原生PaaS,如何让金融业数字化开出“繁花”?

科技云报道原创。 在中国金融业数字化转型的历史长卷中,过去十年无疑是一部磅礴的史诗。 2017年,南京银行第一次将传统线下金融业务搬到了线上。那一年,它的互联网金融信贷业务实现了过去10年的业务总额。 2021年,富滇银行通过…

幻兽帕鲁服务器游戏怎么升级版本?

幻兽帕鲁服务器游戏怎么升级版本?自建幻兽帕鲁服务器进入Palworld游戏提示“您正尝试加入的比赛正在运行不兼容的游戏版本,请尝试升级游戏版本”什么原因?这是由于你的客户端和幻兽帕鲁服务器版本不匹配,如何解决?更新…

配置IPv6静态路由

1、静态路由简介 静态路由是一种需要管理员手工配置的特殊路由。 静态路由在不同网络环境中有不同的目的: 当网络结构比较简单时,只需配置静态路由就可以使网络正常工作。 在复杂网络环境中,配置静态路由可以改进网络的性能,并…