MIT线性代数笔记-第17讲-正交矩阵,Schmidt正交化

17.正交矩阵, S c h m i d t Schmidt Schmidt正交化

  1. “标准”经常表示单位长度

  2. 标准正交基:由两两正交的单位向量组成的基

    将标准正交基中的元素记作 q ⃗ 1 , q ⃗ 2 , ⋯ , q ⃗ n \vec{q}_1 , \vec{q}_2 , \cdots , \vec{q}_n q 1,q 2,,q n

    q ⃗ i T q ⃗ j = { 0 , i ≠ j 1 , i = j \vec{q}_i^T \vec{q}_j = \left\{\begin{matrix} 0 , i \ne j \\ 1,i = j \end{matrix}\right. q iTq j={0,i=j1,i=j

    Q = [ ∣ ⋯ ∣ q ⃗ 1 ⋯ q ⃗ n ∣ ⋯ ∣ ] Q = \begin{bmatrix} | & \cdots & | \\ \vec{q}_1 & \cdots & \vec{q}_n \\ | & \cdots & | \end{bmatrix} Q= q 1q n (此时 Q Q Q不一定要是方阵),则 Q T Q = [ − q ⃗ 1 T − ⋯ ⋯ ⋯ − q ⃗ n T − ] [ ∣ ⋯ ∣ q ⃗ 1 ⋯ q ⃗ n ∣ ⋯ ∣ ] = I Q^T Q = \begin{bmatrix} - & \vec{q}_1^T & - \\ \cdots & \cdots & \cdots \\ - & \vec{q}_n^T & - \end{bmatrix} \begin{bmatrix} | & \cdots & | \\ \vec{q}_1 & \cdots & \vec{q}_n \\ | & \cdots & | \end{bmatrix} = I QTQ= q 1Tq nT q 1q n =I

  3. 正交矩阵:乘上转置得到 I I I的方阵,记作 Q Q Q

    证明正交矩阵各列构成一个标准正交基:

    ​    当 i ≠ j i \ne j i=j时, Q T Q Q^T Q QTQ的第 i i i行第 j j j列等于 0 0 0,即 q ⃗ i T q ⃗ j = 0 \vec{q}_i^T \vec{q}_j = 0 q iTq j=0,说明 Q Q Q各列两两正交,又 Q T Q Q^T Q QTQ的第 i i i行第 i i i列为 1 1 1,即 q ⃗ i T q ⃗ i = q ⃗ i 2 = 1 \vec{q}_i^T \vec{q}_i = \vec{q}_i^2 = 1 q iTq i=q i2=1,可得 ∣ q ⃗ i ∣ = 1 |\vec{q}_i| = 1 q i=1,说明 Q Q Q各列均为单位向量,因而 Q Q Q各列为两两正交的单位向量

    由定义可知 Q T = Q − 1 Q^T = Q^{-1} QT=Q1,所以 Q Q T = Q Q − 1 = I Q Q^T = Q Q^{-1} = I QQT=QQ1=I,因而 Q T , Q − 1 Q^T , Q^{-1} QT,Q1也为正交矩阵

    Q Q Q的列空间的投影矩阵 P = Q ( Q T Q ) − 1 Q T = Q Q T = I P = Q (Q^T Q)^{-1} Q^T = Q Q^T = I P=Q(QTQ)1QT=QQT=I

    因而求解 Q x ⃗ = b ⃗ Q \vec{x} = \vec{b} Qx =b 时,由 Q T b ⃗ = Q T Q x ⃗ Q^T \vec{b} = Q^T Q \vec{x} QTb =QTQx 可得 x ⃗ = Q T b ⃗ \vec{x} = Q^T \vec{b} x =QTb

    所以 b ⃗ \vec{b} b 在第 i i i个基方向上的投影 x i = r o w i o f Q ⋅ b ⃗ = q ⃗ i T b ⃗ x_i = row\ i\ of\ Q \cdot \vec{b} = \vec{q}_i^T \vec{b} xi=row i of Qb =q iTb

    • 哈达玛矩阵:由 1 , − 1 1 , -1 1,1构成的满足右乘其转置得到 n I nI nI的方阵,记作 H H H

      易证 n I nI nI中的 n n n H H H的阶数

      与正交矩阵同理可知 H H H各列两两正交,进而由定义可知 1 n H \dfrac{1}{\sqrt{n}} H n 1H为正交矩阵

      哈达玛矩阵的阶数只能为 1 , 2 1 , 2 1,2 4 4 4的倍数(但并不是所有 4 4 4的倍数都可以) (暂时不会证明) \color{OrangeRed}(暂时不会证明) (暂时不会证明)

  4. 格拉姆-施密特正交化

    • 考虑两个线性无关向量 a ⃗ , b ⃗ \vec{a} , \vec{b} a ,b ,先将二者正交化得到两个新向量 A ⃗ , B ⃗ \vec{A} , \vec{B} A ,B

      正交化具体方法为:令其中一个保持不变,如 A ⃗ = a ⃗ \vec{A} = \vec{a} A =a ,将另一个向量投影到前者后用后者减去投影,得到的就是后者的正交化,如 B ⃗ = b ⃗ − A ⃗ A ⃗ T b ⃗ A ⃗ T A ⃗ \vec{B} = \vec{b} - \vec{A} \dfrac{\vec{A}^T \vec{b}}{\vec{A}^T \vec{A}} B =b A A TA A Tb

      再将正交化得到的向量标准化,即每个向量除以自己的模,得到一个标准正交基 { q ⃗ 1 , q ⃗ 2 } \{ \vec{q}_1 , \vec{q}_2 \} {q 1,q 2},如 q ⃗ 1 = A ⃗ ∣ A ⃗ ∣ , q ⃗ 2 = B ⃗ ∣ B ⃗ ∣ \vec{q}_1 = \dfrac{\vec{A}}{|\vec{A}|} , \vec{q}_2 = \dfrac{\vec{B}}{|\vec{B}|} q 1=A A ,q 2=B B

    • 考虑将线性无关向量拓展至三个—— a ⃗ , b ⃗ , c ⃗ \vec{a} , \vec{b} , \vec{c} a ,b ,c

      正交化得到的向量增加一个 C ⃗ \vec{C} C

      C ⃗ ′ = c ⃗ − A ⃗ A ⃗ T c ⃗ A ⃗ T A ⃗ \vec{C}^{'} = \vec{c} - \vec{A} \dfrac{\vec{A}^T \vec{c}}{\vec{A}^T \vec{A}} C =c A A TA A Tc ,有 ( C ⃗ ′ ) T A ⃗ = 0 (\vec{C}^{'})^T \vec{A} = 0 (C )TA =0

      C ⃗ ′ ′ = C ⃗ ′ − B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ \vec{C}^{''} = \vec{C}^{'} - \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}} C ′′=C B B TB B TC ,有 C ′ ′ ⃗ T B ⃗ = 0 \vec{C^{''}}^T \vec{B} = 0 C′′ TB =0

      其中 B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}} B B TB B TC B ⃗ \vec{B} B 共线,而 B ⃗ \vec{B} B A ⃗ \vec{A} A 正交,因而 B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}} B B TB B TC A ⃗ \vec{A} A 正交,又 ( C ⃗ ′ ) T A ⃗ = 0 (\vec{C}^{'})^T \vec{A} = 0 (C )TA =0,所以 ( C ⃗ ′ ′ ) T A ⃗ = ( C ⃗ ′ − B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ ) A ⃗ = 0 (\vec{C}^{''})^T \vec{A} = (\vec{C}^{'} - \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}}) \vec{A} = 0 (C ′′)TA =(C B B TB B TC )A =0,即 C ⃗ ′ ′ \vec{C}^{''} C ′′ A ⃗ \vec{A} A 正交,所以 C ⃗ ′ ′ \vec{C}^{''} C ′′ B ⃗ , A ⃗ \vec{B} , \vec{A} B ,A 均正交, C ⃗ = C ⃗ ′ ′ \vec{C} = \vec{C}^{''} C =C ′′

      下证 B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ = B ⃗ B ⃗ T c ⃗ B ⃗ T B ⃗ \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}} = \vec{B} \dfrac{\vec{B}^T \vec{c}}{\vec{B}^T \vec{B}} B B TB B TC =B B TB B Tc

      等式左右分别表示 C ⃗ ′ , c ⃗ \vec{C}^{'} , \vec{c} C ,c B ⃗ \vec{B} B 的投影, c ⃗ \vec{c} c 的投影由 C ⃗ ′ , A ⃗ A ⃗ T c ⃗ A ⃗ T A ⃗ \vec{C}^{'} , \vec{A} \dfrac{\vec{A}^T \vec{c}}{\vec{A}^T \vec{A}} C ,A A TA A Tc 的投影两部分组成,而后者与 A ⃗ \vec{A} A 共线,即与 B ⃗ \vec{B} B 正交,投影为 0 ⃗ \vec{0} 0 ,所以 c ⃗ \vec{c} c 的投影与 C ⃗ ′ \vec{C}^{'} C 的投影一致,即 B ⃗ B ⃗ T C ⃗ ′ B ⃗ T B ⃗ = B ⃗ B ⃗ T c ⃗ B ⃗ T B ⃗ \vec{B} \dfrac{\vec{B}^T \vec{C}^{'}}{\vec{B}^T \vec{B}} = \vec{B} \dfrac{\vec{B}^T \vec{c}}{\vec{B}^T \vec{B}} B B TB B TC =B B TB B Tc

      ∴ C ⃗ = c ⃗ − A ⃗ A ⃗ T c ⃗ A ⃗ T A ⃗ − B ⃗ B ⃗ T c ⃗ B ⃗ T B ⃗ \therefore \vec{C} = \vec{c} - \vec{A} \dfrac{\vec{A}^T \vec{c}}{\vec{A}^T \vec{A}} - \vec{B} \dfrac{\vec{B}^T \vec{c}}{\vec{B}^T \vec{B}} C =c A A TA A Tc B B TB B Tc (如果 A ⃗ , B ⃗ \vec{A} , \vec{B} A ,B 不正交,是无法这样得到与两者都正交的分向量的)

      当然最后还是要除以模

    • 将最初和最终的向量分别按列摆放组合为两个矩阵 A = [ ∣ ∣ ⋯ a ⃗ 1 a ⃗ 2 ⋯ ∣ ∣ ⋯ ] , Q = [ ∣ ∣ ⋯ q ⃗ 1 q ⃗ 2 ⋯ ∣ ∣ ⋯ ] A = \begin{bmatrix} | & | & \cdots \\ \vec{a}_1 & \vec{a}_2 & \cdots \\ | & | & \cdots \end{bmatrix} , Q = \begin{bmatrix} | & | & \cdots \\ \vec{q}_1 & \vec{q}_2 & \cdots \\ | & | & \cdots \end{bmatrix} A= a 1a 2 ,Q= q 1q 2

      A = Q R A = QR A=QR,有 Q T A = Q T Q R Q^T A = Q^T QR QTA=QTQR,可得 R = Q T A R = Q^T A R=QTA

      思考得到每个 q ⃗ \vec{q} q 的过程不难发现 a ⃗ 1 , a ⃗ 2 , ⋯ , a ⃗ i \vec{a}_1 , \vec{a}_2 , \cdots , \vec{a}_i a 1,a 2,,a i q ⃗ 1 , q ⃗ 2 , ⋯ , q ⃗ i \vec{q}_1 , \vec{q}_2 , \cdots , \vec{q}_i q 1,q 2,,q i线性组合所能得到的空间是一致的,而 q i q_{i} qi之后的向量会正交于 q ⃗ 1 , q ⃗ 2 , ⋯ , q ⃗ i \vec{q}_1 , \vec{q}_2 , \cdots , \vec{q}_i q 1,q 2,,q i中的每一个,即正交于这个空间中的每一个向量,因此想要对 q ⃗ \vec{q} q 线性组合得到 a ⃗ i \vec{a}_i a i,无法使用到 q ⃗ i \vec{q}_i q i之后的 q ⃗ \vec{q} q ,即它们在 R R R中的对应参数为 0 0 0,所以 R R R为一个上三角阵,且其对角线上的元素为对应的未标准化的新向量的模


