用于超声装置中的EMI减小的多信道实时相位调制

摘要

用于减小由于使用阵列中的多个超声换能器而引起的电磁干扰(特别是在在磁共振成像装置内部的人体的内部)的方法。驱动换能器的电信号相对于彼此在相位上偏移以便实现由这样的信号和换能器引起的电场和磁场的最大偏移。相位偏移被动态地调整以对驱动振幅和频率的变化作出响应以便保持电磁干扰的最佳减小。

权利要求

1.一种减小由超声换能器集合引起的电磁干扰的方法，所述超声换能器集合是热治疗设备中的超声换能器阵列的一部分，所述集合包括N个换能器，所述集合中的每个换能器对应于作用的信道，并且被用成一振幅、频率和相位角的驱动信号电驱动，所述驱动信号的所述频率对于所述集合中的所有换能器都是相同的，所述方法包括确定和设置每个驱动信号的相位角Θ1，Θ2，...，ΘN，这样的相位角的确定和设置包括：确定每个换能器的相应的驱动信号的振幅A1，A2，...，AN，每个振幅是非负实数；

确定所述振幅Am中的一个是否大于除了Am之外的所有振幅的总和A1+A2+…+Am‑1+Am+1+…+AN；

如果Am大于或等于这样的其他的振幅的总和，则对于所有的不等于m的i，设置Θm＝

180°，并且设置Θi＝0°，其中，i＝1,2,…,N，Θi是第i个驱动信号的相位角；

如果Am小于这样的其他的振幅的总和：

确定矢量 这样的矢量包括N个元素，每个这样的元素的值要么是+1，要么是‑1，以使得 和包括所有的振幅的矢量[A1，A2，...，AN]的标量积是非负的，并且在幅值上不大于包括N个元素的任何其他可能的矢量的标量积的幅值，并且矢量[A1，A2，...，AN]包括所有的振幅；

对于i＝1,2,…,N，定义包括N个元素的矢量 以使得每个元素Oi等于Pi和Ai的乘积；

确定 的第一个正元素Oa，Oa不小于 的任何其他元素，其中a是 的第一个正元素的位置的索引；

确定 的第二个正元素Ob，Ob不小于 的除Oa之外的任何其他元素，其中b是 的第二个正元素的位置的索引；

将量γ定义为 的除Oa和Θb之外的任何元素的总和的绝对值；以及将Θa和Θb设置为：

‑1

并且对于所有的Θi，除了Θa和Θb之外，设置Θi＝cos Oi；

其中，

Θa是作用的信道a中驱动信号的相位角，Θb是作用的信道b中驱动信号的相位角。

2.如权利要求1所述的方法，其中所述超声换能器阵列被部署在磁共振成像设备的内部。

3.如权利要求2所述的方法，其中所述超声换能器阵列是出于向人类患者施加适形热治疗的目的而被部署的。

4.如权利要求1所述的方法，其中所述振幅和频率每隔一段时间被更新，并且这样的相位角的确定和设置按每个这样的间隔、基于由这样的更新引起的新的振幅和频率被重复。

5.如权利要求4所述的方法，其中更新所述振幅和频率包括将所述阵列中的一个或更多个超声换能器添加到所述集合和/或从所述集合移除一个或更多个超声换能器。

6.如权利要求4所述的方法，其中与所述振幅和频率的更新相结合的、所述驱动信号的相位角的设置与所述驱动信号中的所述新的振幅和频率的实现同步地在所述驱动信号中被实现。

7.如权利要求1所述的方法，其中所述驱动信号是正弦信号。

8.如权利要求1所述的方法，其中所述矢量 的确定通过计算所述矢量[A1，A2，...，AN]的相应的标量积并且比较所得结果来实现，其中每一个可能的矢量的集合中的每个包括N个元素，所述N个元素中的每个要么是1，要么是‑1，或者具有这样的集合的子集。

9.如权利要求1所述的方法，其中所述超声换能器阵列是线性阵列。

10.如权利要求1所述的方法，其中所述超声换能器阵列是聚焦阵列。

说明书

用于超声装置中的EMI减小的多信道实时相位调制

技术领域

[0001] 本申请一般地涉及超声治疗系统，并且更具体地涉及用于通过以下方式改进超声治疗的有效性的方法，即，通过对发送到超声换能器的电信号使用相位调制来减小由这样的换能器引起的电磁干扰(EMI)。

背景技术

[0002] 超声换能器已经被用于超声治疗系统中以实现病变组织和其他组织的治疗加热。操作来形成超声能量的射束的超声换能器阵列在受影响的组织区域或处置体积
(treatment volume)中使声音转换为热能、随后在处置体积中引起有益的温度升高。

[0003] 在图像引导的超声治疗系统中，患者和超声治疗设备一般被设置在成像体积(诸如磁共振成像(MRI)设备)中，这使得可以引导敷抹器放置，并且另外使得可以通过提供可从其计算温度图的实时数据来监视组织的处置效果。临床操作者然后可以监视处置体积或病变组织内的治疗的进展，并且可以基于来自结果的输入和处置的进展来对超声功率信号做出手动或自动改变。通过适当地监视加热效果，超声治疗系统可以被用来处置有害细胞并且在使对于健康组织的损害最小化的同时可控地摧毁肿瘤。

