下面是将一个矩阵化成对角矩阵的matlab代码
其中,B是原矩阵,P是变换矩阵,H是得到的以B的特征值为主对角线的对角矩阵。
我现在的问题是:目的是将原矩阵变成以其特征值为斜对角线的变换矩阵,谢谢大家。
(matlab、maple实现都可以)
syms b c d
B=[c*b/d,-b;c,-d]
B =
[ (b*c)/d, -b]
[ c, -d]
>> [P H]=jordan(B)
P =
[ d/c, b/d]
[ 1, 1]
H =
[ 0, 0]
[ 0, (- d^2 + b*c)/d] |