Ultrasonic focusing method and system based on quantum whale mechanism

Abstract

The invention relates to an ultrasonic focusing method and system based on a quantum whale mechanism, and belongs to the technical field of focused ultrasound. The method mainly comprises the steps that firstly, an error curve table is constructed, and a target optimization function is designed; initializing the quantum position, calculating the fitness of the quantum position, and determining a global quantum position; the quantum position of each whale is updated through iteration, the fitness function value is recalculated, and then the global quantum position is updated; and after the maximum number of iterations is reached, selecting an optimal quantum position, and mapping the optimal quantum position as a global optimal position so as to determine an azimuth weighting coefficient and an axial weighting coefficient of the target. Compared with a traditional table look-up method only using the maximum wave beam and the second maximum wave beam, the method has the advantages that the weighting coefficients in the azimuth direction and the axial direction can be more efficiently obtained through the quantum white whale mechanism, and the information of the five wave beams is fully utilized, so that the focusing precision is remarkably improved.

Description

Technical Field

[0001] The present invention relates to the field of ultrasonic focusing, and in particular to an ultrasonic focusing method and system based on a quantum white whale mechanism. Background Art

[0002] The phased array ultrasound focusing system uses a large number of individually controlled small transducer units arranged in an array. Each transducer unit is managed by an independent delay controller. By precisely controlling the phase of the ultrasound waves emitted by each transducer unit, focused beams in different directions can be synthesized. These ultrasound waves form a strong focus point at a specific location based on the principle of interference, while weaker side lobes are formed elsewhere. At the same aperture and operating frequency, the response speed, focus switching rate, multi-focus generation capability, resolution, and versatility of the phased array ultrasound system are superior to traditional mechanical scanning ultrasound equipment.

[0003] As a key equipment in modern ultrasound technology, it is crucial for the phased array ultrasound focusing system to accurately locate, focus and treat the target area. In particular, ultrasound focusing, that is, determining the accuracy of the ultrasound beam pointing, has a decisive influence on the performance indicators of the entire phased array ultrasound system, especially in applications that require high-precision guidance, such as medical imaging or focused ultrasound surgery (FUS).

[0004] In the related art, phased array ultrasonic focusing systems mostly use methods based on sum and difference beam amplitude ratio or adjacent beam comparison to perform focusing. Although these methods have the advantages of real-time and low computational complexity, their focusing accuracy may be affected by the beam broadening effect and do not make full use of all available beam information, so they may not provide ideal accuracy in some applications. Summary of the invention

[0005] In view of the shortcomings of the prior art, the present invention provides an ultrasonic focusing method and system based on the quantum white whale mechanism, which solves the problem of unsatisfactory multi-beam focusing accuracy of ultrasonic phased array due to factors such as beam broadening effect and low beam information utilization.

[0006] To achieve the above objectives, the present invention is implemented through the following technical solutions:

[0007] In a first aspect, an embodiment of the present application provides an ultrasonic focusing method based on a quantum white whale mechanism, the method comprising: determining the distribution information of a plurality of ultrasonic array elements in an ultrasonic phased array, calculating the azimuthal beam width and the axial beam width of the ultrasonic phased array in the uv space, and determining the beam overlap information to establish an error mapping table; acquiring transient signal data of a received signal, and constructing a target optimization function based on the error mapping table; initializing and setting parameters of the quantum phase delay of each white whale in a preset quantum white whale mechanism, calculating the fitness values ​​of all white whales through the constructed fitness function, and obtaining a global quantum phase delay; based on The Beluga optimization mechanism iteratively updates the quantum phase delay of each beluga to obtain the target quantum phase delay; maps the target quantum phase delay to the physical delay, calculates the fitness function values ​​of all beluga whales after the update, and updates the global quantum phase delay; when the maximum number of iterations is reached, determines the current global quantum phase delay as the optimal quantum phase delay, and transforms the optimal quantum phase delay mapping into the global optimal phase delay; based on the global optimal phase delay, determines the target azimuth weighting coefficient and the target axial weighting coefficient, and weights the azimuth focusing results and axial focusing results respectively to obtain accurate focusing information.

