线性代数|机器学习-P4正交矩阵中的标准正交向量

2024-06-03 10:44

本文主要是介绍线性代数|机器学习-P4正交矩阵中的标准正交向量,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

文章目录

  • 1. 正交矩阵Q
    • 1.1 矩阵Q的推导
    • 1.2 |Qx|=|x|
  • 2. 常见正交矩阵
    • 2.1 旋转矩阵
    • 2.2 镜像矩阵
    • 2.3 Householder矩阵
    • 2.4 Hadamard矩阵
    • 2.4 小波矩阵
    • 2.5 傅里叶级数矩阵

1. 正交矩阵Q

1.1 矩阵Q的推导

方阵A正交的充要条件是A的行(列)向量组是单位正交向量组. 我们定义正交矩阵Q表示如下:
Q = [ q 1 q 2 ⋯ q n ] ; q i T q j = { 1 , i = j 0 , i ≠ j \begin{equation} Q=\begin{bmatrix}q_1&q_2&\cdots&q_n\end{bmatrix};q_i^Tq_j=\left\{ \begin{aligned} 1&,i=j \\ 0&,i\neq j \\ \end{aligned} \right. \end{equation} Q=[q1q2qn]qiTqj={10,i=j,i=j

  • 对矩阵Q展开可得:
    Q T Q = [ q 1 T q 2 T ⋮ q n T ] [ q 1 q 2 ⋯ q n ] = [ q 1 T q 1 q 1 T q 2 ⋯ q 1 T q n q 2 T q 1 q 2 T q 2 ⋯ q 2 T q n ⋮ ⋮ ⋮ ⋮ q n T q 1 q n T q 2 ⋯ q n T q n ] = [ 1 0 ⋯ 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋮ 0 0 ⋯ 1 ] \begin{equation} Q^TQ=\begin{bmatrix}q_1^T\\\\q_2^T\\\\\vdots \\\\q_n^T\end{bmatrix}\begin{bmatrix}q_1&q_2&\cdots &q_n\end{bmatrix}=\begin{bmatrix} q_1^Tq_1&q_1^Tq_2&\cdots&q_1^Tq_n\\\\ q_2^Tq_1&q_2^Tq_2&\cdots&q_2^Tq_n\\\\ \vdots&\vdots&\vdots&\vdots\\\\ q_n^Tq_1&q_n^Tq_2&\cdots&q_n^Tq_n \end{bmatrix}=\begin{bmatrix} 1&0&\cdots&0\\\\ 0&1&\cdots&0\\\\ \vdots&\vdots&\vdots&\vdots\\\\ 0&0&\cdots&1 \end{bmatrix} \end{equation} QTQ= q1Tq2TqnT [q1q2qn]= q1Tq1q2Tq1qnTq1q1Tq2q2Tq2qnTq2q1Tqnq2TqnqnTqn = 100010001

1.2 |Qx|=|x|

首先我们知道对于一个正交矩阵Q来说,满足: Q T Q = I Q^TQ=I QTQ=I

  • 两边分别乘以 x T , x x^T,x xT,x
    x T Q T Q x = x T x → ( Q x ) T ( Q x ) = x T x → ∣ Q x ∣ 2 = ∣ x ∣ 2 → ∣ Q x ∣ = ∣ x ∣ \begin{equation} x^TQ^TQx=x^Tx\rightarrow (Qx)^T(Qx)=x^Tx\rightarrow |Qx|^2=|x|^2\rightarrow |Qx|=|x| \end{equation} xTQTQx=xTx(Qx)T(Qx)=xTxQx2=x2Qx=x
  • 小结:也就是说,对于任意向量x来说,我们设计一个正交矩阵Q,在不断地左乘矩阵Q的时候,因为 ∣ Q x ∣ = ∣ x ∣ |Qx|=|x| Qx=x的原因,所以矩阵大小值不会变,这点在计算编程中至关重要,这样的话,我们就可以保证矩阵相乘的时候值不会溢出。

2. 常见正交矩阵

对称矩阵和正交矩阵的特征向量组为正交单位向量,所以我们为了方便后续的快速计算,我们希望直接找正交矩阵和对称矩阵。

2.1 旋转矩阵

  • 我们将向量在保持大小不变的情况下,逆时针旋转 θ \theta θ, 我们分别看看向量 x 1 = [ 0 , 1 ] T , x 2 = [ 1 , 0 ] T x_1=[0,1]^T,x_2=[1,0]^T x1=[0,1]T,x2=[1,0]T变化结果,如图所述:
    A = [ 1 0 ] → A 1 = [ cos ⁡ θ sin ⁡ θ ] ; B = [ 0 1 ] → B 1 = [ − sin ⁡ θ cos ⁡ θ ] \begin{equation} A=\begin{bmatrix} 1\\\\0 \end{bmatrix}\rightarrow A_1=\begin{bmatrix} \cos{\theta}\\\\\sin{\theta} \end{bmatrix};B=\begin{bmatrix} 0\\\\1 \end{bmatrix}\rightarrow B_1=\begin{bmatrix} -\sin{\theta}\\\\\cos{\theta} \end{bmatrix} \end{equation} A= 10 A1= cosθsinθ ;B= 01 B1= sinθcosθ
  • 可得旋转矩阵Q表示如下:
    Q = [ cos ⁡ θ − sin ⁡ θ sin ⁡ θ cos ⁡ θ ] \begin{equation} Q=\begin{bmatrix} \cos{\theta}&-\sin{\theta}\\\\ \sin{\theta}&\cos{\theta} \end{bmatrix} \end{equation} Q= cosθsinθsinθcosθ
    在这里插入图片描述

2.2 镜像矩阵

  • 我们将向量在保持大小不变的情况下,沿着 1 2 θ \frac{1}{2}\theta 21θ镜像, 我们分别看看向量 x 1 = [ 0 , 1 ] T , x 2 = [ 1 , 0 ] T x_1=[0,1]^T,x_2=[1,0]^T x1=[0,1]T,x2=[1,0]T变化结果,如图所述:
    A = [ 1 0 ] → A 1 = [ cos ⁡ θ sin ⁡ θ ] ; B = [ 0 1 ] → B 1 = [ sin ⁡ θ − cos ⁡ θ ] \begin{equation} A=\begin{bmatrix} 1\\\\0 \end{bmatrix}\rightarrow A_1=\begin{bmatrix} \cos{\theta}\\\\\sin{\theta} \end{bmatrix};B=\begin{bmatrix} 0\\\\1 \end{bmatrix}\rightarrow B_1=\begin{bmatrix} \sin{\theta}\\\\-\cos{\theta} \end{bmatrix} \end{equation} A= 10 A1= cosθsinθ ;B= 01 B1= sinθcosθ
  • 可得旋转矩阵Q表示如下:
    R e f l e c t i o n = [ cos ⁡ θ sin ⁡ θ sin ⁡ θ − cos ⁡ θ ] \begin{equation} Reflection=\begin{bmatrix} \cos{\theta}&\sin{\theta}\\\\ \sin{\theta}&-\cos{\theta} \end{bmatrix} \end{equation} Reflection= cosθsinθsinθcosθ
    在这里插入图片描述

2.3 Householder矩阵

householder变换的作用是将向量x的第一项值为|x|,其他值不变,相当于将x通过镜面反射得到y

  • x:为我们输入的向量x

  • y:为经过householder变换后向量y

  • w: x-y=2w,等腰三角形底边相等

  • u:定义跟w平行的单位向量u,|u|=1;
    在这里插入图片描述

  • 根据 u T x u^Tx uTx公式可得,|u|=1
    u T x = ∣ x ∣ ∣ u T ∣ cos ⁡ α = ∣ x ∣ cos ⁡ α = ∣ w ∣ \begin{equation} u^Tx=|x||u^T|\cos{\alpha}=|x|\cos{\alpha}=|w| \end{equation} uTx=x∣∣uTcosα=xcosα=w

  • 向量w,x,y之间的关系如下:
    y + 2 w = x , 2 w = 2 ∣ w ∣ u = 2 u T x u \begin{equation} y+2w=x,2w=2|w|u=2u^Txu \end{equation} y+2w=x,2w=2∣wu=2uTxu

  • 因为 u T x u^Tx uTx为一个数,可以任意放位置
    y + 2 u u T x = x → y = ( I − 2 u u T ) x \begin{equation} y+2uu^Tx=x\rightarrow y=(I-2uu^T)x \end{equation} y+2uuTx=xy=(I2uuT)x

  • 这里我们可以将H定义为如下:
    y = H x ; H = I − 2 u u T ; y = ( I − 2 u u T ) x \begin{equation} y=Hx;H=I-2uu^T;y=(I-2uu^T)x \end{equation} y=Hx;H=I2uuT;y=(I2uuT)x

  • Householder矩阵,对称性,正交性。
    H = I − 2 u u T → H T = I − 2 u u T = H , H T H = I − 2 u u T − 2 u u T + 4 u u T u u T = I \begin{equation} H=I-2uu^T\rightarrow H^T=I-2uu^T=H,H^TH=I-2uu^T-2uu^T+4uu^Tuu^T=I \end{equation} H=I2uuTHT=I2uuT=H,HTH=I2uuT2uuT+4uuTuuT=I

2.4 Hadamard矩阵

哈达玛(Hadamard)矩阵;

  • H 2 , H 4 H_2,H_4 H2,H4
    H 2 = [ 1 1 1 − 1 ] ; H 4 = [ H 2 H 2 H 2 − H 2 ] = [ 1 1 1 1 1 − 1 1 − 1 1 1 − 1 − 1 1 − 1 − 1 1 ] \begin{equation} H_2=\begin{bmatrix} 1&1\\\\ 1&-1 \end{bmatrix};H_4=\begin{bmatrix} H_2&H_2\\\\ H_2&-H_2 \end{bmatrix}= \begin{bmatrix} 1&1&1&1\\\\ 1&-1&1&-1\\\\ 1&1&-1&-1\\\\ 1&-1&-1&1\\\\ \end{bmatrix} \end{equation} H2= 1111 ;H4= H2H2H2H2 = 1111111111111111

2.4 小波矩阵

我们知道小波矩阵也是一个正交矩阵
在这里插入图片描述

2.5 傅里叶级数矩阵

我们这里要引入复数i,这里先定义,后面再详细展开
F 4 = [ 1 1 1 1 1 i i 2 i 3 1 i 2 i 4 i 6 1 i 3 i 6 i 9 ] \begin{equation} F_4=\begin{bmatrix} 1&1&1&1\\\\ 1&i&i^2&i^3\\\\ 1&i^2&i^4&i^6\\\\ 1&i^3&i^6&i^9 \end{bmatrix} \end{equation} F4= 11111ii2i31i2i4i61i3i6i9

  • 注意,这里第二列和第四列之间正交需要取共轭复数,具体如下:
    f 1 = [ 1 i i 2 i 3 ] = [ 1 i − 1 − i ] → f 1 H = [ 1 − i − 1 i ] \begin{equation} f_1=\begin{bmatrix} 1\\\\i\\\\i^2\\\\i^3 \end{bmatrix}=\begin{bmatrix} 1\\\\i\\\\-1\\\\-i \end{bmatrix}\rightarrow f_1^H=\begin{bmatrix}1&-i&-1&i\end{bmatrix} \end{equation} f1= 1ii2i3 = 1i1i f1H=[1i1i]
    f 4 = [ 1 i 3 i 6 i 9 ] = [ 1 − i 1 − i ] → f 1 H f 4 = [ 1 − i − 1 i ] [ 1 − i 1 − i ] = 1 + i 2 − 1 − i 2 = 0 \begin{equation} f_4=\begin{bmatrix} 1\\\\i^3\\\\i^6\\\\i^9 \end{bmatrix}=\begin{bmatrix} 1\\\\-i\\\\1\\\\-i \end{bmatrix}\rightarrow f_1^Hf_4=\begin{bmatrix}1&-i&-1&i\end{bmatrix}\begin{bmatrix} 1\\\\-i\\\\1\\\\-i \end{bmatrix}=1+i^2-1-i^2=0 \end{equation} f4= 1i3i6i9 = 1i1i f1Hf4=[1i1i] 1i1i =1+i21i2=0
  • 同理,可以逐个计算得到 F 4 F_4 F4的列向量相互正交。

这篇关于线性代数|机器学习-P4正交矩阵中的标准正交向量的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1026759

相关文章

51单片机学习记录———定时器

文章目录 前言一、定时器介绍二、STC89C52定时器资源三、定时器框图四、定时器模式五、定时器相关寄存器六、定时器练习 前言 一个学习嵌入式的小白~ 有问题评论区或私信指出~ 提示:以下是本篇文章正文内容,下面案例可供参考 一、定时器介绍 定时器介绍:51单片机的定时器属于单片机的内部资源,其电路的连接和运转均在单片机内部完成。 定时器作用: 1.用于计数系统,可

问题:第一次世界大战的起止时间是 #其他#学习方法#微信

问题:第一次世界大战的起止时间是 A.1913 ~1918 年 B.1913 ~1918 年 C.1914 ~1918 年 D.1914 ~1919 年 参考答案如图所示

[word] word设置上标快捷键 #学习方法#其他#媒体

word设置上标快捷键 办公中,少不了使用word,这个是大家必备的软件,今天给大家分享word设置上标快捷键,希望在办公中能帮到您! 1、添加上标 在录入一些公式,或者是化学产品时,需要添加上标内容,按下快捷键Ctrl+shift++就能将需要的内容设置为上标符号。 word设置上标快捷键的方法就是以上内容了,需要的小伙伴都可以试一试呢!

AssetBundle学习笔记

AssetBundle是unity自定义的资源格式,通过调用引擎的资源打包接口对资源进行打包成.assetbundle格式的资源包。本文介绍了AssetBundle的生成,使用,加载,卸载以及Unity资源更新的一个基本步骤。 目录 1.定义: 2.AssetBundle的生成: 1)设置AssetBundle包的属性——通过编辑器界面 补充:分组策略 2)调用引擎接口API

Javascript高级程序设计(第四版)--学习记录之变量、内存

原始值与引用值 原始值:简单的数据即基础数据类型,按值访问。 引用值:由多个值构成的对象即复杂数据类型,按引用访问。 动态属性 对于引用值而言,可以随时添加、修改和删除其属性和方法。 let person = new Object();person.name = 'Jason';person.age = 42;console.log(person.name,person.age);//'J

大学湖北中医药大学法医学试题及答案,分享几个实用搜题和学习工具 #微信#学习方法#职场发展

今天分享拥有拍照搜题、文字搜题、语音搜题、多重搜题等搜题模式,可以快速查找问题解析,加深对题目答案的理解。 1.快练题 这是一个网站 找题的网站海量题库,在线搜题,快速刷题~为您提供百万优质题库,直接搜索题库名称,支持多种刷题模式:顺序练习、语音听题、本地搜题、顺序阅读、模拟考试、组卷考试、赶快下载吧! 2.彩虹搜题 这是个老公众号了 支持手写输入,截图搜题,详细步骤,解题必备

《offer来了》第二章学习笔记

1.集合 Java四种集合:List、Queue、Set和Map 1.1.List:可重复 有序的Collection ArrayList: 基于数组实现,增删慢,查询快,线程不安全 Vector: 基于数组实现,增删慢,查询快,线程安全 LinkedList: 基于双向链实现,增删快,查询慢,线程不安全 1.2.Queue:队列 ArrayBlockingQueue:

硬件基础知识——自学习梳理

计算机存储分为闪存和永久性存储。 硬盘(永久存储)主要分为机械磁盘和固态硬盘。 机械磁盘主要靠磁颗粒的正负极方向来存储0或1,且机械磁盘没有使用寿命。 固态硬盘就有使用寿命了,大概支持30w次的读写操作。 闪存使用的是电容进行存储,断电数据就没了。 器件之间传输bit数据在总线上是一个一个传输的,因为通过电压传输(电流不稳定),但是电压属于电势能,所以可以叠加互相干扰,这也就是硬盘,U盘

人工智能机器学习算法总结神经网络算法(前向及反向传播)

1.定义,意义和优缺点 定义: 神经网络算法是一种模仿人类大脑神经元之间连接方式的机器学习算法。通过多层神经元的组合和激活函数的非线性转换,神经网络能够学习数据的特征和模式,实现对复杂数据的建模和预测。(我们可以借助人类的神经元模型来更好的帮助我们理解该算法的本质,不过这里需要说明的是,虽然名字是神经网络,并且结构等等也是借鉴了神经网络,但其原型以及算法本质上还和生物层面的神经网络运行原理存在

Python应用开发——30天学习Streamlit Python包进行APP的构建(9)

st.area_chart 显示区域图。 这是围绕 st.altair_chart 的语法糖。主要区别在于该命令使用数据自身的列和指数来计算图表的 Altair 规格。因此,在许多 "只需绘制此图 "的情况下,该命令更易于使用,但可定制性较差。 如果 st.area_chart 无法正确猜测数据规格,请尝试使用 st.altair_chart 指定所需的图表。 Function signa