证明正定矩阵的充要条件:全部顺序主子式大于0

2023-11-16

定理: f = x T A x f = x^TAx f=xTAx 正定的充要条件是 A A A 的全部顺序主子式大于零。

必要性:即 f = x T A x f = x^TAx f=xTAx 正定 ⇒ \Rightarrow A A A 的全部顺序主子式大于零。

首先,由于 f = x T A x f = x^TAx f=xTAx 正定,即对称矩阵 A A A 是正定矩阵(有人会问为什么 A A A 是对称矩阵,这是研究二次型的必然要求,将二次多项式转换成 f = x T A x f = x^TAx f=xTAx的形式的时候,交叉项的系数可以全部是对称的),由于实对称矩阵必可正交对角化,所以矩阵 A A A 可以表示为:
A = Q T [ λ 1 ⋯ 0 ⋯ ⋮ λ 2 0 ⋯ 0 ⋮ ⋱ ⋯ 0 λ n ] Q A = Q^T \begin{bmatrix} \lambda_1 & \cdots& 0 & \cdots& \\ \vdots&\lambda_2&0 &\cdots& \\ 0 & & && \\ \vdots&&\ddots &\cdots& \\ 0&& &\lambda_n& \end{bmatrix} Q A=QT λ100λ200λn Q

所以 f = x T A x = x T Q T ∧ Q x f = x^TAx = x^TQ^T \wedge Q x f=xTAx=xTQTQx, 令 y = Q x y = Qx y=Qx f = λ 1 y 1 2 + ⋯ + λ n y n 2 f = \lambda_1 y_1^2+\cdots+\lambda_n y_n^2 f=λ1y12++λnyn2,对任意的 y = ( y 1 , . . . , y n ) y = (y_1,...,y_n) y=(y1,...,yn),都有 f f f 大于 0,由于 Q Q Q是正交矩阵,所以 x = ( x 1 , . . . , x n ) x= (x_1, ..., x_n) x=(x1,...,xn) 能够铺满整个n维空间,则 y = Q x = ( y 1 , . . . , y n ) y = Qx = (y_1,...,y_n) y=Qx=(y1,...,yn) 也能铺满整个n维空间。所以对任意非零 y y y 都有 f > 0 f >0 f>0 ,取 y i ≠ 0 , y j = 0 ( j ≠ i ) y_i\neq 0,y_j = 0(j\neq i) yi=0yj=0(j=i),能够得到 λ i y i 2 > 0 \lambda_iy_i^2>0 λiyi2>0,从而 λ i > 0 \lambda_i>0 λi>0,即正定矩阵的所有特征值全大于0,自然,行列式的值 ∣ A ∣ = λ 1 λ 2 . . . λ n > 0 |A|=\lambda_1\lambda_2...\lambda_n>0 A=λ1λ2...λn>0

A k A_k Ak 为 矩阵 A A A 左上角的 k k k 阶矩阵, k = 1 , 2 , ⋯   , n − 1 k = 1, 2, \cdots, n-1 k=1,2,,n1,由 f = x T A x f = x^TAx f=xTAx 对任意 x ≠ 0 x\neq0 x=0 成立,取 x = ( x 1 , x 2 , ⋯   , x k , 0 , ⋯   , 0 ) x=(x_1,x_2,\cdots,x_k,0,\cdots,0) x=(x1,x2,,xk,0,,0),并记 x ′ = ( x 1 , x 2 , ⋯   , x k ) x'=(x_1, x_2, \cdots, x_k) x=(x1,x2,,xk),则 x   T A x = x ′   T A k x ′ > 0 x\ ^TAx = x'\ ^TA_k x'>0 x TAx=x TAkx>0,我们可以得到 A k A_k Ak也是正定矩阵,自然有 ∣ A k ∣ > 0 |A_k|>0 Ak>0。从而,我们有正定矩阵 A A A 的顺序主子式全部大于0.

充分性: A A A 的全部顺序主子式大于零 ⇒ \Rightarrow f = x T A x f = x^TAx f=xTAx 正定 。

对于 1 阶矩阵,任意1 阶矩阵 A = [ a ] A = [a] A=[a] 的行列式为该矩阵的元素值,该矩阵的顺序主子式只有1个,即 ∣ A ∣ = a |A|=a A=a,该矩阵的行列式大于 0 时,自然有对任意 x x x x T A x = a x 2 > 0 x^TAx = a x^2>0 xTAx=ax2>0

