【数值计算方法(黄明游)】解线性代数方程组的迭代法(一):向量、矩阵范数与谱半径【理论到程序】

文章目录

  • 一、向量、矩阵范数与谱半径
    • 1、向量范数
    • 2、矩阵范数
      • a. 矩阵的范数
      • b. 常见的矩阵范数
        • 相容范数
        • 算子范数
    • 3、谱半径
    • 4、知识点总结
      • 1. 向量范数
      • 2. 矩阵范数
      • 3. 谱半径
    • 5、计算例题

  注意:速读可直接跳转至“4、知识点总结”及“5、计算例题”部分

一、向量、矩阵范数与谱半径

  当涉及到线性代数和矩阵理论时,向量、矩阵范数以及谱半径是非常重要的概念,下面将详细介绍这些内容:

1、向量范数

a. 定义及性质

  考虑一个 n n n 维向量 x x x,定义一个实值函数 N ( x ) N(x) N(x),记作 N ( x ) = ∥ x ∥ N(x) = \|x\| N(x)=x。如果 N ( x ) N(x) N(x) 满足以下条件,那么它就是 x x x 上的一个向量范数(或向量模):

  1. 非负性: N ( x ) ≥ 0 N(x) \geq 0 N(x)0,且 N ( x ) = 0 N(x) = 0 N(x)=0当且仅当 x x x 是零向量。

∥ x ∥ ≥ 0 \|x\| \geq 0 x0 ∥ x ∥ = 0  当且仅当  x = 0 \|x\| = 0 \text{ 当且仅当 } x = \mathbf{0} x=0 当且仅当 x=0

  1. 齐次性: 对于任意实数 α \alpha α(或复数),有 N ( α x ) = ∣ α ∣ ⋅ N ( x ) N(\alpha x) = |\alpha| \cdot N(x) N(αx)=αN(x)

∥ α x ∥ = ∣ α ∣ ⋅ ∥ x ∥ \| \alpha x \| = |\alpha| \cdot \|x\| αx=αx

  1. 三角不等式: 对于任意向量 x x x y y y,有 N ( x + y ) ≤ N ( x ) + N ( y ) N(x + y) \leq N(x) + N(y) N(x+y)N(x)+N(y)

    ∥ x + y ∥ ≤ ∥ x ∥ + ∥ y ∥ \|x + y\| \leq \|x\| + \|y\| x+yx+y

补充解释
  • 非负性: 范数是非负的,即它不会为负值。当且仅当向量是零向量时,范数为零。

  • 齐次性: 范数在缩放(乘以常数)下保持一致,即放大或缩小向量会按比例影响其范数。

  • 三角不等式: 范数的三角不等式表示通过两边之和的方式度量两个向量之间的距离。它确保了向量空间中的“三角形”不会变得扭曲。

范数差

  由上述三角不等式可推导出: ∥ x − y ∥ ≥ ∣ ∥ x ∥ − ∥ y ∥ ∣ \|x - y\| \geq |\|x\| - \|y\|| xy∣∥xy∥∣

  • 推导过程
    • 根据向量范数的三角不等式,对于任意向量 x x x y y y,有: ∥ x − y ∥ ≤ ∥ x ∥ + ∥ y ∥ \|x - y\| \leq \|x\| + \|y\| xyx+y 其中

b. 常见的向量范数

l 1 l_1 l1 l 2 l_2 l2 l ∞ l_\infty l 范数

  对于一个 n n n维向量 x = ( x 1 , x 2 , … , x n ) x = (x_1, x_2, \ldots, x_n) x=(x1,x2,,xn)

  1. l 1 l_1 l1 范数:
    ∥ x ∥ 1 = ∑ i = 1 n ∣ x i ∣ \|x\|_1 = \sum_{i=1}^{n} |x_i| x1=i=1nxi

  2. l 2 l_2 l2 范数:
    ∥ x ∥ 2 = ∑ i = 1 n x i 2 \|x\|_2 = \sqrt{\sum_{i=1}^{n} x_i^2} x2=i=1nxi2

  3. l ∞ l_\infty l 范数:
    ∥ x ∥ ∞ = max ⁡ 1 ≤ i ≤ n ∣ x i ∣ \|x\|_\infty = \max_{1 \leq i \leq n} |x_i| x=1inmaxxi

