CN114217254A - 一种高线性度梯度线圈设计方法 - Google Patents

一种高线性度梯度线圈设计方法 Download PDF

Info

Publication number
CN114217254A
CN114217254A CN202111521122.1A CN202111521122A CN114217254A CN 114217254 A CN114217254 A CN 114217254A CN 202111521122 A CN202111521122 A CN 202111521122A CN 114217254 A CN114217254 A CN 114217254A
Authority
CN
China
Prior art keywords
current density
gradient coil
function
magnetic field
current
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111521122.1A
Other languages
English (en)
Inventor
周向阳
王维乾
赵风文
谢晓旋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202111521122.1A priority Critical patent/CN114217254A/zh
Publication of CN114217254A publication Critical patent/CN114217254A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开了一种高线性度梯度线圈设计方法，提出了两个新的傅里叶级数展开的电流密度函数，引入流函数；选定N个目标区域内的离散目标点，依据毕奥‑萨伐尔定律推导目标场点的磁场表达式，与电流密度函数之间构建关系矩阵；通过设定目标场点的取值，选取流函数的最小曲率作为约束线圈的罚函数，设定罚函数的约束因子，通过正则化方法求解电流密度表达式中的待定系数，得到圆柱表面的电流密度分布情况；使用流函数对电流密度进行离散，生成线圈的绕组图。本发明通过改进目标场法设计梯度线圈，极大的提高了梯度场的线性度；并且在计算前预先设定了圆柱线圈的尺寸，克服了对尺寸的限制；引入曲率约束和正则化方法解决解的病态问题，使梯度线圈便于加工制作。

Description

