From 43cf01ff8494f8381c276cb3861682b8f90cd810 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Wed, 11 Nov 2020 15:34:38 +0100
Subject: [PATCH] Mean cross jump length comparison with Stefano works.
---
prob_laplace.m | 8 +++-----
run_jump_lengths.m | 3 +++
2 files changed, 6 insertions(+), 5 deletions(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index 79cb44f..2d1b5a3 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -185,9 +185,8 @@ for i = 1:100%length(T_mov.t)
% print([num2str(i),'.png'],'-dpng')
shg; pause();
end
-%%
-run_jump_lengths(-5, -0.15432893, 7/3, 1, 4, '73_0154_10', 10^-6, 4, 0.5);
-%% chi 7/3
+
+%% Mean cross jump length chi 7/3
run_jump_lengths(-5, 1, 7/3, 1, 4, '73_1_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, 1, 7/3, 1, 4, '73_1_10', 10^-6, 4, 0.5);
run_jump_lengths(-5, 0, 7/3, 1, 4, '73_6_01', 10^-6, 3, 0.5);
@@ -202,8 +201,6 @@ run_jump_lengths(-5, -0.15432893, 7/3, 1, 4, '73_25_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, -0.15432893, 7/3, 1, 4, '73_25_10', 10^-6, 4, 0.5);
run_jump_lengths(-5, -0.16211677, 7/3, 1, 4, '73_30_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, -0.16211677, 7/3, 1, 4, '73_30_10', 10^-6, 4, 0.5);
-
-replace(num2str(round(chi,2)),'.', '')
%% chi 7.7/3
run_jump_lengths(-5, 1, 7.7/3, 1, 4, '773_1_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, 1, 7.7/3, 1, 4, '773_1_10', 10^-6, 4, 0.5);
@@ -219,6 +216,7 @@ run_jump_lengths(-5, 0, 7.7/3, 1, 4, '773_25_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, 0, 7.7/3, 1, 4, '773_25_10', 10^-6, 4, 0.5);
run_jump_lengths(-5, -0.00756985349, 7.7/3, 1, 4, '773_30_01', 10^-6, 3, 0.5);
run_jump_lengths(-5, -0.00756985349, 7.7/3, 1, 4, '773_30_10', 10^-6, 4, 0.5);
+
%% %%%%%%%%%%%%%%%%% STEADY STATE JUMP LENGTH DISTRIBUTION %%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 0, 1, 10, 7, 0,...
diff --git a/run_jump_lengths.m b/run_jump_lengths.m
index c1fa91f..4a55a2d 100644
--- a/run_jump_lengths.m
+++ b/run_jump_lengths.m
@@ -30,4 +30,7 @@ toc
% should sum to one
N = sum(p)/length(p)*l_max;
m = sum(ls.*p/N)/length(p)*l_max;
+T{1}.calc_real_params();
+D_ratio = T{1}.real_params.D_out/T{1}.real_params.D_in;
save(name);
+disp([D_ratio, m]);
\ No newline at end of file
--
GitLab