14.1 矩阵幂级数

2023-11-02

矩阵的幂

  现在讨论下矩阵的n次方的问题,比如下面的矩阵:
A 1 = ( 2 1 − 1 1 2 − 1 − 1 − 1 2 ) A 2 = ( 6 5 − 5 5 6 − 5 − 5 − 5 6 ) A 3 = ( 22 21 − 21 21 22 − 21 − 21 − 21 22 )   A^ 1 = \begin{pmatrix}2 & 1 & -1\\ 1 & 2 & -1\\ -1 & -1 & 2\\ \end{pmatrix}\\ A^ 2 = \begin{pmatrix}6 & 5 & -5\\ 5 & 6 & -5\\ -5 & -5 & 6\\ \end{pmatrix}\\ A^ 3 = \begin{pmatrix}22 & 21 & -21\\ 21 & 22 & -21\\ -21 & -21 & 22\\ \end{pmatrix}\\\ A1= 211121112 A2= 655565556 A3= 222121212221212122  
  后面的我就不算下去了,越算越大。再看看这个矩阵:
A 1 = ( 0.1 0.2 0 0.5 0.5 0.5 0 0.2 0.5 ) A 2 = ( 0.11 0.12 0.1 0.3 0.45 0.5 0.1 0.2 0.35 ) A 3 = ( 0.071 0.102 0.11 0.255 0.385 0.475 0.11 0.19 0.275 ) A 4 = ( 0.058 0.087 0.106 0.218 0.338 0.43 0.106 0.172 0.232 ) ⋮ A ∞ = ( 0 0 0 0 0 0 0 0 0 ) A^ 1 = \begin{pmatrix}0.1 & 0.2 & 0\\ 0.5 & 0.5 & 0.5\\ 0 & 0.2 & 0.5\\ \end{pmatrix} \\ A^ 2 = \begin{pmatrix}0.11 & 0.12 & 0.1\\ 0.3 & 0.45 & 0.5\\ 0.1 & 0.2 & 0.35\\ \end{pmatrix} \\ A^ 3 = \begin{pmatrix}0.071 & 0.102 & 0.11\\ 0.255 & 0.385 & 0.475\\ 0.11 & 0.19 & 0.275\\ \end{pmatrix} \\ A^ 4 = \begin{pmatrix}0.058 & 0.087 & 0.106\\ 0.218 & 0.338 & 0.43\\ 0.106 & 0.172 & 0.232\\ \end{pmatrix} \\ \vdots\\ A^ {\infty} = \begin{pmatrix}0 & 0 & 0\\ 0 & 0 & 0\\ 0 & 0 & 0\\ \end{pmatrix} A1= 0.10.500.20.50.200.50.5 A2= 0.110.30.10.120.450.20.10.50.35 A3= 0.0710.2550.110.1020.3850.190.110.4750.275 A4= 0.0580.2180.1060.0870.3380.1720.1060.430.232 A= 000000000
  就可以看到,数字越来越小了。为什么有的矩阵的幂越来越大,有的却越来越小呢?这是由什么决定的呢?

矩阵幂的极限

  矩阵的幂的极限由什么确定?是由矩阵的谱半径 ρ ( A ) \rho(A) ρ(A)决定的。那么什么是矩阵的谱半径?谱半径的定义是矩阵所有的特征值的模长的最大值。那么为什么会是这样?如果学习了约当标准型就很容易理解了,对于任意一个矩阵,都可以求出它的约当标准型,也就是 A = P − 1 J P A=P^{-1}JP A=P1JP,那么就可以这样计算了:
lim ⁡ n → ∞ A n = lim ⁡ n → ∞ ( P − 1 J n P ) = P − 1 ( lim ⁡ n → ∞ J n ) P \lim_{n\to\infty}A^n=\lim_{n\to\infty}(P^{-1}J^{n}P)=P^{-1}(\lim_{n\to\infty}J^{n})P nlimAn=nlim(P1JnP)=P1(nlimJn)P
  而 lim ⁡ n → ∞ J n \lim_{n\to\infty}J^{n} limnJn这个极限就取决于约当标准型的对角线,如果都小于1,那么就存在极限且极限为0.约当标准型的对角线上就是矩阵的所有特征值。

