线性代数 --- 特征值与特征向量(下)

特征值与特征向量

Eigen Values & Eigen Vectors

Part III:如何求解特征向量与特征值

The Key Equation

\mathbf{det(A-\lambda I)=0}

对于一般矩阵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 _{2}=-1

这里我们顺便基于矩阵的行列式和迹来验证一下两个特征根:

Trace=0+0=0        \lambda _{1}+\lambda _{2}=-1+1=0

det(A)=0-1=-1        \lambda _{1}*\lambda _{2}=-1*1=-1

特例 I:多个重复/相同的特征值会引发特征向量的短缺

令矩阵A为: 

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

先求特征值:

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

这会得到两个重复的特征值:

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

带入其中一个特征值,并求特征向量:

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

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

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

令自由变量x1=1,得到:

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

         现在问题来了,如果你代入第二个特征向量(还是3),并重复上述操作,你是无法找到一个与第一个特征向量线性无关的另一个特征向量了。最多你只是令自由变量x1为一个其他数,然后得到一个特征向量,但最终你得到的这个特征向量与之前得到的那个是线性相关的。

        这个例子说明了一个矩阵的特征向量的线性相关性即取决于矩阵的维数也取决于是否存在重复的特征值。对于一个 2x2 的方阵,如果它有两个不同的特征值,那么就能得到两个线性无关的特征向量。但如果是两个重复的特征值,就只能得到一个特征向量或两个线性相关的特征向量。

特例 II:特征根为复数的情况

已知旋转矩阵Q,求Q的特征向量和特征值:

Q=\begin{bmatrix} 0 & -1\\ 1& 0 \end{bmatrix}

        这个旋转矩阵Q的作用是把任意向量旋转90°,这里我们会遇到问题。根据特征向量的意义,当矩阵Q作用于特征向量x后,向量Qx的方向应当与x相同。但这个要求明显与旋转矩阵Q的作用向左。换句话说,我们需要找出一个向量,他旋转90°仍在原处。

        其次,根据特征值的性质,我们同样也会发现类似的"Bug"。首先,这是一个2x2矩阵,因此有两个特征值\lambda _{1},\lambda _{2}。其次,我们还知道矩阵的迹等于两个特征值的和,矩阵的行列式等于两个特征值的乘积:

Trace=\lambda _{1}+\lambda _{2}=0

Det=\lambda _{1}\cdot \lambda _{2}=1

从第一个等式上说\lambda _{1},\lambda _{2}必定是一正一负,而如果这样的话他们的乘积必定是负数,这和第二个等式相违背。

现在,我们按照一般的方法去求解: 

det(Q-\lambda I)=\begin{vmatrix} -\lambda &-1 \\ 1& -\lambda \end{vmatrix}=\lambda ^{2}+1=0\Rightarrow \lambda ^{2}=-1

最终得到两个复根:

\begin{matrix} \lambda _{1}=i\\ \lambda _{2}=-i \end{matrix}

这两个复数的特征值正好满足和为0,且乘积为1。

带入第一个特征根i:

(Q-\lambda I)x=0\Rightarrow (Q-iI)x=0

Q-iI=\begin{bmatrix} -i &-1 \\ 1& -i \end{bmatrix}

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

令自由变量x2=1,得到x1=1/-i=i:

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

带入第二个特征根-i:

(Q-\lambda I)x=0\Rightarrow (Q+iI)x=0

Q+iI=\begin{bmatrix} i &-1 \\ 1& i \end{bmatrix}

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

令自由变量x2=1,得到x1=1/i=-i:

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


 (全文完)

作者 --- 松下J27 

 参考文献(鸣谢):

1,Introduction to Linear Algebra,Fifth Edition - Gilbert Strang

2,麻省理工Gilbert Strang教授线代大师-线性代数(全)_哔哩哔哩_bilibili

(配图与本文无关) 

版权声明:文中的部分图片,文字或者其他素材,可能来自很多不同的网站和说明,在此没法一一列出,如有侵权,请告知,立即删除。欢迎大家转载,但是,如果有人引用或者COPY我的文章,必须在你的文章中注明你所使用的图片或者文字来自于我的文章,否则,侵权必究。 ----松下J27

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

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