[0004] 在男性中证实使用MRI引导的经尿道超声治疗系统处置诸如前列腺癌的疾病的工作已经被进行。参见例如Chopra等人，“MRI‑compatible transurethral ultrasound system for the treatment of localized prostate cancer using rotational control”,Med Phys 35(4):1346‑1357,2008年。还参见美国公开2007/0239062；美国专利6,589,174，“Technique and apparatus for ultrasound therapy”，2003年；美国专利7,
771,418，“Treatment ofdiseased tissue using controlled ultrasonic heating”，
2010年；美国专利8,998,889，“System and method for control andmonitoring ofconformal thermal therapy”，2015年；美国专利9,707,413，“Ultrasonic therapy applicator”，2017年。这样的系统(包括本申请人所进行的或者为本申请人进行的累积的公布的且获得专利权的工作，所有这些工作特此通过引用并入)教导了对病变的前列腺使用经尿道超声能量以在病变的组织中达到期望的目标温度以实现临床结果，所述临床结果通常是前列腺中的病变组织细胞的坏死。实时的处置的MRI引导和温度监视使得能够控制供给超声治疗换能器的功率，而且还控制沿着在病变的前列腺附近插入到患者的尿道中的细长敷抹器轴向设置的这样的换能器阵列的旋转。

[0005] 如本领域中技术人员已知的，超声换能器被构造和操作为获取电功率并且在一般被称为转导的过程中从换能器元件的表面生成超声能量波。转导的性质和程度取决于用于构造换能器的材料、换能器几何形状以及输入到换能器的电输入。构造超声换能器所用的常见的材料是以几种形式出现的压电换能器晶体材料(锆钛酸铅，PZT)。

[0006] 在以上引述的来源中所公开的用于超声热治疗的系统中，通常在频带4‑4.5MHz和13‑14.4MHz中产生射频(RF)电输出，并且使用该输出来驱动将电能量转换为声学压力(例如，超声)的多达10个压电元件。所述元件可以在相对较高的功率下被驱动，分别地，在低频带和高频带处在大约4W和2W的功率下被驱动，并且已知的是，PZT材料可以产生非线性响应。这种情况的影响是产生以谐波形式的电磁干扰(EMI)，所述谐波可能污染MRI图像，并且当换能器被用在MRI设备内部时，可能干扰MRI测温，这使其操作依赖于RF电磁信号。这样的EMI是由在换能器元件中产生的电场、以及由于流入和流出换能器元件的电流产生的磁场引起的，因为如本领域中众所周知的，时变的电场和磁场可以产生电磁辐射。

[0007] 因此，需要通过减小EMI来改进超声热治疗的准确性和有效性。当EMI的来源被部署在人体内部的狭窄的空间中(例如，男性尿道和前列腺中)时，传统的屏蔽方法不总是可行的，或者具有有限的功效。因此，需要不增大部署的设备的主体、因此没有伤害患者的风险、减小EMI的其他方法。本公开是针对用于通过调整发送到超声换能器阵列中的不同换能器的电信号的相位角来减小由于该阵列引起的EMI的方法。

发明内容

[0008] 本文中的公开是针对用于通过使用相位调制技术减小由于超声换能器阵列和连接到这样的换能器的电线引起的电磁干扰来改进治疗过程和其他过程中的超声的使用的有效性的方法。所提出的方法为发送到相应的换能器的信号确定并且实现相位偏移角，以便通过阵列中的相应的换能器产生的电场和电流的相互偏移来实现EMI的最佳减小程度。

[0009] 实施方案是针对一种减小由于作为热治疗设备中的超声换能器阵列的一部分的超声换能器集合引起的电磁干扰的方法，所述集合包括N个换能器，所述集合中的每个换能器对应于作用的信道，并且被用成一振幅、频率和相位角的驱动信号电驱动，所述驱动信号的所述频率对于所述集合中的所有换能器都是相同的，所述方法包括确定和设置每个驱动信号的相位角Θ1，Θ2，....，ΘN，这样的相位角的确定和设置包括：确定每个换能器的相应的驱动信号的振幅A1，A2，...，AN，每个振幅是非负实数；确定所述振幅Am中的一个是否大于除了Am之外的所有振幅的总和A1+A2+...Am‑1+Am‑1+Am+1+...+AN；如果Am大于或等于这样的其他的振幅的总和，则对于所有的不等于m的i，设置Θm＝180°，并且设置Θi＝0°；如果Am小于这样的其他的振幅的总和：确定矢量 这样的矢量包括N个元素，每个这样的元素要么是1，要么是‑1，以使得 和包括所有的振幅的矢量[A1，A2，...，AN]的标量积是非负的，并且在幅值上不大于包括N个元素的任何其他的可能的矢量的标量积的幅值，每个这样的元素要么是1，要么是‑1，并且矢量[A1，A2，...，AN]包括所有的振幅；对于i＝1,2,…,N，定义包括N个元素的矢量 以使得每个元素Oi等于Pi和Ai的乘积；确定 Oa的不小于 的任何其他元素的第一个正元素；确定 Ob的不小于 的除Oa之外的任何其他元素的第二个正元素；将量γ定义为 的除Oa和Ob之外的任何元素的总和的绝对值；将Θa和Θb设置为：

[0010]

[0011]

