An electromagnetic simulation method and an electronic device for an Instrument Landing System

Abstract

The present invention discloses an electromagnetic simulation method and an electronic device for an instrument landing system, including the following steps: Step S100: Obtain the actual simulation area of the airport and the on-site environment of the actual simulation area, and construct an airport electromagnetic environment model; Step S200: Construct a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna; Step S300: Use the alternating direction implicit waveguide equation to construct an implicit-explicit step mathematical model for the far-source electromagnetic field, forming a tridiagonal equation; Step S400: Solve the finite-difference time-domain mathematical model and the alternating direction waveguide equation through a cloud computing platform to obtain the electric field value of the electromagnetic environment model; Step S500: Obtain the modulation difference through the electric field value. The simulation method of the present invention reduces the computational complexity, reduces the bandwidth pressure by reducing the data transmission between nodes, and enables ultra-large-scale electromagnetic simulation to operate efficiently with extremely high accuracy on the cloud platform.

Description

Technical Field

[0001] The present invention relates to the field of radio equipment, and in particular to an electromagnetic simulation method and an electronic device for an instrument landing system. Background Art

[0002] In civil aviation equipment, the instrument landing system is the most widely used precision approach and landing guidance system for aircraft. Its function is to achieve course guidance and glide path guidance through two radio signals transmitted from the ground, establish a virtual path pointing from the runway to the air, and the aircraft determines its relative position with this path through on-board receiving equipment, so that the aircraft flies towards the runway in the correct direction and smoothly descends in altitude, and finally achieves a safe landing. In the case of poor weather conditions and very low runway visibility, the instrument landing system can help the pilot guide the aircraft to approach and land. The complex terrain environment between the runway and the aircraft has a great impact on the wireless communication between the instrument landing system and the aircraft. Using large-scale electromagnetic simulation technology to predict the performance of the instrument landing system is of great significance for ensuring the safe flight of civil aviation aircraft.

[0003] With the development of computer hardware technology, large-scale electromagnetic simulation has become a hot topic in recent years. Large-scale electromagnetic simulation can be used to simulate the propagation of radio waves in a real environment to test the interference of complex media in the environment on electromagnetic waves. Using the results obtained from electromagnetic simulation, the environmental settings between the source and the target can be optimized to ensure the normal propagation of electromagnetic waves. The currently mainly adopted technology is to parallelize the finite-difference time-domain equations using a graphics processing unit (GPU) to obtain a higher speed. However, due to the limitation of computer memory, for electromagnetic simulations with a length and width of more than one hundred meters, this method has great limitations, and may even be unable to run due to the overly large simulation scale. Summary of the Invention

[0004] Aiming at the problems existing in the prior art when performing ultra-large-scale electromagnetic simulation calculations, the embodiments of the present invention provide an electromagnetic simulation method applied to an instrument landing system, which reduces the calculation complexity, reduces the data transmission between nodes, reduces the bandwidth pressure, and can perform electromagnetic simulations with a simulation range scale of more than one hundred meters in length and width.

[0005] A second object of the present invention is to provide an electronic device that executes the above-mentioned electromagnetic simulation method applied to an instrument landing system, which reduces the calculation complexity compared with the traditional calculation method, reduces the calculation complexity, reduces the data transmission between nodes, reduces the bandwidth pressure, and can perform electromagnetic simulations with a simulation range scale of more than one hundred meters in length and width.

[0006] The first object of the present invention is achieved by adopting the following technical solutions:

[0007] An electromagnetic simulation method for an instrument landing system, comprising the following steps:

[0008] Step S100: Obtain the actual simulation area of the airport and the on-site environment of the actual simulation area, and construct an airport electromagnetic environment model;

[0009] Step S200: Construct a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna;

[0010] Step S300: Use the alternating direction implicit waveguide equation to construct an implicit-explicit step mathematical model for the far-source electromagnetic field, forming a tridiagonal equation;

[0011] Step S400: Solve the finite-difference time-domain mathematical model and the alternating direction waveguide equation through a cloud computing platform to obtain the electric field value of the electromagnetic environment model;

[0012] Step S500: Obtain the modulation difference through the electric field value.

[0013] Further, the step S200 constructs a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna, including the following steps:

[0014] Step S210: Construct Maxwell's equations:

[0015]

[0016]

[0017]