相关文章

Git win与linux换行符转换的问题

转载:warning: in the working copy of ‘package-lock.json‘, LF will be replaced by CRLF the next time Git_warning: in the working copy of package.json, lf-CSDN博客 warning: in the working copy of ‘package-lock.json‘, LF will be replaced by CRL…

“人工智能+”成国策,天翼云如何打造大模型时代的驱动力样板?

文 | 智能相对论 作者 | 叶远风 两会政府工作报告中明确指出要制定支持数字经济高质量发展政策,开展“人工智能”行动,这意味着,“人工智能”已经成为国策。 在这种大背景下,大模型的重要性无疑被抬到了新的高度。 此时&#…

数据结构和算法:栈与队列

栈 栈 &#xff08;stack&#xff09;是一种遵循先入后出逻辑的线性数据结构 把堆叠元素的顶部称为“栈顶”&#xff0c;底部称为“栈底”。 将把元素添加到栈顶的操作叫作“入栈”&#xff0c;删除栈顶元素的操作叫作“出栈”。 栈的常用操作 /* 初始化栈 */ stack<int&g…

小明的背包——01背包问题

经典版 题目链接&#xff1a;1.小明的背包1 - 蓝桥云课 (lanqiao.cn) 01背包问题中&#xff0c;每种物品只有两种状态&#xff0c;即拿或不拿。设状态dp[i][j]max(dp[i-1][j],dp[i-1][j-w]v)&#xff1b;如果不拿物品i&#xff0c;那么最大价值就是dp[i-1][j]&#xff0c;如果…

Excel第26享:模糊查找之Hlookup函数与通配符的嵌套

1、需求描述 如下图所示&#xff0c;现第一行有三个参考值&#xff1a;人S、羊E、猪3&#xff0c;在第三行有5个字&#xff1a;马、牛、人、羊、猪&#xff0c;每个字如果出现在第一行的三个参考值中&#xff0c;就返回该单元格的数值。如&#xff0c;人&#xff0c;就返回“人…

Windows 应用商城无法打开,解决方法

Windows 应用商城无法打开&#xff0c;也会导致有些软件无法使用&#xff0c;比如动态主题无法更新。 解决方&#xff1a; 1、按windowsR键&#xff0c;输入inetcpl.cpl。确定。 2、Internet属性&#xff0c;高级&#xff0c;勾选“使用TLS1.1",“使用TLS1.2"。确定…

docker镜像ssh服务

基于commit命令实现 首先我们是基于Ubuntu:18.04版本做ssh服务&#xff0c;拉取镜像 [rootmaster ~]# docker pull ubuntu:18.04 [rootmaster ~]# docker images REPOSITORY TAG IMAGE ID CREATED SIZE ubuntu 18.04 f9a80a5…

RabbitMQ 模拟实现【三】:存储设计

文章目录 数据库设计SQLite配置数据库实现 数据库关于哈希表等复杂类的存储启动数据库 文件设计消息持久化消息属性格式核心方法消息序列化消息文件回收 统一硬盘存储管理内存存储管理线程安全数据结构实现 数据库设计 数据库主要存储交换机、队列、绑定 SQLite 此处考虑的是…

GPT-4.5 Turbo意外曝光,最快明天发布?OpenAI终于要放大招了!

大家好&#xff0c;我是木易&#xff0c;一个持续关注AI领域的互联网技术产品经理&#xff0c;国内Top2本科&#xff0c;美国Top10 CS研究生&#xff0c;MBA。我坚信AI是普通人变强的“外挂”&#xff0c;所以创建了“AI信息Gap”这个公众号&#xff0c;专注于分享AI全维度知识…

基于springboot实现数据资产管理系统项目【项目源码+论文说明】计算机毕业设计

基于springboot实现数据资产管理系统演示 摘要 固定资产管理系统主要是完成对系统用户管理、资产信息管理、资产变更管理、资产用途管理、资产类别管理和资产增减管理。因为利用本系统管理员可以直接录入信息&#xff0c;修改信息&#xff0c;删除信息&#xff0c;并且若在录入…

使用docker-compose部署MySQL三主六从半同步集群(MMM架构)

