Dipole array jammer modeling method, device and system

Abstract

The application provides a dipole array interference source modeling method, device and system, and belongs to the field of interference source modeling, which comprises the following steps: detecting the magnetic field amplitude of a printed circuit board; obtaining initial vector parameters and an initial vector parameter set; performing linear transformation and integration on the initial vector parameter set by using optimal initial vector parameters in the initial vector parameter set to obtain integrated vector parameters and an experiment set; determining the initial difference between the detected magnetic field amplitude and the dipole array magnetic field amplitude corresponding to the initial vector parameter set, and determining the experimental difference between the detected magnetic field amplitude and the dipole array magnetic field amplitude corresponding to the experiment set; determining the current optimal vector parameters according to the initial difference and the experimental difference; obtaining the vector parameters of the optimal dipole array; and establishing a printed circuit board interference source model. According to the method of the application, the radiation source model can be easily simulated, the simulation time is shortened, and the system occupation is reduced.

Description

Technical Field

[0001] This invention relates to the field of interference source modeling technology, specifically to a dipole array interference source modeling method, a dipole array interference source modeling device, and a dipole array interference source modeling system. Background Technology

[0002] Chips play an indispensable role in modern electronic products. As operating frequencies continue to expand and chip integration and complexity increase, the electromagnetic interference radiated by chips to nearby circuits cannot be ignored. In this context, anticipating such interference in advance is extremely useful.

[0003] In existing technologies, integrated circuits are mostly quite complex, making it very difficult to model interference sources; moreover, some interference source models are quite complex, requiring a lot of time and consuming a large amount of system resources during simulation analysis. Summary of the Invention

[0004] To address the technical problem that modeling interference sources is very difficult in existing technologies, this invention provides a method, device, and system for modeling interference sources using a dipole array. This method can quickly and easily establish an equivalent model of a radiation source that is easy to simulate and analyze, thereby shortening simulation time and reducing system usage.

[0005] To achieve the above objectives, the first aspect of the present invention provides a method for modeling interference sources of a dipole array, comprising the following steps: S101: detecting the magnetic field amplitude within a set range on a printed circuit board; S102: setting the number of dipoles in the dipole array to n; S103: initializing the vector parameters of the dipole array to obtain initial vector parameters and an initial vector parameter set; wherein the initial vector parameter set includes initial vector parameters; S104: using the optimal initial vector parameters in the initial vector parameter set to perform a linear transformation and integration on the initial vector parameter set to obtain integrated vector parameters and an experimental set; wherein the experimental set includes integrated vector parameters; S105: determining the initial difference between the magnetic field amplitude of the dipole array corresponding to the initial vector parameter set and the magnetic field amplitude within the set range on the printed circuit board, and determining the experimental difference between the magnetic field amplitude of the dipole array corresponding to the experimental set and the magnetic field amplitude within the set range on the printed circuit board; S106: based on the relationship between the initial difference and the experimental difference, from the integrated vector parameters... S107: Using the current optimal vector parameters and the current set of vector parameters as initial vector parameters, and using the current set of vector parameters as the initial set of vector parameters, repeat S104-S106 a set number of times; S108: Increase the number of dipoles in the dipole array by s to form a new dipole array, repeat S103-S107 until the number of dipoles in the new dipole array reaches a set number N, determine the current new dipole array as the optimal dipole array, and determine the vector parameters of the optimal dipole array; wherein, the initial vector parameters and the initial set of vector parameters in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106; wherein, s≥1, N>n+s; S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0006] Furthermore, the initial vector parameters and the set of initial vector parameters are obtained in the following way:

[0007] d i,j,0 =d j,L +r j (d j,U -d j,L );

[0008]

[0009] Where, d i,j,G Let i be the initial vector parameters, i = 1, 2, 3, ..., N P j = 1, 2, ..., D, where D is the number of parameters in the dipole vector parameters; D Gd is the initial vector parameter set; j,U d represents the maximum value of the j-th parameter in the dipole vector parameters; j,L Let r be the minimum value of the j-th parameter in the dipole vector parameters; j is a random number uniformly distributed from [0,1].

[0010] Furthermore, the step of using the optimal initial vector parameter in the initial vector parameter set to perform a linear transformation and integration of the initial vector parameter set to obtain integrated vector parameters and an experimental set includes: using the optimal initial vector parameter in the initial vector parameter set and two randomly selected different initial vector parameters to perform a linear transformation of the initial vector parameter set to obtain linearly changed vector parameters; and integrating some initial vector parameters and some linearly changed vector parameters to obtain integrated vector parameters and an experimental set.

[0011] Furthermore, the parameters of the linear transformation vector are determined in the following way:

[0012] v i,G =d i,G +F1(d best,G -d i,G )+F2(d r1,G -d r2,G );

[0013] Among them, v i,G d is a linearly changing vector parameter; i,G d is the initial vector parameter; best,G d represents the optimal initial vector parameters. r1,G and d r2,G Two different initial vector parameters are randomly selected; r1 and r2 are from [1, N] P The mutually exclusive integers randomly selected from [i] are not equal to i; F1 and F2 are mutation factors.

[0014] Furthermore, the integrated vector parameters are obtained in the following way:

[0015] ;

[0016] Where ui,j,G are the integrated vector parameters, j=1,2,...,k; r j C is a random number uniformly distributed from [0,1]. R ∈[0,1] is the crossover factor; r(i) is a random integer uniformly distributed from [0,k].

[0017] Furthermore, the experimental set was obtained in the following way:

[0018]

[0019] Among them, DG,i This is the experimental set.

[0020] Furthermore, the initial difference and experimental difference were determined in the following ways:

[0021]

[0022]

[0023] Where, f(D) G f(D) represents the initial difference; G,i The difference is the experimental value. and The magnitude of the dipole array magnetic field corresponding to the initial vector parameter set; and To detect the amplitude of the magnetic field; and M represents the amplitude of the dipole array magnetic field corresponding to the experimental set; M is the number of near-field scanning points when detecting the magnetic field amplitude within a set range on the printed circuit board.

[0024] Furthermore, determining the current optimal vector parameter and the current set of vector parameters from the integrated vector parameters and the initial vector parameters based on the relationship between the initial difference and the experimental difference includes: determining the integrated vector parameter as the current optimal vector parameter when the initial difference is greater than or equal to the experimental difference; and determining the initial vector parameter as the current optimal vector parameter when the initial difference is less than the experimental difference.

[0025] The current set of vector parameters is obtained based on the current optimal vector parameters:

[0026] ;

[0027] Where DG+1 is the current vector parameter set.

[0028] A second aspect of the present invention provides a modeling device for interference sources of a dipole array, the device comprising: a detection unit for performing step S101: detecting the magnetic field amplitude within a set range on a printed circuit board;

[0029] The modeling unit performs the following steps: S102: Set the number of dipoles in the dipole array to n; S103: Initialize the vector parameters of the dipole array to obtain initial vector parameters and an initial vector parameter set; wherein the initial vector parameter set includes the initial vector parameters; S104: Perform linear transformation and integration on the initial vector parameter set using the optimal initial vector parameters in the initial vector parameter set to obtain integrated vector parameters and an experimental set; wherein the experimental set includes the integrated vector parameters; S105: Determine the initial difference between the magnetic field amplitude of the dipole array corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the magnetic field amplitude of the dipole array corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board; S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter from the integrated vector parameters and the initial vector parameters. Parameter set; where the current vector parameter set includes the current optimal vector parameter; S107: Use the current optimal vector parameter as the initial vector parameter, and use the current vector parameter set as the initial vector parameter set, repeat S104-S106 a set number of times; S108: Increase the number of dipoles in the dipole array by s to form a new dipole array, repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N, determine the current new dipole array as the optimal dipole array, and determine the vector parameters of the optimal dipole array; where the initial vector parameter and initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameter determined in the previous execution of S106; where s≥1, N>n+s; S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0030] A third aspect of the present invention provides a dipole array interference source modeling system, including the dipole array interference source modeling device described above.

[0031] The present invention has at least the following technical effects through the technical solution provided by the present invention:

[0032] The dipole array interference source modeling method of the present invention first detects the magnetic field amplitude within a set range on the printed circuit board, and then performs interference source modeling: First, a dipole array is set up, and the vector parameters of the dipoles are initialized to obtain initial vector parameters and an initial vector parameter set. Then, the optimal initial vector parameters in the initial vector parameter set are used to perform a linear transformation and integration on the initial vector parameter set to obtain integrated vector parameters and an experimental set. The initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range on the printed circuit board is determined, as well as the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range on the printed circuit board is determined. Based on the initial difference... The relationship between the experimental and experimental differences is determined from the integrated vector parameters and the initial vector parameters to identify the current optimal vector parameters and the current set of vector parameters. The current optimal vector parameters are used as the initial vector parameters, and the current set of vector parameters is used as the initial set of vector parameters. The linear transformation, integration, and difference calculation are restarted, and the iteration is repeated. After repeating the process a set number of times, the number of dipoles in the dipole array is increased by s to form a new dipole array. The newly formed dipole array is then initialized, linearly transformed, integrated, and the difference calculation is repeated until the number of dipoles in the new dipole array reaches the set number N. The current new dipole array is then identified as the optimal dipole array, and the vector parameters of the optimal dipole array are determined.

