From 47a56000e7c6d248087da92412aedbf60da6cc1b Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 18 Nov 2011 00:15:33 -0500 Subject: small modifications --- faraday_and_psr/builder_tasks.m | 3 +-- faraday_and_psr/pp_tasks.m | 59 ++++++++++++++++++++++++----------------- 2 files changed, 36 insertions(+), 26 deletions(-) (limited to 'faraday_and_psr') diff --git a/faraday_and_psr/builder_tasks.m b/faraday_and_psr/builder_tasks.m index 78eea9e..3d81057 100755 --- a/faraday_and_psr/builder_tasks.m +++ b/faraday_and_psr/builder_tasks.m @@ -14,14 +14,13 @@ gmg=.7; % gyro magnetic ration for ground level %[psr_rad]=psr_vs_detuning(Ep, psi_el, B_field, theta, phi) % phi is angle between linear polarization and axis x -%phi=pi/4; 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) psi_el=0*3/180*pi; -Ep=sqrt(10.1); +Ep=sqrt(100); % we are going to sweep B so parts related to sweep paramer(s) zeeman_splitting=+0.5; diff --git a/faraday_and_psr/pp_tasks.m b/faraday_and_psr/pp_tasks.m index 82bcdac..a3919ba 100644 --- a/faraday_and_psr/pp_tasks.m +++ b/faraday_and_psr/pp_tasks.m @@ -29,37 +29,48 @@ 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); -plot(B_field, real(xi_x), B_field, real(xi_y)); -legend('x','y'); -title('xi-real'); -figure(2); -plot(B_field, imag(xi_x), B_field, imag(xi_y)); -legend('x','y'); -title('xi-imag'); +%figure(1); +%plot(B_field, real(xi_x), B_field, real(xi_y)); +%legend('x','y'); +%title('xi-real'); -figure(3); -plot(B_field, real(xi_left), B_field, real(xi_right)); -legend('left','right'); -title('xi-real'); +%figure(2); +%plot(B_field, imag(xi_x), B_field, imag(xi_y)); +%legend('x','y'); +%title('xi-imag'); -figure(4); -plot(B_field, imag(xi_left), B_field, imag(xi_right)); -legend('left','right'); -title('xi-imag'); +%figure(3); +%plot(B_field, real(xi_left), B_field, real(xi_right)); +%legend('left','right'); +%title('xi-real'); -figure(5); -plot(B_field, abs(Ep_out.x).^2, B_field, abs(Ep_out.y).^2 ); -legend('x','y'); -title('Signal out'); +%figure(4); +%plot(B_field, imag(xi_left), B_field, imag(xi_right)); +%legend('left','right'); +%title('xi-imag'); -figure(6); -plot(B_field, abs(Ep_out.left).^2, B_field, abs(Ep_out.right).^2 ); -legend('left','right'); -title('Signal out'); +%figure(5); +%plot(B_field, abs(Ep_out.x).^2, B_field, abs(Ep_out.y).^2 ); +%legend('x','y'); +%title('Signal out'); + +%figure(6); +%plot(B_field, abs(Ep_out.left).^2, B_field, abs(Ep_out.right).^2 ); +%legend('left','right'); +%title('Signal out'); figure(7); +plot(B_field, abs(Ep_out.x).^2 ); +legend('x'); +title('Signal out'); + +figure(8); plot(B_field, abs(Ep_out.y).^2 ); legend('y'); title('Signal out'); + +figure(9); +plot(B_field, abs(Ep_out.x).^2-abs(Ep_out.y).^2 ); +legend('x-y'); +title('Differential x-y out'); -- cgit v1.2.3