[0012] 并且对于所有的Θi，除了Θa和Θb之外，设置Θi＝cos‑1Oi。

[0013] 另一个实施方案是针对一种用于减小操作电驱动的超声热治疗设备时的电磁干扰的方法，所述方法包括：相对于指定的处置区域定位包括所述设备的超声阵列的所述设备；在耦合到所述治疗设备的基于计算机的主机单元中，确定共同的驱动频率，并且确定用其驱动所述阵列的多个换能器元件中的每个的振幅；用相应的电压源产生的相应的驱动信号驱动所述多个换能器元件中的每个，用于每个元件的所述驱动信号包括所述共同的驱动频率、所述振幅和相应的相位角；以及修改至少一个驱动信号的相位角以减小所述热治疗设备的净电磁输出。

[0014] 处置区域可以是患者的身体内的内腔、孔口或其他自然体积或人造体积，包括尿道、直肠或其他器官或腔体。在这样的情况下，处置在内部(例如，经尿道)被递送。在其他情况下，处置从外部被递送给身体，并且超声能量被引导到身体中，例如，通过患者的皮肤、外部的器官和组织层。处置设备因此可以包括成直线(沿着直线装置的行或轴线)的阵列，或者可以是具有弯曲的、有轮廓的或其他的几何布置的几何聚焦的阵列。

附图说明

[0015] 为了更充分地理解本发明的性质和优点，结合附图参照以下对优选实施方案的详细描述，在附图中：

[0016] 图1图示说明用于向患者提供图像引导的超声治疗的示例性系统；

[0017] 图2图示说明示例性超声阵列；

[0018] 图3是图示说明两换能器阵列的简化模型的示意性电路图，通过该电路图，展示了构成本文中所公开的方法的基础的某些原理；

[0019] 图4是两音调混音模型的示意图，该图进一步展示了构成本文中所描述的方法的基础的原理中的一些；

[0020] 图5是示出在使用本文中所公开的方法的情况和不使用本文中所公开的方法的情况这两种情况下被驱动的两个换能器的结果的说明性绘图，该图展示了通过使用这样的方法获得的结果；

[0021] 图6A、6B、6C和6D是用几何展示本文中所公开的方法的某些方面的矢量图；以及[0022] 图7A、7B、7C和7D图示说明示出根据实施方案的、本文中所公开的方法的操作的流程图。

具体实施方式

[0023] 本文中的公开是针对用于使用相位调制技术来减小由于用于超声热治疗的设备中的超声换能器的多元素阵列引起的EMI的方法。

[0024] 图1图示说明用于向患者提供图像引导的超声治疗的示例性系统10。简化的图示示出了主计算机100，诸如便携式PC、工作站、或具有处理器、存储器并且耦合到某个输入/输出设备的其他处理装置。主计算机100可以包括显示器，并且可以支持用户界面110以便利热治疗处置过程的控制和观察。

[0025] 主计算机100适于通过计算机接口连接器120耦合到其他系统和组件。连接装置120将数据和信息传载到主计算机100并且从主计算机100传载数据和信息，并且可以包括标准的或专用的电布线连接电缆，诸如串联连接电缆等。此外，连接装置120可以如无线通信领域中的技术人员已知的那样无线地实现，并且可以进一步通过多次连接、通过网络、或者通过另一合适的方法来实现。

[0026] 在一些实施方案中，主计算机100通过连接装置120耦合到功率控制单元130。功率控制单元130可以被实现为独立的硬件设备，但是可以被实现为主计算机100的一部分，例如，通过被构建到容纳这样的硬件组件的计算机或服务器系统中的特殊卡中。

[0027] 功率控制单元130可以专门至少包括适于处理机器或程序指令的处理器，这些指令可以从系统10的另一组件提供给处理器，并且可以被存储在功率控制单元130中的存储器装置上。包括模拟和/或数字电路系统的电路系统可以在功率控制单元130内被操作以便确定供给超声治疗设备150中的一个或更多个超声治疗换能器元件的输出功率。

[0028] 在一些实施方案中，功率控制单元130可以将受控的电驱动信号递送给超声治疗设备150中的多个超声换能器元件(例如，PZT阵列元件)。驱动信号可以被控制为将程控的功率量递送给治疗设备150的每个元件或一组元件。驱动信号也可以被控制为向治疗设备150的所述超声换能器提供确定的驱动电压、电流、振幅、波形或频率。如本文中所讨论的，驱动信号的相对相位也可以被控制为例如减小EMI。这样的电驱动信号通过合适的导线、电缆或总线140从功率控制单元130被传载到超声治疗设备150。适当的插件接口或连接器可以被包括以便配合与它们的相关联的组件往来的连接器或总线的各种端部。

[0029] 在操作中，超声治疗设备150包括部分155，该部分155被插入到患者的身体的一部分中，以将合适剂量的超声能量递送给患者的身体的病变的区域中的组织。

[0030] 患者和超声治疗设备150一般被设置在成像体积160(诸如磁共振成像(MRI)设备)中，成像体积160可以将患者的相关部分(例如，处置体积)的实时图像提供给主计算机100或显示器和用户界面110。在一些实施方案中，热治疗的实时监视被执行以使得临床操作者可以监视处置体积或病变的组织内的治疗的进展。可以基于来自结果的输入和处置的进展来对来自功率控制单元130的功率信号做出手动或自动改变。

