schur

语法

schur(obj, [sort])

参数

obj 是一个方阵。

sort 是一个字符串,表示根据所选参数对高特征值进行排序。它的可选参数可以是 'lhp'(特征值为小于零的实数),'rhp'(特征值为大于零的实数),'iuc'(特征值的绝对值<=1.0),或 'ouc'(特征值的绝对值>1.0)。
注: 若没有指定该参数,表示不进行排序。

详情

计算矩阵的 Schur(舒尔)分解。

假设输入方阵为 A:

  • 如果不指定 sort 参数,返回两个矩阵:T 矩阵(矩阵 A 的 schur 形式,正交矩阵)与酉矩阵 Z(对于实数矩阵而言,Z 的转置矩阵与其逆矩阵相等),使得 A = ZTZ-1 。

  • 如果指定 sort 参数,返回结果还包含一个整数,表示满足排序条件的特征值的数量。

例子

m=matrix([[0,0,1],[2,1,0],[2,2,1]]);
T,Z=schur(m)
T;
#0 #1 #2
2.658967 1.424405 -1.929334
0 -0.329484 -0.490637
0 1.311789 -0.329484
Z
#0 #1 #2
0.727116 -0.601562 0.330796
0.528394 0.798019 0.289768
0.438294 0.035904 -0.898114
T,Z,s=schur(m, 'lhp');
T;
#0 #1 #2
-0.329484 1.570974 2.251318
-0.40969 -0.329484 -0.092398
0 0 2.658967
Z
#0 #1 #2
0.703818 -0.632169 0.324042
0.509043 0.766983 0.390655
-0.495495 -0.109999 0.861618
s
# output
2
T,Z,s=schur(m, 'rhp');

s;
# output
1

m=matrix([[0,0,9],[-2,1,0],[2,2,1]]);
T,Z,s=schur(m, 'iuc');
s;
# output
0

T,Z,s=schur(m, 'ouc');
s;
# output
1