一种高线性度梯度线圈设计方法
技术领域
本发明涉及一种高线性度梯度线圈设计方法，设计的径向梯度线圈线性度误差小于0.9％，轴向梯度线圈线性度误差小于0.5％，可应用于核磁共振成像系统、SERF原子磁强计测量系统等需要高线性度梯度磁场的领域。
背景技术
梯度线圈是磁共振系统最重要的部件之一，梯度线圈加电流时可以在成像区域产生梯度变化的磁场。梯度线圈的作用主要是对空间进行编码，区分不同位置的核磁共振信号。此外，在基于SERF原子磁强计测量成像领域，为实现原子自旋的SERF态，常采用被动磁屏蔽包括磁屏蔽桶或磁屏蔽房，使其免受电磁噪声和地磁场的干扰。然而在实际应用中，磁屏蔽桶的一端必须打开，外部磁场进入磁屏桶中，因此需要主动磁补偿减少内部剩余磁场，需要匀场线圈与梯度线圈，补偿轴向磁场与梯度磁场，实现测量所需要的高稳定大均匀区的磁场环境。
梯度线圈设计问题属于电磁场逆问题，主要思想是根据给定的磁场分布，确定能够产生相应的磁场且具有一定结构的梯度线圈。设计方法可分为两类：分离绕线法和分布绕线法。分离绕线法是先选择线圈的基本结构，然后通过解析方法或数值优化方法调整线圈的尺寸和位置。分离绕线法直接，结构简单，但是很难找到满足条件的参数，计算复杂。分布绕线法需要计算连续的电流分布，然后将电流密度离散，得到线圈的具体绕线形式。目标场法是分布绕线方法中最经典的方法，然而原始的目标场法假设线圈无限长，需要额外的函数来限制线圈的固定长度；将电流密度展开成傅里叶级数，在求解过程中会产生病态问题，因此要添加正则化约束，解决积分方程病态问题。在磁屏蔽桶中，大多只考虑了采用鞍形线圈和亥姆霍兹线圈，产生匀强磁场，但是这样距离零磁环境还有一定距离，因此需要找到一种高线性度梯度线圈设计方法，设计高线性度梯度线圈应用于磁屏蔽桶中。
发明内容
本发明的技术解决问题：克服现有技术的不足，提供一种高线性度梯度线圈的设计方法，解决磁屏蔽桶单端开口轴向磁场与梯度磁场较大的问题，通过改进目标场法设计梯度线圈，极大的提高了梯度场的线性度；并且在计算前预先设定了圆柱线圈的尺寸，克服了对尺寸的限制；引入曲率约束和正则化方法解决解的病态问题，使梯度线圈便于加工制作。
本发明的技术解决方案是：一种高线性度梯度线圈设计方法，包括以下步骤：
(1)预先确定梯度线圈包括径向梯度线圈与轴向梯度线圈所在圆柱骨架的尺寸、梯度线圈的目标区域范围，给出圆柱表面上连续分布的面电流的电流密度表达式；
(2)构建双傅里叶级数用于限制圆柱表面上连续分布的面电流的电流密度得到新的径向梯度线圈与轴向梯度线圈的电流密度，并引入流函数；
(3)由毕奥-萨伐尔定律确定空间中选定的目标区域内N个离散目标场点的磁场表达式；
(4)设定梯度线圈目标区域N个场点的目标磁场值，对流函数施加最小曲率约束条件作为梯度线圈罚函数，选定惩罚因子，通过Tikhonov正则化方法求解电流密度表达式中待定系数，得到圆柱表面的电流密度分布情况；
(5)基于电流连续性方程，采用流函数对得到的电流密度进行离散化处理，得到梯度线圈的实际电流走线形状即梯度线圈的绕组图，满足均匀性要求，从而完成梯度线圈的设计工作。
所述步骤(1)具体实现如下：
步骤(11)设定梯度线圈所在圆柱骨架尺寸，即半径为R和高度为2H；确定梯度线圈的目标区域范围为与圆柱线圈骨架同几何中心的球体，另外，目标区域代形状也可以是立方体或圆柱体；
步骤(12)圆柱表面上连续分布的面电流的电流密度J(r′)沿径向电流密度分量为Jρ，沿轴向电流密度分量为
Figure BDA0003408577700000031
在圆柱坐标系下，径向电流密度分量为零，产生梯度场，圆柱表面上连续分布的面电流的电流密度J(r′)是方位角θ与高度z的函数，J(r′)表达式为：
J(r′)＝Jθ(θ,z)eθ+Jz(θ,z)ez
其中，r′为圆柱表面上任一坐标点(R,θ,z)，r为待求磁场坐标点
Figure BDA0003408577700000032
Jθ(θ,z)为沿着方位角的电流密度分量，Jz(θ,z)为沿着高度方向的电流密度分量；eθ与ez为沿着方位角与沿着高度方向的单位向量。
所述步骤(2)具体实现如下：
步骤(21)利用三角基函数，将电流密度展开成傅里叶级数，从而将梯度线圈上的电流分布限制在有限半径内，并且范围以外的电流密度都为零；
步骤(22)对于径向梯度线圈的电流密度J(r′)表示为：
Figure BDA0003408577700000033
对于轴向梯度线圈，提出一种新的电流密度函数，电流密度J(r′)表示为：
Figure BDA0003408577700000034
其中，Pm为各阶三角函数的展开系数即电流密度系数，M为展开式阶数；
步骤(23)电流密度满足稳流条件：
Figure BDA0003408577700000035
在线圈上引入流函数I(θ,z)则其满足：
Figure BDA0003408577700000041
根据流函数满足的关系式以及径向梯度线圈的电流密度与轴向梯度线圈的电流密度，得流函数在径向和轴向的表达式为：
Figure BDA0003408577700000042
Figure BDA0003408577700000043
所述步骤(3)具体实现如下:
步骤(31)根据毕奥-萨伐尔定律，电流密度J(r′)在任意点r′处的感应磁场B可表示为：
Figure BDA0003408577700000044
其中，μ0为真空磁导率，S为整个圆柱表面积；
步骤(32)在目标区域选取N个离散分布的目标场点，然后将步骤(22)中得到的径向梯度线圈和轴向梯度线圈的新电流密度傅里叶级数展开式，分别代入步骤(31)感应磁场表达式中，得到目标场点所产生的感应磁场在x轴和z轴分量的表达式：
Figure BDA0003408577700000045
Figure BDA0003408577700000046
其中，定义横向线圈的电流密度傅里叶级数代入后关于
Figure BDA0003408577700000047
的函数为
Figure BDA0003408577700000048
纵向线圈的电流密度傅里叶级数代入后关于
Figure BDA0003408577700000049
的函数为
Figure BDA00034085777000000410
所述步骤(4)具体实现如下：
步骤(41)设定目标区域N个目标场点的目标磁场值，得到电流密度级数的系数，目标点的感应磁场表达式的矩阵形式：
Figure BDA0003408577700000051
其中，B′v,i为沿z轴N个磁场点目标磁场值；
步骤(42)求解目标点的感应磁场表达式矩阵时，有多个满足条件的数值解，这是高度病态。因此，利用Tikhonov正则化方法修正待定系数求解所述矩阵中的元素，克服待定系数求解的病态问题，首先，定义计算磁场与目标磁场之差的平方和为误差函数E，表达式为：
Figure BDA0003408577700000052
然后引入罚函数F和罚函数因子λ，对于误差函数E施加限定约束，使用称为残余误差的新误差E′代替E，将误差函数修改为E′：
Figure BDA0003408577700000053
步骤(43)为了便于优化梯度线圈导线的平滑程度，使线圈电流走线形状规则，选择流函数的最小曲率为罚函数，在x轴和z轴数学定义为：
Figure BDA0003408577700000054
将步骤(23)中的Ix(θ,z)与Iz(θ,z)代入罚函数数学定义中，得到惩罚函数的数学表达式：
Figure BDA0003408577700000055
对步骤(42)中误差函数E′各项系数求导，将Pm的取值设为使得这些导数为零的值，得到E′最小值，该过程表示为：
Figure BDA0003408577700000061
即E′＝0可以求解出电流密度表达式的待定系数。
所述步骤(5)具体实现如下：
步骤(51)将电流密度系数Pm代入到流函数I(θ,z)中；
步骤(52)选取Imax表示平面内流函数的最大值，Imin表示平面内流函数的最小值，梯度线圈离散后，影响电流匝数的参数设为N，获得I(θ,z)的等值线为：
{I(θ,z)＝Imin+(i+1/2)I0,i＝0,1,2,…,N-1}
其中，I0＝Imax-Imin/N，得到的流函数等值线分布代表了梯度线圈离散化后导线的分布情况，每匝导线通入的电流均为I0；
步骤(53)依据梯度线圈离散后导线的分布情况，计算流函数离散后的实际磁场均匀度。其中，磁场均匀度为任意目标点感应磁场大小与中心感应磁场的差与中心磁场的比值的绝对值，若离散后实际磁场均匀度不满足均匀性要求，则改变电流匝数的控制参数N，直到满足均匀性要求，则完成梯度线圈的设计。
本发明与现有技术相比的优点在于：
(1)本发明作为梯度线圈设计的逆方法，提出了两个新的傅里叶级数展开的电流密度函数，通过预先限定圆柱线圈的尺寸及选取目标场点的坐标位置，再计算得出梯度线圈的实际电流分布，克服了尺寸的约束，有效的控制了圆柱型线圈的长径比；
(2)本发明引入流函数的曲率作为惩罚函数，采用Tikhonov正则化求解，克服了待定系数求解的病态问题；
(3)本发明改进了目标场法，设计的梯度线圈线性度好，横线梯度线圈的场线性度是Golay线圈的2.5倍，纵向梯度线圈的场线性度是Maxwell线圈的1.2倍。
附图说明
图1为本发明中的圆柱形梯度线圈中场点、源点的坐标关系示意图；
图2为本发明中高线性度梯度线圈设计方法的流程图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
图1为本发明中圆柱形梯度线圈中场点、源点的坐标示意图。圆柱骨架半径为R，高度为2H，r′为圆柱表面上任一坐标点，r为待求磁场坐标点，θ为圆柱坐标系中方位角，TSV表示目标区域。
如图2所示，本发明的一种高线性度梯度线圈方法，具体包括以下步骤：
(一)梯度线圈包括径向梯度线圈与轴向梯度线圈所在圆柱骨架的尺寸，目标区域以及所述圆柱表面的电流密度表达式确定的过程，该过程分为以下2个步骤：
步骤(1)设定梯度线圈所在圆柱骨架尺寸，即半径为R和高度为2H；确定梯度线圈的目标区域范围为与圆柱线圈骨架同几何中心的球体，另外，目标区域代形状也可以是立方体或圆柱体；
步骤(2)圆柱表面连续分布的面电流的电流密度J(r′)沿径向电流密度分量为Jρ，沿轴向电流密度分量为
Figure BDA0003408577700000071
在圆柱坐标系下，径向电流密度分量为零，产生梯度场，圆柱表面上的连续电流密度J(r′)是方位角θ与高度z的函数，因此电流密度J(r′)可以表示为：
J(r′)＝Jθ(θ,z)eθ+Jz(θ,z)ez
其中，r′为圆柱表面上任一坐标点(R,θ,z)，r为待求磁场坐标点
Figure BDA0003408577700000072
Jθ(θ,z)为沿着方位角的电流密度分量，Jz(θ,z)为沿着高度方向的电流密度分量；eθ与ez为沿着方位角与沿着高度方向的单位向量。
(二)构建双傅里叶级数用于限制圆柱表面上连续分布的面电流的电流密度得到新的径向梯度线圈与轴向梯度线圈的电流密度，并引入流函数的过程，该过程分为以下3个步骤：
步骤(1)利用三角基函数，将电流密度展开成傅里叶级数，从而将线圈上的电流分布限制在有限半径内，并且范围以外的电流密度都为零；
步骤(2)对于径向梯度线圈的电流密度J(r′)表示为：
Figure BDA0003408577700000081
对于轴向梯度线圈，提出一种新的电流密度函数，电流密度J(r′)表示为：
Figure BDA0003408577700000082
其中，Pm为各阶三角函数的展开系数即电流密度系数，M为展开式阶数；步骤(3)电流密度满足稳流条件：
Figure BDA0003408577700000083
在线圈上引入流函数I(θ,z)则其满足：
Figure BDA0003408577700000084
根据流函数满足的关系式以及径向梯度线圈的电流密度与轴向梯度线圈的电流密度，得流函数在径向和轴向的表达式为：
Figure BDA0003408577700000085
Figure BDA0003408577700000086
(三)通过毕奥-萨伐尔定律确定空间中选定的目标区域内N个离散目标场点的磁场表达式的过程，该过程分为以下2个步骤：
步骤(1)根据毕奥-萨伐尔定律，在任意点r′处，磁感应矢量B与电流密度J(r′)可表示为：
Figure BDA0003408577700000087
其中，μ0为真空磁导率，S为整个圆柱表面积。
步骤(2)在目标区域选取N个离散分布的目标场点，然后将(二)步骤(2)中得到的径向梯度线圈和轴向梯度线圈的新电流密度傅里叶级数展开式，分别代入(三)步骤(1)感应磁场表达式中，得到目标场点所产生的感应磁场在x轴和z轴分量的表达式：
Figure BDA0003408577700000091
Figure BDA0003408577700000092
其中，定义横向线圈的电流密度傅里叶级数代入后关于
Figure BDA0003408577700000093
的函数为
Figure BDA0003408577700000094
纵向线圈的电流密度傅里叶级数代入后关于
Figure BDA0003408577700000095
的函数为
Figure BDA0003408577700000096
(四)设定梯度线圈目标区域N个场点的目标磁场值，对流函数施加最小曲率约束条件作为梯度线圈罚函数，选定惩罚因子，通过Tikhonov正则化方法求解电流密度表达式中待定系数，得到圆柱表面的电流密度分布情况的过程，该过程分为以下3个步骤：
步骤(1)设定目标区域N个目标场点的目标磁场值，得到电流密度级数的系数，目标点的感应磁场表达式可以写成矩阵的形式：
Figure BDA0003408577700000097
其中，B′v,i为沿z轴N个磁场点目标磁场值。
步骤(2)求解目标点的感应磁场表达式矩阵时，有多个满足条件的数值解，这是高度病态。因此，利用Tikhonov正则化方法修正待定系数求解所述矩阵中的元素，克服待定系数求解的病态问题，首先，定义计算磁场与目标磁场之差的平方和为误差函数E，表达式为：
Figure BDA0003408577700000098
然后引入罚函数F和罚函数因子λ，对于误差函数E施加限定约束，使用称为残余误差的新误差E′代替E，将误差函数修改为：
Figure BDA0003408577700000101
步骤(3)为了便于优化梯度线圈导线的平滑程度，使线圈电流走线形状规则，选择流函数的最小曲率为罚函数，在x轴和z轴数学定义为：
Figure BDA0003408577700000102
将(二)步骤(3)中的Ix(θ,z)与Iz(θ,z)代入上述表达式，得到惩罚函数的数学表达式：
Figure BDA0003408577700000103
对(四)步骤(2)中误差函数E′各项系数求导，将Pm的取值设为使得这些导数为零的值，得到E′最小值，该过程表示为：
Figure BDA0003408577700000104
即E′＝0可以求解出电流密度表达式的待定系数。
(五)基于电流连续性方程，采用流函数对得到的电流密度进行离散化处理，得到梯度线圈的实际电流走线形状即梯度线圈的绕组图，满足均匀性要求，从而完成梯度线圈设计工作的过程，该过程分为以下3个步骤：
步骤(1)将电流密度系数Pm代入到流函数I(θ,z)中；
步骤(2)选取Imax表示平面内流函数的最大值，Imin表示平面内流函数的最小值，将梯度线圈离散后，影响电流匝数的参数设为N，获得I(θ,z)的等值线为：
{I(θ,z)＝Imin+(i+1/2)I0,i＝0,1,2,…,N-1}
其中，I0＝Imax-Imin/N，得到的流函数等值线分布代表了线圈离散化后导线的分布情况，每匝导线通入的电流均为I0。
步骤(3)依据梯度线圈离散后导线的分布情况，计算流函数离散后的实际磁场均匀度。其中，磁场均匀度为任意目标点感应磁场大小与中心感应磁场的差与中心磁场的比值的绝对值，若离散后实际磁场均匀度不满足均匀性要求，则改变电流匝数的控制参数N，直到满足均匀性要求，则完成梯度线圈的设计。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (6)

