-
Notifications
You must be signed in to change notification settings - Fork 0
/
Ag_experimental.m
59 lines (56 loc) · 1.17 KB
/
Ag_experimental.m
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
function eps=Ag_experimental(wavelength_um)
Ag_data=...
[0.1879 1.07 1.212;
0.1916 1.1 1.232;
0.1953 1.12 1.255;
0.1993 1.14 1.277;
0.2033 1.15 1.296;
0.2073 1.18 1.312;
0.2119 1.2 1.325;
0.2164 1.22 1.336;
0.2214 1.25 1.342;
0.2262 1.26 1.344;
0.2313 1.28 1.357;
0.2371 1.28 1.367;
0.2426 1.3 1.378;
0.249 1.31 1.389;
0.2551 1.33 1.393;
0.2616 1.35 1.387;
0.2689 1.38 1.372;
0.2761 1.41 1.331;
0.2844 1.41 1.264;
0.2924 1.39 1.161;
0.3009 1.34 0.964;
0.3107 1.13 0.616;
0.3204 0.81 0.392;
0.3315 0.17 0.829;
0.3425 0.14 1.142;
0.3542 0.1 1.419;
0.3679 0.07 1.657;
0.3815 0.05 1.864;
0.3974 0.05 2.07;
0.4133 0.05 2.275;
0.4305 0.04 2.462;
0.4509 0.04 2.657;
0.4714 0.05 2.869;
0.4959 0.05 3.093;
0.5209 0.05 3.324;
0.5486 0.06 3.586;
0.5821 0.05 3.858;
0.6168 0.06 4.152;
0.6595 0.05 4.483;
0.7045 0.04 4.838;
0.756 0.03 5.242;
0.8211 0.04 5.727;
0.892 0.04 6.312;
0.984 0.04 6.992;
1.088 0.04 7.795;
1.216 0.09 8.828;
1.393 0.13 10.1;
1.61 0.15 11.85;
1.937 0.24 14.08];
buf=(Ag_data(:,2)+1j*Ag_data(:,3)).^2;
Ag_data(:,4)=real(buf);
Ag_data(:,5)=imag(buf);
eps=interp1(Ag_data(:,1),Ag_data(:,4), wavelength_um,'spline') -1j*interp1(Ag_data(:,1),Ag_data(:,5), wavelength_um,'spline');
end