矩阵奇异值

news2024/9/30 13:49:23

一、ATA

任给一个矩阵A,都有: ATA 为一个对称矩阵

例子:A为一个m×n的矩阵,A的转置为一个n×m的矩阵

对称矩阵的重要性质如下

① 对称矩阵的特征值全为实数(实数特征根)

② 任意一个n阶对称矩阵,一定有代数重数=几何重数

若对称矩阵A的某个特征值的重数=k,则对应的几何重数=k
即一个特征根的重数若为k,则:
该特征根可找到对应线性无关的特征向量个数一定也为k

③ 任意一个n阶对称矩阵,一定有n个线性无关的特征值向量

如果无重根,即n个特征值对应n个线性无关的特征向量
倘若有k重根,也该重根必有k个线性无关的特征向量
倘若有k个重根,可理解为k个相同的特征值而已(总的特征值个数还是n)
也就是:一定有n个特征值,一定有n个线性无关的特征值向量

④ 任意一个n阶对称矩阵,一定能相似对角化

任意一个n阶对称矩阵,一定有n个线性无关的特征值向量
所以对称矩阵一定满足对角化条件,即:
任意一个n阶对称矩阵A,都有 A=PDP−1
任意一个n阶对称矩阵A,都有相似矩阵为对角矩阵的矩阵
对角矩阵D上,主对角线上的元素即为A或D的特征值

⑤ 任意一个n阶对称矩阵,一定能正交对角化

不同特征值对应的特征向量是彼此正交的
彼此正交是指:两个向量的乘积=0,即两个向量相互垂直
也就是: A=PDP−1 中的 P 为正交矩阵

二、ATA 的特征值与奇异值

由于 ATA 相乘后得到的一定是一个对称矩阵(n×n)

那么 ATA 一定有:

n个实数特征值,记为: λ1,λ2,...,λn

n个线性无关的特征向量,记为: v1→,v2→,...,vn→

也就是: v1→,v2→,...,vn→ 为n个相互垂直的特征向量

那么此时称: λi 为 ATA 的奇异值

记作: σi=λi

三、ATA 奇异值的非负性

结论:ATA的特征值 λi≥0

从而可得:ATA的奇异值 λi≥0

即: σi≥0

推导过程如下:

矩阵A左乘特征向量 vi→ 得到一个新的向量

那么这个新的向量模的平方记为: ||Avi→||2

因为两个相同的向量点乘就等于该向量模的平方

a→·a→=||a→||×||a→||×cos(0)=||a→||×||a→||×1=||a→||2

那么: ||Avi→||2=(Avi→)(Avi→)

将两个向量的点乘转为两个矩阵相乘,即为:

(Avi→)(Avi→)=(Avi)T(Avi)

将转置打开得: (Avi)T(Avi)=viTATAvi

由矩阵乘法结合律可得: viTATAvi=viT(ATAvi)

特征值的定义为: Axi→=λxi→

因为 ATA 的特征值为: λi

于是: ATAvi=λivi

从而: viT(ATAvi)=viTλivi

因为: λi 为一个常数,满足矩阵乘法交换律

于是: viTλivi=λiviTvi

接着: viTvi 转为两个向量的点乘 vi→·vi→=||vi||2

那么: λiviTvi=λi||vi||2

此时我们将:特征向量进行归一化处理

也就是将:v1→,v2→,...,vn→ 转为n个相互垂直的单位特征向量

记为: u1→,u2→,..,un→

那么将 vi 替换成 ui 后为: λi||ui||2

因为: ui 为相互垂直的单位特征向量

所以:每个 ui 的模长为1,于是: ||ui||2=1

那么: λi||ui||2=λi

最终: ||Avi→||2=λi

由于: ||Avi→||2 为非负数(一个向量的模长的平方为非负数)

所以: λi 也为非负数,即 λi≥0

四、换个角度看ATA的奇异值

由上面可知: σi=λi=||Avi→||2=||Avi→||

那么可以说:奇异值 σi 就是这个向量 Avi→ 的模长