假设:对任意的 n - 1 阶矩阵,结论“ A A A 的全部顺序主子式大于零 ⇒ \Rightarrow f = x T A x f = x^TAx f=xTAx 正定” 成立。

对于 n 阶矩阵,

x T A x = ( x 1 , x 2 , . . . , x n ) [ a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ a n 1 a n 2 ⋯ a n n ] ( x 1 x 2 ⋮ x n ) = a 11 x 1 2 + a 12 x 1 x 2 + a 13 x 1 x 3 + ⋯ + a 1 n x 1 x n a 21 x 2 x 1 + a 22 x 2 2 + a 23 x 2 x 3 + ⋯ + a 2 n x 2 x n ⋮ a n 1 x n x 1 + a n 2 x n x 2 + a n 3 x n x 3 + ⋯ + a n n x n 2 \begin{align*} x^TAx &= (x_1,x_2,...,x_n) \begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix} \left( \begin{array}{rrrr} x_1\\ x_2\\ \vdots \\ x_n \end{array} \right) \\ &= \begin{align*} a_{11}x_1^2 + a_{12}x_1x_2 + a_{13}x_1x_3+\cdots+a_{1n}x_1x_n\\ a_{21}x_2x_1 + a_{22}x_2^2 + a_{23}x_2x_3+\cdots+a_{2n}x_2x_n\\ \vdots \\ a_{n1}x_nx_1 + a_{n2}x_nx_2 + a_{n3}x_nx_3+\cdots + a_{nn}x_n^2 \end{align*} \end{align*} xTAx=(x1,x2,...,xn) a11a21an1a12a22an2a1na2nann x1x2xn =a11x12+a12x1x2+a13x1x3++a1nx1xna21x2x1+a22x22+a23x2x3++a2nx2xnan1xnx1+an2xnx2+an3xnx3++annxn2
x T A X x^TAX xTAX 表示为:

x T A x = 1 a 11 ( a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n ) 2 + ∑ i = 2 n ∑ j = 2 n b i j x i x j \begin{align*} x^TAx &= \frac{1}{a_{11}} ( a_{11}x_1 + a_{12}x_2+\cdots+a_{1n}x_n)^2+\sum_{i=2}^n\sum_{j=2}^n b_{ij} x_ix_j \end{align*} xTAx=a111(a11x1+a12x2++a1nxn)2+i=2nj=2nbijxixj

其中, b i j = a i j − a 1 i a 1 j a 11 b_{ij}= a_{ij} - \frac{a_{1i}a_{1j}}{a_{11}} bij=aija11a1ia1j,由于 A A A 是实对称矩阵,所以 a i j = a j i a_{ij} = a_{ji} aij=aji,所以 b j i = a j i − a 1 j a 1 i a 11 = a i j − a 1 i a 1 j a 11 = b i j b_{ji} = a_{ji} - \frac{a_{1j}a_{1i}}{a_{11}} = a_{ij} - \frac{a_{1i}a_{1j}}{a_{11}} =b_{ij} bji=ajia11a1ja1i=aija11a1ia1j=bij.

∑ i = 2 n ∑ j = 2 n b i j x i x j = ( x 2 , x 3 , ⋯   , x n ) [ b 22 b 23 ⋯ b 2 n b 32 b 33 ⋯ b 3 n ⋮ ⋱ ⋮ b n 2 b n 3 ⋯ b n n ] ( x 2 x 3 ⋮ x n ) \sum_{i=2}^n\sum_{j=2}^n b_{ij} x_ix_j = (x_2, x_3,\cdots,x_n) \begin{bmatrix} b_{22} & b_{23} & \cdots & b_{2n} \\ b_{32} & b_{33} & \cdots & b_{3n} \\ \vdots & & \ddots & \vdots \\ b_{n2} & b_{n3} & \cdots & b_{nn} \end{bmatrix} \begin{pmatrix} x_2\\ x_3\\ \vdots\\ x_n \end{pmatrix} i=2nj=2nbijxixj=(x2,x3,,xn) b22b32bn2b23b33bn3b2nb3nbnn x2x3xn
∑ i = 2 n ∑ j = 2 n b i j x i x j = x ′ T B x ′ \sum_{i=2}^n\sum_{j=2}^n b_{ij} x_ix_j = x'^TBx' i=2nj=2nbijxixj=xTBx,则 B 也是实对称矩阵。

记矩阵 A A A 左上角的 k 阶矩阵为 A k , k = 2 , 3 , ⋯   , n A_k,k=2,3,\cdots,n Ak,k=2,3,,n,则
∣ A k ∣ = ∣ a 11 a 12 a 13 ⋯ a 1 k a 21 a 22 a 23 ⋯ a 2 k a 31 a 32 a 33 ⋯ a 3 k ⋮ ⋱ ⋮ a k 1 a k 2 a k 3 ⋯ a k k ∣ = ∣ a 11 a 12 a 13 ⋯ a 1 k 0 a 22 − a 12 a 21 a 11 a 23 − a 13 a 21 a 11 ⋯ a 2 k − a 1 k a 21 a 11 0 a 32 − a 12 a 31 a 11 a 33 − a 13 a 31 a 11 ⋯ a 3 k − a 1 k a 31 a 11 ⋮ ⋱ 0 a k 2 − a 12 a k 1 a 11 a k 3 − a 13 a k 1 a 11 ⋯ a k k − a 1 k a k 1 a 11 ∣ = ∣ a 11 a 12 a 13 ⋯ a 1 k 0 b 22 b 23 ⋯ b 2 k 0 b 32 b 33 ⋯ b 3 k ⋮ ⋱ 0 b k 2 b k 3 ⋯ b k k ∣ = a 11 ∣ b 22 b 23 ⋯ b 2 k b 32 b 33 ⋯ b 3 k ⋮ ⋱ b k 2 b k 3 ⋯ b k k ∣ \begin{align*} |A_k| &= \begin{vmatrix} a_{11} & a_{12} & a_{13} & \cdots & a_{1k}\\ a_{21} & a_{22} & a_{23} & \cdots & a_{2k}\\ a_{31} & a_{32} & a_{33} & \cdots & a_{3k}\\ \vdots & & & \ddots & \vdots \\ a_{k1} & a_{k2} & a_{k3} & \cdots & a_{kk} \end{vmatrix}\\ & = \begin{vmatrix} a_{11} & a_{12} & a_{13} & \cdots & a_{1k}\\ 0 & a_{22}-\frac{a_{12}a_{21}}{a_{11}} & a_{23} -\frac{a_{13}a_{21}}{a_{11}}& \cdots & a_{2k}-\frac{a_{1k}a_{21}}{a_{11}}\\ 0 & a_{32}-\frac{a_{12}a_{31}}{a_{11}} & a_{33} -\frac{a_{13}a_{31}}{a_{11}}& \cdots & a_{3k}-\frac{a_{1k}a_{31}}{a_{11}}\\ \vdots & & & \ddots & \\ 0 & a_{k2}-\frac{a_{12}a_{k1}}{a_{11}} & a_{k3} -\frac{a_{13}a_{k1}}{a_{11}}& \cdots & a_{kk}-\frac{a_{1k}a_{k1}}{a_{11}}\\ \end{vmatrix}\\ &= \begin{vmatrix} a_{11} & a_{12} & a_{13} & \cdots & a_{1k}\\ 0 & b_{22} & b_{23} & \cdots & b_{2k} \\ 0 & b_{32} & b_{33} & \cdots & b_{3k} \\ \vdots & & & \ddots & \\ 0 & b_{k2} & b_{k3} & \cdots & b_{kk}\\ \end{vmatrix} \\ &= a_{11}\begin{vmatrix} b_{22} & b_{23} & \cdots & b_{2k} \\ b_{32} & b_{33} & \cdots & b_{3k} \\ \vdots & & \ddots & \\ b_{k2} & b_{k3} & \cdots & b_{kk}\\ \end{vmatrix} \\ \end{align*} Ak= a11a21a31ak1a12a22a32ak2a13a23a33ak3a1ka2ka3kakk = a11000a12a22a11a12a21a32a11a12a31ak2a11a12ak1a13a23a11a13a21a33a11a13a31ak3a11a13ak1a1ka2ka11a1ka21a3ka11a1ka31akka11a1kak1 = a11000a12b22b32bk2a13b23b33bk3a1kb2kb3kbkk =a11 b22b32bk2b23b33bk3b2kb3kbkk
A A A 的全部顺序主子式大于0, 有 a 11 > 0 , ∣ A k ∣ > 0 a_{11}>0, |A_k|>0 a11>0,Ak>0,可得 ∣ B k ∣ = ∣ A k ∣ a 11 > 0 |B_k| = \frac{|A_k|}{a_{11}}>0 Bk=a11Ak>0

