线性代数 --- 特征值与特征向量

特征值与特征向量

        已知任意向量x,现有矩阵A对x进行操作后,得到新的向量Ax。这就好比是自变量x与函数f(x)的关系一样,向量x通过类似“函数”的处理得到了一个新的向量Ax。这个新的向量可能和原向量x方向相同,也可能不同(事实上大多都不同)。此外,新的向量与原向量的长度可能向量,也可能不同。而特征向量(Eigen vector)指的就是那些和原始向量x平行的那些Ax,这是线性代数所研究的两大问题的的另一个部分(在我看来,线性代数的两个主要方向一个是研究垂直,另一个就是这里的平行)。

特征向量与特征值的意义:

        对于矩阵A而言,存在一些特殊的向量x,使得矩阵A作用于向量x之后所得到的Ax任然与原向量x平行(这里的平行可理解为与x的方向相同或正好相反),得到:

Ax=\lambda x

满足上述描述的向量x,就被称为特征向量。其中,系数\lambda是一个数,又叫特征值(Eigen value),可以看成是向量x长度的一个倍数。可以是1,表示Ax=x即,矩阵A没有改变向量x的长度。可以是-5,即改变了x的方向和大小,可以是0.26,也可以是复数。

Tips:当特征值\lambda=0,且x不是零向量时,矩阵A应该满足什么条件呢?

答:若\lambda=0,则\lambdax为零向量,得到Ax=0。向量x属于A的零空间(左零空间)又不是平凡解零向量,则矩阵A必须是奇异矩阵(singular matrix)才能让非零向量x通过线性组合得到零向量。这样一来,我们就得出了一条重要推导,若矩阵A为奇异矩阵,即不可逆矩阵,则矩阵A的特征值包含0。(这里我们再补充一点,对于矩阵A而言,可以有多个特征向量,对应多个特征值)。

特征向量与特征值的性质:

        1,n维度矩阵A有n个特征值。

        2,这n个特征值的乘积等于矩阵A的行列式。

        3,n个特征值的和等于矩阵A主对角线元素的和,这个和被称为Trace(迹)

一些常见矩阵的特征向量与特征值:

1,投影矩阵P(projection matrix)

        已知投影矩阵P可把任意向量投影到n维子空间S上,例如下图中,投影矩阵P作用于向量b上,得到了b在S上的投影Pb。

        现在我们要问的是,对于投影矩阵P而言,他作用在什么向量x上,得到的结果依然与x平行?如果平行,这个新向量的长度是x的几倍呢(即,求\lambda)?根据特征向量的定义,投影矩阵P作用于x上后得到的Px应当平行于Px。那么对投影矩阵P而言,特征向量x必须是n维子空间S内的向量。因为,对子空间S内的任意向量x施加P矩阵,得到投影就是x自己。方向与x相同,且长度不变,即特征值为1。用数学表达式来表示就是:

Px=x,\lambda =1

此外,我们还知道当向量b垂直于子空间S时,b在S上的投影为零向量。因此,我们又发现了S空间之外的特征向量,即所有垂直于S空间的向量x,投影Px与x的方向相同,但长度为0。表达为数学公式就是:

Px=0x\; or\; Px=0 ,\lambda =0

小结这里我们总结一下投影矩阵的特征向量与特征值,对于投影矩阵P而言,他的特征向量为投影矩阵P所投影的子空间内的所有向量和垂直于该子空间的所有向量。

2,置换矩阵(permutation matrix)

先看一个例子,已知置换矩阵A为(交换两个元素的位置):

A=\begin{bmatrix} 0 &1 \\ 1& 0 \end{bmatrix}

同样,还是按照特征向量的定义出发,现有的置换矩阵A作用在哪个向量x上,才能使得新向量Ax平行于x?换句话说,就是要找到一个向量x,交换元素位置后仍然是x或\lambda倍的x? 明显,如果向量x中的两个元素相同,则不论怎么交换两个元素的位置,得到的结果仍是x。

例如:

