diff options
Diffstat (limited to 'psr')
| -rw-r--r-- | psr/field_description.m | 2 | ||||
| -rw-r--r-- | psr/output_psr_results.m | 68 | ||||
| -rw-r--r-- | psr/psr.m | 3 |
3 files changed, 71 insertions, 2 deletions
diff --git a/psr/field_description.m b/psr/field_description.m index e9fa227..51131da 100644 --- a/psr/field_description.m +++ b/psr/field_description.m @@ -1,7 +1,7 @@ 1; %EM field definition -Ed=0.002; %drive +Ed=0.2; %drive Edc=conj(Ed); Ep=0.8*Ed; %probe Epc=conj(Ep); diff --git a/psr/output_psr_results.m b/psr/output_psr_results.m new file mode 100644 index 0000000..c6120fc --- /dev/null +++ b/psr/output_psr_results.m @@ -0,0 +1,68 @@ +1; + + +load '/tmp/xi_vs_detuning.mat' ; + +Er=(1+I*xi_right)*E_field_pos_freq.right; +El=(1+I*xi_left) *E_field_pos_freq.left; + +Ex=(Er+El)/sqrt(2); +Ey=I*(Er-El)/sqrt(2); + +%extra rotation to compensate rotation due to ellipticity +% actually no need for it since x-polarization shifts by positive phase +% and y-pol by negative phase +%el_rot=0*psi_el; +%Ex=cos(el_rot)*Ex-sin(el_rot)*Ey; +%Ey=sin(el_rot)*Ex+cos(el_rot)*Ey; + +Ipos=(abs(Ey).^2)/2; +Ineg=(abs(Ex).^2)/2; + +figure(1); +hold off; +plot(detuning_freq, real(xi_left-xi_right), '-'); +title("differential real xi"); +xlabel("two photon detuning"); + +figure(2); +hold off; +plot(detuning_freq, imag(xi_left-xi_right), '-'); +title("differential imag xi"); +xlabel("two photon detuning"); + +figure(3); +hold off; +plot(detuning_freq, imag(xi_left), '-', detuning_freq, imag(xi_right), '-'); +title("imag xi"); +xlabel("two photon detuning"); + +figure(4); +hold off; +plot(detuning_freq, real(xi_left), '-', detuning_freq, real(xi_right), '-'); +title("real xi"); +xlabel("two photon detuning"); + +figure(5); +hold off; +plot(detuning_freq, (Ipos-Ineg), '-'); +title("BPD signal xi"); +xlabel("two photon detuning"); + +%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; + @@ -29,8 +29,9 @@ field_description; detuning_p=0; N_detun_steps=100; %detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting -detuning_p_min=-40.0; +detuning_p_min=-200.0; detuning_p_max=-detuning_p_min; +detuning_p_max=1000; detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); |