modelvalidate.m

MDPSAS工具箱是马里兰大学开发的

function [Aprederror,MSE,pcnterr,Apred] = modelvalidate(P,A)

% modelvalidate.m Compare wafer profile predictions from polymodel
%                 to a set of experimental data. Called as
%
%                 [Aprederror,MSE,pcnterr,Apred] = modelvalidate(P,A)
%
% INPUT PARAMETERS
%
%                P : a polymodel object
%                A : vector of wafermap objects
%
% OUTPUT PARAMETERS
%
%       Aprederror : vector of wafermap objects
%              MSE : vector of mean squared errors
%          pcnterr : vector of prediction errors, %
%            Apred : vector of predicted wafermap objects

% Copyright (c) by Raymond A. Adomaitis, 1998-2006

% the polymodel parameters

polymodelparam = get(P,'param');

Apred = A;
Aprederror = A;

for i = 1:length(A)

% get parameters from each wafermap object

    Ap   = get(A(i),'param');
    Akey = keys(Ap);

% now reset the parameter values in the polymodel according
% to the wafer param values

    for j = 1:length(Akey)
        P = set(P,'param',getval(Ap,Akey{j}),Akey{j});
    end

    [P,Apred(i)] = residual(P);
    Aprederror(i) = A(i) - Apred(i);

    MSE(i) = sqrt(mean(Aprederror(i)^2));
    pcnterr(i) = mean(abs(Aprederror(i)))/mean(A(i);

end