From 2890c93644aada7b06b467be8cc548f6a9004769 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Tue, 22 Nov 2011 23:05:51 -0500
Subject: added final plot for circular polarizations

---
 .../circular_perpendicular_results/Makefile        | 31 +++++++++++++
 .../circular_perpendicular_results/grand_output.m  | 54 ++++++++++++++++++++++
 2 files changed, 85 insertions(+)
 create mode 100644 squeezing_filter/circular_perpendicular_results/Makefile
 create mode 100644 squeezing_filter/circular_perpendicular_results/grand_output.m

diff --git a/squeezing_filter/circular_perpendicular_results/Makefile b/squeezing_filter/circular_perpendicular_results/Makefile
new file mode 100644
index 0000000..144b263
--- /dev/null
+++ b/squeezing_filter/circular_perpendicular_results/Makefile
@@ -0,0 +1,31 @@
+###  -*- make -*-
+### This makefile can be used to build and run the XMDS examples
+
+
+
+eps_targets  = $(wildcard *.eps)
+pdf_targets  = $(eps_targets:%.eps=%.pdf)	
+png_targets  = $(pdf_targets:%.pdf=%.png)
+
+pdf: $(pdf_targets) 
+
+$(pdf_targets): %.pdf : %.eps
+	cat $< | ps2eps -B > __tt.eps
+	epspdf __tt.eps $@
+	rm -f __tt.eps
+	#ps2eps -B $< | epspdf $< $@
+
+png: pdf $(png_targets) 
+
+$(png_targets): %.png : %.pdf
+	convert -density 300 $< $@
+
+clean: 
+	rm -f $(pdf_targets)
+	rm -f $(png_targets)
+
+real_clean: clean
+	rm -f $(eps_targets)
+
+.PRECIOUS: %.m
+.PHONY: all clean
diff --git a/squeezing_filter/circular_perpendicular_results/grand_output.m b/squeezing_filter/circular_perpendicular_results/grand_output.m
new file mode 100644
index 0000000..b8ac424
--- /dev/null
+++ b/squeezing_filter/circular_perpendicular_results/grand_output.m
@@ -0,0 +1,54 @@
+
+fnames=glob('./*.mat');
+
+Nsteps=length(fnames)
+
+line_colors= [ ...
+	[ 0,   0,   1]; ...
+	[ 1,   0,   0]; ...
+	[ 0,   1,   0]; ...
+	[ 0,   0,   0]; ...
+	[ 0, 0.8,   0]; ...
+	[ 1,   0,   1]; ...
+	[ 0,   0,  .6] ...
+	];
+
+figure(1); hold off;
+legend_str={};
+
+% read the information from resulting files
+for i=1:Nsteps
+	d=load(fnames{i});
+	detunings=d.detunings;
+	B_field=d.B_field(1);
+	xi_left=d.xi_left;
+	xi_right=d.xi_right;
+	xi_x=d.xi_x;
+	xi_y=d.xi_y;
+	Ep  = d.Ep;
+	Ed=d.Ed;
+	Ep_out.x  = d.Ep_out.x;
+	Ep_out.y  = d.Ep_out.y;
+	Ep_out.left  = d.Ep_out.left;
+	Ep_out.right  = d.Ep_out.right;
+
+	B_str=num2str(B_field(1),"%g");
+	B_str=strcat('B= ',B_str, ' G');
+	legend_str{i}=B_str;
+
+	figure(1);
+	plot(detunings, abs(Ep_out.right).^2, 'color', line_colors(i,:) ); hold on
+	legend(legend_str);
+
+endfor
+
+
+figure(1); 
+xlabel('Detuning (MHz)');
+ylabel('Transmission (Arb.Units)');
+title('Signal out right-polarization');
+legend(legend_str);
+legend('location', 'south');
+hold off;
+
+print('-color','-depsc2', '-tight','-S200,120', 'right-polarization_output_vs_detuning.eps')
-- 
cgit v1.2.3