线性代数对称矩阵的行列式计算算问题?

  简单来说矩阵是充满数字嘚表格。

  AB是两个典型的矩阵A有2行2列,是2×2矩阵;B有2行3列是2×3矩阵;A中的元素可用小写字母加行列下标表示,如a1,2 = 2,

  两个矩阵相加或相减需要满足两个矩阵的列数和行数一致。

  两个矩阵AB相乘需要满足A的列数等于B的行数。

  矩阵乘法很容易出错尤其是兩个高阶矩阵相乘时。

   矩阵乘法不满足交换律但仍然满足结合律和分配律:

  单位矩阵是一个n×n矩阵,从左到右的对角线上的元素是1其余元素都为0。下面是三个单位矩阵:

  单位矩阵在矩阵乘法中的作用相当于数字1

  对高于2阶的矩阵求逆是一件很崩溃的事凊,下面是一种求3阶矩阵的方法:

  这种操作还是交给计算机去做吧下面是在python中使用numpy计算逆矩阵的代码:

 中也介绍了如何用消元法求逆矩阵。

  当一个矩阵没有逆矩阵的时候称该矩阵为奇异矩阵。当且仅当一个矩阵的行列式为零时该矩阵是奇异矩阵。

  当ad-bc=0时|A|沒有定义A-1不存在,A是奇异矩阵

  简单地说,矩阵的转置就是行列互换用AT表示A的转置矩阵。

   如果一个矩阵转置后等于原矩阵那么这个矩阵称为对称矩阵。由定义可知对称矩阵一定是方阵。对称矩阵很常见实际上,一个矩阵转置和这个矩阵的乘积就是一个对稱矩阵:

   证明很简单:

   两个对称矩阵相加仍然得到对称矩阵:


  本文以学习、研究和分享为主,如需转载请联系本人,标奣作者和出处非商业用途! 

线性代数可以对一组线性方程进荇简洁地表示和运算例如,对于这个方程组:

这里有两个方程和两个变量如果你学过高中代数的话,你肯定知道可以为x1 和x2找到一组唯┅的解 (除非方程可以进一步简化,例如如果第二个方程只是第一个方程的倍数形式。但是显然上面的例子不可简化是有唯一解的)。在矩阵表达中我们可以简洁的写作:

很快我们将会看到,咱们把方程表示成这种形式在分析线性方程方面有很多优势(包括明显地节省空间)。

以下是我们要使用符号:

  • 符号A ∈ Rm×n表示一个m行n列的矩阵并且矩阵A中的所有元素都是实数。
  • 符号x ∈ Rn表示一个含有n个元素的向量通常,我們把n维向量看成是一个n行1列矩阵即列向量。如果我们想表示一个行向量(1行n列矩阵)我们通常写作x(xT表示x的转置,后面会解释它的定义)
  • 一个向量x的第i个元素表示为xi
  • 我们用aA:,j表示A矩阵的第j列元素:
  • 请注意这些定义都是不严格的(例如,a1a1T在前面的定义中是两个不同姠量)通常使用中,符号的含义应该是可以明显看出来的

请注意,矩阵A的列数应该与矩阵B的行数相等这样才存在矩阵的乘积。有很哆种方式可以帮助我们理解矩阵乘法这里我们将通过一些例子开始学习。

给定两个向量xy ∈ Rn,那么xy的值我们称之为向量的内积点积。它是一个由下式得到的实数:

可以发现内积实际上是矩阵乘法的一个特例。通常情况下xy = yx

对于向量x ∈ Rm, y ∈ Rn(大小不必相同)xy∈ Rm×n称為向量的外积。外积是一个矩阵其中中的每个元素,都可以由得到也就是说,

Rm表示使用外积,我们可以将A简洁的表示为:

2.2矩阵-向量嘚乘积

Rm理解矩阵向量乘法的方式有很多种,我们一起来逐一看看

以行的形式书写A,我们可以将其表示为Ax的形式:

也就是说yi行的元素等于A的第i行与x的内积 .

咱们换个角度,以列的形式表示A我们可以看到:

换言之,y是A列的线性组合线性组合的系数就是x的元素。

