矩阵的正定(positive definite)性质的作用

1. 定义

注意,本文中正定和半正定矩阵不要求是对称或Hermite的。

2. 性质

3. 作用

(1)Ax=b直接法求解

|----------|--------------|
| cholesky | 实对称正定矩阵求解 |
| cholesky | 复共轭对称正定矩阵求解 |
| LDL | 实对称非正定矩阵求解 |
| LDL | 复共轭对称非正定矩阵求解 |
| LDL | 复对称矩阵求解 |
| LU | 实非对称矩阵求解 |
| LU | 复非对称矩阵求解 |

(2)特征值求解

在ARPACK(隐式重启Arnoldi算法)中,对K*x=lambda*M*x该广义特征值问题

M必须得是

|------------------|-------------------------------------------|----------------------|
| Mode | Operator | M |
| Shift | OP = inv[M]*K, and B = M. | 对称-正定 或Hemitian-正定 |
| Shift-and-invert | OP = (inv[K - sigma*M])*M, and B = M. | 对称-半正定 或Hemitian-半正定 |

注释:

  1. OP:operator,表示Arnoldi过程中与向量作用的算子,用户需要提供矩阵向量乘积w ← OPv
  2. M-inner product: <x ,y > =
  3. M-orthogonal: x, y称为M-orthogonal若<x ,y> = 0
  4. B: 用来定义M-inner product中的矩阵,用户需要提供矩阵向量乘积w ← Mv

Slepc有提及,若M不是正定也不是半正定的话,可以用EPS_GHIEP求解。

特征值中,正定或半正定性质对于 M 矩阵来说是一个优良属性,因为它确保了问题的物理可解性和数值计算的稳定性。例如,在结构动力学中,M 作为质量矩阵时,其正定性意味着系统的质量分布是非负的,这是物理上合理的。正定或半正定的 M 矩阵也有助于保证广义特征值问题解的良好性质,如确保所有特征值是实数且特征向量是良定义的。

然而,在某些情况下,M 矩阵可能不是正定或半正定的,这并不意味着广义特征值问题就无法求解。这些情况下,问题可能更加复杂,需要特别的数值方法来处理可能出现的数值不稳定性或解的不确定性。

相关推荐
我是大聪明.3 小时前
CUDA矩阵乘法优化:共享内存分块与Warp级执行机制深度解析
人工智能·深度学习·线性代数·机器学习·矩阵
做cv的小昊9 小时前
【TJU】研究生应用统计学课程笔记(6)——第二章 参数估计(2.4 区间估计)
人工智能·笔记·线性代数·算法·机器学习·数学建模·概率论
EnCi Zheng11 小时前
02a-什么是矩阵
线性代数·矩阵
AI科技星12 小时前
《全域数学》第一部:数术本源·第二卷《算术原本》之十四附录(二)全域数学体系下三大数论猜想的本源推演与哲学阐释【乖乖数学】
人工智能·线性代数·机器学习·量子计算·agi
有为少年1 天前
从概率估计到“LLM 训练是有损压缩”
人工智能·线性代数·机器学习·计算机视觉·矩阵
风落无尘1 天前
第二章《概率与生存》完整学习资料
人工智能·矩阵·概率论
大江东去浪淘尽千古风流人物1 天前
【UV-SLAM】eLSD/LBD 数据维度 UV-SLAM吸收借鉴
数据库·线性代数·oracle·矩阵·uv·augmented reality
风落无尘1 天前
《智能重生:从垃圾堆到AI工程师》——第三章 矩阵与防线
人工智能·线性代数·矩阵
玛丽莲茼蒿1 天前
Leetcode hot100 螺旋矩阵【中等】
算法·leetcode·矩阵
个微管理2 天前
小红书新规深度拆解:从被封到破局,2026年矩阵号生存手册
大数据·人工智能·矩阵