正定矩阵 Positive definite matrix

2024-02-11 21:20

本文主要是介绍正定矩阵 Positive definite matrix,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!


Positive definite matrix

  说道正定矩阵,这名字又能糊住很多人。心里默默念叨“我X,这是个什么矩阵,干嘛用的,名字这么叼”


  首先,本来可能很浅显的数学概念,从原来的英文描述翻译成中文之后会变得更加抽闲!

  我看了MIT的公开课视频才知道所谓的正定矩阵是 positive definite matrix。

下面是授课视频的link

   http://v.163.com/movie/2010/11/3/P/M6V0BQC4M_M6V2B5J3P.html


blog符号说明: V’ 表示矩阵V的转置

                           x^2 表示x的平方

-------------------------------------------------------------------------------------------------------------------------------------------






首先回归到矩阵的本质,矩阵是什么,就是多元函数的系数一种很巧妙的表示方法!无他!

对于这些系数的探讨和研究直接可以分析出函数的特性。


首先由下面的问题引出positive definite matrix的定义


对于一个系数矩阵

A = [a b 

       c  d];

列向量V = [x 

                  y];

V' A V = [ax + cy, bx +dy] * V

          = (ax+cy)*x + (bx+dy)*y

          = ax^2 + dy^2 + (c+b)x*y;

想想,对于任意的x和y变量,我们能够找到适当的a b c d值,是不是可以把V'AV的值确定为大于等于0的?!

比方A = [1 0

      0  1];

这样V' AV = x^2 + y^2; 任意的xy,该式子的值都大于等于0;

如果除开V = [0 ;0]这种情况,其余的地方都是大于0的。

这就是正定矩阵的来源!


正定矩阵的定义就是,对于V不等于零向量的任意向量,有矩阵A使得,V' AV > 0成立,我们把这样的矩阵A称作正定矩阵。


关于positive definite matrix的判定:

1. 矩阵特征值均大于0

2. 行列式的值大于0.


举例,下图中的待确定的矩阵,第二行第二列的元素未知,如果我们填入适当的数值,应该填写多大的数值可以使得该矩阵正定呢?

显然问题很简单,大于18即可。

如果是刚好等于18呢?我们把这种情况称为半正定——positive semi-definite matrix。



恩,既然是说函数V'AV的值在除开零向量出,任意的V均使得函数大于0,二元函数我们可以左图检验。


比方说这里矩阵A = [ 2 6

                                  6 7];

很明显,这里的A是非正定的,对于变量V = [x1 x2],


观察该函数的图像,这个函数的图像像个马鞍,可以观察到,函数是有小于0的部分的。很好的印证了非正定矩阵。


我们只需稍作A的调整,把7改成大于18的数值(比方说20),即可得到正定矩阵,函数的图像除开0点之外,对于所有V的取值,都大于0.

下面是函数图像,Z变量即函数值。我们可以直观的看出函数最小不小于0。



到现在应该明白了,正定矩阵干嘛用的?使得V' AV大于0。进一步的,我们可以利用函数求导的性质探究函数极值的情况。





这篇关于正定矩阵 Positive definite matrix的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/700838

相关文章

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

[论文笔记]LLM.int8(): 8-bit Matrix Multiplication for Transformers at Scale

引言 今天带来第一篇量化论文LLM.int8(): 8-bit Matrix Multiplication for Transformers at Scale笔记。 为了简单,下文中以翻译的口吻记录,比如替换"作者"为"我们"。 大语言模型已被广泛采用,但推理时需要大量的GPU内存。我们开发了一种Int8矩阵乘法的过程,用于Transformer中的前馈和注意力投影层,这可以将推理所需

线性代数|机器学习-P35距离矩阵和普鲁克问题

文章目录 1. 距离矩阵2. 正交普鲁克问题3. 实例说明 1. 距离矩阵 假设有三个点 x 1 , x 2 , x 3 x_1,x_2,x_3 x1​,x2​,x3​,三个点距离如下: ∣ ∣ x 1 − x 2 ∣ ∣ 2 = 1 , ∣ ∣ x 2 − x 3 ∣ ∣ 2 = 1 , ∣ ∣ x 1 − x 3 ∣ ∣ 2 = 6 \begin{equation} ||x

【线性代数】正定矩阵,二次型函数

本文主要介绍正定矩阵,二次型函数,及其相关的解析证明过程和各个过程的可视化几何解释(深蓝色字体)。 非常喜欢清华大学张颢老师说过的一段话:如果你不能用可视化的方式看到事情的结果,那么你就很难对这个事情有认知,认知就是直觉,解析的东西可以让你理解,但未必能让你形成直觉,因为他太反直觉了。 正定矩阵 定义 给定一个大小为 n×n 的实对称矩阵 A ,若对于任意长度为 n 的非零向量 ,有 恒成

python科学计算:NumPy 线性代数与矩阵操作

1 NumPy 中的矩阵与数组 在 NumPy 中,矩阵实际上是一种特殊的二维数组,因此几乎所有数组的操作都可以应用到矩阵上。不过,矩阵运算与一般的数组运算存在一定的区别,尤其是在点积、乘法等操作中。 1.1 创建矩阵 矩阵可以通过 NumPy 的 array() 函数创建。矩阵的形状可以通过 shape 属性来访问。 import numpy as np# 创建一个 2x3 矩阵mat

【UVA】10003-Cutting Sticks(动态规划、矩阵链乘)

一道动态规划题,不过似乎可以用回溯水过去,回溯的话效率很烂的。 13988658 10003 Cutting Sticks Accepted C++ 1.882 2014-08-04 09:26:49 AC代码: #include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include

算法练习题17——leetcode54螺旋矩阵

题目描述 给你一个 m 行 n 列的矩阵 matrix ,请按照 顺时针螺旋顺序 ,返回矩阵中的所有元素。  代码 import java.util.*;class Solution {public List<Integer> spiralOrder(int[][] matrix) {// 用于存储螺旋顺序遍历的结果List<Integer> result = new ArrayList

线性代数 第六讲 特征值和特征向量_相似对角化_实对称矩阵_重点题型总结详细解析

文章目录 1.特征值和特征向量1.1 特征值和特征向量的定义1.2 特征值和特征向量的求法1.3 特征值特征向量的主要结论 2.相似2.1 相似的定义2.2 相似的性质2.3 相似的结论 3.相似对角化4.实对称矩阵4.1 实对称矩阵的基本性质4.2 施密特正交化 5.重难点题型总结5.1 判断矩阵能否相似对角化5.2 已知两个矩阵相似,求某个矩阵中的未知参数5.3 相似时,求可逆矩阵P,使

最大子矩阵和问题归纳总结

一,最大子矩阵问题: 给定一个n*n(0< n <=100)的矩阵,请找到此矩阵的一个子矩阵,并且此子矩阵的各个元素的和最大,输出这个最大的值。 Example: 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2 其中左上角的子矩阵: 9 2 -4 1 -1 8 此子矩阵的值为9+2+(-4)+1+(-1)+8=15。 二,分析 子矩阵是在矩阵