矩阵分析 学习笔记3 多项式矩阵 jordan标准型

多项式矩阵

就是说这个矩阵里面的各个元素都是多项式，多项式的主角是类目大（自变量）。

多项式矩阵的秩

0多项式就是完全0的那种，就一个0，类目大都没有了。

多项式矩阵的秩和带一个类目大进去变成普通矩阵的秩不是一回事哦，毕竟类目大可以取到一个特定的值使行列式的值为0，刚好凑到0，秩就变小了。

单位模阵（可逆多项式矩阵）

单位模阵就是可逆多项式矩阵

假如这个多项式矩阵是可逆的，那么它的行列式肯定是一个固定的值，还不能是0！证明如下：

多项式矩阵初等变换

多项式矩阵的Smith标准形

deg就是次数的意思。

这里介绍一下多项式的知识

多项式知识复习

带余除法

f(x)=g(x)q(x)+r(x)，其中f(x)是被除式子，g(x)是除数，q(x)是我们要求的商，r(x)是余数，它的次数显然要比g(x)要小。那么这种式子要怎么算呢？

就是一直凑，先从f(x)中凑出一个可以整除（刚好多几个x的那种）的式子，然后再从剩下的式子里面继续凑，直到剩下的式子次数比g(x)小为止。

质因式分解

f(x)=q1(x)q2(x)……     。q(x)是不可约的多项式

回到之前的知识。我们可以通过多项式除法，求出a** = a11q(x)+r(x)的q(x)和r(x)来，这样我们就可以通过初等变换把原本多种多样的a**降次，变成r(x)，这样整体的多项式矩阵的次数就降低了，达到了

注意最后的结果要首1化（系数为1）

这一点也很好证明，做法如下：

1.挑选一个次数最低的元素（当然不能是0），用行列变换把他换到a11位置。

2.再使用之前提到的方法，把a**降次

3.用a11把一行一列变成0.剩下的部分组成一个矩阵，都可以被a11整除，又可以重复上述步骤。

举个例子：

第一步就是要找出那个次数最低的（0的次数是无穷哦）然后化简。其实一开始的时候已经很标准了，就是(a+1)^2不听话，所以把他拉上来减，看看能不能有变化。

这里提出一个问题：这个标准形是唯一的吗？（指的是r的大小，不是说式子也要一模一样）

注意这些子式都是不包括0的。

想要证明这一点，我们需要证明：经过一次初等行变换，行列式因子不变。实际上交换两行，各阶子式只相差一个正负号，首项归一就一样了。

这样我们也获得了一个新的求标准型的方法。

用不变因子获得smith标准型

感觉这种方法是最简单的，以后就用这种方法了

数域上矩阵的特征矩阵

特征矩阵很久以前就学过了，为啥忽然讲这个？：

多项式矩阵的等价问题更简单！特征矩阵作为多项式矩阵等价，那就相当于数字矩阵相似。

多项式矩阵的次数

可以看到一个多项式矩阵的秩其实就是看各个元素中次数最大的一个。

和多项式一个道理，两个相乘，次数就是相加。当然要Ar非奇异（满秩）别给我乘没了

和多项式带余除法一个道理，要么余数为0，要么余数次数要比除数要小。

smith标准型的求法上面讲过

n = n1 + n2 ……，左右同时减p就成了。

就比如第一个矩阵，其不变因子就是1，（a+1）^2,a(a+1)^2，给他们做质因式分解（不包括1）得到(a+1)^2,a,(a+1)^2。（别把平方也拆了）这样你会发现，这些因子乘起来就是原来的那个矩阵的行列式。

注意smith标准型对角线上的那些式子都是可以整除上一个式子的。因此（2）实际上只有一种可能。

矩阵可以求出各阶行列式因子，由各阶行列式因子可以求出不变因子，再把不变因子们拆开就得到了初等因子组。

复数域上的矩阵jordan标准型

这里只考虑原本那些不变因子中的一个，相当于第二标准型中的一小块。

完全搞不懂了！！！初等因子组到底干嘛用的？

大概就是这么个流程：看见一个普通矩阵，找出其特征矩阵，然后求smith标准型（两种方法都随便），从对角线的不变因子上我们可以求出初等因子们，从每一个初等因子中我们都可以直接求出诺当标准型当中的一个小块，使用的就是上面这个原则，看它减去几，是几次方，次方决定了小块多大，减去几决定了对角线上的元素。次对角线的元素总是1，假如次数只有1的话那就没有次对角线了。

费老大劲找到了求矩阵相似于一个诺当标准型的方法，拿来干嘛？：

广义特征向量

推理过程：

乘自己对应的向量就是0，乘下一个向量就是非0，这样我们就找到了P（也就是特征向量们）注意这个是发生在诺当标准型的一个小块里面的，因此Ci都是一样的。

举个例子：

接下来就是解方程环节，都是老操作了。正好这个矩阵的两个C还都是一样大，更好算了。

不要一开始就把P2定死，解第三个式子的时候可能会有非齐次方程组特有的问题，要保证下面那几行能等于0！确定了P2有啥限制之后再取。

算完之后可以验证一下特征矩阵是不是满秩的，也可以把P带过去看看能不能把A变成J。P不唯一。

诺当标准型的应用：

我们要求出xi与t之间的关系函数。x = f(t)这种。

这个y是怎么弄出来的？搞不懂啊。接下来要讲别的了，我裂开了。

应该就是说把x做一些变换，弄出一个y来，先求y与t之间的式子，再改回x?把dy/dt的式子求出来（尽量简单）然后就直接用微分方程公式求出y，相应的求出x。

我好像懂了？先是把微分方程写成AX的形式，然后关注这个A，求出其诺当标准型和对应的特征矩阵什么的，然后就得到了AP = PJ。

然后令X = PY（也就是我们想要的把x转化为Y的式子），两边对t求导，也就是dx/dt = P* dy/dt，别忘了dx/dt我们已知，就是Ax,因此就是AX = P * dy/dt，但X是PY啊，也就是APY = P * dy/dt

AP就是PJ，特征矩阵和诺当标准型我们都知道，那么PJY = P * dy/dt。都是左乘一个P，那抵掉了，也就是JY = dy/dt，这样我们就求出了一个单纯关于y的微分方程。由于J是很简单的矩阵，得到的微分方程也很简单，自然可以直接求出来。y解出来了，x也就解出来了。