1.一种基于矩量法的快速分析磁共振射频线圈的方法,包括:
获取输入的磁共振射频线圈的参数信息,根据所述磁共振射频线圈的参数信息建立八叉树数据结构; 选取所述八叉树数据结构的基准层,获取所述基准层的第一阻抗矩阵,所述第一阻抗矩阵存储的数据为所述基准层中非空子组之间的远场弱相互作用对应的阻抗矩阵子块; 对所述基准层的第一阻抗矩阵进行奇异值分解,得到所述基准层的基矩阵和耦合矩阵; 根据所述基准层的基矩阵和耦合矩阵,生成所述基准层的第二阻抗矩阵; 将基准层作为目标层,获取所述目标层的父层,所述父层为与所述目标层邻接的上层,若所述父层不为所述八叉树数据结构的最大层,则根据所述目标层的基矩阵与第一阻抗矩阵,得到所述父层的基矩阵; 获取所述父层的第一阻抗矩阵,所述父层的第一阻抗矩阵存储的数据为所述父层中非空子组之间的远场弱相互作用对应的阻抗矩阵子块,根据所述父层第一阻抗矩阵和基矩阵生成所述父层的耦合矩阵; 根据所述父层的基矩阵和耦合矩阵,生成所述父层的第二阻抗矩阵,并将所述父层作为目标层,迭代执行所述获取所述目标层的父层的步骤。
2.根据权利要求1所述的方法,其特征在于,所述建立八叉树数据结构的步骤之后,还包括: 获取输入的电流和/或电压源和所述八叉树数据结构各层的近场作用阻抗矩阵,计算得到所述八叉树数据结构各层的近场作用激励系数; 根据所述八叉树数据结构各层的近场作用激励系数,计算八叉树数据结构的近场作用总激励系数; 所述得到所述基准层的基矩阵和耦合矩阵的步骤之后,还包括:
在基准层的非空子组中选取任一非空子组作为特定非空子组,获取输入的电压和/或电流源和所述基准层的特定非空子组的基矩阵,计算得到所述基准层的特定非空子组的远场作用激励系数; 根据所述基准层的耦合矩阵,计算得到所述基准层的远场作用激励系数集合,根据所述基准层的远场作用激励系数集合,计算得到所述基准层的远场作用激励系数; 所述将基准层作为目标层,获取所述目标层的父层的步骤之后,还包括: 根据所述目标层的基矩阵与第一阻抗矩阵得到所述目标层的转移矩阵; 根据所述目标层的特定非空子组的远场作用激励系数及所述目标层的转移矩阵,计算得到所述目标层的特定非空子组在所述父层上对应的非空子组的远场作用激励系数; 所述生成所述父层的耦合矩阵的步骤之后,还包括: 根据父层上对应的非空子组的远场作用激励系数和所述父层的耦合矩阵,计算得到所述父层的远场作用激励系数集合,根据所述父层的远场作用激励系数集合,计算得到所述父层的远场作用激励系数; 所述方法还包括:
根据所述八叉树数据结构中各层的远场作用激励系数,计算八叉树数据结构的远场作用总激励系数; 根据所述八叉树数据结构远场作用总激励系数和近场作用总激励系数,计算得到所述磁共振射频线圈的激励系数。
3.根据权利要求1所述的方法,其特征在于,所述根据所述磁共振射频线圈的参数信息建立八叉树数据结构的步骤包括: 建立立方体空间坐标模型,所述立方体空间坐标模型为完全包围所述磁共振射频线圈的最小立方体空间; 将所述立方体空间坐标模型分割为八个大小相等的子立方体空间坐标模型; 根据所述磁共振射频线圈的参数信息对所述子立方体空间坐标模型进行迭代分割,建立八叉树数据结构。
4.根据权利要求1所述的方法,其特征在于,所述选取所述八叉树数据结构的基准层,获取所述基准层的第一阻抗矩阵的步骤还包括: 选取所述八叉树数据结构的基准层;
获取所述基准层的特定非空子组及与其通过远场弱相互作用对应的非空子组集合; 根据所述特定非空子组及与所述特定非空子组对应的远场作用非空子组集合,获取所述基准层的第一阻抗矩阵。
5.根据权利要求1所述的方法,其特征在于,所述对所述基准层的第一阻抗矩阵进行奇异值分解,得到所述基准层的基矩阵和耦合矩阵的步骤包括: 对所述基准层的第一阻抗矩阵进行奇异值分解,得到所述基准层的第一左酉矩阵、第一奇异值矩阵及第一右酉矩阵; 将所述第一左酉矩阵与所述第一奇异值矩阵相乘并按行进行排列,得到所述基准层的中间矩阵; 根据预设的截断误差,对所述中间矩阵进行奇异值分解,得到所述基准层的基矩阵、第二奇异值矩阵及第二右酉矩阵; 根据所述基准层的基矩阵,生成所述基准层的耦合矩阵。
6.根据权利要求5所述的方法,其特征在于,所述根据所述基准层的基矩阵,生成所述基准层的耦合矩阵的步骤包括: 根据公式:
生成所述基准层的耦合矩阵;其中,L表示所述基准层,i表示所述基准层的非空子组,j表示所述非空子组i对应的远场作用非空子组集合中的任一非空子组,集合 表示所述基准层的第二奇异值矩阵, 表示所述基准层的第二右酉矩阵 的共轭转置, 表示所述基准层的第一右酉矩阵, 表示所述基准层的基矩阵 的共轭, 表示所述基准层的耦合矩阵。
7.根据权利要求1所述的方法,其特征在于,所述根据所述基准层的基矩阵和耦合矩阵,生成所述基准层的第二阻抗矩阵的步骤包括: 根据公式:
生成所述基准层的第二阻抗矩阵;其中,Z(L)表示所述基准层的第二阻抗 矩阵,表示所述基准层所有的基矩阵 组成的对角矩阵, 表示所述基准层的耦合矩阵, 表示所述基准层所有的基矩阵组成 的对角矩阵的转置。
8.根据权利要求1所述的方法,其特征在于,所述根据所述目标层的基矩阵与第一阻抗矩阵,得到所述父层的基矩阵的步骤包括: 将所述目标层的基矩阵与所述目标层的第一阻抗矩阵的第一左酉矩阵相乘,得到所述目标层的转移矩阵; 将所述目标层的转移矩阵与所述目标层的基矩阵相乘,得到所述父层的基矩阵。
9.根据权利要求1所述的方法,其特征在于,所述获取所述父层的第一阻抗矩阵,所述父层的第一阻抗矩阵存储的数据为所述父层中非空子组之间的远场弱相互作用对应的阻抗矩阵子块,根据所述父层第一阻抗矩阵和基矩阵生成所述父层的耦合矩阵的步骤包括: 根据公式: 生成所述父层的耦合矩阵;其中,L-1表示所述父层, 表示所述父层的基矩阵, 表示 的共轭转置, 表示所述父层的第一阻抗矩阵。
10.根据权利要求1所述的方法,其特征在于,所述根据所述父层的基矩阵和耦合矩阵,生成所述父层的第二阻抗矩阵步骤包括: 根据公式:
生成所述父层的第二阻抗矩阵;其中,Z(L-1)表示所述父层的第二阻抗矩阵,表示所述父层所有的基矩阵 组成的对角矩阵, 表示所述父层的耦合矩阵, 表示所述父层所有的基矩阵组成 的对角矩阵的转置。
11.一种基于矩量法的快速分析磁共振射频线圈的系统,包括:
建模模块,用于获取输入的磁共振射频线圈的参数信息,根据所述磁共振 射频线圈的参数信息建立八叉树数据结构; 基准层第一阻抗矩阵获取模块,用于选取所述八叉树数据结构的基准层,获取所述基准层的第一阻抗矩阵,所述第一阻抗矩阵存储的数据为所述基准层中非空子组之间的远场弱相互作用对应的阻抗矩阵子块; 基准层计算模块,用于对所述基准层的第一阻抗矩阵进行奇异值分解,得到所述基准层的基矩阵和耦合矩阵; 基准层第二阻抗矩阵生成模块,用于根据所述基准层的基矩阵和耦合矩阵,生成所述基准层的第二阻抗矩阵; 父层第一计算模块,用于将基准层作为目标层,获取所述目标层的父层,所述父层为与所述目标层邻接的上层,若所述父层不为所述八叉树数据结构的最大层,则根据所述目标层的基矩阵与第一阻抗矩阵,得到所述父层的基矩阵; 父层第二计算模块,用于获取所述父层的第一阻抗矩阵,所述父层的第一阻抗矩阵存储的数据为所述父层中非空子组之间的远场弱相互作用对应的阻抗矩阵子块,根据所述父层第一阻抗矩阵和基矩阵生成所述父层的耦合矩阵; 父层第二阻抗矩阵生成模块,用于根据所述父层的基矩阵和耦合矩阵,生成所述父层的第二阻抗矩阵; 迭代模块,用于将所述父层作为目标层,迭代执行所述获取所述目标层的父层。
12.根据权利要求11所述的的系统,其特征在于,所述系统还包括: 基准层特定非空子组的远场作用激励系数获取模块,用于在基准层的非空子组中选取任一非空子组作为特定非空子组,获取输入的电压和/或电流源和基准层的特定非空子组的基矩阵,计算得到所述基准层的特定非空子组的远场作用激励系数; 父层对应的非空子组的远场作用激励系数获取模块,用于将所述基准层作为目标层,获取所述目标层的父层,所述父层为与所述目标层邻接的上层,若获取到,则根据所述目标层的特定非空子组的远场作用激励系数及所述目标层的转移矩阵,计算得到所述目标层的特定非空子组在所述父层上对应的非空子组的远场作用激励系数; 八叉树对应的非空子组的远场作用激励系数获取模块,用于将所述父层作为目标层,迭代执行所述获取所述目标层的父层的步骤,得到所述目标层的特定非空子组在八叉树数据结构各层中对应的非空子组的远场作用激励系数; 远场作用总激励系数获取模块,用于根据所述八叉树数据结构中各层的耦合矩阵,计算得到所述八叉树数据结构中各层的远场作用激励系数集合,根据所述各层的远场作用激励系数集合,计算得到各层的远场作用激励系数,根据所述各层的远场作用激励系数,计算得到八叉树数据结构的远场作用总激励系数; 近场作用总激励系数获取模块,用于获取所述八叉树数据结构中各层的近场作用阻抗矩阵,根据输入的电流和/或电压源,计算得到所述八叉树数据结构中各层的近场作用激励系数,根据所述各层的近场作用激励系数,计算得到八叉树数据结构的近场作用总激励系数; 磁共振射频线圈的激励系数获取模块,用于根据所述八叉树数据结构的远场作用总激励系数和近场作用总激励系数,计算得到所述磁共振射频线圈的激励系数。
13.根据权利要求11所述的系统,其特征在于,所述建模模块包括: 立方体空间坐标模型建立模块,用于建立立方体空间坐标模型,所述立方体空间坐标模型为完全包围所述磁共振射频线圈的最小立方体空间; 分割模块,用于将所述立方体空间坐标模型分割为八个大小相等的子立方体空间坐标模型; 八叉树建立模块,用于根据所述磁共振射频线圈的参数信息对所述子立方体空间坐标模型进行迭代分割,建立八叉树数据结构。
14.根据权利要求11所述的系统,其特征在于,所述基准层第一阻抗矩阵获取模块包括: 基准层选取模块,用于选取所述八叉树数据结构的基准层;
特定非空子组获取模块,用于获取所述基准层的特定非空子组及与其通过远场弱相互作用对应的非空子组集合; 参数获取模块,用于根据所述特定非空子组及与所述特定非空子组对应的 远场作用非空子组集合,获取所述基准层的第一阻抗矩阵。
15.根据权利要求11所述的系统,其特征在于,基准层计算模块包括: 第一奇异值分解模块,用于对所述基准层的第一阻抗矩阵进行奇异值分解,得到所述基准层的第一左酉矩阵、第一奇异值矩阵及第一右酉矩阵; 中间矩阵获取模块,将所述第一左酉矩阵与所述第一奇异值矩阵相乘并按行进行排列,得到所述基准层的中间矩阵; 第二奇异值分解模块,用于根据预设的截断误差,对所述中间矩阵进行奇异值分解,得到所述基准层的基矩阵、第二奇异值矩阵及第二右酉矩阵; 耦合矩阵生成模块,用于根据所述基准层的基矩阵,生成所述基准层的耦合矩阵。
16.根据权利要求15所述的系统,其特征在于,所述耦合矩阵生成模块根据公式: 生成所述基准层的耦合矩阵;其中,L表示所述基准层,i表示所述基准层的非空子组,j表示所述非空子组i对应的远场作用非空子组集合中的任一非空子组,集合 表示所述基准层的第二奇异值矩阵, 表示所述基准层的第二右酉矩阵 的共轭转置, 表示所述基准层的第一右酉矩阵, 表示所述基准层的基矩阵 的共轭, 表示所述基准层的耦合矩阵。
17.根据权利要求11所述的系统,其特征在于,所述基准层第二阻抗矩阵生成模块根据公式: 生成所述基准层的第二阻抗矩阵;其中,Z(L)表示所述基准层的第二阻抗矩阵,表示所述基准层所有的基矩阵 组成的对角矩阵, 表示所述基准层的耦合矩阵, 表示所述基准层所有的基矩阵组成 的对角矩阵的转置。
18.根据权利要求11所述的系统,其特征在于,所述父层第一计算模块包括: 转移矩阵获取模块,用于将所述目标层的基矩阵与所述目标层的第一阻抗矩阵的第一左酉矩阵相乘,得到所述目标层的转移矩阵; 父层基矩阵获取模块,用于将所述目标层的转移矩阵与所述目标层的基矩阵相乘,得到所述父层的基矩阵。
19.根据权利要求11所述的系统,其特征在于,所述父层第二计算模块根据公式: 生成所述父层的耦合矩阵;其中,L-1表示所述父层, 和 表示所述父层的基矩阵, 表示 的共轭转置, 表示所述父层的第一阻抗矩阵。
20.根据权利要求11所述的系统,其特征在于,所述父层第二阻抗矩阵生成模块根据公式: 生成所述父层的第二阻抗矩阵;其中,Z(L-1)表示所述父层的第二阻抗矩阵,表示所述父层所有的基矩阵 组成的对角矩阵, 表示所述父层的耦合矩阵, 表示所述父层所有的基矩阵组成 的对角矩阵的转置。