x=\begin{bmatrix} 1\\ 1 \end{bmatrix}\Rightarrow Ax=\begin{bmatrix} 0 &1 \\ 1& 0 \end{bmatrix}\begin{bmatrix} 1\\ 1 \end{bmatrix}=\begin{bmatrix} 1\\ 1 \end{bmatrix}=x\Rightarrow \lambda =1

因此,对于矩阵A而言,所有两个元素相同的向量x都是A的特征向量,且特征值为1。又因为,矩阵A为二维矩阵,因此,可能还存在一对特征向量与特征值。比如说,如果我们先令特征值为-1,那么就允许特征向量x中的元素大小相同但符号相反,这样一来,交换顺序再乘以-1后,得到的仍然还是x。

例如:

x=\begin{bmatrix} 1\\ -1 \end{bmatrix}\Rightarrow Ax=\begin{bmatrix} 0 &1 \\ 1& 0 \end{bmatrix}\begin{bmatrix} 1\\ -1 \end{bmatrix}=\begin{bmatrix} -1\\ 1 \end{bmatrix}=-x\Rightarrow \lambda =-1

 小结对于置换矩阵而言,那些交换元素位置后依然还是原始向量的向量就是A的特征向量,对应的特征值为1。而那些交换元素相应位置后,只改变了符号但不改变数值的向量,也是A的特征向量,他们所对应的特征值为-1。

对于一般矩阵A,如何找到他的特征值与特征向量?

Step I: Find λ first!

首先,我们有方程:

Ax=\lambda x

但这里有两个未知数,因此我们把上面的方程改写一下:

Ax=\lambda x\Rightarrow Ax=\lambda Ix\Rightarrow (A-\lambda I)x=0

这个齐次方程的解就是矩阵(A-\lambdaI)的零空间,抛开平凡解全0向量不说。要想让矩阵的零空间存在非零向量,则矩阵的A必为奇异矩阵,即不可逆矩阵。同时,结合之前学到的行列式的概念,若一个矩阵是奇异矩阵,则矩阵的行列式为0。这样一来,我们就不用考虑未知数x,也就是特征向量,先求未知数\lambda,也就是特征值。如下:

det(A-\lambda I)=0

这个方程是一个非常重要的方程(Key equation),叫特征值方程(Eigen-value equation)或者叫特征方程(characteristic equation)。

Step II: Substitue λ and solve equations!

 求解完特征值方程后会得到n个\lambda(可能会有相同的\lambda),把这些\lambda代入到经过改写后的方程组中,求解齐次方程组,或者说是求解零空间,得到相应的特征向量。

举个例子来说明上述求解过程:

A=\begin{bmatrix} 0 &1 \\ 1& 0 \end{bmatrix}

 对于上面提到过的置换矩阵A,第一步,我们先求det(A-λI)=0:

det(A-\lambda I)=\begin{vmatrix} -\lambda &1 \\ 1& -\lambda \end{vmatrix} =0

得到:

\lambda ^{2}-1=0\Rightarrow \lambda ^{2}=1

由此得到两个特征根:

\lambda _{1}=1, \lambda _{2}=-1

把第一个特征值λ=1代入改写后的方程:

(A-\lambda I)x=0\Rightarrow (A-I)x=0

其中(A-I)为:

A-I=\begin{bmatrix} -1 &1 \\ 1& -1 \end{bmatrix} 

        对该矩阵消元得到矩阵U,找到主元列和自由列,其中与自由列位置所对应的向量x中的元素就是自由变量。设自由变量为一个任意数,求出相应的特解。

U=\begin{bmatrix} -1 &1 \\ 0& 0 \end{bmatrix}\Rightarrow -x_{1}+x_{2}=0

消元矩阵U的秩为1,存在一个特解。第一列为主元列,第二列为自由列。因此,对应的自由变量为x2。令x2=1, 则x1=1得到一个特解,也就是我们要找的特征向量x:

x=\begin{bmatrix} 1\\ 1 \end{bmatrix},\lambda _{1}=1

相应的代入第二个特征值λ=-1:

(A-\lambda I)x=0\Rightarrow (A+I)x=0