打赏

制作不易,若有帮助,欢迎打赏!
赞赏码

支付宝付款码

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

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

相关文章

基于FactoryBean、实例工厂、静态工厂创建Spring中的复杂对象

😉😉 学习交流群: ✅✅1:这是孙哥suns给大家的福利! ✨✨2:我们免费分享Netty、Dubbo、k8s、Mybatis、Spring...应用和源码级别的视频资料 🥭🥭3:QQ群:583783…

会声会影2024旗舰版系统配置要求及格式支持

会声会影2024旗舰版是一款广受欢迎的视频编辑软件,它的最新版本,会声会影2023,已经发布。在这篇文章中,我们将探讨会声会影2024旗舰版系统配置要求及格式支持 会声会影2024是一款专业的视频剪辑软件,能够帮助用户制作高…

geemap学习笔记016:获取图像的基本属性和描述性信息

前言 遥感数据中通常包含众多信息,例如图像获取的时间、云覆盖量、以及每个波段的最大值最小值等等。 1 导入库并显示地图 import ee import geemapMap geemap.Map() Map2 添加图像数据 centroid ee.Geometry.Point([-122.4439, 37.7538]) #创建一个点坐标lan…

springcloud==openfeign

单独使用 创建一个服务端 import org.springframework.boot.SpringApplication; import org.springframework.boot.autoconfigure.SpringBootApplication; import org.springframework.web.bind.annotation.GetMapping; import org.springframework.web.bind.annotation.Path…

