modifications to fit filter problem

2 files changed, 12 insertions, 39 deletions

diff --git a/squeezing_filter/builder_tasks.m b/squeezing_filter/builder_tasks.m
index 970f5e5..3f06913 100755
--- a/squeezing_filter/builder_tasks.m
+++ b/squeezing_filter/builder_tasks.m
@@ -23,7 +23,7 @@ gmg=.7; % gyro magnetic ration for ground level
 zeeman_splitting=+0.5;
 B_field=zeeman_splitting/gmg;
 
-Nsteps=101;
+Nsteps=100;
 detunings=linspace(detuning_freq_start, detuning_freq_stop, Nsteps);
 
 for i=1:Nsteps
diff --git a/squeezing_filter/pp_tasks.m b/squeezing_filter/pp_tasks.m
index a3919ba..45cc0dd 100644
--- a/squeezing_filter/pp_tasks.m
+++ b/squeezing_filter/pp_tasks.m
@@ -11,66 +11,39 @@ xi_right=zeros(1,Nsteps);
 % read the information from resulting files
 for i=1:Nsteps
 	d=load(fnames{i});
+	detunings(i)=d.detuning_freq;
 	B_field(i)=d.B_field;
 	xi_linear(i)=d.xi_linear;
 	xi_left(i)=d.xi_left;
 	xi_right(i)=d.xi_right;
-	Ep.linear(i)=d.E_field_pos_freq.linear;
-	Ep.left(i)=d.E_field_pos_freq.left;
-	Ep.right(i)=d.E_field_pos_freq.right;
+	Ep.left  = d.E_field_pos_freq.left(2);
+	Ep.right = d.E_field_pos_freq.right(2);
+	Ed=d.Ed;
 endfor
 
-Ep_out.linear=(1-xi_linear).*Ep.linear;
+xi_x=(xi_right+xi_left)/sqrt(2);
+xi_y=(1i*xi_right-1i*xi_left)/sqrt(2);
+
+%Ep_out.linear=(1-xi_linear).*Ep.linear;
 Ep_out.left=(1-xi_left).*Ep.left;
 Ep_out.right=(1-xi_right).*Ep.right;
 
 Ep_out.x= (Ep_out.left + Ep_out.right)/sqrt(2);
 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(3);
-%plot(B_field, real(xi_left), B_field, real(xi_right));
-%legend('left','right');
-%title('xi-real');
-
-%figure(4);
-%plot(B_field, imag(xi_left), B_field, imag(xi_right));
-%legend('left','right');
-%title('xi-imag');
-
-%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 );
+plot(detunings, abs(Ep_out.x).^2 );
 legend('x');
 title('Signal out');
 
 figure(8);
-plot(B_field, abs(Ep_out.y).^2 );
+plot(detunings, 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 );
+plot(detunings, abs(Ep_out.x).^2-abs(Ep_out.y).^2 );
 legend('x-y');
 title('Differential x-y out');