1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
|
x=0.1*10^-6;
w=0.08*10^-6;
Lambda=400e-9:10e-9:800e-9; % en nm
c=3e8;
h= 4.135e-15;
E=h.*c./Lambda;
C = [5030.02, 483.916;1634.30,79.4079]; % en cm-1
C3= 1046.08; % en cm-1
Ek = 1.09969; % en eV
E3 = 1.40985; % en eV
dE = [0.0583148,0.0220161]; % en eV
N = 0.394122;
dN = 1.23084; % en eV-1.
a=x+w;
for w = 1:2
for j=1:2
Alpha= C(j,w)* power(E-Ek+power(-1,j)*dE(w)+abs(E-Ek+power(-1,j)*dE(w))/2,2)./E+C3*power((E-E3)+abs((E-E3)),N+dN*E)./E
end
end
exp (-Alpha.*(x+w)); |
Partager