A+I=\begin{bmatrix} 1 &1 \\ 1& 1 \end{bmatrix}

U=\begin{bmatrix} 1 &1 \\ 0& 0 \end{bmatrix}\Rightarrow x_{1}+x_{2}=0

令自由变量x2=1,则x1=-1:

x=\begin{bmatrix} -1\\ 1 \end{bmatrix},\lambda _{1}=-1

 

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

五、布局布线约束、系统优化参数、时序优化收敛 关键技术点

在实际的工程当中,出现了时序违例的情况如何解决呢? 本章内容将介绍例外约束、布局布线的具体操作,实现系统参数的优化。 **前言:**通过约束时钟,比如基准时钟,和生成时钟,让我们的综合工具知道我们的时序…

2024年阿里Android高级面试题分享,送给正在迷茫的你

前言 很多公司在招人这件事情上都会面临一个问题; “我们的招聘要求又不高,能做项目就行,但为什么就是招不到人?” 很多公司还面临一个问题,招聘的时候这人各方面都不错,但上岗了就是不出活,绩…

2024春招面试,文末有彩蛋

前言 Flutter 作为Google出品的一个新兴的跨平台移动客户端UI开发框架,正在被越来越多的开发者和组织使用,包括阿里的咸鱼、腾讯的微信等。 今天,我主要讲解Flutter中文本组件方面的Widget,包括Text、RichText、TextField&#…

元宇宙线上展厅——如何利用3D展示平台创新吸引客户与展示产品

在数字化转型的浪潮中,元宇宙线上展厅作为一种全新的虚拟展示平台,正以其独特的创新功能和沉浸式体验,吸引着越来越多的企业和用户的目光。 一、元宇宙线上展厅的创新功能 1、沉浸式体验与交互设计 元宇宙线上展厅通过先进的3D建模技术和VR技…

selenium4的相对定位

selenium4相对定位 Selenium 4新增了相对定位器,能帮助用户查找元素附近的其他元素。可用的相对定位器有above、below、toLeftOf、toRightOf、near。在Selenium 4中,find_element方法能够接受一个新方法withTagName,它将返回一个RelativeLoca…

《2023年DDoS攻击现状及趋势报告》

执行概要 与第一季度相比,2023年第二季度的DDoS攻击活动增长了387%。 电信公司遭受的攻击最为频繁,占总攻击量的50%,在2023年上半年发生了超过37,000次攻击;然而,在所有行业中,电信公司的攻击活动在同一时…

mac电脑总卡蓝屏是怎么回事,苹果电脑老卡蓝屏怎么办

电脑老卡蓝屏是比较常见的电脑故障之一,导致这一问题的出现很可能是电脑本身的硬件,或电脑上的驱动安装错误,没法运行,当然也不排除其他的一些因素。虽说电脑蓝屏是电脑几乎都会出现的小毛病,不足以致命,但…

Day19:信息打点-红蓝队自动化项目资产侦察武器库部署企查产权网络空间

目录 各类红蓝队优秀工具项目集合 自动化-武器库部署-F8x 自动化-网络空间-AsamF 自动化-企查信息-ENScan 自动化-综合架构-ARL&Nemo 思维导图 章节知识点 Web:语言/CMS/中间件/数据库/系统/WAF等 系统:操作系统/端口服务/网络环境/防火墙等 应…

一本书讲透ChatGPT,实现从理论到实践的跨越!大模型技术工程师必读

这里写目录标题 前言内容简介作者简介专家推荐读者对象目录直播预告 前言 OpenAI 在 2022 年 11 月推出了人工智能聊天应用—ChatGPT。它具有广泛的应用场景,在多项专业和学术基准测试中表现出的智力水平,不仅接近甚至有时超越了人类的平均水平。这使得 …

垃圾分类网站|基于Springboot框架+java+MYSQL数据库的垃圾分类网站开发设计与实现(可运行源码+数据库+文档)

目录 1.摘 要 2.系统结构设计 3.系统顺序图设计 4.数据库设计 5.系统详细设计 用户前台功能模块 管理员功能模块 垃圾分类管理员功能模块 论文参考 文末获取源码 1.摘 要 本论文主要论述了如何使用JAVA语言开发一个垃圾分类网站 ,本系统将严格按照软件开发…

