diff options
| -rw-r--r-- | squeezing_filter/susceptibility_problem.m | 2 | ||||
| -rwxr-xr-x | squeezing_filter/task_solver.m | 7 | ||||
| -rwxr-xr-x | squeezing_filter/tasks_builder.m | 6 | ||||
| -rw-r--r-- | squeezing_filter/tasks_pp.m | 23 |
4 files changed, 27 insertions, 11 deletions
diff --git a/squeezing_filter/susceptibility_problem.m b/squeezing_filter/susceptibility_problem.m index 56808e9..140204c 100644 --- a/squeezing_filter/susceptibility_problem.m +++ b/squeezing_filter/susceptibility_problem.m @@ -1,4 +1,4 @@ -function [xi_linear, xi_left, xi_right]=susceptibility_problem(detuning_freq, Ep, psi_el, B_field, theta, phi) +function [xi_linear, xi_left, xi_right, E_field_pos_freq, light_positive_freq]=susceptibility_problem(detuning_freq, Ep, psi_el, B_field, theta, phi) % calculates transmission if light polarizations vs B field in the cell % for given laser probe and B fields array % Probe field defined by field strength (Ep) and ellipticity angle (pse_el) diff --git a/squeezing_filter/task_solver.m b/squeezing_filter/task_solver.m index 0ba57ea..c044cf6 100755 --- a/squeezing_filter/task_solver.m +++ b/squeezing_filter/task_solver.m @@ -11,9 +11,10 @@ output_fname=arg_list{2}; load(params_file); -[xi_linear, xi_left, xi_right]=susceptibility_problem(detuning_freq, Ep, psi_el, B_field, theta, phi) +[xi_linear, xi_left, xi_right, E_field_pos_freq, light_positive_freq]=susceptibility_problem(detuning_freq, Ep, psi_el, B_field, theta, phi); -save(output_fname, 'detuning_freq', 'Ep', 'psi_el', 'B_field', 'theta', 'phi', ... - 'xi_linear', 'xi_left', 'xi_right' ... +save(output_fname, 'detuning_freq', 'Ep', 'psi_el', 'B_field', 'theta', 'phi' ... + , 'xi_linear', 'xi_left', 'xi_right' ... + , 'E_field_pos_freq', 'light_positive_freq' ... ); diff --git a/squeezing_filter/tasks_builder.m b/squeezing_filter/tasks_builder.m index 8af439e..1131113 100755 --- a/squeezing_filter/tasks_builder.m +++ b/squeezing_filter/tasks_builder.m @@ -15,7 +15,7 @@ gmg=.7; % gyro magnetic ration for ground level % phi is angle between linear polarization and axis x %phi=pi/4; -phi=pi/2; +phi=0; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') theta=0; % psi_el is the ellipticity parameter (phase difference between left and right polarization) @@ -24,8 +24,8 @@ psi_el=0*3/180*pi; Ep=sqrt(10.0); % we are going to sweep B so parts related to sweep paramer(s) -zeeman_splitting=+0.1; -Nsteps=100; +zeeman_splitting=+6; +Nsteps=1001; B_fields=linspace(-zeeman_splitting/gmg, zeeman_splitting/gmg, Nsteps); for i=1:Nsteps diff --git a/squeezing_filter/tasks_pp.m b/squeezing_filter/tasks_pp.m index edc952a..82bcdac 100644 --- a/squeezing_filter/tasks_pp.m +++ b/squeezing_filter/tasks_pp.m @@ -15,8 +15,18 @@ for i=1:Nsteps xi_linear(i)=d.xi_linear; xi_left(i)=d.xi_left; xi_right(i)=d.xi_right; + Ep.linear(i)=d.E_field_pos_freq.linear; + Ep.left(i)=d.E_field_pos_freq.left; + Ep.right(i)=d.E_field_pos_freq.right; endfor +Ep_out.linear=(1-xi_linear).*Ep.linear; +Ep_out.left=(1-xi_left).*Ep.left; +Ep_out.right=(1-xi_right).*Ep.right; + +Ep_out.x= (Ep_out.left + Ep_out.right)/sqrt(2); +Ep_out.y= 1i*(Ep_out.left - Ep_out.right)/sqrt(2); + xi_x=(xi_right+xi_left)/sqrt(2); xi_y=(1i*xi_right-1i*xi_left)/sqrt(2); figure(1); @@ -40,11 +50,16 @@ legend('left','right'); title('xi-imag'); figure(5); -plot(B_field, abs(xi_x).^2, B_field, abs(xi_y).^2 ); +plot(B_field, abs(Ep_out.x).^2, B_field, abs(Ep_out.y).^2 ); legend('x','y'); -title('absorption'); +title('Signal out'); figure(6); -plot(B_field, abs(xi_left).^2, B_field, abs(xi_right).^2 ); +plot(B_field, abs(Ep_out.left).^2, B_field, abs(Ep_out.right).^2 ); legend('left','right'); -title('absorption'); +title('Signal out'); + +figure(7); +plot(B_field, abs(Ep_out.y).^2 ); +legend('y'); +title('Signal out'); |