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
|
% suppress output of plots on display
set (0, "defaultfigurevisible", "off")
realistic_Rb_and_fields
%% 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, Ip2_out_1); colorbar
tmin=-0.4;
tmax= 0.4;
ylim([tmin,tmax],'manual');
xlabel('z (cm)')
ylabel('t (uS)')
zlabel('I_p2')
title('I_p2')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Ip2.eps')
figure(1)
%set(gca,'fontsize',20);
imagesc(z_1, t_1, Im2_out_1); colorbar
tmin=-0.4;
tmax= 0.4;
ylim([tmin,tmax],'manual');
xlabel('z (cm)')
ylabel('t (uS)')
zlabel('I_m2')
title('I_m2')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Im2.eps')
%% field I4 propagation
figure(3)
imagesc(z_1, t_1, Ip4_out_1); colorbar
tmin=-0.4;
tmax= 0.4;
ylim([tmin,tmax],'manual');
xlabel('z (cm)')
ylabel('t (uS)')
zlabel('I_p4')
title('I_p4')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Ip4.eps')
figure(4)
imagesc(z_1, t_1, Im4_out_1); colorbar
tmin=-0.4;
tmax= 0.4;
ylim([tmin,tmax],'manual');
xlabel('z (cm)')
ylabel('t (uS)')
zlabel('I_m4')
title('I_m4')
print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Im4.eps')
%% fields before and after the cell
figure(5)
%set(gca,'fontsize',30);
plot( ...
t_1,Ip2_out_1(:,1),'.-;I_p2 before;', "linewidth", 4 ...
,t_1,Im2_out_1(:,1),'.-;I_m2 before;', "linewidth", 4 ...
,t_1,Ip2_out_1(:,end), '-;I_p2 after;', "linewidth", 4 ...
,t_1,Im2_out_1(:,end), '-;I_m2 after;', "linewidth", 4 ...
,t_1,Ip4_out_1(:,end), '-;I_p4 after;', "linewidth", 4 ...
,t_1,Im4_out_1(:,end), '-;I_m after;', "linewidth", 4 ...
)
xlabel('t (uS)')
ylabel('I (1/s)^2')
title('Fields before and after cell')
legend('location', 'northeast');
print('-color','-depsc2', 'fields_before_after_cell.eps')
fid_stamp=fopen ("pp_I2.stamp", "a+");
fclose (fid_stamp);
|
