www.gusucode.com > mbctools 工具箱 matlab 源码程序 > mbctools/@mdev_local/pevinit.m

    function TS= pevinit(mdev,ind);
% MDEV_LOCAL/PEVINIT

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




if nargin>1
	TS= mdev.TwoStage{ind};
    sol= ind==1;
else
	if mle_best(mdev)
		TS= mdev.MLE.Model;
		sol= 0;
	else
		sol=1;
		TS= mdev.TwoStage{1};
	end
end

[Xg,Yrf,Sigma]= mledata(mdev,1);

TS= pevinit(TS,Xg,Yrf,Sigma,sol);