[0018] where E is the electric field, H is the magnetic field, D is the electric displacement, σ is the conductivity, μ is the permeability, and ε is the permittivity;

[0019] Step S220: Use the Yee algorithm to solve Maxwell's equations and construct a three-dimensional discrete Maxwell's equations;

[0020] From Maxwell's curl equation:

[0021]

[0022]

[0023] The curl operator operation in the rectangular coordinate system is:

[0024]

[0025]

[0026] According to formula ④, it is deduced that

[0027]

[0028] According to Equation ⑤, it is derived that

[0029]

[0030] Write out the components in each direction of the vector curl equation independently, and the following scalar equations are obtained:

[0031]

[0032]

[0033]

[0034]

[0035]

[0036]

[0037] Step S230: Use the Yee grid to discretize Equation ⑩ - in space and time, and use the second-order accurate central difference approximation of the first-order derivative to divide the electric field and magnetic field equations obtained in Step S210. The discretized Maxwell's equations are:

[0038]

[0039]

[0040]

[0041]

[0042]

[0043]

[0044] Among them, the coefficients CA(m), CB(m), CP(m), and CQ(m) respectively represent the medium parameters at the positions of each field quantity, and m is the array subscript corresponding to each field quantity; the calculation formula for the medium parameters is as follows:

[0045]

[0046]

[0047]

[0048]

[0049] Further, in step S300, an implicit-explicit step mathematical model is constructed for the far-source electromagnetic field by using the alternating direction implicit pipeline equation to form a tridiagonal equation, including the following steps:

[0050] Step 310: Construct the Helmholtz equation:

[0051]

[0052] where n is the reflection coefficient, and the forward pipeline equation is solved to obtain:

[0053]

[0054] Define the following operations:

[0055]

[0056]

[0057] Using the Crank-Nicolson decomposition method, the formula is decomposed into:

[0058]

[0059] In the formula add the term

[0060]

[0061] on both the left and right sides to decompose the single-step iteration formula into implicit and explicit steps:

[0062]

[0063]

[0064] Further, in step S400, the cloud computing platform is used to solve the finite-difference time-domain mathematical model in step S200 and the alternating direction implicit pipeline equation in step S300 to obtain the electric field value of the electromagnetic environment model, including the following steps:

[0065] Step S410: Use openMPI to perform parallel computing on the discrete finite-difference time-domain equation in step S230 ;

[0066] Step S420: Use the Thomas equation to solve the diagonal equation and expand the implicit-explicit step mathematical model in step S300 into a general form:

[0067]

[0068] Perform row transformation on the coefficient matrix:

[0069]

[0070]

[0071] Formula is transformed into:

[0072]

[0073] Solve the equation according to the above steps to obtain Substitute its result into the formula Calculate and obtain according to the above method and display steps The steps to calculate the solution vector are as follows:

[0074]

[0075] Furthermore, the step S500: obtaining the modulation degree difference through the electric field value includes the following steps:

[0076] Step S510: Initialize the transmission signal in the following manner:

[0077] E CSB (t) = A c (1 + m1sinω 90 t + m2sinω 150 t)cos(ω0t + φ c )F c (θ)e -iαr ;

[0078] E SBO (t) = A s (n1sinω 90 t - n2sinω 150 t)cos(ω0t + φ s )F s (θ)e -iαr ;

[0079] Formula and formula are both the antenna transmission signals of the localizer and the glide slope, where E CSB represents the electric field strength value of the carrier and sideband signals, E SBO represents the electric field strength value of the sideband signals, A c and A srepresents the amplitude value, m1 and m2 represent the modulation depths of the CSB signal, n1 and n2 represent the modulation depths of the SBO signal, ω0 represents the carrier frequency, F c (θ) and F s (θ) represent the directivity coefficients of the CSB antenna and the SBO antenna respectively, α = 2π / λ0, λ0 is the wavelength, r represents the distance from the antenna to the receiver, L is the size of the antenna;

[0080] Step S520: After the receiver receives the CSB signal and the SBO signal, calculate the modulation depth difference through the following formula:

[0081]

[0082] where E CSB represents the electric field amplitude of the CSB signal, E SBO represents the electric field amplitude of the SBO signal, φ represents the phase difference between these two signals.

[0083] Furthermore, step S100 obtains the actual simulation area of the airport and the on-site environment of the actual simulation area, and constructs an airport electromagnetic environment model, which specifically includes the following steps:

