Optimization of Multiresonant Wireless Power Transfer Network Based on Generalized Coupled Matrix

Magnetic coupling resonant wireless power transfer network (MCRWPTN) system can realize wireless power transfer for some electrical equipment real-time and high efficiency in a certain spatial scale, which resolves the contradiction between power transfer efficiency and the power transfer distance of the wireless power transfer. A fully coupled resonant energy transfer model for multirelay coils and ports is established. A dynamic adaptive impedance matching control based on fully coupling matrix and particle swarm optimization algorithm based on annealing is developed for the MCRWPTN. Furthermore, as an example, the network which has twenty nodes is analyzed, and the best transmission coefficient which has the highest power transfer efficiency is found using the optimization algorithm, and the coupling constraints are considered simultaneously. Finally, the effectiveness of the proposed method is proved by the simulation results.


Introduction
Since wireless power transfer was firstly proposed by Nikola Tesla 100 years ago [1], much attempts have been devoted to this field.But these attempts did not get much breakthrough for several decades.Until 2006, a MIT research team found that efficient mid-range wireless power transfer could be achieved by utilizing near field magnetic coupling between two identical resonators [2,3].This progress excited renewed interest in this subject.
The magnetic coupling resonant wireless power transmission network (MCRWPTN) is to achieve wireless energy transmission through the resonant coils.In order to improve the transmission efficiency and power and adapt to the needs of the system, the MCRWPT has a variety of advanced multiresonator transmission structures, including relay resonator structure, multitransmit coil, or transmit coil array structure, multireceive coil, or multiple load structure.With the increase of the complexity of the multiresonator structure, a series of new design problems have been introduced in the MCRWPT system, which is a hot research direction in the field of MCRWPT technology.
From the mutual inductance formula, the transmission distance can be increased by increasing the size of the coil, but it is not feasible to apply the method to the limit of the coil geometry.Reducing the load of the transmitter and receiver can also increase the transmission distance of [4], but the actual application of the load cannot be arbitrarily changed; this method also has a lot of limitations.The increase of relay resonator is an effective way to improve the transmission distance.
The concept of relay coil was put forward in [5] based on the induction of electromagnetically induced transparency in quantum interference.By using the coupled mode theory, the energy transfer efficiency of the system was greatly improved by adding a coil with the same resonant frequency between the transmitting coil and the receiving coil.Further study in [6][7][8][9][10] showed that the higher transmission efficiency and power level can be achieved at a larger distance by increasing the relay resonator.
The effect of the addition of 1 relay coil on the transmission efficiency of the MCRWPT system with four coils was studied in [7], and the coaxial relay coil placed on the system to enhance the transmission efficiency was greater, but the vertical display had wider practical application prospects.The [8] proposed MCRWPT array resonator based on  resonator through the relay placed level could achieve  times the diameter of the resonator of wireless power transmission (WPT) distance, when using the 10-array resonator, the efficiency of energy transmission was still close to 85%.An efficient WPT system with a relay coil was designed in [10].Based on the analysis and simulation of the circuit model, the appropriate operating frequency was selected to achieve a transmission power of 6.6 kW at a distance of 200 mm and transmission efficiency was 95.57%.
In the practical application of radio energy transmission, such as charging for mobile electronic devices, there was a large demand for the supply of multiple devices simultaneously.With the development of WPT system, more and more attention had been paid to the research of multireceiving coil and multiple load structure.
The basic transmission characteristics such as transmission efficiency and power of multireceiving coil in [11][12][13] were studied by using circuit model or coupled mode theory.The research on the multireceiver coil structure in WPT was carried out in [11], which verified the feasibility of transmitting the power to a plurality of receiving coils.In a single transmitting coil to the two receiving coils power supply research in [12], the results show that although the transmission efficiency of each single coil is relatively not high, it overall can achieve higher transmission efficiency of the system.Based on the coupled mode theory, [13] was used to study the transmission of multiple load coils.With the same conclusion as in [12], with the two receiving coils, the overall transmission efficiency of multiple receive coils is higher than that of a single receiver coil, but as the number of receivers increases, the overall efficiency of the system tends to saturate and the efficiency of each individual receive coil becomes very low.In [14], the cross-coupling of multireceiving coils was considered, and the optimal design of multiload reception was realized by impedance matching of the transmitter.The paper [15] presented a novel analysis method using band-pass filter model, which made the circuit model with multireceiving coil system more simple and easier to analyze, but its use conditions were relatively harsh.
In [16], the concept of time-sharing transmission in communication system was used to study the multiload charging in MCRWPT system.It was found that time-sharing wireless transmission can achieve efficient and uniform power transmission.The literature [17] was further extended to  resonant coils.The simulation results showed that when there were  resonant coils in the system and the frequency splitting occurred, there would be  efficiency peak points.When there were odd resonant coils, the system transmission efficiency could reach the peak at the original resonant frequency of the coil, and when the number of resonant coils was even, the transmission efficiency was shifted at the original resonant frequency.
The impedance matching method on the load side can simultaneously realize the compensation of the frequency splitting phenomenon and the control of the power distribution.Considering the limitation of the frequency tracking method, the impedance matching is the main control strategy in the multicoil system.The impedance matching strategy for power allocation was studied in [15,18,19].In the actual system, the coupling cannot be ignored when the number of receiving coils is very close [15,18], considering the cross-coupling between multiple receiving coils.The impedance matching and power division method utilizing impedance inverters only at receiver sides was proposed.The mathematical equations in the proposed method were then generalized for arbitrary number of receivers and arbitrary number of repeaters.The coupling of the driving coil and the multiple receiving coils was studied in [19], and it was pointed out that the coupling will change the optimal impedance matching condition given the design method of matching impedance at this time.
Using the scattering parameters as the objective function of the fully coupling matrix extraction and combining the particle swarm optimization algorithm with simulated annealing technique to optimize impedance matching are the goal of this work.The control strategy must overcome some difficulties such as the highly nonlinear and coupled dynamics more over the dynamic, complex, and unstructured environments which may cause unpredictable disturbances to the impedance matching.So dealing with some states in each step will improve matching speed and accuracy.Using the coupling analysis and extraction condition of fully coupling matrix to design an optimal impedance matching strategy for MCRWPTN, it has been shown that the resulting optimization has excellent performance, as demonstrated by simulations.
The paper is outlined as follows.The model of multirelay MCRWPTN is established, and the coupling relation between multiple coils is analyzed in Section 2. The generalized coupling matrix of the system is extracted and an optimal design method is given in Sections 3 and 4. Simulations using FEKO and MATLAB and experiments are performed to validate the new method in Section 5. Finally, Section 6 presents some conclusions.