[0031] 处置系统组件对于系统10中的控制组件的反馈和耦合可以被用来确保最佳的射频(RF)功率信号被提供给病变的组织的处置中所用的超声阵列155的每个元件。一些示例包括在男性患者中使用MRI引导的超声治疗应用的前列腺癌肿瘤的处置。

[0032] RF功率控制单元130可以包括单独的具有单个的处理器、放大器、滤波器和其他组件的电路卡以实现输出到超声处置设备150的超声阵列155的元件的期望的驱动功率。可替换地，单个处理器可以被用来控制对于每个阵列元件的各种功率信道的行为。

[0033] 图2图示说明超声阵列200，超声阵列200与超声阵列155可以是相同的。超声阵列200中的每个超声换能器205经由至少一个电线220被单独的电压源210驱动。因此，用于每个换能器205的每个驱动信号的电压和相对相位可以被控制。电压源210被包括在RF功率控制单元230中，RF功率控制单元230与RF功率控制单元130可以是相同的。注意，RF功率控制单元230可以包括附加的组件，诸如如以上关于RF功率控制单元130讨论的电路、处理器、放大器、滤波器和其他组件。

[0034] 每个换能器205具有2个电端子，其中第一个经由被称为驱动线的电线被电连接到其相应的电压源210的一个端子，并且其中第二个与其他所有的这样的第二换能器端子共同地，经由被称为返回线或地线的电线被电连接到电压源210的其他端子。驱动线(阵列200中的每个超声换能器一个)和公共返回线连接到敷抹器的近端，并且延伸到电压源210。当所述设备被用于热治疗时，在典型的布置中，敷抹器将被插入到患者的身体的内部以在身体的接收治疗的一部分的附近，例如，被经尿道插入到病变的男性前列腺的附近，此时患者在MRI机器的内部。通常被一起绑在单个护套的内部的驱动线和返回线从敷抹器延伸到MRI机器的外部，在MRI机器中，安置了电压源，连同用于控制所述设备的各种装置，诸如RF功率控制单元230。

[0035] 换能器(每个由单独的电压源驱动)的线性阵列的使用使得可以更精确地控制对于正被处置的组织的超声能量的应用，从而与换能器全都在相同的电压下被驱动时的情况相比，提高热治疗的有效性。在对于每个换能器使用单独的电压源的情况下，每个换能器的电压振幅可以根据处置计划，响应于反馈被动态地调整，以在任何给定的时间，保持针对每个换能器的位置和方向输出的超声的最佳水平。用于使用超声控制和监视热治疗的系统和方法在例如标题为“用于超声治疗系统的RF功率控制器(RF Power Controller for Ultrasound Therapy System)”的美国专利申请公开No.2011/0270366和标题为“用于适形热治疗的控制和监视的系统和方法(System and Method for Control and Monitoring of Conformal Thermal Therapy)”的美国专利No.8,998,889中被描述，这些专利和申请特此通过引用并入。

[0036] 每个电压源以给定的电压振幅(任何给定时间时的该振幅由处置计划以及系统的控制算法和机制确定)将给定频率(被称为驱动频率)的交流(AC)信号(通常是正弦信号)发送给换能器中的一个。AC电压信号在换能器内引起振荡电场，这继而借助于压电感应在换能器中引起机械振荡；这些机械振荡将超声声学能量传送到正被处置的组织中，其中这样的能量最终被转换为热能以实现预期的治疗效果。

[0037] 如以上所指出的，MRI设备内部的EMI可以产生于由流过驱动线和返回线的电流产生的磁场、以及在换能器中产生的电场。当多个换能器被单独的电压源以共同的驱动频率驱动时，可以通过使电压信号彼此异相来减小来自这些源的EMI。这样，通过叠加原理，方向彼此相反的电场分量彼此将会叠加和偏移，如方向彼此相反的电流将会叠加和偏移的那样。当多个超声换能器被部署在线性阵列上彼此邻近时，换能器产生的电场将在空间上部分重叠。类似地，彼此紧邻的导线(诸如驱动线)产生的磁场也将在空间上重叠。因此在两种情况下，方向相反的场将会偏移，从而减小EMI。并且返回线中的电流是来自所有的换能器的电流的叠加，所以这样的电流将直接偏移到它们在相反方向上的程度，从而减小由这样的电流引起的任何的因此而产生的磁场。

[0038] 本文中所公开的相位调制方法背后的原理可以通过考虑包括用共同的频率驱动的2个超声换能器的阵列的简单情况来证实。图3是示出2个换能器301和302的简化的示意图，每个换能器被表示为电容器，每个换能器被连接到正弦电压源311和312中的一个，并且由该电压源供电，被驱动线321和322连接，并且被共享的返回线324连接。换能器301和302内的和周围的电场分别用场线331和332示出。

[0039] 来源于源311和312的正弦电压信号可以被描述为时间的函数如下：

[0040] V1(t)＝A cos(ω0t+φ1)

[0041] V2(t)＝B cos(ω0t+φ2)

[0042] 这里，A和B是相应信号的电压振幅，ω0是基于2π弧度或360度的整个周期，以每秒弧度为单位或者可替换地以每秒度数为单位测量的共同的驱动频率，φ1和φ2是以弧度为单位或者可替换地以度数为单位测量的相位角或相位偏移。

