线性代数 --- Gram-Schmidt, 格拉姆-施密特正交化(下)

news2024/10/1 15:11:06

Gram-Schmidt正交化过程

        到目前为止,我们都是在反复强调“对于无解的方程组Ax=b而言,如果矩阵A是标准正交矩阵的话,就怎么怎么好了。。。。”。因为,不论是求投影还是计算最小二乘的正规方程,他们都包含了A^{T}A。当A为标准正交矩阵Q时Q^{T}Q=I,同时可以直接求出最优解\hat{x}=Q^{T}b。从现在开始,就会详细介绍如何用Gram Schmidt的方法去构造一组彼此正交的标准正交基。

Gram-Schmidt process:

        现有一个方程组Ax=b,矩阵A由三个线性无关的列向量a,b,c组成。首先,我们基于a,b,c去构造三个相互正交的向量A,B,C。构造好以后,我们分别让A,B,C除以他们各自的长度,最终得到我们想要的一组标准正交基q_{1},q_{2},q_{3}。令q_{1},q_{2},q_{3}分别为矩阵的三个列向量,得到标准正交矩阵Q。

第一步:令向量A等于向量a,得到第一个向量A,确定了第一个方向。

A=a

第二步:因为我要构造的一组正交基是相互正交的,因此,我们的第二个向量B必须垂直于向量A。我们令b减去b在A上的投影向量p_{b},得到b在垂直于A方向上的另一个分量B。事实上,B就是垂直于投影向量p_{b}的误差向量e_{b}

B=b-p_{b}=b-\frac{A^{T}b}{A^{T}A}A

 第三步:用向量c减去c在A和B所张成的子空间(平面)上的投影p_{c},得到垂直于该子空间的另一个分量C。即,垂直于投影向量p_{c}的误差向量e_{c}。新的向量C即垂直于A也垂直于B。

C=c-p_{c}=c-\frac{A^{T}c}{A^{T}A}A-\frac{B^{T}c}{B^{T}B}B

第四步:分别对彼此正交的A,B,C进行归一化,得到向量长度都为1的一组标准正交基q_{1},q_{2},q_{3}

q_{1}=\frac{A}{\left \| A \right \|},\; q_{2}=\frac{B}{\left \| B \right \|},\; q_{3}=\frac{C}{\left \| C \right \|}

        如果还有d的话,则需要从d减去d在已经构造好的向量A,B,C三个方向上的投影(或者说是减去d在A,B,C所张成的空间上的投影),得到垂直于向量A,B,C的另一个分量D。 

        Gram Schmidt正交化过程的核心思想就是:不断的用新的已知向量,减去该向量在已经构造好的向量上的投影分量,得到我们要找的垂直分量。

即:old_vector - projection = new_vector

Example: 

最后,我们给出一个Gram Schmidt正交化计算过程的例子


 (全文完)

作者 --- 松下J27

参考文献(鸣谢):

1,Introduction to Linear Algebra,Fifth Edition - Gilbert Strang

2,线性代数及其应用,候自新,南开大学出版社 1990

(配图与本文无关)

版权声明:文中的部分图片,文字或者其他素材,可能来自很多不同的网站和说明,在此没法一一列出,如有侵权,请告知,立即删除。欢迎大家转载,但是,如果有人引用或者COPY我的文章,必须在你的文章中注明你所使用的图片或者文字来自于我的文章,否则,侵权必究。 ----松下J27

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

yolov4论文解读

数据层面上的数据增强 四张照片拼接成一张进行训练 相当于增大了batch-size,更适合于单GPU。 Mosaic data augmentation 马赛克数据增强 self-adversarial training(SAT) 自我对抗训练 DropBlock Label Smoothing 损失函数 由IOU改进到CIOU 网络结构 CSPNet&…

Win10 WLAN驱动正常但仍然不显示无线网络解决办法

Win10 WLAN驱动正常但仍然不显示无线网络解决办法 写作背景过程解决方案结尾 写作背景 本菜鸡重置了电脑的网络,然后重新启动后 WLAN 不见了,连不了 WIFI 了,很疑惑,后来经过一番搜索找到了问题所在,写下本篇文章以记…

Spark/Flink广播实现作业配置动态更新

前言 在实时计算作业中,往往需要动态改变一些配置,举几个栗子: 实时日志ETL服务,需要在日志的格式、字段发生变化时保证正常解析;实时NLP服务,需要及时识别新添加的领域词与停用词;实时风控服…

访问学者J1签证面签的七个问题

作为访问学者,申请J1签证面签时可能会遇到一些常见问题。下面知识人网小编将介绍七个访问学者面签可能遇到的问题,并提供相应的答案。 问题一:您将在美国进行何种类型的学术研究? 答案:我将在美国从事学术研究&#x…

普冉PY32L020单片机简介,主频最高48MHZ

PY32L020单片机是一颗32 位 ARM Cortex-M0内核,宽电压工作范围的 MCU。这颗MCU的价格跟八位单片机相差不大,性价比可以说是非常的高了。来看看PY32L020的配置吧。 PY32L020单片机产品特性: 内核: — 32 位 ARM Cortex - M0 — 最…

飞浆AI studio人工智能课程学习(2)-Prompt优化思路|十个技巧高效优化Prompt|迭代法|Trick法|通用法|工具辅助

文章目录 优化思路上节课的例子问题分析思路解析 Prompt优化技巧Prompt优化原理 十个技巧高效优化Prompt迭代法Trick法工具法通用技巧│定基础通用技巧│做强调需求强调怎么做? 通用技巧│提预设Trick法│戴高帽原理 Trick法│说好话以基础计算为例: Trick法│给提示…

