基于傅里叶变换的磁化率张 tensor 快速磁场扰动计算法：利用张量形式实现麦克斯韦方程组傅里叶变换的 MATLAB 算法

简介：
本研究提出了一种高效的MATLAB算法，通过傅里叶变换和张量方法快速求解磁化率张量在复杂磁场扰动下的响应，基于麦克斯韦方程组。 该算法采用张量形式实现麦克斯韦方程的傅里叶变换，能够从前向计算已知磁化率分布的三维圆柱体中的感应场。在MRI中，通过傅里叶变换可以将麦克斯韦方程组（Marques 和 Bowtell, 2005 年和 Salomir 等人, 2003 年）应用于从磁化率到磁场图的转换。尽管人体组织中的大多数结构具有各向同性的磁化率，但一些生物分子如大脑中髓鞘被认为表现出磁化率各向异性，并可以用二阶张量来描述这种特性。适应了这些磁化率张量公式的麦克斯韦方程傅里叶变换在Liu 2010年的研究（关于磁化率张量成像）的公式9中导出。 此例程将磁化率的张量形式融入到三维实心圆柱体场畸变前向计算之中。函数 `[A]=sus_tensor(PHI,THETA,B0,sus, AA,BB,CC)` 的输入参数包括在适用球坐标下的 `PHI` 和 `THETA`，外加磁场与圆柱之间的角度范围为 0 至 π/2 ，以及磁化率和其它相关的张量系数。