K-L变换

(一)K-L变换的基本原理

① 对D维随机向量x\in R^D,用一个完备的正交归一向量系u_j(j=1,2,...\infty )来展开:

                                                  x=\sum^{\infty }_{j=1}c_ju_j \rightarrow c_j=u^T_jx

将x表示成u_j的线性组合,其中:

                                                         u^T_iu_j=\left\{\begin{matrix} 1, &i=j \\ 0,&i \neq j \end{matrix}\right.

② 当只有有限的d项来逼近x,即\hat x=\sum_{j=1}^d c_ju_j。则与原向量的均方误差表示为:

                                  e=E[(x-\hat x)^T(x-\hat x)]=\sum_{j=d+1}^\infty u_j^TE(xx^T)u_j

若记\psi =E(xx^T),则e可以写成下式:

                                                             e = \sum_{j=d+1}^\infty u_j^T\psi u_j

③ 要最小化这一均方误差,就是优化下列问题:

                                                     \left\{\begin{matrix} min &e = \sum_{j=d+1}^\infty u_j^T\psi u_j \\ s.t.& u_j^Tu_j=1,\forall j \end{matrix}\right.

得到无约束的目标函数为:

                                                 g(u) = e-\sum_{j=d+1}^\infty \lambda_j(u_j^Tu_j-1)

④对各个向量求偏导,令其为0,即:

                                         \frac{\partial}{\partial u_j}g(u)=0\rightarrow (\psi -\lambda_jI)u_j=0\rightarrow \psi u_j=\lambda_ju_j

⑤ 因此均方误差记作:

                               e = \sum_{j=d+1}^\infty u_j^T\psi u_j=\sum_{j=d+1}^\infty u_j^T\lambda_j u_j=\sum_{j=d+1}^\infty \lambda_j u_j^Tu_j=\sum_{j=d+1}^\infty \lambda_j

把矩阵\psi的特征值从大到小进行排列,选择前d个特征值对应的特征向量会使误差最小。

(二) 用于监督模式识别的K-L变换

算法步骤:

① 计算总类内离散度矩阵S_w

                                             S_w=\sum_{i=1}^cp_i\Sigma _i,\Sigma _i=E[(x-u_i)(x-u_i)^T]

其中各类w_i的先验概率为p_i,均值为u_i,协方差矩阵为\Sigma_i

②用S_w作为产生矩阵进行K-L变换,求解S_w的特征值\lambda_i和对应的特征向量v_i,得到一组新的特征y_i = v_i^Tx

③计算新特征的分类性能指标,并按大小进行排序。即:

                                        J(y_i)=\frac{v_i^TS_bv_i}{\lambda_i},J(y_1) \geq J(y_2) \geq...\geq J(y_D)

其中,S_b=\sum_{i=1}^cp(w_i)(u_i-u)^T(u_i-u)u_iu分别是第i类的均值和总体均值。

④ 满足③中性能指标的前d个新特征y_1,y_2...y_d(当有J(y_1) \geq J(y_2) \geq...\geq J(y_d)...\geq J(y_D)),构成新的特征(y_i = v_i^Tx)。而相应的v_i构成了特征变换矩阵V=[v_1 v_2...v_d]

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值