[0008] According to the first aspect of the embodiment of the present application, in an ultrasonic phased array, a plurality of ultrasonic array elements are respectively distributed in azimuth and axial directions. The aforementioned determination of the distribution information of a plurality of ultrasonic array elements in the ultrasonic phased array, calculation of the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and determination of the beam overlap information to establish an error mapping table may specifically include the following steps: determining the horizontal spacing of adjacent ultrasonic array elements distributed in the azimuth direction, and the axial spacing of adjacent ultrasonic array elements distributed in the axial direction; based on the horizontal spacing and the axial spacing, calculating the azimuth width and axial beam width of the ultrasonic phased array: mapping the azimuth width and the axial beam width to the uv space to obtain the azimuth beam width and the axial beam width; based on the azimuth beam width and the axial beam width, determining the error functions corresponding to the azimuth and axial directions, and establishing an error mapping table.

[0009] According to the first aspect of the embodiment of the present application, the aforementioned initialization and parameter setting of the quantum phase delay of each beluga in the preset quantum beluga mechanism, the fitness values ​​of all beluga are calculated by constructing a fitness function to obtain a global quantum phase delay, which may specifically include the following steps: setting the number of beluga to Q, the maximum number of iterations to G, and in the g-th iteration, the quantum phase delay of the q-th beluga in the 2-dimensional search space is When g=1, each dimension of the quantum phase delay of all beluga whales in the first generation is initialized to a uniform random number between [0, 1]; according to the preset target mapping rule, each dimension of the quantum phase delay of all beluga whales is mapped into the search space range to obtain the physical delay of the qth beluga whale in the gth generation; the fitness value of each beluga whale is calculated according to the constructed fitness function, and the quantum phase delay corresponding to the optimal fitness value is determined as the global quantum phase delay corresponding to the gth generation

[0010] According to the first aspect of the embodiment of the present application, the aforementioned target mapping rule satisfies the expression:

[0011]

[0012] In the formula, is the physical delay of the qth beluga whale of the gth generation, is the quantum phase delay of the qth beluga whale of the gth generation, u o，max and u o，min are the upper and lower limits of the o-th dimension search space, q=1, 2,…, Q, 0=1, 2; q and q correspond to the azimuthal weighting coefficient β and the axial weighting coefficient α respectively.

[0013] According to the first aspect of the embodiment of the present application, the aforementioned iterative update of the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain the target quantum phase delay may specifically include the following steps: in the iterative update process, for the qth beluga whale of the gth generation, a balance factor is generated; based on the relationship between the balance factor and the preset critical threshold, the target quantum beam phase delay of the oth dimension of the qth beluga whale of the g+1th generation is determined; according to the initial quantum beam phase delay, the quantum phase delay is updated using a simulated quantum revolving gate to obtain the corresponding target quantum phase delay:

[0014] Among them, when the balance factor is less than the preset critical threshold, the o-dimensional target quantum beam phase delay of the q-th beluga whale of the g+1th generation satisfies the expression: in, is the o-th dimensional quantum phase delay of the x-th beluga whale of the g-th generation, x is a random integer between [1, Q] except q, and is a random number uniformly distributed between [0, 1];

[0015] When the balance factor is greater than or equal to the preset critical threshold, the target quantum beam phase delay of the o-th dimension of the q-th beluga whale of the g+1th generation satisfies the expression:

[0016]

[0017] in, For the gth generation The oth-dimensional quantum phase delay of a white whale, is a random integer between [1, Q] except q, and is a random number uniformly distributed between [0, 1]. is the Levy flight coefficient, and is a random number that obeys a Gaussian distribution with a mean of 0 and a variance of 1, μ is a control parameter between [0, 2], and Γ represents the Gamma function.

[0018] According to the first aspect of the embodiment of the present application, the balance factor is in is a random number uniformly distributed between [0, 1];

[0019] When the balance factor is less than the preset critical value, the target quantum phase delay is Where, q = 1, 2, ..., Q, o = 1, 2;

[0021] When the balance factor is greater than or equal to the preset critical threshold, the target quantum phase delay is in, is the o-th dimension global optimal quantum phase delay until the g-th generation; the target quantum phase delay is