谱半径与范数

  计算特征值是十分麻烦的,所以计算谱半径也是十分麻烦的。那么实际中有没有好的办法呢?可以利用矩阵的范数来判断。我们知道对于任意矩阵范数,有以下定理:
ρ ( A ) = inf ⁡ α ∥ A ∥ α \rho(A)=\inf_{\alpha}\parallel A\parallel_{\alpha} ρ(A)=αinfAα
  也就是说,所有矩阵的所有范数都要大于等于谱半径:
ρ ( A ) ≤ ∥ A ∥ α \rho(A)\le \parallel A\parallel_{\alpha} ρ(A)≤∥Aα
  所以实际应用中,用这个就可以快速地判断一个矩阵幂的极限是否存在。比如列和范数(1-范数),行和范数(无穷范数)就是非常好的判断方法。第一步肉眼扫一下,如果有模长和小于1的行或列,根据上述不等式,那么谱半径一定小于1.所以矩阵幂的极限就一定收敛。

矩阵幂级数

  矩阵幂级数和数学分析里的幂级数定义差不多,不过是把变量x换成了矩阵A而已。如果谱半径小于1,那么矩阵幂级数存在,并且有以下公式:
lim ⁡ n → ∞ ∑ i = 0 n A i = ( E − A ) − 1 \lim_{n\to\infty}\sum_{i=0}^nA^i=(E-A)^{-1} nlimi=0nAi=(EA)1
  证明我就不证明了。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

14.1 矩阵幂级数 的相关文章

  • 「第二篇」全国一等奖,经验帖。

    点击上方 大鱼机器人 选择 置顶 星标公众号 福利干货 第一时间送达 阅读文本大概需要 6 分钟 0 前言 本文作者 谢斌 曾经获得2017年控制题 板球控制系统 全国一等奖 他之前有写过几篇关于比赛的文章 大家可以点击阅读 全国一等奖 他
  • 如何使用html制作网页

    如何使用html制作网页 一 html简介 1 1概念 HTML即HyperText Mark up Language 意思是超文本标记语言 HTML不是一种编程语言 而是一种标记语言 超文本指的是超链接 标记指的是标签 是一种用来制作网页

