光栅衍射C语言程序,基于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，运行结果