若尔当型,Jordon Form

文章目录

    • 一、相似
    • 二、若尔当型
      • 1.1 认识若尔当型
      • 1.2 凯莱-哈密顿定理 (Cayley-Hamilton Theorem)

一、相似

Every matrix C=B−1ABC = B^{-1}ABC=B1AB has the same eigenvalues as A. These C’s are similar to A.

任意一个矩阵C,满足 C=B−1ABC = B^{-1}ABC=B1AB 都和A有着相同的特征值。C 和 A 是相似的。

相似具有:

  • 自反性
  • 对称性
  • 传递性

All the matrices A=B−1CBA = B^{-1}CBA=B1CB are “similar”,They all share the eigenvalues of C.

所有满足 A=B−1CBA = B^{-1}CBA=B1CB 的矩阵都是相似的,他们的特征值相同。

证明
如果Cx=λx,则BCB−1有特征值同样为λ的特征向量Bx事实上,(BCB−1)(Bx)=λ(Bx)同理可证A的特征值都是C的特征值\begin{align} & 如果 Cx = \lambda x, 则 BCB^{-1} 有特征值同样为\lambda 的特征向量 Bx \\ & 事实上,(BCB^{-1})(Bx) = \lambda (Bx) \\ & 同理可证 A 的特征值都是C的特征值 \end{align} 如果Cx=λx,BCB1有特征值同样为λ的特征向量Bx事实上,(BCB1)(Bx)=λ(Bx)同理可证A的特征值都是C的特征值
如果矩阵A 和 B 均可相似对角化,那么我们可以轻松判别二者是否相似。

那么不可对角化的情况呢?

二、若尔当型

1.1 认识若尔当型

(Jordan form) If A has s independent eigenvectors, it is similar to a matrix J that has s Jordan blocks J1,...,JsJ_1, ..., J_sJ1,...,Js on its diagonal. B contains “generalized eigenvectors”:

若尔当型(Jordan form):如果A 有s 个线性无关的特征向量,那么 A 相似于矩阵J,J的对角线上有 s 个 若尔当块J1,...,JsJ_1, ..., J_sJ1,...,Js,矩阵B包含了**“广义特征向量”**。

广义特征向量定义如下:

一个向量 vvv 是对应于特征值 λ\lambdaλ 的广义特征向量,如果它满足: (A−λI)kv=0(A - \lambda I)^k v = 0(AλI)kv=0 对于某个正整数 k≥1k \ge 1k1 成立。

k>1k>1k>1 时,我们得到了一系列向量,它们构成一个若尔当链 (Jordan Chain)

一个长度为 kkk 的若尔当链由向量 v_1,v_2,…,v_k{v\_1, v\_2, \dots, v\_k}v_1,v_2,,v_k 构成,它们满足: (A−λI)vk=vk−1(A - \lambda I)v_k = v_{k-1}(AλI)vk=vk1

  • 所有的特征向量都是广义特征向量,但并非所有广义特征向量都是特征向量。
  • k 是使该零空间维度达到最大并稳定的最小整数(也就是 λ 对应的最大若尔当块的阶数)。
  • 特征空间完全包含在广义特征空间之中。
  • 如何求解vkv_kvk =>只需找到$k => v_k = N(A - \lambda I)^{k} \setminus N(A - \lambda I)^{k - 1} $

A 的若尔当型:
Jordan form of AB−1AB=[J1⋱Js]=J\text { Jordan form of } A \quad B^{-1} A B=\left[\begin{array}{lll} J_{1} & & \\ & \ddots & \\ & & J_{s} \end{array}\right]=J  Jordan form of AB1AB=J1Js=J
每一个若尔当块 JiJ_iJi 主对角线上都是特征值 λi\lambda_iλi,主对角线上方恰有1个1.
Jordan block in JJi=[λi1⋅⋅⋅1λi]. \text { Jordan block in } J \quad J_{i}=\left[\begin{array}{cccc} \lambda_{i} & 1 & & \\ & \cdot & \cdot & \\ & & \cdot & 1 \\ & & & \lambda_{i} \end{array}\right] \text {. }  Jordan block in JJi=λi11λi
一个 k×kk \times kk×k 的若尔当块 J_k(λ)J\_k(\lambda)J_k(λ) 就对应一个长度为 kkk若尔当链。将这些若尔当链的所有向量集合在一起,就构成了若尔当基 (Jordan Basis),也就是相似变换矩阵 BBB 的列向量。

一些性质:

  • 每一个矩阵都相似于一个若尔当型。
  • 若尔当块的数目等于线性无关特征向量的个数(几何重数)