探索直播美颜SDK背后的算法:如何实现高效的美颜处理?

直播中,美颜功能更是成为了吸引用户的一大利器,为了实现这一目标,各大直播平台纷纷推出了美颜功能,而直播美颜SDK背后的算法成为了实现这一功能的关键。 一、美颜算法的重要性 美颜算法在直播美颜SDK中扮演着至关重要的角色。它不…

isNaN和Number.isNaN()的区别

一句话概括: isNaN()会先尝试转换为数字,如果无法转换为数字则返回true,否则返回false Number.isNaN():直接检查一个值是否为NaN 示例如下: 对于isNaN() NaN直接就返回true "abc"是字符串且无法转换为数…

连锁经营如何降低财务成本和税务风险

连锁经营的财务是一个比较复杂的体系。连锁经营通过规模化运作,连锁企业可以享受采购、生产和销售方面的经济规模优势,从而降低采购成本、生产成本和运营成本。可以通过统一管理和监控各个门店的财务状况,实现资源的最优配置,减少…

如何单独设置cPanel的PHP扩展

我们在上周遇到购买Hostease的Linux虚拟主机客户网站页面需要使用mb_strlen函数。像这种需要特定PHP函数的设置需求,我们是可以单独在cPanel面板进行设置。 步骤 1:登录到 cPanel 打开您的 Web 浏览器,登录您的 cPanel 控制面板登录页面。 步…

Linux - 进程控制

1、进程创建 1.1、fork函数初识 在linux中fork函数时非常重要的函数&#xff0c;它从已存在进程中创建一个新进程。新进程为子进程&#xff0c;而原进程为父进程&#xff1b; #include <unistd.h> pid_t fork(void); 返回值&#xff1a;自进程中返回0&#xff0c;父进…

EC600模块通过AT指令接入阿里云物联网平台并发布属性

摘要&#xff1a;本文介绍一下如何通过EC600模块的AT指令&#xff0c;将设备属性值发送到阿里云物联网平台的方法。 这个模块供电可以是 5-16V 和电脑通过USB串口连接&#xff0c;4线即可。未来集成到自己的系统中的时候&#xff0c;可以直接发送指令即可。 使用的软件是FreeAT…

Claude3荣登榜首,亚马逊云科技为您提供先行体验!

Claude3荣登榜首&#xff0c;亚马逊云科技为您提供先行体验&#xff01; 个人简介前言抢先体验关于Amazon BedrockAmazon Bedrock 的功能 Claude3体验教程登录Amazon Bedrock试用体验管理权限详细操作步骤1.提交应用场景详细信息2.请求模型的访问权限3.请求成功&#xff0c;开始…

智慧合同管理平台:企业合同管理的数字化革新

在当今的商业环境中&#xff0c;合同管理已经成为企业运营中不可或缺的一环。传统的合同管理方式不仅耗时耗力&#xff0c;而且容易出错。为了解决这个问题&#xff0c;越来越多的企业开始转向使用智慧合同管理平台。本文将以道本科技智慧合同管理平台为例&#xff0c;详细介绍…

中航证券2024 AI智算时代已至,算力芯片加速升级

一、AI处史上最长繁荣期&#xff0c;算力国产化需求迫切 AI正处史上最长繁荣大周期 人工智能从1956 年被正式提出以来&#xff0c;经历了数十年的 发展历程。人工智能诞生初期&#xff0c;其研究主要分为三个流 派&#xff0c;即逻辑演绎、归纳统计和类脑计算。 人工智能研究的…

JasperReport指定自定义字体文件,解决中文不显示问题支持

原文链接&#xff1a;【JasperReports笔记05】JasperReport指定自定义字体文件&#xff0c;解决中文不显示问题支持Android_jasperreports导出pdf不显示中文-CSDN博客 目录 一、自定义字体文件 1.1、创建字体配置文件 1.2、创建fonts.xml字体文件 1.3、在Jasper Studio中添…