www.pudn.com > trackingdemos.zip > plot_res_prob.m

% plot_res_prob.m 
 
function fig = plot_res_prob() 
r = linspace(.1, 4, 40); 
b = linspace(.1, 4, 40); 
 
for i=1:length(r) 
    for j=1:length(b) 
        if r(i)>3 | b(j)>3 
            prob(i,j)=1; 
        elseif r(i)>1 & b(j)>1 
            prob(i,j) = (r(i)-1)*(b(j)-1)/4; 
        else 
            prob(i,j) = 0; 
        end 
    end 
end 
fig = figure; 
meshc(r, b, prob); 
xlabel('\Deltar/C_r'); 
ylabel('\Deltab/C_b'); 
zlabel('Resolution Probability');