基于matlab的超短脉冲经反射型体光栅后的光束传输特性分析源程序

《基于matlab的超短脉冲经反射型体光栅后的光束传输特性分析源程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、在用Matlab模拟中所用的源程序代码如下:1・光栅厚度对衍射光强的影响clear;closeall%****参数定义区*******c二3*10,8;%真空光速delta_tau=50*10^(-15);%超短脉冲宽度lambdaO=l.5*1(T(-6);%入射脉冲中心波长omega0=2*pi*c/lambda0;%中心角频率T=delta_tau./(2*log(2))0.5;%入射高斯光束常数TLambda=l*10-6;%体光栅周期d=0.1*10=3;%体光栅厚度nl=0.0001;%折射率调制幅值omega=1.l*l(T15:0.0001*l(T15:1.4*1(T15;a

2、O二4.9048;al二0.117680;a2=0.0475;a3=0.027169;%角频率变化区间n_omega=sqrt(a0+(al.Iomega.2*10-12)./((2*pi*c)"2-a2*omega・2*10-12)-a3*(2*pi*c./omega).八2*10八-12);%色散导致折射率变化n_omegaO=sqrt(a0+(a1*omega02*10-12)/((2*pi*c)2-a2*omega02*10^-12)-a3*(2*pi*c/omeg80)"2*l(T-12);%中心频率对应折射率thetl=acos(lambdaO/(2*Lambda*nomegdO

3、));%入射角mu_omega=exp(-1*(T“2*(omega—omegaO)•”2)/4).2%频谱函数costhet2_omega=sqrt(l-(sin(thetl)./nomega).”2);%色散后的入射角phi_omega=omega*nl*d./(2*c.*costhet2_omega);xi_omega=pi2*c*d./(Lambda2.*costhet2_omega).*(1./(omega.*nomega)-1/(omega0*nomega。));eta_omega=(1+(1-(xiomega./phi_omega).”2)./sinh(sqrt(phi_ome

4、ga.八2-%取最大值%归一化频谱效率etaO=max(eta_omega);etal=eta_omega/etaO;1s_omega=mu_omega.*etal;plot(omega,muomega,')，holdonomega=1.1*10^15:0.00001*10^15:1.4*10^15;%角频率变化区间aO二4.9048;al二0.117680;a2=0.0475;a3=0.027169;plot(omega,Isomega,'--k'),hold%****参数定义区*******onc=3*10'8;%真空光速delta_tau=50*10^(-15);%超短脉冲宽度lamb

5、da0=l.5*1(T(-6);%入射脉冲中心波长omega0=2*pi*c/lambda0;%中心角频率T=delta_tau./(2*log(2)厂0.5;%入射高斯光束常数TLambda=l*10-6;%体光栅周期d=0.3*10=3;%体光栅厚度nl=0.0001;%折射率调制幅值n_omega=sqrt(a0+(al.Iomega.2*10-12)./((2*pi*c)"2-a2*omega・2*10-12)-a3*(2*pi*c./omega).八2*10八-12);%色散导致折射率变化n_omega0=sqrt(a0+(a1*omega02*10-12)/((2*pi*c)2-

6、a2*omega02*10^-12)-33*(2*pi*c/omegaO)"2*l(T-12);%中心频率对应折射率thetl=acos(lambdaO/(2*Lambda*nomegdO));%入射角mu_omega=exp(-1*(T2*(omega-omegaO).”2)/4).“2;%频谱函数costhet2_omega=sqrt(l-(sin(thetl)./nomega).八2);%色散后的入射角phi_omega=omega*nl*cl・/(2*c.*costhet2_omega);xi_omega=pi^2*c*d./(Lambda2.*costhet2_omega).*(1

7、・/(omega.*nomega)-1/(omegaO*n_omegaO));eta_omega=(1+(1-(xiomega・/phiomega).2)./sinh(sqrt(phiomega・2~xi_omega."2)).”2).=1;%频谱衍射效率etaO=max(etaomega);%取最大值etal=eta_omega/etaO;%归一化频谱效率1s_omega=mu_omega・*etal;pl