Spectrum

计算矩阵的谱，作为来自乘积 AT*A 的本征值集合。

vector matrix::Spectrum()

返回值

矩阵的谱作为矩阵本征值的向量。

举例

double MatrixCondSpectral(matrix& a)

  {

   double norm=0.0;

   vector v=a.Spectrum();



   if(v.Size()>0)

     {

      double max_norm=v[0];

      double min_norm=v[0];

      for(ulong i=1; i<v.Size(); i++)

        {

         double real=MathAbs(v[i]);

         if(max_norm<real)

            max_norm=real;

         if(min_norm>real)

            min_norm=real;

        }

      max_norm=MathSqrt(max_norm);

      min_norm=MathSqrt(min_norm);

      if(min_norm>0.0)

         norm=max_norm/min_norm;

     }



   return(norm);

  }