1.一种高线性度梯度线圈设计方法，其特征在于，包括以下步骤：
(1)预先确定梯度线圈包括径向梯度线圈与轴向梯度线圈所在圆柱骨架的尺寸、梯度线圈的目标区域范围，给出圆柱表面上连续分布的面电流的电流密度表达式；
(2)构建双傅里叶级数用于限制圆柱表面上连续分布的面电流的电流密度得到新的径向梯度线圈与轴向梯度线圈的电流密度，并引入流函数；
(3)由毕奥-萨伐尔定律确定空间中选定的目标区域内N个离散目标场点的磁场表达式；
(4)设定梯度线圈目标区域N个场点的目标磁场值，对流函数施加最小曲率约束条件作为梯度线圈罚函数，选定惩罚因子，通过Tikhonov正则化方法求解电流密度表达式中待定系数，得到圆柱表面的电流密度分布情况；
(5)基于电流连续性方程，采用流函数对得到的电流密度进行离散化处理，得到梯度线圈的实际电流走线形状即梯度线圈的绕组图，满足均匀性要求，从而完成梯度线圈的设计工作。
2.根据权利要求1所述的高线性度梯度线圈设计方法，其特征在于：所述步骤(1)具体实现如下：
步骤(11)设定梯度线圈所在圆柱骨架尺寸，即半径为R和高度为2H；确定梯度线圈的目标区域范围为与圆柱线圈骨架同几何中心的球体，另外，目标区域代形状也可以是立方体或圆柱体；
步骤(12)圆柱表面上连续分布的面电流的电流密度J(r′)沿径向电流密度分量为Jρ，沿轴向电流密度分量为
Figure FDA0003408577690000011
在圆柱坐标系下，径向电流密度分量为零，产生梯度场，圆柱表面上连续分布的面电流的电流密度J(r′)是方位角θ与高度z的函数，J(r′)表达式为：
J(r′)＝Jθ(θ,z)eθ+Jz(θ,z)ez
其中，r′为圆柱表面上任一坐标点(R,θ,z)，r为待求磁场坐标点
Figure FDA0003408577690000021
Jθ(θ,z)为沿着方位角的电流密度分量，Jz(θ,z)为沿着高度方向的电流密度分量；eθ与ez为沿着方位角与沿着高度方向的单位向量。
3.根据权利要求1所述的高线性度梯度线圈设计方法，其特征在于，所述步骤(2)具体实现如下：
步骤(21)利用三角基函数，将电流密度展开成傅里叶级数，从而将梯度线圈上的电流分布限制在有限半径内，并且范围以外的电流密度都为零；
步骤(22)对于径向梯度线圈的电流密度J(r′)表示为：
Figure FDA0003408577690000022
对于轴向梯度线圈，提出一种新的电流密度函数，电流密度J(r′)表示为：
Figure FDA0003408577690000023
其中，Pm为各阶三角函数的展开系数即电流密度系数，M为展开式阶数；
步骤(23)电流密度满足稳流条件：
Figure FDA0003408577690000024
在线圈上引入流函数I(θ,z)则其满足：
Figure FDA0003408577690000025
根据流函数满足的关系式以及径向梯度线圈的电流密度与轴向梯度线圈的电流密度，得流函数在径向和轴向的表达式为：
Figure FDA0003408577690000026
Figure FDA0003408577690000031
4.根据权利要求1所述的高线性度梯度线圈设计方法，其特征在于，所述步骤(3)具体实现如下：
步骤(31)根据毕奥-萨伐尔定律，电流密度J(r′)在任意点r′处的感应磁场B可表示为：
Figure FDA0003408577690000032
其中，μ0为真空磁导率，S为整个圆柱表面积；
步骤(32)在目标区域选取N个离散分布的目标场点，然后将步骤(22)中得到的径向梯度线圈和轴向梯度线圈的新电流密度傅里叶级数展开式，分别代入步骤(31)感应磁场表达式中，得到目标场点所产生的感应磁场在x轴和z轴分量的表达式：
Figure FDA0003408577690000033
Figure FDA0003408577690000034
其中，定义横向线圈的电流密度傅里叶级数代入后关于
Figure FDA0003408577690000035
的函数为
Figure FDA0003408577690000036
纵向线圈的电流密度傅里叶级数代入后关于
Figure FDA0003408577690000037
的函数为
Figure FDA0003408577690000038
5.根据权利要求3所述的高线性度梯度线圈设计方法，其特征在于，所述步骤(4)具体实现如下：
步骤(41)设定目标区域N个目标场点的目标磁场值，得到电流密度级数的系数，目标点的感应磁场表达式的矩阵形式：
Figure FDA0003408577690000039
其中，Bv′,i为沿z轴N个磁场点目标磁场值；
步骤(42)求解目标点的感应磁场表达式矩阵时，有多个满足条件的数值解，这为高度病态，因此利用Tikhonov正则化方法修正待定系数求解所述矩阵中的元素，克服待定系数求解的高病态问题，首先，定义计算磁场与目标磁场之差的平方和为误差函数E，表达式为：
Figure FDA0003408577690000041
然后引入罚函数F和罚函数因子λ，对于误差函数E施加限定约束，使用称为残余误差的新误差E′代替E，将误差函数修改为E′：
Figure FDA0003408577690000042
步骤(43)为了便于优化梯度线圈导线的平滑程度，使线圈电流走线形状规则，选择流函数的最小曲率为罚函数，在x轴和z轴数学定义为：
Figure FDA0003408577690000043
将步骤(23)中的Ix(θ,z)与Iz(θ,z)代入罚函数数学定义中，得到惩罚函数的数学表达式：
Figure FDA0003408577690000044
对步骤(42)中误差函数E′各项系数求导，将Pm的取值设为使得这些导数为零的值，得到E′最小值，该过程表示为：
Figure FDA0003408577690000045
即E′＝0求解出电流密度表达式的待定系数。
6.根据权利要求1所述的一种高线性度梯度线圈设计方法，其特征在于，所述步骤(5)具体实现如下：
步骤(51)将电流密度系数Pm代入到流函数I(θ,z)中；
步骤(52)选取Imax表示平面内流函数的最大值，Imin表示平面内流函数的最小值，梯度线圈离散后，影响电流匝数的参数设为N，获得I(θ,z)的等值线为：
{I(θ,z)＝Imin+(i+1/2)I0,i＝0,1,2,…,N-1}
其中，I0＝Imax-Imin/N，得到的流函数等值线分布代表了梯度线圈离散化后导线的分布情况，每匝导线通入的电流均为I0；
步骤(53)依据梯度线圈离散后导线的分布情况，计算流函数离散后的实际磁场均匀度，其中，磁场均匀度为任意目标点感应磁场大小与中心感应磁场的差与中心磁场的比值的绝对值，若离散后实际磁场均匀度不满足均匀性要求，则改变电流匝数的控制参数N，直到满足均匀性要求，则完成梯度线圈的设计。
CN202111521122.1A 2021-12-13 2021-12-13 一种高线性度梯度线圈设计方法 Pending CN114217254A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111521122.1A CN114217254A (zh) 2021-12-13 2021-12-13 一种高线性度梯度线圈设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111521122.1A CN114217254A (zh) 2021-12-13 2021-12-13 一种高线性度梯度线圈设计方法

