﻿功能:矩阵的舒分解(Schur分解)

格式:Schur(a);Schur(a,loopn);Schur(a,loopn,erro)

说明:a为方阵变量.我们知道方阵的Schur分解是把矩阵a分解成a=qt*t*q的形式,其中qt与q是互为转置的矩阵,且q为正交矩阵,t为上三角矩阵.t的对角线元素为a的特征值.执行本函数返回的是t矩阵的值.loopn为循环控制变量,默认为500,erro为误差控制变量,其值在1-16之间,默认为10.

例子:

a =
[ 11.0551527287137   71.0685545443876   55.0673381216206   5.49302664841201   53.9487880905852
  27.3912818298635   41.6194105248989   37.0397641030325   71.2605744466468   67.5347192527422
  38.2129086825125   85.7217742994995   34.9285981780517   84.9352803942446   24.7666571404630
  0.22585745911386   1.09301931275661   44.8080903593488   35.0322951725835   32.6288567076572
  8.50495319278210   68.5980763140126   77.6573113061755   18.3509931985061   54.1977917562229 ]

schur(a)//执行命令
ans =
[ 212.101244342624   131.207813493537   -29.993394254844   86.7739152800068   66.7197301337804
  0.00000000000000   -30.384921202575   -34.862466866375   26.5082695418676   -49.382254721654
  0.00000000000000   43.6486410274716   -31.880037713496   20.2918940149537   -39.235420223935
  0.00000000000000   0.00000000000000   0.00000000000000   -1.3475031695787   -14.270724953314
  0.00000000000000   0.00000000000000   0.00000000000000   0.00000000000000   28.3444661034974 ]

//注意:目前本函数还有问题,有待改善.其结果仅供参考.

//完成时间 By 2011/09/24 17:00:00