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

    function [B0,minb,maxb,OK]= initial(ts,x,y);
%INITIAL

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



nk= length(ts.knots);

k= sort(ts.knots);

xs= sort(x);

r= xs(end)-xs(1);
tol= r/100;

maxb= xs(end-nk)- tol;
minb= xs(nk) + tol;

B0= double(ts);

bk= linspace(minb,maxb,nk+2);

B0(1:nk)= bk(2:end-1);

OK=1;