5.6 奇异值

news2024/11/18 19:39:49

对单位圆的乘法

  首先我们在单位圆上遍历所有的点,作为二维向量,来研究某个矩阵乘以这些向量得到的结果,我们选三种矩阵,秩为0的矩阵,秩为1的矩阵和秩为2的矩阵。
  秩为0的矩阵就一个,也就是0矩阵,毫无疑问,它把所有向量变成0向量,没什么可以多说的了。
  我们任选一个秩为1的矩阵,就选这个吧:
( 1 1 1 1 ) \begin{pmatrix} 1 & 1\\ 1 & 1 \end{pmatrix} (1111)
  然后我们再遍历单位圆上的所有向量,毫无以为可以用以下形式表示任意单位向量:
( cos ⁡ θ sin ⁡ θ   ) \begin{pmatrix} \cos \theta \\ \sin \theta \ \end{pmatrix} (cosθsinθ )
  那么矩阵乘以这些单位向量就变成了以下向量:
( cos ⁡ θ + sin ⁡ θ cos ⁡ θ + sin ⁡ θ ) \begin{pmatrix} \cos \theta+\sin \theta \\ \cos \theta+\sin \theta \end{pmatrix} (cosθ+sinθcosθ+sinθ)
  最终它把所有单位向量都变到这根直线上:
在这里插入图片描述
  再找个秩为2的矩阵,就选个简单点的:
( 1 1 1 0 ) \begin{pmatrix} 1 & 1\\ 1 & 0 \end{pmatrix} (1110)
  那么矩阵乘以这些单位向量就变成了以下向量:
( cos ⁡ θ + sin ⁡ θ cos ⁡ θ ) \begin{pmatrix} \cos \theta+\sin \theta \\ \cos \theta \end{pmatrix} (cosθ+sinθcosθ)
  它把单位向量变到这个斜放的椭圆上:
在这里插入图片描述

矩阵2-范数的定义

  回忆下矩阵2-范数的定义:
∥ A ∥ 2 = max ⁡ ∥ A x ∥ 2 ∥ x ∥ 2 \parallel A \parallel_2=\max \frac{\parallel Ax \parallel_2}{\parallel x \parallel_2} A2=maxx2Ax2
  现在这些向量都是单位向量,2-范数都是1。所以上式k可以变成:
∥ A ∥ 2 = max ⁡ ∥ A x ∥ 2 \parallel A \parallel_2=\max \parallel Ax \parallel_2 A2=maxAx2
  那就是在上图的绿色图形中找到离原点最远的那个点的长度就是2-范数了。以第一个秩为1的矩阵为例子,它绿色曲线上2-范数最大的向量是:
( ± 2 ± 2 ) \begin{pmatrix} \pm \sqrt2 \\ \pm \sqrt2 \end{pmatrix} (±2 ±2 )
  所以它的2-范数是2.
  下面这个矩阵呢?它那个绿色曲线最圆的点是什么呢?直接硬算是不好算的,因为是斜的,所以需要把坐标变一下,这根绿色的曲线,按 y = 5 − 1 2 x y=\cfrac{\sqrt5-1}2x y=25 1x这根直线轴对称。如图:
在这里插入图片描述
  解以下方程组:
x 2 − 2 x y + 2 y 2 = 1 y = 5 − 1 2 x ⇒ x 2 − 2 x 5 − 1 2 x + 2 ( 5 − 1 2 x ) 2 = 1 ⇒ x 2 − ( 5 − 1 ) x 2 + 2 6 − 2 5 4 x 2 = 1 ⇒ x 2 − ( 5 − 1 ) x 2 + ( 3 − 5 ) x 2 = 1 ⇒ 5 x 2 − 2 5 x 2 = 1 ⇒ x 2 = 5 + 2 5 5 x = ± 5 + 2 5 5 y 2 = 3 − 5 2 × 5 + 2 5 5 y = ± 5 + 5 10 x^2-2xy+2y^2=1\\ y=\cfrac{\sqrt5-1}2x\\ \Rightarrow x^2-2x\cfrac{\sqrt5-1}2x+2(\cfrac{\sqrt5-1}2x)^2=1\\ \Rightarrow x^2-(\sqrt5-1)x^2+2\cfrac{6-2\sqrt5}4x^2=1\\ \Rightarrow x^2-(\sqrt5-1)x^2+(3-\sqrt5)x^2=1\\ \Rightarrow 5x^2-2\sqrt5x^2=1\\ \Rightarrow x^2=\frac{5+2\sqrt5}{5}\\ x=\pm \sqrt{\frac{5+2\sqrt5}{5}}\\ y^2=\frac{3-\sqrt5}2\times\frac{5+2\sqrt5}{5}\\ y=\pm \sqrt{\frac{5+\sqrt5}{10}} x22xy+2y2=1y=25 1xx22x25 1x+2(25 1x)2=1x2(5 1)x2+24625 x2=1x2(5 1)x2+(35 )x2=15x225 x2=1x2=55+25 x=±55+25 y2=235 ×55+25 y=±105+5
  所以绿色曲线上2-范数最大的向量是:
