beifen.txt

来自「一个激光谐振腔的稳态过程模拟程序。通过设置波长」· 文本 代码 · 共 48 行

for(i=0;i<M;i++)
{
	buffer_cos=0;
	buffer_sin=0;
	for(m=0;m<M;m++)
	{
		buffer_cos+=cos(c*(u2[i].x-u1[m].x)*(u2[i].x-u1[m].x)
			-u1[m].P)*u1[m].A*step;
		buffer_sin-=sin(c*(u2[i].x-u1[m].x)*(u2[i].x-u1[m].x)
			-u1[m].P)*u1[m].A*step;
	}
	u2[i].x=-a+2*a*i/M;
	u2[i].GetAmplitude(buffer_cos,buffer_sin);
	u2[i].GetPhrase(buffer_sin,buffer_cos);

}

//归一化
for(i=0;i<M;i++)
{
	u2[i].A/=u2[M/2].A;
	u2[i].P-=u2[M/2].P;
}


/////////////u2---->u1的衍射//////////
for(i=0;i<M;i++)
{
	buffer_cos=0;
	buffer_sin=0;
	for(m=0;m<M;m++)
	{
		buffer_cos+=cos(c*(u2[m].x-u1[i].x)*(u2[m].x-u1[i].x)
			-u2[m].P)*u2[m].A*step;
		buffer_sin+=sin(c*(u2[m].x-u1[i].x)*(u2[m].x-u1[i].x)
			-u2[m].P)*u2[m].A*step;
	}
	u1[i].x=-a+2*a*i/M;
	u1[i].GetAmplitude(buffer_sin,buffer_cos);
	u1[i].GetPhrase(buffer_sin,buffer_cos);
}

//归一化
for(i=0;i<M;i++)
{
	u1[i].A/=u1[M/2].A;
	u1[i].P-=u1[M/2].P;
}