[0043] 注意，给定频率ω的正弦信号可以被认为是在复平面中、以每秒ω弧度(或度数)或者每一个2π弧度(或360度)一个整周期的速率逆时针旋转的相位矢量或相量；矢量的长度或幅值是信号的振幅，并且任何给定时间的信号等于相量到水平(实)轴(即，相量的实数部分)上的投影。信号的相位偏移或相位角是在时间t＝0、被指定为周期的起点的任意时间(该时间对于正被考虑的所有相量是相同的)、从轴到相量逆时针测量的、相量的相对于正实轴的角度位置。当共同的频率的多个AC信号被表达为相量时，矢量算术可以被用来计算这样的信号的组合效果。本文中所公开的方法使用这样的计算技术。

[0044] 返回到图3，流过换能器301和302的电流分别指定为I1和I2，其中返回线324中的电流被指定为IT，其中IT＝I1+I2。如果相位角φ1和φ2被设置为使得V1和V2为相反的相位，诸如通过设置φ1＝0°和φ2＝180°，并且A＝B，即，信号的幅值相等，则V1＝‑V2并且I1＝‑I2，并且和IT＝I1+I2将被相加为零，即，流过换能器的净电流为零，并且返回线中的电流为零。

[0045] 该情形在图3中被描绘，在图3中，电场线331和332被示为指向相反的方向；换能器中的电场的方向将取决于施加的电压的意义。假定，换能器301和302在物理上彼此并联，因为它们将在线性换能器阵列上。可以看出，来自相应的换能器301和302的电场在换能器之间的空间333中重叠。按照叠加原理，同一空间中的这样的重叠的场将会彼此偏移到它们指向相反的方向的程度。这种场抵消使由这样的场引起的RF发射减小。另外，来自返回线的辐射减小，因为返回线传载更小的电流，因此产生更小的磁场。此外，由于这样的电流的方向相反，驱动线中的电流产生的磁场将是彼此反向的，因此这样的磁场也将彼此叠加和偏移，进一步减小RF发射。

[0046] 超声换能器中所用的压电材料(诸如PZT)已经被观察到响应于施加的信号，表现出非线性行为。在两个作用的元件的简单的情况下，两音调非线性混合模型可以被用来理解因此而产生的谐波和互调内容：

[0047]

[0048] 在以上方程中，V1和V2表示施加于两个被驱动(在这种情况下，以ω0的共同的频率操作)的元件的信号。功率电平自由变化。V0是包含新的谐波和互调内容的“混合器”的输出，这里用PZT传递函数的幂级数逼近表示。图4以示意性的形式示出混合模型。输入信号410和420被馈送到混合器400中，混合器400输出信号430。为了简化，特定于压电的模型参数已经被忽视。在A＝B并且信号为相反相位(例如，φ1＝0°并且φ2＝180°)的情况下，那么V1＝‑V2并且求和项取消，从而消除了谐波内容。

[0049] 在两个元件的情况下对信号进行相位调制的效果可以在图5中被看见，图5描绘了基于当以相同的振幅和频率(可替代地同相的和异相180度)驱动2个元件时获得的实验结果的说明性绘图500。水平轴510表示频率，垂直轴520表示组合信号的功率。迹线530示出同相(φ1＝φ2)驱动信号的组合，而迹线540示出异相(φ1+180°＝φ2)驱动信号的组合。峰550表示以被驱动的压电元件的自然频率的整数倍出现的谐波。可以容易地看出，当驱动信号异相时，峰低得多，从而证实了由彼此异相的驱动信号造成的谐波内容大幅减少。

[0050] RF电压信号可以是提前的或推迟的，也就是说，它们可以使它们的相位角被调整，而不极大地影响对于目标的超声声学能量的递送。换能器产生的超声能量取决于驱动该换能器的电压信号的振幅和频率，而它不取决于该信号的相位；并且这样的能量在取决于换能器的位置和几何形状以及周围组织、而不取决于该信号的相位的方向上被传送。因此，本文中所公开的方法(其涉及选择并且实现用于驱动超声换能器阵列的各种信号的不同的相位角)对于减小来自这样的信号的EMI是有效的，同时保持这样的信号对于驱动换能器以实现期望的治疗目标的有效性。

[0051] 以上所讨论的两元件模型可以容易地被扩展到多个元件，其中这样的元件的一个或更多个子集被以共同的频率驱动。分配相位角的简单的方法将是将共同的频率的交替的元件设置为0/180/0/180…度。然而，该方法只有在所有元件上输出的功率相同时才是最佳的。在许多超声应用中，可取的是，以相同的频率、但是不同的功率电平驱动阵列上的不同的换能器元件，并且在过程期间，响应于反馈动态地调整这样的功率电平以便实现该过程的治疗目标或其他目标。就在换能器元件之间并且随着时间变化的功率电平而言，需要更好的算法。

