【正文】 old on。 endlegend(39。d=39。,39。d=39。,39。d=39。,39。d=39。)。title(39。不同光柵厚度下的波長選擇性39。)。xlabel(39。波長變化（mm）39。)。ylabel(39。歸一化的衍射效率39。)。grid onendfor i=1:4 %不同光柵線對下的波長選擇性 switch i case 1 kai=1/800。 case 2 kai=1/1000。 case 3 kai=1/1200。 case 4 kai=1/1400。 end d=1。 lambda=lambda0+deltalambda。 theta=asin((lambda0*n)/(2*kai))。 nu=(pi*n1*d)./(lambda.*cos(theta))。 xi=(deltalambda*pi*d)./(2*n*(kai^2)*cos(theta))。 eta=(sin(sqrt(nu.^2+xi.^2)).^2)./(1+(xi./nu).^2)。 a=(pi*n1*d)./(lambda0*cos(theta))。 b=(sin(a))^2。 c=eta./b。 figure(4) subplot(2,2,i)。 plot(deltalambda,c,39。r39。)。 switch i case 1 title(39。kai=1/800 mm39。)。 case 2 title(39。kai=1/1000 mm39。)。 case 3 title(39。kai=1/1200 mm39。)。 case 4 title(39。kai=1/1400 mm39。)。 end axis([5*10^6,5*10^6,0,1])。 xlabel(39。波長變化（mm）39。)。 ylabel(39。歸一化的衍射效率39。)。 grid onfigure(8) switch i case 1 plot(deltalambda,c)。 hold on。 case 2 plot(deltalambda,c,39。39。)。 hold on。 case 3 plot(deltalambda,c,39。:39。)。 hold on。 case 4 plot(deltalambda,c,39。.39。)。 hold on。 endlegend(39。kai=1/800 mm39。,39。kai=1/1000 mm39。,39。kai=1/1200 mm39。,39。kai=1/1400 mm39。)。title(39。不同光柵線對下的波長選擇性39。)。xlabel(39。波長變化（mm）39。)。ylabel(39。歸一化的衍射效率39。)。grid onend（3）單色發(fā)散光束經(jīng)透射型布拉格體光柵的特性matlab程序clc。clear。syms t。d=1。n=。n1=。%介質(zhì)的折射率調(diào)制度lambda=*10^3。c=linspace(0,0,1000)。b=linspace(10^5,10^1,1000)。for num=1:3 %不同光柵線對 switch num case 1 kai=1/500。 case 2 kai=1/1000。 case 3 kai=1/2000。 end theta0=asin((lambda)/(2*kai*n))。 for i=1:1000 f=@(t)(((sin(sqrt((pi*n1*d./(lambda*cos(t))).^2+(t*pi*d*cos(theta0)./(kai*cos(t))).^2)).^2)./(1+((t*pi*d*cos(theta0)./(kai*cos(t)))./((pi*n1*d)./(lambda*cos(t)))).^2)).*(exp(2*((t)./b(i)).^2)))。 c(i)=quadl(f,1,1)。%積分 end c=c./b。 c=c*sqrt(2/pi)。 c=c/max(c)。%歸一化 figure(9) switch num case 1 semilogx(b,c)。 hold on。 case 2 semilogx(b,c,39。39。)。 hold on。 case 3 semilogx(b,c,39。.39。)。 hold on。 end endtitle(39。單色發(fā)散光束經(jīng)透射型體光柵衍射效率與發(fā)散角關(guān)系圖39。)。xlabel(39。發(fā)散角（rad）39。)。ylabel(39。歸一化的衍射效率39。)。legend(39。500/mm39。,39。1000/mm39。,39。2000/mm39。)。（4）多色平面波經(jīng)透射型布拉格體光柵的特性matlab程序clc。clear。syms t。n=。d=1*10^6。n1=。%介質(zhì)的折射率調(diào)制度lambda0=*10^3。m=lambda0*(11*10^12)。c=linspace(0,0,1000)。w=linspace(10^2,10^2,1000)。for num=1:3 switch num case 1 kai=1/500*1000000。 case 2 kai=1/1000*1000000。 case 3 kai=1/2000*1000000。 end theta=asin((lambda0)/(2*kai*n))。 for i=1:1000 f=@(t)(((sin(sqrt(((pi*n1*d)./(t.*cos(theta))).^2+(((tlambda0)*pi*d)./(2*n*kai^2*cos(theta))).^2))).^2)./(1+((((tlambda0)*pi*d)./(2*n*kai^2*cos(theta)))./((pi*n1*d)./(t.*cos(theta)))).^2).*(exp(2*((tlambda0)./w(i)).^2)))。 c(i)=quadl(f,lambda0m,lambda0+m)。 end c=c./w。 c=c/max(c)。 figure(10) switch num case 1 semilogx(w,c)。 hold on。 case 2 semilogx(w,c,39。39。)。 hold on。 case 3 semilogx(w,c,39。.39。)。 hold on。 end endtitle(39。多色平面波衍射效率與角譜寬度關(guān)系圖39。)。xlabel(39。角譜寬度（nm）39。)。ylabel(39。歸一化的衍射效率39。)。legend(39。500/mm39。,39。1000/mm39。,39。2000/mm39。