diff options
| -rw-r--r-- | compass_circ.m | 6 | ||||
| -rw-r--r-- | compass_circ_output_results.m | 14 |
2 files changed, 15 insertions, 5 deletions
diff --git a/compass_circ.m b/compass_circ.m index 106b895..d80d179 100644 --- a/compass_circ.m +++ b/compass_circ.m @@ -133,12 +133,12 @@ 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); % strangely parcell is slower than cellfun 20 seconds vs 29 -%total_relative_transmission_vs_phi=parcellfun(2, @total_relative_transmission, problems_cell_array); -total_relative_transmission_vs_phi=cellfun(@total_relative_transmission, problems_cell_array); +%total_relative_transmission_vs_theta=parcellfun(2, @total_relative_transmission, problems_cell_array); +total_relative_transmission_vs_theta=cellfun(@total_relative_transmission, problems_cell_array); %save 'xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; problem_cntr--; -save '/tmp/total_relative_transmission_vs_phi.mat' detuning_freq total_relative_transmission_vs_phi thetas problem_cntr; +save '/tmp/total_relative_transmission_vs_theta.mat' detuning_freq total_relative_transmission_vs_theta thetas problem_cntr; compass_circ_output_results; diff --git a/compass_circ_output_results.m b/compass_circ_output_results.m index 627ead7..b2c7a04 100644 --- a/compass_circ_output_results.m +++ b/compass_circ_output_results.m @@ -1,7 +1,7 @@ 1; -load '/tmp/total_relative_transmission_vs_phi.mat' ; +load '/tmp/total_relative_transmission_vs_theta.mat' ; % let's create sideband transmission vs angle vectors % 1st of all we need to create matrix instead of a vector @@ -11,7 +11,7 @@ load '/tmp/total_relative_transmission_vs_phi.mat' ; N_detunings=length(detuning_freq); N_angles=length(thetas); -transmission_matrix=reshape(total_relative_transmission_vs_phi, N_detunings, N_angles); +transmission_matrix=reshape(total_relative_transmission_vs_theta, N_detunings, N_angles); % the last sideband is not in two-photon resonance % we use it as a reference for background transmission @@ -54,4 +54,14 @@ legend(labels); hold off; +figure(2) +j=4; k= 3; +j=2; k= 7; +% plotting parametric line of a sideband amplitude (j) vs another one (k) +plot(zoom_factor*(transmission_matrix(j,:)), zoom_factor*(transmission_matrix(k,:))) +xlabel( labels{j}); +ylabel( labels{k}); +title('One sideband amplitude vs another for different angles theta, between B-field and light propagation direction'); + + % vim: ts=2:sw=2:fdm=indent |