27正定矩阵

news2024/10/5 19:15:12

这一节进入正定矩阵的内容,什么叫做正定矩阵?为什么我们对矩阵正定这么感兴趣?

PS:这一节将前面所有的概念都融合在一起:主元、行列式、特征值、不稳定性

一、正定矩阵的判断方法

为了说明问题,我们先考虑二阶矩阵 A A A
A = [ a b b c ] A=\begin{bmatrix}a&b\\b&c\\\end{bmatrix} A=[abbc]

  • 特征值判断法 λ 1 > 0 \lambda_1>0 λ1>0 λ 2 > 0 \lambda_2>0 λ2>0
  • 行列式判断法 a > 0 a>0 a>0 a c − b 2 > 0 ac-b^2>0 acb2>0
  • 主元判断法 a > 0 a>0 a>0 a c − b 2 a > 0 \frac{ac-b^2}{a}>0 aacb2>0
  • 判断式 x T A x > 0 x^TAx>0 xTAx>0

PS:线性代数范围内,正定矩阵需要是对称矩阵。

例子:
[ 2 6 6 ? ] \begin{bmatrix}2&6\\6&?\end{bmatrix} [266?]
填入什么数字会使得其为正定矩阵?

  • 填入数字18,该矩阵行列式恰好为0,矩阵称为半正定;
  • 填入大于18的数字,矩阵行列式顺序主子式均为正数

OK!为什么正定对于我们线性代数非常重要?因为它与工程联系非常密切,切入点就是最后一个判断正定的条件 x T A x > 0 x^TAx>0 xTAx>0

对于半正定矩阵 A = [ 2 6 6 18 ] A=\begin{bmatrix}2&6\\6&18\end{bmatrix} A=[26618],我们写出其计算式 x T A x x^TAx xTAx
x T A x = [ x 1 x 2 ] [ 2 6 6 18 ] [ x 1 x 2 ] = 2 x 1 2 + 12 x 1 x 2 + 18 x 2 2 x^TAx=\begin{bmatrix}x_1&x_2\end{bmatrix}\begin{bmatrix}2&6\\6&18\end{bmatrix}\begin{bmatrix}x_1\\x_2\end{bmatrix}=2x_1^2+12x_1x_2+18x_2^2 xTAx=[x1x2][26618][x1x2]=2x12+12x1x2+18x22
对上述的式子进行配方:
z = a x 1 2 + 2 b x 1 x 2 + c x 2 2 = 2 x 1 2 + 12 x 1 x 2 + 18 x 2 2 = 2 ( x 1 + 3 x 2 ) 2 \begin{aligned} z&=ax_1^2+2bx_1x_2+cx_2^2\\ &=2x_1^2+12x_1x_2+18x_2^2=2(x_1+3x_2)^2 \end{aligned} z=ax12+2bx1x2+cx22=2x12+12x1x2+18x22=2(x1+3x2)2
我们知道 A x Ax Ax 是线性的,但是左乘了一个 x T x^T xT就变成了二阶的,这种形式称为二次型形式 (Quadratic From),它是“纯” 二次形式的,没有线性部分,常数项和3、4或者其他次方。给出正定矩阵的定义,对于任何 x 1 , x 2 x_1,x_2 x1,x2所有值都大于零。

无论 x 1 , x 2 x_1,x_2 x1,x2如何取值,对应的结果均大于零。如何研究这个曲线的形状,一种可行的方法就是截面法,通过固定其中的一个值,然后研究其形状,最后在将固定的值进行遍历就能大概知道这个曲线的形状了。

通过配方,固定 x 1 x_1 x1 可以看出 z = 2 ( x 0 + 3 x 2 ) 2 z=2(x_0+3x_2)^2 z=2(x0+3x2)2,它是一个抛物线且顶点坐标为:
x 2 = − 1 3 x 0 x_2=-\frac{1}{3}x_0 x2=31x0
从与 x x x 垂直的截面来看,它首先是一个抛物线,随着 x 1 x_1 x1 的增加,对称轴坐标 x 2 x_2 x2逐渐减少。从 X X X轴上截面看就好像一个不断平移的抛物线。同理,如果固定 x 2 x_2 x2 研究这个界面,也是一个不断平移的抛物面。
在这里插入图片描述

二、正定矩阵与空间曲线草图

考察一个对称矩阵:
A = [ 2 6 6 ? ] A=\begin{bmatrix}2&6\\6&?\end{bmatrix} A=[266?]
在问号中填入不同的数字,矩阵可能满足正定,也可能是非正定的,甚至有可能是半正定。现在的问题是,正定性与其对应的二次型曲线有什么关系?满足正定是怎么样得到?不满足的又是何种情况?