浅析函数防抖节流

防抖和节流都是前端开发中常用的优化性能的技术。 一、定义 防抖: 防抖指的是在事件触发后,在规定的时间内若再次触发,则重新计时,直到规定时间内没有再次触发事件,才执行事件处理。这样可以避免在短时间内频繁地触发…

Redis实战命令

实战命令 单值缓存 set key value get key 对象缓存 (1)set user:1 value(json格式) (2)mset user:1:name junfeng user:1:age 18 mget user:1:name user:1:age 分布式锁 分布式锁解决了什么问题? 分布式锁解…

2 线、3 线和 4 线 RTD 配置之间有什么区别?

电阻温度检测器 (RTD) 是温度传感器的一种,由于其准确性、可重复性和稳定性而广泛应用于各种工业应用。这些设备通过感测材料温度变化时电阻的变化来测量温度。 RTD 探头有多种配置,包括 2 线、3 线和 4 线型号。这些类型之间存在显着差异,在…

智能优化算法应用:基于回溯搜索算法无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用:基于回溯搜索算法无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用:基于回溯搜索算法无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.回溯搜索算法4.实验参数设定5.算法结果6.参考…

ubuntu改window任务栏

经常在ubuntu和win之间切换,任务栏的布局不统一会让人很别扭,个人很喜欢win任务栏的不折叠图标功能,而ubuntu没有,又很喜欢的ubuntu的多工作空间,效率比副屏还高,还可以自定义切换工作空间的快捷键。鱼和熊…