Publications (1)

Publication Number Publication Date
CN114217254A true CN114217254A (zh) 2022-03-22

Family

ID=80701437

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111521122.1A Pending CN114217254A (zh) 2021-12-13 2021-12-13 一种高线性度梯度线圈设计方法

Country Status (1)

Country Link
CN (1) CN114217254A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116203492A (zh) * 2023-04-28 2023-06-02 山东奥新医疗科技有限公司 一种梯度线圈线性度测试装置

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101996273A (zh) * 2010-11-29 2011-03-30 浙江大学 Mri系统梯度线圈的有限差分设计方法
CN102096053A (zh) * 2010-12-02 2011-06-15 苏州安科医疗系统有限公司 一种核磁共振系统中的梯度线圈设计法
CN104007407A (zh) * 2014-06-13 2014-08-27 河海大学 结合泛函与模拟退火算法的梯度线圈设计方法
CN105718677A (zh) * 2016-01-22 2016-06-29 中国科学院电工研究所 自屏蔽超导核磁共振成像系统梯度线圈设计方法
CN109946630A (zh) * 2019-03-27 2019-06-28 中国计量大学 最小化最大温度值的盘式梯度线圈设计方法
CN110162926A (zh) * 2019-06-05 2019-08-23 江苏麦格迈医学科技有限公司 一种用于低场核磁共振双平面型梯度线圈设计的方法
CN110568390A (zh) * 2019-10-14 2019-12-13 惠仁望都医疗设备科技有限公司 一种双平面匀场线圈的设计方法
CN110703170A (zh) * 2019-09-27 2020-01-17 河北惠仁医疗设备科技有限公司 一种乳腺专用mri梯度线圈的设计方法
CN110780245A (zh) * 2019-11-29 2020-02-11 中国科学院电工研究所 用于平面超导磁共振系统的高屏蔽梯度线圈设计方法及其梯度线圈
CN112597617A (zh) * 2019-09-16 2021-04-02 河海大学 一种梯度场线圈优化方法
CN112684392A (zh) * 2019-10-19 2021-04-20 惠仁望都医疗设备科技有限公司 一种非平面梯度线圈的设计方法
CN113030811A (zh) * 2021-05-19 2021-06-25 中国科学院精密测量科学与技术创新研究院 一种圆筒型匀场线圈的设计方法
CN113033028A (zh) * 2021-05-19 2021-06-25 中国科学院精密测量科学与技术创新研究院 一种磁共振Halbach型匀场线圈的设计方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101996273A (zh) * 2010-11-29 2011-03-30 浙江大学 Mri系统梯度线圈的有限差分设计方法
CN102096053A (zh) * 2010-12-02 2011-06-15 苏州安科医疗系统有限公司 一种核磁共振系统中的梯度线圈设计法
CN104007407A (zh) * 2014-06-13 2014-08-27 河海大学 结合泛函与模拟退火算法的梯度线圈设计方法
CN105718677A (zh) * 2016-01-22 2016-06-29 中国科学院电工研究所 自屏蔽超导核磁共振成像系统梯度线圈设计方法
CN109946630A (zh) * 2019-03-27 2019-06-28 中国计量大学 最小化最大温度值的盘式梯度线圈设计方法
CN110162926A (zh) * 2019-06-05 2019-08-23 江苏麦格迈医学科技有限公司 一种用于低场核磁共振双平面型梯度线圈设计的方法
CN112597617A (zh) * 2019-09-16 2021-04-02 河海大学 一种梯度场线圈优化方法
CN110703170A (zh) * 2019-09-27 2020-01-17 河北惠仁医疗设备科技有限公司 一种乳腺专用mri梯度线圈的设计方法
CN110568390A (zh) * 2019-10-14 2019-12-13 惠仁望都医疗设备科技有限公司 一种双平面匀场线圈的设计方法
CN112684392A (zh) * 2019-10-19 2021-04-20 惠仁望都医疗设备科技有限公司 一种非平面梯度线圈的设计方法
CN110780245A (zh) * 2019-11-29 2020-02-11 中国科学院电工研究所 用于平面超导磁共振系统的高屏蔽梯度线圈设计方法及其梯度线圈
CN113030811A (zh) * 2021-05-19 2021-06-25 中国科学院精密测量科学与技术创新研究院 一种圆筒型匀场线圈的设计方法
CN113033028A (zh) * 2021-05-19 2021-06-25 中国科学院精密测量科学与技术创新研究院 一种磁共振Halbach型匀场线圈的设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
FENGWEN ZHAO: "Design of Highly Linear Gradient Field Coils Based on an Improved Target-Field Method", 《IEEE SENSORS JOURNAL》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116203492A (zh) * 2023-04-28 2023-06-02 山东奥新医疗科技有限公司 一种梯度线圈线性度测试装置
CN116203492B (zh) * 2023-04-28 2023-07-14 山东奥新医疗科技有限公司 一种梯度线圈线性度测试装置