性质
  • 非负性:

    ∥ x ∥ 1 , ∥ x ∥ 2 , ∥ x ∥ ∞ ≥ 0 \|x\|_1, \|x\|_2, \|x\|_\infty \geq 0 x1,x2,x0

  • 齐次性: 对于每个 x x x 和标量 α \alpha α,这三种范数都满足齐次性,即
    ∥ α x ∥ 1 = ∣ α ∣ ⋅ ∥ x ∥ 1 \|\alpha x\|_1 = |\alpha| \cdot \|x\|_1 αx1=αx1 ∥ α x ∥ 2 = ∣ α ∣ ⋅ ∥ x ∥ 2 \|\alpha x\|_2 = |\alpha| \cdot \|x\|_2 αx2=αx2 ∥ α x ∥ ∞ = ∣ α ∣ ⋅ ∥ x ∥ ∞ \|\alpha x\|_\infty = |\alpha| \cdot \|x\|_\infty αx=αx

  • 三角不等式: 对于每对向量 x x x y y y,这三种范数都满足三角不等式:
    ∥ x + y ∥ 1 ≤ ∥ x ∥ 1 + ∥ y ∥ 1 \|x + y\|_1 \leq \|x\|_1 + \|y\|_1 x+y1x1+y1 ∥ x + y ∥ 2 ≤ ∥ x ∥ 2 + ∥ y ∥ 2 \|x + y\|_2 \leq \|x\|_2 + \|y\|_2 x+y2x2+y2 ∥ x + y ∥ ∞ ≤ ∥ x ∥ ∞ + ∥ y ∥ ∞ \|x + y\|_\infty \leq \|x\|_\infty + \|y\|_\infty x+yx+y

关系
  • l 1 l_1 l1 范数、 l 2 l_2 l2 范数、 l ∞ l_\infty l 范数之间存在关系:
    ∥ x ∥ ∞ ≤ ∥ x ∥ 2 ≤ n ∥ x ∥ ∞ \|x\|_\infty \leq \|x\|_2 \leq \sqrt{n}\|x\|_\infty xx2n x ∥ x ∥ ∞ ≤ ∥ x ∥ 1 ≤ n ∥ x ∥ ∞ \|x\|_\infty \leq \|x\|_1 \leq n\|x\|_\infty xx1nx

2、矩阵范数

a. 矩阵的范数

  矩阵的范数是定义在矩阵空间上的实值函数,用于度量矩阵的大小或度量。对于一个矩阵 A A A,矩阵范数通常表示为 N ( A ) N(A) N(A) ∣ ∣ A ∣ ∣ ||A|| ∣∣A∣∣,满足以下条件:

  1. 非负性(Non-negativity):对于任意矩阵 A A A,有 N ( A ) ≥ 0 N(A) \geq 0 N(A)0,且等号成立当且仅当 A A A 是零矩阵。

  2. 齐次性(Homogeneity):对于任意标量 k k k 和矩阵 A A A,有 N ( k A ) = ∣ k ∣ ⋅ N ( A ) N(kA) = |k| \cdot N(A) N(kA)=kN(A)

  3. 三角不等式(Triangle Inequality):对于任意两个矩阵 A A A B B B,有 N ( A + B ) ≤ N ( A ) + N ( B ) N(A + B) \leq N(A) + N(B) N(A+B)N(A)+N(B)

b. 常见的矩阵范数

相容范数
  • 对于任意两个矩阵 A A A B B B,有 ∣ ∣ A B ∣ ∣ ≤ ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ B ∣ ∣ ||AB|| \leq ||A|| \cdot ||B|| ∣∣AB∣∣∣∣A∣∣∣∣B∣∣,这被称为相容性质。
  • 对于任意矩阵 A A A 和向量 x x x,有 ∣ ∣ A x ∣ ∣ ≤ ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ x ∣ ∣ ||Ax|| \leq ||A|| \cdot ||x|| ∣∣Ax∣∣∣∣A∣∣∣∣x∣∣,这也是相容性质。
算子范数

在这里插入图片描述

