[论文解读] Fast and Accurate Optimization of Metasurfaces with Gradient Descent and the Woodbury Matrix Identity

Jordan Budhu, Anthony Grbic|arXiv (Cornell University)|Jul 7, 2021

Advanced Antenna and Metasurface Technologies被引用 4

一句话总结

本文提出了一种基于梯度的快速且精确的优化方法，用于有限尺寸超表面，采用矩量法（Method of Moments）和伍德伯里矩阵恒等式（Woodbury Matrix Identity）。通过在梯度计算过程中利用伍德伯里恒等式高效计算扰动矩阵的逆，该方法相较于标准数值微分实现了高达26.5倍的加速，从而实现了对无源、无损耗宽角反射超表面的高效设计。

ABSTRACT

A fast metasurface optimization strategy for finite-size metasurfaces modeled using integral equations is presented. The metasurfaces considered are constructed from finite patterned metallic claddings supported by grounded dielectric spacers. Integral equations are used to model the response of the metasurface to a known excitation and solved by Method of Moments. An accelerated gradient descent optimization algorithm is presented that enables the direct optimization of such metasurfaces. The gradient is normally calculated by solving the method of moments problem N+1 times where N is the number of homogenized elements in the metasurface. Since the calculation of each component of the N-dimensional gradient involves perturbing the moment method impedance matrix along one element of its diagonal and inverting the result, this numerical gradient calculation can be accelerated using the Woodbury Matrix Identity. The Woodbury Matrix Identity allows the inverse of the perturbed impedance matrix to be computed at a low cost by forming a rank-r correction to the inverse of the unperturbed impedance matrix. Timing diagrams show up to a 26.5 times improvement in algorithm times when the acceleration technique is applied. An example of a passive and lossless wide-angle reflecting metasurface designed using the accelerated optimization technique is reported.

研究动机与目标

• 开发一种基于积分方程建模的有限尺寸超表面计算高效优化框架。
• 克服在超表面设计中数值梯度计算带来的高计算成本。
• 利用矩阵恒等式加速在梯度评估过程中扰动阻抗矩阵的求逆运算。
• 实现在高精度下、运行时间显著减少的超表面大规模优化。
• 在无源、无损耗宽角反射超表面设计中验证该方法。

提出的方法

• 采用矩量法（MoM）结合积分方程对超表面进行建模，以描述其电磁响应。
• 通过有限差分法计算梯度，需对N个单元分别进行N+1次MoM阻抗矩阵求逆。
• 应用伍德伯里矩阵恒等式，将扰动阻抗矩阵的逆表示为对原始逆矩阵的低秩修正。
• 该方法将每个梯度分量的计算成本从O(N³)降低至每个单元O(N²)，显著加速优化过程。
• 将加速后的梯度下降算法应用于优化超表面单元结构，以实现所需的波前调控。
• 在一种无源、无损耗的金属图案在介质基底上的宽角反射超表面结构上验证了该方法。

实验结果

研究问题

• RQ1 伍德伯里矩阵恒等式是否可有效用于加速超表面优化中的梯度计算？
• RQ2 基于矩阵恒等式的求逆方法在梯度驱动的超表面优化中可实现多大程度的加速？
• RQ3 该加速方法能否在复杂波前调控下保持对无源、无损耗超表面设计的高精度？
• RQ4 该方法在计算效率和收敛速度方面的性能极限是什么？
• RQ5 该方法能否实现具有宽角响应特性的有限尺寸超表面的实际设计？

主要发现

• 使用伍德伯里矩阵恒等式将每个梯度分量的计算成本从O(N³)降低至O(N²)，从而实现更快的优化。
• 加速算法相较于标准数值梯度计算，运行时间最高可提升26.5倍。
• 该方法保持了高精度，成功实现了对无源、无损耗宽角反射超表面的设计。
• 时序图显示，该方法在不同尺寸和配置的超表面中均表现出一致的性能提升。
• 该方法具有良好的可扩展性，适用于具有复杂单元结构的大规模超表面优化。
• 该方法可直接对有限尺寸超表面进行优化，无需依赖均质化近似。