向量与矩阵范数_矩阵范数与谱半径的关系

向量与矩阵范数_矩阵范数与谱半径的关系范数(norm),是具有“长度”概念的函数。在线性代数、泛函分析及相关的数学领域,范函是一个函数,其为矢量空间内的所有矢量赋予非零的正长度或大小。半范数反而可以为非零的矢量赋予零长度。举一个简单的例子,在二维的欧氏几何空间R就可定义欧氏范数。在这个矢量空间中的元素常常在笛卡儿坐标系统中被画成一个从原点出发的带有箭头的有向线段。每一个矢量的欧氏范数就是有向线段的长度。其中定义范数的矢量空间就是赋范矢


Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

范数(norm),是具有“长度”概念的函数。在线性代数、泛函分析及相关的数学领域,范函是一个函数,其为矢量空间内的所有矢量赋予非零的正长度或大小。半范数反而可以为非零的矢量赋予零长度。

举一个简单的例子,在二维的欧氏几何空间 R就可定义欧氏范数。在这个矢量空间中的元素常常在笛卡儿坐标系统中被画成一个从原点出发的带有箭头的有向线段。每一个矢量的欧氏范数就是有向线段的长度。

其中定义范数的矢量空间就是赋范矢量空间。同样,其中定义半范数的矢量空间就是赋半范矢量空间。


向量范数

概念

