www.gusucode.com > mbcmodels 工具箱 matlab 源码程序 > mbcmodels/@localpspline/eval.m

    function y=eval(s,x)
%EVAL evaluate model

%  Copyright 2000-2008 The MathWorks, Inc. and Ford Global Technologies, Inc.



y=zeros(size(x));

if isfinite(s.knot)
   m= x <= s.knot;
   y(m) = polyval_mex(s.polylow, ( x(m)-s.knot )) ;
   y(~m) = polyval_mex(s.polyhigh, ( x(~m)-s.knot ));
else
   y(:)=NaN;
end