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

文章目录

  • 一、向量、矩阵范数与谱半径
    • 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]

特征值为 λ 1 = 11 + 2 10 \lambda_1 = 11+2\sqrt{10} λ1=11+210 , λ 2 = 11 − 2 10 \lambda_2 = 11-2 \sqrt{10} λ2=11210

∥ A ∥ 2 = λ max = 11 + 2 10 = 4.162277 \|A\|_2 = \sqrt{\lambda_{\text{max}}} = \sqrt{11+2\sqrt{10}} =4.162277 A2=λmax =11+210 =4.162277

  1. 谱半径:

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

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

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

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

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

数字图像处理(实践篇)二十二 使用opencv进行人脸、眼睛、嘴的检测

目录 1 xml文件 2 涉及的函数 3 实践 使用opencv进行人脸、眼睛、嘴的检测。 1 xml文件 方法① 下载 地址:https://github.com/opencv/opencv/tree/master/data/haarcascades 点击haarcascade_frontalface_default.xml文件 对着Raw右键,选择“链接…

和鲸科技携手深圳数据交易所,“数据+数据开发者生态”赋能人工智能产业发展

信息化时代,数据驱动决策的重要性日益凸显。通过利用数据可以深入了解市场需求、客户行为、竞争态势等关键信息,从而制定更为有效的战略和决策。围绕推动数据要素产业发展,近日,深圳数据交易所(以下简称“深数所”&…

ALTERNET STUDIO 9.1 Crack

ALTERNET STUDIO 9.1 发布 宣布 AlterNET Studio 9.1 版本今天上线。AlterNET Studio 9.0 是一个中期更新,重点是改进我们所有的组件库。 以下是 AlterNET Studio 9.1 的发布亮点: Roslyn C# 和 Visual Basic 解析器现在支持代码修复/代码重构。 代码修复…

JVM虚拟机系统性学习-对象存活判断算法、对象引用类型和垃圾清除算法

垃圾回收 在 JVM 中需要对没有被引用的对象,也就是垃圾对象进行垃圾回收 对象存活判断算法 判断对象存活有两种方式:引用计数法、可达性分析算法 引用计数法 引用计数法通过记录每个对象被引用的次数,例如对象 A 被引用 1 次&#xff0c…

docker mysql8 设置不区分大小写

docker安装Mysql8.0的坑之lower_case_table_names_docker mysql lower_case_table_names-CSDN博客https://blog.csdn.net/p793049488/article/details/108365929 docker run ‐di ‐‐nametensquare_mysql ‐p 33306:3306 ‐e MYSQL_ROOT_PASSWORD123456 mysql

广东佛山开房屋租赁发票

我是20223年12月办理的,给大家做个参考。 一、准备材料 (如果非房东本人办理,还需要房东签份授权书,多复印几份或者直接签多份,不然会被税务局收走) 废话不多说,直接上图。 二、线上预约 附个…

变量与方法面试题

char 型变量中能不能存储一个中文汉字,为什么? char 类型可以存储一个中文汉字,因为 Java 中使用的编码是 Unicode(不选择任何特定的编码,直接使用字符在字符集中的编号,这是统一的唯一方法)&a…

内存cache大量使用问题导致应用异常问题

概述 28s应用崩溃查看内存使用有大量cache。 分析 查看free 信息平时的确存在大量cache使用的情况查看dmes信息发现filesendserver崩溃 崩溃信息为系统调用 查看到page allocation failure:order 5 同时也看到系统内存使用情况 查看到系统实际还有部分内存为空闲内存&am…

如何部署自己的服务渲染页面为Pdf文档

前言 相信大家都觉得官方发布的文档生成模块https://docs.mendix.com/appstore/modules/document-generation/很有用,它能把Mendix页面像素级导出到Pdf文件中,这对于归档等业务非常有价值。但部署依赖公有云提供的渲染服务,而中国本土用户对…

<软考高项备考>《论文专题 - 4 论文题目和框架》

1 大数据分析论文的题目 2 论文写作的基础 2.1 10大知识域、49个过程、5大过程组 领域启动计划执行监控收尾整合管理1.制定项目章程2.制订项目管理计划3.指导与管理项目工作4.管理项目知识5.监控项目工作6.实施整体变更控制7.结束项目或阶段范围管理1.规划范围管理2.收集需求…

mysql数据库学习笔记(1)

今天开始学mysql数据库,为什么要学这个呢,因为数据库可结构化存储大量的数据信息,方便用户进行有效的检索和访问。数据库可有效地保持数据信息的一致性、完整性、降低数据冗余。数据库可满足应用的共享和安全方面的要求,把数据放在…

《Spring Cloud Alibaba 从入门到实战》分布式消息(事件)驱动

分布式消息(事件)驱动 1、简介 事件驱动架构(Event-driven 架构,简称 EDA)是软件设计领域内的一套程序设计模型。 这套模型的意义是所有的操作通过事件的发送/接收来完成。 传统软件设计 举个例子,比如一个订单的创建在传统软…

Vuex快速上手

一、Vuex 概述 目标:明确Vuex是什么,应用场景以及优势 1.是什么 Vuex 是一个 Vue 的 状态管理工具,状态就是数据。 大白话:Vuex 是一个插件,可以帮我们管理 Vue 通用的数据 (多组件共享的数据)。例如:购…

数字化手术麻醉临床信息系统源码,全方位实时数据调阅及采集,自动生成各种手术医疗文书

医院手术麻醉管理系统源码,自主版权,二次开发方便快捷 技术架构:PHP、 js 、mysql、laravel、vue2 手术麻醉临床信息管理系统是数字化手段应用于手术过程中的重要组成部分,用数字形式获取并存储手术相关信息,既便捷又…

【深度学习】注意力机制(二)

本文介绍一些注意力机制的实现,包括EA/MHSA/SK/DA/EPSA。 【深度学习】注意力机制(一) 【深度学习】注意力机制(三) 目录 一、EA(External Attention) 二、Multi Head Self Attention 三、…

使用wire重构商品微服务

一.wire简介 Wire 是一个轻巧的Golang依赖注入工具。它由Go Cloud团队开发,通过自动生成代码的方式在编译期完成依赖注入。 依赖注入是保持软件 “低耦合、易维护” 的重要设计准则之一。 此准则被广泛应用在各种开发平台之中,有很多与之相关的优秀工…

Vue3封装一个轮播图组件

先看效果 编写组件代码 CarouselChart.vue <template><div classimg-box><el-button clickpreviousImages v-ifprops.showBtn>←</el-button><div classimg><div styledisplay: flex;gap: 20px idmove><imgclassimg-item v-for(item…

C语言--clock()时间函数【详细介绍】

一.clock()时间函数介绍 在 C/C 中&#xff0c;clock() 函数通常用于处理和测量程序运行时间&#xff08;时钟时间&#xff09;。它是一种数据类型&#xff0c;表示 CPU 执行指定任务所耗费的“时钟计数”数量&#xff0c;单位为“时钟周期”。 这个函数通常包含在 time.h 头文…

每日一题,头歌平台c语言题目

任务描述 题目描述:输入一个字符串&#xff0c;输出反序后的字符串。 相关知识&#xff08;略&#xff09; 编程要求 请仔细阅读右侧代码&#xff0c;结合相关知识&#xff0c;在Begin-End区域内进行代码补充。 输入 一行字符 输出 逆序后的字符串 测试说明 样例输入&…

大语言模型有什么意义?亚马逊训练自己的大语言模型有什么用?

近年来&#xff0c;大语言模型的崭露头角引起了广泛的关注&#xff0c;成为科技领域的一项重要突破。而在这个领域的巅峰之上&#xff0c;亚马逊云科技一直致力于推动人工智能的发展。那么&#xff0c;作为一家全球科技巨头&#xff0c;亚马逊为何会如此注重大语言模型的研发与…