Mercurial > hg > map
view multithreshold 1.46/compare.m @ 25:d2c4c07df02c
path problem solved
| author | Ray Meddis <rmeddis@essex.ac.uk> |
|---|---|
| date | Fri, 17 Jun 2011 14:30:28 +0100 |
| parents | f233164f4c86 |
| children |
line wrap: on
line source
f=[100 200 500 1000 2000 5000 10000]; levels=-10:10:90; without=[ 1.44E-10 8.24E-10 1.15E-09 1.43E-09 1.08E-09 4.75E-10 1.43E-09 4.56E-10 2.61E-09 3.63E-09 4.52E-09 3.37E-09 1.50E-09 4.52E-09 1.44E-09 8.24E-09 1.15E-08 1.44E-08 1.06E-08 4.75E-09 1.44E-08 4.56E-09 2.61E-08 4.45E-08 5.16E-08 3.73E-08 1.15E-08 5.16E-08 1.44E-08 7.05E-08 9.57E-08 9.34E-08 7.44E-08 1.98E-08 9.34E-08 4.56E-08 1.28E-07 1.46E-07 1.34E-07 1.04E-07 4.35E-08 1.34E-07 1.06E-07 1.76E-07 2.29E-07 1.99E-07 1.47E-07 5.80E-08 1.99E-07 1.72E-07 2.48E-07 3.65E-07 3.94E-07 3.39E-07 1.31E-07 3.94E-07 2.54E-07 3.59E-07 7.88E-07 9.97E-07 9.07E-07 3.51E-07 9.97E-07 3.82E-07 6.00E-07 2.13E-06 2.85E-06 2.66E-06 1.03E-06 2.85E-06 5.81E-07 1.37E-06 6.30E-06 8.65E-06 8.17E-06 3.30E-06 8.65E-06 ]; with=[ 2.60E-11 1.48E-10 2.09E-10 2.61E-10 1.93E-10 8.54E-11 2.61E-10 8.21E-11 4.69E-10 6.61E-10 8.25E-10 6.10E-10 2.70E-10 8.25E-10 2.60E-10 1.48E-09 2.09E-09 2.61E-09 1.93E-09 8.54E-10 2.61E-09 8.21E-10 4.69E-09 6.61E-09 8.25E-09 6.10E-09 2.70E-09 8.25E-09 2.60E-09 1.48E-08 2.35E-08 3.10E-08 1.96E-08 8.54E-09 3.10E-08 8.21E-09 4.56E-08 7.17E-08 7.64E-08 5.59E-08 1.33E-08 7.64E-08 2.60E-08 1.02E-07 1.20E-07 1.16E-07 9.13E-08 3.13E-08 1.16E-07 7.49E-08 1.51E-07 1.79E-07 1.70E-07 1.25E-07 4.88E-08 1.70E-07 1.38E-07 2.10E-07 2.90E-07 2.70E-07 2.24E-07 8.63E-08 2.70E-07 2.10E-07 3.02E-07 5.08E-07 6.24E-07 5.58E-07 2.15E-07 6.24E-07 3.13E-07 4.33E-07 1.29E-06 1.70E-06 1.57E-06 6.05E-07 1.70E-06 ]; figure(7), clf criterion = 1e-7; c1=contour(f, levels, with,[criterion criterion]) hold on c2=contour(f, levels, without,[criterion criterion]) title(['criterion= ' num2str(criterion)]) set(gca,'Xscale','log')