由此可得实对称矩阵 B B B 的所有顺序主子式大于0。

根据假设,对于任意n-1阶矩阵,全部顺序主子式大于0能够得到该矩阵正定,所以 ∑ i = 2 n ∑ j = 2 n b i j x i x j = x ′ T B x ′ \sum_{i=2}^n\sum_{j=2}^n b_{ij} x_ix_j = x'^TBx' i=2nj=2nbijxixj=xTBx 对于任意 x ′ ≠ 0 x'\neq 0 x=0,都大于0. 则:

x T A x = 1 a 11 ( a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n ) 2 + ∑ i = 2 n ∑ j = 2 n b i j x i x j > 0 \begin{align*} x^TAx &= \frac{1}{a_{11}} ( a_{11}x_1 + a_{12}x_2+\cdots+a_{1n}x_n)^2+\sum_{i=2}^n\sum_{j=2}^n b_{ij} x_ix_j \\ &>0 \end{align*} xTAx=a111(a11x1+a12x2++a1nxn)2+i=2nj=2nbijxixj>0

由此我们可以得到,在假设对于n-1阶矩阵结论成立的情况下,能够推出对于n阶矩阵结论仍然成立,又对于1阶矩阵结论显然是成立的。所以我们就能得到:

A A A 的全部顺序主子式大于零 ⇒ \Rightarrow f = x T A x f = x^TAx f=xTAx 正定 。

必要性和充分性都得到了证明,所以我们就证明了最开始的定理。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

证明正定矩阵的充要条件:全部顺序主子式大于0 的相关文章

  • vue添加水印踩坑

    介绍 前景 app页面添加水印展示 技术实现 watermark dom 完整代码 vue watermark 实现效果 功能描述 添加 删除 更新水印 引入 方式一 推荐 方便拓展 在index html引入相关文件 方式二 npm包引入