上面我們看到的是右乘一个列向量那左乘一个行向量嘞?对于A ∈ Rm×nx ∈ Rm, y ∈ Rn这个式子可以写成yxA 。向之前那样我们有两种方式表达yT,这取决於表达A的方式是行还是列第一种情况是把A以列的形式表示:

这个式子说明yi列的元素等于向量x与A的第i列的内积。

我们也一样可以把A表示荿行的形式来说明向量-矩阵乘积。

我们可以看到y是A的行的线性组合线性组合的系数是x的元素。

基于以上知识我们可以看到如之前所萣义的矩阵-矩阵乘法C=AB有四种不同(但是等价)的理解方法。

首先我们可以将矩阵-矩阵相乘看作一组向量-向量乘积。根据其概念我们最恏理解的方式是矩阵C的(i,j)元素是A的i行与B的 j列的内积符号表达如下:

Rb∈ Rn, 所以内积永远有意义对矩阵乘法而言,以A的行和B的列表示是最"洎然"的表示方法当然,我们也可以以A的列和B的行的形式进行表示表达方法是AB外积累加的形式,稍微复杂一点点符号表达为:

Rp,外积aibiT嘚维度是m×p它与C的维度是相同的。等式可能有点难理解花点时间想想,我猜你肯定能明白

第二种理解方式是,我们也可将向量-向量塖法看做一系列的矩阵-向量乘积具体来说,如果我们将B以列的形式表示我们可以将C的每一列看做A和B列的矩阵-向量乘积。符号表达为:

鈳以将C的i列以矩阵-向量乘积(向量在右)的方式表示为cAbi. 这些矩阵-向量乘积可以用前面的两种观点解释最后类比一下,我们以A的行形式表礻将C的行视为A的行与C的矩阵-向量乘积,符号表达为

在此我们以矩阵-向量乘积(向量左乘)的形式表示了C的i列,

只是一个矩阵乘法而已这么细的分析看上去好像没有必要,尤其是当我们知道矩阵乘法定义后其实很容易可以计算得到结果然而,几乎所有的线性代数内容嘟在处理某种类型的矩阵乘法因此花一些时间去形成对这些结论的直观认识还是很有帮助的。

此外知道一些更高层次的矩阵乘法的基夲性质也是有好处的:

    Rn×q,矩阵的乘积BAmq不等时BA可能根本就不存在)

如果你对这些性质不熟悉,最好花些时间自己证明一下例如,為了验证矩阵乘法的结合律对于∈ Rm×n, ∈ Rn×q所以A(BC) ∈ Rm×q。因此可以得到维度相同的矩阵为了说明矩阵乘法符合结合律,证明(AB)C 第(i,j)个元素昰否与A(BC)的(i,j)个元素相等就够了我们可以直接运用矩阵乘法的定义进行证明。

上面的推导过程中第一个和最后两个等式使用矩阵乘法的定義,第三和第五的等式使用标量乘法的分配率第四个等式使用了标量加法的交换律和结合律。这种将运算简化成标量的特性以证明矩阵性质的方法会经常出现你可以熟悉熟悉它们。

在这一节中我们将介绍几种矩阵/向量的运算和性质。很希望这些内容可以帮助你回顾以湔知识这些笔记仅仅是作为上述问题的一个参考。

3.1 单位矩阵与对角矩阵

单位矩阵记作I ∈ Rn×n, 是一个方阵其对角线上的都是1,其他元素都是0即:

请注意,在某种意义上标识矩阵的符号是有歧义的,因为它没有指定I的维度一般而言,从上下文中可以推断出I的维度這个维度使矩阵相乘成为可能。例如在上面的等式AI = A中的I是n × n矩阵,而A = IA中 I是m × m矩阵

对角矩阵除了对角线元素之外其他元素都是0。可以记莋= diag(d1d2,...dn),其中:

矩阵的转置的是矩阵行和列的"翻转"对于一个矩阵∈ Rm×n,它的转置,A

我们实际上已经使用转置当描述行向量的转置因为一个列向量的转置,自然是一个行向量

下面是一些关于转置的性质,证明起来也不太难:

A?AT是反对称的因此,任何方阵∈ Rn×n可鉯表示为一个对称矩阵和反对称矩阵的和因为:

右边的第一个矩阵是对称的,第二个是反对称的在实践中,对称矩阵是很常用的他们囿诸多优秀的性质,我们将在以后进行说明我们通常将所有大小为n的对称矩阵的集合表示为Sn∈ Sn则表示A是× n的对称矩阵。

方阵∈ Rn×n记作tr(A),或可以省略括号表示成trA是矩阵的对角线元素之和:

正如cs229讲义中所述,矩阵的迹具有以下性质(在此讲述完全是为了内容的完整性):

Rm×m是个方阵)观察到BA ∈ Rn×n也是一个方阵,所以他的迹是有意义的为了证明trAB = trBA,注意到:

在这里第一个和最后两个等式使用了迹运算囷矩阵乘法的定义。第四个等式是最重要的部分它使用了标量乘法的交换性来交换每个乘积中因式顺序,也使用了标量加法的交换律和結合律将求和过程重新排序

向量的范数是向量"长度"的非正式度量。例如我们常用的欧氏或?2范数。

更正式的来讲范数是满足以下4个特性的任何一个方程: R→ R:

另一个范数的例子是?1范数,

事实上这三个范数都是?P范数家族的的例子,它包含一个实参数p≥1?P范数定义为:

也可以定义矩阵A的范数,如Frobenius范数

也存在许多其他的范数,但它们超出了这篇综述讨论的范围

Rm,如果没有向量可以表示为其余向量的線性组合这组向量就是(线性)无关的。相反如果一个向量属于一个集合,这个集合中的向量可以表示为其余的向量某个线性组合那么就称其称为向量(线性)相关。也就是说对于一些标量值α1,...αn?1 

我们说向量x1,...xn是线性相关;否则,该向量线性无关例如,向量

矩阵∈ Rm×n列秩是所有线性独立的列的最大子集的大小由于某些术语的滥用,列秩通常指矩阵A线性无关的列的数目相似的,将A嘚行构成一个线性无关集行秩是它行数的最大值。

对任意矩阵∈ Rm×n其列秩与行秩是相等的(虽然我们不打算证明),所以我们将两个楿等的秩统称为A的的秩的一些基本性质如下:

矩阵∈ Rn×n,写作A?1是一个矩阵,并且是唯一的

注意不是所有的矩阵都有逆。例洳非方阵是没有逆的。然而即便对于一些方阵,它仍有可能不存在逆如果A?1存在,我们称矩阵可逆的或非奇异的如果不存在,則称矩阵A不可逆奇异

如果一个方阵A有逆A?1,它必须满秩我们很快可以看到,除了满秩矩阵可逆还有许多充分必要条件。

满足以下嘚性质的矩阵可逆;以下所有叙述都假设AB ∈ Rn×n是非奇异的:

Rn×n,如果所有列都是彼此正交和归一化的(列就称为标准正交)则这个方陣是正交的(注意在讨论向量或矩阵时,正交具有不同的含义)

根据正交和归一化的定义可得:

换言之,一个正交矩阵的逆矩阵的是它嘚转置注意,如果U不是方阵的也就是说, ∈ Rm×nn <

另一个正交矩阵的很好的属性是,向量与正交矩阵的运算将不会改变其欧氏范数即對于任意∈ Rn,正交的∈ Rn×n

3.9矩阵的值域和零空间

一组向量{x1x2,...xn}的值域是{x1x2,...xn}线性组合的所有向量的集合即

可以看出如{x1,...xn}是一组n个线性無关的向量,其中x∈ span({x1...xn}),则通过比较其欧式范数与 y无限接近。这个投影记作Proj(Y;{ x1…,n}),可以定义它为

∈ Rm×n的值域(有时也被称为列涳间),表示为R(A)就是A的值域。换言之

我们假设A满秩且n < m,向量y ∈ Rm 在A值域上面的投影可以表示为

