Gyromagnetic medium equivalent magnetic conductivity tensor calculation method and system and electronic equipment

Abstract

The invention relates to the technical field of electromagnetic wave regulation and control, and particularly discloses a gyromagnetic medium equivalent magnetic conductivity tensor calculation method and system and electronic device.The gyromagnetic medium equivalent magnetic conductivity tensor calculation method comprises the following steps that a composite structure composed of at least two gyromagnetic media is constructed, and the composite structure is distributed in a layered mode in the x-axis direction, the thickness of each layer of gyromagnetic medium is respectively d1 and d2, determining an electromagnetic field boundary condition at the interface of the composite structure, establishing a magnetic conductivity tensor equation of the two gyromagnetic media based on the boundary condition, and solving an equivalent magnetic conductivity tensor of the composite structure through tensor operation by combining material proportionality coefficients f1 and f2. According to the method for calculating the equivalent magnetic conductivity tensor of the gyromagnetic medium, the magnetic conductivity tensor is directly processed through the electromagnetic field boundary condition, anisotropic characteristics are compatible, the nonreciprocal effect of the gyromagnetic medium is accurately represented, and an explicit equivalent magnetic conductivity formula is exported based on the material thickness proportionality coefficient and tensor operation.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of electromagnetic wave regulation, and in particular to a method and system for calculating the equivalent permeability tensor of a gyromagnetic medium, and an electronic device. BACKGROUND

[0002] Gyromagnetic medium plays a key role in electromagnetic wave regulation (such as isolator and circulator design) due to its permeability tensor characteristics (including off-diagonal elements). Traditional theory considers that the permeability is constant at optical frequencies. Although the emergence of metamaterials has achieved the regulation of diagonal elements, the efficient calculation and regulation of off-diagonal elements (gyromagnetic characteristics) still face bottlenecks, which are core challenges in the field of optoelectronics.

[0003] Existing methods cannot directly handle the off-diagonal components of the permeability tensor when calculating the equivalent permeability tensor, resulting in the lack of analysis of gyromagnetic characteristics (such as non-reciprocity). In the calculation process, the continuity conditions of the normal magnetic induction intensity (Bx) and the tangential magnetic field intensity (Hy, Hz) at the material interface are not fully utilized, and the model deviates from the actual field distribution. SUMMARY

[0004] The present application aims to at least partially solve one of the technical problems in the related art. To this end, the present application aims to provide a method and system for calculating the equivalent permeability tensor of a gyromagnetic medium, and an electronic device, to efficiently and accurately calculate the equivalent permeability tensor of a gyromagnetic medium containing off-diagonal elements.

[0005] To achieve the above-mentioned purpose, the first aspect of the present application provides a method for calculating the equivalent permeability tensor of a gyromagnetic medium, which comprises the following steps:

[0006] A composite structure composed of at least two gyromagnetic media is constructed, which is distributed in a layered manner along the x-axis direction, and the thickness of each layer of gyromagnetic medium is d1, d2, respectively, and the total thickness d = d1 + d2, the proportionality coefficient of the first material is defined as The proportionality coefficient of the second material is

[0007] The electromagnetic field boundary conditions at the interface of the composite structure are determined: the normal magnetic induction intensity B x is continuous: B 1x = B 2x = B x ; the tangential magnetic field intensity H y and H z are continuous: H 1y = H2 y = H y , H 1z = H 2z = H z ;

[0008] Based on the boundary conditions, the permeability tensor equations of two gyromagnetic media are established, and the permeability tensors are respectively:

[0009] Wherein, wherein u1 and u2 are the permeability parameters of two gyromagnetic media respectively, and i is the imaginary unit;

[0010] Combined with the material proportion coefficients f1 and f2, the equivalent permeability tensor μ of the composite structure is solved by tensor operation, and the expression is:

[0011]

[0012] In some embodiments of the present application, the material proportion coefficients f1 and f2 satisfy:

[0013] f1+f2=1.

