矩陣的范數和逆矩陣的范數的關系

文章目錄

矩陣的范數和逆矩陣的范數之間的關系是線性代數中的一個重要主題,尤其是在數值分析和矩陣理論中。

前提條件

  • 我們討論的是方陣 A A A 的情況,因為只有方陣才可能存在逆矩陣。
  • 矩陣 A A A可逆的(即 det ? ( A ) ≠ 0 \det(A) \neq 0 det(A)=0),這意味著 A ? 1 A^{-1} A?1 存在。

關鍵結論

對于一個給定的矩陣范數 ∥ ? ∥ \|\cdot\| ?,如果 A A A 是可逆的,則有以下不等式成立:

∥ A ? 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A?1A1?

其中, ∥ A ∥ \|A\| A 表示矩陣 A A A 的范數,而 ∥ A ? 1 ∥ \|A^{-1}\| A?1 則表示其逆矩陣的范數。注意這里的范數可以是任何相容的矩陣范數,如誘導范數(induced norm)、Frobenius范數等。

推導過程簡述

考慮向量 x x x 和單位向量 y = A x / ∥ A x ∥ y = Ax/\|Ax\| y=Ax/∥Ax,我們有:

1 = ∥ y ∥ = ∥ A x ∥ A x ∥ ∥ ≤ ∥ A ∥ ? ∥ x ∥ A x ∥ ∥ 1 = \|y\| = \left\|\frac{Ax}{\|Ax\|}\right\| \leq \|A\| \cdot \left\|\frac{x}{\|Ax\|}\right\| 1=y= ?AxAx? ?A? ?Axx? ?

由此可以得出:

1 ∥ A ∥ ≤ ∥ x ∥ ∥ A x ∥ \frac{1}{\|A\|} \leq \frac{\|x\|}{\|Ax\|} A1?Axx?

對所有非零向量 x x x 成立。特別地,取 x = A ? 1 z x = A^{-1}z x=A?1z 對于任意非零向量 z z z,我們得到:

1 ∥ A ∥ ≤ ∥ A ? 1 z ∥ ∥ z ∥ \frac{1}{\|A\|} \leq \frac{\|A^{-1}z\|}{\|z\|} A1?zA?1z?

從而得出:

∥ A ? 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A?1A1?

注意事項

  • 這個關系表明了逆矩陣的范數至少為原矩陣范數的倒數。然而,在實際應用中, ∥ A ? 1 ∥ \|A^{-1}\| A?1 可能遠大于 1 / ∥ A ∥ 1/\|A\| 1/∥A,特別是在矩陣接近奇異的情況下(即當 det ? ( A ) \det(A) det(A) 非常接近于0時)。
  • 矩陣條件數(condition number)定義為 κ ( A ) = ∥ A ∥ ? ∥ A ? 1 ∥ \kappa(A) = \|A\| \cdot \|A^{-1}\| κ(A)=A?A?1,它衡量了矩陣求解線性系統時的敏感度。條件數越大,意味著矩陣越“病態”,數值計算的結果可能越不穩定。

示例說明

假設有一個 2 × 2 2 \times 2 2×2 矩陣 A A A

A = [ 4 7 2 5 ] A = \begin{bmatrix} 4 & 7 \\ 2 & 5 \end{bmatrix} A=[42?75?]

其逆矩陣為:

A ? 1 = [ 5 / 6 ? 7 / 6 ? 1 / 3 2 / 3 ] A^{-1} = \begin{bmatrix} 5/6 & -7/6 \\ -1/3 & 2/3 \end{bmatrix} A?1=[5/6?1/3??7/62/3?]

如果我們使用譜范數(spectral norm),即最大奇異值作為范數,則可以通過計算各自的奇異值來估計范數,并驗證上述不等式是否成立。

📘 譜范數定義

矩陣 A A A譜范數定義為它的最大奇異值:

∥ A ∥ 2 = σ max ? ( A ) \|A\|_2 = \sigma_{\max}(A) A2?=σmax?(A)