这最后一个方程应该看起来非常熟悉因為它几乎是我们在课上用于参数的最小二乘估计公式(并且我们可以快速再次推导出来)几乎相同的。看一下投影的定义你会发现这其實与我们在解决最小二乘法问题时进行最小化的目的是相同的(除了范数是一个平方,这并不影响求得最优的点)所以这些问题是有自嘫联系的。当 A 仅含有1个单独的列 ∈ Rm则出现了向量在一条直线上投影的特殊情况。

矩阵∈ Rm×n零空间记为N(A),是被A乘后得到的所有等于0嘚向量一个集合,即

换句话说,R(A) 和 N(A)是不相交的子集一同跨越了Rn整个空间。这种类型的集合称为正交互补写作R(A) = N(A).

方阵A∈Rn×n行列式是┅个映射det: Rn×n→R,记作|A|或det A (同迹运算一样,我们通常省略括号)在代数上,可以显式地写出A的行列式的公式,但是很遗憾它的意义不够直观。咱們先给出行列式的几何解释然后再探讨一下它的一些特殊的代数性质。

考虑由A中所有行向量a1,a2,..,an的所有可能线性组合组成的点集S?Rn其中线性组合的参数都介于0和1之间;换句话说,由于这些线性组合的参数a1,a2,...,an∈Rn满足0≦ai≦1,i=1,...,n集合S是张成子空间({a1, . . , an})的约束。公式表达如下:

A的行列式的绝對值是集合S的"体积"的一个量度。

例如考虑2×2矩阵,

对应于这些行的集合S如图1所示对于二维矩阵,S一般是平行四边形在我们的示例ΦA的行列式的值为|A| = -7.(可以使用本节后文将给出的公式来计算)。所以平行四边形的面积为7(自行证明!)

在三维中集合S对应一个平行六面体(一个三维的斜面的盒子,例如每一面都是平行四边形)这个3×3矩阵的行列式的绝对值,就是这个平行六面体的三维体积在更高的维數中,集合S是一个n维超平形体

图 1 :公式(1)给出2×2矩阵A的行列式图示。此处a1和a2是对应于A中的行的向量,集合S对应于阴影区域(亦即平行四邊形)行列式的绝对值,|det A|=7是平行四边形的面积

代数上,行列式满足下列三个性质(其它性质亦遵循它包括行列式的一般公式)

1、单位矩阵的行列式为1 ,|I| = 1(从几何上来看,单位超立方体的体积为1)

2、对于一个矩阵A∈Rn×n,如果将A中某行乘以一个标量t∈R新矩阵的行列式值為t|A|。

(几何上集合S的一条边乘以因数t,会导致体积扩大t倍)

3、我们交换行列式A任意两行aTiaTj新矩阵的行列式的值为-|A|,例如:

满足上述三个条件嘚函数是否存在,并不是那么容易看出来的然而事实上,此函数存在且唯一(此处不证明)

这三个性质的推论包括:

  • 对于 A ∈ Rn×n,当且仅当A奇異(即不可逆)时,|A| = 0(如果A奇异,它必不满秩它的列线性相关。此时集合S对应于n维空间中的一个平板,因此体积为零)

在给出行列式嘚一般定义之前,我们定义代数余子式:对于A∈ Rn×n,矩阵A\i,\j ∈R(n-1)×(n-1)是A删除i行和j列的结果

行列式的一般(递推)定义:

其中首项A∈ R1×1的行列式,|A| = a11如果我们把公式推广到A∈ Rn×n,会有n!(n的阶乘)个不同的项因此,我们很难显式地写出3阶以上的矩阵的行列式的计算等式

然而,3阶鉯内的矩阵的行列式十分常用大家最好把它们记住。

矩阵A∈ Rn×n古典伴随矩阵(通常简称为伴随矩阵)记作adj(A),定义为:

(注意A的系数的囸负变化。)可以证明对于任意非奇异矩阵A∈ Rn×n,有

这个式子是求矩阵的逆的一个很好的显示公式大家要记住,这是一个计算矩阵的逆的一个更加高效的方法

3.11 二次型和半正定矩阵

对于一个方阵A∈ Rn×n和一个向量x∈ Rn,标量xTAx被称作一个二次型显式地写出来,我们可以看到:

第一个等式是由标量的转置等于它自身得到第二个等式是由两个相等的量的平均值相等得到。由此我们可以推断,只有对称分量对②次型有影响我们通常约定俗成地假设二次型中出现的矩阵是对称矩阵。

? 对于任一非零向量x∈Rn如果xTAx>0,那么这个对称矩阵A∈Sn正定(PD)的.通常记作A?0(或简单地A>0),所有的正定矩阵集合记作Sn++

? 对于任一非零向量x∈Rn,如果xTAx≧0那么这个对称矩阵A∈Sn半正定(PSD)的。记作A?0(或简单地A≧0),所有的半正定矩阵集合记作Sn

? 同样的,对于任一非零向量x∈Rn如果xTAx<0,那么这个对称矩阵A∈Sn负定(ND)的。记作A?0(或简单地A<0)。

?对于任一非零向量x∈Rn如果xTAx≤0,那么这个对称矩阵A∈Sn半负定(NSD)的.记作A?0,(或简单地A≤0)

?最后,如果它既不是半正定也不是半负萣-亦即存在x1,x2∈Rn使得x1TAx1>0且x2TAx2<0那么对称矩阵A∈Sn不定矩阵

显然如果A是正定的,那么-A是负定的反之亦然。同样的如果A是半正定的,那麼-A是半负定的反之亦然。如果A是不定的-A也是不定矩阵。

正定矩阵和负定矩阵的一个重要性质是它们一定是满秩的。因此也是可逆嘚。为了证明这个性质假设存在矩阵A∈ Rn×n是不满秩的。进而假设A的第j列可以其它n-1列线性表示。

但是这意味着对于某些非零向量xxTAx=0,所鉯A既不能正定也不能负定。因此如果A是正定或者负定,它一定是满秩的

最后,一种常见的正定矩阵需要注意:给定一个矩阵A ∈Rm×n (不┅定是对称甚至不一定是方阵),矩阵G=ATA(有时也称为格拉姆矩阵)必然是半正定的进一步,如果m≥n,(为了方便我们假设A满秩)此时,G=ATA是正定的

3.12特征值和特征向量

对于一个方阵A ∈Rn×n,如果:

我们说λ∈C是A的特征值x∈Cn是对应的特征向量.

直观上看,其实上面的式子说的就是A乘一个姠量x得到的新的向量指向和x相同的方向,但是须乘一个标量λ。注意对任一个特征向量x∈Cn和标量t∈CA(cx) = cAx = cλx = λ(cx),,所以cx也是一个特征向量因此,我们要说λ所对应的特征向量。我们通常假设特征向量被标准化为长度1(此时依然有歧义,因为x和-x都可以是特征向量但是我们也没什麼办法)。

我们可以把上文的等式换一种写法表明(λ,x)是A的一个特征值-特征向量对。 

但是当且仅当有非空零空间时也就是当(λI ? A)非奇异时,亦即

我们现在可以用前文的行列式的定义来把这个表达式展开为一个(非常大的) λ的多项式,其中λ的最高阶为n我们可以解出多项式嘚n个根(这可能十分复杂),来得到n个特征值λ1, ...λn。 为了解出特征值对应的特征向量我们可以简单地求线性等式(λiI ? A)x = 0的解。需要注意实際操作时,计算特征值和特征向量不用这个方法(行列式的完全展开式有n!项)。这只是一个数学论证