OK!假如问号中的数字是20,对称矩阵变成:
A = [ 2 6 6 20 ] A=\begin{bmatrix}2&6\\6&20\end{bmatrix} A=[26620]
对于二阶矩阵最简单的方式就是通过从左上到左下行列式,显然分别是2和4都是大于零的,所以矩阵 A A A 是正定的。因为我们需要研究的是对应二次型的空间曲线,所以先将其写成二次型:
x T A x = [ x 1 x 2 ] [ 2 6 6 20 ] [ x 1 x 2 ] = 2 x 1 2 + 12 x 1 x 2 + 20 x 2 2 x^TAx=\begin{bmatrix}x_1&x_2\end{bmatrix}\begin{bmatrix}2&6\\6&20\end{bmatrix}\begin{bmatrix}x_1\\x_2\end{bmatrix}=2x_1^2+12x_1x_2+20x_2^2 xTAx=[x1x2][26620][x1x2]=2x12+12x1x2+20x22
为了在直角坐标系中表示,将方程写成我们熟悉的形式:
z = 2 x 2 + 12 x y + 20 y 2 z=2x^2+12xy+20y^2 z=2x2+12xy+20y2
不妨从特殊点着手,首先方程过 ( 0 , 0 , 0 ) (0,0,0) (0,0,0) 点、 ( 1 , 0 , 2 ) (1,0,2) (1,0,2)
在这里插入图片描述
用一个平行于 X O Y XOY XOY 且过原点平面去截取这个曲面:
z = 2 x 2 + 12 x y + 20 y 2 0 = x 0 x \begin{aligned} z&=2x^2+12xy+20y^2\\ 0&=x_0x \end{aligned} z0=2x2+12xy+20y2=x0x
联立上述式子,有:
z = 20 y 2 z=20y^2 z=20y2
在这里插入图片描述
再用 x − y = 0 x-y=0 xy=0 平面去截取这个曲面:
z = 2 x 2 + 12 x y + 20 y 2 x − y = 0 \begin{aligned} z&=2x^2+12xy+20y^2\\ x-y&=0 \end{aligned} zxy=2x2+12xy+20y2=0
联立得:
z = 34 y 2 z=34y^2 z=34y2
在这里插入图片描述
不过,现代软件可以很轻松帮你绘制这样的图形:
**加粗样式**
所谓正定,就是这个曲线所有的高度值都是大于零的。

再来看一个非正定的矩阵对应的二次曲线:
[ 2 6 6 2 ] \begin{bmatrix}2&6\\6&2\end{bmatrix} [2662]

在这里插入图片描述
有部分的值是负数,所以不是正定的。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/363452.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

低代码平台和无代码平台哪个更适合开发企业管理系统?

编者按:本文分析了开发企业管理系统所需要的平台特性,并根据这些特点和低代码无代码的优劣比较,得出低代码平台更适合开发企业管理系统。关键词:私有化部署,可视化设计,源码交付,数据集成&#…

研究人员在野外发现大量的信息窃取者 “Stealc “的样本

一个名为Stealc的新信息窃取者正在暗网上做广告,它可能成为其他同类恶意软件的一个值得竞争的对象。 "SEKOIA在周一的一份报告中说:"威胁行为者将Stealc作为一个功能齐全、随时可用的窃取者,其开发依赖于Vidar、Raccoon、Mars和Re…

LeetCode 73. 矩阵置零

LeetCode 73. 矩阵置零 难度:middle\color{orange}{middle}middle 题目描述 给定一个 KaTeX parse error: Double subscript at position 3: _m_̲ x _n_ 的矩阵,如果一个元素为 0 ,则将其所在行和列的所有元素都设为 0 。请使用 原地 算法…

Power Query 再谈日期表

Power Query 再谈日期表 需要附件的朋友可以点赞关注后,在公众号(焦棚子)后台回复 日期表 三个字即可领取。 一、背景 在 Power BI 或者 Power Pivot 做业务分析中,必不可少的就是日期表。我们继续优化调制出符合国人习惯日期表。 本次更新&#xff…

harbor 仓库迁移升级

harbor 仓库迁移升级 harbor仓库安装数据传输仓库切换版本 v1.8.0 v2.3.5 harbor仓库安装 环境准备:安装docker详见:docker 的介绍和部署,并下载docker-compose详见:docker 三剑客compose。 现有支持的安装harbor仓库的方式有两…

犹太教、基督教、伊斯兰教的区别与联系

一、犹太教、基督教、伊斯兰教的简明关系图二、犹太教、基督教、伊斯兰教的主要区别注:弥赛亚(希伯莱语)就是基督(希腊语),意思是“救世主”。注:伊斯兰教的观点是:穆罕默德不是伊斯…

华为OD机试 - 优雅数组(Python)【2023-Q1 新题】

