Cpp五条

这个屌丝很懒，什么也没留下！

热门标签

全息图算法：加权Gerchberg-Saxton算法（GSW）_gsw算法

作者：Cpp五条 | 2024-04-20 04:43:17

踩

gsw算法

背景：GS（Gerchberg-Saxton）在计算全息图、产生光阱方面，对比于直接算法，包括随机相位掩膜法(RM)、棱镜透镜叠加(S)算法，随机叠加(SR)算法，具有光能利用率高的优点，但是和GAA算法一样，都具有一个缺点，即光能的分布上，均匀性不够好。从该缺陷出发，便产生了加权的GS算法(GSW)算法。

本文只讨论纯相位型空间光调制器。设入射在空间光调制器(SLM)的光波为均匀平面波，设振幅为 $\left | u \right |$ ,则从SLM出射后的光波振幅为 $u_{j}=\left | u \right |exp(i\varphi _{j})$ ,其中， $\varphi _{j}$ 为SLM上第j个像素产生的相移。使用标量衍射理论，得到从SLM上第i行第j列像素的光传播到焦面(focal plane)上第m个光阱处的光振幅为 $\nu _{m}$

$\nu _{m}=\frac{exp(ik(2f+z_{m})))\cdot d^{2}}{i\lambda f}\sum_{j=1}^{N}\left | u \right |exp(i(\varphi_{j}-\Delta _{j}^{m} )))$

$\Delta_{j}^{m}=\frac{k z_{m}}{2f^2}(x_{j}^2+y_{j}^2)+\frac{k}{f}(x_jx_m+y_jy_m)$

其中，k为波矢，有 $k=\frac{2\pi}{\lambda}$ ,λ为波长， $\Delta_j^m$ 实际上就表示的是位于 $(x_m,y_m.z_m)$ 处的光阱经过菲涅耳变换后得到的相位信息，右端包含了其空间频率信息，左端包含了轴向的信息，是一个透镜项。f为傅里叶变换透镜的焦距，d为像素边长，N为像素个数。引入无量纲量 $V_m$

$V_m=\sum_{j=1}^{N}\frac{1}{N}exp(i\varphi_j-\Delta_j^m)$

其物理意义通过 $I_m=\left | V_m \right |^{2}$ 来表示， $I_{m}$ 表示第m个光阱能量与总能量之比。计算全息的算法的意义在于对于给定的m个位于 $(x_m,y_m.z_m)$ 的光阱，或者给定的 $\Delta_{j}^{m}$ (两者是等效的)，给出最优的 $\varphi_{j}$ ，使得各个光阱的强度最接近于设计值。对于GS算法、GAA算法这两种迭代算法，其以 $\sum\left | V_m \right |$ 达到最大值为目的，这也是导致其计算产生的光阱强度不均匀的间接原因之一。

首先将GS算法的思路作为切入，为了使以 $\varphi_{j}$ 为自变量的函数 $\sum\left | V_m \right |$ 达到最大， $\sum\left | V_m \right |$ 对

$\varphi_{j}$ 求偏微分得到

$\frac{\partial }{\partial \varphi_j}\sum\left | V_m \right |=Re\left \{ \frac{iexp(i\varphi_j))}{N}\sum_mexp(-i\Delta_j^m)\frac{V_m^*}{\left | V_m \right |} \right \}=0$

最终得到（原文中用了Hessian矩阵，博主这不太明白QAQ，直接略过了）

$\varphi_j=arg\left [ \sum_mexp(i\Delta_j^m)V_m/\left | V_m \right | \right ]$

而GSW算法，是在该表达式中引入了一个隐式的权重 $w_m$ ,表达式变为

$\varphi_j=arg\left [ \sum_mexp(i\Delta_j^m)w_mV_m/\left | V_m \right | \right ]$

而 $w_m$ 在初次迭代时为1，在第k次迭代中， $w_m^k=w_m^{k-1}\frac{\left \langle \left | V_m^k-1 \right | \right \rangle}{\left | V_m^k-1 \right | }$ ，其中 $\left \langle .. \right \rangle$ 表示求平均。

在对算法基本原理讲解后，大家可能还是一头雾水的(反正博主是一头雾水的QAQ)，github上找到一个python编写的GSW算法源码后，给了我很大启发。如下是我绘制的算法框图，简单明了的讲解了这篇文章想要表达的一个思路。

GSW算法的流程图如图

输入图像灰度值作为振幅输入，并对图像每一像素赋予随机相位，作为输入图像。

其与GS算法不同的是，使用了随迭代次数改变的权重系数，即，输入图像中光阱的能量与每次迭代得到的复原图像中光阱能量之比值乘上一次迭代中的权重系数，作为本次迭代的权重。权重的引入，使得全息图还原得到的复原图的光阱能量会朝设计值靠近。

实验：如图所示为设计的图像

如下图为计算得到的计算全息图

下图为使用全息图还原得到的还原图

下图为还原图与原图的相关系数随迭代次数的图像

本篇文章中可能还有很多问题，欢迎大家在评论区讨论

参考文献:

[1]Di Leonardo Roberto,Ianni Francesca,Ruocco Giancarlo. Computer generation of optimal holograms for optical trap arrays.[J]. Optics express,2007,15(4).

声明：本文内容由网友自发贡献，不代表【wpsshop博客】立场，版权归原作者所有，本站不承担相应法律责任。如您发现有侵权的内容，请联系我们。转载请注明出处：https://www.wpsshop.cn/w/Cpp五条/article/detail/455425