[0084] Step S110: According to the on-site environment of the actual simulation area of the airport, use modeling software to construct a three-dimensional model of the actual simulation area;

[0085] Step S120: Perform regional grid division on the three-dimensional model; Step S130: Perform height modeling and electromagnetic parameter modeling on each grid to construct a regional electromagnetic parameter model.

[0086] Furthermore, in step S120, the grid size is 1 / 20 of the wavelength of the ILS signal.

[0087] Furthermore, the grid size is: the finite-difference time-domain calculation area is 0.05m * 0.05m * 0.05m, and the alternating direction implicit pipeline equation calculation area is 0.05m * 1.0m * 0.05m.

[0088] Furthermore, in step S110, use CAD software to construct a three-dimensional model of the actual simulation area; in step S130, call the API interface provided by CAD software for programming analysis, and model the relative permittivity, relative conductivity, and relative permeability of each grid to construct a regional electromagnetic parameter model.

[0089] The second object of the present invention is achieved through the following technical solutions:

[0090] An electronic device includes a plurality of processors, a memory, and a computer program stored on the memory and executable on the processors. When the processors execute the computer program, an electromagnetic simulation method for an instrument landing system as described in any one of the above is implemented.

[0091] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0092] The present invention provides an electromagnetic field simulation method and an electronic device for an instrument landing system. By obtaining the actual simulation area of the airport and considering the on-site environment such as terrain and building distribution when constructing the three-dimensional model of the simulation area, the constructed electromagnetic environment model is very close to the actual electromagnetic environment of the airport. At the same time, by constructing an alternating direction implicit waveguide equation to construct an implicit-explicit step mathematical model, the computational complexity is reduced, the data transmission between nodes is reduced, and the bandwidth pressure is reduced. Thus, the parallel hybrid time-domain finite difference / waveguide equation model can be processed through a cloud computing platform, so as to quickly and accurately obtain the electric field value, and thus quickly calculate the accurate modulation difference, realizing large-scale electromagnetic simulation with a simulation range scale of more than one hundred meters in length, width, and height applicable to the instrument panel landing system. BRIEF DESCRIPTION OF THE DRAWINGS

[0093] Figure 1 It is a schematic flowchart of an electromagnetic simulation method for an instrument landing system according to Embodiment 1 provided by the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0094] Next, in combination with the drawings and the specific embodiments, the present invention will be further described. It should be noted that, on the premise of no conflict, the following described embodiments or technical features can be arbitrarily combined to form new embodiments.

[0095] Embodiment 1

[0096] As Figure 1 shown, an electromagnetic method for an instrument landing system provided by an embodiment of the present invention includes the following steps:

[0097] Step S100: Obtain the actual simulation area of the airport and the on-site environment of the actual simulation area, and construct an electromagnetic environment model of the airport, which specifically includes the following steps:

[0098] Step S110: According to the on-site environment of the actual simulation area of the airport, use modeling software to construct a three-dimensional model of the actual simulation area; among them, various existing three-dimensional modeling software can be used for three-dimensional modeling, such as CAD, UG, etc. Preferably, CAD software is used because CAD software has been widely used in various infrastructure and engineering applications, providing more modeling resources. For some airports, the previously built airport CAD three-dimensional model diagram can be directly imported.

[0099] Step S120: Perform mesh division on the 3D model; for the near-source part area, the mesh size needs to be adjusted according to the boundary conditions of the obstacle. For the far-source area where the obstacle is relatively continuous, the mesh size can be appropriately increased to reduce the computational complexity. In a preferred implementation, the mesh size is specifically 1 / 20 of the blind landing signal wavelength, that is, the finite-difference time-domain calculation area is 0.05m * 0.05m * 0.05m, and the alternating direction implicit pipeline equation calculation area is 0.05m * 1.0m * 0.05m.

[0100] Step S130: Perform height modeling and electromagnetic parameter modeling on each mesh to form a regional electromagnetic parameter model. When using CAD software for modeling, call the API interface provided by the CAD software for programming analysis, and model the relative permittivity, relative conductivity, and relative permeability of each mesh to construct a regional electromagnetic parameter model.

[0101] Step S200: Construct a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna; specifically including the following steps:

[0102] Step S210: Construct Maxwell's equations:

[0103]

[0104]

[0105]

[0106] Among them, E is the electric field, H is the magnetic field, D is the electric displacement, σ is the conductivity, μ is the permeability, and ε is the permittivity;

