NumPy 教程

NumPy 首页 NumPy 简介 NumPy 环境

数组

NumPy Ndarray 对象 NumPy 数据类型

创建和操作数组

NumPy 数组创建例程 NumPy 数组操作 NumPy 基于现有数据创建数组 NumPy 基于数值范围创建数组 NumPy 数组迭代 NumPy 重塑数组 NumPy 连接数组 NumPy 堆叠数组 NumPy 拆分数组 NumPy 展平数组 NumPy 转置数组

索引与切片

NumPy 索引 &切片 NumPy 索引 NumPy 切片 NumPy 高级索引 NumPy 高级索引 NumPy 字段访问 NumPy 使用布尔数组切片

数组属性与操作

NumPy 数组属性 NumPy 数组形状 NumPy 数组大小 NumPy 数组步长 NumPy 数组元素大小 NumPy 广播 NumPy 算术运算 NumPy 数组加法 NumPy 数组减法 NumPy 数组乘法 NumPy 数组除法

高级数组运算

NumPy 交换数组的轴 NumPy 字节交换 NumPy 副本和视图 NumPy 逐元素数组比较 NumPy 过滤数组 NumPy 连接数组 NumPy 排序、搜索& 计数函数 NumPy 搜索数组 NumPy 数组并集 NumPy 查找唯一行 NumPy 创建日期时间数组 NumPy 二元运算符 NumPy 字符串函数 NumPy 矩阵库 NumPy 线性代数 NumPy Matplotlib NumPy 使用 Matplotlib 绘制直方图

排序和高级操作

NumPy 数组排序 NumPy 沿轴排序 NumPy 使用花式索引进行排序 NumPy 结构化数组 NumPy 创建结构化数组 NumPy 操作结构化数组 NumPy 记录数组 NumPy 加载数组 NumPy 保存数组 NumPy 将值附加到数组 NumPy 交换列数组 NumPy 将轴插入数组

处理缺失数据

NumPy 处理缺失数据 NumPy 识别缺失值 NumPy 移除缺失数据 NumPy 插补缺失值数据

性能优化

NumPy 使用数组进行性能优化

线性代数

NumPy 线性代数 NumPy 矩阵库 NumPy 矩阵加法 NumPy 矩阵减法 NumPy 矩阵乘法 NumPy 逐元素矩阵运算 NumPy 点积 NumPy 矩阵求逆 NumPy 行列式计算 NumPy 特征值 NumPy 特征向量 NumPy 奇异值分解 NumPy 求解线性方程 NumPy 矩阵范数

元素级矩阵运算

NumPy 总和 NumPy 平均值 NumPy 中位数 NumPy 最小值 NumPy 最大值

集合运算

NumPy 唯一元素 NumPy 交集 NumPy 并集 NumPy 差集

随机数生成

NumPy 随机数生成器 NumPy 排列和重排 NumPy 均匀分布 NumPy 正态分布 NumPy 二项分布 NumPy 泊松分布 NumPy 指数分布 NumPy 瑞利分布 NumPy 逻辑分布 NumPy 帕累托分布 NumPy 使用 Seaborn 可视化分布 NumPy 多项分布 NumPy 卡方分布 NumPy Zipf 分布

文件输入 &输出

NumPy 使用 NumPy 进行 I/O NumPy 从文件读取数据 NumPy 将数据写入文件 NumPy 支持的文件格式

数学函数

NumPy 数学函数 NumPy 三角函数 NumPy 指数函数 NumPy 对数函数 NumPy 双曲函数 NumPy 舍入函数

傅里叶变换

NumPy 离散傅里叶变换 (DFT) NumPy 快速傅里叶变换 (FFT) NumPy 逆傅里叶变换 NumPy 傅里叶级数和变换 NumPy 信号处理应用 NumPy 卷积

多项式

NumPy 多项式表示 NumPy 多项式运算 NumPy 求多项式的根 NumPy 求多项式的根

统计

NumPy 统计函数 NumPy 描述性统计

日期时间函数

NumPy 日期和时间基础知识 NumPy 表示日期和时间 NumPy 日期和时间运算 NumPy 使用日期时间进行索引 NumPy 时区处理 NumPy 时间序列分析 NumPy 处理时间增量 NumPy 闰秒处理 NumPy 矢量化日期时间运算

ufunc

NumPy ufunc 简介 NumPy 创建通用函数 (ufunc) NumPy 算术通用函数 (ufunc) NumPy 小数舍入 ufunc NumPy 对数通用函数(ufunc) NumPy 求和通用函数 (ufunc) NumPy 乘积通用函数 (ufunc) NumPy 差分通用函数 (ufunc) NumPy 寻找最小公倍数 (LCM) NumPy 寻找最大公约数 (GCD) NumPy 三角函数 (ufunc) NumPy 双曲线 (ufunc) NumPy 集合运算(ufunc)

实用资源

NumPy 快速指南 NumPy 备忘单


NumPy - 矩阵范数

什么是矩阵范数?

矩阵范数是一个将非负数赋给矩阵的函数。它衡量矩阵的大小。

