Magnetic Stimulation Method with Controllable Induced Field Direction

Abstract

The disclosure discloses a magnetic stimulation method with a controllable induced field direction, and belongs to the technical field of noninvasive neural regulation. The method includes the following steps: S100, calculating currents i.sub.1j, i.sub.2j, i.sub.3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point P.sub.t, the same below; S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E.sub.1, E.sub.2, E.sub.3; S300, calculating currents that may generate the electric field components E.sub.1, E.sub.2, E.sub.3 at the target point, I.sub.1j=E.sub.1i.sub.1j, I.sub.2j=E.sub.2i.sub.2j, I.sub.3j=E.sub.3i.sub.3j; S400, superimposing the currents of the three energization modes to obtain a resultant current I.sub.j=I.sub.1j+I.sub.2j+I.sub.3j=E.sub.1i.sub.1j+E.sub.2i.sub.2j+E.sub.3i.sub.3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.

Claims

1. A magnetic stimulation method with a controllable induced field direction, the magnetic stimulation method with the controllable induced field direction comprising the following steps: S100, calculating currents i.sub.1j, i.sub.2j, i.sub.3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point P.sub.t, the same below; S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E.sub.1, E.sub.2, E.sub.3; S300, calculating currents that may generate the electric field components E.sub.1, E.sub.2, E.sub.3 at the target point, I.sub.1j—E.sub.1i.sub.1j, I.sub.2j=E.sub.2i.sub.2j, I.sub.3j=E.sub.3i.sub.3j; S400, superimposing the currents of three energization modes to obtain a resultant current I.sub.j=I.sub.1j+I.sub.2j+I.sub.3j=E.sub.1i.sub.1j+E.sub.2i.sub.2j+E.sub.3i.sub.3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.

2. The magnetic stimulation method with the controllable induced field direction according to claim 1, wherein before S100, the method further comprises: S000, constructing a magnetic stimulation coil group.

3. The magnetic stimulation method with the controllable induced field direction according to claim 2, wherein S000 specifically comprises the following steps: S001, constructing a three-dimensional space; S002, placing the target point P.sub.t in the three-dimensional space; setting q lead wires, end points of ends of which are intersected, above the target point P.sub.t, wherein all the lead wires form an included angle with an xoy plane and form included angles θ.sub.1, θ.sub.2, . . . θ.sub.q with the forward direction of the z axis; S003, closing the other ends of the q lead wires with the intersected end points by line segments or polygonal connecting lead wires to form a magnetic stimulation coil group which has n coils.

4. The magnetic stimulation method with the controllable induced field direction according to claim 1, wherein in S100, a current is made into the magnetic stimulation coil group by the following energization modes: in a first energization mode, currents in all the coils are I.sub.11, I.sub.12, I.sub.13, . . . , I.sub.1j, . . . , I.sub.1n, thus generating an induced electric field E.sub.1; in a second energization mode, currents in all the coils are I.sub.21, I.sub.22, I.sub.23, . . . , I.sub.2j, . . . , I.sub.2n, thus generating an induced electric field E.sub.2; in a third energization mode, currents in all the coils are I.sub.31, I.sub.32, I.sub.33, . . . , I.sub.3j, . . . , I.sub.3n, thus generating an induced electric field E.sub.3.

5. The magnetic stimulation method with the controllable induced field direction according to claim 4, wherein in S200, specifically: the coil group in each energization mode respectively generates combined induced electric fields E.sub.1(a.sub.1,b.sub.1,c.sub.1), E.sub.2(a.sub.2,b.sub.2,c.sub.2), E.sub.3(a.sub.3,b.sub.3,c.sub.3) at the target point P.sub.t below coil intersections; three induced electric fields are not coplanar, that is, E.sub.1, E.sub.2, E.sub.3 are linearly independent, wherein $.Math.\begin{array}{ccc}{a}_1& b_1& c_1\\ a_2& b_2& c_2\\ a_3& b_3& c_3\end{array}.Math.\ne 0.$

