WIP, adding test cases.

ternary_frap.m

@@ -165,7 +165,7 @@ surf(X,Y,log(Z));
 set(hAx,'xscale','log');
 set(hAx,'yscale','log');
 h = colorbar('northoutside', 'XTickLabel',{'10^{-6}','10^{-4}','10^{-2}','1','10^2','10^4','10^6','10^8'}, ...
-            'YTick',[-6:2:8]);
+            'YTick',-6:2:8);
 ylabel(h, [char(964) '_{in}/' char(964) '_{out}']);
 % set(h,'YTick',[-6:2:8], '10);
 colormap(c)
@@ -180,5 +180,28 @@ nu = 10^-6;
 chi = 7/3;
 b = nu^(1/3)*sqrt(chi/(chi-2));
 e = sqrt(3/8*(chi-2));
-T = Ternary_model(2, 'FRAP', [-1, b, 0.5, e, 1, 0.5], linspace(0.0, 400, 5000));
-T.solve_tern_frap();
\ No newline at end of file
+t = linspace(0.0, 40, 500);
+T = Ternary_model(2, 'FRAP', [-1, b, 0.5, e, 1, 0.5], t);
+T.solve_tern_frap();
+%%
+x_q = linspace(0, 0.9, 100);
+v_q = interp1(T.x, T.sol', x_q);
+v_q = v_q/max(v_q(:));
+% plot(sum(v_q, 1))
+x_fit = [10, 2];
+model = 'simple_drop';
+f_min = @(x_fit) intensity_fit(v_q', t', x_fit, model, sum(v_q, 1));
+opt = optimset('MaxFunEvals', 2000, 'PlotFcns',@optimplotfval);
+for i = 1:1
+    x_fit = fminsearch(f_min, x_fit, opt);
+end
+%%
+x_fit = [8, 2];
+intensity_fit(v_q', t', x_fit, model, sum(v_q, 1))
+[~, v_fit] = to_minimize(v_q', t', x_fit, model);
+figure; hold on;
+plot(v_q);
+plot(v_fit')
+figure; hold on;
+plot(sum(v_q));
+plot(sum(v_fit, 2));
\ No newline at end of file