[0014] In some embodiments of the present application, the normal magnetic induction intensity continuous in the boundary conditions is specifically:

[0015] B 1x =B x =μ1H 1x -iμ1'H y

[0016] B 2x =B x =μ2H 2x -iμ'2H y ;

[0017] Solving:

[0018] In some embodiments of the present application, the transverse component of the equivalent permeability tensor is solved by linear superposition, and specifically:

[0019] μ xx =f1μ 1xx +f2μ 2xx

[0020] Wherein, wherein μ 1xx and μ 2xx are the permeability components of two gyromagnetic media in the x-axis direction, and μ 1xx =0, μ 2xx =0.

[0021] In some embodiments of the present application, the longitudinal component of the equivalent permeability tensor is solved by harmonic mean, and specifically:

[0022]

[0023] In some embodiments of the present application, the method assumes that the external magnetic field is applied along the z-axis direction, and the magnetic field component H z is continuous.

[0024] In some embodiments of the present application, the method can be extended to a composite structure containing N layers of gyromagnetic media, wherein the material proportionality coefficient of the kth layer is and the total thickness

[0025] In some embodiments of the present application, the method further combines nonlinear magnetic response analysis, including the magnetic hysteresis effect of ferromagnetic materials, by introducing a nonlinear correction term Δu to correct the equivalent permeability tensor.

[0026] To achieve the above-mentioned purpose, the second aspect of the present application proposes a gyromagnetic medium equivalent permeability tensor calculation system, comprising:

[0027] Model construction module: for defining a layered composite structure composed of at least two gyromagnetic media, inputting the thickness d1, d2 of each layer and calculating the proportionality coefficient

[0028] Boundary condition module: for setting the electromagnetic field continuity condition at the interface, i.e. the normal B x is continuous and the tangential H y is continuous; z

[0029] Tensor operation module: based on the boundary conditions and the proportionality coefficient, the equivalent permeability tensor is derived through analytical formula, and the tensor matrix containing off-diagonal components is output.

[0030] To achieve the above-mentioned purpose, the third aspect of the present application proposes an electronic device, comprising a memory, a processor and a computer program stored in the memory, wherein the computer program is executed by the processor to realize the above-mentioned gyromagnetic medium equivalent permeability tensor calculation method.

[0031] The gyromagnetic medium equivalent permeability tensor calculation method, system and electronic device of the embodiments of the present application directly process the permeability tensor through the electromagnetic field boundary conditions, are compatible with anisotropic characteristics, accurately characterize the non-reciprocal effect of gyromagnetic media, and based on the material thickness proportionality coefficient and tensor operation, derive an explicit equivalent permeability formula, which significantly improves the calculation efficiency. BRIEF DESCRIPTION OF DRAWINGS

[0032] Figure 1 is a flowchart of a gyromagnetic medium equivalent permeability tensor calculation method according to an embodiment of the present application;

[0033] Figure 2 is a schematic diagram of the composite structure constructed according to an embodiment of the present application;​

[0034] Figure 3 Fig. 1 is a structural schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0035] Embodiments of the present application are described in detail below with reference to examples illustrated in the accompanying drawings, in which the same or similar components are denoted by the same or similar reference numerals throughout. The embodiments described below by reference to the drawings are exemplary and are intended to explain the present application, and should not be understood as limiting the present application.

[0036] A method, system and electronic device for calculating the equivalent permeability tensor of a gyromagnetic medium are described below with reference to the accompanying drawings.

[0037] Figure 1 Fig. 1 is a structural schematic diagram of an electronic device according to an embodiment of the present application.

[0038] As shown in Fig. 1, a method for calculating the equivalent permeability tensor of a gyromagnetic medium includes the following steps: Figure 1

[0039] 1. Composite structure construction and parameter definition