[0052] 本文中所公开的方法是针对确定并且设置驱动超声换能器阵列的正弦电压信号集合的相位角，以使得对于给定驱动频率的每个信号子集，这样的子集中的所有信号的矢量和最小化。在计算矢量和时，给定频率的每个电压信号被表达为相量，也就是说，被表达为复平面中的具有一长度、或等于该信号的振幅的幅值、并且成从正实轴逆时针测量的、等于这样的信号的相位角或相位偏移的角度的矢量。也就是说，诸如V1(t)＝A1cos(ω0t+φ1)的信号在相量表示法中将被表示为长度为A1的、从正实(即，指向右的)轴逆时针指向角度φ1的矢量。可替换地，作为矢量，相量可以用其分量(即，其实数部分和虚数部分)被表示为[A1 cosφ1，A1 sinφ1]、或者被表示为单个复数 其中2
j ＝‑1。当被这样表示时，相量可以被乘以 以获得作为时间的函数的电压信号：
。因为给
定频率的所有的相量都被乘以相同的因子以获得因此而产生的信号，所以给定频率的两个或更多个信号的组合效果可以通过简单地将这样的信号的相量相加来计算。AC信号的这样的表示在电领域中是众所周知的。

[0053] 给定频率的所有信号的矢量和将是表示正被发送到阵列中的正被以该频率驱动的所有换能器的净信号。通过最小化该净信号，驱动线和返回线中的净电流(其电流与这样的线传载的电压信号成比例)也将最小化。所述线中的净电流的减小将使由这样的电流引起的磁场减小，从而将使RF发射减小。同样地，施加于设置在线性阵列上、因此物理上彼此并联的超声换能器元件的净电压的减小将导致这样的换能器产生的电场的偏移更大，从而来自这样的场的RF发射减小。

[0054] 对于正被以给定的共同的频率驱动的给定的换能器子集，本文中所公开的算法可以被用来优化发送到集合中的每个换能器的信号的相应的相位角。设N为子集中的换能器的数量，并且设A1、A2等一直到AN为将被发送到集合中的换能器的信号的相应的振幅(在电压上，或者可替换地在功率上)；这样的振幅将通过用于这样的目的的任何手段来确定和选择，以便实现超声过程的治疗目标和/或其他目标。注意，所有的Ai都是非负数，并且对于所有的作用的信道，一般都将是正的。(不作用的信道将用零振幅表示)。目标是找到相位角θ1,θ2,…,θN的集合以使得所有的相量的矢量和最小化。

[0055] 所述算法继续进行如下：定义由振幅组成的N矢量(也就是说，N个标量量的有序集合)

[0056]

[0057] 作为第一遍或粗略的逼近，我们仅考虑0度和180度的相位角、以及这样的角度的什么组合将使总信号最小化。180度的相位角等同于将信号乘以‑1，从而使该级的计算简化。计算振幅的和并且确定最大振幅：

[0058]

[0059] Amax＝max[A1，A2，…，AN]

[0060] 确定 是否成立，即，中的最大的元素的幅值是否大于或等于 中的所有元素的和的一半，或者换句话说，最大的振幅是否大于或等于所有其他的振幅的和。注意，当N＝2时，情况将总是如此。可以通过相位调制实现的信号的最大偏移于是通过将具有最大的振幅的信号的相位设置为180度、并且将所有其他的信号的相位设置为0度来获得。
如果 则信号将完全偏移，并且如果 则信号将不会完全偏移，而
是将偏移到对于这样的振幅可能的最大程度。

[0061] 尤其是当数量较多的(例如，10个)换能器被以给定频率驱动时，情况通常可以是，最大的振幅将不大于所有其他的振幅的和，即， 在这种情况下，下一个步骤是确定将得到信号的最大偏移的、0度和180度相位角(或者，分别地，等于1和‑1)的组合。为了该目的，将“相位偏移矢量”定义为由1和‑1组成的N矢量。(注意，该矢量不同于以上所讨论N的相位矢量或“相量”。)存在2 个可能的这样的相位偏移矢量，这些矢量表示1和‑1的所有的可能的置换。对于给定的相位偏移矢量 由施加用 表示的相位角(0度和180度)而产生的净振幅或残留振幅可以被计算为 和 的标量积，也被称为内积或点积。该乘积是通过将
2个矢量中的N对对应元素相乘、然后求取所得的N个乘积的和、得到标量结果而计算得到的：

[0062]

[0063] 所述算法中的下一个步骤是确定将得到具有最小幅值或绝对值的残留的R的相位偏移矢量

[0064]

[0065] 这里，pk是所有的可能的2N个相位偏移矢量。例如，如果N＝3，则 至 将如下(注意，这里的次序没有关系)：

[0066]

[0067]

[0068]

[0069]

[0070]

[0071]

[0072]

[0073]

[0074] 可以通过获取 与所有的可能的 的标量积并且比较结果、或者通过本领域中已知的其他方法(诸如通过使用查找表，或者使用其成本函数根据某个惩罚度量计算残差的迭代最小化算法)来确定。注意，按照对称性，对于每个 将存在其中所有元素都被保留的另一个 例如，在以上列表中 所以2个 将得到R个相等幅值、但是相N‑1
反符号的值。因此，因为我们仅对残差的幅值感兴趣，所以最多2 个可能性需要被检查以确定Rmin。

[0075] 设 ＝得到Rmin的 中的一个：

[0076]