[0107] Step S220: Use the Yee algorithm to solve the Maxwell's equations and construct a three-dimensional discrete Maxwell's equations:

[0108] For calculation on a computer, the continuous equation needs to be discretized. First, use the Yee algorithm to solve the Maxwell's equations and construct a three-dimensional discrete Maxwell's equations:

[0109] From Maxwell's curl equation:

[0110]

[0111]

[0112] The curl operator operation in the Cartesian coordinate system can be written as:

[0113]

[0114]

[0115] According to formula ④, it can be deduced that

[0116]

[0117] According to formula ⑤, it can be deduced that

[0118]

[0119] Writing out the components in each direction in the vector curl equation independently, six scalar equations can be obtained:

[0120]

[0121]

[0122]

[0123]

[0124]

[0125]

[0126] After that, using the Yee grid to discretize the above formulas ⑩ - in space and time, and using the second-order accurate central difference approximation formula of the first derivative to divide the electric field and magnetic field formulas obtained in step S210, so that the sampling points in space of the six components of the electromagnetic field in the Yee grid are placed on the edges and the center points of the surfaces of the cube, and the electromagnetic field propagates through the coupling of the electric field and the magnetic field. The discretized Maxwell's equations are:

[0127]

[0128]

[0129]

[0130]

[0131]

[0132]

[0133] Among them, the coefficients CA(m), CB(m), CP(m), CQ(m) respectively represent the medium parameters at the positions of each field quantity, and m is the array subscript corresponding to each field quantity. The calculation formulas of the medium parameters are as follows:

[0134]

[0135]

[0136]

[0137]

[0138] Step S300: Construct an alternating direction implicit pipeline equation to construct an implicit-explicit step mathematical model, forming a tridiagonal equation. Specifically, it further includes the following steps:

[0139] Step S310:: Construct the Helmholtz equation in free space:

[0140]

[0141] where n is the reflection coefficient, and the forward pipeline equation is solved as:

[0142]

[0143] To simplify the formula form, the following operations are defined:

[0144]

[0145]

[0146] Using the Crank-Nicolson decomposition method, the formula is decomposed into:

[0147]

[0148] In the formula add the term

[0149]

[0150] to both sides of the single-step iteration formula and decompose it into implicit and explicit steps:

[0151] Implicit step:

[0152]

[0153] Explicit step:

[0154]

[0155] It can be seen that the iteration formula is transformed into the form of a tridiagonal equation, which is easy to solve in a parallel system.

[0156] Step S400: Solve the finite-difference time-domain mathematical model in S200 and the alternating direction implicit pipeline equation in S300 through a cloud computing platform to obtain the electric field value of the electromagnetic environment model. Specifically, it includes the following steps:

[0157] Step S410: Use openMPI to perform parallel computing on the discrete finite-difference time-domain equation in Step S230 to solve and obtain the near-source field electric field value.

[0158] Step S420: Use the near-source field electric field value as the far-field initial value and solve the diagonal equation using the Thomas equation Specifically, it includes

[0159] Expand the implicit-explicit step mathematical model in Step S300 into a general form:

[0160]

[0161] Perform row transformation on the coefficient matrix:

[0162]

[0163]

[0164] Formula is transformed into:

[0165]

[0166] Solve the equation according to the above steps to obtain Substitute its result into the formula Calculate and obtain according to the above method for the display step The steps for calculating the solution vector are as follows:

[0167]

[0168] Step S500: Obtain the modulation difference through the electric field value. Specifically, it includes the following steps:

[0169] Step S510: Initialize the transmitted signal in the following manner:

[0170] E CSB (t) = A c (1 + m1sinω 90 t + m2sinω 150 t)cos(ω0t + φ c )F c (θ)e -iαr ;

[0171] E SBO E(t)=A s (n1sinω 90 t - n2sinω 150 t)cos(ω0t + φ s )F s (θ)e -iαr ;

[0172] Formula and formula are both the antenna emission signals of the localizer and the glide path. Among them, E CSB represents the electric field intensity value of the carrier and sideband signals, and E SBO represents the electric field intensity value of the sideband signal. A c and A s represent amplitude values, m1 and m2 represent the modulation degrees of the CSB signal, n1 and n2 represent the modulation degrees of the SBO signal, ω0 represents the carrier frequency, F c (θ) and F s (θ) represent the directivity coefficients of the CSB antenna and the SBO antenna respectively, α = 2π / λ0, λ0 is the wavelength, r represents the distance from the antenna to the receiver, L is the size of the antenna;

