耦合波理论计算体全息衍射效率matlab程序.pdf

体全息光栅衍射效率计算程序

耦合波理论

半高全宽计算

程序如下：

clear

formatlong

lam0=0.532;%um,%记录

k0=2*pi/lam0;

n=1.68;

k=n*k0;

d=40;%um厚度

n0=n;%曝光后的平均折射率

n1=0.003;%折射率调制度

alpha0=0;%平均折射率虚部

alpha1=0;%折射率虚部调制度

couple_mode=prism;%couple_mode=air;

beta=[-503];

beta1=beta(1);

beta2=beta(2);%左侧入射为beta1，右侧为beta2+z朝右+x朝上

switchcouple_mode

caseair

beta1_media=asind(sind(beta1)/n);

beta2_media=asind(sind(beta2)/n);

k1=k*[cosd(beta1_media)sind(beta1_media)];%[zx]

k2=k*[cosd(180+beta2_media)sind(180+beta2_media)];

K=k1-k2;

caseprism

beta1_media=beta1;

beta2_media=beta2;

k1=k*[cosd(beta1_media)sind(beta1_media)];

k2=k*[cosd(180+beta2_media)sind(180+beta2_media)];

K=k1-k2;

end

norm_K=norm(K);

lam_p0=0.535;%衍射过程

kp0=2*pi/lam_p0;

kp=n0*kp0;

uni_norm_vec=[10];

theta_pin=(-30:0.1:30);

theta_p=asind(sind(theta_pin)/n0);%介质外部入射到内部

kp_vec=kp*[cosd(theta_p)sind(theta_p)];

forii=1:length(theta_pin)

kp_vec_ii=[kp_vec(ii,1),kp_vec(ii,2)];

kappa_TE=n1/n0/2*kp-1j*alpha1/2;%TEwave%--coupledwavetheory

kappa=kappa_TE;

zeta=(2*dot(kp_vec_ii,K)-norm_K^2)/(2*kp);

衍射光方向余弦，以为起始轴

cs=dot((kp_vec_ii-K),uni_norm_vec)/kp;%+z

入射光方向余弦，以为起始轴

cr=dot(kp_vec_ii,uni_norm_vec)/kp;%+z

gamma1=-1/2*(alpha0/cs+alpha0/cr+1j*zeta/cs)...

+1/2*sqrt((alpha0/cr-alpha0/cs-1j*zeta/cs)^2-4*kappa^2/(cr*cs));

gamma2=-1/2*(alpha0/cs+alpha0/cr+1j*zeta/cs)...

-1/2*sqrt((alpha0/cr-alpha0/cs-1j*zeta/cs)^2-4*kappa^2/(cr*cs));

E_refl_1order=-1j*kappa/(alpha0+1j*zeta+...

cs*(gamma1*exp(gamma2*d)-gamma2*exp(gamma1*d))/(exp(gamma2*d)-

exp(gamma1*d))...

);

E_refl_0order=cs*(gamma1-gamma2)*...

((alpha0+1j*zeta+cs*gamma1)*exp(-gamma1*d)-...

(alpha0+1j*zeta+cs*gamma2)*exp(-gamma2*d)...

)^-1;

eta_refl_1order(ii,1)=real(abs(cs)/cr*E_refl_1order*conj(E_refl_1order));%前面的因子是由

于入射光和衍射光角度不同带来的因子，与投影面积相关

eta_refl_0order(ii,1)=E_refl_0order*con