[0022] According to the first aspect of the embodiment of the present application, the aforementioned mapping of the target quantum phase delay to a physical delay, calculating the fitness function values ​​of all beluga whales after the update, and updating the global quantum phase delay can specifically include the following steps: mapping the target quantum phase delay to a physical delay, and calculating the fitness function values ​​of all beluga whales; determining the o-th dimensional quantum beam phase delay of the q-th beluga whale of the g+1th generation based on the fitness function value of the q-th beluga whale: based on the o-th dimensional quantum beam phase delay of the q-th beluga whale of the g+1th generation, using a simulated quantum rotating gate to update the quantum phase delay to obtain the quantum phase delay of the g+1th generation; updating the global quantum phase delay based on the quantum phase delay of the g+1th generation.

[0023] In the second aspect, an embodiment of the present application provides an ultrasonic focusing system based on a quantum white whale mechanism, the system comprising a calculation module, an optimization function construction module, a parameter configuration module, an iteration module, a first mapping module, a second mapping module and a weighting module; wherein the calculation module is used to determine the distribution information of multiple ultrasonic array elements in an ultrasonic phased array, calculate the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and determine the beam overlap information to establish an error mapping table; the optimization function construction module is used to obtain transient signal data of a received signal, and construct a target optimization function based on the error mapping table; the parameter configuration module is used to initialize and set parameters of the quantum phase delay of each white whale in the preset quantum white whale mechanism, and calculate all through the constructed fitness function. The fitness value of the beluga whale is used to obtain the global quantum phase delay; the iteration module is used to iteratively update the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain the target quantum phase delay; the first mapping module is used to map the target quantum phase delay to the physical delay, calculate the fitness function values ​​of all beluga whales after the update, and update the global quantum phase delay; the second mapping module is used to determine the current global quantum phase delay as the optimal quantum phase delay when the maximum number of iterations is reached, and transform the optimal quantum phase delay mapping into the global optimal phase delay; the weighting module is used to determine the target azimuth weighting coefficient and the target axial weighting coefficient based on the global optimal phase delay, and weight the azimuth focusing results and the axial focusing results respectively to obtain accurate focusing information.

[0024] The present invention provides an ultrasonic focusing method and system based on the quantum white whale mechanism. Compared with the prior art, it has the following beneficial effects:

[0025] The present application calculates the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and establishes an error mapping table in the uv space to cope with different situations where the central beam points; linear interpolation is used to reduce the quantization error when looking up the table, and weighted processing is performed on the azimuth focusing result and the axial focusing result respectively from the azimuth and axial directions; the present application is based on the quantum white whale mechanism, and through initialization of quantum phase delay, parameter setting, calculation of fitness values ​​of all white whales, iterative update of quantum phase delay of each white whale, phase delay mapping and other processing, finally when the maximum number of iterations is reached, the optimal quantum phase delay is determined and mapped to the global optimal phase delay, and the global optimal phase delay is used to characterize the estimated values ​​of the azimuth weighting coefficient and the axial weighting coefficient. The present application quickly obtains the target azimuth weighting coefficient and the target axial weighting coefficient through the quantum white whale mechanism. Compared with the traditional method of looking up the table using only the largest beam and the second largest beam in the azimuth and axial directions, the method of the present application fully utilizes the information of the five beams and has higher focusing accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the related technologies, the drawings required for use in the embodiments or the related technical descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0027] Figure 1 It is a schematic flow chart of an ultrasonic focusing method based on the quantum white whale mechanism provided by the present invention;

[0028] Figure 2 yes Figure 1 An example flow chart of S110;

[0029] Figure 3 It is a schematic diagram of the structure of an ultrasonic focusing system based on the quantum white whale mechanism provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0030] In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in this application will be clearly and completely described below in conjunction with the drawings in this application. Obviously, the described embodiments are part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0031] According to an embodiment of the first aspect of the present application, Figure 1 As shown, a schematic flow chart of an ultrasonic focusing method based on the quantum white whale mechanism provided in an embodiment of the present application, the ultrasonic focusing method may include the following steps S110-S170.

[0032] S110, determining distribution information of multiple ultrasonic array elements in the phased array, calculating the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and determining beam overlap information to establish an error mapping table.

[0033] S120, acquiring transient signal data of a received signal, and constructing a target optimization function based on an error mapping table.

[0034] S130, initializing and setting parameters of the quantum phase delay of each beluga whale in the preset quantum beluga mechanism, calculating the fitness values ​​of all beluga whales through the constructed fitness function, and obtaining the global quantum phase delay.