6. The magnetic stimulation method with the controllable induced field direction according to claim 5, wherein in S300, specifically: obtaining corresponding unit-direction vector electric fields, e.sub.1=E.sub.1/|E.sub.1|, e.sub.2=E.sub.2/|E.sub.2|, e.sub.3=E.sub.3/|E.sub.3|, and generating currents needing to be made into the unit-direction vector electric fields at the target point,
i.sub.11=I.sub.11/|E.sub.1|,i.sub.12=I.sub.12/|E.sub.1|,i.sub.13=I.sub.13/|E.sub.1|, . . . ,i.sub.1j=I.sub.1j/|E.sub.1|, . . . ,i.sub.1n=I.sub.1n/|E.sub.1|;
i.sub.21=I.sub.21/|E.sub.2|,i.sub.22=I.sub.22/|E.sub.2|,i.sub.23=I.sub.23/|E.sub.2|, . . . ,i.sub.2j=I.sub.2j/|E.sub.2|, . . . ,i.sub.2n=I.sub.2n/|E.sub.2|;
i.sub.31=I.sub.31/|E.sub.3|,i.sub.32=I.sub.32/|E.sub.3|,i.sub.33=I.sub.33/|E.sub.3|, . . . ,i.sub.3j=I.sub.3j/|E.sub.3|, . . . ,i.sub.3n=I.sub.3n/|E.sub.3|;

7. The magnetic stimulation method with the controllable induced field direction according to claim 6, wherein in S400, specifically: according to the superposition principle of fields, three directional fundamental vector electric fields at the target point are used to combine a unit electric field e in arbitrary direction in a space,
e=λ.sub.1e.sub.1+λ.sub.2e.sub.2+λ.sub.3e.sub.3,
0≤λ.sub.i≤1 correspondingly, the currents in all the coils are also superposed, $i_1={\lambda }_1{i}_{11}+{\lambda }_2{i}_{21}+{\lambda }_3{i}_{31}{i}_2={\lambda }_1{i}_{12}+{\lambda }_2{i}_{22}+{\lambda }_3{i}_{32}{i}_3={\lambda }_1{i}_{13}+{\lambda }_2{i}_{23}+{\lambda }_3{i}_{33}.Math.i_j={\lambda }_1{i}_{1j}+{\lambda }_2{i}_{2j}+{\lambda }_3{i}_{3j}.Math.i_n={\lambda }_1{i}_n+{\lambda }_2{i}_{2n}+{\lambda }_3{i}_{3n}$ a ratio i.sub.1:i.sub.2:i.sub.3: . . . :i.sub.j: . . . :i.sub.n of the currents made into the coils in the coil group is adjusted, that is, the direction of the electric field e at the target point is controlled; the sizes ki.sub.1, ki.sub.2, ki.sub.3, . . . :ki.sub.j, . . . :ki.sub.n of single currents are overall adjusted, that is, the size ke of the electric field at the target point is adjusted.

Description

BRIEF DESCRIPTION OF FIGURES

[0030] FIG. 1A-FIG. 1C are projections of a coil center structure on an xoy plane;

[0031] FIG. 2 is a diagram of a coil center structure, which is composed of q wires, q=1, 2, 3, . . . n;

[0032] FIG. 3 is an example diagram of making a current;

[0033] FIG. 4 illustrates magnetic stimulation coil groups in four rectangular forms;

[0034] FIG. 5A-FIG. 5D illustrate an implementation method of controlling a direction of an induced electric field, where FIG. 5A is an energization mode 1; FIG. 5B is an energization mode 2; FIG. 5C is an energization mode 3; FIG. 5D is an energization mode 4;

[0035] FIG. 6 is a diagram of electric field distribution on a section under various energization modes;

[0036] FIG. 7A-FIG. 7D illustrate electric field distributions on a section when a resultant current is made and a size of a field value of an electric field of a target point, where FIG. 7A is E; FIG. 7B is an x-directional component of E; FIG. 7C is a y-directional component of E; FIG. 7D is a z-directional component of E;

[0037] FIG. 8A-FIG. 8E are example diagrams of five magnetic stimulation coil groups, where FIG. 8A illustrates four quadrilateral forms; FIG. 8B illustrates four triangular forms; FIG. 8C illustrates three triangular forms; FIG. 8D illustrates five triangular forms; FIG. 8E illustrates three hexagonal forms.

DETAILED DESCRIPTION

[0038] The technical solutions in the examples of the present disclosure will be clearly and completely described below in conjunction with the accompanying drawings in the examples of the present disclosure. Apparently, the described examples are only a part of the examples of the present disclosure, rather than all the examples. Based on the examples in the present disclosure, all other examples obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of the present disclosure.

[0039] Referring to FIG. 1A-FIG. 1C, the disclosure provides a magnetic stimulation method with a controllable induced field direction, the magnetic stimulation method with the controllable induced field direction including the following steps:

[0040] S100, calculating currents i.sub.1j, i.sub.2j, i.sub.3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point P.sub.t, the same below;

[0041] S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E.sub.1, E.sub.2, E.sub.3;

[0042] S300, calculating currents that may generate the electric field components E.sub.1, E.sub.2, E.sub.3 at the target point, I.sub.1j=E.sub.1i.sub.1j, I.sub.2j=E.sub.2i.sub.2j, I.sub.3j=E.sub.3i.sub.3j;

[0043] S400, superimposing the currents of the three energization modes to obtain a resultant current I.sub.j=I.sub.1j+I.sub.2j+I.sub.3j=E.sub.1i.sub.1j+E.sub.2i.sub.2j+E.sub.3i.sub.3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.

[0044] Further, before S100, the method further includes:

[0045] S000, constructing a magnetic stimulation coil group.

[0046] Further, in S000 specifically includes the following steps:

[0047] S001, constructing a three-dimensional space;

[0048] S002, placing the target point P.sub.t in the three-dimensional space; setting q lead wires, the end points of ends of which are intersected, above the target point P.sub.t, where all the lead wires form an included angle with an xoy plane and form included angles θ.sub.1, θ.sub.2, . . . θ.sub.q with the forward direction of the z axis;

[0049] S003, closing the other ends of the q lead wires with the intersected end points by line segments or polygonal connecting lead wires to form a magnetic stimulation coil group which has n coils.

[0050] Specifically, the physical principle of a magnetic intervention technology is Maxwell's equations. A time-varying electromagnetic field generated by a coil pulse current can be regarded as a magnetic quasi-static field. An induced electric field in an air domain has the following relationship with the current:

[00003]$E=-\frac{\partial A}{\partial t}A=\frac{{\mu }_0}{4\pi }\int \frac{\mathrm{Idl}}{r}E=-\frac{{\mu }_0}{4\pi }\frac{\partial I}{\partial t}\int \frac{\mathrm{dl}}{r}$

[0051] where A is a vector magnetic potential; E is the induced electric field; IdI is a current element; r is a distance from the current element to a point to be determined; and μ.sub.0 is a vacuum permeability.

[0052] It can be known from the above formula that the induced electric field is parallel to the current element IdI, and has an opposite direction to that of the current element IdI. Based on the relationship between the current and the direction of an electric field and the superposition principle of fields, it is easy to construct a coil group structure that is capable of controlling arbitrary direction of an electric field in a space above a target point P.sub.t through intersecting currents. In order to obtain an electric field in arbitrary direction in a three-dimensional space, three non-coplanar fundamental vectors are required, which provide electric field components in the x, y, and z directions of the Cartesian coordinate system.

[0053] Further, in S100, a current is made into the magnetic stimulation coil group by the following energization modes:

[0054] in a first energization mode, currents in all the coils are I.sub.11, I.sub.12, I.sub.13, . . . , I.sub.11, . . . , I.sub.1n, thus generating an induced electric field E.sub.1;

[0055] in a second energization mode, currents in all the coils are I.sub.21, I.sub.22, I.sub.23, . . . , I.sub.2j, . . . , I.sub.2n, thus generating an induced electric field E.sub.2;

[0056] in a third energization mode, currents in all the coils are I.sub.31, I.sub.32, I.sub.33, . . . , I.sub.3j, . . . , I.sub.3n, thus generating an induced electric field E.sub.3.

[0057] Specifically, an electric field in arbitrary direction is constructed on a two-dimensional plane. FIG. 1 shows a projection of a center structure of a coil group on an xoy plane. An intersection of lead wires is located directly above the target point. The lead wire structures alone or that overlap into a straight line cannot bring electric field components in the x and y directions, and it is impossible to construct, on the plane, a coil structure that controls the direction of the electric field by a current. The lead wire structures that radiate 2, 3, . . . , q lead wires from the center of the target point can bring components in the x and y directions on the plane, thus combining an electric field in arbitrary direction in the xoy plane.

[0058] In order to enable the coil to be able to generate an electric field component in the z direction for the electric field in arbitrary direction in the three-dimensional space, the lead wires need to form a certain included angle with the xoy plane. Therefore, the coil structure is a three-dimensional structure rather than a planar structure. As shown in FIG. 2, included angles between the lead wires and the forward direction of the z axis are θ.sub.1, θ.sub.2, . . . θ.sub.q, respectively.