而奇異值是矩陣 A T A A^T A ATA 的特征值的平方根。因此,我們可以按以下步驟計算:

  1. 計算 A T A A^T A ATA
  2. 求其特征值
  3. 取最大特征值的平方根即為譜范數

? 步驟一:計算 A T A A^T A ATA

A T = [ 4 2 7 5 ] , A T A = [ 4 2 7 5 ] [ 4 7 2 5 ] = [ 4 2 + 2 2 4 ? 7 + 2 ? 5 7 ? 4 + 5 ? 2 7 2 + 5 2 ] = [ 20 38 38 74 ] A^T = \begin{bmatrix} 4 & 2 \\ 7 & 5 \end{bmatrix}, \quad A^T A = \begin{bmatrix} 4 & 2 \\ 7 & 5 \end{bmatrix} \begin{bmatrix} 4 & 7 \\ 2 & 5 \end{bmatrix} = \begin{bmatrix} 4^2 + 2^2 & 4\cdot7 + 2\cdot5 \\ 7\cdot4 + 5\cdot2 & 7^2 + 5^2 \end{bmatrix} = \begin{bmatrix} 20 & 38 \\ 38 & 74 \end{bmatrix} AT=[47?25?],ATA=[47?25?][42?75?]=[42+227?4+5?2?4?7+2?572+52?]=[2038?3874?]


? 步驟二:求 A T A A^T A ATA 的特征值

M = A T A = [ 20 38 38 74 ] M = A^T A = \begin{bmatrix} 20 & 38 \\ 38 & 74 \end{bmatrix} M=ATA=[2038?3874?]

特征方程為:

det ? ( M ? λ I ) = ∣ 20 ? λ 38 38 74 ? λ ∣ = ( 20 ? λ ) ( 74 ? λ ) ? 3 8 2 \det(M - \lambda I) = \begin{vmatrix} 20 - \lambda & 38 \\ 38 & 74 - \lambda \end{vmatrix} = (20 - \lambda)(74 - \lambda) - 38^2 det(M?λI)= ?20?λ38?3874?λ? ?=(20?λ)(74?λ)?382

展開:

( 20 ? λ ) ( 74 ? λ ) = 1480 ? 94 λ + λ 2 3 8 2 = 1444 (20 - \lambda)(74 - \lambda) = 1480 - 94\lambda + \lambda^2 \\ 38^2 = 1444 (20?λ)(74?λ)=1480?94λ+λ2382=1444

所以特征方程為:

λ 2 ? 94 λ + 36 = 0 \lambda^2 - 94\lambda + 36 = 0 λ2?94λ+36=0

使用求根公式:

λ = 94 ± 9 4 2 ? 4 ? 1 ? 36 2 = 94 ± 8836 ? 144 2 = 94 ± 8692 2 \lambda = \frac{94 \pm \sqrt{94^2 - 4 \cdot 1 \cdot 36}}{2} = \frac{94 \pm \sqrt{8836 - 144}}{2} = \frac{94 \pm \sqrt{8692}}{2} λ=294±942?4?1?36 ??=294±8836?144 ??=294±8692 ??

8692 ≈ 93.23 ? λ 1 ≈ 94 + 93.23 2 = 93.615 , λ 2 ≈ 94 ? 93.23 2 = 0.385 \sqrt{8692} \approx 93.23 \Rightarrow \lambda_1 \approx \frac{94 + 93.23}{2} = 93.615, \quad \lambda_2 \approx \frac{94 - 93.23}{2} = 0.385 8692 ?93.23?λ1?294+93.23?=93.615,λ2?294?93.23?=0.385


? 步驟三:取最大特征值的平方根

∥ A ∥ 2 = λ max ? = 93.615 ≈ 9.68 \|A\|_2 = \sqrt{\lambda_{\max}} = \sqrt{93.615} \approx 9.68 A2?=λmax? ?=93.615 ?9.68


? 對 A ? 1 A^{-1} A?1 做同樣的操作

