pipeijiao.m

来自「BBO I类匹配角计算 通过修改相应的参数可以计算其他负单轴晶体的I类匹配角」· M 代码 · 共 13 行

for i=1:1:1000
    lamuda(i)=0.4+0.001*i;
    n2o_w=2.7359+0.01878/(lamuda(i)^2-0.01822)-0.01354*lamuda(i)^2;
    n2o_2w=2.7359+0.01878/((lamuda(i)/2)^2-0.01822)-0.01354*(lamuda(i)/2)^2;
    n2e_w=2.3753+0.01224/(lamuda(i)^2-0.01667)-0.01516*lamuda(i)^2;
    n2e_2w=2.3753+0.01224/((lamuda(i)/2)^2-0.01667)-0.01516*(lamuda(i)/2)^2;
    theta(i)=asin((n2e_2w/n2o_w)*(n2o_2w-n2o_w)/(n2o_2w-n2e_2w))*180/pi;
end ;   
plot(lamuda,theta,'rx')
title('BBO晶体倍频I类匹配角随波长的变化')
xlabel('入射光的波长（um）')
ylabel('匹配角大小(\circ)')