首页/文章/ 详情

数学笔记:相似矩阵 · 奇异值分解 · 伪逆

2月前浏览886


一、相似矩阵(Similar Matrices)

1. 定义

设      和      都是      的方阵,若存在可逆矩阵     ,使得

 

则称      和     互为相似矩阵,记作     。


2. 直观含义

相似矩阵的本质是:

同一个线性变换,在不同基下所对应的矩阵表示。

  • 在张量分析中,一个        型线性变换张量,在不同坐标系(基)下写成的矩阵就是相似矩阵。
  • 相似 ≠ 数值相等,而是“几何本质相同,观测视角(基)不同”。

3. 相似矩阵的性质

既然相似矩阵描述同一个线性变换,那么以下量必然保持不变:

性质      
说明      
反身性      
       
对称性      
若         ,则              
传递性      
若         ,则              
特征值      
       
行列式      
            
(特征值之积)      
秩      
            
(无关特征向量数)      
迹      
            
(特征值之和)      
可逆性      
若          可逆,则          可逆,且              

只要认识到“相似矩阵 = 同一线性变换的不同矩阵表示”,上述性质都几乎是显然的。


4. 若尔当标准型(Jordan Canonical Form)

经典对角化公式:

 

此时     ,     为对角矩阵。

但并非所有矩阵都能对角化,例如:

 

该矩阵特征值     ,但代数重数 > 几何重数,线性无关特征向量不足。

若尔当指出:

任意方阵       ,即使不可对角化,也一定存在可逆矩阵       ,使得

 

其中        为若尔当标准型。

  •       由若干 若尔当块(Jordan Block) 组成;
  • 若尔当块特点:
    1. 主对角线元素相同(同一特征值);
    2. 次对角线(若存在)全为 1。
若尔当标准型  

二、奇异值分解(SVD)

在讨论伪逆之前,需要先理解矩阵的四个基本子空间的关系:

四个基本子空间  

1. 从可逆方阵到一般矩阵

  • 当        为 可逆方阵 时:

       
       

    表示行空间中的        经        映射到列空间;逆映射由        给出。

  • 当       不可逆 / 非方阵 时:

    • 零空间中的分量对输出无贡献;
    • 任何        都可分解为:      
    • 有效信息仅来自行空间分量       ,且      
      在列空间中唯一对应。

行空间与列空间维数相同(等于     ),在      作用下一一对应。


2. SVD 的构造思路

在行空间中选取一组标准正交基     ,其在      下的像:

 

可标准化为列空间的一组标准正交基。

写成矩阵形式:

 

即:

 

两边右乘     (注意      正交,    ),得 奇异值分解:

 

其中:

  •      :列空间的标准正交基(       或       );
  •      :行空间的标准正交基(       或       );
  •      :奇异值对角(或半对角)矩阵。

三、伪逆(Moore–Penrose Pseudoinverse)

1. 动机

即使      不可逆或是非方阵,行空间与列空间之间仍存在一一对应,因此可以定义一种“广义的逆”     ,把列空间向量映射回行空间。

由 SVD:

 

逆向操作自然是:

 

其中      将      中非零奇异值取倒数,并保持维度对齐。


2. 伪逆的作用

  • 左乘伪逆:

       
       

    效果是:保留向量在          行空间的投影,零空间分量被清零。

  • 右乘伪逆:

       
       

    效果是:保留向量在          列空间的投影,这正是最小二乘中的投影矩阵。


3. 特殊情形

列满秩(    ,左逆)

 

此时:

 

行满秩(    ,右逆)

 

4. 小结

  • 列满秩:有左逆,左乘得单位阵,右乘得列空间投影矩阵;
  • 行满秩:有右逆,右乘得单位阵,左乘得行空间投影矩阵;
  • 一般情况:伪逆        给出行空间 / 列空间的最优投影映射,是最小二乘解的统一表达式。

📌 一句话总结
相似矩阵是“换基不改变换”,SVD 是“任意矩阵的结构拆解”,伪逆是在不可逆世界中,对逆运算的最合理推广。