[0173] Step S520: When the receiver receives the CSB signal and the SBO signal, for the calculated electric field, the modulation degree difference can be calculated in the following way:

[0174]

[0175] where E CSB represents the electric field amplitude of the CSB signal, and E SBO represents the electric field amplitude of the SBO signal, and φ represents the phase difference between these two signals.

[0176] The present invention provides an electromagnetic field simulation method applied to the instrument landing system. By obtaining the actual simulation area of the airport and considering the field environment such as terrain and building distribution when constructing the three-dimensional model of the simulation area, an airport electromagnetic environment model is constructed. The constructed electromagnetic environment model is very close to the actual electromagnetic environment of the airport. At the same time, by constructing the alternating direction implicit waveguide equation to construct the implicit-explicit step mathematical model, the computational complexity is reduced, the data transmission between nodes is reduced, and the bandwidth pressure is reduced. Thus, the parallel hybrid time-domain finite difference / waveguide equation model can be processed through the cloud computing platform, so as to quickly and accurately obtain the electric field value, and thus quickly obtain the accurate modulation degree difference, realizing a large-scale electromagnetic simulation with a simulation range scale of more than one hundred meters in length, width and height applicable to the instrument panel landing system.

[0177] The electromagnetic simulation method for instrument landing system provided by Embodiment 1 has the following remarkable beneficial effects:

[0178] 1) By combining finite-difference time-domain and alternating direction implicit pipeline equation to construct a mathematical model, the simulation accuracy is ensured while the simulation efficiency is improved.

[0179] 2) When solving the pipeline equation, only three groups of data, namely r, H, and t, are transmitted between nodes, which can greatly reduce the bandwidth pressure on the cloud platform;

[0180] 3) Compared with the existing algorithms, this algorithm has the characteristics of higher efficiency and more accurate results;

[0181] 4) The experimental results show that in a certain application scenario, this algorithm is 3.4 times faster than the traditional FDTD algorithm.

[0182] Embodiment 2

[0183] Based on the same inventive concept, the present invention also provides an electronic device, including a memory, a plurality of processors, and a program stored in the memory, and the program is configured to be executed by the processors. When the plurality of processors execute the program, the steps of the above electromagnetic simulation method for instrument landing system are implemented.

[0184] The device in this embodiment and the method in the foregoing embodiment are two aspects based on the same inventive concept. The implementation process of the method has been described in detail above. Therefore, those skilled in the art can clearly understand the device and the implementation process in this embodiment according to the foregoing description. For the sake of brevity of the specification, it will not be repeated here.

[0185] The above are the preferred embodiments of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements are also regarded as the protection scope of the present invention.

Claims