[0059] The lead wires are closed with polygonal connecting line segments, so that a coil group structure with a controllable induced field direction is obtained. The number q of the intersecting lead wires above the target point corresponds to the number n of coil groups.

[0060] The implementation principle of the magnetic stimulation method with the controllable induced field direction is as follows: three energization modes are used to make a current into the coil groups. In a first energization mode, the currents in all coils are I.sub.11, I.sub.12, I.sub.13, . . . , I.sub.1j, . . . , I.sub.1n; in a second energization mode, the currents in all the coils are I.sub.21, I.sub.22, I.sub.23, . . . , I.sub.2j, . . . , I.sub.2n; in a third energization mode, the currents in all the coils is I.sub.31, I.sub.32, I.sub.33, . . . , I.sub.3j, . . . , I.sub.3n (FIG. 3 is taken as an example, m=1, 2, 3 represent the energization modes).

[0061] Further, in S200, specifically:

[0062] the coil group in each energization mode respectively generates combined induced electric fields E.sub.1(a.sub.1,b.sub.1,c.sub.1), E.sub.2(a.sub.2,b.sub.2,c.sub.2), E.sub.3(a.sub.3,b.sub.3,c.sub.3) at the target point P.sub.t below coil intersections; the three induced electric fields are not coplanar, that is, E.sub.1, E.sub.2, E.sub.3 are linearly independent, where

[00004]$.Math.\begin{array}{ccc}{a}_1& b_1& c_1\\ a_2& b_2& c_2\\ a_3& b_3& c_3\end{array}.Math.\ne 0.$

[0063] Further, in S300, specifically:

[0064] obtaining corresponding unit-direction vector electric fields, e.sub.1=E.sub.1/|E.sub.1|, e.sub.2=E.sub.2/|E.sub.2|, e.sub.3=E.sub.3/|E.sub.3|, and generating currents needing to be made into the unit-direction vector electric fields at the target point,

i.sub.11=I.sub.11/|E.sub.1|,i.sub.12=I.sub.12/|E.sub.1|,i.sub.13=i.sub.13/|E.sub.1|, . . . ,i.sub.1j=I.sub.1j/|E.sub.1|, . . . ,i.sub.1n=I.sub.1n/|E.sub.1|;

i.sub.21=I.sub.21/|E.sub.2|,i.sub.22=I.sub.22/|E.sub.2|,i.sub.23=i.sub.23/|E.sub.2|, . . . ,i.sub.2j=I.sub.2j/|E.sub.2|, . . . ,i.sub.2n=I.sub.2n/|E.sub.2|;

i.sub.31=I.sub.31/|E.sub.3|,i.sub.32=I.sub.32/|E.sub.3|,i.sub.33=i.sub.33/|E.sub.3|, . . . ,i.sub.3j=I.sub.3j/|E.sub.3|, . . . ,i.sub.3n=I.sub.3n/|E.sub.3|;

[0065] Further, in S400, specifically: according to the superposition principle of fields, three directional fundamental vector electric fields at the target point are used to combine a unit electric field e in arbitrary direction in a space,

e=λ.sub.1e.sub.1+λ.sub.2e.sub.2+λ.sub.3e.sub.3,

0≤λ.sub.i≤1

[0066] correspondingly, the currents in all the coils are also superposed,

[00005]$i_1={\lambda }_1{i}_{11}+{\lambda }_2{i}_{21}+{\lambda }_3{i}_{31}{i}_2={\lambda }_1{i}_{12}+{\lambda }_2{i}_{22}+{\lambda }_3{i}_{32}{i}_3={\lambda }_1{i}_{13}+{\lambda }_2{i}_{23}+{\lambda }_3{i}_{33}.Math.i_j={\lambda }_1{i}_{1j}+{\lambda }_2{i}_{2j}+{\lambda }_3{i}_{3j}.Math.i_n={\lambda }_1{i}_n+{\lambda }_2{i}_{2n}+{\lambda }_3{i}_{3n}$

[0067] a ratio i.sub.1:i.sub.2:i.sub.3: . . . :i.sub.j: . . . :i.sub.n of the currents made into the coils in the coil group is adjusted, that is, the direction of the electric field e at the target point is controlled; the sizes ki.sub.1, ki.sub.2, ki.sub.3, . . . :ki.sub.j, . . . :ki.sub.n of single currents are overall adjusted, that is, the size ke of the electric field at the target point is adjusted.