[0077] 可能存在多于一个的这样的 但是任何一个都可以被选择。按照对称性，对于其与 的内积为负的每个可能的 将还存在其与 的内积为正并且幅值相同的另一个 即，符号反过来的第一个 不失一般性，如果Rmin>0，则我们可以任意地选择与 得到正的内积的 在一些实施方案中，如果存在得到相同的Rmin的多个可能的相位偏移矢量，则基于其他准则(诸如设法使在彼此相反相位上被驱动的多对物理上相邻的换能器的数量最大化)，从这些可能的候选中选择相位偏移矢量

[0078] 如果Rmin＝0，则电压信号的整个偏移可以通过根据 设置信号的相位角来实现；也就是说，如果 那么如果Pi＝1，则来自信道i的每个信号(具有振幅
Ai)被设置为具有相位0度，或者如果Pi＝‑1，则被设置为具有相位180度。

[0079] 如果Rmin>0，则这样(即，根据粗略逼近步骤)设置相位角将不会得到电压信号的整个偏移，而是将使净信号减小到振幅Rmin中的一个。在这种情况下，所述算法继续通过获取电压信号中的两个并且进一步调整它们的相位角以便实现自所有信号的整个偏移。

[0080] 我们将N矢量 定义为 与 的逐元素相乘，其中 是如上所述那样选择的：

[0081]

[0082] 假定 已经如上所述被选择为使得 即， 的元素的和是正的，则设α是的最大的正元素，其中a是其位置的索引，并且设β是 的第二大的正元素，其中b是其位置的索引：

[0083]

[0084]

[0085] 注意，对于Rmin>0，在 中将总是存在至少2个正的元素。(如果情况并非如此，则所述一个正的元素将大于所有的其他的元素的幅值的和，这与如以上所讨论的的情况是不一致的。)设γ为 中的其余的元素的和：

[0086] 注意，Rmin＝α+β+γ。

[0087] 注意，γ不总是负的。(如果γ是正的，则Rmin将不会是最小的可能的残差，因为它可以通过将Pa和Pb这二者的符号反过来而变得更小。)还注意，γ的幅值将小于α和β的和(如果γ不是，则我们将具有Rmin<0)。γ的幅值将大于α和β的差：

[0088] γ＜0

[0089] α‑β＜|γ|＜|α+β|

[0090] 量γ表示来自除了α和β所表示的那些信道之外的所有信道的电压信号的因此而产生的信号的相量。下一个步骤是对这两个信道设置相位角以使得它们偏离所述因此而产生的信号。这些角度可以从几何上通过形成具有长度α、β和|γ|的边的三角形而被找到。这样的三角形可以被形成，因为三个量满足以上“三角形不等式”。

[0091] 余弦定理陈述了，对于具有与长度c的边相对的角度C的任何三角形，其中其他边的长度为a和b：

[0092] c2＝a2+b2‑2ab cos C

[0093] 通过重新布置位，余弦定理被用来分别计算用于α和β所表示的信道的相位角A和B：

[0094]

[0095]

[0096] 信道a被分配相位角‑B，并且信道b被分配相位角A，其余的信道在粗略逼近步骤中被分配，即，分别地，根据 的对应元素是+1、还是‑1，要么被分配0度，要么被分配180度。

[0097] 上面描述的相同的方法也可以被用于Rmin<0。在这种情况下，计算是相同的，除了α和β被选择为 的(在幅值上)最大的负元素，并且γ是正的；并且信道a和b的相位角分别被分配180°‑B和180°+A。

[0098] 选择相位角的方法可以用数值和图形示例来例示说明。考虑N＝4的情况，即，4个作用的信道，其中振幅分别为1、1、3和2，即， 可以确定最小的可能的残差的幅值为1。 将产生 其中
(其他置换也将产生Rmin＝‑1或Rmin＝1，但是1对于计算
是足够的。)因此，粗略逼近步骤得到[180°,180°,0°,180°]的相位角。这在图6A中被用矢量描绘，其中 的元素在顶部被示出，负元素指向左，单个正元素指向右。‑1的残差(R_min)被描绘在下面。

[0099] 因为残差为负的，所以α和β被选为负方向上的两个最大的矢量，其中长度为α＝‑2和β＝‑1，分别地，对应于信道a＝4和b＝1。图6B示出在这2个矢量被移除之后的剩余的矢量，其中新的残差γ＝+2被描绘在下面。

[0100] 下一个步骤是调整矢量α和β的相位以便使γ偏移，即，以使得3个矢量将在矢量上相加为零。应用余弦定理，

[0101]

[0102]

[0103] 因此我们将信道a的相位设置为180°‑B＝180°‑28.96°＝151.04°，并且将信道b的相位设置为180°+A＝180°+75.52°＝255.52°。在图6C中，矢量α和β被示为被放置成这些角度，其中γ被放置成0°。角度A和B在由矢量α、β和γ形成的三角形中被指示。可以看出，α离指向右的(零相位)位置180°‑B，并且β离这样的位置180°+A。

[0104] 因此，通过所述方法得到的最终的相位角为[255.52°,180°,0°,155.04°]。在图6D中，单个的信道所表示的相量被示为具有这些相位角。4个相量的矢量和为零，从而在驱动线和返回线中得到零的净电流以及由换能器产生的电场的最大偏移。