[0040] As shown in Fig. 2, the geometric model: a composite structure is constructed along the x-axis direction, which contains at least two gyromagnetic media (denoted as medium 1 and medium 2), and the thicknesses of the two media are d1 and d2, respectively, and the total thickness d = d1 + d2. Figure 2 Proportionality coefficient: the proportionality coefficient of the first material is defined as f1 = d1 / d, and the proportionality coefficient of the second material is defined as f2 = d2 / d.

[0041] f1 + f2 = 1, representing the volume ratio of the two media in the composite structure.

[0042] Parameter definition supplement:

[0043] d1, d2: physical thickness of each layer of gyromagnetic medium, unit: meter (m);

[0044] f1, f2: dimensionless coefficient, used to quantify the material ratio, support dynamic adjustment of the anisotropy degree of equivalent permeability.

[0045] 2. Electromagnetic field boundary condition is clear

[0046] Boundary condition refinement:

[0047] Normal continuity: the normal magnetic induction intensity B x is continuous at the interface, i.e. B 1x n1 = B 2x n2 = B x ​​​The normal direction of medium 1 and medium 2

[0048] B 1x = B x = μ1H 1x -iμ1'H y

[0049] The magnetic induction intensity is equal). Wherein, B 2x = B x = μ2H 2x -iμ'2H y , which embodies the coupling effect of the off-diagonal component of the permeability tensor to the normal magnetic field.

[0050] Tangential continuity: the tangential magnetic field intensity H y and H z are continuous, that is, H 1y = H2 y = H y , H 1z = H 2z = H z , the continuity of the tangential magnetic field ensures the consistency of the field distribution across the medium, which is the key constraint for deriving the equivalent permeability.

[0051] 3. Permeability tensor equation

[0052] Tensor form is specified:

[0053] The permeability tensors of medium 1 and medium 2 are respectively:

[0054]

[0055] Wherein, u1 and u2 are the permeability parameters of the two kinds of paramagnetic medium respectively (characterizing isotropic magnetic response), μ1', μ2' are off-diagonal permeability components (characterizing paramagnetic effect, such as Faraday rotation); The existence of off-diagonal element ±iμ' embodies the coupling between magnetic field components (such as the interaction of H x and H y ), i is the imaginary unit.

[0056] 4. Solution of equivalent permeability tensor

[0057] Using the continuity of B x , H y and H z , the equation group is established:

[0058] B 1x = B x = μ1H 1x -iμ1'H y

[0059] B 2x = Bx = μ2H 2x - iμ'2H y

[0060] Solving:

[0061]

[0062] Substitute:

[0063]

[0064] Solving:

[0065]

[0066] B y = f1B 1y + f2B 2y = f1(iμ1'H 1x + μ1H 1y ) + f2(iμ'2H 2x + μ2H 2y )

[0067] = (f1iμ1'H 1x + f2iμ'2H 2x ) + (f1μ1 + f2μ2)H y

[0068] Combining the material proportionality coefficients f1, f2, solve for the equivalent permeability components. From the system of equations: B 1x = B x = μ1H 1x - iμ1'H y

[0069] B 2x = B x = μ2H 2x - iμ'2H y

[0070] Resulting in:

[0071] μ1H 1x - μ2H 2x = i(μ1' - μ'2)H y

[0072] Combining f1H 1x + f2H 2x = H x

[0073] Solving:

[0074]

[0075] Substitute:

[0076]

[0077] B z = f1B 1z + f2B 2z = f1μ1H z + f2μ2H z = (f1μ1+ f2μ2)H z

[0078] Finally, the equivalent permeability tensor μ is obtained, and its expression is:

[0079]

[0080] In some embodiments of the present application, the transverse component of the equivalent permeability tensor is solved by linear superposition, specifically:

[0081] μ xx = f1μ 1xx + f2μ 2xx