( ± 5 + 2 5 5 ± 5 + 5 10 ) \begin{pmatrix} \pm \sqrt{\frac{5+2\sqrt5}{5}} \\ \pm \sqrt{\frac{5+\sqrt5}{10}} \end{pmatrix} ±55+25 ±105+5
  那么它的二范数就是:
x 2 + y 2 = 5 + 2 5 5 + 5 + 5 10 = 15 + 5 5 10 = 3 + 5 2 x^2+y^2=\frac{5+2\sqrt5}{5}+\frac{5+\sqrt5}{10}\\=\frac{15+5\sqrt5}{10}=\frac{3+\sqrt5}{2} x2+y2=55+25 +105+5 =1015+55 =23+5

奇异值

  好了,但是这样求2-范数也太麻烦了吧,要解二元二次方程组,还要用到微积分或代数几何知识?其实大可不必。奇异值就是用来做这个事情的,奇异值的定义是 A H A A^HA AHA的所有特征值。以上例中的秩为2的矩阵来说,它有两个奇异值:
σ 1 = 3 + 5 2 , σ 2 = 3 − 5 2 \sigma_1=\frac{3+\sqrt5}{2},\sigma_2=\frac{3-\sqrt5}{2} σ1=23+5 ,σ2=235
  这恰好是对单位圆拉伸的最大与最小长度,而且这两个拉伸方向是成正交的。对于高维的矩阵,就分别是球体的 n n n个方向的拉伸倍数,这几个方向中有最大拉伸倍数也有最小拉伸倍数。所以通过二维的变换就可以更容易理解奇异值。

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

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

相关文章

【页面设计】02.CSS

CSS&#xff1a;Cascading Style Sheets 层叠样式表1.不是编程语言。2.告诉浏览器如何指定样式、布局等一、基本格式效果&#xff1a;二、三种方式添加CSS1.外部样式表&#xff08;1&#xff09;CSS保存在.css文件中&#xff08;2&#xff09;在HTML的<head>中使用<li…

棕黑色L-CQDs/TiO2 离子液体修饰/500nm粒径氮硫掺杂碳量子点N,S-CQDs/ZnO的制备过程

棕黑色L-CQDs/TiO2 离子液体修饰/500nm粒径氮硫掺杂碳量子点N,S-CQDs/ZnO的制备过程 今天小编分享L-CQDs/TiO2 离子液体修饰碳量子点的制备过程&#xff0c;一起看看吧: L-CQDs/TiO2 离子液体修饰碳量子点的制备过程&#xff1a; 通过水热法制备L-CQDs/TiO2复合催化剂。将50 …

Harbor安装对接Containerd

使用docker-compose安装harbor 先决条件&#xff1a; 安装docker安装docker-compose 安装参考 下载并解压 wget -c https://github.com/goharbor/harbor/releases/download/v2.3.5/harbor-offline-installer-v2.3.5.tgztar -zxvf harbor-offline-installer-v2.3.5.tgz cd harbo…

Aurora、Chip2chip、Ethernet(一)

摘要&#xff1a;之前的文章对aurora、chip2chip以及Ethernet这三个IP都进行介绍、仿真和使用说明。但是在实际使用中一定没有那么简单&#xff0c;在复杂联合使用的情况下&#xff0c;肯定会碰到各种各样的问题。此系列文章主要说明如何解决联合使用情况下碰到的一系列问题。 …

操作系统(day02)

“指令”就是处理器&#xff08;CPU&#xff09;能识别、执行的最基本的命令也可以叫做机器指令 两种指令、两种处理器状态、两种程序 两种指令 特权指令 如内存清零指令&#xff0c;不允许用户程序使用非特权指令 如普通的运算指令 既然有两种指令&#xff0c;且特权指令不…

IB学习者培养目标-知识渊博

“We explore concepts, ideas and issues that have local and global significance. In so doing, we acquire in-depth knowledge and develop understanding across a broad and balanced range of disciplines.” -IB definition of the attribute Knowledgeable“Being …

C++ —— 容器适配器和仿函数

目录 1.什么是容器适配器 2.stack的模拟实现 3.queue的模拟实现 4.deque概述 5.priority_queue的模拟实现 5.1仿函数 5.2模拟实现 6.反向迭代器 1.什么是容器适配器 在已有的容器(vector、list)的基础上适配出其他的容器。就类似于手机、笔记本电脑的电源适配器&…

一些lc周赛

