码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 矩阵分析与应用


    本文学习来源是《矩阵分析与应用》

    矩阵的二次型

    任意一个正方矩阵A的二次型x^{H}Ax是一个实标量,以实矩阵为例:
     

    x^{T}Ax =\begin{bmatrix} x_{1} ,x_{2} , x_{3} \end{bmatrix}\begin{bmatrix} 1 & 4&2 \\ -1& 7&5 \\ -1&6 & 3 \end{bmatrix}\begin{bmatrix} x_{1}\\ x_{2}\\ x_{3} \end{bmatrix} \\ =x_{1}^{2}-x_{2}x_{1}-x_{3}x_{1}+4x_{1}x_{2}+7x_{2}^{2}+6x_{3}x_{2}+2x_{1}x_{3}+5x_{2}x_{3}+3x_{3}^{2}\\=x_{1}^{2}+7x_{2}^{2}+3x_{3}^{2}+3x_{1}x_{2}+x_{1}x_{3}+11x_{2}x_{3}

    这是变量x的二次型函数,我们称x^{T}Ax为矩阵A的二次型。

    推而广之,若x=[x_{1},x_{2},.. .,x_{n}]^{T},且n\times n矩阵A的元素为a_{ij},则二次型

    x^{T}Ax=\sum_{i=1}^{n}\sum_{j=1}^{n}x_{i}x_{j}a_{ij} \\=\sum_{i=1}^{n}a_{ii} x_{i}^{2}+\sum_{i=1,i\neq j}^{n}\sum_{j=1}^{n}a_{ij}x_{i}x_{j}\\=\sum_{i=1}^{n}a_{ii} x_{i}^{2}+\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}(a_{ij}+a_{ji})x_{i}x_{j}

    若两个矩阵A,B满足对角线相等,且满足a_{ij}+a_{ji}=b_{ij}+b_{ji},则这两个矩阵的二次型就相等。

    如:

    A=\begin{bmatrix} 1 & 4&2 \\ -1& 7&5 \\ -1&6 & 3 \end{bmatrix}B=\begin{bmatrix} 1 & -1&-1 \\ 4& 7&6 \\ 2&5 & 3 \end{bmatrix}C=\begin{bmatrix} 1 & 114&52\\ -111& 7&2 \\ -51&9 & 3 \end{bmatrix}

    具有相同的二次型。即:

    x^{T}Ax=x^{T}Bx=x^{T}Cx\\=x_{1}^{2}+7x_{2}^{2}+3x_{3}^{2}+3x_{1}x_{2}+x_{1}x_{3}+11x_{2}x_{3}

    即对于任何一个二次型函数

    f\left ( x_{1},x_{2},.. .,x_{n} \right )=\sum_{i=1}^{n}\alpha _{ii}x_{i}^{2}+\sum_{i=1,i \neq j}^{n}\sum_{j=1}^{n}\alpha _{ij}x_{i}x_{j}

    存在许多矩阵A,它们的二次型x^{T}Ax=f\left ( x_{1},x_{2},.. .,x_{n} \right )相同,但是,只有一个唯一的对称矩阵A满足x^{T}Ax=f\left ( x_{1},x_{2},.. .,x_{n} \right ),其元素为a_{ii}=\alpha _{ii}和a_{ij}=a_{ji}=\frac{1}{2}(\alpha _{ij}+\alpha _{ji}),其中,i=1,2,.. .,n,j=1,2,.. .,n,i\neq j。因此在讨论矩阵A的二次型时,通常都假定A为实对称矩阵或复共轭对称矩阵。

    一个复共轭对称矩阵A称为:

    • 正定矩阵,若二次型x^{H}Ax> 0,\forall x \neq 0
    • 半正定矩阵,若二次型x^{H}Ax\geqslant 0,\forall x \neq 0
    • 负定矩阵,若二次型x^{H}Ax< 0,\forall x \neq 0
    • 半负定矩阵,若二次型x^{H}Ax\leqslant 0,\forall x \neq 0
    • 不定矩阵,若二次型x^{H}Ax既可能取正值,也可能取负值。

    例如:实对称矩阵

    R=\begin{bmatrix} 3 & -1&0 \\ -1 & 3& -1\\ 0&-1 & 3 \end{bmatrix}

    是正定的,因为二次型x^{H}Rx=2x_{1}^{2}+x_{2}^{2}+2x_{3}^{2}+\left ( x_{1}-x_{2} \right )^{2}+\left ( x_{2}-x_{3} \right )^{2}> 0,仅x_{1}=x_{2}=x_{3}=0时等式为0。

    这里继续学习一个非负矩阵的概念。

    设矩阵A_{m \times n}的元素为a_{ij},若

    a_{ij}\geqslant 0, \forall i=1,2,.. .,m,j=1,2,.. .,n

    即A的所有元素都是非负的,则称A为非负矩阵。

    若A的所有元素都是正的,则称A为正矩阵。

    非负矩阵和正矩阵分别用符号A\geqslant 0和A> 0简记,以下是正定矩阵和正矩阵之间的区别:

    正定矩阵一定是正方的矩阵,但是正矩阵可以是非正方的。

    正定矩阵A的定义为:x^{H}Ax> 0,\forall x \neq 0,正矩阵的定义为:a_{ij} > 0

    正定矩阵的符号为:x^{H}Ax>0,正矩阵的符号为:A>0

    矩阵的迹

    n*n矩阵A的对角元素之和称为A的迹,非正方矩阵没有迹。迹记作tr(A),即:
    tr(A)=a_{11}+a_{22}+.. .+a_{nn}=\sum_{i=1}^{n}a_{ii}

    以下是关于迹的等式:

    1. 若A和B均为n*n的矩阵,则tr(A\pm B)=tr(A) \pm tr(B)
    2. 若c是一个复常数或者实常数,则tr(cA)=ctr(A)
    3. 若A和B均为n*n的矩阵,且c_{1}和c_{2}为常数,则tr(c_{1}A\pm c_{2}B)=c_{1}tr(A) \pm c_{2}tr(B)
    4. 矩阵A的转置,复数共轭的迹为:tr(A^{T})=tr(A),tr(A^{*})=[tr(A)]^{*}
    5. 若A为m*n矩阵,B为n*m矩阵,则tr(AB)=tr(BA)
    6. 若A,B均是一个m*m矩阵,且B非奇异,则tr(BAB^{-1})=tr(B^{-1}AB)=tr(A)
    7. 若A是一个m*n的矩阵,则tr(A^{H}A)=0\Leftrightarrow A=O_{m \times n}
    8. 迹等于特征值之和,即tr(A)=\lambda _{1}+\lambda _{2}+.. .+\lambda _{n}
    9. tr(A^{k})=\sum_{i=1}^{n}\lambda _{i}^{k}
  • 相关阅读:
    Unity --- 场景/场景管理类 与 异步场景加载法
    SwiftUI SQLite教程之带有历史的搜索栏List App (教程含完整代码)
    【前端知识】Three 学习日志(四)—— 相机控件
    【python运维脚本实践】python实践篇之使用Python处理有序文件数据的多线程实例
    Flutter dio Http请求之Cookie管理
    82.(cesium之家)cesium点在3d模型上运动
    ChatGPT推出新“朗读”功能 支持多语言与声音;使用大型语言模型增强分类数据集
    步进电机驱动板在电机运行中起到什么作用?
    CPU使用率和负载区别及分析
    MyBatis一对多关联查询
  • 原文地址:https://blog.csdn.net/weixin_46007132/article/details/125591473
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号