完整瑞利-索末菲模型的 Matlab 实现，比较适合搞光学模拟的你参考。代码清晰，逻辑直接，用起来挺顺手。适合做圆形孔径的近场衍射计算，不用担心孔径尺寸或传播距离的限制，场景覆盖得蛮广的。

Gilllen 那篇 AJP 老论文里提到的近场传播计算方法，这段代码就是对它的一个还原。像那种需要沿 z 轴观察不同传播距离下的强度分布场景——比如检测器放在(0,0,z)点，z 从 1000a 一路减小——这个模型就挺派得上用场。

代码默认圆孔径半径是100 µm，波长设定为10 µm，入射是平面波，电场幅度是1。这些数值在实验或仿真中也挺常见的。核心其实就是个积分计算，轴上强度的结果也蛮稳定。

哦对了，这玩意比较适合搞Matlab的同学，代码里没花活，也没有太多调参门槛。响应也快，适合拿来做教学演示或者论文复现用。

顺便推荐几个相关资源，如果你想深入搞衍射建模，可以看看：

如果你在做光学模拟项目，是近场类的，可以考虑把这段代码塞进工具库里，省时省力还靠谱。