也可以说:矩阵A左乘 ATA 特征向量 vi 后

得到该新向量的模长 (||Avi→||) 即为 ATA 的奇异值

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

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

相关文章

思科dhcp的配置

以路由器为例 让pc3 自动获取ip地址并获取的网段为172.16.4.100-172.16.4.200 配置如下: R1(config)#interface GigabitEthernet0/2 R1(config)#ip address 172.16.4.254 255.255.255.0 R1(config)# no shutdown R1(config)#ip dhcp pool 4_pool //创建dhcp地址池…

实际有库存却提示可用量不足保存不了杂发单

财务要统计研发费用,成本的金额。研发人员没有足够的意识配合。开立请购单时兴之所致,任性自由。想弄一个项目号就弄一个。不开心就没有项目号啦。哪管他人死活。 U9的逻辑,请购单如果带入项目号(客制化的功能)&#x…

c语言200例 067

大家好,欢迎来到无限大的频道 今天给大家带来的是c语言200例 题目要求: 设计一个共用体类型,使其成员包含多种数据类型,根据不同的数据类型,输出不同的结果 要设计一个共用体(union)类型&…

如何判断主机字节序

测试代码: #include <stdio.h> void byteorder() {union{short value;char union_bytes[sizeof(short)];//union_bytes数组}test;test.value 0x0102;if((test.union_bytes[0] 1) && (test.union_bytes[1]2)){printf("big endian\n");}else if((test…

初识Java反序列化漏洞

目录 为什么需要序列化&#xff1f; 序列化与反序列化基础案例 Serializable 接口 序列化对象 反序列化对象 Java 反序列化漏洞 readObject() 序列化&#xff1a;将对象的状态信息转换为可以存储或传输的形式的过程&#xff0c;即将对象转换为字节序列。反序列化&#x…

如何突破科技服务领域的客户管理困境?

在知识产权与科技服务领域&#xff0c;企业面临着独特的客户管理需求和挑战&#xff0c;这些挑战不仅要求高度的专业性和精细化操作&#xff0c;还涉及复杂的法律流程、数据保密性以及不断变化的客户需求。传统的客户管理方式&#xff0c;如纸质档案、简单的电子表格或人工处理…

[教程]Crystal源码下载及编译

描述&#xff1a; 随着 Crystal Source 代码的更新&#xff0c;用于构建源代码和编译它们的指南已经过时&#xff0c;这导致了很多混淆和寻求帮助。 本指南将是一个完整的分步指南&#xff0c;从下载 Visual Studio 到启动到您的服务器。 此外&#xff0c;请确保下载此存储库中…

如何使用 Python 读取数据量庞大的 excel 文件

使用 pandas.read_excel 读取大文件时&#xff0c;的确会遇到性能瓶颈&#xff0c;特别是对于10万行20列这种规模的 .xlsx 文件&#xff0c;常规的 pandas 方法可能会比较慢。 要提高读取速度&#xff0c;关键是找到更高效的方式处理 Excel 文件&#xff0c;特别是在 Python 的…

Docker版MKVtoolnix的安装及中文显示

本文是应网友 kkkhi 要求折腾的&#xff0c;只研究了 MKVtoolnix 的安装及中文显示&#xff0c;未涉及到软件的使用&#xff1b; 什么是 MKVtoolnix &#xff1f; MKVToolnix 是一款功能强大的多媒体处理工具&#xff0c;用于在 Linux、其他 Unix 系统和 Windows 上创建、修改和…

uniapp/vue项目 import 导入文件时提示Module is not installed,‘@/views/xxx‘路径无法追踪

文章目录 背景解决方案1.IDE配置2.alias&#xff08;别名&#xff09;配置webpackvue-clivite 3.检查 jsconfig.json 或 tsconfig.json 写在最后 前往闪闪の小窝以获得更好的阅读和评论体验 背景 Vue3在我自学Vue的时候看过一点&#xff0c;实操过一点&#xff0c;但是太久没用…

