diff options
Diffstat (limited to 'compass_lin_extrema_vs_theta_output_results.m')
| -rw-r--r-- | compass_lin_extrema_vs_theta_output_results.m | 25 |
1 files changed, 15 insertions, 10 deletions
diff --git a/compass_lin_extrema_vs_theta_output_results.m b/compass_lin_extrema_vs_theta_output_results.m index 89ad2f6..44b20df 100644 --- a/compass_lin_extrema_vs_theta_output_results.m +++ b/compass_lin_extrema_vs_theta_output_results.m @@ -20,6 +20,14 @@ background_transmission=repmat( background_vector , N_detunings, 1); transmission_matrix=-background_transmission+transmission_matrix; +% Resonances amplitude +am3 = (transmission_matrix(1,:)); % a_{-3} +am2 = (transmission_matrix(2,:)); % a_{-2} +am1 = (transmission_matrix(3,:)); % a_{-1} +a0 = (transmission_matrix(4,:)); % a0 +ap1 = (transmission_matrix(5,:)); % a_1 +ap2 = (transmission_matrix(6,:)); % a_2 +ap3 = (transmission_matrix(7,:)); % a_3 line_colors= [ ... @@ -71,18 +79,15 @@ 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. -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 = ap3./(a0+ap3); +y_combo = ap2./(a0+ap2); +eps=2e-5; +ind= abs(a0+ap3)<eps; x_combo(ind)=1; -y_combo(ind)=0; +%y_combo(ind)=0; %x_combo = peak_1; %y_combo = peak_2; -ind=thetas<=pi/2; +ind=(thetas<=pi/2); surface([x_combo(ind);x_combo(ind)], ... [y_combo(ind);y_combo(ind)], ... @@ -100,7 +105,7 @@ surface([x_combo(ind);x_combo(ind)], ... %xlabel( labels{j}); %ylabel( labels{k}); xlabel('Sideband ratio: a_3/(a_0+a_3)'); -ylabel('Sideband ratio: a_0/(a_0+a_3)'); +ylabel('Sideband ratio: a_2/(a_0+a_2)'); %clabel('\theta'); cbh=colorbar; colormap('jet'); |