典型相关分析(CCA)附算法应用与程序文件.doc

《典型相关分析(CCA)附算法应用与程序文件.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、......典型相关分析专业技术资料......摘要利用典型相关分析的思想,提出了解决了当两组特征矢量构成的总体协方差矩阵奇异时,典型投影矢量集的求解问题,使之适合于高维小样本的情形,推广了典型相关分析的适用范围.首先,探讨了将典型分析用于模式识别的理论构架,给出了其合理的描述.即先抽取同一模式的两组特征矢量,建立描述两组特征矢量之间相关性的判据准则函数,然后依此准则求取两组典型投影矢量集,通过给定的特征融合策略抽取组合的典型相关特征并用于分类.最后,从理论上进一步剖析了该方法之所以能有效地用于识别的内在本质.该方法巧妙地将两组特征矢量之间的相关性特征作为有效判别信息,既达到了信息融合之目

2、的,又消除了特征之间的信息冗余,为两组特征融合用于分类识别提出了新的思路.专业技术资料......一、典型相关分析发展的背景随着计算机技术的发展,信息融合技术已成为一种新兴的数据处理技术,并已取得了可喜的进展.信息融合的3个层次像素级、特征级、决策级。特征融合，对同一模式所抽取的不同特征矢量总是反映模式的不同特征的有效鉴别信息，抽取同一模式的两组特征矢量，这在一定程度上消除了由于主客观因素带来的冗余信息，对分类识别无疑具有重要的意义典型相关分析(CanoniealComponentAnalysis:CCA)是一种处理两组随机变量之间相互关系的统计方法。它的意义在于:用典型相关变量之间的关系

3、来刻画原来两组变量之间的关系!实现数据的融合和降维!降低计算复杂程度。二、典型相关分析的基本思像CCA的目的是寻找两组投影方向,使两个随机向量投影后的相关Txcc...c性达到最大。具体讲,设有两组零均值随机变量1，2，pTydd...d和1，2，qTCCA首先要找到一对投影方向1和1,使得投影v11yT和u11x之间具有最大的相关性,u1和v1为第一对典型变量;同理,寻找第二对投影方向2和2,得到第二对典型变量u2和v2,使其与第一对典型变量不相关,且u2和v2之间又具有最大相关性。这样下去,直到x与y的典型变量提取完毕为止。从而x与y之间的相关性分析,只需通

4、过分析少数几对典型变量的关系即可达到目的。专业技术资料......三、CCA算法详解(1)(2)(3)考虑到:的极值只与和的方向有关,而与它们的大小无关,为了得到唯一解不失一般性,加入限制条件:TTSS1(4)xxyy问题变为在约束条件式下,求使准则函数式取最大值的典型投影矢量对和求解上述优化问题,可定义拉格朗日函数:(5)分别对和求导数,并令为零,得到：(6)专业技术资料......(7)(8)(9)再对上式两端分别左乘和得：(10)(11)(12)(13)T对H进行奇异值分解：HUSV令TTTTTT分别将1x，2x，…dx与1x，2x，…dx看做是

5、变换后的特征分量：(14)(15)（16）投影后的组合特征用于分类，其中变换矩阵为：专业技术资料......（17）四、典型相关分析应用实例欲研究儿童形态与肺通气功能的关系，测得某小学40名8~12岁健康儿童（身高X1，体重X2，胸围X3）与肺通气功能（肺活量Y1，静息通气Y2和每分钟最大通气量Y3），分析儿童形态和肺通气指标的相关性，确定典型变量的对数。x1=[140.6,135.7,140.2,152.1,132.2,147.1,147.5,130.6,154.9,142.4,136.5,162,148.9,136.3,159.5,165.9,134.5,152.5,138.2,144

6、.2];x2=[43.7,39.5,48,52,36,45,47,38,48,42,38,58,42,33,49,55,41,53,35.5,42];x3=[77,63,75,88,62,78,76,61,87,74,69,95,80,68,87,93,61,83,66,76];y1=[2.6,2,2.6,2.8,2.1,2.8,3.1,2,2.9,2.33,1.98,3.29,2.7,2.4,2.98,3.1,2.25,2.96,2.13,2.52];y2=[7,7,6.1,10.1,7.4,9.25,8.78,5.31,10.6,11.1,7.77,3.35,10.1,7.8,11.7

7、7,13.14,8.75,6.6,6.62,5.59];y3=[108,91,101,112,97,92,95,77,80,76,49,58,82,76,88,110,75,71,105,82];（1）仿真结果分析结：（实验平台：Matlab2014，程序见附录）R1=0.9282R2=0.5302R3=0.0081专业技术资料......R1=0.9282R2==0.5302R3=0.0081专业技术资料......(