随机推荐

  • java byte[] 学习总结

    最近在学习netty 突然发现自己对字符数组是那么的陌生 吓死宝宝了 然后各种学习 然后测试 终于会用一些了 下线的都是本人的学习笔记 byte表是字符 一个字节 8位 可以组成2 8 256中不同数字 byte存值范围 128 127 1
  • pytorch基本使用_02

    import numpy as np import torch 从numpy引入tensor a np array 2 3 3 print torch from numpy a tensor 2 0000 3 3000 dtype torc
  • java线上CPU100%如何排查

    定位耗费CPU的进程 top c 就可以显示进程列表 然后输入P 按照cpu使用率排序 你会看到类似下面的东西 2 定位耗费CPU的线程 top Hp 1500 就是输入那个进程id就好了 然后输入P 按照cpu使用率排序 你会看到类型下面
  • 单片机c语言数码管显示0到9,单片机如何让8个数码管同时流水显示0到9,大家帮我看看!...

    按你的要求修改如下 include reg52 h 此文件中定义了单片机的一些特殊功能寄存器 typedef unsigned int u16 对数据类型进行声明定义 typedef unsigned char u8 sbit LSA P2
  • Java 网络编程UDP协议之发送数据和接收数据的详解

    博主前些天发现了一个巨牛的人工智能学习网站 通俗易懂 风趣幽默 忍不住也分享一下给大家 点击跳转到网站 UDP协议 用户数据报协议 User Datagram Protocol UDP是无连接通信协议 即在数据传输时 数据的发送端和接收端不
  • 《信号与系统》4.10.2工频干扰的滤除

    平台 版本 Multisim14 1 参考书籍 信号与系统 4 10 2工频干扰的滤除 工程上 滤除工频干扰比较常用的电路是无源双T陷波滤波器 图示双T的无源陷波滤波器电路 陷波器是某一小频率范围内的带阻滤波器 陷波器的一个常见的应用是滤除
  • Seaborn入门详细教程

    作者 luanhz 来源 小数志 Seaborn入门详细教程 导读 今天我们来介绍 seaborn 这是一个基于matplotlib进行高级封装的可视化库 相比之下 绘制图表更为集成化 绘图风格具有更高的定制性 教程目录 01 初始seab
  • 一文带你了解序列化与反序列化基本原理与操作

    文章目录 一 什么是序列化与反序列化 二 为什么我们需要序列化与反序列化 三 步骤说明 四 注意说明 五 代码说明 六 序列化与反序列化原理 一 什么是序列化与反序列化 序列化是指将对象转换为字节序列的过程 以便于存储或传输 在序列化过程中
  • torch的交叉熵损失函数(cross_entropy)计算(含python代码)

    1 调用 首先 torch的交叉熵损失函数调用方式为 torch nn functional cross entropy input target weight None size average None ignore index 100
  • idea 配置log4j

    1 导入log4j jar包 放在lib目录下 右键jar包 createlibrary 好像是 还有一种方法 2 配置log4j properties文件 设置输出信息 DEBUG级别和ERROR级别 log4j rootLogger d
  • 【图像检测】基于计算机视觉和滤波实现热红外图像温度检测系统含Matlab源码

    1 简介 读入给定的热红外图像 完成彩色图像灰度化 对步骤1中的灰色图像加噪 并采用不同的方法对不同的加噪图像去噪 使得去噪效果更佳 基于灰度图像中标准温度对照带 提取图像中各像素的温度值 2 部分代码 function varargout
  • Keras-多输入多输出【多任务】

    模型结果设计 代码 from keras import Input Model from keras layers import Dense Concatenate import numpy as np from keras utils i
  • Reid 论文文章

    Reid 论文探索 Harmonious Attention Network for Person Re Identification 着重在空间注意力和通道注意力上 结构设计比较巧妙 通过global和local的两级将特征细化 参考链接
  • [QT_027]Qt学习之按钮类控件(QPushButton、QToolButton、QRadioButton、QCheckBox、QCommandLinkButton)

    本文转自 Qt编程指南 作者 奇先生 Qt编程指南 Qt新手教程 Qt Programming Guide 5 1 按钮类的控件 本节介绍图形程序里常见的按钮 包括普通的按钮按钮 之前都用过多次了 还有单选按钮 复选框 命令链接按钮等 本节
  • 【python】jupyter notebook的快捷键

    Jupyter Notebook 提供了许多快捷键以方便操作 这些快捷键分为两类 命令模式快捷键和编辑模式快捷键 命令模式快捷键 在命令模式下 细胞边框是蓝色的 这些快捷键主要用于操作单元格 Enter 进入编辑模式 Shift Enter
  • NoSQL数据库如何选型

    NoSQL数据库在体系结构和功能上各不相同 整体来说 NoSQL数据库选型需要从数据模型和存储特性两方面综合考虑 一 数据模型 NoSQL数据模型包括键值 宽列 图形 文档等 与关系型数据库使用的数据结构不同 具有高扩展性 弱一致性 适合分
  • Ozone作为Hadoop FileSystem的配置使用

    文章目录 前言 Ozone FileSystem的Hadoop兼容性文件系统实现原理 Ozone FileSystem的配置 Ozone FileSystem的使用 引用 前言 Ozone作为同样隶属于Hadoop大数据生态圈的一个系统 尽
  • Dynamics CRM 365 电子邮箱配置

    Dynamics CRM 365 如何配置电子邮箱服务 电子邮件服务器配置 电子邮件配置设置 邮箱 注意 天天在踩坑 坑坑不一样 启用CRM的电子邮箱通知 1 先配置电子邮件服务器文件 2 在电子邮件配置中设置好相关属性配置 3 在邮箱中设
  • NPOI从数据库中调取数据直接导出到EXCEL中

    https blog csdn net weixin 30725315 article details 98769248 一 关于NPOI NPOI是POI项目的 NET版本 是由 Tony Qu http tonyqus cnblogs
  • 证明正定矩阵的充要条件:全部顺序主子式大于0

    定理 f x T A x f x TAx f xTAx 正定的充要条件是