Problem Formulation
Wireless power transfer network (WPTN) is an energycentric network, which includes an exciting coil, a transmitting coil, a relay coil, a receiving coil, and a pickup coil.The generalized coupling model of the single energy transmission link composed of  relay energy nodes in the WPTN network is shown in Figure 1.
An internal resonance network is composed of a transmitting, a relay, and a receiving coil, and the resonant frequency and the structure of the transmitting coil, the relay coil, and the receiving coil are the same.The exciting coil and the pickup coil are connected to the power and load ports, respectively, and form an external resonant network with the internal resonant network.The magnetic coupling of the internal resonant network forms a contactless energy transmission link, and multiple transmission links constitute the entire WPTN.
In Figure 1,   ,   , and   ( = , ) are power and load port parameters;   ,   , and   are the parameters of each coupling node.  is the mutual inductance between the port and the node;  ,+1 is the mutual inductance between adjacent nodes.The  parameter equation of the two-port system can be expressed as

Relay resonant network
where [] is the system equivalent impedance matrix and [] is a fully coupling matrix.For systems with  ( = 1, . . ., ) resonant networks and  ports, the order of the equivalent impedance matrix [] is  + 2.  is the normalized frequency factor;  is the unit matrix. is the port matrix, except for  11 =   and  +2,+2 =   and other elements are zero.The impedance matrix of the system is shown in (3).Matrix [] is derived from the loop equation formulation of the source and load, external and inner-resonator couplings, and resonator impedances.
Without loss of generality, we can assume that all nodes are fixed and stationary randomly distributed within the network, and each node has a unique number.All nodes have the same structure and circuit topology, and the coupling between nonadjacent nodes is negligible.Each node has the same function and status equivalence as the source node, relay node, and destination node.
Generalize matrix [] such that the resultant matrix can be normalized in terms of couplings and port impedances only.The resultant matrix is the normalized impedance matrix [ Z], as shown in Normalized variables are used in representing the impedance matrix as follows.  =  0   /  is the unloaded  factor of resonator ,  = (/ 0 −  0 /) is the normalized frequency factor, which reflects the relationship between the resonant network and the frequency in the frequency response of MCRWPT.The coupling coefficient of the internal resonator is defined as  ,+1 =  ,+1 /√   +1 , and the external coupling between resonator  and port  for  = 1, and  = , is defined as  1 =  1 /√ 0    1 and   =   /√ 0     .The normalized reactance of port  for  = ,  is   =   /  and   =   /  .
This normalized impedance matrix will be referred to as [ Z]. Substituting in these variables allows us to make the decomposition of the matrix seen below: Using algebraic manipulation, [ Z] is separated into constitutive components into the form given by For the purposes of WPT optimization, the variable  simply reflects the frequency dependency of the MRCWPT frequency response.The unit matrix [] of the resonant network is  *  identity matrix, the port normalized impedance matrix [] is zero except for the power and load port elements, the coupling matrix [] is a reciprocal matrix of  * , and the other nondiagonal elements represent the coupling between the two resonators.If it is zero, the corresponding two resonators are not coupled.

