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

    function J= CalcJacob(ts,x);
% TRUNCPS/CALCJACOB

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




nk= length(ts.knots);

J= zeros(length(x),numParams(ts.xreglinear)+nk);

X= x2fx(ts,x);

J(:,nk+1:end)= X(:,Terms(ts.xreglinear));

m= ts.order;
p= double(ts.xreglinear);
for i= 1:length(ts.knots);
   rhs= x>ts.knots(i);
	if any(rhs)
		J(rhs,i)= J(rhs,i) - (m-1)*p(m+i)*(x(rhs,:)-ts.knots(i)).^(m-2);
	end
end