diff options
| -rw-r--r-- | compass_lin_extrema_vs_theta_output_results.m | 48 |
1 files changed, 41 insertions, 7 deletions
diff --git a/compass_lin_extrema_vs_theta_output_results.m b/compass_lin_extrema_vs_theta_output_results.m index cfb101b..eb151b8 100644 --- a/compass_lin_extrema_vs_theta_output_results.m +++ b/compass_lin_extrema_vs_theta_output_results.m @@ -33,7 +33,9 @@ line_colors= [ ... ]; figure(1); +if (exist('keep','var') && keep) clf(); +end hold off; labels={}; for i=1:N_detunings-1 @@ -58,21 +60,39 @@ print('compass_lin_sidebands_vs_theta.png') hold off; figure(2) +if (exist('keep','var') && keep) clf(); +end +% resonance -3 -2 -1 0 1 2 3 +% index 1 2 3 4 5 6 7 j=4; k= 3; -j=2; k= 7; +j=6; k= 7; +j=4; k= 7; % plotting parametric line of a sideband amplitude (j) vs another one (k) z = zeros(1,N_angles); col = thetas; % This is the color, vary with x in this case. -surface([zoom_factor*(transmission_matrix(j,:));zoom_factor*(transmission_matrix(j,:))], ... - [zoom_factor*(transmission_matrix(k,:));zoom_factor*(transmission_matrix(k,:))], ... +peak_1 = (transmission_matrix(j,:)); +peak_2 = (transmission_matrix(k,:)); +peak_3 = (transmission_matrix(6,:)); +x_combo = peak_2./(peak_1+peak_2); +y_combo = peak_1./(peak_1+peak_2); +eps=1e-5; +ind= abs(peak_1+peak_2)<eps; +x_combo(ind)=1; +y_combo(ind)=0; +%x_combo = peak_1; +%y_combo = peak_2; +surface([x_combo;x_combo], ... + [y_combo;y_combo], ... [z;z],[col;col],... 'facecol','no',... 'edgecol','interp',... 'linew',2); %plot(zoom_factor*(transmission_matrix(j,:)), zoom_factor*(transmission_matrix(k,:))) -xlabel( labels{j}); -ylabel( labels{k}); +%xlabel( labels{j}); +%ylabel( labels{k}); +xlabel('Sideband ratio: a_3/(a_0+a_3)'); +ylabel('Sideband ratio: a_0/(a_0+a_3)'); %clabel('\theta'); colorbar; set(gca, 'ZTick', [0,pi/4, pi/2, 3*pi/4, pi]) @@ -80,8 +100,22 @@ set(gca, 'ZTickLabel', {'0', '\pi/4', '\pi/2', '3\pi/4', '\pi'}) %colorbar('Ticks',[0, pi/2, pi], 'TicksLabels', ['0', '\pi/2', '\pi']) %ctickslabel(['0', '\pi/2', '\pi']) -title('One sideband amplitude vs another for different angles theta'); -print('compass_lin_sidebands_combo_vs_theta.png') +title('Phase diagram for different angles theta'); +print('compass_lin_sidebands_phase_diagram_vs_theta.png') + +figure(3) +if (exist('keep','var') && keep) +clf(); +end + +plot(thetas, (x_combo)) +set(gca, 'XTick', [0,pi/4, pi/2, 3*pi/4, pi]) +set(gca, 'XTickLabel', {'0', '\pi/4', '\pi/2', '3\pi/4', '\pi'}) +title("Sideband ratio: a_3/(a_0+a_3)"); +xlabel('Angle \theta, between B-field and light propagation direction'); +ylabel('Ratio'); +ylim([0:1]); +print('compass_lin_sidebands_ratio_vs_theta.png') % vim: ts=2:sw=2:fdm=indent |