gr_algorithm.m

主成分分析和偏最小二乘SquaresPrincipal成分分析（ PCA ）和偏最小二乘（ PLS ）

function gr=gr_algorithm(modelinfo);
% Main computational algorithm for global reconstruction

% last modified 28.09.05

coeffs=modelinfo.coeffs;
funcs=modelinfo.funcs;
left_part=modelinfo.left_part;

crd_number=length(modelinfo.x.name);
for i=1:crd_number
    assignin('caller',sprintf('x%d',i),modelinfo.x.data(:,i));
end
A=zeros(length(modelinfo.x.data(:,1)),length(funcs));
for i=1:length(funcs)
    A(:,i)=evalin('caller',funcs{i});
end
for i=1:crd_number
    evalin('caller',['clear ' sprintf('x%d',i)]);
end
B=modelinfo.x.data(:,str2num(left_part(2:end)));

C=inv(A'*A)*(A'*B);
err=B-A*C;
clear A;
gr.err=err;
gr.gm=modelinfo.gm;
gr.coeffs_txt=coeffs;
gr.funcs_txt=funcs;
gr.left_part=left_part;
gr.coeffs=C;
gr.B=B;
gr.x=modelinfo.x;

% f=figure('Units','characters','Name','Global model approximation error','NumberTitle','off','color',[1 1 1]);
% subplot(311);
% plot(get_time(modelinfo.x(1)),err),grid on
% xlim([min(get_time(modelinfo.x(1))) max(get_time(modelinfo.x(1)))]);
% ylabel('model error');
 disp_result_gr_model(gr);