Similar Documents

Publication Publication Date Title
CN109765510B (zh) 一种带有圆角的径向超导匀场线圈设计方法
JP5695241B2 (ja) 傾斜磁場コイルおよび核磁気共鳴撮像装置
CN106990373B (zh) 一种磁共振系统的解耦轴向匀场线圈设计方法
Faugeras An overview of the numerical methods for tokamak plasma equilibrium computation implemented in the NICE code
Tomasi Stream function optimization for gradient coil design
CN112711876B (zh) 一种可降低磁屏蔽耦合效应的偶极均匀磁场线圈设计方法
CN113343551B (zh) 一种基于磁屏蔽筒耦合的圆柱形和球形匀强磁场线圈设计方法
CN108872896A (zh) 一种双平面磁共振成像系统梯度线圈的设计方法
CN104198969B (zh) 一种梯度线圈设计方法
CN110780245B (zh) 用于平面超导磁共振系统的屏蔽梯度线圈设计方法及其梯度线圈
CN102096053A (zh) 一种核磁共振系统中的梯度线圈设计法
CN110568390A (zh) 一种双平面匀场线圈的设计方法
CN114217254A (zh) 一种高线性度梯度线圈设计方法
CN101852843A (zh) 一种超导磁体外磁屏蔽线圈的优化设计算法
CN109407022B (zh) 磁场漂移的控制方法、装置和存储介质
JP4873803B2 (ja) 磁気共鳴のための非対称ゾーナルシムコイル
CN114200368B (zh) 一种磁屏蔽筒内轴向匀场线圈的设计方法
CN112231880A (zh) 一种超导磁体线圈设计方法
CN112444762B (zh) 一种八棱柱形径向匀场线圈设计方法
CN112597617A (zh) 一种梯度场线圈优化方法
US8581586B2 (en) Method of calculating coil pattern, and gradient magnetic field coil
Hu et al. Shim coil set for NMR using a novel target field method based on trigonometric series
Hu et al. Optimization magnetic resonance imaging shim coil using second derivative discretized stream function
Zheng et al. Gradient coils design with regularization method for superconducting magnetic resonance imaging
CN116070297A (zh) 一种圆柱径向匀强磁场线圈设计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20220322

RJ01 Rejection of invention patent application after publication