diff options
| -rw-r--r-- | delta_2photon_noise.xmds | 10 | ||||
| -rw-r--r-- | pp.m | 13 |
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); @@ -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'} ) |