diff options
Diffstat (limited to 'compass_lin_output_results.m')
| -rw-r--r-- | compass_lin_output_results.m | 57 |
1 files changed, 57 insertions, 0 deletions
diff --git a/compass_lin_output_results.m b/compass_lin_output_results.m new file mode 100644 index 0000000..dd02f41 --- /dev/null +++ b/compass_lin_output_results.m @@ -0,0 +1,57 @@ +1; + + +load '/tmp/total_relative_transmission_vs_phi.mat' ; + +% let's create sideband transmission vs angle vectors +% 1st of all we need to create matrix instead of a vector +% with rows corresponding to absorption for each sideband +% and columns to each phi angle + +N_detunings=length(detuning_freq); +N_angles=length(phis); + +transmission_matrix=reshape(total_relative_transmission_vs_phi, N_detunings, N_angles); + +% the last sideband is not in two-photon resonance +% we use it as a reference for background transmission +background_vector=transmission_matrix(N_detunings,:); +background_transmission=repmat( background_vector , N_detunings, 1); + +transmission_matrix=-background_transmission+transmission_matrix; + + + +line_colors= [ ... + [ 0, 0, 1]; ... + [ 1, 0, 0]; ... + [ 0, 1, 0]; ... + [ 0, 0, 0]; ... + [ 0, 0.8, 0]; ... + [ 1, 0, 1]; ... + [ 0, 0, .6] ... + ]; + +figure(1); +clf(); +hold off; +labels={}; +for i=1:N_detunings-1 + %we will skip the very last row since it the reference transmission + + zoom_factor=1; + %plot_style=strcat("-", num2str(i)); + %plot( phis, zoom_factor*(transmission_matrix(i,:)), plot_style); + labels = {labels{:}, strcat("Sideband ", num2str(i-4))}; + line( phis, zoom_factor*(transmission_matrix(i,:)), "color", line_colors(i,:) ); + hold on; +endfor + +title("relative transmission"); +xlabel("angle"); +ylabel("amplitude"); +legend(labels); + +hold off; + +% vim: ts=2:sw=2:fdm=indent |