1. An electromagnetic simulation method for an instrument landing system, characterized in that It includes the following steps: Step S100: Obtain the actual simulation area of the airport and the on-site environment of the actual simulation area, and construct an airport electromagnetic environment model; Step S200: Construct a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna; Step S300: Use the alternating direction implicit waveguide equation to construct an implicit-explicit step mathematical model for the far-source electromagnetic field, forming a tridiagonal equation; Step S400: Solve the finite-difference time-domain mathematical model and the alternating direction waveguide equation through a cloud computing platform to obtain the electric field value of the electromagnetic environment model; Step S500: Obtain the modulation difference through the electric field value, The step S200 constructs a finite-difference time-domain mathematical model according to the terrain within the preset range of the antenna, including the following steps: Step S210: Construct Maxwell's equations: ；① ；② ；③ Among them, is the electric field, is the magnetic field, is the electric displacement, is the conductivity, is the magnetic permeability, is the permittivity; Step S220: Use the Yee algorithm to solve Maxwell's equations and construct a three-dimensional discrete Maxwell's equations; From Maxwell's curl equation: ；④ ；⑤ The curl operator operation in the Cartesian coordinate system is: ；⑥ ；⑦ According to formula ④, it is deduced that ；⑧ According to formula ⑤, it is deduced that ；⑨ Write out each direction component in the vector curl equation independently to obtain the following scalar equations: ；⑩ ；⑪ ；⑫ ；⑬ ；⑭ ；⑮ Step S230: Discretize formulas ⑩-⑮ in space and time using a Yee grid, and use the second-order accurate central difference approximation of the first derivative to divide the electric and magnetic field formulas obtained in step S210. The discretized Maxwell's equations are: ；⑯ ；⑰ ；⑱ ；⑲ ；⑳ ；㉑ Among them, the coefficient respectively represents the medium parameters at the positions of each field quantity, m is the array subscript corresponding to each field quantity; the calculation formula of the medium parameters is as follows: ；㉒ ；㉓ ；㉔ ㉕。 2. The method according to claim 1, characterized in that, The step S300 uses the alternating direction implicit waveguide equation to construct an implicit-explicit step mathematical model for the far-source electromagnetic field, forming a tridiagonal equation, including the following steps: Step 310: Construct the Helmholtz equation: ；㉖ where n is the reflection coefficient, and the forward waveguide equation is solved to obtain: ；㉗ Define the following operations: ； ； Using the Crank-Nicolson decomposition method, formula ㉗ is decomposed into: ；㉘ Add the term on both sides of formula ㉑ ； Decompose the single-step iteration formula into implicit and explicit steps: ；㉙ ㉚。 3. The method according to claim 2, wherein The step S400 solves the finite-difference time-domain mathematical model in step S200 and the alternating direction implicit waveguide equation in step S300 through a cloud computing platform to obtain the electric field value of the electromagnetic environment model, including the following steps: Step S410: Use openMPI to perform parallel computing on the discrete finite-difference time-domain equations ⑯-㉑ in step S230; Step S420: Use the Thomas equation to solve the diagonal equation, and expand the implicit-explicit step mathematical model in step S300 into a general form: ；㉛ Perform row transformation on the coefficient matrix: ；㉜ ；㉝ Formula ㉛ is reduced to: ；㉞ Solving the equation ㉙ according to the above steps will obtain , substituting its result into formula ㉚ and calculating the display steps according to the above method will obtain , the steps for calculating the solution vector are as follows: ㉟。 4. The method according to claim 3, wherein Step S500 obtains the modulation difference through the electric field value, including the following steps: Step S510: Initialize the transmitted signal in the following manner: ；㊱ ；㊲ Both Equation ㊱ and Equation ㊲ are the antenna emission signals of the localizer and glide path, where represents the electric field strength value of the carrier and sideband signals, represents the electric field strength value of the sideband signal, and represent the amplitude value, and represent the modulation depth of the CSB signal, and represent the modulation depth of the SBO signal, represents the carrier frequency, and respectively represent the directivity coefficients of the CSB antenna and the SBO antenna, , is the wavelength, represents the distance from the antenna to the receiver, , is the antenna size; Step S520: When the receiver receives the CSB signal and the SBO signal, calculate the modulation difference through the following formula: ；㊳ wherein represents the electric field strength value of the SBO signal, represents the electric field strength value of the CSB signal, represents the electric field strength of the SBO signal, represents the phase difference between these two signals.

5. The method according to claim 1 or 4, characterized in that: The step S100 obtains the actual simulation area of the airport and the on-site environment of the actual simulation area, and constructs an airport electromagnetic environment model, specifically including the following steps: Step S110: According to the on-site environment of the actual simulation area of the airport, use modeling software to construct a three-dimensional model of the actual simulation area; Step S120: Perform regional mesh division on the three-dimensional model; Step S130: Perform height modeling and electromagnetic parameter modeling on each mesh to construct a regional electromagnetic parameter model.

6. The method according to claim 5, wherein In the step S120, the mesh size is 1 / 20 of the wavelength of the ILS signal.

7. The method according to claim 6, wherein The mesh size is: the finite-difference time-domain calculation area is 0.05m * 0.05m * 0.05m, and the alternating direction implicit pipeline equation calculation area is 0.05m * 1.0m * 0.05m.

8. The method according to claim 7, characterized in that, In step S110, use CAD software to construct a three-dimensional model of the actual simulation area; in the step S130, call the API interface provided by CAD software for programming analysis, model the relative permittivity, relative conductivity, and relative permeability of each mesh, and construct a regional electromagnetic parameter model.

9. An electronic device, characterized in that, It includes multiple processors, a memory, and a computer program stored on the memory and executable on the processors. When the processors execute the computer program, they implement the electromagnetic simulation method for an instrument landing system according to any one of claims 1 to 8.

Images