2023-11-28-直播单细胞图表美化-seurat数据结构 featureplot dotplot vlnplot

单细胞常见的可视化方式有DimPlot,FeaturePlot ,DotPlot ,VlnPlot 和 DoHeatmap几种 ,Seurat中均可以很简单的实现,但是文献中的图大多会精美很多。 之前 跟SCI学umap图| ggplot2 绘制umap图,坐标位置 &am…

[SaaS] 广告创意中stable-diffusion的应用

深度对谈:广告创意领域中 AIGC 的应用这个领域非常快速发展,所以你应该保持好奇心,不断尝试新事物,不断挑战自己。https://mp.weixin.qq.com/s/ux9iEABNois3y4wwyaDzAQ我对AIGC领域应用调研,除了MaaS服务之外&#xff…

安全防控 | AIRIOT智能安防管理解决方案

现代社会对安全和便捷性的需求越来越高,特别是在大型商业园区、住宅社区和办公大楼等场所。传统的安防系统往往存在一些痛点: 通行效率问题:传统门禁系统通常导致人员排队等待,降低了通行效率。车辆通行管理不当会导致交通拥堵和停车问题。 …

Elasticsearch 线上实战问题及解决方案探讨

1、reindex相关问题 1.1 问题描述 我有 1tb 的一个大索引若干,要迁移到另外一个新集群去,有没有好办法?reindex好像会中断...... reindex 是不是就算设置了频率也会莫名的中断,而且没地方查到错误?1000多万的数据&…