通常,矩阵范数用于量化矩阵的大小,并且在涉及矩阵方程的问题中起着重要作用,例如求解线性方程组或执行矩阵因式分解。

矩阵范数的常见类型

矩阵范数有几种类型,但最常用的是 -

  • Frobenius 范数
  • 1-范数
  • 无穷范数
  • 2-范数(谱范数)

Frobenius 范数

Frobenius 范数是最简单、最常用的矩阵范数之一。它被定义为矩阵元素绝对值平方和的平方根。数学上,它表示为 −

‖A‖F = √(i=1 j=1 |aij|2)

其中 A 为矩阵,aij 为矩阵元素。Frobenius 范数等价于矩阵作为向量时的 L2 范数。

1-范数

矩阵的 1-范数(也称为最大列和范数)定义为列和的绝对值之和。从数学上讲,它表示为 −

‖A‖1 = maxj i=1 |aij|

简单来说,1范数是矩阵任意列元素绝对值之和的最大值。

无穷范数

矩阵的无穷范数(也称为最大行和范数)定义为行和的绝对值之和的最大值。从数学上讲,它表示为 −

‖A‖ = maxi j=1 |aij|

无穷范数给出了矩阵中任意一行元素绝对值之和的最大值。

2-范数(谱范数)

矩阵的2-范数(也称为谱范数)定义为矩阵的最大奇异值。它衡量矩阵应用于向量时的最大拉伸因子。2-范数由下式给出:-

‖A‖2 = max(A)

其中,max(A)是矩阵A的最大奇异值。在这种情况下,2-范数与矩阵的奇异值相关,可以使用奇异值分解 (SVD) 计算。

NumPy 中的矩阵范数

NumPy 提供了用于计算各种矩阵范数的函数。numpy.linalg.norm() 函数可用于计算大多数常见的矩阵范数。让我们探索如何将此函数用于不同类型的矩阵范数。

使用 NumPy 计算 Frobenius 范数

要使用 NumPy 计算 Frobenius 范数,我们使用带有参数 ord='fro'numpy.linalg.norm() 函数。

示例

以下示例中,矩阵 A 的弗罗贝尼乌斯范数是通过对矩阵中所有元素的平方和取平方根来计算的 -

import numpy as np

# 定义矩阵 A
A = np.array([[1, 2], [3, 4]])

# 计算矩阵的弗罗贝尼乌斯范数
frobenius_norm = np.linalg.norm(A, 'fro')

print("A 的弗罗贝尼乌斯范数:", frobenius_norm)

以下是得到的输出 -

A 的弗罗贝尼乌斯范数: 5.477225575051661

使用 NumPy 计算 1-范数

为了计算 1-范数,我们使用 numpy.linalg.norm() 函数,并传入参数 ord=1。矩阵的 1-范数是矩阵任意列元素绝对值之和的最大值。

示例

在本例中,列和分别为 4 和 6,因此 1-范数为 6 −

import numpy as np

# 定义矩阵 A
A = np.array([[1, 2], [3, 4]])

# 计算矩阵的 1-范数
one_norm = np.linalg.norm(A, 1)

print("A 的 1-范数:", one_norm)

以下是得到的输出 −

A 的 1-范数:6.0

使用 NumPy 计算无穷范数

为了计算无穷范数,我们使用 numpy.linalg.norm() 函数,并传入参数 ord=np.inf。矩阵的无穷范数是矩阵任意行元素绝对值之和的最大值。

示例

本例中,行和为 3 和 7,因此无穷范数为 7 −

import numpy as np

# 定义矩阵 A
A = np.array([[1, 2], [3, 4]])

# 计算矩阵的无穷范数
infinity_norm = np.linalg.norm(A, np.inf)

print("A 的无穷范数:", infinity_norm)

以下是得到的输出 −

A 的无穷范数: 7.0

使用 NumPy 计算 2-范数

为了计算 2-范数(谱范数),我们使用 numpy.linalg.norm() 函数,并传入参数 ord=2。矩阵的 2-范数(谱范数)是矩阵的最大奇异值,它衡量矩阵应用于向量时的最大拉伸因子。

示例

以下是在 NumPy 中计算 2-范数的示例 -

import numpy as np

# 定义矩阵 A
A = np.array([[1, 2], [3, 4]])

# 计算矩阵的 2-范数
two_norm = np.linalg.norm(A, 2)

print("A 的 2-范数(谱范数):", two_norm)

以下是得到的输出 -

2-范数A 的谱范数:5.464985704219043

矩阵范数的应用

矩阵范数在数值分析、机器学习、优化等领域有着广泛的实际应用 -

  • 数值稳定性:矩阵范数用于分析数值算法的稳定性,尤其是在求解线性系统或执行矩阵分解时。
  • 机器学习:在机器学习中,矩阵范数通常用于正则化模型并防止过拟合。例如,L2 正则化使用 Frobenius 范数来惩罚模型中的大权重。
  • 优化:矩阵范数用于度量优化问题中与期望解的误差或偏差。
  • 信号处理:在信号处理中,矩阵范数用于度量信号和滤波器的"能量"或幅度。