www.gusucode.com > mbcmodels 工具箱 matlab 源码程序 > mbcmodels/@localtruncps/x2fx.m
function X= x2fx(ts,x) % TRUNCPS/X2FX % Copyright 2000-2004 The MathWorks, Inc. and Ford Global Technologies, Inc. X= zeros(length(x),ts.order+length(ts.knots)); X(:,ts.order)= 1; if ts.order>1 X(:,ts.order-1)= x; for i=ts.order-2:-1:1 X(:,i)= x.*X(:,i+1); end end m= ts.order; for i=1:length(ts.knots) rhs= x>ts.knots(i); X(rhs,i+m)= (x(rhs)-ts.knots(i)).^(m-1); end