相位恢复算法

诺兮

Mr.rui 发表于 2025-8-18 15:51
你好，我这个传统的激光全息光学，我的系统是用硬X射线同轴法得到的相衬图对其进行相位恢复，光源是微焦 ...

我边看你的代码边写哈，可能会乱，会列出一些我觉得有些奇怪的地方，希望能对你有些启发，不一定对；

1.A你写的是均匀不是随机分布，注释里又写的随机，怪欸；
2.输入x-ray source应该不是均匀输入吧，看你的迭代过程，正向传播过程应该是object到detector（在我这的理解，第一个傅里叶变换是object变成频谱，第二个反傅里叶变换是频谱经过传播后在detector上的强度分布，通过对detector的强度约束进行迭代），如果是这样理解的话，Ei处设置的振幅分布（虽然你是直接写了一个相位，但是这实际默认的振幅分布是均匀分布，相当于初始振幅分布给的是ones(N,N);)和你的实验模型并不一致，你可以在实验系统中source前加透镜让其变为平行光，可能更符合你写的。论文中没有加可能是距离比较大导致输入光近似认为平行光？我只是这样猜测。
3.关于流程图，正向传播部分你做的是对的，通过物面振幅替换掉传播后的振幅分布，反向传播的振幅分布也是对的，用均匀振幅分布替换反向传播计算得到的振幅分布，但是反向传播后的相位分布原理图中的那个箭头我并不清楚相位分布如何替换的（我接触到的绝大部分是直接保留上一次的相位分布，虽然也有类似双限制相位算法分信号区与非信号区有对相位进行不同的约束，但是真的很少），你可以着重研究一下反向传播中相位约束条件这部分；
4.你最后的结果图没给代码，我提醒一下你得到的结果貌似是Ei这个物面时域，需要你正向传播走一次再进行振幅归一化（abs后再/max(max()))才能得到正确的仿真结果；