首页 > 编程知识 正文

矩阵顺序主子式计算,矩阵主子式怎么求

时间:2023-05-03 22:09:15 阅读:221722 作者:3691

1.名词介绍:

特征多项式:在数域P上的某个方阵A,则行列式即为其特征多项式,简记为;

逆序数:1,2,...n的某个排序相对自然序列相反的个数;记为:

主子式:某个行列式的部分列和部分行:行号和列号一样交线处的元素按相对位置不变交出来的方阵的行列式:

按照子阵的表达方式:即行列交出来的子阵,其行列式记为;

2.下面介绍两个引理:

引理1:行列式乘积表达式:

这个引理的证明可参看《高等代数》帅气的八宝粥,震动的小笼包 著,第四版的56页

引理2:行列式求导:

对于某个行列式求导等于对其每一行求导并且将其还原到原矩阵并求和:

也即:

这个引理很显然,利用引理1,和乘积求导的公式,并把每次被求导的元素按照位置不变的带回原来的行列式就能得到引理2

定理1:

对于一个实数域上n阶方阵A,其特征多项式被表示为:

,其中

证明:利用多项式sdjc展开式可以得到以下关系:

利用行列式求导法:

1由于都在主对角线上,按行降阶:

即表示成了所有n-1阶主子阵的特征多项式之和

对于任意一个n-2阶主子阵,一定真包含且只包含在两个n-1阶主子阵中:

即只真包含在与

中,因此对于要除以2才是

依次类推:

包含且只包含在n-k个k+1阶主子阵中,故要除以才是:

定理1证毕;

证明到这里我们可以看到方阵的迹就是这里k取1

利用定理1我们可以得到关于方阵主子式的又一个性质:这在最优化的海森矩阵判断正定性里有部分应用:

定理2:

对于实对称矩阵A是半正定的充分必要条件是A的一切主子式大于等于0:

实对称,取其主子式:

则存在一个可逆阵其中B,D是除去的一定规格的矩阵

事实上:对于A的行与列做对等的第二类初等变换(交换)就可得到该矩阵:

这里运用到一个常用的求实对称矩阵合同变换矩阵的方法:

构造:

矩阵,这里的表示可以取任意n阶方阵,每次交换两行的同时交换两列:最终的得到A下方位置的矩阵就是C:

事实上,这里根据初等矩阵的理论:每次对M做一对对等的初等行变换:相当于给M左乘一个初等阵

结果为:

在此结果上做一次对等的列变换:等价于右乘一个初等阵::结果为:

依次类推可得经过n次该变换的矩阵为:

即根据初等矩阵的性质:即可证明经过m次该对等的行列初等变换的和A合同;

必要性:

按照交换排列的方法:就比如选择排序,操作M,使得被选中的行列在上或左,得到的的前k行k列交出来的矩阵是排序之前的

,注意这里的形式相同,表达着不同的含义:前者是m经过m次对等初等行列变换后原本的A位置的矩阵;而是原本A矩阵的k阶主子式;

现在证明这个k阶主子阵也是半正定的:

取向量:由于合同变换前后的矩阵正定性相同,即也是半正定的故对于分块:

即也为半正定,故其行列式大于等于零:如若不然则存在负的特征根;注意这个反之不然;

充分性:

因为A的任意主子式都大于或等于0,故考虑的正定性:这里可取任意小

判断正定性的方法其一就是其充要条件:所有顺序主子式都大于等于0,那么这里的的任意一个k阶顺序主子式即:

的展开式为(***),而同样考虑的特征多项式

其展开式为

化为首一多项式为:

由定理1:

故:

故在时恒正

即推导出任意阶的顺序主子式大于0,即正定;

故考虑非零实向量:对fkdsg成立,然而其一定可以看成

关于的有限阶多项式,一定在R上连续,令从正向趋于0可得

这就证明了A的半正定性

定理2证毕

版权声明:该文观点仅代表作者本人。处理文章:请发送邮件至 三1五14八八95#扣扣.com 举报,一经查实,本站将立刻删除。