如何求解若尔当型

无需求解相似变换矩阵B

只需要知道每个特征值对应的若尔当块的数量大小

步骤如下

  1. 求解所有特征值,以及其代数重数(代数重数是指特征值作为特征方程根的重数)。
  2. 对每一个特征值 λi\lambda_iλi 单独分析,假设特征值 λ\lambdaλ 的代数重数是 mmm。这意味着与 λ\lambdaλ 相关的所有若尔当块的阶数之和必须为 mmm关键在于确定块的数量和各自的大小:
    1. 若尔当块的数量: 对应于特征值 λ\lambdaλ 的若尔当块的数量,等于该特征值的几何重数 (Geometric Multiplicity)(即线性无关特征向量的个数)
    2. k阶若尔当块的数量:rank((A−λI)k−1)+rank((A−λI)k−1)−2rank((A−λI)k)rank((A - \lambda I)^{k - 1}) + rank((A - \lambda I)^{k - 1}) - 2rank((A - \lambda I)^{k})rank((AλI)k1)+rank((AλI)k1)2rank((AλI)k)

至此,我们了解了若尔当型,并且得出 A 和 B 相似的充要条件是 A和B 有着相同的若尔当标准型

1.2 凯莱-哈密顿定理 (Cayley-Hamilton Theorem)

凯莱-哈密顿定理 (Cayley-Hamilton Theorem):每一个n阶方阵 A,都满足其自身的特征方程。
记p(λ)=det(A−λI)展开:p(λ)=∑i=0nciλi定理称:p(A)=∑i=0nciAi=0\begin{align} & 记 p(\lambda) = det(A - \lambda I) \\ & 展开:p(\lambda) = \sum_{i=0}^{n}c_i \lambda^i \\ & 定理称:p(A) = \sum_{i=0}^{n}c_i A^i = 0 \end{align} p(λ)=det(AλI)展开:p(λ)=i=0nciλi定理称:p(A)=i=0nciAi=0
证明
设A=PJP−1(若尔当型)Ak=PJkP−1p(A)=∑i=0nciAi=∑i=0nciPJiP−1=P(∑i=0nciJi)P−1=Pp(J)P−1由于,p(λ)=det(A−λI)=det(J−λI)(相似矩阵具有相同的特征方程)可写成:p(λ)=(λ−λ1)a1(λ−λ2)a2⋯(λ−λk)ak⋯,αk为代数重数则有:p(Jk)=(Jk−λ1I)a1⋯(Jk−λkI)ak⋯\begin{align} 设 A &= PJP^{-1}(若尔当型) \\ A^k &= PJ^kP^{-1} \\ p(A) &= \sum_{i=0}^{n}c_iA^i \\ &= \sum_{i=0}^{n}c_iPJ^{i}P^{-1} \\ &= P(\sum_{i=0}^{n} c_iJ^i)P^{-1} \\ &= Pp(J)P^{-1} \\ 由于,p(\lambda)&= det(A - \lambda I) = det(J - \lambda I)(相似矩阵具有相同的特征方程) \\ 可写成&:p(\lambda) = (\lambda - \lambda_1)^{a_1} (\lambda - \lambda_2)^{a_2} \cdots (\lambda - \lambda_k)^{a_k} \cdots, \alpha_k 为代数重数 \\ 则有&: p(J_k) = (J_k - \lambda_1 I)^{a_1} \cdots (J_k - \lambda_k I)^{a_k} \cdots \\ \end{align} AAkp(A)由于,p(λ)可写成则有=PJP1(若尔当型)=PJkP1=i=0nciAi=i=0nciPJiP1=P(i=0nciJi)P1=Pp(J)P1=det(AλI)=det(JλI)(相似矩阵具有相同的特征方程):p(λ)=(λλ1)a1(λλ2)a2(λλk)ak,αk为代数重数:p(Jk)=(Jkλ1I)a1(JkλkI)ak

