前言
许多波导结构，例如布拉格光栅，包含着沿传播方向（z方向）的扰动。标准的单向BPM算法不能计算由前向和后向传输场之间的这些扰动引起的耦合。为了适应这些情况，在BeamPROP中实现了几种双向BPM算法。这些技术考虑了两个方向上的波传播，并提供了确定稳态场的方法。此功能目前仅适用于2D仿真，并且最适用于但不限于index contrast（扰动）小且满足近轴条件（前向和后向波沿着或几乎沿着z轴传播）的情况。
仿真部分
1、创建布局
2、绘制波导
光栅结构由三段定义：输入段，光栅段和输出段。
3、设置光栅中折射率周期性变化
查看设定好的折射率分布
4、设定路径，监控和发射场
5、运行仿真
6、扫描波长进行仿真
最后的透射谱和反射谱如下图
但文档里给出的是下图，运行给出的示例出来的谱，跟上面两个一样，出不来下面的这么好看的图。(还是用MATLAB做吧，，，，，)
参考文档
beamprop.pdf