预览加载中，请您耐心等待几秒...

matlab主成分分析法.pdf

预览
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲，该文档总共11页，到这已经超出免费预览范围，如果喜欢就直接下载吧～

下载提示 文本预览

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

2005 §10.Matlab 1. Matlab() Matlab 1.1 L r11r12r1p  rrLr R=21222p MMMM  L rp1rp2rpp 1 3.5.3rijij=12…pxixj n ∑(xki−xi)(xkj−xj) k=1 rij= nn 22 ∑(xki−xi)∑(xkj−xj) k=1k=1 2 Rrij=rji 58 2005 λI−R=0Jacobi λ=L i(i,2,1,p)L λ1≥λ2≥,≥λp≥0 p L2 λiei(i=,2,1,p)ei=1∑eij=1 j=1 eijeij zi λ i=L p(i,2,1,p) ∑λk k=1 i ∑λk k=1=L p(i,2,1,p) ∑λk k=1 L 85里95%λ1,λ2,,λm …mmp l=p(z,x)=λei,(j=,2,1L,p) ijijiij3 59 2005 3.5.2 L z11z12z1m  zzLz Z=21222m MMMM  zzLz n1n2nm4 2. Matlab MatlabMatlab 2.1 Cwprint.m Cwstd.mCwfac.mCwscore.m 2.2 Cwstd.m里里 Cwfac.m里里 85% 60 2005 Cwscore.m里里 Cwprint.m里里 3. 3.1cwstd.m %cwstd.m, functionstd=cwstd(vector) cwsum=sum(vector,1);% [a,b]=size(vector);%,a,b fori=1:a forj=1:b std(i,j)=vector(i,j)/cwsum(j); end end 3.2cwfac.m %cwfac.m functionresult=cwfac(vector); fprintf(':\n') std=CORRCOEF(vector)% fprintf('(vec)(val)\n') [vec,val]=eig(std)%(val)(vec) newval=diag(val); [y,i]=sort(newval);%yi fprintf('\n') forz=1:length(y) newy(z)=y(length(y)+1-z); end f