Fully Coupling Matrix Extraction
From (1) we can see that the current vector can be expressed as Thus, the scattering parameter  of the entire coupling circuit can be expressed as We can get from ( 7) that where the matrix [] is built from the coupling matrix [], the source matrix [p], and the resonator matrix [].It can be shown that, from inverting the [] matrix, the scattering parameters can be found.
The objective function  is uniquely determined by the zero position, the pole position, and the in-band return loss.

𝐾 =
where   is artificially specified quantity and   is a function of   .  and   are the zero and pole of the characteristic function   (), respectively.The latter two are used to determine the in-band return loss.When  = 0, the target is satisfied and the optimization is finished.Set the following expression: 1 () can be obtained as substituting the first pole  1 into (14), and the values of  0 ,  1 , V 0 , V 1 are obtained by comparing the coefficients, and then determine  1 (),  1 ().
2 () can be obtained as substituting the first pole  2 , and so on, until the  pole is replaced.From formulas ( 16) and ( 17), the molecules of   () can be obtained by iteration.
Solving the roots of   () = 0 can obtain the  zeros of the characteristic function.The parameters of the objective function can be determined after substituting the zero point, the pole, and the return loss into (12)

Impedance Matching Optimization
In order to improve the extraction speed of generalized coupling matrix, a particle swarm optimization algorithm based on simulated annealing is presented.Fitness function can be expressed as follows: Firstly, the initial population is generated randomly, and a new group of individuals is generated by the basic particle swarm optimization algorithm.
where  is the best location of the particle itself and  is the best location for all particles in the population.
Then, the crossover and the Gaussian mutation operation are carried out independently.Simulated annealing is performed for each individual, and the result is used as an individual in the next generation.
In each evolution, crossover operation selects a specified number of particles into a pool according to crossing probability, which are randomly crossed to produce the same number of children, instead of parents to keep the number of particles unchanged.The position of the child particle is calculated by the weighted sums of the parent particles.
where  is the -dimension position vector and child  () and parent  (),  = 1, 2, respectively, are the location of the child or parent particles. is a uniform random number vector of -dimension, with each component in [0, 1] value.The speed of the child particles is given by the following formula: where V is the velocity vector of the -dimension.
In each evolution, the specified number of particles is selected by mutation operation according to the mutation probability, which was mutated by Gaussian mutation operator, and the mutated particles are used instead of the original particles.mutation () =  * (1 + Gaussian ()) . ( The implementation of the algorithm is composed of two parts.Firstly, the evolution of the basic particle swarm optimization algorithm (global search) is used to produce a better group, and then the crossover and mutation operations are applied, and the particles are further optimized by simulated annealing (local search).The evolution process is iterated until the termination condition is satisfied.The algorithm flow is as follows: (1) Initialization parameters.
(2) The population of randomly generated  particles.
(3) Using (24) to manipulate the particles in the population.(6) If the current optimal individual satisfies the convergence condition, the evolution process ends successfully and returns the global optimal solution.
(7) If the number of evolutions is less than the maximum number of times of evolution, the annealing temperature of the population is modified, and then return to step (3).

Simulation Results
In order to optimize the extraction of generalized coupling matrix parameters fast convergence, the initial value should be set close to the actual value.The parameters of the MCRWPTN system are shown in  When the number of ports is  = 2, the number of coupling relay coils is  = 20, and the system scattering parameter curve is shown in Figures 2-4.These three curves represent the specific distribution of the scattering parameters of the total transmission network and the two transmission subnetworks of the MCRWPTN, respectively, which can reflect the law of the scattering parameters of the system with frequency, where  21 is forward transmission coefficient and  11 is reflection coefficient of port 1.They can respond to the transmission characteristics of the system, the higher  21 is better, and the opposite of  11 is as low as possible.
As can be seen from Figure 2, the frequency division of the system has low frequency resonance and high frequency resonance.At the resonant frequency,  21 reaches the maximum, and  11 = 0.1 is the smallest.So the maximum power transmission efficiency is obtained.The response curve  obtained after optimization is in good agreement with the S parameters obtained from the simulation.From Figures 3 and  4, the existence of external coupling network makes the port impedance match well.
Figure 5 shows the simulated  11 and  22 responses of the complex load impedance accommodated response measured at the interresonator coupling location  ST = 0.01.As shown, the resonance dip of  22 response approaches the center of the Smith chart.This indicates that the output impedance of the impedance matching optimized system is the complex conjugate of   =   +   .
Figure 6 is the modulus value of the reflection coefficient obtained by different algorithms.By comparing the three different algorithms of reflection coefficient magnitude, the particle swarm optimization combined with simulated annealing algorithm is a new type of impedance matching algorithm, and the reflection coefficient tends to zero to achieve impedance matching finally, which is better than the improved particle swarm algorithm in speed.The relationship between the number of coupled nodes and the number of iterations is shown in Figure 7.As the number of nodes increases, the algorithm is superior to other algorithms and is easy to extend to multirelay MCRWPTN.

Conclusion
A fully coupling matrix for multirelay MCRWPTN is proposed and theory to impedance matching optimization algorithm is developed.This combination will emphasize not only the energy transfer but also the matching speed and scalability.The optimization algorithm combines the ability of the simulated annealing algorithm to jump out of the local optimal solution and the global optimization ability of the particle swarm optimization algorithm, which avoids the shortcomings of falling into the local extreme point and improves the algorithm convergence speed and accuracy.Theoretical results are supported by numerical simulations with impedance matching of the proposed fully coupled matrix and optimization algorithm.

Figure 1 :
Figure 1: Schematic of a wireless power transfer route.

Figure 3 :
Figure 3: The curves of the scattering parameters of external coupling network of power port.

Figure 5 :
Figure 5: Simulation validation of complex load impedance accommodation.

Figure 6 :Figure 7 :
Figure 6: The modulus value of the reflection coefficient obtained by different algorithms.
. The relation between generalized coupling matrix [] and objective function  can be established by the above derivation, and the optimization model of the coupling matrix extraction is shown in formula (22).

Table 1 :
The system parameters configuration.

Table 1 .
The optimized and simulated  parameters are shown in Figures2-4.