function res = deterministic_vs_prob rand('seed',42); n = 10; x0 = rand(n,1); A = generate_connected(10); A = graph_circular(10); x0 = ((1:10)' / 10) + 0.2*rand(n,1); % epsilons = [0.2 0.1 0.05 0.025 0.0125 0.0064] epsilons = 0.02; % colors = {'r','g','b','k','m'}; for i=1:numel(epsilons) epsilon=epsilons(i); for strategy=[6 1 2 3 5] res{i,strategy} = cons00(A, x0, epsilon, strategy, 2000); end % p = plot(log_or_small(rprob.mse, -30), '.'); % set(p, 'color', colors{mod(i, numel(colors))+1}); % set(p, 'markersize', 2); % pause(0.1) % res{i,5} = rdet; % res{i,2} = rprob; end % for i=1:numel(epsilons) % legends{i} = sprintf('\\epsilon = %f',epsilons(i)); % end % legend(legends); function r = log_or_small(k, min) for i=1:numel(k) r(i) = log(max(k(i),exp(min))); end