华为OD机试300题大纲 参加华为od机试,一定要注意不要完全背诵代码,需要理解之后模仿写出,通过率才会高。 华为 OD 清单查看地址:blog.csdn.net/hihell/category_12199275.html 华为OD详细说明:https://dream.blog.csdn.net/article/details/128980730 优雅数组 | 华为…

S5P6818_系统篇(2)源码编译及烧录

源码获取 源码获取和操作流程 1.下载liunux下的系统制作脚本,可以烧录系统和构建镜像 git clone https://github.com/friendlyarm/sd-fuse_s5p6818.git 如果出现git错误可使用如下方法: git config --global http.sslverify false 2.阅读该工具rea…

Linux文件系统操作与磁盘管理

查看磁盘和目录的容量 使用 df 命令查看磁盘的容量 df在实验楼的环境中你将看到如下的输出内容: 但在实际的物理主机上会更像这样: 物理主机上的 /dev/sda2 是对应着主机硬盘的分区,后面的数字表示分区号,数字前面的字母 a 表示…

Magisk Yyds.Msu定制方案 原理 ROOT隐藏 风控绕过

介绍 与市面上所有方案都不一样,基于最新Magisk代码分支,通过修改Magisk的挂载方式,解决了root隐藏的问题,可绕过如momo独立进程,挂载检测等各种风控,并因为个人项目原因,默认把Zygisk禁用&…

【手把手一起学习】(二) Altium Designer 20 工程创建

1 Altium Designer 20 工程创建 手把手一起学习,从零开始制作一款属于自己的STM32核心板,规范AD使用,养成良好习惯 1.1 新建工程 首先,创建工程,操作如图所示: 选择默认工程,输入工程名字和路…

吊炸天的1500亿,到底有多魔性

超1500亿的数据安全规模 数据安全是近几年来最热门的行业话题,数字化转型使得数据作为新型生产要素的重要性日益凸显,数据经济已成为经济发展的核心驱动力,数字经济的蓬勃催生了巨大的数据安全需求。 临近春节前夕,在工信部、国…

Java API 文档搜索引擎

1. 认识搜索引擎:在搜狗搜索的搜索结果页中, 包含了若干条结果, 每一个结果包含了图标, 标题, 描述, 展示URL等搜索引擎的本质:输入一个查询词, 得到若干个搜索结果, 每个搜索结果包含了标题, 描述, 展示URL和点击URL2. 搜索引擎思路:2.1 搜索的核心思路:当前我们有很多的网页(…

日志收集笔记(架构设计、Log4j2项目初始化、Lombok)

1 架构设计 ELK 技术栈架构设计图: 从左往右看, Beats:主要是使用 Filebeat,用于收集日志,将收集后的日志数据发送给 Kafka,充当 Kafka 的生产者Kafka:高性能消息队列,主要起缓冲…

GitHub标星30K+的Java面试八股文长啥样?

2023年的互联网行业竞争越来越严峻,面试也是越来越难,一直以来我都想整理一套完美的面试宝典,奈何难抽出时间,这套1000道的Java面试手册我整理了整整1个月,上传到Git上目前star数达到了30K 一、32 道 MySQL 面试题 1&…

论文阅读:RA-Depth: Resolution Adaptive Self-Supervised Monocular Depth Estimation

中文标题: 分辨自适应自监督单目深度估计 提出问题 传统的无监督深度估计方法如果在固定分辨率上训练,则迁移到其他分辨率上存在严重退化。 创新点 提出一种分辨率自适应的无监督深度学习框架(RA-Depth)来学习场景深度尺度不变性。具体地…

postman入门到精通之【接口知识准备】(一)

postman入门到精通之【接口知识准备】(一) 目录:导读 前言 接口测试概念 接口测试 接口测试的原理 常用接口测试工具 接口测试基础知识 接口的定义 接口的分类 HTTP接口 Web Service接口 RESTful接口 HTTP请求 统一资源定位符&…

使用终端工具给你的电脑发送弹窗提醒

大家好,我是良许。 现在人手一部智能手机,这些智能手机都有个非常实用的功能,那就是弹窗提醒。当我们收到短信,或者微信信息时,手机就会弹窗显示信息的大致内容。有了这个功能你就不会错过重要信息了。 电脑上也有类…

Android 流量统计

Android 流量统计最近项目上有一个应用流量统计的功能需要实现,在此总结一下 流量统计架构 在Android9.0之前,流量监控是基于xt_qtaguid模块的,通过读取/proc/net/xt_qtaguid/stats文件内容进行解析获取对应流量数据。 Android9.0之后&…

51单片机学习笔记_12 LCD1602 原理及其模块化代码

LCD1602 liquid crystal display 液晶显示屏,一种字符型液晶显示模块,可以显示 16*2 个字符,每个字符是 5*7 点阵。 P0 P2 会和数码管、LED 一定程度上冲突。 地。 Vcc。 调对比度的。 RS:数据指令端。1代表 DB 是数据&#x…