具体而言,常用的算子范数是 p p p范数,其中 p p p 是一个实数。

  • p = ∞ p = \infty p= 时,算子范数被定义为矩阵行的绝对值之和的最大值。即,
    ∣ ∣ A ∣ ∣ ∞ = max ⁡ 1 ≤ i ≤ n ∑ j = 1 n ∣ a i j ∣ ||A||_\infty = \max_{1 \leq i \leq n} \sum_{j=1}^n |a_{ij}| ∣∣A=1inmaxj=1naij
  • p = 1 p = 1 p=1 时,算子范数被定义为矩阵列的绝对值之和的最大值。即,
    ∣ ∣ A ∣ ∣ 1 = max ⁡ 1 ≤ j ≤ n ∑ i = 1 n ∣ a i j ∣ ||A||_1 = \max_{1 \leq j \leq n} \sum_{i=1}^n |a_{ij}| ∣∣A1=1jnmaxi=1naij
  • p = 2 p = 2 p=2 时,算子范数被定义为 A A A 的谱半径。谱半径是矩阵的特征值的按模最大值,表示为 p ( A ) = max ⁡ ∣ λ ∣ p(A) = \max |\lambda| p(A)=maxλ其中 λ \lambda λ A A A 的特征值。

3、谱半径

  待完善……

4、知识点总结

1. 向量范数

  • l 1 l_1 l1 范数(曼哈顿范数)
    ∣ ∣ x ∣ ∣ 1 = ∑ i = 1 n ∣ x i ∣ ||x||_1 = \sum_{i=1}^{n} |x_i| ∣∣x1=i=1nxi

  • l 2 l_2 l2 范数(欧几里得范数)
    ∣ ∣ x ∣ ∣ 2 = ∑ i = 1 n x i 2 ||x||_2 = \sqrt{\sum_{i=1}^{n} x_i^2} ∣∣x2=i=1nxi2

  • l ∞ l_\infty l 范数(无穷范数)
    ∣ ∣ x ∣ ∣ ∞ = max ⁡ 1 ≤ i ≤ n ∣ x i ∣ ||x||_\infty = \max_{1 \leq i \leq n} |x_i| ∣∣x=1inmaxxi

2. 矩阵范数

  • 弗罗贝尼乌斯范数(矩阵中每项数的平方和的开方值)
    ∣ ∣ A ∣ ∣ F = ∑ i = 1 n ∑ j = 1 n ∣ a i j ∣ 2 ||A||_F = \sqrt{\sum_{i=1}^{n} \sum_{j=1}^{n} |a_{ij}|^2} ∣∣AF=i=1nj=1naij2
  • 算子范数
    • 行和范数:当 p = ∞ p = \infty p= 时,算子范数被定义为矩阵中各行元素按绝对值求和所得的最大和数,即,
      ∣ ∣ A ∣ ∣ ∞ = max ⁡ 1 ≤ i ≤ n ∑ j = 1 n ∣ a i j ∣ ||A||_\infty = \max_{1 \leq i \leq n} \sum_{j=1}^n |a_{ij}| ∣∣A=1inmaxj=1naij
    • 列和范数:当 p = 1 p = 1 p=1 时,算子范数被定义为
      矩阵列的绝对值之和的最大值。即,
      ∣ ∣ A ∣ ∣ 1 = max ⁡ 1 ≤ j ≤ n ∑ i = 1 n ∣ a i j ∣ ||A||_1 = \max_{1 \leq j \leq n} \sum_{i=1}^n |a_{ij}| ∣∣A1=1jnmaxi=1naij
    • p = 2 p = 2 p=2 时,算子范数即 A A A 的谱半径,谱半径是矩阵的特征值的按模最大值
      ∣ ∣ A ∣ ∣ 2 = λ max ( A T A ) = p ( A ) = max ⁡ ∣ λ ∣ ||A||_2 = \sqrt{\lambda_{\text{max}}(A^TA)} = p(A) = \max |\lambda| ∣∣A2=λmax(ATA) =p(A)=maxλ

3. 谱半径

  谱半径是矩阵的特征值按模最大的那个值,对于一个 n × n n \times n n×n 的矩阵 A A A,其谱半径 p ( A ) p(A) p(A) 定义为:

p ( A ) = max ⁡ { ∣ λ ∣   ∣   λ  是  A  的特征值 } p(A) = \max \{|\lambda| \ | \ \lambda \text{ 是 } A \text{ 的特征值}\} p(A)=max{λ  λ  A 的特征值}

5、计算例题

对于矩阵 A = [ 2 1 − 1 4 ] A = \begin{bmatrix} 2 & 1 \\ -1 & 4 \end{bmatrix} A=[2114]计算其各种范数:

∥ A ∥ 1 = max ⁡ j ∑ i ∣ a i j ∣ = max ⁡ { 3 , 5 } = 5 \|A\|_1 = \max_j \sum_i |a_{ij}| = \max\{3, 5\} = 5 A1=jmaxiaij=max{3,5}=5

