diff Copy_of_multithreshold 1.46/compareMicrophonic.m @ 28:02aa9826efe0

mainly multiThreshold
author Ray Meddis <rmeddis@essex.ac.uk>
date Fri, 01 Jul 2011 12:59:47 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Copy_of_multithreshold 1.46/compareMicrophonic.m	Fri Jul 01 12:59:47 2011 +0100
@@ -0,0 +1,43 @@
+% 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')