矩阵的范数和逆矩阵的范数的关系

文章目录

矩阵的范数和逆矩阵的范数之间的关系是线性代数中的一个重要主题,尤其是在数值分析和矩阵理论中。

前提条件

  • 我们讨论的是方阵 A A A 的情况,因为只有方阵才可能存在逆矩阵。
  • 矩阵 A A A可逆的(即 det ⁡ ( A ) ≠ 0 \det(A) \neq 0 det(A)=0),这意味着 A − 1 A^{-1} A1 存在。

关键结论

对于一个给定的矩阵范数 ∥ ⋅ ∥ \|\cdot\| ,如果 A A A 是可逆的,则有以下不等式成立:

∥ A − 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A1A1

其中, ∥ A ∥ \|A\| A 表示矩阵 A A A 的范数,而 ∥ A − 1 ∥ \|A^{-1}\| A1 则表示其逆矩阵的范数。注意这里的范数可以是任何相容的矩阵范数,如诱导范数(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\|} A1Axx

对所有非零向量 x x x 成立。特别地,取 x = A − 1 z x = A^{-1}z x=A1z 对于任意非零向量 z z z,我们得到:

1 ∥ A ∥ ≤ ∥ A − 1 z ∥ ∥ z ∥ \frac{1}{\|A\|} \leq \frac{\|A^{-1}z\|}{\|z\|} A1zA1z

从而得出:

∥ A − 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A1A1

注意事项

  • 这个关系表明了逆矩阵的范数至少为原矩阵范数的倒数。然而,在实际应用中, ∥ A − 1 ∥ \|A^{-1}\| A1 可能远大于 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)=AA1,它衡量了矩阵求解线性系统时的敏感度。条件数越大,意味着矩阵越“病态”,数值计算的结果可能越不稳定。

示例说明

假设有一个 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=[4275]

其逆矩阵为:

A − 1 = [ 5 / 6 − 7 / 6 − 1 / 3 2 / 3 ] A^{-1} = \begin{bmatrix} 5/6 & -7/6 \\ -1/3 & 2/3 \end{bmatrix} A1=[5/61/37/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=[4725],ATA=[4725][4275]=[42+2274+5247+2572+52]=[20383874]


✅ 步骤二:求 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=[20383874]

特征方程为:

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λ383874λ =(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λ)=148094λ+λ2382=1444

所以特征方程为:

λ 2 − 94 λ + 36 = 0 \lambda^2 - 94\lambda + 36 = 0 λ294λ+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±9424136 =294±8836144 =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λ1294+93.23=93.615,λ229493.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} A1 做同样的操作

先写出:
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} A1=[65316732](A1)T=[65673132]

计算 ( A − 1 ) T A − 1 (A^{-1})^T A^{-1} (A1)TA1

( 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} (A1)TA1=[65673132][65316732]=[(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} 363592=3635368=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} (A1)TA1=[3629364336433665]

接下来求这个矩阵的最大特征值:

特征方程:
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λ364336433665λ])=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)λ+(3622965432)=0

计算常数项:

  • 29 ⋅ 65 = 1885 29 \cdot 65 = 1885 2965=1885
  • 4 3 2 = 1849 43^2 = 1849 432=1849
  • 所以分子为 1885 − 1849 = 36 1885 - 1849 = 36 18851849=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 λ23694λ+129636=λ23694λ+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)2364 22.611±6.8170.111 =22.611±6.706 22.611±2.590

得到两个特征值:

  • λ 1 ≈ 5.201 2 = 2.6005 \lambda_1 \approx \frac{5.201}{2} = 2.6005 λ125.201=2.6005
  • λ 2 ≈ 0.021 2 = 0.0105 \lambda_2 \approx \frac{0.021}{2} = 0.0105 λ220.021=0.0105

所以:

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


✅ 最终结果总结
矩阵谱范数
A A A ∣ A ∣ 2 ≈ 9.68 |A|_2 \approx 9.68 A29.68
A − 1 A^{-1} A1 ∣ A − 1 ∣ 2 ≈ 1.613 |A^{-1}|_2 \approx 1.613 A121.613

📌 验证关系: ∥ A − 1 ∥ ≥ 1 ∥ A ∥ \|A^{-1}\| \geq \frac{1}{\|A\|} A1A1

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}\| A19.6810.1033<1.613=A1

✅ 成立!

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