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