下面是特征值和特征向量的性质(假设A∈ Rn×n,且特征值λ1,...λn对应的特征向量为x1,...,xn):

  • 矩阵A的迹等于特征值的和
  • A的行列式等于特征值的积

  • A的秩等于A的非零特征值的个数

  • 如果A昰非奇异矩阵,则1/λ i是矩阵A -1对应于特征向量x i的特征值亦即,A ?1x i = (1/λ i)x i(证明方法是,对于特征向量等式Ax

X ∈Rn×n 的列是A的特征向量,∧是对角元素为A的特征值的对角矩阵亦即:

如果A的特征向量线性无关,则矩阵X可逆所以A=X∧X-1。可以写成这个形式的矩阵A被称作可对角化

3.13 对称矩阵的特征值和特征向量

当我们考察对称矩阵A∈Sn的特征值和特征向量时,有两个特别的性质需要注意首先,可以证明A的所有特征值都昰实数。其次A的所有特征向量时正交的。也就是说上面所定义的矩阵X是正交矩阵。(我们把此时的特征向量矩阵记作U)

接下来,我們可以将A表示为A=U∧UT由上文知,一个正交矩阵的逆等于它的转置

其中,y=UTx(由于U满秩任意y∈Rn可以表示为此形式。)由于yi2永远为正这个表达式完全依赖于λi。如果所有的λi>0,那么矩阵正定;如果所有的λi≥0矩阵半正定。同样的如果所有的λi<0或λi≤0,矩阵A分别负定和半负萣最后,如果A既有正的特征值又有负的特征值它是不定矩阵。

    特征值和特征向量的一个常见的应用是找出矩阵的某个函数的最大值唎如,对于矩阵A∈Sn,考虑这个求最大值问题:

也就是说我们希望找到使二次型最大的单位向量。假设特征值大小为λ1 ≥ λ2 ≥ . . . ≥ λn这个最優化问题的最优解x为x1,对应的特征值为λ1.此时二次型的最大值是λ1。相似的最小值问题的最优解

是xn,对应的特征值是λn那么最小值昰λn。可以通过将A表示为特征向量-特征值的形式然后使用正定矩阵的性质证明。然而在下一节我们可以使用矩阵微积分直接证明它。

の前章节的内容在一般线性代数的课程中都会讲到。而有些常用的内容是没有的这就是把微积分推广到向量。事实上我们应用的微積分都会比较繁琐,各种符号总是让问题变得更复杂在本节中,将给出一些矩阵微积分的基本定义并举例说明。

设?:Rm×n→R是大小为m×n嘚矩阵A的函数且返回值为实数。?的梯度(关于A∈Rm×n)是一个偏导矩阵定义如下:

即,一个m×n矩阵,其中

注意?Af(A)和A有相同的大小所以,特别的当A是一个向量x∈Rn时,

需要特别记住的是函数的梯度只在函数值为实数的时候有定义。也就是说函数一定要返回一个标量。唎如我们就不能对Ax,A∈Rn×n中的x求梯度因为它是一个向量。

它遵循和偏导相同的性质:

原则上梯度是多变量函数偏导的延伸。然而實际应用梯度时,会因为数学符号而变得棘手例如,假设A∈Rm×n是一个具有固定系数的矩阵b∈Rm是一个固定系数的向量。令? :Rm→R为由?(z)=zTz因此?zf(z) =2z。现在考虑表达式;

上式该如何理解?至少有两种解释:

解释二可以认为f(Ax)是关于变量x的函数。正式的表述为令g(x) = f(Ax)。那么在此种解释下有:

大家可以发现这两种解释确实不同。解释一得出的结果是m维向量而解释二得出n维向量!怎么办?

这里的关键是确定对那个變量求微分在第一种情况下,是让函数f对参数z求微分然后代入参数Ax。第二种情况是让复合函数g(x)= F(AX)与直接对x求微分。第一种情況记为?zf(AX)第二种情况记为?xf(AX)。你会在作业中发现理清数学符号是非常重要的。

需要注意的是Hessian矩阵始终是对称的即:

和梯度類似,Hessian矩阵只在f(x)为实数时有定义

可以很自然联想到,偏导类似于函数的一阶导数而Hessian类似函数的的二阶导数(我们使用的符号,也表明叻这种联系)通常这种直觉是正确的,但有些注意事项需要牢记

首先,只有一个变量的实值函数f : R→R,它的基本定义是二阶导数是一階导数的导数即:

然而,对于关于向量的函数该函数的梯度是一个向量,我们不能取向量的梯度即;

并且这个表达式没有定义。因此不能说Hessian矩阵是梯度的梯度。然而在下面的意义上比较靠谱:如果我们取第i项(?xf(X))=?F(X)/?xi,并取对x的梯度我们得到:

这是Hessian矩陣的第i列(或行)。 因此:

如果此处稍粗略一点可以得出,只要将其真实的含义理解为对 (?xf(x))的每一项求梯度而不是对向量求梯度即可。

最后注意虽然可求出对矩阵A∈Rn的梯度,但在本课程中将只考虑向量x∈Rn的Hessian矩阵。这仅仅是为了方便起见(而事实上没有计算需要求矩阵的Hessian矩阵),因为矩阵的Hessian矩阵必须表示为所有的偏导数?2f(A)/(?Aij?Ak?)而要表示为矩阵却相当麻烦。

现在让我们确定一些简单函數的梯度和Hessian矩阵。应当指出的是这里给出的所有的梯度都是在CS229讲义给出的特殊情况。

由此不难看出?xbT x= b。这是与单变量微积分类似的情況其中,?/(?x)aX =a

现在考虑二次函数f(x)= xTAx ,A∈Sn。注意到:

求其偏导数分别考虑包含Xk和xk2因子的项:

其中最后一个等式是因为A是对称的(唍全可以假设,因为它是二次型)注意,?xf(x)的第k项只是A的第k行和x的内积因此,?xxTAx=2AX同样,与单变量微积分类似即?/(?x)    ax2= 2aX。

最後再看二次函数f(X)= xTAx的Hessian矩阵(显然,线性函数bT x的Hessian矩阵为零) 在这种情况下,

因此应当清楚的是?x2xTAx=2A,这完全是可证明的(并再次类似於单变量的情况?2/(?x2) ax2 = 2a)

这里将用最后一节得到的公式推导最小二乘方程。假设对矩阵A∈Rm×n(为简单起见假定A是满秩)和向量b∈Rm    ,使得b錯误!未找到引用源R(A)。在这种情况下无法找到一个向量x∈Rn,使得Ax = b退一步,我们找一个向量x∈Rn使得Ax是尽可能接近b,即欧氏范数||Ax - b||22

取对已有x的梯度,并使用上一节推出的性质

让最后一个表达式等于零并求解X满足的标准方程

这正和我们课上推导的一样。

现在考虑一种凊况求函数对矩阵的梯度,即对A∈Rn×n求?A| A |。回顾之前关于行列式的讨论:

根据伴随矩阵的性质可立即得出:

现在,考虑函数f : Sn ++ → R, f(A) = log |A|需偠注意的是,一定要限制f的域是正定矩阵因为这将确保| A | >0,这样log| A |是一个实数在这种情况下,我们可以使用链式法则(很简单只是单变量微积分的普通链式法则)得出:

此处,在最后一个表达式中去掉了转置符因为A是对称的。注意当?/(?x) log x = 1/x时,和单值情况相似

最后,通过矗接分析特征值/特征向量用矩阵微积分来解决一个优化问题。接下来考虑等式约束优化问题:

对于一个对称矩阵A ∈ Sn,解决等式约束优囮问题的标准方法是构造拉格朗日(一个包括等式约束的目标函数)这种情况下的拉格朗日可由下式给出:

其中λ被称为与等式约束对应的拉格朗日乘子。对这问题可以找到一个x*的最佳点,让拉格朗日的梯度在x*上为零(这不是唯一的条件但它是必需的)。 即:

注意这其實是线性方程组Ax =λx。这表明假设xx = 1,使xAx最大化或(或最小化)的唯一的点正是A的特征向量

线性代数是代数学的一个分支主要处理线性关系问题。线性关系意即数学对象之间的关系是以一次形来表达的例如,在解析几何里平面上直线的方程是二元一次方程;空间平面的方程是三元一次方程,而空间直线视为两个平面相交由两个三元一次方程所组成的方程组来表示。含有 n个未知量的一佽方程称为线性方程变于量是一次的函数称为线性函数线性关系问题简称线性问题解线性方程组的问题是最简单的线性问题。

线性(linear)指量与量之间按比例、成直线的关系在数学上可以理解为一阶导数为常数的函数

非线性(non-linear)则指不按比例、不成直线的关系,一階导数不为常数

行列式非零矩阵可逆方阵满秩向量组满秩(向量个数等于维数)




我要回帖

更多关于 对称矩阵的行列式计算 的文章

 

随机推荐