[0035] S140, based on the white whale optimization mechanism, iteratively update the quantum phase delay of each white whale to obtain the target quantum phase delay

[0036] S150, mapping the target quantum phase delay to the physical delay, calculating the fitness function values ​​of all beluga whales after the update, and updating the global quantum phase delay.

[0037] S160, when the maximum number of iterations is reached, determine the current global quantum phase delay as the optimal quantum phase delay, and transform the optimal quantum phase delay mapping into the global optimal phase delay. It can be understood that in the process of iterative updating, it is determined whether the current corresponding algebra g reaches the maximum number of iterations G. If not, set g=g+1 and return to S140; if reached, terminate the loop iteration.

[0038] S170, based on the global optimal phase delay, determine the target azimuth weighting coefficient and the target axial weighting coefficient, and perform weighted processing on the azimuth focusing result and the axial focusing result respectively to obtain accurate focusing information.

[0039] The above is a specific implementation method of an ultrasonic focusing method based on the quantum white whale mechanism provided in an embodiment of the present application. The present application calculates the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and establishes an error mapping table in the uv space to cope with different situations where the central beam points in different directions; when looking up the table, linear interpolation is used to reduce the quantization error and weighted processing is performed on the azimuth focusing result and the axial focusing result respectively from the axial direction.

[0040] This application is based on the quantum white whale mechanism. By initializing quantum phase delay, setting parameters, calculating the fitness values ​​of all white whales, iteratively updating the quantum phase delay of each white whale, and mapping phase delay, the optimal quantum phase delay is finally determined and mapped to the global optimal phase delay when the maximum number of iterations is reached. The global optimal phase delay can characterize the estimated values ​​of the azimuth weighting coefficient and the axial weighting coefficient. This application uses the quantum white whale mechanism to quickly obtain the target azimuth weighting coefficient and the target axial weighting coefficient. Compared with the traditional method of only using the largest beam and the second largest beam to look up the table in the azimuth and axial directions, the method of this application fully utilizes the information of the five beams, is not affected by the beam broadening effect, and can solve the designed azimuth and axial weighting coefficients with high precision, and quickly and accurately obtain the estimated results.

[0041] In some embodiments, Figure 2 As shown, the above-mentioned determination of the distribution information of multiple ultrasonic array elements in the ultrasonic phased array, calculation of the azimuth beam width and the axial beam width of the ultrasonic phased array in the uv space, and determination of the beam overlap information to establish an error mapping table, that is, the above-mentioned S110 may specifically include the following steps:

[0042] S210: Determine a horizontal spacing between adjacent ultrasonic array elements along the azimuth direction and an axial spacing between adjacent ultrasonic array elements along the axial direction.

[0043] S220. Based on the horizontal spacing and the axial spacing, calculate the azimuth width and the axial beam width of the ultrasonic phased array.

[0044] S230 , mapping the azimuth beam width and the axial beam width to the uv space to obtain the azimuth beam width and the axial beam width.

[0045] S240: Based on the azimuth beam width and the axial beam width, determine the error functions corresponding to the azimuth direction and the axial direction and establish an error mapping table.

[0046] In the embodiment of the present application, it can be understood that, assuming that the ultrasonic phased array is M 1 ×M 2 dimensional and contains M ultrasonic array elements. The number of ultrasonic array elements distributed along the azimuth direction is M 1 , the number of ultrasonic array elements distributed along the axial direction is M 2 , M=M 1 ×M 2 , then the normal beam width and axial beam width of the ultrasonic phased array are θ 3dB =

[0047] 0.886λM 1 d 1 ,and Where λ is the signal wavelength, d 1 is the spacing between adjacent ultrasonic array elements distributed along the azimuth direction, d 2 is the spacing between adjacent ultrasonic array elements distributed along the axial direction; then the normal beam width θ 3dB . and axial beam width Convert to uv space to get azimuth beam width and axial beam width, that is, u and v dimensional beam width and

[0048] It should be noted that, assuming that the azimuth beam patterns of two adjacent beams with 3dB overlap are and in, represents the axial delay of the two beams, sinc represents the normalized sinc function; assuming that the axial beam patterns of the other two adjacent beams with 3dB overlap are and in, The error functions in the azimuth and axial directions are and According to the error function, a normalized error and an error mapping table are established in the uv space.

