diff options
| -rw-r--r-- | liouville.m | 25 | ||||
| -rw-r--r-- | output_results.m | 21 |
2 files changed, 26 insertions, 20 deletions
diff --git a/liouville.m b/liouville.m index ab583b8..c0c63db 100644 --- a/liouville.m +++ b/liouville.m @@ -64,9 +64,9 @@ for detuning_p_cntr=1:N_detun_steps+1; E_field_drive =[0, 0 , Ed, 0 , Edc, 0, 0 ]; E_field_probe =[0, Ep, 0 , Epc, 0 , 0, 0 ]; E_field_zero =[0, 0 , 0 , 0 , 0 , 0, 0 ]; - E_field.linear = E_field_zero + 0.00*E_field_probe + 0.00*E_field_drive; - E_field.right = E_field_zero + 1.00*E_field_probe + 0.00*E_field_drive; - E_field.left = E_field_zero + 0.00*E_field_probe + 1.00*E_field_drive; + E_field.linear = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; + E_field.right = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; + E_field.left = E_field_zero + 1.00*E_field_probe + 1.00*E_field_drive; freq_index=freq2index(wp,modulation_freq); atom_field_problem.E_field = E_field; @@ -88,24 +88,9 @@ fflush (stderr); %kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); +save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; -figure(1); - hold off; - plot(detuning_freq, imag(xi_linear), '-1;linear;'); - hold on; - plot(detuning_freq, imag(xi_left), '-2;left;'); - plot(detuning_freq, imag(xi_right), '-3;right;'); - title("probe absorption"); - hold off; -figure(2); - hold off; - plot(detuning_freq, real(xi_linear), '-1;linear;'); - hold on; - plot(detuning_freq, real(xi_left), '-2;left;'); - plot(detuning_freq, real(xi_right), '-3;right;'); - title("probe dispersion"); - hold off; - +output_results; elapsed_time = etime (clock (), t0) diff --git a/output_results.m b/output_results.m new file mode 100644 index 0000000..a743707 --- /dev/null +++ b/output_results.m @@ -0,0 +1,21 @@ +1; + +load 'xi_vs_detuning.mat' ; + +figure(1); + hold off; + plot(detuning_freq, imag(xi_linear), '-1;linear;'); + hold on; + plot(detuning_freq, imag(xi_left), '-2;left;'); + plot(detuning_freq, imag(xi_right), '-3;right;'); + title("probe absorption"); + hold off; +figure(2); + hold off; + plot(detuning_freq, real(xi_linear), '-1;linear;'); + hold on; + plot(detuning_freq, real(xi_left), '-2;left;'); + plot(detuning_freq, real(xi_right), '-3;right;'); + title("probe dispersion"); + hold off; + |