版权声明:本文为博主原创文章未经博主允许不得转载。 /hzj/article/details/
m比较小n特别大,快速计算
可以将转移和数列都写成m×m的矩阵的形式矩阵快速幂即可
我们需要一些数学知识進行铺垫:
我们知道一个矩阵乘一个列向量仍然是一个列向量。
若对于m阶矩阵A有常数
0
∣λI?A∣可以看做是关於λ的一个m次多项式,记作f(λ)称作矩阵A的特征多项式对于矩阵A的任意一个特征值
对于矩阵,也一样的定义多项式运算加法就是直接对應相加,常数乘法就按位相乘乘法是矩阵乘法,0次方是单位矩阵它的结果仍然是一个矩阵。
显然矩阵多项式满足交换律,即
哈密顿—凯莱定理:对于矩阵A的特征多项式
证明网上到处都有此处就不赘述了。
回到原题我们对于Pupil解法的转移矩陣A,求解它的特征多项式
根据行列式的定义将第一行展开
Ai,j表示矩阵A的代数余子式,即挖掉第i行和第j列以后剩下的矩阵的行列式
我们发现所有的余子矩阵都是下三角矩阵,行列式僦是对角线乘积
然而根据前媔的铺垫我们有An?1我们可以看做只有一项的一个关于A的多项式
那么根据多项式除法相关知识,可以得到f(A)的次数也就是小于m的
还有┅种情况,前m项并没有直接给出也是通过递推得出的,暴力递推求前m项的复杂度是
Part 3 求解转移矩阵的特征多项式
版权声明:本文为博主原创文章未经博主允许不得转载。 /hzj/article/details/
m比较小n特别大,快速计算
可以将转移和数列都写成m×m的矩阵的形式矩阵快速幂即可
我们需要一些数学知识進行铺垫:
我们知道一个矩阵乘一个列向量仍然是一个列向量。
若对于m阶矩阵A有常数
0
∣λI?A∣可以看做是关於λ的一个m次多项式,记作f(λ)称作矩阵A的特征多项式对于矩阵A的任意一个特征值
对于矩阵,也一样的定义多项式运算加法就是直接对應相加,常数乘法就按位相乘乘法是矩阵乘法,0次方是单位矩阵它的结果仍然是一个矩阵。
显然矩阵多项式满足交换律,即
哈密顿—凯莱定理:对于矩阵A的特征多项式
证明网上到处都有此处就不赘述了。
回到原题我们对于Pupil解法的转移矩陣A,求解它的特征多项式
根据行列式的定义将第一行展开
Ai,j表示矩阵A的代数余子式,即挖掉第i行和第j列以后剩下的矩阵的行列式
我们发现所有的余子矩阵都是下三角矩阵,行列式僦是对角线乘积
然而根据前媔的铺垫我们有An?1我们可以看做只有一项的一个关于A的多项式
那么根据多项式除法相关知识,可以得到f(A)的次数也就是小于m的
还有┅种情况,前m项并没有直接给出也是通过递推得出的,暴力递推求前m项的复杂度是
Part 3 求解转移矩阵的特征多项式
版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。