∥ A ∥ ∞ = max ⁡ i ∑ j ∣ a i j ∣ = max ⁡ { 3 , 5 } = 5 \|A\|_\infty = \max_i \sum_j |a_{ij}| = \max\{3, 5\} = 5 A=imaxjaij=max{3,5}=5

∥ A ∥ 2 = λ max ( A T A ) \|A\|_2 = \sqrt{\lambda_{\text{max}}(A^TA)} A2=λmax(ATA)

计算 A T A A^TA ATA 的特征值,找到最大特征值 λ max \lambda_{\text{max}} λmax

A T A = [ 5 − 2 − 2 17 ] A^TA = \begin{bmatrix} 5 & -2 \\ -2 & 17 \end{bmatrix} ATA=[52217]

特征值为 $\lambda = $。

∥ A ∥ 2 = λ max = \|A\|_2 = \sqrt{\lambda_{\text{max}}} = \sqrt{} A2=λmax =

  1. 谱半径:

    p ( A ) = max ⁡ { ∣ λ ∣ } p(A) = \max \{|\lambda|\} p(A)=max{λ}

    A A A 求特征值,找到最大的绝对值。

  • 1范数:5
  • ∞范数:5
  • 2范数:
  • 谱半径:

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

Dexie 查询sql速度优化

Dexie查询速度慢的原因主要一个优化点是复杂查询下的count执行。 以下摘自Dexie官方文档:https://dexie.org/docs/Collection/Collection.count() If executed on simple queries, the native IndexedDB ObjectStore count() method will be called (fast execution…

PPT插件-好用的插件-字距快速设置-大珩助手

字距快速设置 包含两端对齐、段首缩进、取消缩进、字间距、行间距、段后距 段首缩进 每次缩进两个字符,可对选中的文字、选中的多个文本对象两个层级操作 取消缩进 将缩进取消,可对选中的文字、选中的多个文本对象两个层级操作 字间距 预设了常用…

网页设计--第6次课后作业

