主子式大于等于零的矩阵是半正定矩阵的证明方法之二
命题
主子式大于等于零的实对称矩阵是半正定矩阵。
即:对于任意的 n∈N+,n \in \mathbb N ^+, 对于任意一个实对称矩阵 An×n,A_{n \times n}, 若它的的任意一个主子式 ≥0,\ge 0, 则 X⊺AX≥0,∀X∈RnX ^ \intercal A X \ge 0, \forall X \in \mathbb R ^n 。
证明
n=1n = 1 时, 对于任意一个实对称矩阵 A1×1,A_{1 \times 1}, a11=|A|≥0,a_{11} = \left \vert A \right \vert \ge 0, 因此 a11x2≥0,∀x∈Ra_{11} x ^2 \ge 0, \forall x \in \mathbb R 。命题成立。
假设 nn 时命题成立。则 n+1n + 1 时,
记 A=(A1α⊺αan+1,n+1),A = \begin{pmatrix} A_1 & \alpha \\ \alpha ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix},
An+1×n+1A_{n + 1 \times n + 1} 的任意一个主子式 ≥0,\ge 0, 因此 an+1,n+1≥0a_{n + 1, n + 1} \ge 0 。
若 an+1,n+1=0,a_{n + 1, n + 1} = 0, 则对于任意的 i∈N,1≤i<n+1,i \in \mathbb N, 1 \le i \lt n + 1,
∣∣∣ai,ian+1,iai,n+1an+1,n+1∣∣∣=ai,ian+1,n+1−ai,n+1an+1,i=−ai,n+12, \begin{vmatrix} a_{i,i} & a_{i,n + 1} \\ a_{n + 1,i} & a_{n + 1,n + 1} \end{vmatrix} = a_{i,i} a_{n + 1,n + 1} - a_{i,n + 1} a_{n + 1,i} = -{a_{i,n + 1} }^2,
由于 ∣∣∣ai,ian+1,iai,n+1an+1,n+1∣∣∣≥0, \begin{vmatrix} a_{i,i} & a_{i,n + 1} \\ a_{n + 1,i} & a_{n + 1,n + 1} \end{vmatrix} \ge 0, 因此 ai,n+1=an+1,i=0, a_{i,n + 1} = a_{n + 1, i} =0, 于是 A=(A10⃗ ⊺0⃗ 0),A = \begin{pmatrix} A_1 & \vec 0 \\ {\vec 0} ^ {\intercal} & 0 \end{pmatrix},
令 X=(X1xn+1)∈Rn+1,X = \begin{pmatrix} X_1 \\ x_{n + 1} \end{pmatrix} \in \mathbb R ^{n + 1}, 则 X⊺AX=(X1⊺xn+1)(A10⃗ ⊺0⃗ 0)(X1xn+1)=X1⊺A1X1,X ^ {\intercal} A X = \begin{pmatrix} {X_1} ^ {\intercal} & x_{n + 1} \end{pmatrix} \begin{pmatrix} A_1 & \vec 0 \\ {\vec 0} ^ {\intercal} & 0 \end{pmatrix} \begin{pmatrix} X_1 \\ x_{n + 1} \end{pmatrix} = {X_1} ^ {\intercal} A_1 {X_1} ,
由归纳假设,X1⊺A1X1≥0, {X_1} ^ {\intercal} A_1 {X_1} \ge 0, 因此命题成立。
若 an+1,n+1>0,a_{n + 1, n + 1} \gt 0, 则
令 P=(E−an+1,n+1−1α⊺0⃗ 1),P = \begin{pmatrix} E & \vec 0 \\ - {a_{n + 1, n + 1} }^{-1} \alpha ^ {\intercal} & 1 \end{pmatrix},
则 P⊺APP ^ {\intercal} AP
=(E0⃗ ⊺−an+1,n+1−1α1)(A1α⊺αan+1,n+1)(E−an+1,n+1−1α⊺0⃗ 1)= \begin{pmatrix} E & - {a_{n + 1, n + 1} }^{-1} \alpha \\ \vec 0 ^ {\intercal} & 1 \end{pmatrix} \begin{pmatrix} A_1 & \alpha \\ \alpha ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} \begin{pmatrix} E & \vec 0 \\ - {a_{n + 1, n + 1} }^{-1} \alpha ^ {\intercal} & 1 \end{pmatrix}
=(A1−an+1,n+1−1αα⊺α⊺0⃗ an+1,n+1)(E−an+1,n+1−1α⊺0⃗ 1)= \begin{pmatrix} A_1- {a_{n + 1, n + 1} }^{-1} \alpha \alpha ^ {\intercal} & \vec 0 \\ \alpha ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} \begin{pmatrix} E & \vec 0 \\ - {a_{n + 1, n + 1} }^{-1} \alpha ^ {\intercal} & 1 \end{pmatrix}
=(A1−an+1,n+1−1αα⊺0⃗ ⊺0⃗ an+1,n+1)= \begin{pmatrix} A_1- {a_{n + 1, n + 1} }^{-1} \alpha \alpha ^ {\intercal} & \vec 0 \\ \vec 0 ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix}
令 B=A1−an+1,n+1−1αα⊺,B = A_1- {a_{n + 1, n + 1} }^{-1} \alpha \alpha ^ {\intercal} , 则 P⊺AP=(B0⃗ ⊺0⃗ an+1,n+1)P ^ {\intercal} AP = \begin{pmatrix} B & \vec 0 \\ \vec 0 ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix}
对于 BB 的任意一个主子式 B′k×k,B'_{k \times k}, 设 BB 中与 B′B' 对应的行号为 n1,⋯,nk,{n}_1, \cdots, {n}_k,
令 A1′=⎛⎝⎜⎜an1,n1⋮ank,n1⋯⋱⋯an1,nk⋮ank,nk⎞⎠⎟⎟k×k,α′=⎛⎝⎜⎜an1,n+1⋮ank,n+1⎞⎠⎟⎟k×1,{A_1} ' = \begin{pmatrix} a_{{n}_1, {n}_1} & \cdots & a_{{n}_1,{n}_k}\\ \vdots & \ddots & \vdots \\ a_{{n}_k, {n}_1} & \cdots & a_{{n}_k, {n}_k} \end{pmatrix} _{k \times k}, \alpha' = \begin{pmatrix} a_{{n}_1, n + 1} \\ \vdots \\ a_{{n}_k, n + 1} \end{pmatrix} _{k \times 1},
则 B′=A′1−an+1,n+1−1α′α′⊺,B' = A_1' - {a_{n + 1, n + 1} }^{-1} \alpha' \alpha' ^ {\intercal} ,
令 A′=(A′1α′⊺α′an+1,n+1),Q=(Ek×k−an+1,n+1−1α′⊺0⃗ 1),A' =\begin{pmatrix} A_1' & \alpha' \\ \alpha' ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} , Q = \begin{pmatrix} E_{k \times k} & \vec 0 \\ - {a_{n + 1, n + 1} }^{-1} \alpha' ^ {\intercal} & 1 \end{pmatrix}, 则
Q⊺A′QQ ^ {\intercal} A' Q
=(A′1−an+1,n+1−1α′α′⊺0⃗ ⊺0⃗ an+1,n+1)=(B′0⃗ ⊺0⃗ an+1,n+1)= \begin{pmatrix} A_1'- {a_{n + 1, n + 1} }^{-1} \alpha' \alpha' ^ \intercal & \vec 0 \\ \vec 0 ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} = \begin{pmatrix} B' & \vec 0 \\ \vec 0 ^ \intercal & a_{n + 1,n + 1} \end{pmatrix}
由于 A′=(A′1α′⊺α′an+1,n+1)=⎛⎝⎜⎜⎜⎜⎜an1,n1⋮ank,n1an+1,n1⋯⋱⋯⋯an1,nk⋮ank,nkan+1,nkan1,n+1⋮ank,n+1an+1,n+1⎞⎠⎟⎟⎟⎟⎟(k+1)×(k+1)A' = \begin{pmatrix} A_1' & \alpha' \\ \alpha' ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} = \begin{pmatrix} a_{{n}_1, {n}_1} & \cdots & a_{{n}_1,{n}_k} & a_{n_1, n + 1}\\ \vdots & \ddots & \vdots & \vdots\\ a_{{n}_k, {n}_1} & \cdots & a_{{n}_k, {n}_k} & a_{n_k, n + 1} \\ a_{{n + 1}, n_1} & \cdots & a_{{n + 1}, n_k} & a_{n + 1, n + 1} \end{pmatrix} _{(k + 1) \times (k + 1)}
是 AA 的一个主子式,因此 |A′|≥0,\vert A' \vert \ge 0,
又 |A′|=|Q|2|A′|=|Q⊺A′Q|=∣∣∣B′0⃗ ⊺0⃗ an+1,n+1∣∣∣=an+1,n+1|B′|,\vert A' \vert = \vert Q \vert ^2 \vert A' \vert = \vert Q ^ {\intercal} A' Q \vert = \begin{vmatrix} B' & \vec 0 \\ \vec 0 ^ \intercal & a_{n + 1,n + 1} \end{vmatrix} = a_{n + 1,n + 1} \vert B' \vert ,
因此 |B′|≥0, \vert B' \vert \ge 0, 故 BB 的任意一个主子式 ≥0, \ge 0, 由归纳假设,X⊺BX≥0,∀X∈RnX ^ \intercal B X \ge 0, \forall X \in \mathbb R ^{n}
因此对于任意的 X=(X1xn+1)∈Rn+1,X = \begin{pmatrix} X_1 \\ x_{n + 1} \end{pmatrix} \in \mathbb R ^{n + 1},
X⊺(B0⃗ ⊺0⃗ an+1,n+1)X=(X1⊺xn+1)(B0⃗ ⊺0⃗ an+1,n+1)(X1xn+1)=X1⊺BX1+an+1,n+1xn+12≥0,X ^ {\intercal} \begin{pmatrix} B & \vec 0 \\ \vec 0 ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix} X = \begin{pmatrix} {X_1} ^ {\intercal} & x_{n + 1} \end{pmatrix} \begin{pmatrix} B & \vec 0 \\ {\vec 0} ^ {\intercal} & a_{n + 1, n + 1} \end{pmatrix} \begin{pmatrix} X_1 \\ x_{n + 1} \end{pmatrix} = {X_1} ^ {\intercal} B {X_1} + a_{n + 1,n + 1} {x_{n + 1}} ^2 \ge 0,
由于 P⊺AP=(B0⃗ ⊺0⃗ an+1,n+1),P ^ {\intercal} AP = \begin{pmatrix} B & \vec 0 \\ \vec 0 ^ {\intercal} & a_{n + 1,n + 1} \end{pmatrix}, 因此 A<script type="math/tex" id="MathJax-Element-61">A</script> 也是半正定矩阵。
因此命题成立。
主子式大于等于零的矩阵是半正定矩阵的证明方法之二相关推荐
- 半正定矩阵的判定方法_线性代数30——正定矩阵和最小值
我们经常在判定一个函数是否有最小值时使用正定矩阵,正定矩阵和最小值有什么关系呢? 1 判断正定矩阵 给出一个矩阵: 有4个途径可以判定该矩阵是否是正定矩阵(注意这个矩阵的4个元素中有2个b,这是因为正 ...
- 判定(半)正定矩阵的特殊大于(等于)简写符号
判定(半)正定矩阵的特殊大于(等于)简写符号 这些符号的意义是: (1)特殊的大于号: M 特殊大于 0 (2)特殊的大于等于号:以此类推按照半正定矩阵的意义理解特殊符号的含义 posted on ...
- ker矩阵是什么意思_如何理解正定矩阵和半正定矩阵
乍看正定和半正定会被吓得虎躯一震,因为名字取得不知所以,所以老是很排斥去理解这个东西是干嘛用的,下面根据自己和结合别人的观点解释一下什么是正定矩阵(positive definite, PD) 和半正 ...
- 正定矩阵、负定矩阵、半正定矩阵、半负定矩阵
正定矩阵.负定矩阵.半正定矩阵.半负定矩阵 载▼ 1.正定矩阵 一个n×n的实对称矩阵M是正定的,当且仅当对于所有的非零实系数向量z,都有zTMz > 0.其中zT表示z的转置. 2.负 ...
- 线性代数学习笔记8-4:正定矩阵、二次型的几何意义、配方法与消元法的联系、最小二乘法与半正定矩阵A^T A
正定矩阵Positive definite matrice 之前说过,正定矩阵是一类特殊的对称矩阵: 正定矩阵满足对称矩阵的特性(特征值为实数并且拥有一套正交特征向量.正 / 负主元的数目等于正 / ...
- r语言中正定矩阵由于误差不正定_浅谈「正定矩阵」和「半正定矩阵」
在众多的机器学习模型中,线性代数的身影无处不在,当然,我们也会时常碰到线性代数中的正定矩阵和半正定矩阵.例如,多元正态分布的协方差矩阵要求是半正定的. --------------×--------- ...
- 「正定矩阵」和「半正定矩阵」
在众多的机器学习模型中,线性代数的身影无处不在,当然,我们也会时常碰到线性代数中的正定矩阵和半正定矩阵.例如,多元正态分布的协方差矩阵要求是半正定的. --------------×--------- ...
- 矩阵正定 matlab,matlab 错误使用 chol 矩阵必须为正定矩阵
警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇异值或缩放错误.结果可能不准确.RCOND = NaN. 警告: 矩阵为奇异值.接近奇 ...
- matlab矩阵半张量积,矩阵的半张量积_一个便捷的新工具.pdf
您所在位置:网站首页 > 海量文档  > 高等教育 > 微积分 矩阵的半张量积_一个便捷的新工具.pdf11页 本文档 ...
- Java(79):Character.isDigit() 判断输入的数据是否大于等于零整数
/*** @param str* @function 判断输入的数据是否是大于等于零的整数*/public static boolean isNumeric(String str) {for (int ...
最新文章
- Product Backlog:终极任务清单
- Python中的测试工具
- CUDA学习----sp, sm, thread, block, grid, warp概念
- Android EditText密码框的可见性设置
- 有关EnableQ在线问卷调查引擎的记者访谈记录
- html:(2):制作第一个网页和html和css的关系
- LeetCode OJ 160. Intersection of Two Linked Lists
- 去广州见了我大学老师标哥
- sql 生成一列1到10的数字_SQL 打印矩阵(三)
- 水培设备行业调研报告 - 市场现状分析与发展前景预测
- day17【前台】支付案例
- C++实现LRU(Least-Recently Used)缓存算法
- java 快速排序算法简单_Java 快速排序算法的简单说明及实现
- 主板声卡坏了会不会有电流声_不同领域适使用不同类型声卡,千万别选错了~~~...
- 五笔字根表识别码图_五笔输入法口诀(五笔字根表快速记忆图)
- 标书制作详细教程(零基础速成,助力公司中标)
- sql 触发器(Insert,update,delete)
- Sloth演示程序及源代码发布
- [日常] Go语言圣经前言
- 鼓捣中兴的TCM SSX44B
热门文章
- 杀软自己做 编写autorun病毒免疫工具
- linux tab 缩进,codemirror TAB 缩进问题记录
- 数据结构实验5-递归
- python哈夫曼编码注意_[Python]哈夫曼编码
- 使用Git上传文件至Github中
- 1060 Are They Equal (25 分)科学计数法,stl中string的各种函数用法
- nyoj Splits the string 749
- android arraymap排序,内存优化之ArrayMap、SparseArray、SparseIntArray
- 【POJ2826】An Easy Problem?!(线段相交+分情况讨论+精度)
- 【codeup22562】最长回文子串(dp基础题)