Gram-Schmidt正交化方法(过目不忘)
标准正交基的Gram-Schmidt求解方法(过目不忘)
- 一、定理
- 二、几何意义
- 1、两向量垂直
- 2、定义垂直操作
- 3、α2\alpha_2α2和α3\alpha_3α3同时关于β1\beta_1β1垂直变换
- 4、β3′\beta_3'β3′关于β2\beta_2β2进行垂直变换
一、定理
如果α1,…,αn\alpha_1,\dots,\alpha_nα1,…,αn是n维Euclid空间V的一个基(基肯定是线性无关的),若令
β1=α1\beta_1 = \alpha_1 β1=α1β2=α2−<α2⋅β1><β1⋅β1>β1\beta_2 = \alpha_2 - \frac{<\alpha_2·\beta_1>}{<\beta_1·\beta_1>}\beta_1 β2=α2−<β1⋅β1><α2⋅β1>β1⋮\vdots ⋮βn=α2−<α2⋅β1><β1⋅β1>β1−,⋯−<αn⋅βn−1><βn−1⋅βn−1>βn−1,\beta_n = \alpha_2 - \frac{<\alpha_2·\beta_1>}{<\beta_1·\beta_1>}\beta_1 - ,\dots - \frac{<\alpha_n·\beta_{n-1}>}{<\beta_{n-1}·\beta_{n-1}>}\beta_{n-1}, βn=α2−<β1⋅β1><α2⋅β1>β1−,⋯−<βn−1⋅βn−1><αn⋅βn−1>βn−1,
则β1,…,βn\beta_1,\dots,\beta_nβ1,…,βn就是V的一个正交基.若再令
ei=βi∣∣βi∣∣,i=1,2,…,n,e_i = \frac{\beta_i}{||\beta_i||},\quad i = 1,2,\dots,n, ei=∣∣βi∣∣βi,i=1,2,…,n,就得到了V的标准正交基e1,e2,…,ene_1,e_2,\dots,e_ne1,e2,…,en.
\quad虽然上面的公式已经十分有规律了,时间长了还难免会把小角标记错。书上也有推导的原理(非常简单),但是数学公式还是太抽象了;那我们可以想一下3维Euclid空间中该公式的意义,然后再将其推广更高纬度的空间中去。
二、几何意义
\quad简单推导一下。
1、两向量垂直
如果α1,α2,α3\alpha_1,\alpha_2,\alpha_3α1,α2,α3是3维Euclid空间VVV的一个基,β1,β2,β3\beta_1,\beta_2,\beta_3β1,β2,β3是待求的正交基,设
β1=α1,\beta_1 = \alpha_1, β1=α1,β2=α2−kβ1,由于正交求得k=<α2⋅β1><β1⋅β1>.\beta_2 = \alpha_2 - k\beta_1,由于正交求得k = \frac{<\alpha_2·\beta_1>}{<\beta_1·\beta_1>}. β2=α2−kβ1,由于正交求得k=<β1⋅β1><α2⋅β1>.
我们考虑一下β2\beta_2β2是怎么得到的:在α2\alpha_2α2的基础上减了kβ1k\beta_1kβ1,这个值是α2\alpha_2α2在α1\alpha_1α1上的投影长度。这个操作之后确实就垂直了,如下图所示:
2、定义垂直操作
我们定义下面操作
βi=αi−<αi⋅βj><βj⋅βj>βj\beta_i = \alpha_i -\frac{<\alpha_i·\beta_j>}{<\beta_j·\beta_j>}\beta_j βi=αi−<βj⋅βj><αi⋅βj>βj
为αi\alpha_iαi关于βj\beta_jβj的垂直变换。
3、α2\alpha_2α2和α3\alpha_3α3同时关于β1\beta_1β1垂直变换
在α2\alpha2α2关于β1\beta_1β1垂直变换得到β2\beta_2β2后,我们对α3\alpha_3α3也进行关于β1\beta_1β1垂直变换得到β3′\beta_3'β3′,有
β3′=α3−<α3⋅β1><β1⋅β1>β1\beta_3' = \alpha_3 -\frac{<\alpha_3·\beta_1>}{<\beta_1·\beta_1>}\beta_1 β3′=α3−<β1⋅β1><α3⋅β1>β1
此时假设空间直角坐标系的x,z轴与β1、β2\beta_1、\beta_2β1、β2重合;β3′\beta_3'β3′在平面xOyxOyxOy,垂直β1\beta_1β1,不垂直β2\beta_2β2。那么在三维空间中绘制出来的效果如下:
4、β3′\beta_3'β3′关于β2\beta_2β2进行垂直变换
明显我们离完成正交基的求解只有一步了:只需将β3′\beta_3'β3′垂直变换后得到β3\beta_3β3使其垂直β2\beta_2β2。
β3=β3′−<β3′⋅β2><β2⋅β2>β2,\beta_3 = \beta_3' - \frac{<\beta_3'·\beta_2>}{<\beta_2·\beta_2>}\beta_2, β3=β3′−<β2⋅β2><β3′⋅β2>β2,由于<β3′⋅β2>=<(α3−<α3⋅β1><β1⋅β1>β1)⋅β2>=<α3⋅β2>,由于<\beta_3'·\beta_2> = <(\alpha_3 -\frac{<\alpha_3·\beta_1>}{<\beta_1·\beta_1>}\beta_1)·\beta_2> = <\alpha_3·\beta_2>, 由于<β3′⋅β2>=<(α3−<β1⋅β1><α3⋅β1>β1)⋅β2>=<α3⋅β2>,β3=α3−<α3⋅β1><β1⋅β1>β1−<α3⋅β2><β2⋅β2>β2\beta_3= \alpha_3 - \frac{<\alpha_3·\beta_1>}{<\beta_1·\beta_1>}\beta_1 - \frac{<\alpha_3·\beta_2>}{<\beta_2·\beta_2>}\beta_2 β3=α3−<β1⋅β1><α3⋅β1>β1−<β2⋅β2><α3⋅β2>β2\quad至此我们已经得到了三维空间中的正交基β1,β2,β3\beta_1,\beta_2,\beta_3β1,β2,β3了,观察β3的结构\beta_3的结构β3的结构每一项都相当于一次关于已知的正交基的垂直转动。想像空间中三个起点连接的向量,在空间中随机分布,我们要用手把它们掰成相互垂直关系。
\quad使用相同的方式我们就能够得到n维空间中的正交向基了,单位化之后得到标准正交基。
Gram-Schmidt正交化方法(过目不忘)相关推荐
- 通俗理解“Schmidt正交化”和“正交矩阵” 此博文包含图片 (2015-05-19 09:50:47) 施密特正交化在空间上是不断建立垂直于原次维空间的新向量的过程。 如图β2垂直于β1(1维)
通俗理解"Schmidt正交化"和"正交矩阵" 施密特正交化在空间上是不断建立垂直于原次维空间的新向量的过程. 如图β2垂直于β1(1维)构建新2维,β3垂直于 ...
- 线性空间标准正交基的构造——Schmidt正交化
证明: 两两正交的一组向量必然线性无关 线性空间的标准正交基定义 (定义摘自<矩阵论>,程云鹏,1999年第2版,p88) 下面的定理证明了欧式空间一定存在标准正交基,其证明过程可以看作把 ...
- 线性代数学习笔记——第七十讲——格拉姆—施密特(Gram-Schmidt)正交化方法
1. 格拉姆-施密特(Gram-Schmidt)正交化方法详解 2. 格拉姆-施密特(Gram-Schmidt)正交化示例
- 【线性代数】施密特正交化方法——Python实现
文章目录 思想 Python代码 思想 施密特正交化方法: 将n维子空间中的任意一组基向量变换成标准正交向量. 假设有两个向量a⃗\vec{a}a和b⃗\vec{b}b,若要使两向量正交,则a⃗\ve ...
- schmidt正交化
设r1,r2,r3线性无关 第一步,正交化 第二步,单位化 证明施密特正交化
- 格拉姆--施密特(Gram-Schmidt)正交化方法笔记
前段时间,身边的小伙伴问过我关于线性代数的一些知识,其中包含了特征向量的正交化的内容.遥想当初在学习线性代数的时候,只知道施密特正交化可以用来解决这个问题,但是老师猝不及防地甩过来一大段公式(可见下式 ...
- 正交矩阵和一种简单基础的Gram-Schmidt斯密特正交化方法
定义 Q为方阵且: QTQ=I,QT=Q−1Q^{T} Q=I, Q^{T}=Q^{-1} QTQ=I,QT=Q−1 则Q为正交矩阵. 组成Q的向量有: qiTqj={0i≠j(ortho) 1i=j ...
- 矩阵论复习-过渡矩阵、生成子空间、表示矩阵、度量矩阵、酉空间、内积
一小部分矩阵论的整理复习,这个由于公式输入的太麻烦了,所以就弄了一点.后面直接看着书复习的. 矩阵论复习 线性空间 基与维数 基变换公式.过渡矩阵.坐标变换公式 线性子空间 齐次方程组的解空间 特征子 ...
- 矩阵分析——QR分解
Gram-Schmidt正交化 在提到矩阵的QR分解前,必须要提到Gram–Schmidt方法,理论上QR分解是由Gram–Schmidt正交化推出来的.那么Gram–Schmidt正交化究竟是什么. ...
最新文章
- 深入浅出统计学(十三)假设检验
- matlab uitree update,elementUI tree 懒加载 更新节点
- .NET仓储模式高级用例
- vs2005制作安装包(自动安装.net framework 2.0)(小寿原创)
- python 中 feedparser的简单用法
- [转]VS2010中水晶报表安装应用及实例
- SQL多字段模糊查询
- Linux下使用shell实现上传linux下某个目录下所有文件到ftp
- (王道408考研数据结构)第三章栈和队列-第三节2:栈的应用之递归
- OpenCV4每日一练day11:单目位姿估计
- 品方牌可以利用视号频做什么
- zedboard连接pmod oledrgb
- 快速创建UI控件的 方法 ,值得总结1
- 【旧文章搬运】深入分析Win7的对象引用跟踪机制
- Arduino下的marlin固件配置基础
- SATA,SAS,SSD 读写性能测试结果
- 全国计算机四级薪资,全国计算机四级通过率有多少
- 迷宫小游戏c语言代码,C语言编写的迷宫小游戏-源代码
- 域名申请步骤史上最全
- 怎么写网站的需求文档