专利名称：：一种用于磁共振成像超导磁体设计的离散优化方法
技术领域：
：本发明属于应用超导
技术领域：
，尤其涉及磁共振成像(MRI)超导磁体设计。
背景技术：
：磁共振成像(MRI)是一种利用核磁共振现象进行成像的技术。主磁体是MRI系统的重要部件，其作用是在成像区域内产生一个具有特定强度和一定均匀度的背景磁场。其电磁设计指标有-(1)磁场强度^，指成像区域(DSV)中心点磁场值。(2)磁场均匀度，其一种定义方法为一丑血-。—"^"，其中，//o表示磁场均匀度，5dx和5di分别表示成像区域内磁场的最大值和最小值。(3)杂散场范围，通常指0.5mT(豪特斯拉)磁场等值线所占用空间范围。中高场MRI系统中的主磁体通常采用超导磁体。超导磁体是使用超导线材(例如NbTi多芯超导线)在磁体骨架绕制一个或多个超导线圈来实现的。超导磁体充磁后，以恒流模式在低温条件下进行工作。由于超导线成本较高，所以也将超导线用线量作为设计指标。超导磁体的一种通常设计方法是，预先选定工作电流和线圈基本结构，将成像区域内的磁场用谐波方法展开，以线圈尺寸和位置作为自变量，通过主磁场强度和低阶谐波分量抵消关系建立非线性方程组，求解该方程完成设计。另一种常用设计方法是，在特定的线圈结构下，建立优化问题模型，将线圈尺寸和位置作为优化变量，以前述设计指标为目标函数，采用一定的优化算法求解。在这两种方法的最后阶段，由于实际加工的需要，需要考虑超导线截面尺寸，对线圈尺寸进行离散化处理，对位置和半径进行取整。取整和离散后，磁体的均匀度等指标通常会以一定幅度下降，使设计结果偏离最优解。此外，由于实际加工时可能存在误差，在设计最后需要分析可能误差对磁场指标的影响，即进行鲁棒性分析，选择最差情况下仍能满足设计要求的解。但这样的后果是难以实现磁场指标和鲁棒性同时达到最优。本发明中通过采用离散变量优化设计，避免了取整和离散误差，实现了最优设计。通过在优化问题中考虑误差鲁棒性，实现了强鲁棒性的最优设计。
发明内容本发明的目的在于提供一种制磁共振成像(MRI)超导磁体设计方法，其步骤如下步骤(l)，设定如下优化变量，并输入计算机线圈层数和每层匝数，处理为离散变量；线圈半径和位置在不考虑加工精度时，处理为连续变量；在考虑加工精度时，也将线圈半径和位置变量处理为离散变量。所述线圈至少包含主线圈和屏蔽线圈。步骤(2)，按下式计算使用所述计算机采用模拟退火算法进行所述离散优化方法时用的目标函数/，<formula>formulaseeoriginaldocumentpage8</formula>其中，w，^，…^为预先指定的权系数，是设定值，B为中心磁场强度，//为磁场均匀度，S为杂散场范围，Z为超导线用线量，W为制作误差鲁棒性系数。步骤(2.1)，将磁共振成像超导磁体分解为包括电流同向或电流反向的圆柱型线圈在内的线圈组合，使用螺线管内部的Legendre展开法，按下式计算中心磁场强度B和磁场均匀度H:首先，对于每个圆柱型线圈，计算前10阶偶次Legendre系数。以第z个线圈为例，其前10阶偶次Legendre系数M。,，M2,，…，M,。,，如下式所示<formula>formulaseeoriginaldocumentpage8</formula>M4,=*M6,:^^[C/2(8+12C2,+15C2,2_70C2,3+315C2,4)-C3,3/2(8+12C4,+15C4,2-70C4,3+315C4,4)]Ms=~^[C,,2(16+24C2,+30C2,2+35C2,3+315C2,8'396"h-2079C2,5+3003C2,6)-C3,3/2(16+24C4,+30C4,2+35C4,3+315C4,4-2079C4,5+3003C4,6)]M'0'=n^[Ci'3/2(128+192C2'+240C2'2+280C2'3+315C2,4-2772C2,5+42042C2,6-128700C2,7+109395C2,8)-C/2(128+192C4,+240C4,2+280Q,3+315C4,4—2772Q,5+42042C4,6—128700C4i7+109395C4,8)]其中，<formula>formulaseeoriginaldocumentpage9</formula><formula>formulaseeoriginaldocumentpage9</formula>4为第1个线圈的内径，4,为第Z个线圈的外径，A为第l'个线圈的半长。然后，按步骤(2.1)所述方法计算所有圆柱型线圈的所述前10阶中各阶Legendre系数。接着，计算所述磁共振成像超导磁体的各阶Legendre系数M0,M2,…M10。计算时，记各线圈电流密度幅值为lJiI=J，记半径最小的线圈为1号线圈，将其他各线圈的各阶Legendre系数进行电流和内径按其折算求和，如下<formula>formulaseeoriginaldocumentpage9</formula>其中，Mj=0,2,…10为磁体Legendre系数，A11为1号线圈内径，A11为第z个线圈内径，为第i个线圈的7阶Legendre系数，SIGNJ表示第i个线圈电流方向符号，K为线圈个数。再按下式计算中心磁场强度万，再按下式计算磁场均匀度//，其中，^为空气磁导率，J为电流密度幅值，4为最内侧线圈内径，加权系数^入,…,^为设定值。步骤(2.2)，把所述磁共振成像超导磁体分解成包含电流同向或电流反向的圆柱型线圈在内的线圈组合后，按下述方法使用所述螺线管外部的Legendre展开法计算所述杂散场范围首先，对每个圆柱型线圈，计算前6阶奇次Legendre展开系数。以第i个线圈为例，其前6阶奇次Legendre系数A^，A^，，如下式所示iV"(l一",)肌<formula>formulaseeoriginaldocumentpage10</formula>其中，<formula>formulaseeoriginaldocumentpage10</formula>4，4,和丄,含义同前。然后，按步骤(2.2)所述方法计算所有圆柱型线圈的各阶Legendre系数。接着，计算所有圆柱型线圈的各阶Legendre系数的代数和，用A^,A^,A^表示。计算时，记各线圈电流密度幅值为^卜/,记半径最小的线圈为1号线圈，将其他各线圈的各阶Legendre系数进行电流和内径按其折算求和，如下:a:<formula>formulaseeoriginaldocumentpage11</formula>其中，jv"7=1,3,5为磁体Legendre系数，A^,为第！个线圈的_/阶Legendre系数，其他含义同前。再计算杂散场范围S:<formula>formulaseeoriginaldocumentpage11</formula>其中，/,，^/5为加权系数，为设定值。步骤(2.3)，按下式计算超导线圈用线量丄其中，F为线圈总体积，&,为导线截面积。步骤(2.4)，按下式计算磁体鲁棒性系数及。步骤(2.4.1)，按步骤(2.1)所述方法，计算优化变量无制作误差时的磁场均匀度//'。步骤(2.4.2)，按步骤(2.1)所述方法，计算优化变量存在制作误差时的磁场均匀度//"。步骤(2.4.3)，按下式计算制作误差鲁棒性系数，i=|I。步骤(3)，用模拟退火算法对所述目标函数进行优化，其步骤如下步骤(3.1)，向计算机输入模拟退火算法初始参数，包括初始温度，终止温度，温度下降速度，同一温度下最大搜索次数，搜索步长，变量初值。将变量初值置为当前解和最优解，并计算其目标函数值。步骤(3.2)，在当前解的邻域内搜索新解。邻域定义为解空间中，以当前解为中心，其分量中一个或几个在搜索步长设定的幅度(如原变量值的10%)内随机变化的区域。步骤(3.3)，计算机计算新解对应的目标函数值。步骤(3.4)，计算机计算新解对应的约束，如果超过约束则对目标函数采用加罚函数方法惩罚。约束条件为尺寸约束，如果新解的各分量中存在超过设定值的情况，则视为超过约束。在需要对特定性能指标进行限制时，也可将其处理为约束条件，使用与目标函数中对应项相同的方法进行计算，如果超过设定值则进行惩罚。使用加罚函数方法惩罚的方法如下<formula>formulaseeoriginaldocumentpage12</formula>其中，乂为原目标函数值，/2为惩罚后目标函数值，A为惩罚量。步骤(3.5)，按照Metropolis准则判断是否接受新解为当前解，并记录最优解:Metropolis准则如下<formula>formulaseeoriginaldocumentpage12</formula>其中，x,为新解，；c为当前解，r为当前温度，ra^w]为间的随机数。步骤(3.6)，如果同一温度下最大搜索次数达到设定值，则温度下降，否则返回(3.2)。步骤(3.7)，如果温度达到终止温度则算法停止，输出最优解，否则返回(3.2)。本发明中的设计方法，通过加权法设置目标函数，综合考虑上述各项优化指标；通过采用离散变量，避免了通常方法中的取整和离散误差；通过在设计目标中考虑误差鲁棒性指标，实现了强鲁棒性的最优设计。使用本设计方法进行超导磁体设计，可以综合考虑多个设计指标，设计结果具有便于加工制作，误差鲁棒性强的优点。图1是应用本发明设计的一种有源屏蔽超导磁体结构。图2是应用本发明设计的一种有源屏蔽超导磁体的优化变量设置图。图3是应用本发明设计的一种模拟退火启发式算法流程图。图4是应用本发明设计的一种有源屏蔽超导磁体的成像区域磁场均匀度分布。图5是应用本发明设计的一种有源屏蔽超导磁体的0.5mT等值线。表1是应用本发明设计的一种有源屏蔽超导磁体的优化变量。表2是应用本发明设计的一种有源屏蔽超导磁体的设计结果。具体实施例方式下面结合附图来说明一下本发明的原理和具体的实施方式。图1是本发明实施例中所选用磁体结构。线圈关于Z轴对称，并关于Z-O平面对称。磁体内侧为主线圈，采用外凹口形式，外侧为两对屏蔽线圈，屏蔽线圈与主线圈反向串连，磁体内孔半径为0.15m。选择SuperConC54S43型超导线，电流设为/=240A。图2是本发明实施例中优化模型变量设置图。图示为轴截面第一象限。变量X呵xl，x2,x3,x4,x5,x6，x7，x8,x9,x10]，变量各分量含义如下表所示表l<table>tableseeoriginaldocumentpage13</column></row><table>其中xl，x2，x3,x4，x5，x7，x9为离散变量，x6，x8，x10为连续变量。预先对各变量对磁场指标//造成的影响进行人工分析，选择yl=^x8和=^(xl-x2)作为误差鲁棒性系数校验对象，其中^为超导线直径，制作误差取0.5mm。约束条件为磁体半长小于0.5m，屏蔽线圈外径小于0.5m。图3是本发明实施例中用于求解离散优化问题的模拟退火算法流程图，其基本步骤如下(1)设置优化变量，包括各线圈层数、每层匝数、半径和位置。(2)建立前述优化目标函数/。(3)向计算机输入模拟退火算法初始参数，包括初始温度，终止温度，温度下降速度，同一温度下最大搜索次数，搜索步长(原值的10%)，变量初值。将变量初值置为当前解和最优解，并计算其目标函数值。(4)计算机在当前解的邻域内搜索新解。邻域定义为解空间中，以当前解为中心，其分量在搜索步长设定内随机变化的区域。(5)计算机计算新解对应的目标函数值。(6)计算机计算新解对应的约束，如果超过约束则对目标函数采用加罚函数方法惩罚。约束条件为所述尺寸约束，如果新解的各分量中存在超过设定值的情况，则视为超过约束。使用加罚函数方法惩罚的方法如下力为原目标函数值，/2为惩罚后目标函数值，A为惩罚量。(7)按照Metropolis准则判断是否接受新解为当前解，并记录最优解；Metropolis准则如下，<formula>formulaseeoriginaldocumentpage14</formula>其中，、为新解，x为当前解，r为当前温度。ra"《。u为间的随机数(8)如果同一温度下最大搜索次数达到设定值，则温度下降，否则返回(4)。(9)如果温度达到终止温度则算法停止，输出最优解，否则返回(4)继续执行。(10)结束优化算法。优化算法结束后，进行如下后处理，计算最终性能参数使用单积分法计算中心磁场强度&。如果5o-7.00T，则调节电流/,使其相等；使用单积分法计算磁场均匀度使用单积分法计算0.5mT等值线范围；使用单积分法计算线圈内最大磁场强度Smax。将Bmax带入超电线材的Sc—/c特性曲线，求得&=5,时的失超电流化，得到超导线利用率/Zfc;验算所有变量的磁场均匀度//o鲁棒性系数。本发明实施例中的设计结果如卜.表所示表2<table>tableseeoriginaldocumentpage14</column></row><table>权利要求1、一种用于磁共振成像超导磁体设计的离散优化设计方法，其特征在于，依次含有以下步骤步骤(1)，设定如下优化变量，并输入计算机线圈层数和每层匝数，处理为离散变量；线圈半径和位置在不考虑加工精度时，处理为连续变量；在考虑加工精度时，也将线圈半径和位置变量处理为离散变量；所述线圈至少包含主线圈和屏蔽线圈；步骤(2)，按下式计算使用所述计算机采用模拟退火算法进行所述离散优化方法时用的目标函数f，f＝w1·B+w2·H+w3·S+w4·L+w5·R，其中，w1，w2，…w5为预先指定的权系数，是设定值，B为中心磁场强度，H为磁场均匀度，S为杂散场范围，L为超导线用线量，R为制作误差鲁棒性系数；步骤(2.1)，将磁共振成像超导磁体分解为包括电流同向或电流反向的圆柱型线圈在内的线圈组合，使用螺线管内部的Legendre展开法，按下式计算中心磁场强度B和磁场均匀度H首先，对于每个圆柱型线圈，计算前10阶偶次Legendre系数；以第i个线圈为例，其前10阶偶次Legendre系数M0i，M2i，…，M10i，如下式所示<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>β</mi><mi>i</mi></msub><mi>ln</mi><mfrac><mrow><msub><mi>α</mi><mi>i</mi></msub><mo>+</mo><msup><mrow><mo>(</mo><msup><msub><mi>α</mi><mi>i</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow><mrow><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></mfrac></mrow>]]></math></maths><mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>β</mi><mi>i</mi></msub></mrow></mfrac><mrow><mo>(</mo><msup><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>-</mo><msup><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>24</mn><msup><msub><mi>β</mi><mi>i</mi></msub><mn>3</mn></msup></mrow></mfrac><mo>[</mo><msup><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>+</mo><mn>3</mn><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>15</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><msup><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>+</mo><mn>3</mn><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>15</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>)</mo></mrow><mo>]</mo></mrow>]]></math></maths><mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>6</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>240</mn><msup><msub><mi>β</mi><mi>i</mi></msub><mn>5</mn></msup></mrow></mfrac><mo>[</mo><msup><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>8</mn><mo>+</mo><mn>12</mn><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>15</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>-</mo><mn>70</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>4</mn></msup><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0005"num="0005"><math><![CDATA[<mrow><mo>-</mo><msup><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>8</mn><mo>+</mo><mn>12</mn><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>15</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>-</mo><mn>70</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>4</mn></msup><mo>)</mo></mrow><mo>]</mo></mrow>]]></math></maths><mathsid="math0006"num="0006"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>8</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>396</mn><msup><msub><mi>β</mi><mi>i</mi></msub><mn>7</mn></msup></mrow></mfrac><mo>[</mo><msup><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>16</mn><mo>+</mo><mn>24</mn><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>30</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><mn>35</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>4</mn></msup></mrow></mrow>]]></math></maths><mathsid="math0007"num="0007"><math><![CDATA[<mrow><mo>-</mo><mn>2079</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>5</mn></msup><mo>+</mo><mn>3003</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>6</mn></msup><mo>)</mo></mrow>]]></math></maths><mathsid="math0008"num="0008"><math><![CDATA[<mrow><mo>-</mo><msup><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>16</mn><mo>+</mo><mn>24</mn><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>30</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><mn>35</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>4</mn></msup></mrow></mrow>]]></math></maths><mathsid="math0009"num="0009"><math><![CDATA[<mrow><mo>-</mo><mn>2079</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>5</mn></msup><mo>+</mo><mn>3003</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>6</mn></msup><mo>)</mo><mo>]</mo></mrow>]]></math></maths><mathsid="math0010"num="0010"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mn>10</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>11520</mn><msup><msub><mi>β</mi><mi>i</mi></msub><mn>9</mn></msup></mrow></mfrac><mo>[</mo><msup><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>128</mn><mo>+</mo><mn>192</mn><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>240</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><mn>280</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo></mrow></mrow>]]></math></maths><mathsid="math0011"num="0011"><math><![CDATA[<mrow><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>4</mn></msup><mo>-</mo><mn>2772</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>5</mn></msup><mo>+</mo><mn>42042</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>6</mn></msup><mo>-</mo><mn>128700</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>7</mn></msup><mo>+</mo><mn>109395</mn><msup><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mn>8</mn></msup><mo>)</mo></mrow>]]></math></maths><mathsid="math0012"num="0012"><math><![CDATA[<mrow><mo>-</mo><msup><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>128</mn><mo>+</mo><mn>192</mn><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>+</mo><mn>240</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><mn>280</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>3</mn></msup><mo>+</mo></mrow></mrow>]]></math></maths><mathsid="math0013"num="0013"><math><![CDATA[<mrow><mn>315</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>4</mn></msup><mo>-</mo><mn>2772</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>5</mn></msup><mo>+</mo><mn>42042</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>6</mn></msup><mo>-</mo><mn>128700</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>7</mn></msup><mo>+</mo><mn>109395</mn><msup><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mn>8</mn></msup><mo>)</mo><mo>]</mo></mrow>]]></math></maths>其中，<mathsid="math0014"num="0014"><math><![CDATA[<mrow><msub><mi>α</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>A</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msub><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mfrac></mrow>]]></math>id="icf0014"file="A2008101147360003C12.tif"wi="13"he="10"top="148"left="90"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0015"num="0015"><math><![CDATA[<mrow><msub><mi>β</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>L</mi><mi>i</mi></msub><msub><mi>A</mi><mn>1</mn></msub></mfrac></mrow>]]></math>id="icf0015"file="A2008101147360003C13.tif"wi="12"he="10"top="148"left="109"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0016"num="0016"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow>]]></math>id="icf0016"file="A2008101147360003C14.tif"wi="20"he="10"top="165"left="53"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0017"num="0017"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mrow><mn>1</mn><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow>]]></math>id="icf0017"file="A2008101147360003C15.tif"wi="20"he="10"top="164"left="78"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0018"num="0018"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><msup><msub><mi>α</mi><mi>i</mi></msub><mn>2</mn></msup><mrow><msup><msub><mi>α</mi><mi>i</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow>]]></math>id="icf0018"file="A2008101147360003C16.tif"wi="24"he="11"top="164"left="104"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0019"num="0019"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mn>4</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mrow><msup><msub><mi>α</mi><mi>i</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup></mrow></mfrac></mrow>]]></math>id="icf0019"file="A2008101147360003C17.tif"wi="24"he="10"top="164"left="133"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>A1i为第i个线圈的内径，A2i为第i个线圈的外径，Li为第i个线圈的半长；然后，按步骤(2.1)所述方法计算所有圆柱型线圈的所述前10阶中各阶Legendre系数；接着，计算所述磁共振成像超导磁体的各阶Legendre系数M0，M2，…M10；计算时，记各线圈电流密度幅值为|Ji|＝J，记半径最小的线圈为1号线圈，将其他各线圈的各阶Legendre系数进行电流和内径按其折算求和，如下<mathsid="math0020"num="0020"><math><![CDATA[<mrow><msub><mi>M</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>sign</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>M</mi><mi>ji</mi></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>A</mi><mn>11</mn></msub><msub><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mfrac><mo>)</mo></mrow><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow>]]></math></maths>其中，Mj，j＝0，2，…10为磁体Legendre系数，A11为1号线圈内径，A1i为第i个线圈内径，Mji为第i个线圈的j阶Legendre系数，sign(Ji)表示第i个线圈电流方向符号，K为线圈个数；再按下式计算中心磁场强度B，B＝μ0JA1M0；再按下式计算磁场均匀度H，H＝k2M22+k4M42+…+k10M102，其中，μ0为空气磁导率，J为电流密度幅值，A1为最内侧线圈内径，加权系数k2，k4，…，k10为设定值；步骤(2.2)，把所述磁共振成像超导磁体分解成包含电流同向或电流反向的圆柱型线圈在内的线圈组合后，按下述方法使用所述螺线管外部的Legendre展开法计算所述杂散场范围S首先，对每个圆柱型线圈，计算前6阶奇次Legendre展开系数；以第i个线圈为例，其前6阶奇次Legendre系数N1i，N3i，N5i，如下式所示<mathsid="math0021"num="0021"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>β</mi><mi>i</mi></msub><msub><mi>D</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mrow>]]></math></maths><mathsid="math0022"num="0022"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mn>30</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>β</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mrow><mn>9</mn><mi>D</mi></mrow><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>-</mo><mn>20</mn><msub><mi>D</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mo>)</mo></mrow></mrow>]]></math></maths><mathsid="math0023"num="0023"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mn>5</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>56</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>β</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mrow><mn>15</mn><mi>D</mi></mrow><mrow><mn>3</mn><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mn>84</mn><mi>D</mi></mrow><mrow><mn>2</mn><mi>i</mi></mrow></msub><msup><msub><mi>β</mi><mi>i</mi></msub><mn>2</mn></msup><mo>+</mo><mn>56</mn><msub><mi>D</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><msup><msub><mi>β</mi><mi>i</mi></msub><mn>4</mn></msup><mo>)</mo></mrow></mrow>]]></math></maths>其中，<mathsid="math0024"num="0024"><math><![CDATA[<mrow><msub><mi>α</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>A</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msub><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mfrac></mrow>]]></math>id="icf0024"file="A2008101147360004C4.tif"wi="13"he="10"top="195"left="88"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0025"num="0025"><math><![CDATA[<mrow><msub><mi>β</mi><mi>i</mi></msub><mo>=</mo><mfrac><mi>L</mi><msub><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mfrac></mrow>]]></math>id="icf0025"file="A2008101147360004C5.tif"wi="13"he="10"top="195"left="107"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0026"num="0026"><math><![CDATA[<mrow><msub><mi>D</mi><msub><mn>1</mn><mi>i</mi></msub></msub><mo>=</mo><mn>1</mn><mo>+</mo><msub><mi>α</mi><mi>i</mi></msub><mo>+</mo><msup><msub><mi>α</mi><mi>i</mi></msub><mn>2</mn></msup></mrow>]]></math>id="icf0026"file="A2008101147360004C6.tif"wi="26"he="5"top="212"left="66"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>D2i＝1+αi+αi2+αi3+αi4D3i＝1+αi+αi2+αi3+αi4+αi5+αi6A1i，A2i和Li含义同前；然后，按步骤(2.2)所述方法计算所有圆柱型线圈的各阶Legendre系数；接着，计算所有圆柱型线圈的各阶Legendre系数的代数和，用N1，N3，N5表示；计算时，记各线圈电流密度幅值为|Ji|＝J，记半径最小的线圈为1号线圈，将其他各线圈的各阶Legendre系数进行电流和内径按其折算求和，如下<mathsid="math0027"num="0027"><math><![CDATA[<mrow><msub><mi>N</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>sign</mi><mrow><mo>(</mo><msub><mi>J</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>N</mi><mi>ji</mi></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>A</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><msub><mi>A</mi><mn>11</mn></msub></mfrac><mo>)</mo></mrow><mrow><mi>j</mi><mo>+</mo><mn>3</mn></mrow></msup><mo>,</mo></mrow>]]></math></maths>其中，Nj，j＝1，3，5为磁体Legendre系数，Nji为第i个线圈的j阶Legendre系数，其他含义同前；再计算杂散场范围SS＝l1N12+l3N32+l5N52，其中，l1，l3，l5为加权系数，为设定值；步骤(2.3)，按下式计算超导线圈用线量LL＝V/Scr，其中，V为线圈总体积，Scr为导线截面积；步骤(2.4)，按下式计算磁体鲁棒性系数R；步骤(2.4.1)，按步骤(2.1)所述方法，计算优化变量无制作误差时的磁场均匀度H′；步骤(2.4.2)，按步骤(2.1)所述方法，计算优化变量存在制作误差时的磁场均匀度H″；步骤(2.4.3)，按下式计算制作误差鲁棒性系数R，R＝|H′-H″|；步骤(3)，用模拟退火算法对所述目标函数进行优化，其步骤如下步骤(3.1)，向计算机输入模拟退火算法初始参数，包括初始温度，终止温度，温度下降速度，同一温度下最大搜索次数，搜索步长，变量初值；将变量初值置为当前解和最优解，并计算其目标函数值；步骤(3.2)，在当前解的邻域内搜索新解；邻域定义为解空间中，以当前解为中心，其分量中一个或几个在搜索步长设定的幅度(如原变量值的10％)内随机变化的区域；步骤(3.3)，计算机计算新解对应的目标函数值；步骤(3.4)，计算机计算新解对应的约束，如果超过约束则对目标函数采用加罚函数方法惩罚；约束条件为尺寸约束，如果新解的各分量中存在超过设定值的情况，则视为超过约束；在需要对特定性能指标进行限制时，也可将其处理为约束条件，使用与目标函数中对应项相同的方法进行计算，如果超过设定值则进行惩罚；使用加罚函数方法惩罚的方法如下f2＝f1+Δ，其中，f1为原目标函数值，f2为惩罚后目标函数值，Δ为惩罚量；步骤(3.5)，按照Metropolis准则判断是否接受新解为当前解，并记录最优解；Metropolis准则如下其中，x1为新解，x2为当前解，T为当前温度，rand为间的随机数；步骤(3.6)，如果同一温度下最大搜索次数达到设定值，则温度下降，否则返回(3.2)；步骤(3.7)，如果温度达到终止温度则算法停止，输出最优解，否则返回(3.2)。2、根据权利要求1所述的一种磁共振成像超导磁体设计的离散化优化方法，其特征在于，在计算磁场均匀度的鲁棒性系数时，优化变量的变化量为0.5mm或lmm。全文摘要一种用于磁共振成像(MRI)超导磁体的离散优化设计方法，属于超导磁体设计技术，其特征在于，通过加权法设置目标函数，综合考虑磁场强度，磁场均匀度，杂散场范围，用线量，制作误差等指标。采用离散变量作为优化变量，避免了通常方法中的取整和离散误差。通过在设计目标中考虑误差鲁棒性指标，实现了强鲁棒性的最优设计。使用本设计方法进行超导磁体设计，可综合考虑多个设计指标，设计结果具有便于加工制作，误差鲁棒性强的优点。文档编号G06F17/50GK101339580SQ20081011473公开日2009年1月7日申请日期2008年6月11日优先权日2008年6月11日发明者梁嘉祺,蒋晓华,伟郭申请人:清华大学