若爾當型,Jordon Form

文章目錄

    • 一、相似
    • 二、若爾當型
      • 1.1 認識若爾當型
      • 1.2 凱萊-哈密頓定理 (Cayley-Hamilton Theorem)

一、相似

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

任意一個矩陣C,滿足 C=B?1ABC = B^{-1}ABC=B?1AB 都和A有著相同的特征值。C 和 A 是相似的。

相似具有:

  • 自反性
  • 對稱性
  • 傳遞性

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

所有滿足 A=B?1CBA = B^{-1}CBA=B?1CB 的矩陣都是相似的,他們的特征值相同。

證明
如果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,BCB?1有特征值同樣為λ的特征向量Bx事實上,(BCB?1)(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?=vk?1?

  • 所有的特征向量都是廣義特征向量,但并非所有廣義特征向量都是特征向量。
  • 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?AB?1AB=?J1????Js???=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?=?λi??1?????1λ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)k?1)+rank((A?λI)k?1)?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=0n?ci?λi定理稱:p(A)=i=0n?ci?Ai=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(λ)可寫成則有?=PJP?1(若爾當型)=PJkP?1=i=0n?ci?Ai=i=0n?ci?PJiP?1=P(i=0n?ci?Ji)P?1=Pp(J)P?1=det(A?λI)=det(J?λI)(相似矩陣具有相同的特征方程):p(λ)=(λ?λ1?)a1?(λ?λ2?)a2??(λ?λk?)ak??,αk?為代數重數:p(Jk?)=(Jk??λ1?I)a1??(Jk??λk?I)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??λk?I)Jk??λk?I=?00?10????010??這是一個冪零矩陣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(λ)=λ2?6λ+5A2?6A+5I=05I=6A?A25A?1=6I?AA?1=56I?A???

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/diannao/95750.shtml
繁體地址,請注明出處:http://hk.pswp.cn/diannao/95750.shtml
英文地址,請注明出處:http://en.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…