矩阵的正定(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.hqwc.cn/news/457606.html

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

相关文章

[晓理紫]AI专属会议截稿时间订阅

AI专属会议截稿时间订阅 关注{晓理紫|小李子}&#xff0c;每日更新最新AI专属会议信息&#xff0c;如感兴趣&#xff0c;请转发给有需要的同学&#xff0c;谢谢支持&#xff01;&#xff01; 如果你感觉对你有所帮助&#xff0c;请关注我&#xff0c;每日准时为你推送最新AI专属…

pycharm deployment 灰色 一直无法点击

我的development的配置如下&#xff0c;我看了很多教程一直不知道为什么一直是灰色的&#xff0c; 文件夹配置&#xff1a; 如果你这里 Autodect&#xff0c;那么你Mapping 的文件夹应该是应该省略这个前缀的&#xff0c;例如我下面&#xff0c;我应该将本地文件夹映射到/home…

【经典项目】Java实现打地鼠小游戏(附源码)

一、游戏回顾 打地鼠游戏是一款简单而有趣的反应游戏。游戏中&#xff0c;你需要在地洞中出现的地鼠出现时迅速点击它们&#xff0c;以获得分数。以下是一般的打地鼠游戏玩法介绍&#xff1a; 准备阶段&#xff1a;游戏开始时&#xff0c;你会看到一个由多个地洞组成的游戏界面…

用python编写爬虫,爬取房产信息

题目 报告要求 工程报告链接放在这里 https://download.csdn.net/download/Samature/88816284使用 1.安装jupyter notebook 2.用jupyter notebook打开工程里的ipynb文件&#xff0c;再run all就行 注意事项 可能遇到的bug 暂无&#xff0c;有的话私信我

C:\Users\ShuYixiao>mysql ‘mysql‘ 不是内部或外部命令,也不是可运行的程序 或批处理文件

这个错误信息表明 mysql 命令在你的系统中没有被识别。这通常意味着 MySQL 没有安装在你的电脑上&#xff0c;或者它的可执行文件路径没有添加到系统的环境变量中。以下是一些解决这个问题的步骤&#xff1a; 确认 MySQL 是否已安装&#xff1a; 如果你还没有安装 MySQL&#x…

Oracle Vagrant Box 扩展根文件系统

需求 默认的Oracle Database 19c Vagrant Box的磁盘为34GB。 最近在做数据库升级实验&#xff0c;加之导入AWR dump数据&#xff0c;导致空间不够。 因此需要对磁盘进行扩容。 扩容方法1&#xff1a;预先扩容 此方法参考文档Vagrant, how to specify the disk size?。 指…

【ADI PLL 】时序小结

原文地址&#xff1a;https://ez.analog.com/cn/support-reference-library/rf/f/forum/95507/thread ADI的图一如既往的需要“推敲” ADI的所有锁相环产品控制接口均为三线串行控制接口。如上所示。要注意的是&#xff1a;在ADI的PLL产品中&#xff0c;大多数的时序图如图上所…

2024Node.js零基础教程(小白友好型),nodejs新手到高手,(六)NodeJS入门——http模块

047_http模块_获取请求行和请求头 hello&#xff0c;大家好&#xff0c;那第二节我们来介绍一下如何在这个服务当中来提取 HTT 请求报文的相关内容。首先先说一下关于报文的提取的方法&#xff0c;我在这个文档当中都已经记录好了&#xff0c;方便大家后续做一个快速的查阅。 …

腾讯云游戏联机服务器配置价格表,4核16G/8核32G/4核32G/16核64G

2024年更新腾讯云游戏联机服务器配置价格表&#xff0c;可用于搭建幻兽帕鲁、雾锁王国等游戏服务器&#xff0c;游戏服务器配置可选4核16G12M、8核32G22M、4核32G10M、16核64G35M、4核16G14M等配置&#xff0c;可以选择轻量应用服务器和云服务器CVM内存型MA3或标准型SA2实例&am…

Leetcode刷题笔记题解(C++):590. N 叉树的后序遍历

思路&#xff1a;类似于二叉树的排序&#xff0c;这里需要将子树进行依次递归遍历&#xff0c;前序遍历也与之类似 /* // Definition for a Node. class Node { public:int val;vector<Node*> children;Node() {}Node(int _val) {val _val;}Node(int _val, vector<N…

mysql 多数据源

依赖 <dependencies><!--mysql连接--><dependency><groupId>mysql</groupId><artifactId>mysql-connector-java</artifactId><scope>runtime</scope></dependency><!--多数据源--><dependency><g…

14.0 Zookeeper环球锁实现原理

全局锁是控制全局系统之间同步访问共享资源的一种方式。 下面介绍zookeeper如何实现全民锁&#xff0c;讲解他锁和共享锁两类全民锁。 排他锁 排他锁&#xff08;Exclusive Locks&#xff09;&#xff0c;又被称为写锁或独占锁&#xff0c;如果事务T1对数据对象O1加上排他锁…