小红书数据分析:如何用ChatGPT输出爆文笔记

ChatGPT的热度依旧不减,随着技术升级,越来越多更高级的玩法被发掘。今天我们就来聊聊,如何用ChatGPT写出小红书风格的文章。 首先,小红书笔记制作分为两个步骤: 1、找选题 2、写小红书风格的笔记 我们用例子说话&a…

全国自考本科通过率仅7%,为什么还有这么多人报考?

根据教育部官网公布的《2021年全国教育事业发展统计公报》得知:2021年,全国高等教育自学考试学历教育报考625.78万人次,取得毕业证书48.94万人。也就是说2021年我国自考平均通过率大概在7%左右。 自考通过率为什么这么低? ①自考…

Android平台外部编码数据(H264/H265/AAC/PCMA/PCMU)实时预览播放技术实现

开发背景 好多开发者可能疑惑,外部数据实时预览播放,到底有什么用? 是的,一般场景是用不到的,我们在开发这块前几年已经开发了非常稳定的RTMP、RTSP直播播放模块,不过也遇到这样的场景,部分设…

“ChatGPT之父”勇闯币圈!数十亿人的空投计划,只需交出你的虹膜?

最近,Worldcoin(世界币)热度持续提升,这个由OpenAI创始人SamAltman亲自操持的加密项目,让沉寂已久的币圈开始躁动起来。 虽然Worldcoin并非最新项目,但颇具乌托邦色彩的理念,独特的虹膜识别机制…

10几个国产免费ChatGPT网页版(内附体验网址)

文章目录 前言1. AI文本工具站效率工具自媒体创作工具代码工具 2.道和顺ChatIC3.星期五4.文心一言5.讯飞星火认知大模型6.通义千问7.商汤-日日新8.Moss9.ChatGLM10. 360智脑写在最后 前言 随着ChatGPT迅速走红,国内各大企业纷纷发力认知大模型领域。经过一段时间的酝酿&#x…

剩余参数、扩展运算符、将伪数组转换伪数组

剩余参数 剩余参数语法允许我们将一个不定数量的参数表示为一个数组 语法: array.forEach(function(currentValue,index,arr){}) currentValue:数组当前项的值index:数组当前项的索引arr:数组对象本身 示例 <!DOCTYPE html> <html lang"en"><h…

实力认证!开源网安入选《2023年中国网络安全市场全景图》

近日&#xff0c;数说安全正式发布《2023年中国网络安全市场全景图》&#xff08;以下简称全景图&#xff09;。开源网安凭借强大的技术能力与服务经验上榜本次全景图开发安全/DevSecOps、应用安全测试&#xff08;AST&#xff09;、软件成分分析&#xff08;SCA&#xff09;、…

野火STM32电机系列(三)Cubemx配置CAN通信

CAN接口: PI9 PB9 1.配置CAN 通信参数 由于F4的 CAN外设挂载在APB1上&#xff0c;时钟配置后APB1的时钟速率为42MHz&#xff0c;目标通信速率为1000KHz&#xff0c;由公式&#xff1a; BaudRate 1/NominalBitTime NominalBitTime 1tq tBS1 tBS2 设置参数如下&#xff1a;…

pip 安装 pytorch

一.使用pip安装 pytorch pytorch博客教程 最好先创建一个虚拟环境&#xff0c;因为如果在同一环境&#xff0c;升级某一个包可能会导致另一个包无法使用&#xff0c;比如 a的包是1.0&#xff0c;b的包是1.0依赖a的1.0&#xff0c;然后a升级为1.1&#xff0c;b可能无法使用 创…

Typora 含多图片笔记快速上传到 CSDN 上发表

Typora 含多图片笔记快速上传到 CSDN 上发表 适用场景解决方案具体步骤 适用场景 在 Typora 里面记笔记,上传的图片是本地保存的,如果要将笔记上传到 CSDN 上发表的话,图片得一张一张地拖拽非常麻烦 解决方案 Typora PicGo Gitee 具体步骤 先安装 PicGo 和 node.js, 创建…

开源之夏 2023 | 欢迎报名Rust相关项目

开源之夏是中国科学院软件研究所联合openEuler发起的开源软件供应链点亮计划系列暑期活动&#xff0c;旨在鼓励在校学生积极参与开源软件的开发维护&#xff0c;促进优秀开源软件社区的蓬勃发展。活动联合各大开源社区&#xff0c;针对重要开源软件的开发与维护提供项目&#x…

java+springboot 做日志链路追踪

一、 为什么要做日志链路追踪 日志链路追踪&#xff08;Log Path Tracing&#xff09;是Spring Boot项目的一项关键功能&#xff0c;它通过将日志消息的源头与其对应的请求或响应路径相关联&#xff0c;实现对日志数据的可视化跟踪。随着项目规模的扩大和复杂性的增加&#xf…

Padding, Spacer, Initializer 的使用

1. Padding 的使用 1.1 样式一 1) 实现 func testText1()-> some View{Text("Hello, World!").background(Color.yellow) // 背景颜色//.padding() // 默认间距.padding(.all, 10) // 所有的间距.padding(.leading, 20) // 开始的间距.ba…

真题详解(数字签名算法)-软件设计(七十八)

真题详解(有限自动机)-软件设计&#xff08;七十七)https://blog.csdn.net/ke1ying/article/details/130748759 可用于数字签名算法的是_____。 答案&#xff1a;非对称RSA 移植性&#xff1a;易安装、易替换、适应性。 UML状态图转换不正确的是______。 活动可以在转换时执…