我们来仔细分析其中的一个关键因子:(Jk−λkI)。Jk−λkI=(0100⋱⋱100)这是一个幂零矩阵Nm×m,满足Nm=0而m≤代数重数=α故p(Jk)=0,定理得证\begin{align} &我们来仔细分析其中的一个关键因子:(J_k - \lambda_k I)。\\ &J_k - \lambda_k I = \begin{pmatrix} 0 & 1 & & \mathbf{0} \\ & 0 & \ddots & \\ & & \ddots & 1 \\ \mathbf{0} & & & 0 \end{pmatrix}\\ & 这是一个幂零矩阵N_{m\times m}, 满足 N^{m} = 0 \\ & 而 m \le 代数重数 = \alpha \\ & 故p(J_k) = 0, 定理得证 \end{align} 我们来仔细分析其中的一个关键因子:(JkλkI)JkλkI=0010010这是一个幂零矩阵Nm×m,满足Nm=0m代数重数=αp(Jk)=0,定理得证
一个比较实用的应用:求解矩阵的逆
p(λ)=λ2−6λ+5A2−6A+5I=05I=6A−A25A−1=6I−AA−1=6I−A5\begin{align} & p(\lambda) = \lambda^2 -6\lambda + 5 \\ & A^2−6A+5I=0 \\ & 5I = 6A - A^2 \\ & 5A^{-1} = 6I - A \\ & A^{-1} = \frac{6I - A}{5} \end{align} p(λ)=λ26λ+5A26A+5I=05I=6AA25A1=6IAA1=56IA

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

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

相关文章

解决uni-app微信小程序编译报错:unexpected character `1`

问题原因在uni-app微信小程序开发中&#xff0c;当template模板中包含英文符号<或>时&#xff0c;微信小程序的编译器会将这些符号误认为是HTML标签的开闭符号&#xff0c;从而导致类似unexpected character 1的编译错误。错误示例<view class"plan-bmi">…

[Linux] RAID存储技术

目录 RAID实现方式 RAID 0 RAID 1 RAID 5 RAID 10 管理RAID0 创建RAID 查看RAID 格式化和挂载 删除RAID 管理RAID1 创建RAID 查看RAID 格式化和挂载 增加热备盘 模拟故障 删除故障磁盘 删除RAID 管理RAID5 创建RAID 查看RAID md5设备划分分区 RAID实现方…

程序设计|C语言教学——C语言基础4:进阶

一、预处理指令预处理指令在编译前执行&#xff0c;除了#include&#xff0c;还有以下常用指令&#xff1a;1. #define 宏定义无参宏&#xff1a;定义常量或代码片段&#xff0c;编译时直接替换&#xff08;无类型检查&#xff09;。#define PI 3.1415926 // 定义常量 #define…

数据结构之heap算法

文章目录前言1. heap结构概述2. push_heap3. pop_heap4. sort_heap5. make_heap前言 heap这种数据结构&#xff0c;允许用户以任何次序将任何数据放入该结构中&#xff0c;但是最后取出数据的时候一定是权值最高&#xff08;或者最低&#xff09;的元素。主要和实现有关&#x…

MCU 软件断点调试注意事项!!!

——为什么调试器会在运行中改我的Flash程序&#xff1f;调试单片机时&#xff0c;很多人都有这样的疑问&#xff1a;明明我在调试前刷进去的固件是好的&#xff0c;为什么加了一个断点之后&#xff0c;调试器居然去改了 Flash&#xff1f; 如果我拔掉调试器&#xff0c;这个固…

启发式合并 + 莫队 恋恋的心跳大冒险

题目来源&#xff1a;2025 Wuhan University of Technology Programming Contest 比赛链接&#xff1a;Dashboard - 2025 Wuhan University of Technology Programming Contest - Codeforces 题目大意&#xff1a; Solution&#xff1a; 首先肯定要预处理出以每个节点为起点…

JCTools 无锁并发队列基础:ConcurrentCircularArrayQueue

ConcurrentCircularArrayQueue ConcurrentCircularArrayQueue 是一个抽象类&#xff0c;它为基于数组的并发循环队列提供了基础功能。从其命名可以看出几个关键特性&#xff1a;​​Concurrent​​&#xff1a;常指无锁并发。​​Circular Array​​&#xff1a;内部使用循环数…

力扣(LeetCode) ——622. 设计循环队列(C语言)

题目&#xff1a;622. 设计循环队列示例1&#xff1a; MyCircularQueue circularQueue new MyCircularQueue(3); // 设置长度为 3 circularQueue.enQueue(1); // 返回 true circularQueue.enQueue(2); // 返回 true circularQueue.enQueue(3); // 返回 true circularQueue.…

在JVM跑JavaScript脚本 | Oracle GraalJS 简介与实践

这是2024年初的 GraalVM 系列博文&#xff0c;当时写了大纲&#xff0c;知道一年半后的现在才得以完成发布&#x1f604; 1、概述 实话说&#xff0c;标题的场景为小众需求&#xff0c;日常开发基本用不到&#xff0c;我是最近在做一个低代码轮子玩具 app-meta 需要实现 FaaS&…

基于 EC 数据与大模型技术实现天气预报:从数据到上线的全栈方法

