Bug fix to EDM SBS

Fitting is now more robust. + Update to plotting in main_0 script
tsipkens · Jan 20, 2024 · 8770102 · 8770102
1 parent 10cbc9e
commit 8770102
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/+pp/edm_sbs.m b/+pp/edm_sbs.m
@@ -103,11 +103,11 @@
     ome = -0.8515; % Omega parameter in sigmoid function
     a = 0.9966;
     sigmoid = @(x) a ./ (1 + ...
-        exp(((log(x(1)) - log(dp_bin)) ./ log(x(2)) - bet) ./ ome));
+        exp(((log(x(1)) - log(dp_bin)) ./ x(2) - bet) ./ ome));
         % x(1) = dpg, x(2) = spg
 
     opts = optimset('Display','off');
-    x0 = [25, 1.05];
+    x0 = [dp_bin(find(Sa < 0.5, 1)), log(1.5)];
     x1 = lsqnonlin(@(x) (sigmoid(x) - Sa) ./ 100, x0, [], [100, 1.1], opts);
     x1 = real(x1);
     x1(x1<0) = NaN;

diff --git a/main_0.m b/main_0.m
@@ -35,7 +35,7 @@
 
 
 figure(4);
-loglog([Aggs_pcm.dp], [Aggs_edm.dp], '.');
+loglog([Aggs_pcm.dp], [Aggs_edm.dp], 'o');
 hold on;
 loglog(xlim,xlim);
 hold off;