[0068] As shown in FIG. 8, the coil group can have various forms. A quadrilateral form including four coils in an air domain is taken as an example to describe the implementation of the magnetic stimulation method for controlling the direction of an induced electric field. The schematic diagram of a coil group structure is shown in FIG. 4. The coil group is composed of 4 identical rectangles of a “book-shaped” structure, and edge lengths of the rectangles are a and b. Coils 1 and 2 are coplanar; coils 3 and 4 are coplanar and two coplanar coilsform a wing. An included angle between the plane of each wing coil and a negative direction of the z axis is α. The coil group is of a symmetrical structure, that is, an included angle between two wings is 2α.

[0069] In energization mode 1, the currents made into the coil group are I.sub.11, I.sub.12, I.sub.13, I.sub.14, and the magnitudes of the currents are equal: I.sub.11=I.sub.12=I.sub.13=I.sub.14; the current directions are clockwise, clockwise, counterclockwise, and counterclockwise when seen from the forward direction of the z axis. In energization modes 2 and 3, the magnitudes of the currents made into the coil group are also the same: I.sub.21=I.sub.22=I.sub.23=I.sub.24, I.sub.31=I.sub.32=I.sub.33=I.sub.34, and the current directions are as shown in FIG. 5.

[0070] In the air domain, the direction of the induced electric field is opposite to that of the current. FIG. 6 shows an electric field distribution on a certain section below the coils under the three energization modes. The magnitude of the current is 1 A; the frequency is 2500 Hz; the side lengths a and b are both 70 mm; the included angle is 60°. Table 1 shows electric field values at the target point. The three energization modes only bring electric fields in the x, y, and z directions at the target point.

TABLE-US-00001 Energization I.sub.m1, I.sub.m2, I.sub.m3, I.sub.m4 E.sub.x E.sub.y E.sub.z mode (A) (V/m) (V/m) (V/m) 1 1, 1, −1, −1 −1.50 × 10.sup.−3 0 0 2 1, −1, 1, −1 0 2.44 × 10.sup.−4 0 3 1, −1, −1, 1 0 0 8.61 × 10.sup.−4 [0071] Table 1 The induced electric field values of the target point (0, 0, −0.08) under all the energization modes

[0072] Note: m=1, 2, 3 are the energization modes; the current is clockwise positive and counterclockwise negative when seen from the forward direction of the z axis, and the same below.

[0073] Therefore, the current required to generate a unit fundamental vector electric field at the target can be calculated, as shown in Table 2.

TABLE-US-00002 Energization e.sub.x e.sub.y e.sub.z i.sub.m1, i.sub.m2, i.sub.m3, i.sub.m4 mode (V/m) (V/m) (V/m) (A) 1 1 0 0 −0.67 × 10.sup.3, −0.67 × 10.sup.3, 0.67 × 10.sup.3, 0.67 × 10.sup.3 2 0 1 0 4.10 × 10.sup.3, −4.10 × 10.sup.3, 4.10 × 10.sup.3, −4.10 × 10.sup.3 3 0 0 1 1.16 × 10.sup.3, −1.16 × 10.sup.3, −1.16 × 10.sup.3, 1.16 × 10.sup.3 [0074] Table 2 Current required to generate a unit fundamental vector electric field at the target point (0, 0, −0.08)

[0075] A vector electric field E(2, 1, −1) V/m generated at the target point is taken as an example. Its components in a fundamental vector direction are E.sub.1=2, E.sub.2=1, E.sub.3=−1, I.sub.j=2i.sub.1j+i.sub.2j−i.sub.3j (j=1, 2, 3, 4); a resultant current needing to be made to the coil group is (I.sub.1, I.sub.2, I.sub.3, I.sub.4)=(1.61×10.sup.3, −4.27×10.sup.3, 6.59×10.sup.3, −3.93×10.sup.3) A.

[0076] FIG. 7A˜D shows an electric field distribution of a section including a target point when a resultant current is made. An electric field value generated by the target point (0, 0, −0.08) is consistent with an expected value.

[0077] In a computational domain with a complex medium and a complex boundary, a finite element calculation method can be used to calculate an induced electric field, or actual measurement can be carried out to calibrate a fundamental vector of an electric field and a current made into it.

[0078] The above implementation examples are only used to help understand the method of the disclosure and its core idea. For those skilled in the art, according to the idea of the disclosure, several improvements and modifications can be made in the specific implementation and application scope. These improvements and modifications shall also fall within the protection scope of the disclosure.