基于matcom的光栅衍射仿真程序

1，核心代码

double lam=500e-9;//波长
double N=2;//缝数
double a=2e-4;//光栅常数
double D=1;//屏距
double d=1e-3;//缝宽
double ym=2*lam*D/a;//x轴观察范围
double xs=ym;//y轴观察范围
int n=500;//取样点数

Mm ys=linspace(-ym,ym,n);
Mm B;
B=zeros(10,n);

//计算光强分布
for(int i=1;i<=n;i++)
{
 double sinphi=ys.r(i)/D;
 double alpha=3.1415926*a*sinphi/lam;
 double beta=3.1415926*d*sinphi/lam;
 for(int j=1;j<=10;j++)
  B(j,i)=(sin(alpha)/alpha)*(sin(alpha)/alpha)*(sin(N*beta)/sin(beta))*(sin(N*beta)/sin(beta));
    
}

CRect m_client; 
//绘制曲线
winaxes(GetDlgItem(IDC_LINE)->m_hWnd );
GetDlgItem(IDC_LINE)->GetClientRect(m_client);
axesposition(0,0,m_client.Width(),m_client.Height());
Mm m_line=plot((CL(ys),B/4));
set(m_line,TM("color"),TM("blue"));

//绘制图像
winaxes(GetDlgItem(IDC_IMAGE)->m_hWnd );
GetDlgItem(IDC_IMAGE)->GetClientRect(m_client);
axesposition(0,0,m_client.Width(),m_client.Height());
imagesc((ys,xs,CL(B)));

//颜色渲染
Mm color;
int colorIndex=64;
color=zeros(colorIndex,3);
for(int k=1;k<=colorIndex;k++)
 color.r(k,2)=(double)k/colorIndex;
colormap(color);

drawnow();
exitM();

2，运行结果