[0049] In some embodiments, in the process of acquiring the transient signal data of the received signal and constructing the target optimization function based on the error mapping table, that is, in the aforementioned S120, the maximum number of transient signals can be defined as K, assuming that the beam pointing direction of the central receiving beam of the kth transient signal of the ultrasonic phased array is Convert it to uv space to get the uv beam pointing as Then the uv beam pointing directions of two adjacent beams in azimuth are and The uv beam directions of two adjacent axial beams are and Then the target echo signals of the central receiving beam and the azimuth axial receiving beam are and Among them, x(k)=a(u,v)s(k)+n(k) is the kth transient signal data received by the ultrasonic phased array,

[0051] θ k and represents the phase delay of the kth transient signal target, Ultrasonic phased array Line The phase delay coordinates of the ultrasonic array element where the column is located, s(k) is the received signal, n(k) is an M×1 dimensional complex Gaussian noise vector, k=1, 2, …, K, j is a complex unit, T is a transpose, and H is a conjugate transpose.

[0052] Then, the errors between the central beam and the axially adjacent beams are calculated respectively, and we get and According to the axial error mapping table, the error corresponding to and Then calculate and Then the axial delay of the kth transient signal target is estimated to be Among them, α is the axial weighting coefficient, which calculates the error between the central beam and the adjacent beam in azimuth and

[0053] According to the azimuth error mapping table, the error corresponding and Then calculate and

[0054] Then the azimuth delay of the kth transient signal target is estimated as Among them, β is the azimuth weighting coefficient, and then the target optimization function is constructed after obtaining all transient signal data. The target optimization function satisfies the expression: Here, argmin means finding the variable with the minimum function value. Satisfies the expression: Satisfies the expression:

[0055] In some embodiments, the aforementioned initialization and parameter setting of the quantum phase delay of each beluga whale in the preset quantum beluga mechanism, the fitness values ​​of all beluga whales are calculated by the constructed fitness function to obtain the global quantum phase delay, that is, the aforementioned S130 may specifically include the following steps:

[0056] S310, set the number of white whales to Q, the maximum number of iterations to G, and in the g-th iteration, the quantum phase delay of the q-th white whale in the 2-dimensional search space is When g=1, each dimension of the quantum phase delay of all beluga whales in the first generation is initialized to a uniform random number between [0, 1].

[0057] S320. According to a preset target mapping rule, each dimension of the quantum phase delay of all beluga whales is mapped into the search space to obtain the physical delay of the qth beluga whale of the gth generation.

[0058] S330, calculating the fitness value of each beluga whale according to the constructed fitness function, and determining that the quantum phase delay corresponding to the optimal fitness value is the global quantum phase delay corresponding to the g-th generation

[0059] In one example, the fitness function of the qth beluga whale of the gth generation satisfies the expression:

[0060]

[0061] The target mapping rule satisfies the expression:

[0062]

[0063] In the formula, is the physical delay of the qth beluga whale of the gth generation, is the quantum phase delay of the qth beluga whale of the gth generation, u o，max and u o，min are the upper and lower limits of the o-th dimension search space, q=1, 2,…,Q, o=1, 2; and They correspond to the azimuth weighting coefficient B and the axial weighting coefficient a respectively.

[0064] In some embodiments, the aforementioned iterative update of the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain the target quantum phase delay, that is, the aforementioned S140 may specifically include the following steps:

[0065] S410. During the iterative update process, a balance factor is generated for the qth beluga whale of the gth generation.

[0066] S420. Based on the relationship between the balance factor and the preset critical threshold, determine the o-th dimension target quantum beam phase delay of the q-th beluga whale of the g+1-th generation.

[0067] S430, according to the initial quantum beam phase delay, using a simulated quantum rotating gate to update the quantum phase delay, to obtain a corresponding target quantum phase delay.

[0068] Among them, when the balance factor is less than the preset critical threshold, the o-dimensional target quantum beam phase delay of the q-th beluga whale of the g+1th generation satisfies the expression:

[0069] in, is the o-th dimensional quantum phase delay of the x-th beluga whale of the g-th generation, x is a random integer between [1, Q] except q, and It is not a random number uniformly distributed between [0, 1]. When the balance factor is greater than or equal to the preset critical threshold, the target quantum beam phase delay of the o-th dimension of the q-th beluga whale of the g+1th generation satisfies the expression:

[0070]

[0071] in, For the gth generation The oth-dimensional quantum phase delay of a white whale, is a random integer between [1, Q] except q, and is a random number uniformly distributed between [0, 1]. is the Levy flight coefficient, and is a random number that obeys a Gaussian distribution with a mean of 0 and a variance of 1, μ is a control parameter between [0, 2], and Γ represents the Gamma function.

[0072] In the embodiment of the present application, it can be understood that the global optimal quantum phase delay is updated to the g+1th generation. When the fitness function value of the qth beluga whale is not improved, the o-th dimension quantum beam phase delay of the qth beluga whale of the g+1th generation is determined as in and is a random number uniformly distributed between [0, 1]. For the gth generation The oth-dimensional quantum phase delay of a white whale, is a random integer between [1, Q] except q. Based on the o-dimensional quantum beam phase delay of the q-th beluga whale in the g+1th generation, the target quantum phase delay of the g+1th generation can be obtained by updating the quantum phase delay using a simulated quantum rotating gate; the target quantum phase delay of the g+1th generation satisfies the expression:

[0073] In some embodiments, the balance factor is in, is a random number uniformly distributed between [0, 1]; when the balance factor is less than the preset critical threshold, the target quantum phase delay is Where q = 1, 2, ..., Q, o = 1, 2; when the balance factor is greater than or equal to the preset critical threshold, the target quantum phase delay is in, is the o-th dimension global optimal quantum phase delay until the g-th generation; the target quantum phase delay is

[0075] In the embodiments of the present application, it can be understood that the present application executes a corresponding update strategy based on the relationship between the balance factor and the preset critical threshold to update the quantum beam phase delay to obtain the target quantum beam phase delay; then based on the White Whale optimization mechanism, a simulated quantum rotating gate is used to update the quantum phase delay to obtain the target quantum phase delay; the preset critical threshold can be 0.5.

[0076] In one example, the aforementioned mapping of the target quantum phase delay to the physical delay, calculating the fitness function values ​​of all white whales after the update, and updating the global quantum phase delay, that is, the aforementioned S150 may specifically include the following steps:

[0077] S510, mapping the target quantum phase delay to the physical delay, and calculating the fitness function values ​​of all beluga whales.

[0078] S520. Determine the o-th dimensional quantum beam phase delay of the q-th beluga whale of the g+1-th generation based on the fitness function value of the q-th beluga whale.

[0079] S530. Based on the o-th dimensional quantum beam phase delay of the q-th beluga whale of the g+1-th generation, the quantum phase delay is updated using a simulated quantum rotating gate to obtain the quantum phase delay of the g+1-th generation.

[0080] S540, updating the global quantum phase delay based on the quantum phase delay of the g+1th generation.

[0081] In some embodiments, Figure 3 As shown, the present application provides a phased array radar multi-beam angle measurement system based on the quantum white whale mechanism, which may include the following modules:

[0082] The calculation module 610 is used to determine the distribution information of multiple ultrasonic array elements in the phased array radar, calculate the azimuth beam width and axial beam width of the phased array radar in the uv space, and determine the beam overlap information to establish an error mapping table;

[0083] An optimization function construction module 620, used to obtain transient signal data of a received signal and construct a target optimization function based on an error mapping table;

[0084] The parameter configuration module 630 is used to initialize and set parameters of the quantum phase delay of each beluga whale in the preset quantum beluga mechanism, calculate the fitness value of all beluga whales through the constructed fitness function, and obtain the global quantum phase delay;

[0085] Iteration module 640, used for iteratively updating the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain a target quantum phase delay;

[0086] The first mapping module 650 is used to map the target quantum phase delay to the physical delay, calculate the fitness function values ​​of all white whales after the update, and update the global quantum phase delay;

[0087] A second mapping module 660 is used to determine that the current global quantum phase delay is the optimal quantum phase delay when the maximum number of iterations is reached, and transform the optimal quantum phase delay mapping into the global optimal phase delay;

[0088] The weighting module 670 is used to determine the target azimuth weighting coefficient and the target axial weighting coefficient based on the global optimal phase delay, and perform weighted processing on the azimuth angle measurement result and the axial angle measurement result respectively to obtain accurate angle measurement information.