文章目录 &#x1f50a;博主介绍&#x1f964;本文内容部署MySQL三主六从半同步集群一主二从同步集群规划需要安装docker和docker-compose命令形式安装安装docker安装docker-compose 宝塔面板形式安装 部署node1节点的master1docker-compose.yaml文件my.cnf文件授权启动 部署no…

ACM中Java输入输出

ACM中Java输入输出 最初写算法时&#xff0c;是用Scanner的。因为当时接触的测试数据基本都是以算法的复杂度为主&#xff0c;但是后面遇到大量的输入数据时。发现Scanner远远不能满足条件。下面列出几种常用的输入输出方式。(输出统一用printwriter&#xff0c;系统的system.o…

海豚调度系列之:认识海豚调度

海豚调度系列之&#xff1a;认识海豚调度 一、海豚调度二、特性三、建议配置四、名次解释 一、海豚调度 Apache DolphinScheduler 是一个分布式易扩展的可视化DAG工作流任务调度开源系统。适用于企业级场景&#xff0c;提供了一个可视化操作任务、工作流和全生命周期数据处理过…

Prompt Engineering(提示工程)

Prompt 工程简介 在近年来&#xff0c;大模型&#xff08;Large Model&#xff09;如GPT、BERT等在自然语言处理领域取得了巨大的成功。这些模型通过海量数据的训练&#xff0c;具备了强大的语言理解和生成能力。然而&#xff0c;要想充分发挥这些大模型的潜力&#xff0c;仅仅…

2024/3/14打卡k倍区间(8届蓝桥杯)——前缀和+优化***

题目 给定一个长度为 N 的数列&#xff0c;A1,A2,…AN&#xff0c;如果其中一段连续的子序列 Ai,Ai1,…Aj 之和是 K 的倍数&#xff0c;我们就称这个区间 [i,j] 是 K 倍区间。 你能求出数列中总共有多少个 K 倍区间吗&#xff1f; 输入格式 第一行包含两个整数 N 和 K。 以下 N…

O2OA(翱途)开发平台系统安全-用户登录IP限制

O2OA(翱途)开发平台[下称O2OA开发平台或者O2OA]支持对指定的用户设置可以连接的客户端计算机的IP地址&#xff0c;以避免用户在不安全的环境下访问系统。本篇主要介绍如何开启O2OA用户登录IP限制。 一、先决条件&#xff1a; 1、O2Server服务器正常运行&#xff0c;系统安装部…

0基础学习VR全景平台篇第145篇:图层控件功能

大家好&#xff0c;欢迎观看蛙色VR官方——后台使用系列课程&#xff01;这期&#xff0c;我们将为大家介绍如何使用图层控件功能。 一.如何使用图层控件功能&#xff1f; 进入作品编辑页面&#xff0c;点击左边的控件后就可以在右边进行相应设置。 二.图层控件有哪些功能&am…

跨境电商SaaS独立站的真面目...(网站建站)

跨境电商独立站自外贸交易开始&#xff0c;就一直存在&#xff0c;接触过电商的朋友应该都听过&#xff0c;但大部分人仅仅只是停留在听过的阶段&#xff0c;并没有真正的去了解它&#xff1b;独立站&#xff0c;顾名思义就是一个独立的网站&#xff0c;不依附任何平台&#xf…

基于Java+SpringBoot+vue+element实现物流管理系统

基于JavaSpringBootvueelement实现物流管理系统 博主介绍&#xff1a;多年java开发经验&#xff0c;专注Java开发、定制、远程、文档编写指导等,csdn特邀作者、专注于Java技术领域 ** 作者主页 央顺技术团队** 欢迎点赞 收藏 ⭐留言 文末获取源码联系方式 文章目录 基于JavaSpr…

OSCP靶场--Exfiltrated

OSCP靶场–Exfiltrated 考点(1.cms 站点地图插入php反弹shell 2. CVE-2021-4034提权 3.root定时任务提权[CVE-2021-22204]) 1.nmap扫描 ┌──(root㉿kali)-[~/Desktop] └─# nmap -sV -sC -p- 192.168.155.163 --min-rate 2500 Starting Nmap 7.92 ( https://nmap.org ) a…
最新文章