view multithreshold 1.46/compareMicrophonic.m @ 25:d2c4c07df02c

path problem solved

% these data are generated by the testBM program


f=[250   500  1000  2000  4000  8000];
levels=-10:10:90;

with=[
3.00036e-010	2.9476e-010	3.67369e-010	2.69424e-010	2.61262e-010	7.18365e-011
9.48796e-010	9.32111e-010	1.16172e-009	8.51994e-010	8.26184e-010	2.27167e-010
3.00036e-009	2.9476e-009	3.67369e-009	2.69424e-009	2.61262e-009	7.18365e-010
9.48796e-009	9.32126e-009	1.16172e-008	8.51994e-009	8.26184e-009	2.27167e-009
2.98964e-008	3.30801e-008	4.3799e-008	2.85983e-008	2.12715e-008	7.18365e-009
8.08142e-008	1.0005e-007	1.0854e-007	8.23884e-008	6.41103e-008	1.57936e-008
1.4489e-007	1.70294e-007	1.67989e-007	1.37031e-007	9.77912e-008	3.74565e-008
1.9871e-007	2.56245e-007	2.50429e-007	1.82704e-007	1.2266e-007	5.73054e-008
2.80266e-007	4.15464e-007	4.42073e-007	3.65116e-007	3.24109e-007	6.7619e-008
4.49308e-007	8.26643e-007	1.0501e-006	9.22625e-007	9.46576e-007	2.13766e-007
8.73737e-007	2.1498e-006	2.9006e-006	2.61691e-006	2.90337e-006	6.36048e-007
];

without=[
1.85179e-009	1.84842e-009	2.29365e-009	1.68885e-009	1.64802e-009	4.53194e-010
5.85587e-009	5.84521e-009	7.25317e-009	5.34061e-009	5.2115e-009	1.43313e-009
1.85179e-008	1.84842e-008	2.39525e-008	1.67348e-008	1.41015e-008	4.53194e-009
5.70059e-008	6.97074e-008	8.049e-008	5.93827e-008	4.75289e-008	1.43313e-008
1.21363e-007	1.42557e-007	1.42791e-007	1.14796e-007	8.59727e-008	2.31062e-008
1.7631e-007	2.1607e-007	2.05558e-007	1.61584e-007	1.09552e-007	5.05732e-008
2.41242e-007	3.44485e-007	3.25688e-007	2.57527e-007	2.17035e-007	6.31655e-008
3.64591e-007	5.90911e-007	7.23077e-007	6.17822e-007	6.09472e-007	1.38217e-007
6.52753e-007	1.44434e-006	1.89854e-006	1.69908e-006	1.84634e-006	4.09033e-007
1.43691e-006	4.03513e-006	5.54045e-006	5.04855e-006	5.73398e-006	1.24307e-006
3.75649e-006	1.20961e-005	1.69372e-005	1.55783e-005	1.7989e-005	3.87293e-006
];


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')