6285. 执行 K 次操作后的最大分数(327 贪心 优先队列模拟) Math.ceil(val) 向上取整函数 public long maxKelements(int[] nums, int k) {PriorityQueue<Integer> queuenew PriorityQueue<>((a,b)->(b-a));for(int n:nums){queue.add(n);}long sum0;for(int i0…

JavaScript高级 ES6新特性

ES6~ES13新特性1. ECMA新描述概念1. 概念区别回顾2. 词法环境3. 环境记录4. 内存图的表示2. let、const的使用1. 基础的使用2. 作用域提升3. 暂时性死区 (TDZ)4. window 添加属性的区别5. 块级作用域的使用6. var、let、const的选择3. 模板字符串的详解4. ES6函数的增强用法1. …

VSCODE 系列(七)格式化工具clang-format

文章目录一、VS Code中使用生成.clang-format文件VS Code设置参考一、VS Code中使用 VS Code 中自带clang-format.exe 生成.clang-format文件 使用命令 .\clang-format.exe -stylellvm -dump-config > .clang-format或者新建.clang-format文件&#xff0c;将自己的配置…

vscode firefox xdebug 安装及配置

一、安装 vscode扩展中招xdebug直接安装。 firefox扩展中找xdebug直接安装。 xdebug下载&#xff0c;以window为例。 根据配置选下载内容。 设置成和ide相同的 。 二、配置 文档地址&#xff1a;Xdebug: Documentation 所有配置说明&#xff1a;Xdebug: Documentation All …

.mp4 文件转化成 .bag 文件并在 rviz 中显示

文章目录一、Python实现.mp4和.bag相互转化1、.mp4转.bag验证是否转换成功&#xff1a;使用 rosplay2、.bag转.mp4二、rviz 读取 *.bag 数据包并显示1、查看bag数据包的基本信息2、rviz 显示信息一、Python实现.mp4和.bag相互转化 1、.mp4转.bag # -*- coding: utf-8 -*- ##i…

Webpack 的 Chunk,想怎么分就怎么分

想必大家都用过 webpack&#xff0c;也或多或少了解它的原理&#xff0c;但是不知道大家有没有写过 Webpack 的插件呢&#xff1f; 今天我们就一起来写一个划分 Chunk 的 webpack 插件吧&#xff0c;写完后你会发现想怎么分 Chunk 都可以&#xff01; 首先我们简单了解下 web…

图像配准:基于 OpenCV 的高效实现

在这篇文章中&#xff0c;我将对图像配准进行一个简单概述&#xff0c;展示一个最小的 OpenCV 实现&#xff0c;并展示一个可以使配准过程更加高效的简单技巧。什么是图像配准图像配准被定义为将不同成像设备或传感器在不同时间和角度拍摄的两幅或多幅图像&#xff0c;或来自同…

什么牌子的护眼灯最好推荐?盘点口碑好的护眼灯品牌

护眼灯是目前大部分家庭都在使用的灯具之一&#xff0c;利用光源起到保护视力的效果&#xff0c;预防近视&#xff0c;可谓是现代生活中伟大的发明&#xff0c;今天由小编来列出优秀的护眼灯品牌&#xff0c;并详细的介绍&#xff0c;告诉大家哪个护眼灯品牌好。① 南卡护眼台…

【微信小程序-原生开发】实用教程07 - Grid 宫格导航,详情页,侧边导航(含自定义页面顶部导航文字)

开始前&#xff0c;请先完成成员页的开发&#xff0c;详见 【微信小程序-原生开发】实用教程 06-轮播图、分类页签 tab 、成员列表&#xff08;含Tdesign升级&#xff0c;切换调试基础库&#xff0c;设置全局样式&#xff0c;配置组件按需注入&#xff0c;添加图片素材&#x…

Canal快速入门

Canal 一、Canal 入门 1.1、什么是 Canal ​ 阿里巴巴 B2B 公司&#xff0c;因为业务的特性&#xff0c;卖家主要集中在国内&#xff0c;买家主要集中在国外&#xff0c;所以衍生出了同步杭州和美国异地机房的需求&#xff0c;从 2010 年开始&#xff0c;阿里系公司开始逐步…

PaddleSeg图像分割预测pyinstaller打包报错:No module named ‘framework_pb2‘,问题解决

报错 在使用PaddlePaddle的FastDeploy部署图像分割时&#xff0c;需要验证目标电脑环境&#xff0c;所以先将预测代码predict.py用pyinstaller打包来试试&#xff0c;指令&#xff1a; pyinstaller -D predict.py 打包完成&#xff0c;成功生成exe 运行时报错如下&#xff1a;…

项目五linux 内网完整渗透测试实例笔记

1.DDD4 靶场介绍本靶场存在三个 flag 把下载到的虚拟机环境导入到虚拟机&#xff0c;本靶场需要把网络环境配置好。1.1.网络示意图2. 信息收集2.1.主机发现sudo netdiscover -i eth0 -r 192.168.1.0/242.2.nmap 主机发现nmap -sn 192.168.1.0/242.3.masscan 端口探测sudo massc…

Python批量改文件名

对以下路径中的文件名批量修改。 一、读取指定路径中的文件名 #导入标准库 import os#读取文件名 filesDir "路径……" fileNameList os.listdir(filesDir)#输出路径中的所有文件 for filename in fileNameList:print(filename)二、正则表达式提取需要保留的部分 …