来源:我的博士日记
科普
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-07-30
最近编辑:2月前
此生君子意逍遥
博士 签名征集中
获赞 62粉丝 133文章 155课程 0
点赞
收藏
作者推荐

VPSC耦合元胞自动机模拟金属热压缩的再结晶行为

文章名称《A polycrystal plasticity based discontinuous dynamic recrystallization simulation method and its application to copper》DOI:10.1016/j.ijplas.2017.01.001传统动态再结晶模型往往更关注宏观流动应力或平均晶粒尺寸,例如通过临界应变、Avrami 方程或经验型晶粒尺寸方程描述软化行为。这类方法虽然计算方便,但很难回答一个关键问题:不同取向晶粒在变形过程中如何累积位错?哪些晶粒更容易发生再结晶?新生晶粒又如何影响后续织构演化?这篇文章正是针对这一问题,提出了一种基于 VPSC 的 DDRX 模拟方法。核心流程框架如下:模型的基本思想是:将每个晶粒看作嵌入均匀等效介质中的黏塑性夹杂,通过 VPSC 计算晶粒尺度的应力、应变率、滑移活动和取向演化。在此基础上,作者引入位错密度硬化模型,用 Kocks-Mecking 型方程描述位错累积与动态回复过程。当某个晶粒的位错密度达到临界值后,该晶粒具备动态再结晶形核条件。随后,模型通过形核概率判断是否生成新的 DRX 晶粒,并赋予新晶粒较低的初始位错密度和新的初始取向。文章中的另一个重要部分是晶粒长大。由于 VPSC 是平均场模型,不包含真实晶粒邻接关系,因此作者采用平均场晶界迁移假设:晶粒长大驱动力来自单个晶粒位错密度与整体平均位错密度之间的差值。低位错密度的新生 DRX 晶粒倾向于长大,而高位错密度的变形母晶逐渐被消耗。通过不断更新晶粒半径和体积分数权重,模型可以同时描述 DRX 体积分数、平均晶粒尺寸和织构演化。作者将该模型应用于 OFHC 铜的热变形模拟。结果表明,模型能够较好地预测流动应力曲线中的峰值和后续软化行为,也能反映动态再结晶导致的晶粒细化过程。更重要的是,模型能够进一步分析织构演化:塑性变形会使晶粒取向向特定稳定取向旋转,而动态再结晶新晶粒的不断形成会削弱原有织构强度,使整体织构趋于分散。这一点体现了多晶塑性模型相比传统经验型 DRX 模型的优势。这篇文章给我们的启示主要有三点。第一,动态再结晶建模不应只停留在宏观应力软化层面,而应尽可能联系到晶粒尺度的位错密度、取向演化和体积分数变化。第二,VPSC 这类平均场晶体塑性模型虽然没有空间分辨能力,但计算效率高,适合进行多参数、多温度、多应变率条件下的组织演化预测。第三,对于 FCC 金属材料,尤其是铜、铝合金以及镍基高温合金,位错密度驱动的 DDRX 模型具有较强的可扩展性。总体来看,这篇文章的最大贡献不是提出某一个复杂公式,而是建立了一条清晰的建模路线:以 VPSC 为力学骨架,以位错密度为内部变量,以临界位错密度控制形核,以储能差驱动晶粒长大,最终实现热变形过程中力学响应与微观组织演化的耦合预测。对于从事晶体塑性、动态再结晶和高温合金热加工模拟的研究者来说,这是一篇非常值得精读和复现的文章。使用作者提供的完整公式和计算流程。我们可以使用最新的vpsc8.0为基础,并在此上面完成二次开发,加入再结晶影响。在加入时可以将再结晶形核和生长做成相对独立的模块,这样可以不破坏原始结构。使用原始模型和参数,以及对应的温度(775K)和晶粒尺寸(13μm)进行简单的压缩模拟(ND方向压缩100%)。初始300个随机晶粒考虑DRX:初始随机织构:开启再结晶后的热变形模拟织构:开启再结晶后的热变形晶粒数量演化:开启再结晶后的热变形应力应变曲线:不考虑DRX:热变形模拟织构:热变形模拟应力应变曲线:来源:我的博士日记

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