blob: f17d0d110ce77e2d5fc50d0fd50d9ae035013984 (
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
|
basis_transformation;
fnames=glob('results/*.mat');
Nsteps=length(fnames)
B_field=zeros(1,Nsteps);
xi_linear=zeros(1,Nsteps);
xi_left=zeros(1,Nsteps);
xi_right=zeros(1,Nsteps);
% read the information from resulting files
for i=1:Nsteps
d=load(fnames{i});
B_field(i)=d.B_field;
xi_linear(i)=d.xi_linear;
xi_left(i)=d.xi_left;
xi_right(i)=d.xi_right;
endfor
figure(1);
plot(B_field, real((xi_right+xi_left))/sqrt(2));
title('phase of linear x-polarization');
figure(2);
plot(B_field, imag((xi_right+xi_left))/sqrt(2));
title('absorption of linear x-polarization');
figure(3);
plot(B_field, real((1i*xi_right-1i*xi_left))/sqrt(2));
title('phase of linear y-polarization');
figure(4);
plot(B_field, imag((1i*xi_right-1i*xi_left))/sqrt(2));
title('absorption of linear y-polarization');
|
