1 files changed, 50 insertions, 0 deletions
diff --git a/squeezing_filter/tasks_pp.m b/squeezing_filter/tasks_pp.m
new file mode 100644
index 0000000..edc952a
--- /dev/null
+++ b/squeezing_filter/tasks_pp.m
@@ -0,0 +1,50 @@
+basis_transformation;
+
+fnames=glob('results/*.mat');
+
+Nsteps=length(fnames)
+B_field=zeros(1,Nsteps);
+xi_linear=zeros(1,Nsteps);
+xi_left=zeros(1,Nsteps);
+xi_right=zeros(1,Nsteps);
+
+% read the information from resulting files
+for i=1:Nsteps
+	d=load(fnames{i});
+	B_field(i)=d.B_field;
+	xi_linear(i)=d.xi_linear;
+	xi_left(i)=d.xi_left;
+	xi_right(i)=d.xi_right;
+endfor
+
+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(xi_x).^2, B_field, abs(xi_y).^2 );
+legend('x','y');
+title('absorption');
+
+figure(6);
+plot(B_field, abs(xi_left).^2, B_field, abs(xi_right).^2 );
+legend('left','right');
+title('absorption');