-
Notifications
You must be signed in to change notification settings - Fork 3
/
input_ARcoating.m
73 lines (47 loc) · 1.7 KB
/
input_ARcoating.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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%% AR coating structure %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
nL=1;
nR=3;
l1=3e-6;
n1=1;
l3=5e-6;
n3=3;
n2=1.732+0.0i; % sqrt(n3)
lambda0=850e-9; % Central wavelength
l2=lambda0/(4*abs(n2)); % thickness at lambda/4
%alpha2=4*pi*imag(n2)./lambda;
layer=[
l1 n1
l2 n2
l3 n3
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%% Formula computation %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r12 = (n1-n2) / (n1+n2);
r23 = (n2-n3) / (n2+n3);
t12 = 2*n1 / (n1+n2);
t23 = 2*n2 / (n2+n3);
D1 = zeros(2,2,length(lambda));
D1(1,1,:) = exp(+1i*2*pi*n2*l2./lambda) ;
D1(2,2,:) = exp(-1i*2*pi*n2*l2./lambda) ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
P1(1,1,:) = (1./t12) ;
P1(1,2,:) = (1./t12) .* r12 ;
P1(2,1,:) = (1./t12) .* r12 ;
P1(2,2,:) = (1./t12) ;
P2(1,1,:) = (1./t23) ;
P2(1,2,:) = (1./t23) .* r23 ;
P2(2,1,:) = (1./t23) .* r23 ;
P2(2,2,:) = (1./t23) ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j=1:length(lambda)
S(:,:,j)=P1*D1(:,:,j)*P2;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Rf = ( abs( S(1,2,:)./S(2,2,:) ) ).^2;
Rf = reshape(Rf,1,length(lambda));
Tf = abs( 1./S(2,2,:) ).^2 ;
Tf = real(n3./n1).*reshape(Tf,1,length(lambda));
Af=1-exp(-4*pi*imag(n2)./lambda * l2);