css 下拉框展示:当hover的时候展示下拉框 z-index的用法解释

代码如下&#xff1a; <template><div class"outer"><div class"left"></div><div class"aTest2"><div class"box">显示方框</div><div class"aTest3"></div></…

2024年度全球前2%顶尖科学家榜单发布!大陆10687人上榜,占比仅4.92%,附大陆Top 1000榜单!

9月16日&#xff0c;美国斯坦福大学的John Ioannidis教授团队新近在爱思唯尔(Elsevier)出版集团官网上发布了2024年全球各学科前2%科学家榜单(第七版)&#xff0c;通过统计每个科学家发表的论文数及其引用次数&#xff0c;采用6种指标建模进行综合评分(模型和指标附后)&#xf…

MIC麦克风工作原理:【图片+公式计算讲解】

Mic一个声电传换装置&#xff0c;广泛应用于电子产品中&#xff0c;最常见的就是手机。我们通常也叫它送话器&#xff0c;今天我们就来好好学习一下Mic的电路和它的工作原理。 1&#xff1a;电路结构图 1&#xff1a;C&#xff1a;Mic的核心部件&#xff0c;是一个可以膜片震动…

Linux之Docker虚拟化部署

上传docker安装包 解压安装包 将解压后的docker文件夹移动到/usr/local/文件夹下 docker 启动命令/usr/local/docker/dockerd 但是启动报错&#xff0c;意思是没有docker用户组 创建docker用户组&#xff0c;执行完会生成套接字文件 将套接字文件加入该用户组管理 第二个错误原…

【Lcode 随笔】C语言版看了不后悔系列持续更新中。。。

文章目录 题目一&#xff1a;最长回文子串题目描述&#xff1a;示例输入与输出&#xff1a;题目分析&#xff1a;解题思路&#xff1a;示例代码&#xff1a;深入剖析&#xff1a; 题目二&#xff1a;合并K个有序链表题目描述&#xff1a;示例输入与输出&#xff1a;题目分析&am…

k8s中,ingress的实现原理,及其架构。

图片来源&#xff1a;自己画的 图片来源&#xff1a;k8s官网 首先&#xff0c;什么是ingress? 是服务还是控制器&#xff1f; 都不精确 ingress是一个api资源 service和deployment也是api资源。 这几个相互协作&#xff0c;组建成一个对外提供服务的架构。 ingress提供的…

java学习-idea编辑器基础使用设置

首先打开电脑中的idea编辑器&#xff0c;点击头部&#xff1a;File按钮 → Settings… 打开设置界面&#xff1b; 设置idea的主题 设置idea代码注释的字体颜色 设置idea编辑器的字体和字体大小 设置idea通过提示回车自动导入包 设置idea输入忽略大小写进行提示

计算机网络面试题——第二篇

1. TCP拆包和粘包 现象 粘包&#xff1a;指在TCP传输中&#xff0c;发送方的多个数据包在接收方被合并在一个包接收&#xff0c;导致多条消息数据粘在一起&#xff0c;接收方无法正确区分这些消息的边界。拆包&#xff1a;指的是发送方的一个数据包在接收方被分成了多个包接收…

CHARLS数据库系列教程(4)--多模型效应分析、Per SD、P for trend及限制立方样条图绘制

CHARLS 是一项具备中国大陆 45 岁及以上人群代表性的追踪调查&#xff0c;旨在建设一个高质量的公共微观数据库&#xff0c;采集的信息涵盖社会经济状况和健康状况等多维度的信息&#xff0c;以满足老龄科学研究的需要。 为利用国际上最佳的数据采集方式&#xff0c;并确保研究…

shinyproxy部署streamlit记录

shinyproxy部署streamlit记录 streamlit 也是构建shinyproxy的后段是docker&#xff0c;所以不但可以部署shiny应用&#xff0c;还可以部署streamlit应用以及其它的应用程序。 部署的过程分两步&#xff0c;第一步是构建streamlit应用的docker镜像&#xff0c;就是构建Docker…