1. 先校准“EC 数据”与“AI 预报”的语境 EC 数据家族(常用) IFS/HRES:确定性全球模式,水平分辨率约 9 km,常用预报范围 10 天; IFS/ENS:51 成员集合预报,提供 15 天概率信息; ERA5:再分析数据,小时级、0.25,可追溯至 1940 年,用作训练/评测黄金基准。 AI 预报…

迭代器模式及优化

迭代器模式&#xff08;Iterator Pattern&#xff09;是一种行为型设计模式&#xff0c;用于提供一种统一的方式遍历聚合对象&#xff08;如集合、容器&#xff09;中的元素&#xff0c;而无需暴露对象的内部实现细节。它将遍历逻辑与聚合对象分离&#xff0c;使得遍历操作可以…

纯Qt手撕gb28181协议/gb28181协议服务端/gb28181协议设备端/gb28181设备模拟器/gb28181虚拟监控设备

一、前言说明 搞完onvif设备模拟器&#xff0c;总想着把28181设备模拟也实现&#xff0c;因为之前已经花了大力气把28181平台软件端实现了&#xff0c;为了实现这个组件&#xff0c;头发掉了一大把&#xff0c;专门把国标文档看了好几遍&#xff0c;逐行阅读&#xff0c;针对需…

【渗透实战】无下载器环境(curl/wget)下玩转 Metasploit 自动利用

1. 背景与问题场景 在渗透测试或漏洞利用中&#xff0c;Metasploit&#xff08;MSF&#xff09;是业界最常用的框架之一。 其许多 RCE&#xff08;远程代码执行&#xff09;模块在落地 payload&#xff08;如 Meterpreter 或反弹 shell&#xff09;时&#xff0c;采用了 CMD St…

jd-hotkey探测热点key

对任意突发性的无法预先感知的热点数据&#xff0c;包括并不限于热点数据&#xff08;如突发大量请求同一个商品&#xff09;、热用户&#xff08;如恶意爬虫刷子&#xff09;、热接口&#xff08;突发海量请求同一个接口&#xff09;等&#xff0c;进行毫秒级精准探测到。然后…

C#WPF实战出真汁07--【系统设置】--菜品类型设置

1、菜品设置介绍 菜品设置跟餐桌设置的功能目的是相同的&#xff0c;包括了新增&#xff0c;删除&#xff0c;编辑&#xff0c;分页&#xff0c;查询&#xff0c;重置&#xff0c;全选&#xff0c;全消&#xff0c;列表功能&#xff0c;实现流程也是布局设计&#xff0c;后台逻…

aave v3 存款与借款利息的计算方式

本文只涉及到利率计算的数学原理&#xff0c;不作源码解析:存款首先我们假设小明在aave里面存了10000usdt&#xff0c;存的时候年化收益率是5%,那么半年后其存款的利息是多少呢?常规的计算方式如下:利息10000*5%*(存款的时长/一年的时长)这么做有什么问题呢&#xff1f;假设现…

Windows MCP.Net:基于.NET的Windows桌面自动化MCP服务器深度解析

&#x1f4cb; 目录 项目概述 技术架构深度解析 核心功能模块详解 代码实现分析 使用场景与实战案例 性能优化与最佳实践 扩展开发指南 总结与展望 项目概述 什么是Windows-MCP.Net&#xff1f; Windows MCP.Net是一个基于.NET 10.0开发的Windows桌面自动化MCP&…

Boost.Asio学习(7):Boost.Beast实现简易http服务器

namespace beast boost::beast;beast::flat_buffer是一个用于 Boost.Asio 和 Boost.Beast 网络读写的缓冲区实现。专为 一次性顺序读取 / 消费 场景设计&#xff0c;比 std::string 或 std::vector 高效&#xff0c;因为它是扁平内存结构&#xff08;contiguous memory&#x…

深入解析JVM内存区域划分:从理论到实践

Java虚拟机&#xff08;JVM&#xff09;是Java程序运行的核心环境&#xff0c;它负责管理内存分配、垃圾回收、字节码执行等关键任务。理解JVM的内存区域划分&#xff0c;对于优化Java应用性能、排查内存问题&#xff08;如OutOfMemoryError、StackOverflowError&#xff09;至…

滑窗|贪心|✅滚动数组

lc17.08pair按身高升序、相同时体重降序排序结果是找体重序列的最长递增子序列长度核心&#xff1a;转化为二维最长递增子序列问题求解vector<int> dp;for (auto& p : hw) {int w p.second;auto it lower_bound(dp.begin(), dp.end(), w);if (it dp.end()) {dp.pu…