Mercurial > hg > trimatlab
changeset 6:0f16c069423b
Added randomise.m to provide random init; also added minimum radius to mt_get_transmat_near
| author | samer |
|---|---|
| date | Mon, 06 Feb 2012 15:55:08 +0000 |
| parents | 80e9103b057f |
| children | 917d5d11aa48 |
| files | mt_get_transmat_near.m private/randomise.m |
| diffstat | 2 files changed, 14 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- 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)));
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/private/randomise.m Mon Feb 06 15:55:08 2012 +0000 @@ -0,0 +1,11 @@ +function s=randomise +% randomise - get entropy from time and reset random number generators + + +m=[5000,700000,1000000,1000000,1000000,10000000]; +rand('state',sum(m.*clock)); +randn('state',sum(m.*clock)); +if nargout>0, + s=getrndstate; +end +