先寫出:
A ? 1 = [ 5 6 ? 7 6 ? 1 3 2 3 ] ? ( A ? 1 ) T = [ 5 6 ? 1 3 ? 7 6 2 3 ] A^{-1} = \begin{bmatrix} \frac{5}{6} & -\frac{7}{6} \\ -\frac{1}{3} & \frac{2}{3} \end{bmatrix} \Rightarrow (A^{-1})^T = \begin{bmatrix} \frac{5}{6} & -\frac{1}{3} \\ -\frac{7}{6} & \frac{2}{3} \end{bmatrix} A?1=[65??31???67?32??]?(A?1)T=[65??67???31?32??]

計算 ( A ? 1 ) T A ? 1 (A^{-1})^T A^{-1} (A?1)TA?1

( A ? 1 ) T A ? 1 = [ 5 6 ? 1 3 ? 7 6 2 3 ] [ 5 6 ? 7 6 ? 1 3 2 3 ] = [ ( 5 6 ) 2 + ( ? 1 3 ) 2 5 6 ? ( ? 7 6 ) + ( ? 1 3 ) ? 2 3 ( ? 7 6 ) ? 5 6 + 2 3 ? ( ? 1 3 ) ( ? 7 6 ) 2 + ( 2 3 ) 2 ] (A^{-1})^T A^{-1} = \begin{bmatrix} \frac{5}{6} & -\frac{1}{3} \\ -\frac{7}{6} & \frac{2}{3} \end{bmatrix} \begin{bmatrix} \frac{5}{6} & -\frac{7}{6} \\ -\frac{1}{3} & \frac{2}{3} \end{bmatrix} = \begin{bmatrix} \left(\frac{5}{6}\right)^2 + \left(-\frac{1}{3}\right)^2 & \frac{5}{6} \cdot \left(-\frac{7}{6}\right) + \left(-\frac{1}{3}\right)\cdot\frac{2}{3} \\ \left(-\frac{7}{6}\right)\cdot\frac{5}{6} + \frac{2}{3}\cdot\left(-\frac{1}{3}\right) & \left(-\frac{7}{6}\right)^2 + \left(\frac{2}{3}\right)^2 \end{bmatrix} (A?1)TA?1=[65??67???31?32??][65??31???67?32??]=[(65?)2+(?31?)2(?67?)?65?+32??(?31?)?65??(?67?)+(?31?)?32?(?67?)2+(32?)2?]

逐項計算:

  • 第 (1,1) 項: 25 36 + 1 9 = 25 36 + 4 36 = 29 36 \frac{25}{36} + \frac{1}{9} = \frac{25}{36} + \frac{4}{36} = \frac{29}{36} 3625?+91?=3625?+364?=3629?
  • 第 (1,2) 項: ? 35 36 ? 2 9 = ? 35 36 ? 8 36 = ? 43 36 -\frac{35}{36} - \frac{2}{9} = -\frac{35}{36} - \frac{8}{36} = -\frac{43}{36} ?3635??92?=?3635??368?=?3643?
  • 第 (2,2) 項: 49 36 + 4 9 = 49 36 + 16 36 = 65 36 \frac{49}{36} + \frac{4}{9} = \frac{49}{36} + \frac{16}{36} = \frac{65}{36} 3649?+94?=3649?+3616?=3665?

所以:
( A ? 1 ) T A ? 1 = [ 29 36 ? 43 36 ? 43 36 65 36 ] (A^{-1})^T A^{-1} = \begin{bmatrix} \frac{29}{36} & -\frac{43}{36} \\ -\frac{43}{36} & \frac{65}{36} \end{bmatrix} (A?1)TA?1=[3629??3643???3643?3665??]

接下來求這個矩陣的最大特征值:

特征方程:
det ? ( [ 29 36 ? λ ? 43 36 ? 43 36 65 36 ? λ ] ) = 0 \det\left( \begin{bmatrix} \frac{29}{36} - \lambda & -\frac{43}{36} \\ -\frac{43}{36} & \frac{65}{36} - \lambda \end{bmatrix} \right) = 0 det([3629??λ?3643???3643?3665??λ?])=0

( 29 36 ? λ ) ( 65 36 ? λ ) ? ( 43 36 ) 2 = 0 \left(\frac{29}{36} - \lambda\right)\left(\frac{65}{36} - \lambda\right) - \left(\frac{43}{36}\right)^2 = 0 (3629??λ)(3665??λ)?(3643?)2=0

展開得:
λ 2 ? ( 94 36 ) λ + ( 29 ? 65 ? 4 3 2 3 6 2 ) = 0 \lambda^2 - \left(\frac{94}{36}\right)\lambda + \left(\frac{29\cdot65 - 43^2}{36^2}\right) = 0 λ2?(3694?)λ+(36229?65?432?)=0

計算常數項:

  • 29 ? 65 = 1885 29 \cdot 65 = 1885 29?65=1885
  • 4 3 2 = 1849 43^2 = 1849 432=1849
  • 所以分子為 1885 ? 1849 = 36 1885 - 1849 = 36 1885?1849=36

代入后:
λ 2 ? 94 36 λ + 36 1296 = λ 2 ? 94 36 λ + 1 36 = 0 \lambda^2 - \frac{94}{36}\lambda + \frac{36}{1296} = \lambda^2 - \frac{94}{36}\lambda + \frac{1}{36} = 0 λ2?3694?λ+129636?=λ2?3694?λ+361?=0

用求根公式解:
λ = 94 36 ± ( 94 36 ) 2 ? 4 36 2 ≈ 2.611 ± 6.817 ? 0.111 2 = 2.611 ± 6.706 2 ≈ 2.611 ± 2.590 2 \lambda = \frac{\frac{94}{36} \pm \sqrt{\left(\frac{94}{36}\right)^2 - \frac{4}{36}}}{2} \approx \frac{2.611 \pm \sqrt{6.817 - 0.111}}{2} = \frac{2.611 \pm \sqrt{6.706}}{2} \approx \frac{2.611 \pm 2.590}{2} λ=23694?±(3694?)2?364? ??22.611±6.817?0.111 ??=22.611±6.706 ??22.611±2.590?

得到兩個特征值:

  • λ 1 ≈ 5.201 2 = 2.6005 \lambda_1 \approx \frac{5.201}{2} = 2.6005 λ1?25.201?=2.6005
  • λ 2 ≈ 0.021 2 = 0.0105 \lambda_2 \approx \frac{0.021}{2} = 0.0105 λ2?20.021?=0.0105

所以:

∥ A ? 1 ∥ 2 = λ max ? = 2.6005 ≈ 1.613 \|A^{-1}\|_2 = \sqrt{\lambda_{\max}} = \sqrt{2.6005} \approx 1.613 A?12?=λmax? ?=2.6005 ?1.613


? 最終結果總結
矩陣譜范數
A A A ∣ A ∣ 2 ≈ 9.68 |A|_2 \approx 9.68 A2?9.68
A ? 1 A^{-1} A?1 ∣ A ? 1 ∣ 2 ≈ 1.613 |A^{-1}|_2 \approx 1.613 A?12?1.613

📌 驗證關系: ∥ A ? 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A?1A1?

1 ∥ A ∥ ≈ 1 9.68 ≈ 0.1033 < 1.613 = ∥ A ? 1 ∥ \frac{1}{\|A\|} \approx \frac{1}{9.68} \approx 0.1033 < 1.613 = \|A^{-1}\| A1?9.681?0.1033<1.613=A?1

? 成立!

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

如若內容造成侵權/違法違規/事實不符,請聯系多彩編程網進行投訴反饋email:809451989@qq.com,一經查實,立即刪除!

相關文章

成像光譜遙感技術中的AI革命:ChatGPT在遙感領域中的應用

課程將最新的人工智能技術與實際的遙感應用相結合&#xff0c;提供不僅是理論上的&#xff0c;而且是適用和可靠的工具和方法。無論你是經驗豐富的研究人員還是好奇的從業者&#xff0c;本課程都將為分析和解釋遙感數據開辟新的、有效的方法&#xff0c;使你的工作更具影響力和…

Debian12 安裝 sippts