[0082] wherein μ 1xx and μ 2xx represent the permeability components of the two gyromagnetic media in the x-axis direction, and μ 1xx = 0, μ 2xx = 0. This is due to the characteristics of gyromagnetic media, which determine that in the x-axis direction, the permeability component exhibits special physical properties (such as being affected by the direction of the external magnetic field and the intrinsic properties of the medium, resulting in zero permeability component in this direction), and the linear superposition can effectively handle the continuity of the non-diagonal components of the tensor. Specifically, this calculation method is based on the contribution of the two media to the transverse magnetic field response in the composite structure, and f1 and f2 are the material thickness proportionality coefficients, which quantify the proportion of the two media in the total thickness, so that the calculation of the transverse component not only considers the characteristics of the medium itself, but also reflects the influence of the structure proportion on the equivalent permeability.

[0083] Specifically, by defining μ 1xx = 0 and μ 2xx = 0 and using linear superposition, the non-diagonal component continuity problem in the permeability tensor of gyromagnetic media can be directly handled, and the characteristics of gyromagnetic media under the action of transverse magnetic field can be accurately described, providing a more accurate model for analyzing the propagation of electromagnetic waves in the composite structure. Based on the linear superposition of material thickness proportionality coefficients, complex numerical calculations or simulations are avoided, improving the calculation efficiency, and at the same time, the formula has a clear physical meaning, which is convenient for theoretical analysis and engineering application.

[0084] In some embodiments of the present application, the longitudinal component of the equivalent permeability tensor is solved by harmonic mean, specifically:

[0085]

[0086] The physical meaning of the formula is that in the longitudinal (z-axis direction), the distribution characteristics of the magnetic field make the permeability of the two gyromagnetic media have a series effect on the equivalent result. The use of harmonic mean can more accurately reflect this series relationship and avoid errors caused by simple arithmetic mean. Specifically, μ 1zz and μ 2zz are the permeability components of the two gyromagnetic media in the z-axis direction, and f1 and f2 still represent the material thickness proportionality coefficient. Through this calculation method, the hindering effect of the two media on the magnetic field in the longitudinal direction can be considered comprehensively, so that the calculation of the equivalent permeability is more in line with the actual physical law.

[0087] Specifically, the use of harmonic mean is more in line with the actual situation of the longitudinal magnetic field distribution, can accurately simulate the comprehensive effect of the two gyromagnetic media on the magnetic field in the longitudinal direction, and can reflect the series effect of the medium permeability better than other average methods (such as arithmetic average), thereby improving the accuracy of the equivalent permeability calculation. This calculation method makes the model better adapt to the combination of gyromagnetic media with different thickness proportions, provides a flexible and accurate method for designing composite structures with specific longitudinal magnetic properties, and expands the application range of the scheme in practical engineering.

[0088] In some embodiments of the present application, it is assumed that the external magnetic field is applied along the z-axis direction, which is based on the typical scenario of gyromagnetic media in practical applications (such as the magnetic field configuration in microwave devices). When the external magnetic field is along the z-axis direction, the permeability tensor characteristics of the gyromagnetic medium will produce a specific response to the external magnetic field, directly affecting the propagation characteristics of the electromagnetic wave. Under this setting, the magnetic field component H z is continuous, which ensures that the magnetic field strength in the z direction remains consistent at the interface of each layer of the composite structure, in line with the continuity principle of physical fields. Specifically, the continuity of H z allows accurate reflection of the joint action of the gyromagnetic media in each layer on the magnetic field in the z direction when calculating the equivalent permeability tensor, avoiding calculation errors caused by magnetic field discontinuity, providing a key boundary condition for accurate solution of the equivalent permeability tensor, and enhancing the practicality and reliability of the method of the present application.

[0089] In some embodiments of the present application, the method of the present application can be flexibly extended to a composite structure containing N layers of gyromagnetic media. For the kth layer of gyromagnetic medium, its material proportionality coefficient is defined as where d k is the thickness of the kth layer, and the total thickness

