当前位置: 首页 > article >正文

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

1. 定义

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

2. 性质

3. 作用

(1)Ax=b直接法求解

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

(2)特征值求解

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

M必须得是

ModeOperatorM
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> =x^HMy
  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 矩阵可能不是正定或半正定的,这并不意味着广义特征值问题就无法求解。这些情况下,问题可能更加复杂,需要特别的数值方法来处理可能出现的数值不稳定性或解的不确定性。


http://www.kler.cn/a/230931.html

相关文章:

  • STM32 硬件随机数发生器(RNG)
  • 设计模式1-访问者模式
  • 国密SM2算法进行数据的加密、签名和验签、解密
  • 图解支付-金融级密钥管理系统:构建支付系统的安全基石
  • Springboot整合JUnit5框架
  • 【学网攻】 第(23)节 -- PPP协议
  • Redis系列——Lua脚本和redis事务的应用
  • swift结算体系
  • 获取旁站 / C 段:第三方网站(附链接)
  • 问题:0xc8前面加(byte) #人工智能#学习方法的原因是因为0xc8大于??????????? 。 #微信#其他#微信
  • 【数据结构与算法】堆 / 堆排序 / TopK问题(Heap)
  • 洛谷C++简单题小练习day9—[AHOI2017]寻找探监点
  • 单片机学习笔记---串口通信(2)
  • 什么是UI设计?
  • Mapbox Vision SDK 介绍
  • 专业排版设计软件:QuarkXPress 2024 for mac中文激活版
  • 多个Nginx虚拟主机部署脚本
  • 【Linux】基于管道进行进程间通信
  • 【RK3288 Android10 C30 支持sim卡拔掉不弹窗,及热插拔】
  • 代码随想录算法训练营第13天—二叉树02 | ● *层序遍历(对应10道题) ● *226.翻转二叉树 ● 101.对称二叉树