随机推荐

  • 如何保护单例不被反射修改?

    public class Safety private static Safety instance new Safety private Safety if instance null throw new RuntimeException
  • hadoop生态系统的详细介绍-详细一点

    前提 日常喜欢看一些微信分享的好文 总结下来 可以作为过滤器吧 节约更多人的时间 在这里引用的是别人的文章 对原文的作者表示感谢 确实写的很好 hadoop生态系统的详细介绍 简介 Hadoop是一个开发和运行处理大规模数据的软件平台 是A
  • Kafka By the sea——kafka的使用场景

    文章目录 消息队列概述 消息队列应用场景 异步处理 应用解耦 流量削锋 日志处理 消息通讯 消息中间件示例 电商系统 日志收集系统 常用消息队列 ActiveMQ Kafka 消息队列概述 消息队列中间件是分布式系统中重要的组件 主要解决应
  • js怎么做延迟函数delay

    const delay ms gt new Promise resolve reject gt setTimeout resolve ms const getData status gt new Promise resolve reject
  • 评分算法_协同过滤推荐算法:评分预测准确性评估

    纸上得来终觉浅 绝知此事要躬行 宋 陆游 冬夜读书示子聿 对于评分预测常用的准确性评测指标是均方根误差RMSE和平均绝对误差MAE RMSE 均方根误差 对大的偏差更敏感 MAE 平均绝对值误差 注意 R 表示数据集合的长度 准确性指标计算
  • PADS Router VX2.7 操作界面以及常用设置

    打开方式 直接双击Router或者从layout中打开 打开Router 右击工具栏 选择自己想要使用的工具 项目浏览器 输出窗口 电子表格 导航窗口 都在右上角 标志工具栏中 如果不小心关掉 点击即可恢复 坐标以及单位 设置在 工具栏中的
  • 聊聊软件测试的那些事

    笔者入行软件测试行业也有两年左右的时间了 这两年中 在工作中也学习 积累了一些知识 但是每每谈及理论 又好像怎么也说不清一些东西的定义 其实很多人认为 知识学习了会用就可以 但软件测试的道路上 打好基础是很重要的 有些东西你知道但无法清晰表
  • 【leetcode每日刷题】214. 最短回文串

    给定一个字符串 s 你可以通过在字符串前面添加字符将其转换为回文串 找到并返回可以用这种方式转换的最短回文串 示例 1 输入 aacecaaa 输出 aaacecaaa 示例 2 输入 abcd 输出 dcbabcd 解法 KMP算法 假设
  • BUCK电路输入电容计算

    输入电容决定了输入电压的纹波 对于Buck变换器的输入端来说 输入电流是不连续的 在开关管导通的时候会有极大的阶跃电流 芯 片 BUCK控制器 时 间 2021 04 27 说 明 适用于稳态和动态负载 在Buck变换器的输入电压最小时 满
  • Qt中自定义结构体的使用

    Qt的自定义结构体Qt是不认识的 下面就直接列出使用方法 第一步 建议把所需的结构体放在一个单独头文件中 防止头文件相互包含 gg 而且还有条件编译的头自动生成 直接向工作添加C 头文件 自己把名字取好就行了 注意 这样会在 pro中 HE
  • 【JavaScript高级】内存管理与闭包:垃圾回收GC、闭包定义、访问和执行过程、内存泄漏

    文章目录 内存管理 垃圾回收GC 引用计数 Reference counting 标记清除 mark Sweep 闭包 定义 闭包的访问和执行过程 内存泄漏 浏览器的优化 参考 内存管理 所有编程语言 在代码的执行过程中都需要给它分配内存
  • 【python】多进程multiprocessing模块、进程池的使用

    multiprocessing中的多进程Process的基本使用 在python中 进程是通过 multiprocessing 多进程模块来管理的 multiprocessing模块提供了一个Process类来创建进程对象 创建子进程 Pr
  • 【使用 BERT 的问答系统】第 2 章 :用于自然语言处理的神经网络

    大家好 我是Sonhhxg 柒 希望你看完之后 能对你有所帮助 不足请指正 共同学习交流 个人主页 Sonhhxg 柒的博客 CSDN博客 欢迎各位 点赞 收藏 留言 系列专栏 机器学习 ML 自然语言处理 NLP 深度学习 DL fore
  • linux脚本定时调用存储过程,LINUX定时执行SHELL脚本实现DB2对存储过程的调用

    需求分析 本地化零件待办数量对应用户统计存入数据表 定时更新 使用linux的crontab定时任务来完成 1 编写存储过程 设置指向的数据库 SET SCHEMA DB2INST1 设置当前的路径 SET CURRENT PATH SYS
  • Java概述

    文章目录 一 Java简介 1 1 Java版本 1 2 Java特点 二 Java运行机制 2 1 Java运行过程 2 2 JDK JRE JVM 三 Java开发环境 3 1 下载 安装JDK 3 2 配置环境变量 四 Java开发规
  • R语言实战-第八章回归

    第八章 回归 简单线性回归 用到基础包中的women数据集 研究身高与体重的关系 head women fit lt lm weight height data women summary fit fitted fit 列出拟合模型的预测值
  • 【深度学习】RetinaFace人脸检测简要介绍

    介绍 Insight Face在2019年提出的最新人脸检测模型 原模型使用了deformable convolution和dense regression loss 当时在 WiderFace 数据集上达到SOTA 基网络有三种结构 基于
  • Java使用base64格式上传图片

    使用蚂蚁金服ui直接返回的是base64格式的图片 通过post方式进行请求 然后在控制器中以字符串的形式进行接收 接收之后进行转图片存储处理 只保存路径到数据库中 base64字节转图片代码 package com utils impor
  • qsort函数实现对任意数据的排序

    学会使用qsort函数排序 qsort介绍 compare函数介绍 不同的数据类型相应的比较函数定义 对数组元素为数字的 数组元素为字符时比较函数定义 结构体数据比较函数定义 qsort介绍 qsort函数是一个库函数 它的作用是对数据进行
  • 14.1 矩阵幂级数

    文章目录 矩阵的幂 矩阵幂的极限 谱半径与范数 矩阵幂级数 矩阵的幂 现在讨论下矩阵的n次方的问题 比如下面的矩阵 A 1