[0090] Taking a three-layer structure (N = 3) as an example, the thicknesses of the three layers are d1, d2, and d3, and the total thickness d = d1 + d2 + d3, then f1 = d1 / d, f2 = d2 / d, and f3 = d3 / d. In calculating the equivalent permeability tensor, the transverse component can be expressed as The longitudinal component is expressed as This extension can adapt to complex multi-layer structure design and meet various engineering needs. For example, in a broadband electromagnetic wave control device, by reasonably designing the thickness and material properties of the multi-layer gyromagnetic medium, precise control of electromagnetic waves of different frequencies can be achieved.

[0091] In some embodiments of the present application, the present application further combines nonlinear magnetic response analysis. For the magnetic hysteresis effect of ferromagnetic materials, the equivalent permeability tensor is modified by introducing a nonlinear correction term Δu. The magnetic hysteresis effect refers to the fact that the magnetization of a ferromagnetic material is related not only to the current magnetic field strength but also to the historical magnetization state, resulting in nonlinear changes in permeability.

[0092] The specific modification method is to add a correction term related to the magnetic field strength to the original equivalent permeability component. For example, for the longitudinal component μ zz , the modified expression is where Δu zz is related to the magnetic field strength H. The relationship can be determined by experimental data or theoretical model fitting. Through this modification, the magnetic properties of ferromagnetic materials in strong magnetic fields or complex magnetic field environments can be more accurately described, and the accuracy of equivalent permeability tensor calculation in nonlinear scenarios can be improved.

[0093] The above method expands and improves the simulation capability of complex magnetic property materials in actual engineering, making the equivalent permeability tensor calculation more consistent with the real physical situation, providing a more accurate theoretical tool for the design of high-reliability gyromagnetic devices, and reducing the performance deviation of the device caused by model simplification.

[0094] Corresponding to the above method embodiments, the present application also proposes a gyromagnetic medium equivalent permeability tensor calculation system, comprising:

[0095] A model construction module: this module is used to define a layered composite structure composed of at least two gyromagnetic media, has a user interaction interface, and can realize the input of parameters such as the thickness d1, d2 of each layer (supports manual input or file import). After input is completed, the system automatically calculates the proportion coefficient (where d = d1 + d2), and displays the structure information in a visual graph (such as a layered structure diagram, labeling the thickness and material type of each layer).

[0096] Meanwhile, the module supports multi-layer structure expansion (up to N layers), and users can flexibly adjust the number of layers and corresponding parameters to meet different complex structure requirements. Users input parameters through the interactive interface, and the system automatically calculates the scale factor and displays the structure, reducing manual calculation and drawing time, improving design efficiency, and the module supports multi-layer structure expansion, meeting the diversified engineering design requirements and enhancing the applicability of the system.

[0097] Boundary condition module: This module is used to set the electromagnetic field continuity conditions at the interface, including:

[0098] Normal B x Continuity: Ensure that the normal magnetic induction intensity B x is consistent at the interface of each layer of the composite structure, that is, B 1x = B 2x = B x (B 1x , B 2x are the normal magnetic induction intensities at the interface of the two adjacent layers of medium), which is based on the magnetic field Gauss law in the Maxwell equation set to ensure that there is no flux source at the interface of the magnetic field.

[0099] Tangential H y , H z Continuity: Ensure that the tangential magnetic field intensity H y , H z is continuous at the interface, that is, H 1y = H 2y = H y , H 1z = H 2z = H z (H 1y , H 2y are the tangential magnetic field intensities at the interface of the two adjacent layers of medium), which is derived from the tangential component of the Ampere loop law at the interface, ensuring that the tangential component of the magnetic field has no abrupt change.

[0100] The module encapsulates these conditions internally, and users can directly select and apply them without deep understanding of complex electromagnetic theory, reducing the use threshold.

[0101] Tensor operation module: Based on the boundary conditions and scale factors, the module derives the equivalent permeability tensor according to the Maxwell equation set and the definition of the permeability tensor of the gyromagnetic medium.

[0102] The specific process is as follows: call the scale factors f1, f2 and the permeability tensor parameters of each layer of the model construction module;

