summaryrefslogtreecommitdiff
path: root/xmds2/Genas_system/pp_fields.m
blob: c87f73cbbc5952b1223deeb3453e338a4b50fcda (plain)
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
Genas_system

Ex=Ex_re_out_1 + 1i*Ex_im_out_1;
Ey=Ey_re_out_1 + 1i*Ey_im_out_1;

Ecw= Ex+1i*Ey;
Eccw=Ex-1i*Ey;

%% field propagation
t_1=t_1*1e9; % time now measured in uS
figure(1)
subplot(1,2,1); imagesc(z_1, t_1, abs(Ecw)); colorbar
xlabel('z')
ylabel('t (nS)')
zlabel('|E_{cw}|')
title('|E_{cw}|')
subplot(1,2,2); imagesc(z_1, t_1, abs(Eccw)); colorbar
xlabel('z')
ylabel('t (nS)')
zlabel('|E_{ccw}|')
title('|E_{ccw}|')

print('fields_map.pdf')

% fields before and after
[Nr,Nc] = size( Ecw );
Ecw_b=Ecw(:,1);
Ecw_a=Ecw(:,Nc);
Eccw_b=Eccw(:,1);
Eccw_a=Eccw(:,Nc);
figure(2)
subplot(2,1,1);
plot(t_1, abs(Ecw_b), '-;E_{cw_{before}};', t_1, abs(Ecw_a), '-;E_{cw_{after}};'); 
xlabel('t (nS)')
ylabel('Rabi frequency (1/s)')
subplot(2,1,2);
plot(t_1, abs(Eccw_b), '-;E_{ccw_{before}};', t_1, abs(Eccw_a), '-;E_{ccw_{after}};')
xlabel('t (nS)')
ylabel('Rabi frequency (1/s)')

print('fields_before_after.pdf')