线性代数笔记23——矩阵的对角化和方幂

2024-02-14 20:58

本文主要是介绍线性代数笔记23——矩阵的对角化和方幂,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

特征值矩阵

  假设A有n个线性无关的特征向量x1,x2……xn,这些特征向量按列组成矩阵S,S称为特征向量矩阵。来看一下A乘以S会得到什么:

  最终得到了S和一个以特征值为对角线的对角矩阵的乘积,这个对角矩阵就是特征值矩阵,用Λ表示:

  没有人关心线性相关的特征向量,上式有意义的前提是S由n个线性无关的特征向量组成,这意味着S可逆,等式两侧可以同时左乘S-1:

  AS=SΛ和S-1AS=Λ就是对角化的两种方法。需要注意的是,并非所有矩阵A都存在n个线性无关的特征向量,这类矩阵不能对角化。

  矩阵对角化还有另一种表达:

  我们已经知道了矩阵的LU分解,A=LU;格拉姆-施密特正交化,A=QR;现在有多了一种对角化分解,A=SΛS-1

矩阵的方幂

  如果A存在特征值和特征向量,即Ax = λx,那么A2的特征值和特征向量是什么?

  这在上一章的示例中出现过,将Ax = λx的等式两侧同时左乘A就可以表示A的特征向量:

  由于λ是标量,所以可以把λ单独提出来:

  现在可以得出结论了,A2的特征向量不变,特征值变成了λ2

  可以用同样的方式看看A2的对角化:

  按照这个思路可以继续计算Ak的对角化,Ak的特征向量不变,Ak的特征值矩阵是A的特征值矩阵的k次方:

  根据上式,如果k→∞,在所有特征值|λi|<0时,Ak→0,当然,前提是A有n个线性无关的特征向量。

对角化的前提

  对角化的前提是A存在n个线性无关的特征向量,问题是怎样判断A存在n个线性无关的特征向量?一个判断方法是:当A的所有特征组互不相同时,A必然存在n个线性无关的特征向量;如果存在重复的特征值就不好说了,需要另行判断。

  n阶单位矩阵的所有特征值都是1,但是它仍然有n个线性无关的特征向量,因此单位矩阵可以对角化:

  再来看三角矩阵。三角矩阵A的各列是线性无关的,意味着它有唯一解,没有n个线性无关的特征向量,比如下面这个:

  先计算A的特征值:

  作为2×2矩阵,A只有一个特征向量,它无法完成对角化。

使用对角化

  给定一个向量u0和一个能够对角化的矩阵A,如果uk+1=Auk,那么u100 = ?

  可以简单的向后推导一下:

  现在可以得到结论,u100=A100 u0,问题是如何求得A100?

  A有n个线性无关的特征向量x1,x2,……,xn,这意味着u0可以看成这些特征向量的线性组合:

  以单位矩阵为例,假设A是3×3的单位矩阵,则A的三个特征向量是:

  这三个特征向量可以通过线性组合成为任意的三维向量。

  现在可以将Au0写成下面的形式:

  由于Ci是标量,所以可以将Ci写到前面:

  x1,x2,……,xn都是A的特征向量,它们以特征值为媒介和A存在关联,Axi = λixi,因此:

  等式两侧同时左乘A:

  同样,可以把比标量Ciλi放到前面:

  无论等式两侧再左乘几个A都将得到类似的结果,因此:

  这就是最终的答案,如果真要计算A100 u0,可以先把u0展开成特征向量的线性组合,求出具体的C值,在使用SΛ100C求解。

综合示例

  

  a,b都是0的时候没什么可算的,主要看ab≠0的情况。C看起来比较别扭,还是用A来说话。先来看一下特征值:

  特征值矩阵和特征向量矩阵是:

  当a=b=-1时:

 


   作者:我是8位的

  出处:http://www.cnblogs.com/bigmonkey

  本文以学习、研究和分享为主,如需转载,请联系本人,标明作者和出处,非商业用途! 

  扫描二维码关注公众号“我是8位的”

这篇关于线性代数笔记23——矩阵的对角化和方幂的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

安卓链接正常显示,ios#符被转义%23导致链接访问404

原因分析: url中含有特殊字符 中文未编码 都有可能导致URL转换失败,所以需要对url编码处理  如下: guard let allowUrl = webUrl.addingPercentEncoding(withAllowedCharacters: .urlQueryAllowed) else {return} 后面发现当url中有#号时,会被误伤转义为%23,导致链接无法访问

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

线性代数|机器学习-P36在图中找聚类

文章目录 1. 常见图结构2. 谱聚类 感觉后面几节课的内容跨越太大,需要补充太多的知识点,教授讲得内容跨越较大,一般一节课的内容是书本上的一章节内容,所以看视频比较吃力,需要先预习课本内容后才能够很好的理解教授讲解的知识点。 1. 常见图结构 假设我们有如下图结构: Adjacency Matrix:行和列表示的是节点的位置,A[i,j]表示的第 i 个节点和第 j 个

论文阅读笔记: Segment Anything

文章目录 Segment Anything摘要引言任务模型数据引擎数据集负责任的人工智能 Segment Anything Model图像编码器提示编码器mask解码器解决歧义损失和训练 Segment Anything 论文地址: https://arxiv.org/abs/2304.02643 代码地址:https://github.com/facebookresear

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

【C++学习笔记 20】C++中的智能指针

智能指针的功能 在上一篇笔记提到了在栈和堆上创建变量的区别,使用new关键字创建变量时,需要搭配delete关键字销毁变量。而智能指针的作用就是调用new分配内存时,不必自己去调用delete,甚至不用调用new。 智能指针实际上就是对原始指针的包装。 unique_ptr 最简单的智能指针,是一种作用域指针,意思是当指针超出该作用域时,会自动调用delete。它名为unique的原因是这个

查看提交历史 —— Git 学习笔记 11

查看提交历史 查看提交历史 不带任何选项的git log-p选项--stat 选项--pretty=oneline选项--pretty=format选项git log常用选项列表参考资料 在提交了若干更新,又或者克隆了某个项目之后,你也许想回顾下提交历史。 完成这个任务最简单而又有效的 工具是 git log 命令。 接下来的例子会用一个用于演示的 simplegit