激光超声成像算法

 1、Fourier-domain backpropagation（傅里叶反向传播）

1.1 Fourier transforms

1.1.1 平面波的波动方程解

        在二维几何结构中，假设波动方程的解为平面波的形式：

其中：

•          表示在时间 位置 和 处的波场，它通常表示在该时刻和位置上波的振幅或物理量的强度；
•          和  分别是沿方向和方向的波数，分别表示各自方向上的空间频率。波数的单位通常为反长度,表示在单位长度内波的相位变化多少。波数可以理解为波的“波长”的倒数。波数越大，波长越短，波的振动周期越密集。

                具体来说：

                ：描述了波在 方向上的传播特性。越大，意味着波在方向上变化得越快，即波的空间振动频率越高。

                ：描述了波在方向上的传播特性。与类似，越大，波在方向上的振动越密集

•         是波的角频率，角频率是波的振动频率，它描述了波在时间上的振荡速度。角频率的单位是弧度每秒（rad/s），它与通常使用的频率  的关系为：，其中的单位是赫兹（Hz）。角频率表示波在单位时间内经过的相位变化量。角频率越高，波在时间上的变化越快，振动周期越短。换句话说，高频率波的振动次数更多，其时间上的波动更为频繁。

        是虚数单位。

        这个方程描述了波随着时间和空间位置 的传播，频率和波数决定了波的空间和时间变化。

 1.1.2 色散关系

        波数 ​ 和角频率之间的关系由色散方程决定：

         是波速；

        色散关系表明，频率 和波数 ​ ​ 不是独立的。这里 被选择为依赖变量，可以通过  ​ 和 计算得到。

1.1.3 波动方程的一般解

        波动方程的最一般解是所有可能解的组合：

        其中：