试用Vue相关指令完成对以下json数据的显示。显示效果如下: 其中:gender1 显示为女,gender2显示为男。价格超过30元,显示“有点小贵”。价格少于等于30元,则显示“价格亲民”。 data: {books: [{"id": "…

【软件安装】VMware安装Centos7虚拟机并且设置静态IP,实现Windows和Centos7网络互相访问

这篇文章,主要介绍VMware安装Centos7虚拟机并且设置静态IP,实现Windows和Centos7网络互相访问。 目录 一、VMware安装Centos7 1.1、下载Centos7镜像 1.2、安装Centos7系统 二、设置静态IP地址 2.1、查看虚拟机网络IP 2.2、禁用NetworkManager服务 …

PandoraFMS 监控软件 SQL注入漏洞复现

0x01 产品简介 Pandora FMS是西班牙Artica公司的一套监控系统。该系统通过可视化的方式监控网络、服务器、虚拟基础架构和应用程序等。 0x02 漏洞概述 Pandora FMS监控软件存在SQL注入漏洞,攻击者通过chart_generator.php 来执行恶意语句,获取数据库敏感信息。 0x03 复现…

机器学习之全面了解回归学习器

我们将和大家一起探讨机器学习与数据科学的主题。 本文主要讨论大家针对回归学习器提出的问题。我将概要介绍,然后探讨以下五个问题: 1. 能否将回归学习器用于时序数据? 2. 该如何缩短训练时间? 3. 该如何解释不同模型的结果和…

拦截 open调用 (进程白名单,文件白名单)

拦截 open 文章目录 拦截 open第一个需求文件结构进程白名单文件白名单 测试代码第一个版本版本二代码演示 增加一个日志记录代码解释 gcc -shared -fPIC -o libintercept.so intercept.c -ldlLD_PRELOAD./libintercept.so ./processA在Linux中,我们可以使用LD_PREL…

RT-DETR优化:Backbone改进 | UniRepLKNet,通用感知大内核卷积网络,RepLK改进版本 | 2023.11

🚀🚀🚀本文改进: UniRepLKNet,通用感知大内核卷积网络,ImageNet-22K预训练,精度和速度SOTA,ImageNet达到88%, COCO达到56.4 box AP,ADE20K达到55.6 mIoU 🚀🚀🚀RT-DETR改进创新专栏:http://t.csdnimg.cn/vuQTz 学姐带你学习YOLOv8,从入门到创新,轻轻松松…

无人机高空巡查+智能视频监控技术,打造森林防火智慧方案

随着冬季的到来,森林防火的警钟再次敲响,由于森林面积广袤,地形复杂,且人员稀少,一旦发生火灾,人员无法及时发现,稍有疏忽就会酿成不可挽救的大祸。无人机高空巡查智能视频监控是一种非常有效的…

库卡LBR_iisy_3_R760协作机器人导入到coppeliasim

1.从库卡官网xpert下载模型 一般载都是这个step文件格式,其他的好像不太好用。coppeliasim导入格式用的是stl,需要用freeCAD打开重新转换一下。下载下来后,很多都是一个整体,在freeCAD导入中,导入选择要不勾选合并。 下载完用CAD …

Day54力扣打卡

打卡记录 出租车的最大盈利(动态规划) 链接 class Solution:def maxTaxiEarnings(self, n: int, rides: List[List[int]]) -> int:d defaultdict(list)for start, end, w in rides:d[end].append((start, end - start w))f [0] * (n 1)for i in…

linux进入emergency mode

问题描述 linux系统进入emergency mode模式 解决方法 查看问题原因 journalctl -xb -p3 使用fsck 不一定是sda2,也可能是其他,我的是/dev/sda6,然后接受所有的option,完毕后重启电脑 fsck /dev/sda2接受所有的选项&#xff…

华为配置Smart Link主备备份示例

定义 Smart Link,又叫做备份链路。一个Smart Link由两个接口组成,其中一个接口作为另一个的备份。Smart Link常用于双上行组网,提供可靠高效的备份和快速的切换机制。 Monitor Link是一种接口联动方案,它通过监控设备的上行接口…

Leetcode—2646.最小化旅行的价格总和【困难】

2023每日刷题&#xff08;五十三&#xff09; Leetcode—2646.最小化旅行的价格总和 算法思想 看灵神的 实现代码 class Solution { public:int minimumTotalPrice(int n, vector<vector<int>>& edges, vector<int>& price, vector<vector&l…

Spring Boot 整合 xxl-job 保姆级教程!

文章目录 介绍使用初始化“调度数据库”配置调度中心配置“执行器项目”调度任务 介绍 首先我们介绍一下什么是xxl-job&#xff0c;根据官方定义&#xff0c;XXL-JOB是一个分布式任务调度平台&#xff0c;其核心设计目标是开发迅速、学习简单、轻量级、易扩展。现已开放源代码…

快速学会绘制Pyqt5中的所有图(下)

Pyqt5相关文章: 快速掌握Pyqt5的三种主窗口 快速掌握Pyqt5的2种弹簧 快速掌握Pyqt5的5种布局 快速弄懂Pyqt5的5种项目视图&#xff08;Item View&#xff09; 快速弄懂Pyqt5的4种项目部件&#xff08;Item Widget&#xff09; 快速掌握Pyqt5的6种按钮 快速掌握Pyqt5的10种容器&…

【Go自学版】02-goroutine

利用时间片分割进程&#xff0c;致使宏观上A,B,C同时执行&#xff08;并发&#xff09; CPU利用率包含了执行和切换&#xff0c;进程/线程的数量越多&#xff0c;切换成本也会增大 最大并行数&#xff1a;GOMAXPROCS work stealing: 偷其他队列的G hand off: 当前G1阻塞&#…

【Pytorch】Fizz Buzz

文章目录 1 数据编码2 网络搭建3 网络配置&#xff0c;训练4 结果预测5 翻车现场 学习参考来自&#xff1a; Fizz Buzz in Tensorflowhttps://github.com/wmn7/ML_Practice/tree/master/2019_06_10Fizz Buzz in Pytorch I need you to print the numbers from 1 to 100, excep…

数字化转型怎么才能做成功?_光点科技

数字化转型对于现代企业来说是一场必要的革命。它不仅仅是技术的更迭&#xff0c;更是企业战略、文化和运营方式全面升级的体现。一个成功的数字化转型能够使企业更具竞争力、更灵活应对市场变化&#xff0c;并最终实现业务增长和效率提升。那么&#xff0c;数字化转型怎么才能…

JVM常见垃圾回收器

串行垃圾回收器 Serial和Serial Old串行垃圾回收器&#xff0c;是指使用单线程进行垃圾回收&#xff0c;堆内存较小&#xff0c;适合个人电脑 Serial作用于新生代&#xff0c;采用复制算法 Serial Old作用于老年代&#xff0c;采用标记-整理算法 垃圾回收时&#xff0c;只有…
最新文章