[0033] The dipole array interference source modeling method provided by the present invention can easily and quickly establish an equivalent radiation source model that is easy to simulate and analyze, shortening the simulation time and reducing system usage.

[0034] Other features and advantages of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0035] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:

[0036] Figure 1 A flowchart of a dipole array interference source modeling method provided in an embodiment of the present invention;

[0037] Figure 2 The equivalent radiation source model diagram of the equivalent dipole array in the dipole array interference source modeling method provided in the embodiments of the present invention is shown.

[0038] Figure 3 This is a schematic diagram of the printed circuit board of the device under test in the dipole array interference source modeling method provided in the embodiments of the present invention;

[0039] Figure 4This is a graph showing the variation of the fitness of the best individual in each generation in the dipole array interference source modeling method provided in this embodiment of the invention.

[0040] Figure 5 The graph shows the change curves of the magnetic field amplitude within a set range on the printed circuit board in the dipole array interference source modeling method provided in the embodiments of the present invention and the magnetic field strength value calculated by the present invention.

[0041] Figure 6 This is a schematic diagram of the apparatus for modeling a dipole array interference source provided in an embodiment of the present invention. Detailed Implementation

[0042] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the scope of the present invention.

[0043] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other.

[0044] In this invention, unless otherwise stated, directional terms such as "upper," "lower," "top," and "bottom" are generally used to describe the relative positions of components in relation to the directions shown in the accompanying drawings or in relation to the vertical, perpendicular, or gravitational directions.

[0045] The present invention will now be described in detail with reference to the accompanying drawings and embodiments.

[0046] Please refer to Figure 1This invention provides a method for modeling interference sources in a dipole array, comprising the following steps: S101: detecting the magnetic field amplitude within a set range on a printed circuit board; S102: setting the number of dipoles in the dipole array to n; S103: initializing the vector parameters of the dipole array to obtain initial vector parameters and an initial vector parameter set; wherein the initial vector parameter set includes initial vector parameters; S104: performing a linear transformation and integration on the initial vector parameter set using the optimal initial vector parameters in the initial vector parameter set to obtain integrated vector parameters and an experimental set; wherein the experimental set includes integrated vector parameters; S105: determining the initial difference between the magnetic field amplitude of the dipole array corresponding to the initial vector parameter set and the magnetic field amplitude within the set range on the printed circuit board, and determining the experimental difference between the magnetic field amplitude of the dipole array corresponding to the experimental set and the magnetic field amplitude within the set range on the printed circuit board; S106: based on the relationship between the initial difference and the experimental difference, from the integrated vector parameters and... The initial vector parameters determine the current optimal vector parameters and the current vector parameter set; the current vector parameter set includes the current optimal vector parameters; S107: Using the current optimal vector parameters as the initial vector parameters and the current vector parameter set as the initial vector parameter set, repeat S104-S106 a set number of times; S108: Increase the number of dipoles in the dipole array by s to form a new dipole array, repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N, determine the current new dipole array as the optimal dipole array, and determine the vector parameters of the optimal dipole array; wherein, the initial vector parameters and the initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106; wherein, s≥1, N>n+s; S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0047] Specifically, in this embodiment of the invention, the magnetic field amplitude within a set range on the printed circuit board is first detected, then the interference source is modeled, and the dipole array is set in the equivalent dipole plane. It is assumed that a probe is used to detect the equivalent dipole array (e.g.,...). Figure 2 As shown, the magnetic field amplitude of the detection probe and the dipole array in the equivalent dipole plane is calculated. The difference between the magnetic field amplitude and the dipole array amplitude within the set range of the printed circuit board is calculated. The dipole vector parameters of the dipole array with the smallest difference are selected. The vector parameters (number, position and torque amplitude) of the optimal dipole array are determined. Then, the interference source model of the printed circuit board is established based on the vector parameters of the optimal dipole array.

[0048] The interference source modeling process is as follows:

[0049] First, the number of dipoles in the dipole array is set to n. In this application, the parameters of each dipole include the x and y coordinates of its position, as well as the real and imaginary parts of Mx and My, forming a vector a = (x, y, Re[M... x ],Re[M y ],Im[M x ],Im[M y The parameters of all dipoles in the dipole array form a vector parameter d. i =(a1,…,a n Next, the vector parameters of the dipole are initialized to obtain initial vector parameters and an initial vector parameter set. Then, the optimal initial vector parameters in the initial vector parameter set are used to perform a linear transformation and integration on the initial vector parameter set to obtain integrated vector parameters and an experimental set. The initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board is determined, as well as the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board is determined. The vector parameters (number, position, and torque amplitude) of the dipole are filtered by the difference between the magnetic field amplitude within the set range of the printed circuit board and the dipole array amplitude, thus obtaining the current optimal vector parameters and the current vector parameter set. The current optimal vector parameters are used as the initial vector parameters, and the current set of vector parameters is used as the initial vector parameter set. This process is repeated iteratively a set number of times. Then, the number of dipoles in the dipole array is increased by `s`, forming a new dipole array. Initialization, linear transformation, integration, and difference calculation are then performed again until the number of dipoles in the new dipole array reaches the set number `N`. This new dipole array is then determined as the optimal dipole array, along with its vector parameters. During the iterative process, the optimal vector parameters obtained before adding dipoles are used to initialize the vector parameters of the dipole array, obtaining the initial vector parameters and the initial vector parameter set. Then, a printed circuit board interference source model is established based on the vector parameters of the optimal dipole array. Finally, the dipole array is used to replace the unknown or complex radiation source model, facilitating subsequent simulation analysis.

[0050] The method provided by this invention enables the simple and rapid establishment of an equivalent model of a radiation source that is easy to simulate and analyze, thereby shortening simulation time and reducing system usage.

[0051] The specific steps are as follows:

[0052] First, execute S101: Detect the magnetic field amplitude within the set range of the printed circuit board.

[0053] Next, execute S102: Set the number of dipoles in the dipole array to n.

[0054] Furthermore, the initial vector parameters and the set of initial vector parameters are obtained in the following way:

[0055] d i,j,0 =d j,L +r j (d j,U -d j,L );

[0056]

[0057] Where, d i,j,G Let i be the initial vector parameters, i = 1, 2, 3, ..., N P j = 1, 2, ..., D, where D is the number of parameters in the dipole vector parameters, and N P For control parameters; D G d is the initial vector parameter set; j,U d represents the maximum value of the j-th parameter in the dipole vector parameters; j,L Let r be the minimum value of the j-th parameter in the dipole vector parameters; j G is a random number uniformly distributed from [0,1]; G is the algebra for the selection.

[0058] Next, execute S103: initialize the vector parameters of the dipole array to obtain the initial vector parameters and the initial vector parameter set; wherein, the initial vector parameter set includes the initial vector parameters.

[0059] Furthermore, the step of using the optimal initial vector parameters from the initial vector parameter set to perform a linear transformation and integration on the initial vector parameter set to obtain integrated vector parameters and an experimental set includes:

[0060] The initial vector parameter set is linearly transformed using the optimal initial vector parameter from the initial vector parameter set and two randomly selected different initial vector parameters to obtain linearly transformed vector parameters; some initial vector parameters and some linearly transformed vector parameters are integrated to obtain integrated vector parameters and experimental set.

[0061] Furthermore, the parameters of the linear transformation vector are determined in the following way:

[0062] v i,G =d i,G +F1(d best,G -d i,G )+F2(d r1,G -d r2,G );

[0063] Among them, v i,G d is a linearly changing vector parameter; i,G d is the initial vector parameter; best,Gd represents the optimal initial vector parameters. r1,G and d r2,G Two different initial vector parameters are randomly selected; r1 and r2 are from [1, N] P The mutually exclusive integers randomly selected from [i] are not equal to i; F1 and F2 are mutation factors.

[0064] F1(d best,G -d i,G It can guarantee the vector parameter v i,G Including the optimal partial genes (vector parameters) from generation G can improve the convergence of the algorithm.

[0065] Next, execute S104: Use the optimal initial vector parameters in the initial vector parameter set to perform a linear transformation and integration on the initial vector parameter set to obtain the integrated vector parameters and the experimental set, wherein the experimental set includes the integrated vector parameters.

[0066] Furthermore, the integrated vector parameters are obtained in the following way:

[0067] ;

[0068] Where ui,j,G are the integrated vector parameters, j=1,2,...,k; r j C is a random number uniformly distributed from [0,1]. R ∈[0,1] is the crossover factor; r(i) is a random integer uniformly distributed from [0,k].

[0069] Furthermore, the experimental set was obtained in the following way:

[0070]

[0071] Among them, D G,i This is the experimental set.

[0072] Next, execute S105: determine the initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board.

[0073] Furthermore, the initial difference and experimental difference were determined in the following ways:

[0074]

[0075]

[0076] Where, f(D) G f(D) represents the initial difference; G,i The difference is the experimental value. and The magnitude of the dipole array magnetic field corresponding to the initial vector parameter set; and To detect the amplitude of the magnetic field; and M represents the amplitude of the dipole array magnetic field corresponding to the experimental set; M is the number of near-field scanning points when detecting the magnetic field amplitude within a set range on the printed circuit board.

[0077]

[0078]

[0079]

[0080]

[0081]

[0082] Where j is the imaginary unit, e is the natural constant; k represents the free space wavenumber. f is the frequency; c is the speed of light in a vacuum; M x and M y This represents the horizontal magnetic dipole moment, with units of A·m. 2 R represents the distance between the near-field sampled data point and each dipole in the dipole array; H x and H y This is the magnetic field amplitude that has already taken into account the influence of the ground plane.

[0083] Next, execute S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter set from the integrated vector parameter and the initial vector parameter; wherein, the current vector parameter set includes the current optimal vector parameter.

[0084] Furthermore, determining the current optimal vector parameter and the current set of vector parameters from the integrated vector parameters and the initial vector parameters based on the relationship between the initial difference and the experimental difference includes: determining the integrated vector parameter as the current optimal vector parameter when the initial difference is greater than or equal to the experimental difference; and determining the initial vector parameter as the current optimal vector parameter when the initial difference is less than the experimental difference.

[0085] The current set of vector parameters is obtained based on the current optimal vector parameters:

[0086] ;

[0087] Where DG+1 is the current vector parameter set.

[0088] Next, execute S107: using the current optimal vector parameters as the initial vector parameters, and the current set of vector parameters as the initial set of vector parameters, repeat S104-S106 a predetermined number of times to obtain the optimal vector parameters for the n dipole arrays, d i =(a1,…,a n ).

[0089] In this embodiment, the vector parameters d of the n dipole arrays can also be obtained when the maximum number of repetitions or the fitness value reaches a set value. i =(a1,…,a n Please refer to this. Figure 4 The more iterations there are, the lower the fitness becomes, but the higher the accuracy of the model.

[0090] Next, execute S108: Increase the number of dipoles in the dipole array by s to form a new dipole array. Repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N. Determine the current new dipole array as the optimal dipole array and determine the vector parameters of the optimal dipole array. In each execution of S103, the initial vector parameters and the initial vector parameter set are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106. Where s≥1, N>n+s.

[0091] Next, execute S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0092] Specifically, in this embodiment of the invention, after obtaining the vector parameters of the optimal dipole array, the vector parameters of the optimal dipole array are imported into the full-wave simulation software HFSS for analysis via a script, which can more accurately predict the radiation emission of the original device under test.

[0093] Example 1

[0094] Please refer to Figure 3 The device under test (DUT) is a 54mm*72mm*1mm printed circuit board with a chip soldered on it. One side has a 2mm wide L-shaped microstrip line, and the other side has a ground plane. The dielectric constant and loss tangent of the intermediate dielectric layer are 4.7 and 0.025, respectively. One end of the microstrip line is a 1V voltage source, and the other end is a 50Ω source. The operating frequency is 1GHz. Step S101 is performed: In a microwave anechoic chamber, a scanner is used to perform a near-field scan at a point 10mm above the microstrip line. The scanning plane is 30mm*50mm, and the scanning interval in the x and y directions is 2mm. A total of 416 H points are acquired. x H y The actual magnetic field amplitude of the near field data is saved to the computer as reference data.

[0095] Then, interference source modeling is performed:

[0096] S102: Set the number of dipoles in the dipole array to n.

[0097] First, assume the number of dipoles n = 1, and the dipole array vector parameters d1 = a1 = (x1, y1, Re[M1], Re[M1], Im[M1], Im[M1]). In this application, the control parameter N... P C R F1 and F2 are selected as 30, 0.5, 0.1 and 0.8 respectively. Those skilled in the art can set them according to the actual situation.

[0098] In this embodiment, when initially setting the number of dipoles in the dipole array, setting it from one dipole array allows for referencing the results calculated from one dipole array, thus reducing the search space for other numbers (n≥2) of dipole arrays and accelerating the algorithm's running speed.

[0099] Those skilled in the art may also set up an array of other numbers of dipoles instead of starting with a single dipole array; this embodiment does not impose any restrictions.

[0100] S103: Initialize the vector parameters of the dipole array to obtain initial vector parameters and an initial vector parameter set; wherein, the initial vector parameter set includes the initial vector parameters. A trial-and-error initial guess is used for the search space of a1.

[0101] Furthermore, the initial vector parameters and the set of initial vector parameters are obtained in the following way:

[0102] d i,j,G =a 1,j,L +r j (a 1,j,U -a 1,j,L );

[0103]

[0104] Where, d i,j,G Let i be the initial vector parameters, i = 1, 2, 3, ..., 30, j = 1, 2, ..., 6, and D be the number of parameters in the dipole vector parameters; D G The initial vector parameter set; a 1,j,L Let a be the maximum value of the j-th parameter in the dipole vector parameter a1. 1,j,U Let a be the minimum value of the j-th parameter in the dipole vector parameter a1, and the first a 1,j,L and a 1,j,U Using trial and error as an initial guess; r jG is a random number uniformly distributed from [0,1]; G is the algebra of the selection (i.e., the number of repetitions).

[0105] S104: Use the optimal initial vector parameters in the initial vector parameter set to perform a linear transformation and integration on the initial vector parameter set to obtain the integrated vector parameters and the experimental set; wherein, the experimental set includes the integrated vector parameters.

[0106] Furthermore, the step of using the optimal initial vector parameter in the initial vector parameter set to perform a linear transformation and integration of the initial vector parameter set to obtain integrated vector parameters includes: using the optimal initial vector parameter in the initial vector parameter set and two randomly selected different initial vector parameters to perform a linear transformation of the initial vector parameter set to obtain linearly changed vector parameters; and integrating some initial vector parameters and some linearly changed vector parameters to obtain integrated vector parameters and an experimental set.

[0107] Furthermore, the parameters of the linear transformation vector are determined in the following manner:

[0108] v i,G =d i,G +F1(d best,G -d i,G )+F2(d r1,G -d r2,G );

[0109] Among them, v i,G d is a linearly changing vector parameter; i,G d is the initial vector parameter; best,G d represents the optimal initial vector parameters. r1,G and d r2,G Two different initial vector parameters are randomly selected; r1 and r2 are from [1, N] P The integers randomly selected from the set are not equal to i; F1 and F2 are mutation factors, F1 = 0.1 and F2 = 0.8.

[0110] F1(d best,G -d i,G It can guarantee the vector parameter v i,G It contains the optimal portion of genes (vector parameters) from generation G.

[0111] Furthermore, the integrated vector parameters are obtained in the following way:

[0112] ;

[0113] Where ui,j,G are the integrated vector parameters, j=1,2,...,k; r j C is a random number uniformly distributed from [0,1]. R∈[0,1] is the crossover factor; r(i) is a random integer uniformly distributed from [0,k].

[0114] Furthermore, the experimental set was obtained in the following way:

[0115]

[0116] Among them, D G,i This is the experimental set.

[0117] S105: Determine the initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board.

[0118] Assuming the number of near-field scan points is M, and the number of equivalent magnetic dipoles is n, where each equivalent magnetic dipole can be represented as a six-dimensional vector, then the fitness function to be optimized is f(D). G It can be defined as the difference between the amplitude of the measured magnetic field and the amplitude of the magnetic field generated by the equivalent magnetic dipole.

[0119] Furthermore, the initial difference and experimental difference were determined in the following ways:

[0120]

[0121]

[0122] Where, f(D) G f(D) represents the initial difference; G,i The difference is the experimental value. and The magnitude of the dipole array magnetic field corresponding to the initial vector parameter set; and To detect the amplitude of the magnetic field; and M represents the amplitude of the dipole array magnetic field corresponding to the experimental set; M is the number of near-field scanning points when detecting the magnetic field amplitude within a set range on the printed circuit board.

[0123]

[0124]

[0125]

[0126]

[0127]

[0128] Where j is the imaginary unit, e is the natural constant; k represents the free space wavenumber. f is the frequency; c is the speed of light in a vacuum; M x and M y This represents the horizontal magnetic dipole moment, with units of A·m. 2 R represents the distance between the near-field sampled data point and each dipole in the dipole array; H x and H y This is the magnetic field amplitude that has already taken into account the influence of the ground plane.

[0129] S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter set from the integrated vector parameter and the initial vector parameter; wherein, the current vector parameter set includes the current optimal vector parameter.

[0130] Compare f(D) G ) and f(D G,i The fitness of ) determines the retention of d. i,G or u i,G .

[0131] Furthermore, determining the current optimal vector parameter and the current set of vector parameters from the integrated vector parameters and the initial vector parameters based on the relationship between the initial difference and the experimental difference includes: determining the integrated vector parameter as the current optimal vector parameter when the initial difference is greater than or equal to the experimental difference; and determining the initial vector parameter as the current optimal vector parameter when the initial difference is less than the experimental difference.

[0132] The current set of vector parameters is obtained based on the current optimal vector parameters:

[0133] ;

[0134] Among them, D G+1 This is the current set of vector parameters.

[0135] S107: Using the current optimal vector parameters as the initial vector parameters and the current set of vector parameters as the initial set of vector parameters, repeat S104-S106 a set number of times, or until the maximum number of repetitions is reached or the fitness value reaches the set value, to obtain the optimal vector parameters d1=a1 of a dipole array.

[0136] S108: Increase the number of dipoles in the dipole array in S107 by 1 to form a new dipole array. Repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N. Determine the dipole array with N dipoles as the optimal dipole array, and determine the vector parameter d of the optimal dipole array. i=(a1,…,a N The initial vector parameters and initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106.

[0137] In this embodiment, the optimal dipole array consists of 6-8 dipoles. Those skilled in the art can set the array according to the actual situation, and this embodiment does not impose any restrictions.

[0138] S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0139] Specifically, in this embodiment of the invention, after obtaining the vector parameters of the optimal dipole array, the vector parameters of the optimal dipole array are imported into the full-wave simulation software HFSS for analysis via a script, which can predict the radiated emission of the original device under test relatively accurately. Please refer to... Figure 5 , Figure 5 The dashed line represents the magnetic field amplitude calculated by the simulation software, while the solid line represents the magnetic field amplitude actually sampled by the printed circuit board. The magnetic field amplitude simulated by this invention is quite close to the actual sampled value.

[0140] Please refer to Figure 6 A second aspect of the present invention provides a dipole array interference source modeling apparatus, the apparatus comprising:

[0141] The detection unit is used to perform step S101: detect the magnetic field amplitude within a set range on the printed circuit board;

[0142] Modeling unit, used to perform the following steps:

[0143] S102: Set the number of dipoles in the dipole array to n;

[0144] S103: Initialize the vector parameters of the dipole array to obtain the initial vector parameters and the initial vector parameter set; wherein, the initial vector parameter set includes the initial vector parameters.

[0145] S104: Using the optimal initial vector parameters in the initial vector parameter set, perform a linear transformation and integration on the initial vector parameter set to obtain the integrated vector parameters and the experimental set; wherein, the experimental set includes the integrated vector parameters;

[0146] S105: Determine the initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board.

[0147] S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter set from the integrated vector parameter and the initial vector parameter; wherein, the current vector parameter set includes the current optimal vector parameter;

[0148] S107: Use the current optimal vector parameters as the initial vector parameters, and use the current set of vector parameters as the initial set of vector parameters, and repeat the execution of S104-S106 a set number of times.

[0149] S108: Increase the number of dipoles in the dipole array by s to form a new dipole array. Repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N. Determine the current new dipole array as the optimal dipole array and determine the vector parameters of the optimal dipole array. The initial vector parameters and the initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106. Where s≥1, N>n+s.

[0150] S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

[0151] A third aspect of the present invention provides a dipole array interference source modeling system, including the dipole array interference source modeling device described above.

[0152] The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings. However, the present invention is not limited to the specific details of the above embodiments. Within the scope of the technical concept of the present invention, various simple modifications can be made to the technical solution of the present invention, and these simple modifications all fall within the protection scope of the present invention.

[0153] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present invention will not describe the various possible combinations separately.

[0154] Furthermore, various different embodiments of the present invention can be combined in any way, as long as they do not violate the spirit of the present invention, they should also be regarded as the content disclosed by the present invention.

Claims

1. A method for modeling interference sources in a dipole array, characterized in that, The dipole array interference source modeling method includes: S101: Detects the magnetic field amplitude within a set range on the printed circuit board; S102: Set the number of dipoles in the dipole array to n; S103: Initialize the vector parameters of the dipole array to obtain the initial vector parameters and the initial vector parameter set; wherein, the initial vector parameter set includes the initial vector parameters. S104: Using the optimal initial vector parameters in the initial vector parameter set, perform a linear transformation and integration on the initial vector parameter set to obtain the integrated vector parameters and the experimental set; wherein, the experimental set includes the integrated vector parameters; S105: Determine the initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board. S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter set from the integrated vector parameter and the initial vector parameter; wherein, the current vector parameter set includes the current optimal vector parameter; S107: Use the current optimal vector parameters as the initial vector parameters, and use the current set of vector parameters as the initial set of vector parameters, and repeat the execution of S104-S106 a set number of times. S108: Increase the number of dipoles in the dipole array by s to form a new dipole array. Repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N. Determine the current new dipole array as the optimal dipole array and determine the vector parameters of the optimal dipole array. The initial vector parameters and the initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106. Where s≥1, N>n+s. S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

2. The dipole array interference source modeling method according to claim 1, characterized in that, The initial vector parameters and the set of initial vector parameters are obtained in the following way: d i,j,0 =d j,L +r j (d j,U -d j,L ); Where, d i,j,G Let i be the initial vector parameters, i = 1, 2, 3, ..., N P j = 1, 2, ..., D, where D is the number of parameters in the dipole vector parameters; D G d is the initial vector parameter set; j,U d represents the maximum value of the j-th parameter in the dipole vector parameters; j,L Let r be the minimum value of the j-th parameter in the dipole vector parameters; j is a random number uniformly distributed from [0,1].

3. The dipole array interference source modeling method according to claim 2, characterized in that, The process of using the optimal initial vector parameters from the initial vector parameter set to perform a linear transformation and integration of the initial vector parameter set to obtain integrated vector parameters and an experimental set includes: By using the optimal initial vector parameter in the initial vector parameter set and two randomly selected different initial vector parameters, a linear transformation is performed on the initial vector parameter set to obtain the linearly transformed vector parameters. By integrating some initial vector parameters and some linearly transformed vector parameters, we obtain integrated vector parameters and an experimental set.

4. The dipole array interference source modeling method according to claim 3, characterized in that, The parameters of the linear transformation vector are determined in the following way: v i,G =d i,G +F1(d best,G -d i,G )+F2(d r1,G -d r2,G ); Among them, v i,G d is a linearly changing vector parameter; i,G d is the initial vector parameter; best,G d represents the optimal initial vector parameters. r1,G and d r2,G Two different initial vector parameters are randomly selected; r1 and r2 are from [1, N] P The mutually exclusive integers randomly selected from [i] are not equal to i; F1 and F2 are mutation factors.

5. The dipole array interference source modeling method according to claim 4, characterized in that, The integrated vector parameters are obtained in the following way: ； Among them, u i,j,G For the integrated vector parameters, j = 1, 2, ..., k; r j C is a random number uniformly distributed from [0,1]. R ∈[0,1] is the crossover factor; r(i) is a random integer uniformly distributed from [0,k].

6. The dipole array interference source modeling method according to claim 5, characterized in that, The experimental set was obtained in the following way: Among them, D G,i This is the experimental set.

7. The dipole array interference source modeling method according to claim 6, characterized in that, The initial difference and experimental difference are determined in the following ways: Where, f(D) G f(D) represents the initial difference; G,i The difference is the experimental value. and The magnitude of the dipole array magnetic field corresponding to the initial vector parameter set; and To detect the amplitude of the magnetic field; and M represents the amplitude of the dipole array magnetic field corresponding to the experimental set; M is the number of near-field scanning points when detecting the magnetic field amplitude within a set range on the printed circuit board.

8. The dipole array interference source modeling method according to claim 7, characterized in that, The step of determining the current optimal vector parameters and the current set of vector parameters from the integrated vector parameters and the initial vector parameters based on the relationship between the initial difference and the experimental difference includes: If the initial difference is greater than or equal to the experimental difference, the integrated vector parameters are determined to be the current optimal vector parameters; If the initial difference is less than the experimental difference, the initial vector parameters are determined to be the current optimal vector parameters; The current set of vector parameters is obtained based on the current optimal vector parameters: ； Where DG+1 is the current vector parameter set.

9. A modeling device for interference sources in a dipole array, characterized in that, The device includes: The detection unit is used to perform step S101: detect the magnetic field amplitude within a set range on the printed circuit board; Modeling unit, used to perform the following steps: S102: Set the number of dipoles in the dipole array to n; S103: Initialize the vector parameters of the dipole array to obtain the initial vector parameters and the initial vector parameter set; wherein, the initial vector parameter set includes the initial vector parameters. S104: Using the optimal initial vector parameters in the initial vector parameter set, perform a linear transformation and integration on the initial vector parameter set to obtain the integrated vector parameters and the experimental set; wherein, the experimental set includes the integrated vector parameters; S105: Determine the initial difference between the dipole array magnetic field amplitude corresponding to the initial vector parameter set and the magnetic field amplitude within the set range of the printed circuit board, and determine the experimental difference between the dipole array magnetic field amplitude corresponding to the experimental set and the magnetic field amplitude within the set range of the printed circuit board. S106: Based on the relationship between the initial difference and the experimental difference, determine the current optimal vector parameter and the current vector parameter set from the integrated vector parameter and the initial vector parameter; wherein, the current vector parameter set includes the current optimal vector parameter; S107: Use the current optimal vector parameters as the initial vector parameters, and use the current set of vector parameters as the initial set of vector parameters, and repeat the execution of S104-S106 a set number of times. S108: Increase the number of dipoles in the dipole array by s to form a new dipole array. Repeat S103-S107 until the number of dipoles in the new dipole array reaches the set number N. Determine the current new dipole array as the optimal dipole array and determine the vector parameters of the optimal dipole array. The initial vector parameters and the initial vector parameter set in each execution of S103 are obtained by initializing the vector parameters of the dipole array according to the optimal vector parameters determined in the previous execution of S106. Where s≥1, N>n+s. S109: Establish a printed circuit board interference source model based on the vector parameters of the optimal dipole array.

10. A modeling system for interference sources in a dipole array, characterized in that, Includes the dipole array interference source modeling device as described in claim 9.

Images