試了試&#xff0c;貌似不復雜&#xff0c;記錄如下&#xff1a; apt-get install -y python3 python3 --version # 3.11.2 apt-get install -y python3-pip pip3 --version # 24.3.1 rm /usr/lib/python3.11/EXTERNALLY-MANAGED cd /usr/src git clone https://github.com/…

VR Panorama 360 PRO Renderer保姆級別教程

總覽: 全景圖及全景視頻錄制插件有兩個 一個是件(以下簡稱VR360插件) 一個是Unity官方的Unity Recorder插件(以下簡稱Recorder插件) 在圖片清晰度上VR 360插件要高于Recorder插件,所以渲染全景圖時,優先使用VR 360插件,當然全景視頻也可以使用VR360插件。 但VR 360插件…

cv610將音頻chn0配置為g711a,chn1配置為 aac編碼,記錄

cv610將音頻chn0配置為g711a,chn1配置為 aac編碼,記錄 工程代碼在文章底部 編譯時放在 sdk的同級目錄 sdk_version: sdk_V010,打了AOV的補丁 aenc可以配置為 chn0=g711a, chn1=aac 設置兩個編碼通道為不同編碼屬性 主要思路為在 ss_mpi_aenc_create_chn時將 chn1配置編碼為…

CAD2018,矩形設計,新增文字,塊新增與打散

一、矩形設計 1.選擇頁面&#xff0c;點擊左鍵&#xff0c;直接輸入【rec】&#xff0c;回車&#xff1b; 2.長按鼠標左鍵&#xff0c;拉出矩形&#xff0c;抬起左鍵。 3. 會生成一個矩形框。 4. ①輸入寬度數值&#xff0c;②輸入逗號切換到高度&#xff0c;③輸入高度。 5.成…

day047-部署我的世界-java版服務器

文章目錄 1. 官方地址2. Ubuntu配置服務端2.1 下載服務端jar包&#xff0c;并上傳2.2 安裝jdk2.3 啟動服務端2.4 設置云安全組 3. 客戶端-我的世界啟動器 1. 官方地址 官方服務端下載地址&#xff1a;[Minecraft 服務器下載 | Minecraft](https://www.minecraft.net/zh-hans/d…

飛算JavaAI:精準切中開發者痛點,專治“AI生成代碼不可用、邏輯混亂”的頑疾

飛算JavaAI&#xff1a;精準切中開發者痛點&#xff0c;專治“AI生成代碼不可用、邏輯混亂”的頑疾 一、前言二、關于飛算JavaAI2.1 飛算JavaAI來源2.2 飛算JavaAI超能力 三、飛算JavaAI我的另一半3.1 Idea安裝配置3.2 Main方法寫個九九乘法表3.3 Main方法寫個冒泡排序3.4 老項…

傘兵 釣魚的肝

題目描述 為了在敵國滲透作戰&#xff0c;指揮官決定&#xff1a;派出傘兵前往敵國&#xff01;然而敵國的風十分強烈&#xff0c;能讓傘兵在同一高度不停轉悠&#xff0c;直到被刮到一個無風區…… 輸入格式 第一行兩個整數 n,m&#xff0c;表示敵國的大小。 以下 n 行&am…

GoFastDFS:輕量級高性能分布式文件存儲解決方案(Linux安裝部署)

一、GoFastDFS概述 GoFastDFS是一個基于HTTP協議的分布式文件存儲系統&#xff0c;采用Go語言開發&#xff0c;具有輕量級、高性能、易部署等特點。它專為互聯網應用設計&#xff0c;特別適合圖片、視頻、文檔等中小文件的存儲與分發場景。 它基于大道至簡的設計理念&#xf…

python 原型污染 perl符號表污染 -- Google 2025 MYTHOS

