1、function[coeff,score,latent,tsquare]=princomp(x,econFlag)%PRINC晦倘佳鳎OMPPrincipalComponentsAnalysis(PCA)fromrawdata.%COEFF=PRINCOMP(X)performsprincipalcomponentsanalysisontheN-by-P%datamatrixX,andreturnstheprincipalcomponentcoefficients,also
3、[COEFF,SCORE]=PRINCOMP(X)returnstheprincipalcomponentscores,i.e.,therepresentationofXintheprincipalcomponentspace.RowsofSCOREcorrespondtoobservations,columnstocomponents.
5、real_T kk,delta;real_T sat(real_T x){if (fabs(x)>delta) return (sign(x)); else return (kk*x); };
7、[~,maxind]=max(abs(coeff),[],1);d=size(coeff,2);colsign=sign(coeff(maxind+(0:p:(d-1)*p)));coeff=bsxfun(@times,coeff,colsign);ifnargout>1score=bsxfun(@times,score,colsign);end
