功能: 对矩阵进行特征分解

格式: 
[U, D] = Evd(A)

A : 为方阵变量

U : 由特征向量组成的正交矩阵
D : 由特征组成的对角矩阵 

说明:
1、本函数主要执行如下分解操作
$$
A = U D U^{-1}
$$

2、在实数域并不是所有矩阵均可以分解
3、有些情况会分解失败
4、U的每一列是对应特征值的特征向量
5、D对角线对应矩阵特征值
6、对称矩阵的特征值可能小于0,因此D对角线元素可能小于0,注意与奇异值分解的得到的对角阵加以区别

例子:
a = 
[ 6.20750296218670    7.80977527508967    6.99137387191475
  9.07738494178158    6.25692807894988    7.32543014331042
  8.83495556136358    2.54731417286550    3.37043793563286 ]

[u,d] = Evd(a)//执行得到如下结果
u =
[-0.61291375327381    0.52696049847470    0.28076480748802
 -0.65802256059664    0.15076871879236    0.48106506294562
 -0.43742775494204   -0.83640984360599   -0.83051040215600 ]

d =
[ 19.5816999958907    0.00000000000000    0.00000000000000
  0.00000000000000   -2.65499076938446    0.00000000000000
  0.00000000000000    0.00000000000000   -1.09184024973680 ]