題目實現了一個Game,分為前后端 part 1 前端存在明顯原型污染 def copy(src, dst):for k, v in src.items():if hasattr(dst, "__getitem__"):if dst.get(k) and type(v) dict:copy(v, dst.get(k))else:dst[k] velif hasattr(dst, k) and type(v) dict:copy(v, ge…

數據結構day4——棧

目錄 一、棧的核心概念 什么是棧&#xff1f; 棧的核心特性 二、棧的基本操作 三、C 語言實現棧的兩種方式 1. 順序棧&#xff08;基于數組實現&#xff09; 實現代碼 順序棧的優缺點 2. 鏈式棧&#xff08;基于鏈表實現&#xff09; 實現代碼 鏈式棧的優缺點 四、…

用戶系統的架構設計與實現策略(二)

一個用戶系統除了基本的用戶業務功能&#xff0c;還應囊括用戶的權限設計及其實現。這本文中我們將探討一下關于用戶權限的設計與實現方法論。 簡介 在構建現代應用系統的過程中&#xff0c;很少有設計決策會像訪問控制機制那樣&#xff0c;對安全性、可擴展性和用戶體驗產生…

深度學習-邏輯回歸

邏輯回歸的目的 邏輯回歸只判斷樣本屬于正類的概率是多大&#xff0c;0-1之間 找到一組最佳的權重&#xff08;w1,w2,w3,…&#xff09; &#xff0c;b&#xff0c;使得模型預測的概率 P(Y1) 盡可能接近樣本的真實標簽&#xff08;1 或 0&#xff09;。 計算過程 前向傳播過程…

對象池模式:減少GC的Kotlin實戰指南

對象池模式通過對象復用機制&#xff0c;將對象生命周期從"創建-銷毀"轉變為"借出-歸還"&#xff0c;顯著減少GC壓力。下面通過完整實例展示其實現細節。 一、對象池工作原理圖解 #mermaid-svg-Edrz4np9hD6DJdNi {font-family:"trebuchet ms",v…

Java接口報錯:Packet for query is too large - 解決方案與架構思考

Java接口報錯&#xff1a;Packet for query is too large - 解決方案與架構思考 背景與技術原理解決方案體系&#xff08;擴展版&#xff09;一、MySQL服務端配置&#xff08;永久生效&#xff09;配置文件修改&#xff08;推薦生產環境&#xff09; 文件路徑參考Linux: /etc/m…

7月2日作業

思維導圖 一、創建一個進程扇 代碼 #include <25041head.h>int main(int argc, const char *argv[]) {pid_t pid;for(int i1;i<4;i){pidfork();if(pid>0){sleep(1);}if(pid0){printf("我是子進程%d:%d,父進程%d\n",i,getpid(),getppid());sleep(1);re…

設計模式(九)

職責鏈模式&#xff08;Chain of Responsibility&#xff09;詳解 一、核心概念 職責鏈模式將請求的發送者和接收者解耦&#xff0c;使多個對象都有機會處理請求。這些對象連接成一條鏈&#xff0c;請求沿著鏈傳遞&#xff0c;直到有一個對象處理它為止。該模式允許動態調整處…

左神算法之Zigzag方式打印矩陣

目錄 Zigzag方式打印矩陣1. 題目2. 解釋3. 思路4. 代碼5. 總結 Zigzag方式打印矩陣 1. 題目 用zigzag的方式打印矩陣&#xff0c;比如下面的矩陣&#xff1a; 0 1 2 3 4 5 6 7 8 9 10 11打印順序為&#xff1a;0 1 4 8 5 2 3 6 9 10 7 11 2. 解釋 Zigzag打印矩陣是指按照…

【前端批量下載圖片,并打包成壓縮包下載】

一、需求說明 我現在有個需求&#xff1a; 1.列表中有個下載按鈕&#xff0c;點擊下載&#xff0c;將列表中所有的圖片打成壓縮包&#xff0c;并下載 2.效果演示點擊查看效果 最終效果&#xff1a; 二、安裝下載插件 實現此功能需要兩個插件&#xff1a;jszip、file-saver …

NV133NV137美光固態閃存NV147NV148

NV133NV137美光固態閃存NV147NV148 美光固態閃存技術矩陣深度解析&#xff1a;NV133至NV148的全面較量 一、性能參數&#xff1a;數據高速公路的“車速”比拼 讀寫速度&#xff1a;從“鄉間小道”到“高鐵動脈” 美光NV系列固態閃存的核心競爭力在于其讀寫速度的躍升。以NV15…