summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2016-03-28 17:16:33 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2016-03-28 17:16:33 -0400
commita8dbe52e26b350fce7e88b6f7c90ee8beb4fa054 (patch)
tree194a6bd45943500e62292f624f5c552d7239105c
parent0e9f79d6a4b7db517510c1b2832724d5953c3eec (diff)
downloadnoisy_eit_xmds-a8dbe52e26b350fce7e88b6f7c90ee8beb4fa054.tar.gz
noisy_eit_xmds-a8dbe52e26b350fce7e88b6f7c90ee8beb4fa054.zip
new changes for noise traces
-rw-r--r--delta_2photon_noise.xmds10
-rw-r--r--pp.m13
2 files changed, 18 insertions, 5 deletions
diff --git a/delta_2photon_noise.xmds b/delta_2photon_noise.xmds
index 7723082..2cb3a52 100644
--- a/delta_2photon_noise.xmds
+++ b/delta_2photon_noise.xmds
@@ -21,7 +21,7 @@
const complex gbc=.001;
const complex Split = .1;
const double tset=5;
- const double alpha = .0;
+ const double alpha = 0.001;
complex Gab, Gca, Gcb, GCB;
complex Ep,Ed,Epc,Edc;
@@ -41,7 +41,7 @@
</transverse_dimensions>
</geometry>
- <driver name="multi-path" paths="2" />
+ <driver name="multi-path" paths="16" />
<noise_vector name="drivingNoise" kind="wiener" type="real" method="dsfmt" >
<components>Eta</components>
@@ -80,10 +80,10 @@
Epc = conj(Ep);
Edc = conj(Ed);
- Gab=g+i*(delta_d+delta_2photon+noise);
+ Gab=g+i*(delta_d+delta_2photon+1*noise);
Gca=g-i*(delta_d);
- Gcb=gbc+i*(Split + delta_2photon+noise);
- GCB=gbc+i*(-Split + delta_2photon+noise);
+ Gcb=gbc+i*(Split + delta_2photon+0*noise);
+ GCB=gbc+i*(-Split + delta_2photon+0*noise);
rba=conj(rab);
diff --git a/pp.m b/pp.m
index 150cb01..f59f19c 100644
--- a/pp.m
+++ b/pp.m
@@ -1,8 +1,21 @@
delta_2photon_noise
+figure(1)
plot(delta_2photon_1, mean_mor_signal_1(end,:) )
hold all
plot(delta_2photon_1, mean_eit_signal_1(end,:) )
legend( {'mor', 'eit'} )
+figure(2)
+plot(delta_2photon_1, stderr_mor_signal_1(end,:) )
+hold all
+plot(delta_2photon_1, stderr_eit_signal_1(end,:) )
+
+legend( {'mor', 'eit'} )
+
+figure(3)
+plot(delta_2photon_1, stderr_eit_signal_1(end,:)./stderr_mor_signal_1(end,:) )
+hold all
+
+legend( {'noise eit/mor'} )