[0089] According to an embodiment of the present application, any multiple modules among the calculation module 610, the optimization function construction module 620, the parameter configuration module 630, the iteration module 640, the first mapping module 650, the second mapping module 660 and the weighting module 670 can be combined into one module for implementation, or any one of the modules can be split into multiple modules. Alternatively, at least part of the functions of one or more of these modules can be combined with at least part of the functions of other modules and implemented in one block.

Claims

1. An ultrasonic focusing method based on quantum white whale mechanism, characterized in that: include: Determine the spatial distribution information of multiple ultrasonic array elements in the ultrasonic phased array, calculate the azimuthal beam width and axial beam width of the ultrasonic phased array in space, and construct an error mapping table based on the beam overlap. The specific steps are: Determine the distribution intervals of the ultrasonic array elements in the azimuth and axial directions, which are the horizontal interval (azimuth) and the axial interval (axial direction) respectively. According to the two intervals, the azimuth width and axial beam width of the ultrasonic phased array are calculated, and the calculated azimuth width and axial beam width are mapped to the uv space to obtain the azimuth beam width and axial beam width. Based on the target width, the error functions in the azimuth and axial directions are further determined, and an error mapping table is established accordingly. The system obtains the transient signal data of the received signal and constructs the target optimization function based on the error mapping table. The system initializes and sets the parameters of the quantum phase delay of each "white whale" in the preset quantum white whale mechanism, and calculates the fitness values ​​of all "white whales" through the constructed fitness function, thereby determining the initial global optimal quantum phase delay. Based on the White Whale optimization mechanism, the system iteratively updates the quantum phase delay of each "White Whale", gradually approaching the target quantum phase delay. During each iteration, the target quantum phase delay is mapped to the physical delay, and the fitness values ​​of all "White Whales" are recalculated to continuously update the global optimal quantum phase delay. After reaching the predetermined maximum number of iterations, the current global optimal quantum phase delay is confirmed as the target quantum phase delay and mapped to the actual physical delay. This process ensures that the system is continuously improved according to the latest fitness evaluation results until the optimal azimuthal and axial weighting coefficients are found. Convert the target quantum phase delay into physical delay and recalculate the fitness values ​​of all "white whales" to update the global optimal quantum phase delay. When the number of iterations reaches the preset maximum value, the system confirms that the current global optimal quantum phase delay is the final optimal quantum phase delay and maps it to the actual global optimal phase delay. Based on the global optimal phase delay, weighting coefficients in the azimuth and axial directions are determined, and then these weighting coefficients are applied to the focusing results in the azimuth and axial directions, respectively, so as to obtain more accurate focusing information.

2. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 1, characterized in that In the ultrasonic phased array, the plurality of ultrasonic array elements are respectively distributed in the azimuth direction and the axial direction; The determining of distribution information of a plurality of ultrasonic array elements in an ultrasonic phased array, calculating an azimuth beam width and an axial beam width of the ultrasonic phased array in a uv space, and determining beam overlap information to establish an error mapping table includes: Determining the horizontal spacing of adjacent ultrasonic array elements along the azimuth direction and the axial spacing along the axial direction; Based on the horizontal interval and the axial interval, the azimuth width and the axial beam width of the ultrasonic phased array are calculated; Mapping the azimuth width and the axial beam width to the uv space to obtain the azimuth beam width and the axial beam width; Based on the azimuth beamwidth and the axial beamwidth, error functions corresponding to the azimuth direction and the axial direction are determined and an error mapping table is established.

3. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 1, characterized in that: The quantum phase delay of each beluga whale in the preset quantum beluga mechanism is initialized and parameterized, and the fitness values ​​of all beluga whales are calculated by constructing a fitness function to obtain a global quantum phase delay, including: Assume the number of white whales is Q, the maximum number of selections is G, and in the g-th iteration, the quantum phase delay of the q-th white whale in the 2D search space is When g = 1, each dimension of the quantum phase delay of all beluga whales in the first generation is initialized to a uniform random number between [0, 1]; According to the preset target mapping rule, each dimension of the quantum phase delay of all beluga whales is mapped into the search space to obtain the physical delay of the qth beluga whale of the gth generation; The fitness value of each beluga whale is calculated according to the constructed fitness function, and the quantum phase delay corresponding to the optimal fitness value is determined as the global quantum phase delay corresponding to the g-th generation.

4. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 3, characterized in that The target mapping rule satisfies the expression: In the formula, is the physical delay of the qth beluga whale of the gth generation, is the quantum phase delay of the qth beluga whale of the gth generation, u o，max and u o，min are the upper and lower limits of the o-th dimension search space, q=1, 2,…,Q, o=1, 2; and They correspond to the azimuth weighting coefficient β and the axial weighting coefficient α respectively.

5. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 3, characterized in that: The iterative updating of the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain the target quantum phase delay includes: During the iterative update process, for the qth beluga whale in the gth generation, a balance factor is generated; Based on the relationship between the balance factor and the preset critical threshold, determining the o-th dimension target quantum beam phase delay of the q-th beluga whale of the g+1-th generation; According to the initial quantum beam phase delay, the quantum phase delay is updated using a simulated quantum rotating gate to obtain a corresponding target quantum phase delay; Wherein, when the balance factor is less than the preset critical threshold, the target quantum beam phase delay of the oth dimension of the qth beluga whale of the g+1th generation satisfies the expression: in, is the o-th dimensional quantum phase delay of the x-th beluga whale of the g-th generation, x is a random integer between [1, Q] except q, and is a random number uniformly distributed between [0, 1]; When the balance factor is greater than or equal to the preset critical threshold, the target quantum beam phase delay of the oth dimension of the qth beluga whale of the g+1th generation satisfies the expression: in, is the o-th dimensional quantum phase delay of the g-th generation beluga whale, is a random integer between [1, Q] except q, and is a random number uniformly distributed between [0, 1]. is the Levy flight coefficient, and is a random number that obeys a Gaussian distribution with a mean of 0 and a variance of 1, μ is a control parameter between [0, 2], and Γ represents the Gamma function.

6. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 5, characterized in that The balance factor is in is a random number uniformly distributed between [0, 1]; When the balance factor is less than the preset critical threshold, the target quantum phase delay is Where, q = 1, 2, ..., Q, o = 1, 2; When the balance factor is greater than or equal to the preset critical threshold, the target quantum phase delay is in, is the global optimal quantum phase delay of the gth dimension until the gth generation; the target quantum phase delay is 7. The ultrasonic focusing method based on the quantum white whale mechanism as claimed in claim 1, characterized in that The mapping of the target quantum phase delay to a physical delay, calculating the fitness function values ​​of all white whales after the update, and updating the global quantum phase delay include: Mapping the target quantum phase delay to a physical delay, and calculating the fitness function values ​​of all beluga whales; Based on the fitness function value of the qth beluga whale, determine the oth-dimensional quantum beam phase delay of the qth beluga whale of the g+1th generation; Based on the o-th dimensional quantum beam phase delay of the q-th beluga whale of the g+1-th generation, the quantum phase delay is updated using a simulated quantum rotating gate to obtain the quantum phase delay of the g+1-th generation; The global quantum phase delay is updated based on the quantum phase delay of the g+1th generation.

8. An ultrasonic focusing system based on the quantum white whale mechanism, characterized in that: include: A calculation module, used to determine the distribution information of multiple ultrasonic array elements in the ultrasonic phased array, calculate the azimuth beam width and axial beam width of the ultrasonic phased array in the uv space, and determine the beam overlap information to establish an error mapping table; An optimization function construction module, used for acquiring transient signal data of a received signal and constructing a target optimization function based on the error mapping table; The parameter configuration module is used to initialize and set the parameters of the quantum phase delay of each beluga whale in the preset quantum beluga mechanism, calculate the fitness value of all beluga whales through the constructed fitness function, and obtain the global quantum phase delay; An iteration module is used to iteratively update the quantum phase delay of each beluga whale based on the beluga optimization mechanism to obtain the target quantum phase delay; A first mapping module is used to map the target quantum phase delay to a physical delay, calculate the fitness function values ​​of all white whales after the update, and update the global quantum phase delay; A second mapping module is used to determine that the current global quantum phase delay is the optimal quantum phase delay when the maximum number of iterations is reached, and transform the optimal quantum phase delay mapping into a global optimal phase delay; The weighting module is used to determine the target azimuth weighting coefficient and the target axial weighting coefficient based on the global optimal phase delay, and perform weighted processing on the azimuth focusing result and the axial focusing result respectively to obtain accurate focusing information.

Images