%闪耀光栅-白光衍射现象
clear;
close all
lambda=[660,610,570,550,460,440,410]*1e-9; %波长
RGB=[1,0,0;1,0.5,0;1,1,0;0,1,0;0,1,1;0,0,1;0.67,0,1]; %七种波长对应的RGB值
Irgb=zeros(150,500,3); %光屏仿真结果图像矩阵,是各个波长的光波衍射结果叠加的结果
Im=zeros(150,500,3); %记录各个波长光衍射结果的RGB值
%定义变量
a=8e-6;N=100;d=4e-5;%光栅常数
I0=100;gamma1=25;i1=gamma1;%闪耀角和入射角相等,自准直入射
gamma=gamma1*pi/180;i=gamma1*pi/180;
for j=1:7
theta=linspace(-0.03*pi,0.03*pi,500);%衍射角度的变化范围
delta=2*pi*d*(sin(theta)+sin(i))/lambda(j);
alpha=pi*a*sin(theta)/lambda(j);
I1=(sinc(alpha/pi)).^2;%单缝衍射因子
I2=(sin(N*delta/2)./sin(delta/2)).^2;%多光束干涉因子
I=I0*I1.*I2;
for k=1:150
Im(k,:,1)=I*RGB(j,1);%把红基色代码计入Iw矩阵红维度
Im(k,:,2)=I*RGB(j,2);%把绿基色代码计入Iw矩阵绿维度
Im(k,:,3)=I*RGB(j,3);%把蓝基色代码计入Iw矩阵蓝维度
end
Irgb=Irgb+Im;%衍射结果叠加
Im=[]; %将Im清零，进行下一个波长的计算
end
Br=1/max(max(max(Irgb)));%调整Irgb矩阵元素的最大值为1的系数
Ib=Irgb*Br;
image(theta,theta,Ib);
xlabel('\gamma /rad');ylabel('y/m');
title(['N=',num2str(N),', 闪耀角\gamma=',num2str(gamma1),'时闪耀光栅的衍射图样']);
