samer@0: % mt_resample - Sample or resample melody triangle transition matrices of given size
samer@0: %
samer@0: % mt_resample ::
samer@0: %    mt_system  ~'initial system'
samer@0: %    natural    ~'size of transition matrices to resample'
samer@0: % -> action mt_system.
samer@0: %
samer@0: % A new set of transition matrices will be sampled and a 3D information space
samer@0: % scatter plot generated in the figure determined by the initial call to
samer@0: % mt_init.
samer@0: 
samer@0: function Sys=mt_resample(Sys,K)
samer@0: 	[TT, II] = Sys.sample_transmats(K);
samer@0: 	figure(Sys.fig); % 'name','info space'); 
samer@0: 	Sys.transmats{K} = TT;
samer@0: 	Sys.info{K}      = II;
samer@0: 	Sys.hScat{K}     = scatc(II,II(:,3),16); 
samer@0: 	axis on; box on; grid off; 
samer@0: 	lc=[0.4,0.4,0.4];
samer@0: 	set(gca,'Color','none');
samer@0: 	set(gca,'XColor',lc);
samer@0: 	set(gca,'YColor',lc);
samer@0: 	set(gca,'ZColor',lc);
samer@0: 	xlabel('entropy rate');
samer@0: 	ylabel('redundancy');
samer@0: 	zlabel('pred-info rate');
samer@0: 	rotate3d on;
samer@0: end