srm_ab_model.m

来自「轨迹聚类工具箱及其相关说明和文章（台风路径分类等）。如有问题可与wyljess@」· M 代码 · 共 48 行

function M = srm_model()
%SRM_AB_MODEL  Example simulation model
%  M = SRM_AB_MODEL()

% Scott J Gaffney   2 February 2003
% Department of Information and Computer Science
% University of California, Irvine.

%% Setup spline basis model
% range of x-axis
b_int = -4;    % beginning of data interval
e_int = 25;    % end of data interval
len_int = e_int-b_int+1; % length of interval

degree = 3;  % degree of polynomial pieces, this can be calculated as Kn-L-1
             % where L is size(M.Mu,1)
Kn = max(ceil(len_int/3),1);  % number of knots
knots = linspace(b_int,e_int,Kn);
M.knots = addendpts(knots,degree+1);
Kn = length(M.knots);  % reset M to (possibly) new value
M.order = degree+1;

% Mu: L-K-D or P-K-D
M.Mu(:,1) = [10.0000    8.0000    9.0000    6.6667    6.0000 ...
    4.0000    3.0000    5.0000    2.0000  2  1  0];
M.Mu(:,2) = [4.0000    1.1111    3.3333    6.6667    6.0000 ...
    13.3333   16.6667   18.8889   10.0000  9  8  7];


% Alpha: K-1
M.Alpha = [.6 .4]';
M.K = length(Alpha);

% R: K-D
M.R = [1.5; .5];

% S: K-D
M.S = [.7; .3];

% Sigma: K-D
M.Sigma(1,1) = .05;
M.Sigma(2,1) = .05;


if (nargout==0)
  assignin('caller','SimModel',M);
  clear M;    % keeps default output from printing on-screen
end