1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
% suppress output of plots on display
set (0, "defaultfigurevisible", "off")
realistic_Rb
%% field I2 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')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps')
desired_x_size=200;
desired_y_size=200;
xskip=ceil(length(z_1)/desired_x_size);
yskip=ceil(length(t_1)/desired_y_size);
map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
%% field I4 propagation
figure(2)
imagesc(z_1, t_1, I4_out_1); colorbar
tmin=-0.4;
tmax= 0.4;
ylim([tmin,tmax],'manual');
xlabel('z (cm)')
ylabel('t (uS)')
zlabel('I_4')
title('I_4')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps')
map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip);
%% fields before and after the cell
figure(3)
%set(gca,'fontsize',30);
plot( ...
t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4 ...
,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ...
,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ...
)
xlabel('t (uS)')
ylabel('I (1/s)^2')
title('Fields 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.eps')
print('-color','-depsc2', 'fields_before_after_cell.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));
I4_max_out=max(I4_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;
I4_out_norm=(I4_out_1(:,end))/I4_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),'.-;I_2 before;', "linewidth", 4, ...
t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ...
,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ...
)
legend('location', 'southeast');
xlim([tmin,tmax],'manual');
ylim([0.5,1],'manual');
xlabel('t (uS)')
ylabel('I normalized')
title('Fields 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', 'fields_before_after_cell_normalized.eps')
print('-color','-depsc2', 'fields_before_after_cell_normalized.eps')
fid_stamp=fopen ("pp_I2.stamp", "a+");
fclose (fid_stamp);
|