[0103] According to the B x , H y , H zContinuity, simultaneous equations are solved. For example, for the transverse component μ xx , by linear superposition μ xx = f1μ 1xx + f2μ 2xx (μ 1xx = 0, μ 2xx = 0).

[0104] For the longitudinal component μ zz , the harmonic mean is adopted

[0105] The final output contains the complete tensor matrix of the off-diagonal component, and supports the visualization and export of the matrix.

[0106] As an example, the system can also add a data verification module to increase the data verification function, to check the rationality of the input parameters such as the thickness d1, d2, etc. (such as whether the thickness is greater than zero, whether the sum of the proportion coefficients is 1, etc.), and if the parameters are abnormal, the system will prompt the user to correct in time, avoiding calculation errors due to input errors.

[0107] As an example, the system can also add a multi-coordinate system support module to expand support for cylindrical coordinate system, spherical coordinate system, etc. in addition to the rectangular coordinate system, to meet the equivalent permeability tensor calculation needs of different shapes of gyromagnetic medium structures, such as cylindrical gyromagnetic resonant cavity, spherical gyromagnetic particles, and further broaden the application field of the system.

[0108] Corresponding to the above embodiments, the application also provides an electronic device.

[0109] As shown in Figure 3 is a structural schematic diagram of an electronic device in the application, the electronic device 200 comprises a processor 201 and a memory 203. Wherein, the processor 201 and the memory 203 are connected, such as connected through a bus 202. Optionally, the electronic device 200 can also comprise a transceiver 204. It should be noted that in actual application, the transceiver 204 is not limited to one, and the structure of the electronic device 200 does not constitute a limitation on the embodiments of the application.

[0110] The processor 201 can be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array) or other programmable logic device, a transistor logic device, a hardware component, or any combination thereof. It can implement or execute various exemplary logical blocks, modules and circuits described in combination with the present disclosure. The processor 201 can also be a combination of computing functions, such as one or more microprocessor combinations, combinations of DSP and microprocessor, etc.

[0111] The bus 202 can include a path for transmitting information between the above-mentioned components. The bus 202 can be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. The bus 202 can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 3 Only one thick line is used in the middle, but it does not mean that there is only one bus or one type of bus.

[0112] The memory 203 is used to store a computer program corresponding to the method for calculating the equivalent permeability tensor of a gyromagnetic medium according to the above-mentioned embodiments of the present application. The computer program is controlled and executed by the processor 201. The processor 201 is used to execute the computer program stored in the memory 203 to realize the content shown in the above-mentioned method embodiments.

[0113] The electronic device 200 includes but is not limited to mobile terminals such as notebook computers, PADs (tablet computers), etc., and fixed terminals such as desktop computers, etc. Figure 3 The electronic device 200 shown is only an example and should not limit the functions and use range of the embodiments of the present application.

[0114] The electronic device 200 of the embodiments of the present application directly processes the permeability tensor based on the electromagnetic field boundary condition, is compatible with the anisotropic characteristics, accurately represents the non-reciprocal effect of the gyromagnetic medium, and derives an explicit equivalent permeability formula based on the material thickness scaling coefficient and tensor operation, so that the calculation efficiency is significantly improved.

[0115] It is to be appreciated that the above description and the examples that follow are intended to be illustrative only and that changes can be made to the description, as represented by the above listed elements, by the steps recited in the flow charts, and by the examples that follow, without departing from the spirit of the application. Accordingly, the scope of the present application is intended to be defined only by the appended claims.

[0116] It should be understood that aspects of the application can be implemented in hardware, software, firmware or a combination thereof. In the above embodiments, various steps or methods can be implemented in software or firmware that is stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any of the following techniques, which are well known in the art of hardware implementation, can be used: a hybrid of the above techniques, a mixture of two or more of the above techniques, or a combination of the above techniques with other techniques not listed above.