[0105] 如果阵列上的超声换能器正被以多于一个的驱动频率驱动，则对于正被以给定的驱动频率驱动的换能器的每个子集执行本文中所公开的方法。因为振幅和驱动频率在超声过程的过程期间被调整和重新校准，所以每次用于信道中的任何一个的任何振幅和/或频率有改变时，重复相位角计算。对于最佳的EMI减小，每次如通过控制算法和/或处置计划确定的、频率和/或振幅有更新时，在振幅和/或频率的改变在发送到换能器的信号中被实现之前，与新的振幅和频率集合相对应的相位角被如本文中所公开的那样计算。然后，当这些信号被改变时，新的频率和振幅被与如本文中所公开的那样确定的新的相位角同时地实现。

[0106] 本文中所公开的方法在图7A‑7D中所描绘的流程图70中被图示说明。所述方法从步骤702开始，其中N矢量 被定义，该矢量包括将被发送到阵列中的相应的超声换能器的、给定驱动频率的N个电压信号的振幅中的每个。这样的振幅根据超声过程的目标并且根据接收的反馈和其他考虑来确定。这些振幅不是通过本文中所公开的方法确定的，而是根据对于这样的方法的输入来确定的。所述方法的输出是矢量 该矢量包括将被用于发送到每个换能器的电压信号的相位角。

[0107] 所述方法继续进行到步骤704，在步骤704中，计算 与每个可能的 矢量(也就是说，其元素要么为+1、要么为‑1的每个可能的N矢量)的点积(或标量积)。然后在步骤706，从步骤704的结果确定Rmin为所计算的在幅值上最小的点积，并且确定和选择矢量 对于该矢量 可能存在得到Rmin(或‑Rmin)的多于一个的可能的 矢量，在这种情况下，这样的可能的 矢量中的一个被任意地选择或者基于其他考虑(诸如具有其相位偏移相反的相邻信道以使电场抵消最大)被选择。注意，为了该流程图的目的，Rmin可以是正的或负的。如本文中其他地方所指出的，在一些实施方案中，Rmin和 可以通过在对 进行所有的可能的置换的情况下不包括计算所有的可能的点积的其他方法来确定。

[0108] 接着是步骤708，在步骤708中，确定最大的振幅(也就是说， 中的最大的元素)是否大于所有的振幅的和的一半。如果是，则在步骤710(参见图7D)，根据 设置相位角，也就是说， 其中反余弦是逐个元素获取的，以使得根据P的对应元素是+1、还是‑1，用于每个信道的相位角要么是0°，要么是180°。在一些实施方案中，步骤708在步骤
704之前发生，以使得如果到达步骤710，步骤704和706不是必要的。在这样的情况下， 将被设置为使得如果信道m是具有最大的振幅的信道，则对于所有的i≠m，θm＝180°，并且θi＝
0°。

[0109] 如果步骤708结果为“否”，则流程图继续进行到占位符A。图7B从占位符A开始，占位符A然后继续进行到步骤712，计算 并且在步骤714，流程图根据Rmin是非负的(“是”分支)、还是负的(“否”分支)，分裂为2个分支。在步骤720和730，将计算相位角中所用的标量量Θ要么设置为0°，要么设置为180°。步骤720和730然后分别继续进行到占位符C和D。

[0110] 图7C从占位符C和D开始。在步骤722和732，α被定义为 中的最大的正的(或负的，视情况而定)条目，其位置的索引为i(α＝Oi)。在步骤524和534，β被定义为 中的第二大的正的(或负的)条目，其位置的索引为j(β＝Oj)。

[0111] 流程图的分支在步骤540处重新会合，在步骤540中，残差γ被定义为移除Oi和Oj之后的、 中的剩余的条目的和。步骤540然后继续进行到占位符E，占位符E在图7D中也被找‑1 2 2 2到。在步骤542，使用余弦定理，从α、β和γ计算角度A和B。具体地说，A＝cos ([β+γ ‑α])/‑1 2 2 2
[2×|β|×|γ|])，并且B＝cos ([γ+α‑β])/[2×|γ|×|α|])。

[0112] 在步骤744，根据粗略逼近步骤，将相位角 的矢量暂时设置为 然后，在步骤746，根据在步骤542中确定的角度，设置信道i和j的相位角(也就是说，矢量α＝Oi和β＝Oj的相位角)，以使得如果Rmin为非负的，则θi被设置为‑B，或者如果Rmin为负的，则θi被设置为180°‑B；并且如果Rmin为非负的，则θj被设置为A，或者如果Rmin为负的，则θj被设置为
180°+A。注意，如果Rmin＝0，则步骤542中的计算将得到A＝B＝0，以使得相位角将全都要么为0°，要么为180°。

[0113] 从步骤746或步骤710，下一个步骤是步骤750，在步骤750中，信道1至N根据矢量被分配相位角，并且这样的相位角在发送到超声换能器的、由 给出的振幅的电压信号中被实现。流程图70可以按每个硬件更新间隔重复(即，返回到步骤502)，即，每次换能器中的任何一个在其被驱动的振幅和/或频率有任何改变时，流程图中所示的方法将被跟随，以便确定将与驱动信号一起使用的最佳相位角，以使得这样的相位角可以与振幅和/或频率的改变同步地被实现到可行的最大程度。

[0114] 本发明不应被认为限于上述特定实施方案。各种修改、等同过程、以及本发明可适用于的许多结构对于本发明针对的领域的技术人员来说在阅览本公开时将是容易明白的。