【axios封装】万字长文,TypeScript实战,封装一个axios - 基础封装篇

目录 前言版本环境变量配置引入的类型1、AxiosIntance: axios实例类型2、InternalAxiosRequestConfig: 高版本下AxiosRequestConfig的拓展类型3、AxiosRequestConfig: 请求体配置参数类型4、AxiosError: 错误对象类型5、AxiosResponse: 完整原始响应体类型 目标效果开始封装骨架…

鸿运主动安全监控云平台存在任意文件读取漏洞 附POC

@[toc] 鸿运主动安全监控云平台存在任意文件读取漏洞 附POC 免责声明:请勿利用文章内的相关技术从事非法测试,由于传播、利用此文所提供的信息或者工具而造成的任何直接或者间接的后果及损失,均由使用者本人负责,所产生的一切不良后果与文章作者无关。该文章仅供学习用途…

蓝桥杯物联网竞赛_STM32L071_5_串口接收发送数据

理论: 串口采取异步通信,即不依赖时钟节拍来接收或发送数据,而是采用互相约定的波特率传输数据。 波特率与单位时间传输的比特数有关,波特率越大传输的数据越多 传输一个比特花费的时间T 1 / 比特率 接受和发送数据的时候需要…

leetcode:455. 分发饼干

一、题目 原题链接:455. 分发饼干 - 力扣(LeetCode) 函数原型:int findContentChildren(int* g, int gSize, int* s, int sSize) 二、思路 将胃口数组和饼干尺寸数组降序排序 遍历两个数组,从胃口数组中找到符合饼干尺…

网络入门---网络的大致了解

目录标题 网络发展的简单认识协议作用的理解协议的本质什么是协议分层网络通信所面对的问题OSI七层模型TCP/IP模型协议报头的理解局域网通信局域网通信基本原理报头的问题局域网的特点跨网的网络链接如何查看mac地址 网络发展的简单认识 通过之前的学习我们知道计算机是给人提…

使用docker-compose优雅部署nacos

查看代码中引入nacos版本 在应用的pom.xml中搜索nacos关键字,找到相关的nacos依赖 点击以来左边的图标,找到依赖管理器中的pom.xml,并全局搜索nacos,即可找到对应的nacos客户端版本 使用docker-compose部署nacos version: 3s…

Jmeter-分布式压测(远程启动服务器,windows)

1 前提条件 JDK已部署,版本一致Jmeter已部署,版本一致多台服务器连接的同一网络(例如:同一wifi)防火墙处于关闭状态(或者对应默认端口处于开放状态)虚拟网络适配器都处于关闭状态查找到每一台服务器的IP 2 主服务器配…