5-8 molecules dohne

bnetlab · Jul 10, 2019 · 78f30ac · 78f30ac
1 parent ee577f5
commit 78f30ac
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 2 deletions.
diff --git a/code/Simulation/plot_phase.m b/code/Simulation/plot_phase.m
@@ -0,0 +1,12 @@
+%plot the phase diagram
+
+T = readtable('dataSimulation8.csv')
+res = table2array(T(:,6:end))
+
+[argvalue, argmax] = max(res');
+
+res2 = reshape(argmax,[20,18])
+imagesc([0.5:0.2:4],[0.1:0.2:4], res2)
+xlabel('\tau')
+ylabel('T')
+saveas(gcf,'phase8.png')
diff --git a/code/Simulation/sim_model_prl.py b/code/Simulation/sim_model_prl.py
@@ -90,8 +90,8 @@ def main_prl(N):
     # run for different parameter combination
 
     d=1; v=1; sigma=1
-    t_range = np.arange(0.1,4,0.1)
-    tau_range = np.arange(0.5,4,0.1)
+    t_range = np.arange(0.1,4.01,0.2)
+    tau_range = np.arange(0.5,4.01,0.2)
     Data = np.zeros((len(t_range)*len(tau_range),5+N-1),dtype=float)
     res = Parallel(n_jobs=16)(delayed(sim)(N,d,v,sigma,i,j) for j in tau_range for i in t_range)    
     res = np.array(res)