diff mt_get_transmat_near.m @ 6:0f16c069423b

Added randomise.m to provide random init; also added minimum radius to mt_get_transmat_near

--- a/mt_get_transmat_near.m	Sun Feb 05 18:20:17 2012 +0000
+++ b/mt_get_transmat_near.m	Mon Feb 06 15:55:08 2012 +0000
@@ -14,17 +14,18 @@
 % The selected transition matrix will be shown using mt_show_transmat,
 
 
-function [T,P,I]=mt_get_transmat_at(Sys,Id,K,X,Y)
+function [T,P,I]=mt_get_transmat_near(Sys,Id,K,X,Y,MinRad)
 	logK=log(K);
 	normpos = Sys.map([X;Y])';
 	target = normpos*logK;
 	tmats = Sys.transmats{K};
 	info = Sys.info{K};
 	L = size(info,1);
+	if nargin<6, MinRad=0.01; end
 
 	% distance from target
 	d2 = sum((info(:,1:2) - repmat(target,L,1)).^2,2);
-	Mask = d2<=max(min(d2)*2,0.01);
+	Mask = d2<=max(min(d2)*2,MinRad);
 
 	M = find(Mask);
 %	J = M(argmax(info(M,3)));