From eaafb6a468b62d637ed9beb2839cfebc2ba717c7 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 30 Dec 2011 17:46:10 -0500 Subject: files renamed to better reflect the processes --- .../fast_light/Makefile | 54 ---------------- .../fast_light/README | 5 -- .../fast_light/map2dat.m | 33 ---------- .../fast_light/plot_fields_propagation_I2.gp | 15 ----- .../fast_light/pp_I2.m | 71 ---------------------- .../fast_to_slow_switch/Makefile | 4 +- .../fast_to_slow_switch/README | 5 ++ .../fast_to_slow_switch/map2dat.m | 34 ++++++++++- .../plot_fields_propagation_I2.gp | 2 +- .../fast_to_slow_switch/pp_I2.m | 21 ++++--- .../pulse_split/Makefile | 56 +++++++++++++++++ .../pulse_split/map2dat.m | 1 + .../pulse_split/plot_fields_propagation_I2.gp | 15 +++++ .../pulse_split/pp_I2.m | 68 +++++++++++++++++++++ 14 files changed, 192 insertions(+), 192 deletions(-) delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README mode change 120000 => 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m (limited to 'xmds2') diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile deleted file mode 100644 index 76d3f5c..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile +++ /dev/null @@ -1,54 +0,0 @@ -### -*- make -*- -### This makefile can be used to build and run the XMDS examples - - -XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil -M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) -GNUPLOT_FILES = $(wildcard *.gp) - -XSIL2GRAPHICS = xsil2graphics - -# fast light -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 -# slow light EIT -#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 - -all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png - -Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run - $< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt - -%.m: %.xsil - $(XSIL2GRAPHICS) $< - -plot: $(M_FILES) $(GNUPLOT_FILES) - octave pp_I2.m - gnuplot plot_fields_propagation_I2.gp - -clean: - rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf - rm -f $(pdf_targets) - rm -f $(eps_targets) - -real_clean: clean - rm -f $(png_targets) - -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 $< $@ - -.PRECIOUS: %.run %.xsil %.m -.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README deleted file mode 100644 index 1de063d..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README +++ /dev/null @@ -1,5 +0,0 @@ -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 - -1.5 cm seems to give largest advance -2.5 cm gives zero delay -10cm slow light again diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m deleted file mode 100644 index 969b6dc..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m +++ /dev/null @@ -1,33 +0,0 @@ -function map2dat(outfile, x,y,z, xskip, yskip) -% saves 3D data in suitable way to be drawn by gnuplot -% x,y - vectors of x,y values -% z map of z values as used by Octave/Matlab -% xskip, yskip - skip paprameters -% only every, xskip, yskip point will be written - - - -Nx=length(x); -Ny=length(y); -Nxs=Nx/xskip; -Nys=Ny/yskip; -points=zeros(1,3*Nxs*Nys); -%points=[]; -tic; -for i=1:Nxs - for k=1:Nys - %points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)]; - points((i-1)*(Nys-1)*3+3*(k-1)+1) = x(i*xskip); - points((i-1)*(Nys-1)*3+3*(k-1)+2) = y(k*yskip); - points((i-1)*(Nys-1)*3+3*(k-1)+3) = z(k*yskip,i*xskip); - end -end -disp('=== points formation complete ==='); -toc -tic; -%points -fd = fopen(outfile, "wt"); -fprintf (fd, "%g %g %g\n", points); -fclose(fd); -disp('=== points saving complete ==='); -toc; diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp deleted file mode 100644 index 8591d87..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I2.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:1.5][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m deleted file mode 100644 index 0d46ec7..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m +++ /dev/null @@ -1,71 +0,0 @@ -Nlevels_no_dopler_with_z_4wm - -%% field propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -%set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -tmin=-0.4; -tmax= 0.4; -ylim([tmin,tmax],'manual'); -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') - -xskip=1; -yskip=10; -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - - -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') - - - -%% fields before and after the cell -figure(2) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'northeast'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.01; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-.05; -tmax=.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... - ) -legend('location', 'southeast'); -xlim([tmin,tmax],'manual'); -xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps') - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile index 95ebb0a..76d3f5c 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile @@ -9,11 +9,9 @@ GNUPLOT_FILES = $(wildcard *.gp) XSIL2GRAPHICS = xsil2graphics # fast light -# PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 # slow light EIT #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 -#Fast light to Slow light switch -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2 all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README new file mode 100644 index 0000000..1de063d --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README @@ -0,0 +1,5 @@ +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 + +1.5 cm seems to give largest advance +2.5 cm gives zero delay +10cm slow light again diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m deleted file mode 120000 index 14fae30..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m +++ /dev/null @@ -1 +0,0 @@ -../map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m new file mode 100644 index 0000000..969b6dc --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m @@ -0,0 +1,33 @@ +function map2dat(outfile, x,y,z, xskip, yskip) +% saves 3D data in suitable way to be drawn by gnuplot +% x,y - vectors of x,y values +% z map of z values as used by Octave/Matlab +% xskip, yskip - skip paprameters +% only every, xskip, yskip point will be written + + + +Nx=length(x); +Ny=length(y); +Nxs=Nx/xskip; +Nys=Ny/yskip; +points=zeros(1,3*Nxs*Nys); +%points=[]; +tic; +for i=1:Nxs + for k=1:Nys + %points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)]; + points((i-1)*(Nys-1)*3+3*(k-1)+1) = x(i*xskip); + points((i-1)*(Nys-1)*3+3*(k-1)+2) = y(k*yskip); + points((i-1)*(Nys-1)*3+3*(k-1)+3) = z(k*yskip,i*xskip); + end +end +disp('=== points formation complete ==='); +toc +tic; +%points +fd = fopen(outfile, "wt"); +fprintf (fd, "%g %g %g\n", points); +fclose(fd); +disp('=== points saving complete ==='); +toc; diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp index a0cd31f..8591d87 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp @@ -12,4 +12,4 @@ set ylabel "t ({/Symbol m}S)" set zlabel "I_2 (1/S)" set nokey #set view map -splot [0:][-0.2:0.4] 'I2.dat' +splot [0:1.5][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m index ecfdbbb..0d46ec7 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m @@ -4,8 +4,11 @@ Nlevels_no_dopler_with_z_4wm z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS figure(1) -set(gca,'fontsize',20); +%set(gca,'fontsize',20); imagesc(z_1, t_1, I2_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); xlabel('z (cm)') ylabel('t (uS)') zlabel('I_2') @@ -17,7 +20,7 @@ map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); -print('-color','fields_propagation_I2.eps') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') @@ -43,26 +46,26 @@ delay_time=t_1(max_pos_after)-t_1(max_pos_before); printf('Second field delay time = %f uS\n',delay_time); %set(gca,'fontsize',40); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','fields_before_after_cell_I2.eps') +%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps') figure(4) I2_max_in=max(I2_out_1(t_good_indx,1)); I2_max_out=max(I2_out_1(t_good_indx,end)); I2_in_norm=(I2_out_1(:,1))/I2_max_in; I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; +tmin=-.05; +tmax=.05; indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region plot( ... t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... ) +legend('location', 'southeast'); xlim([tmin,tmax],'manual'); xlabel('t (uS)') ylabel('I_2') title('I_2 before and after cell normalized') -legend('location', 'northeast'); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','probe_before_after_cell_I2_normalized.eps') +%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps') diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile new file mode 100644 index 0000000..95ebb0a --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile @@ -0,0 +1,56 @@ +### -*- make -*- +### This makefile can be used to build and run the XMDS examples + + +XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil +M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) +GNUPLOT_FILES = $(wildcard *.gp) + +XSIL2GRAPHICS = xsil2graphics + +# fast light +# PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 +# slow light EIT +#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 +#Fast light to Slow light switch +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2 + +all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png + +Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run + $< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt + +%.m: %.xsil + $(XSIL2GRAPHICS) $< + +plot: $(M_FILES) $(GNUPLOT_FILES) + octave pp_I2.m + gnuplot plot_fields_propagation_I2.gp + +clean: + rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf + rm -f $(pdf_targets) + rm -f $(eps_targets) + +real_clean: clean + rm -f $(png_targets) + +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 $< $@ + +.PRECIOUS: %.run %.xsil %.m +.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m new file mode 120000 index 0000000..14fae30 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m @@ -0,0 +1 @@ +../map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp new file mode 100644 index 0000000..a0cd31f --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp @@ -0,0 +1,15 @@ +set terminal postscript portrait enhanced color solid size 5,3.5 +set output 'fields_propagation_I2.eps' +set dgrid3d 100,100 qnorm 4 +set pm3d map +#set contour +set hidden3d +set palette rgb 10,13,31 negative + + +set xlabel "z (cm)" +set ylabel "t ({/Symbol m}S)" +set zlabel "I_2 (1/S)" +set nokey +#set view map +splot [0:][-0.2:0.4] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m new file mode 100644 index 0000000..ecfdbbb --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m @@ -0,0 +1,68 @@ +Nlevels_no_dopler_with_z_4wm + +%% field propagation +z_1=z_1*100; % z in cm +t_1=t_1*1e6; % time now measured in uS +figure(1) +set(gca,'fontsize',20); +imagesc(z_1, t_1, I2_out_1); colorbar +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_2') +title('I_2') + +xskip=1; +yskip=10; +map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); + + + +print('-color','fields_propagation_I2.eps') + + + +%% fields before and after the cell +figure(2) +%set(gca,'fontsize',30); +plot( ... + t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... + t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... + ) +xlabel('t (uS)') +ylabel('I_2 (1/s)^2') +title('I_2 before and after cell') +legend('location', 'northeast'); + +[b, a]=butter(3, 0.05); +I2_out_after=I2_out_1(:,end); +I2_out_after_filtered=filtfilt(b,a,I2_out_after); +settling_time=0.01; %uS +t_good_indx=t_1> min(t_1 + settling_time); +[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); +delay_time=t_1(max_pos_after)-t_1(max_pos_before); +printf('Second field delay time = %f uS\n',delay_time); + +%set(gca,'fontsize',40); +set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','fields_before_after_cell_I2.eps') + +figure(4) +I2_max_in=max(I2_out_1(t_good_indx,1)); +I2_max_out=max(I2_out_1(t_good_indx,end)); +I2_in_norm=(I2_out_1(:,1))/I2_max_in; +I2_out_norm=(I2_out_1(:,end))/I2_max_out; +tmin=-0.05; +tmax=0.05; +indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region +plot( ... + t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... + t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... + ) +xlim([tmin,tmax],'manual'); +xlabel('t (uS)') +ylabel('I_2') +title('I_2 before and after cell normalized') +legend('location', 'northeast'); +set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','probe_before_after_cell_I2_normalized.eps') + -- cgit v1.2.3