[0117] In the description of the present specification, the description of the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" and the like means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Also, the specific features, structures, materials or characteristics described can be combined in any appropriate manner in one or more embodiments or examples.

[0118] Furthermore, the terms "first", "second", etc. are used only for descriptive purposes and do not connote or imply relative importance or a quantity of the indicated technical features. Thus, a feature defined with "first", "second", etc. can include at least one of the features implicitly or explicitly. In the description of the present application, the meaning of "a plurality" is at least two, for example, two, three, etc., unless otherwise specifically defined.

[0119] Although the embodiments of the present application have been shown and described above, it is understood that the above-described embodiments are exemplary and are not to be construed as limiting the present application, and that changes, modifications, substitutions and variations can be made by those skilled in the art without departing from the scope of the present application.

Claims

1. A method for calculating the equivalent permeability tensor of a gyromagnetic medium, characterized in that: The following steps are involved: Construct a composite structure composed of at least two gyromagnetic media. The composite structure is layered along the x-axis. The thickness of each layer of gyromagnetic media is d1 and d2 respectively, and the total thickness is d = d1 + d2. The proportional coefficient of the first material is defined as The proportionality coefficient of the second material is Clarify the electromagnetic field boundary conditions at the interface of the composite structure: normal magnetic induction intensity B x Continuous: B 1x =B 2x =B x ; Tangential magnetic field strength H y and H z Continuous: H 1y =H 2y =H y , H 1z =H 2z =H z ; Based on the boundary conditions, the magnetic permeability tensor equations of the two gyromagnetic media are established, and the magnetic permeability tensors are: Where, u1 and u2 are the magnetic permeability parameters of the two gyromagnetic media, respectively, and i is an imaginary unit; Combined with the material proportional coefficients f1 and f2, the equivalent magnetic permeability tensor of the composite structure is solved by tensor operation. Its expression is:

2. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, characterized in that: The material proportional coefficients f1 and f2 satisfy: f1+f2=1.

3. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, characterized in that: The boundary condition for the continuity of normal magnetic induction intensity is specifically: B 1x =B x =μ1H 1x -iμ1'H y B 2x =B x =μ2H 2x -iμ'2H y ; The solution is:

4. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, wherein: The transverse component of the equivalent permeability tensor is solved by linear superposition, specifically: m xx =f1μ 1xx +f2μ 2xx Among them, μ 1xx and μ 2xx are the magnetic permeability components of the two gyromagnetic media in the x-axis direction, and μ 1xx =0, μ 2xx =0.

5. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, wherein: The longitudinal component of the equivalent permeability tensor is solved by harmonic averaging, specifically:

6. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, characterized in that: The method assumes that the external magnetic field is applied along the z-axis and the magnetic field component H z continuous.

7. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, characterized in that: The method can be extended to a composite structure comprising N layers of gyromagnetic media, where the material ratio coefficient of the kth layer is And the total thickness 8. The method for calculating the equivalent permeability tensor of a gyromagnetic medium according to claim 1, wherein: The method further combines nonlinear magnetic response analysis, including the hysteresis effect of ferromagnetic materials, and modifies the equivalent permeability tensor by introducing a nonlinear correction term Δu.

9. A gyromagnetic medium equivalent permeability tensor calculation system, characterized in that: include: Model building module: used to define a layered composite structure consisting of at least two gyromagnetic media, input the thickness of each layer d1, d2 and calculate the proportional coefficient Boundary condition module: used to set the electromagnetic field continuity condition at the interface, that is, normal B x Continuous and tangential H y 、H z continuous; Tensor operation module: Based on the boundary conditions and the proportional coefficient, the equivalent permeability tensor is derived through an analytical formula, and a tensor matrix containing non-diagonal components is output.

10. An electronic device, characterized in that: The invention comprises a memory, a processor and a computer program stored in the memory, wherein when the computer program is executed by the processor, a method for calculating the equivalent magnetic permeability tensor of a gyromagnetic medium according to any one of claims 1 to 8 is implemented.

Images