parafacc.m

多维数据处理：MATLAB源程序用于处理多维数据

% all data preprocess by prestd
%gc1=gc';%input('gc=');
%gp1=gp';%input('gp=');
%cc1=cc';%input('cc=');
dimx=[1,49,10];
%dimy=[16,29,9];
%[ac2,meanac,stdac,cc1,meancc,stdcc]=prestd(gc1,cc1);
%ap2=trastd(gp1,meanac,stdac);
%[xc,mgc,sgc]=nprocess(cc,dimxc,centc,scalc);
%gx=ac2';
gx=c;
%gy=ap2';
%cc2=cc1';
%factors=parafac(c,dimx,4);
Options(1)=1e-8
Options(2)=1;
%Options(3)=1;
%Options(4)=[];
%Options(6)=[];
%const=[0 0 1];
fixmode=[1 1 1];
fac=1;
factors=parafac(gx,dimx,fac,fixmode);
[Acal Bcal Ccal]=fac2let(factors,dimx,fac);
M=ppp(Bcal,Ccal)';
%b=Acal\cc2;
%Apre=gy/M;
%Y=Apre*b;
%PC1=poststd(Y',meancc,stdcc);
%PC=PC1'
%cp=cp%input('cp=')
%cp=nprocess(cp,dimyc,centc,scalc,mgc,sgc)
%M=PC-cp;
%N=M.^2;
%F=sum(N);
%H=sum(F,2);
%E=PC.^2;
%G=sum(E);
%J=sum(G,2);
%RSD=sqrt(H./J)
%ERSD=sqrt(F./G) %ERSD is every relational standard error