设函数 f:Rn→R ,若 f 满足

  1. 正定性:

    f(x)≥0,∀x∈Rn
    ,等号当且仅当 x=0 时候成了。

    • 正齐次性: f(αx)=|α|×f(x),∀x∈Rn,∀α∈R 。
    • 次可加性(三角不等式): f(x+y)≤f(x)+f(y) 。
    • 例如:

      ∥x∥p=(∑i=1n|xi|p)1p,p∈[1,∞)是Rn上的向量范数。

      常见的范数

      这里以 Cn 空间为例, Rn 空间类似。
      最常用的范数就是p-范数。若 x=[x1,x2,…,xn]T ,那么

      ∥x∥p=(|x1|p+|x2|p+⋯+|xn|p)1p

      可以验证p-范数确实满足范数的定义。其中三角不等式的证明不是平凡的,这个结论通常称为闵可夫斯基(Minkowski)不等式。

      当p取1,2,∞的时候分别是以下几种最简单的情形:
      1-范数: ║x║1=│x1│+│x2│+…+│xn│
      2-范数: ║x║2=(│x1│2+│x2│2+…+│xn│2)1/2
      ∞-范数(最大范数): ║x║∞=max(│x1│,│x2│,…,│xn│)

      其中2-范数就是通常意义下的距离。

      范数的性质

      (1)连续性
      设 f 是

      Rn
      上的任意一个范数,则 f 关于

      x
      的每个分量是连续的。

      (2)等价性
      设 ∥⋅∥s 和 ∥⋅∥t 是 Rn 上的任意两个范数,则存在常数 c1 和 c2 ,使得对任意的 x∈Rn 有

      c1∥x∥s≤∥x∥t≤c2∥x∥s

      (3)不等式

      对于这些常用的范数有以下不等式: ║x║∞≤║x║2≤║x║1≤n12║x║2≤n║x║∞
      另外,若p和q是赫德尔(Hölder)共轭指标,即1/p+1/q=1,那么有赫德尔不等式:

      |<x,y>|=||xH∗y|≤║x║p║y║q

      当p=q=2时就是柯西-许瓦兹(Cauchy-Schwarz)不等式。

      (4)向量序列的收敛性

      limk→∞x(k)=x∗⇐⇒limk→∞∥x(k)−x∗∥=0


      矩阵范数

      概念

      一般来讲矩阵范数除了正定性,齐次性和三角不等式之外,还规定其必须满足相容性: ║XY║≤║X║║Y║ 。

      所以矩阵范数通常也称为相容范数。

      如果 ║⋅║α 是相容范数,且任何满足 ║⋅║β≤║⋅║α 的范数 ║⋅║β 都不是相容范数,那么 ║⋅║α 称为极小范数。对于 n 阶实方阵(或复方阵)全体上的任何一个范数

      ║⋅║
      ,总存在唯一的实数 k>0 ,使得 k×║⋅║ 是极小范数。
      注:如果不考虑相容性,那么矩阵范数和向量范数就没有区别,因为 m×n 矩阵全体和 mn 维向量空间同构。引入相容性主要是为了保持矩阵作为线性算子的特征,这一点和算子范数的相容性一致,并且可以得到Mincowski定理以外的信息。

      常见的矩阵范数

      满足以上设定的矩阵范数可以有多种。由于它们都是定义在 Mm,n(K) 这个有限维向量空间上的范数,所以实质上是等价的。常见的矩阵范数通常是在矩阵的应用中自然定义或诱导的范数。

      向量范数诱导的矩阵范数

      考虑从向量空间 V=Km 映射到 W=Kn 的所有线性映射的构成的空间: Lm,n(K) 。设 V 和

      W
      中分别装备了两个向量范数 ∥⋅∥V 和 ∥⋅∥W ,则可以定义 Lm,n(K) 上的算子范数 ∥⋅∥L:

      ∀A∈Lm,n(K)

      ∥A∥L=max{
      ∥A(x)∥W;x∈V,∥x∥V⩽1}

      而给定了基底后,每个从V映射到W的线性映射都可以用一个m \times n的矩阵来表示,所以同样地可以定义 Mm,n(K) 上的非负映射 ∥⋅∥M:

      ∀A∈Mm,n(K)

      ∥A∥L=max{
      ∥A(x)∥W;x∈V,∥x∥V⩽1}

      可以验证, ∥⋅∥M 满足矩阵范数的定义,因此是一个矩阵范数。这个矩阵范数被称为是由向量空间范数诱导的矩阵范数,可以看作是算子范数在由有限维向量空间之间线性映射组成的空间上的特例。如果 m=n ,所对应的矩阵空间就是n阶方块矩阵空间 Mn(K) 。这时可以验证,诱导范数 ∥⋅∥M 满足一致性条件。

      p-范数诱导的矩阵范数

      当 V 和

      W
      中装备的向量范数都是p-范数的时候,诱导的矩阵范数也称为矩阵的诱导p-范数。具体来说就是:

      ∥A∥p=maxx≠0∥Ax∥p∥x∥p=maxx≠0(∑ni=1|∑mj=1Aijxj|p)1/p(∑mi=1|xi|p)1/p

      在p=1和 p=∞ 的情况下,其范数可以以下方式计算:

      ∥A∥1=max1≤j≤n∑i=1m|aij|∥A∥∞=max1≤i≤m∑j=1n|aij|.


      这些与矩阵的Schatten p-范数不同,也可以用 ∥A∥p 来表示。
      当p = 2(欧几里德范数)且m = n(方阵)时,诱导的矩阵范数就是谱范数。矩阵A的谱范数是A最大的奇异值或半正定矩阵A*A的最大特征值的平方根:

      ∥A∥2=λmax(A∗A)−−−−−−−−−√


      其中 A∗ 代表A的共轭转置。

      矩阵元范数

      这些向量范数将矩阵视为 m×n 向量,并使用类似的向量范数。
      举例说明,使用向量的p-范数,我们得到:

      ∥A∥p=(∑i=1m∑j=1n|aij|p)1/p


      注:不要把矩阵元p-范数与诱导p-范数混淆。

      弗罗贝尼乌斯范数

      对 p=2 ,这称为弗罗贝尼乌斯范数(Frobenius norm)或希尔伯特-施密特范数(Hilbert–Schmidt norm),不过后面这个术语通常只用于希尔伯特空间。这个范数可用不同的方式定义:

      ∥A∥F=∑i=1m∑j=1n|aij|2−−−−−−−−−−⎷=trace(A∗A)−−−−−−−−−√=∑i=1min{
      m,n}
      σ2i
      −−−−−−−−−⎷


      这里 A∗ 表示 A 的共轭转置,

      σi
      是 A 的奇异值,并使用了迹函数。弗罗贝尼乌斯范数与

      Kn
      上欧几里得范数非常类似,来自所有矩阵的空间上一个内积。
      弗罗贝尼乌斯范范数是服从乘法的且在数值线性代数中非常有用。这个范数通常比诱导范数容易计算。

      矩阵范数性质

      矩阵范数的性质

      (1)连续性:设 f 是

      Rn×n
      上的任一矩阵范数,则 f 关于

      A
      的每个分量是连续的。

      (2)等价性:设 ∥⋅∥s 和 ∥⋅∥t 是 Rn×n 上的任意两个矩阵范数,则存在常数 c1 和 c2 ,使得对任意的 A∈Rn×n 有

      c1∥A∥s≤∥A∥t≤c2∥A∥s


      (3)
      任何诱导的矩阵范数都满足此不等式

      ∥A∥≥ρ(A)


      其中 ρ(A) 是 A 的谱半径。事实上,可以证明

      ρ(A)
      是 A 的所有诱导范数的下界。
      此外,我们有


      limr→∞∥Ar∥1/r=ρ(A)
      。
      若A是对称矩阵,则 ρ(A)=∥A∥2

      算子范数的性质

      定理:
      设 ∥⋅∥ 是 Rn 上的任一向量范数,其对应的算子范数也记为 ∥⋅∥ ,则有

      ∥Ax∥≤∥A∥×∥x∥

      定理:
      设 ∥⋅∥ 是任一算子范数,则 ρ(A)≤∥A∥

      定理:
      对任意 ε>0 ,总存在一个算子范数 ∥⋅∥ε ,使得

      ∥⋅∥ε≤ρ(A)+ε

      这个图的计算可能有一些问题
      这里写图片描述

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/191846.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长

相关推荐

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号