samer@0: % mt_init - Initialise Melody Triangle system.
samer@0: %
samer@0: % mt_init :: 
samer@0: %    A0:nonneg  ~'parameter for sampling',
samer@0: %    B0:nonneg  ~'parameter for sampling',
samer@0: %    A1:nonneg  ~'parameter for sampling',
samer@0: %    B0:nonneg  ~'parameter for sampling',
samer@9: %    L:natural  ~'number of transmats to sample',
samer@9: %    Shuffle:book ~'whether or not to shuffle symbols when returning transmat'
samer@0: % -> mt_system.
samer@0: %
samer@0: % Initial system contains no transition matrices - you
samer@0: % must call mt_ensure with a particular value of K
samer@0: % to sample a set of L transition matrices of that size.
samer@0: %
samer@0: % Initial calibration is equivalent to:
samer@0: %    sys=mt_calibrate(sys, 1:3, [0,1,0;0,0,1]);
samer@0: %
samer@0: % The figure for scatter plots is fixed to figure 50 for now.
samer@0: 
samer@11: function Sys=mt_init(A0,B0,A1,B1,L,Shuffle,ErgMethod,Tol)
samer@11: 	if nargin<7, ErgMethod=1; end
samer@11: 	if nargin<6, Shuffle=0; end
samer@11: 	if nargin<8, Tol=0.001; end
samer@11: 	Sys.sample_transmats = @(k,l)sample_transmat_hdp(A0,B0,A1,B1,k,l);
samer@0: 	Sys.transmats = {};
samer@0: 	Sys.info      = {};
samer@0: 	Sys.refpoints = [0,0;1,0;0,1]';
samer@0: 	Sys.fig       = 50;
samer@9: 	Sys.shuffle   = Shuffle;
samer@11: 	Sys.ergmeth   = ErgMethod;
samer@11: 	Sys.L         = L;
samer@11: 	Sys.tol       = Tol;
samer@0: 	Sys = mt_calibrate(Sys);
samer@0: end
samer@11: