A simulation method, system, device and storage medium for aviation radio signal jamming system

Abstract

The present invention provides a simulation method, system, device and storage medium for an aviation radio signal jamming system. The method comprises the following steps: obtaining transmission signal parameters of a target navigation station, a three-dimensional spatial coordinate set of surrounding terrain obstacles and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include a rain intensity parameter and a fog droplet concentration parameter; calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter; analyzing the horizontal azimuth and vertical elevation angles of the surrounding terrain obstacles relative to the target navigation station according to the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set, extracting the main reflection direction therefrom, and generating a directional reflection jamming signal by combining the directional reflection jamming signal with the multipath scattering phase offset modulation; generating a course deviation simulation signal for the aviation radio signal jamming system using the meteorological attenuation factor, the directional reflection jamming signal and the transmission signal parameters; and improving the authenticity and accuracy of the course deviation simulation signal in complex electromagnetic scenarios.

Description

Technical Field

[0001] The present invention relates to the field of electrical communication technology, and in particular to a simulation method, system, device and storage medium for an aviation radio signal jamming system. Background Art

[0002] The increasing complexity and sophistication of aviation navigation systems has placed higher demands on simulation technology for radio signal jamming systems. Especially during the landing guidance phase, the stability of the localizer signal directly impacts the aircraft's navigation accuracy and operational safety. Therefore, high-precision simulation methods are required to evaluate the performance of navigation systems in jamming environments.

[0003] Existing solutions employ simulation methods based on multipath channel modeling and the fusion of a terrain database. These methods collect the target navigation station's transmission parameters, combine them with a digital elevation model to extract terrain reflection path information, and use statistical methods to estimate the arrival angle and amplitude distribution of the multipath signal. Multipath interference signals are generated using real-time signal synthesis technology and superimposed on the original navigation signal to simulate signal deviations that may occur during actual flight. However, this existing solution has significant shortcomings in handling the effects of complex meteorological conditions and terrain coupling. For example, its multipath modeling relies primarily on static terrain data and fails to fully account for the dynamic modulation of signal scattering characteristics by meteorological factors such as rain and fog. This results in deviations from the actual phase characteristics and spatial distribution of the generated interference signal. Furthermore, a lack of refined modeling for identifying the main reflection direction and distributing interference energy makes it difficult to accurately reproduce the effects of strong interference in specific directions, limiting its applicability for high-fidelity aviation interference simulations. Summary of the Invention

[0004] The present invention provides a simulation method, system, device and storage medium for an aviation radio signal jamming system, which are used to solve the problems that the interference signals generated in the prior art deviate from the actual situation in terms of phase characteristics and spatial distribution; and are insufficiently applicable in high-fidelity aviation jamming simulation. This improves the realism and accuracy of the course deviation simulation signal in complex electromagnetic scenarios.

[0005] In a first aspect, the present invention provides a method for simulating an aviation radio signal jamming system, comprising:

[0006] Acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of surrounding terrain obstacles, and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include rain intensity parameters and fog droplet concentration parameters;

[0007] Calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter;

[0008] Analyzing the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set;

[0009] Extracting the main reflection direction from the terrain reflection azimuth vector set, and combining it with the multipath scattering phase offset modulation to generate a directional reflection interference signal;

[0010] Based on the meteorological attenuation factor, the directional reflection interference signal and the transmission signal parameter, a localizer deviation simulation signal for an aviation radio signal jamming system is generated.

[0011] Optionally, calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter includes:

[0012] Based on the rainfall intensity parameter, generating an instantaneous rainfall intensity value sequence, and calculating the fluctuation amplitude of the instantaneous rainfall intensity value sequence;

[0013] Based on the proportional relationship between the fluctuation amplitude and the preset rainfall intensity benchmark, a rainfall intensity persistence coefficient is determined, and combined with a preset atmospheric attenuation constant, a meteorological attenuation factor is generated;

[0014] Analyzing the spatial distribution characteristics of the droplet concentration parameters to identify the three-dimensional position of the droplet concentration peak area;

[0015] The spatial impact weight of the droplet is calculated according to the straight-line distance between the three-dimensional position and the target navigation station, and the multipath scattering phase offset is generated in combination with a preset phase disturbance base value.

[0016] Optionally, according to the three-dimensional space coordinate set, analyzing the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station to generate a terrain reflection azimuth vector set includes:

[0017] Calculating a coordinate difference component set between the coordinates of all obstacles and the coordinates of the target navigation station in the three-dimensional space coordinate set, wherein the coordinate difference component set includes a horizontal plane first direction difference, a horizontal plane second direction difference, and a vertical direction difference;

[0018] determining the horizontal azimuth of the surrounding terrain obstacle according to the ratio of the first horizontal plane direction difference to the second horizontal plane direction difference;

[0019] Calculating the square root of the sum of the squares of the first horizontal plane direction difference and the second horizontal plane direction difference, and using the square root of the sum of the squares as the horizontal projection length;

[0020] determining the vertical elevation angle of the surrounding terrain obstacle based on a proportional relationship between the vertical direction difference and the horizontal projection length;

[0021] The horizontal azimuth angle and the vertical elevation angle are combined into an azimuth elevation angle vector, and the azimuth elevation angle vectors of all surrounding terrain obstacles are collected to form a terrain reflection azimuth vector set.

[0022] Optionally, extracting a main reflection direction from the terrain reflection azimuth vector set and combining it with the multipath scattering phase offset modulation to generate a directional reflection interference signal includes:

[0023] Calculating obstacle distance parameters and size parameters corresponding to each azimuth elevation angle vector in the terrain reflection azimuth vector set;

[0024] Calculating a distance attenuation factor based on the obstacle distance parameter and a reflection gain factor based on the size parameter, and generating a directional strength weight value of each elevation angle vector in combination with the distance attenuation factor;

[0025] Selecting the azimuth elevation angle vector with the largest directional intensity weight value as the main reflection direction, and converting the multipath scattering phase offset into a phase rotation angle;

[0026] generating a basic reflected signal based on the transmit signal parameter, and performing a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal;

[0027] The phase disturbance reflection signal is associated and bound with the spatial orientation of the main reflection direction to generate a directional reflection interference signal.

[0028] Optionally, a distance attenuation factor is calculated according to the obstacle distance parameter, and a reflection gain factor is calculated according to the size parameter. The directional strength weight values ​​of each orientation elevation angle vector are generated in combination with the distance attenuation factor, including:

[0029] Performing an inverse proportional operation on the obstacle distance parameter to generate a distance attenuation factor;

[0030] Comparing the size parameter with a preset size threshold, when the size parameter is less than the preset size threshold, multiplying the size parameter by a unity gain coefficient to obtain a reflection gain factor, and when the size parameter is greater than or equal to the preset size threshold, using a preset fixed gain value as the reflection gain factor;

[0031] Multiplying the distance attenuation factor by the reflection gain factor to generate an original weight coefficient;

[0032] Obtain the obstacle distribution density in the space sector corresponding to each azimuth elevation vector;

[0033] All obstacle distribution densities and the original weight coefficients are weighted to generate directional strength weight values ​​of all orientation elevation angle vectors.

[0034] Optionally, generating a basic reflected signal based on the transmit signal parameter, and performing a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal includes:

[0035] Resolving the carrier frequency value and the initial amplitude value from the transmission signal parameters;

[0036] Determining a reflection intensity correction coefficient based on an obstacle distance parameter and a size parameter corresponding to the main reflection direction;

[0037] Performing a product operation on the initial amplitude value and the reflection intensity correction coefficient to generate a reflection signal amplitude value;

[0038] generating a sinusoidal carrier signal based on the carrier frequency value and the reflected signal amplitude value, and using the sinusoidal carrier signal as a basic reflected signal;

[0039] The phase rotation angle is superimposed on the phase component of the sinusoidal carrier signal to obtain a phase disturbance reflection signal.

[0040] Optionally, generating a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal, and the transmission signal parameter includes:

[0041] Converting the transmission signal parameters into a reference carrier signal, and performing amplitude modulation on the reference carrier signal to generate a reference navigation signal;

[0042] Performing waveform superposition of the directional reflection interference signal and the reference navigation signal to generate a synthetic interference signal;

[0043] Converting the meteorological attenuation factor into a signal attenuation proportional coefficient, and performing a product operation on the signal attenuation proportional coefficient and the amplitude component of the synthetic interference signal to generate an attenuated interference signal;

[0044] Obtaining a standard localizer reference azimuth, and calculating an azimuth deviation angle between the main reflection direction and the standard localizer reference azimuth;

[0045] A signal phase offset is determined according to the azimuth deviation angle, so as to perform offset coding processing on the attenuated interference signal and generate a course deviation simulation signal.

[0046] In a second aspect, the present invention provides a simulation system for an aviation radio signal jamming system, comprising:

[0047] An acquisition module is used to acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of the surrounding terrain obstacles, and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include rain intensity parameters and fog droplet concentration parameters;

[0048] A calculation module, configured to calculate a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter;

[0049] an analysis module, configured to analyze the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set;

[0050] An extraction module, configured to extract a main reflection direction from the terrain reflection azimuth vector set, and generate a directional reflection interference signal by combining the multipath scattering phase offset modulation;

[0051] A generating module is used to generate a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal and the transmission signal parameter.

[0052] In a third aspect, the present invention provides a computing device comprising a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component to implement a simulation method of an aviation radio signal interference system as described in the first aspect above.

[0053] In a fourth aspect, the present invention provides a computer storage medium storing a computer program, which, when executed by a computer, implements a simulation method of an aviation radio signal interference system as described in the first aspect.

[0054] In the present invention, the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of the surrounding terrain obstacles and the rain and fog meteorological parameters are obtained, and the rain and fog meteorological parameters include a rain intensity parameter and a fog droplet concentration parameter; based on the rain intensity parameter and the fog droplet concentration parameter, the meteorological attenuation factor and the multipath scattering phase offset are calculated; according to the three-dimensional spatial coordinate set, the horizontal azimuth and vertical elevation angle of the surrounding terrain obstacles relative to the target navigation station are analyzed to generate a terrain reflection azimuth vector set, and the main reflection direction is extracted therefrom, and a directional reflection interference signal is generated by combining with the multipath scattering phase offset modulation; based on the meteorological attenuation factor, the directional reflection interference signal and the transmission signal parameters, a course deviation simulation signal for an aviation radio signal jamming system is generated. The technical solution provided by the present invention solves the problem of lack of dynamic coupling caused by the isolated processing of meteorological and terrain data in traditional solutions; by quantifying the dynamic modulation effect of rain and fog on signal scattering characteristics, it overcomes the defect of insufficient modeling of the influence of meteorological factors in existing solutions; it realizes dynamic analysis of terrain reflection paths, breaking through the adaptability limitations of static terrain databases to complex fly-around scenarios; through precise identification and phase modulation of the main reflection direction, it reproduces the strong interference effect in a specific direction and solves the problem of interference energy distribution distortion; by combining environmental attenuation and directional interference, it reproduces the distortion process of the actual heading channel signal, meeting the needs of high-precision navigation performance evaluation. Furthermore, based on the terrain reflection azimuth vector set, the reflection energy weight is dynamically calculated through the obstacle distance parameters and size parameters to achieve spatial optimization of the main reflection direction; the multipath scattering phase offset is converted into a physically operable phase rotation angle, and phase perturbation is directly performed on the basic reflection signal; finally, a directional interference signal with azimuth specificity is generated by binding the signal and space direction; the refined distribution of terrain reflection energy is achieved, and the dominant interference path is accurately captured; the abstract phase offset is converted into an operable rotation angle parameter, and phase perturbation is applied at the signal level to reproduce the multipath phase distortion characteristics in rainy and foggy environments; the strong correlation binding of signal and space direction ensures the accurate projection of interference energy in physical direction, meeting the realistic simulation requirements of strong interference scenarios in the heading channel.

[0055] These and other aspects of the present invention will become more readily apparent from the following description of the embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following is a brief introduction to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0057] Figure 1 A flow chart showing a method for simulating an aviation radio signal jamming system provided by the present invention is shown;

[0058] Figure 2 A schematic structural diagram of a simulation system for an aviation radio signal jamming system provided by the present invention is shown;

[0059] Figure 3 A schematic structural diagram of a computing device provided by the present invention is shown. DETAILED DESCRIPTION

[0060] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0061] In some of the processes described in the specification and claims of the present invention and the above-mentioned figures, multiple operations that appear in a specific order are included, but it should be clearly understood that these operations may not be executed in the order in which they appear in this article or may be executed in parallel. The serial numbers of the operations, such as 101, 102, etc., are only used to distinguish between different operations, and the serial numbers themselves do not represent any execution order. In addition, these processes may include more or fewer operations, and these operations may be executed in sequence or in parallel. It should be noted that the descriptions of "first", "second", etc. in this article are used to distinguish different messages, devices, modules, etc., and do not represent the order of precedence, nor do they limit "first" and "second" to be different types.

[0062] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making any creative efforts shall fall within the scope of protection of the present invention.

[0063] In response to the high-precision assessment requirements for navigation system stability in the landing guidance phase of aviation radio signal interference simulation, existing technologies mostly use simulation methods based on the fusion of multipath channel modeling and terrain databases, but they have obvious deficiencies in handling the effects of complex meteorological conditions and terrain coupling. Specifically, traditional solutions mainly rely on static terrain data to extract reflection path information, which makes it difficult to reflect the dynamic impact of meteorological factors such as rain and fog on the multipath scattering characteristics of the signal, resulting in the generated interference signal deviating from the actual environment in phase offset and spatial distribution. At the same time, there is a lack of refined modeling in the identification of the main reflection direction and the distribution of directional interference energy, and it is impossible to accurately reproduce the strong interference effect in a specific direction. In order to solve the above problems, the present invention conducts aviation radio signal interference simulation by comprehensively considering the terrain geometry, meteorological attenuation factors and multipath phase offset. Specifically, based on the acquisition of the target navigation station's transmission parameters, three-dimensional terrain coordinates, and rain and fog meteorological parameters, a terrain reflection azimuth vector set is constructed and the main reflection direction is extracted. The directional reflection interference signal is generated by combining the phase modulation caused by the weather. Finally, the meteorological attenuation factor and the original transmission parameters are integrated to synthesize the course deviation simulation signal, thereby achieving a high-fidelity restoration of the signal distortion characteristics in a complex electromagnetic environment and improving the accuracy and practicality of the interference simulation system. Figure 1 A flowchart of a method for simulating an aviation radio signal jamming system is provided for an embodiment of the present invention. Figure 1 As shown, the method includes:

[0064] Step 101: Acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of the surrounding terrain obstacles, and the rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include the rain intensity parameter and the fog droplet concentration parameter;

[0065] In this step, the target navigation station refers to the localizer in the instrument landing system, including its geographical coordinates, transmission frequency, signal modulation mode and other attributes, which are used to generate the localizer guidance signal. The transmission signal parameters are the set of physical characteristic parameters of the radio signal transmitted by the guidance station, including the carrier frequency value, initial amplitude value, modulation type code, which are used to reconstruct the reference navigation signal. The surrounding terrain obstacles refer to the natural or artificial features around the airport that affect the propagation of radio signals, including mountains, buildings, etc., and their three-dimensional spatial coordinate set contains longitude, latitude, and altitude data. The rain and fog meteorological parameters refer to the set of atmospheric environment data monitored in real time, which are used to quantify the impact of precipitation and fog on radio wave propagation, including the instantaneous value sequence of rainfall intensity and the spatial concentration distribution matrix of droplets. The rainfall intensity parameter refers to the dynamic monitoring value sequence of rainfall per unit time, which reflects the changing characteristics of rainfall intensity, is measured in millimeters per hour, and is used to calculate the degree of signal attenuation. The droplet concentration parameter refers to the distribution density value of droplet particles in unit volume of air, which is measured in grams per cubic meter, reflects the spatial heterogeneity of fog, and is used to calculate the multipath scattering phase offset.

[0066] In this embodiment of the present invention, the transmission signal parameters of the target navigation station are collected in real time through the aviation simulator data interface. The three-dimensional spatial coordinates of surrounding terrain obstacles are simultaneously loaded from a geographic information system. Rain and fog meteorological parameters are also acquired synchronously through a meteorological sensor network. The transmission signal parameters drive the signal generation module, the three-dimensional spatial coordinates are input into the terrain analysis module, and the rain and fog meteorological parameters are transmitted to the meteorological processing module, achieving the coordinated collection and distribution of multi-source data.

[0067] Step 102: Calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter;

[0068] In this embodiment of the present invention, a sequence of instantaneous rainfall intensity values ​​is generated based on rainfall intensity parameters. The difference between the maximum and minimum values ​​in the sequence is calculated to obtain the fluctuation amplitude. The fluctuation amplitude is divided by a preset rainfall intensity reference value to obtain the rainfall intensity persistence coefficient. This coefficient is multiplied by a preset atmospheric attenuation constant to generate a meteorological attenuation factor. Simultaneously, the spatial distribution characteristics of the droplet concentration parameters are analyzed to identify the three-dimensional location of the highest concentration point, and the straight-line distance from this location to the navigation station is calculated. The inverse of the straight-line distance is multiplied by a preset phase perturbation base value to generate the multipath scattering phase offset.

[0069] Step 103: Analyzing the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional space coordinate set to generate a terrain reflection azimuth vector set;

[0070] In an embodiment of the present invention, each obstacle coordinate in the three-dimensional space coordinate set is traversed: the difference between the obstacle coordinate and the navigation station coordinate in the east-west direction (i.e., the first direction difference in the horizontal plane), the difference in the north-south direction (i.e., the second direction difference in the horizontal plane), and the height difference (i.e., the vertical direction difference) are calculated; the first direction difference is divided by the second direction difference to obtain a proportional relationship, and the horizontal azimuth is determined based on this; the square root of the sum of the squares of the first direction difference and the second direction difference is calculated to generate a horizontal projection length; the vertical direction difference is divided by the horizontal projection length to obtain a proportional relationship, and the vertical elevation angle is determined based on this; the horizontal azimuth and the vertical elevation angle are combined to form an azimuth elevation angle vector, and finally all the vectors are collected to form a terrain reflection azimuth vector set.

[0071] Step 104: extracting a main reflection direction from the terrain reflection azimuth vector set, and combining it with the multipath scattering phase offset modulation to generate a directional reflection interference signal;

[0072] In this embodiment of the present invention, obstacle distance and size parameters corresponding to each azimuth elevation vector in the terrain reflection azimuth vector set are calculated. The inverse of the distance parameter is used as the distance attenuation factor, and the size parameter is multiplied by the unity gain coefficient as the reflection gain factor. The two factors are multiplied together to generate the corresponding directional intensity weight value. The vector with the largest weight value is selected as the main reflection direction, and the multipath scattering phase offset is converted into a phase rotation angle. A basic reflection signal is generated based on the transmission signal parameters. A phase rotation operation is performed on the basic reflection signal based on the phase rotation angle to generate a phase-perturbed reflection signal. Finally, this signal is bound to the spatial pointing coordinates of the main reflection direction to generate a directional reflection interference signal.

[0073] Step 105: generating a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal, and the transmission signal parameter;

[0074] In an embodiment of the present invention, a transmission signal parameter is converted into a sinusoidal carrier signal and used as a reference navigation signal; a directional reflection interference signal and a reference navigation signal are waveform-superimposed to generate a synthetic interference signal; a meteorological attenuation factor is converted into an attenuation proportional coefficient, and the coefficient is multiplied by the amplitude component of the synthetic interference signal to generate an attenuated interference signal; a phase offset is calculated based on the azimuth angle between the main reflection direction and the standard localizer reference azimuth, and the attenuated interference signal is phase-shifted and encoded based on the phase offset to generate a localizer deviation simulation signal.

[0075] The embodiments of the present invention achieve accurate quantification of meteorological attenuation factors and multipath scattering phase offsets; break through the angular resolution limitations of traditional digital elevation models; solve the problems of main reflection direction identification deviation and interference signal phase misalignment; and improve the fidelity and reliability of aviation radio interference simulation.

[0076] The present invention provides a specific embodiment, step 102, calculating the meteorological attenuation factor and the multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter, specifically includes the following steps:

[0077] Step 201: generating an instantaneous rainfall intensity value sequence based on the rainfall intensity parameter, and calculating the fluctuation amplitude of the instantaneous rainfall intensity value sequence;

[0078] In this step, the instantaneous rainfall intensity sequence refers to a collection of rainfall intensity data sampled in chronological order. This data reflects the dynamic characteristics of rainfall intensity over time. High-frequency sampling (e.g., once per second) is used to obtain continuous instantaneous rainfall intensity values, providing a data foundation for analyzing rainfall intensity fluctuations. The fluctuation amplitude, the difference between the maximum and minimum values ​​in the instantaneous rainfall intensity sequence, measures the magnitude of rainfall intensity variation over a given period. A larger value indicates more dramatic rainfall intensity fluctuations and a more significant dynamic impact on signal attenuation.

[0079] In an embodiment of the present invention, based on the rainfall intensity parameter, a sequence of instantaneous rainfall intensity values ​​with equal time intervals is generated through time series sampling technology. For example, sampling is performed once per second to form 100 continuous rainfall intensity data points. The fluctuation amplitude is obtained by calculating the difference between the maximum and minimum values ​​in the sequence to characterize the severity of the change in rainfall intensity in the time dimension.

[0080] Step 202: Determine a rainfall intensity persistence coefficient based on a proportional relationship between the fluctuation amplitude and a preset rainfall intensity benchmark, and generate a meteorological attenuation factor in combination with a preset atmospheric attenuation constant;

[0081] In this step, the preset rainfall intensity benchmark refers to a pre-set rainfall intensity reference threshold (e.g., 5 mm / h corresponds to light rain), which is used in proportional calculations with the fluctuation amplitude to quantify the degree of persistence of rainfall intensity changes on signal attenuation. This benchmark value can be adjusted according to different meteorological standards or application scenarios. The rainfall intensity persistence coefficient is a coefficient determined by the proportional relationship between the fluctuation amplitude and the preset rainfall intensity benchmark. It reflects the weight of the impact of the persistence of rainfall intensity changes on signal attenuation. This coefficient, combined with the atmospheric attenuation constant, generates the final meteorological attenuation factor. The preset atmospheric attenuation constant is a constant pre-set based on electromagnetic wave propagation theory (e.g., 0.01 dB / (km·mm / h)). It quantifies the degree of signal attenuation per unit rainfall intensity per unit distance and is one of the basic parameters for calculating the meteorological attenuation factor. The meteorological attenuation factor is a parameter that integrates the characteristics of rainfall intensity fluctuation and atmospheric attenuation. It is used to characterize the degree of attenuation of radio signals in rainy and foggy weather. This factor directly influences the amplitude modulation process of the subsequent interference signal.

[0082] In an embodiment of the present invention, a ratio operation is performed on the fluctuation amplitude and a preset rainfall intensity reference (such as the rainfall intensity threshold corresponding to moderate rain) to obtain a coefficient reflecting the intensity of the continuous effect of rainfall intensity. For example, when the fluctuation amplitude is 1.5 times the reference value, the coefficient is 1.5. This coefficient is then multiplied by a preset atmospheric attenuation constant to generate a meteorological attenuation factor for quantifying the degree of signal attenuation caused by rain and fog.

[0083] Step 203: analyzing the spatial distribution characteristics of the droplet concentration parameters to identify the three-dimensional position of the droplet concentration peak area;

[0084] In this step, spatial distribution characteristics refer to the distribution of droplet concentration in three-dimensional space. These characteristics, described by parameters such as concentration gradient and peak position, are used to identify the region with the greatest impact on signal scattering (i.e., the peak droplet concentration region). The peak droplet concentration region is the three-dimensional coordinate region where the droplet concentration reaches its maximum value. Droplets within this region have the strongest scattering effect on radio signals and serve as a key reference for calculating multipath scattering phase offset.

[0085] In an embodiment of the present invention, the three-dimensional grid distribution data of the droplet concentration parameters is analyzed, and a neighborhood extreme value detection algorithm is used to traverse all grid cells. When the concentration value of a cell is greater than that of eight adjacent cells, it is marked as a peak point. All peak points are gathered to form a droplet concentration peak area, and the three-dimensional position coordinates of its center point are output.

[0086] Step 204: Calculate the droplet spatial impact weight based on the straight-line distance between the three-dimensional position and the target navigation station, and generate a multipath scattering phase offset in combination with a preset phase perturbation base value;

[0087] In this step, the droplet spatial impact weight refers to the weight coefficient calculated based on the straight-line distance between the droplet concentration peak area and the target navigation station. The closer the distance, the greater the weight. It is used to quantify the degree of influence of spatial position on the signal scattering phase offset. The preset phase perturbation base value refers to the pre-set phase offset base value. When combined with the droplet spatial impact weight, the final multipath scattering phase offset is generated, reflecting the basic amplitude of the signal phase change caused by droplet scattering. The multipath scattering phase offset refers to the phase parameter that integrates the spatial distribution of droplets and distance factors. It is used to characterize the phase change of the multipath propagation of radio signals caused by droplet scattering. It is a key parameter for generating phase modulation of directional reflection interference signals.

[0088] In an embodiment of the present invention, the Euclidean distance from the center point of the droplet concentration peak area to the target navigation station is calculated as the straight-line distance; the droplet spatial influence weight is calculated using an inverse proportional function of the distance (such as 1 / distance) (the closer the distance, the greater the weight), and the weight is multiplied by a preset phase perturbation base value (such as π / 2 radians) to generate a multipath scattering phase offset.

[0089] The embodiment of the present invention converts the temporal characteristics of rainfall intensity changes into attenuation factors, making the meteorological attenuation simulation more consistent with the dynamic process of actual rainfall and avoiding the errors of traditional static rain attenuation models; the droplet scattering phase is refined to achieve spatially differentiated simulation of multipath scattering phase characteristics, improving the matching degree between the phase characteristics of the interference signal and the actual environment; the rain intensity attenuation and the droplet phase offset are decoupled and calculated, and then work synergistically, so that the generated interference signal simultaneously reflects the combined effects of amplitude attenuation and phase disturbance, thereby improving the realism and accuracy of the simulation of aviation radio signal interference.

[0090] For example, in the interference simulation scenario of a target navigation station at an airport (coordinate (0,0,0)), the real-time rain intensity parameter is obtained as The continuous sampling data is used to generate an instantaneous rainfall intensity value sequence containing 100 data points by sampling every second. The maximum value of the sequence is calculated to be , the minimum value is , the fluctuation amplitude is ; The default rainfall intensity is , the ratio of the two is 0.8, combined with the preset atmospheric attenuation constant , generating the meteorological attenuation factor as At the same time, the spatial distribution data of the droplet concentration parameters are obtained, and the three-dimensional position of the peak droplet concentration area is analyzed to be (1000, 500, 100) m. The straight-line distance between this area and the navigation station is calculated as , the spatial influence weight of the droplet is Combined with the preset phase perturbation base value radians, generating a multipath scattering phase offset of The calculation of meteorological attenuation factor and multipath scattering phase offset is finally completed, providing key parameters for subsequent interference signal generation.

[0091] The present invention provides a specific embodiment, step 103, based on the three-dimensional space coordinate set, analyzing the horizontal azimuth and vertical elevation angles of the surrounding terrain obstacles relative to the target navigation station to generate a terrain reflection azimuth vector set, specifically comprising the following steps:

[0092] Step 301: Calculate a coordinate difference component set between the coordinates of all obstacles in the three-dimensional space coordinate set and the coordinates of the target navigation station, wherein the coordinate difference component set includes a horizontal first direction difference, a horizontal second direction difference, and a vertical direction difference;

[0093] In this step, the coordinate difference component set refers to the coordinate difference data set between the obstacle and the navigation station in the three-dimensional rectangular coordinate system, including the east-west difference, the north-south difference and the elevation difference, which is used for spatial relative position analysis. The first direction difference in the horizontal plane refers to the coordinate difference between the obstacle and the navigation station in the east-west direction, with the east being the positive direction, the unit is meter, and it reflects the position deviation in the longitude direction. The second direction difference in the horizontal plane refers to the coordinate difference between the obstacle and the navigation station in the north-south direction, with the north being the positive direction, the unit is meter, and it reflects the position deviation in the latitude direction. The vertical direction difference refers to the coordinate difference between the obstacle and the navigation station in the elevation direction, with the upward direction being the positive direction, the unit is meter, and it characterizes the difference in terrain height.

[0094] In the embodiment of the present invention, the three-dimensional coordinates (X0, Y0, Z0) of the target navigation station and the coordinates (X i ,Y i ,Z i ), calculate the first direction difference of the horizontal plane X i =X i -X0 (due east), horizontal plane second direction difference Y i =Y i -Y0 (due north), vertical difference Z i =Z i-Z0 (elevation difference), constituting the coordinate difference component set of the i-th obstacle ( X i , Y i , Z i ).

[0095] Step 302: Determine the horizontal azimuth angle of the surrounding terrain obstacle based on the ratio of the first horizontal plane direction difference to the second horizontal plane direction difference;

[0096] In this step, the horizontal azimuth refers to the angle between the horizontal projection line from the navigation station to the obstacle and the north direction, which ranges from 0° to 360° clockwise.

[0097] In the embodiment of the present invention, the first direction difference of the horizontal plane is taken X i Difference in the second direction from the horizontal plane Y i The algebraic ratio of X i / Y i , calculate the azimuth angle value through the four-quadrant inverse tangent function, and map it to the range of 0°-360° to generate the horizontal azimuth angle i , characterizing the planar orientation relationship of the obstacle relative to the navigation station.

[0098] Step 303: Calculate the square root of the sum of the squares of the horizontal plane first direction difference and the horizontal plane second direction difference, and use the square root of the sum of the squares as the horizontal projection length;

[0099] In this step, the horizontal projection length refers to the straight-line distance between the obstacle projection point on the horizontal plane and the navigation station, in meters.

[0100] In the embodiment of the present invention, the first direction difference of the horizontal plane is calculated X i The square value and the difference between the second direction of the horizontal plane Y i The square root of the sum is calculated to get the horizontal projection length L. i , that is, L i = , reflecting the radial distance of the obstacle on the horizontal plane.

[0101] Step 304: determining the vertical elevation angle of the surrounding terrain obstacle based on the proportional relationship between the vertical direction difference and the horizontal projection length;

[0102] In this step, the vertical elevation angle refers to the angle between the line connecting the navigation platform and the obstacle and the horizontal plane, which is positive upwards and is used to quantify the impact of terrain undulation.

[0103] In the embodiment of the present invention, the vertical difference Z i Divide by the horizontal projection length Lᵢ to get the ratio, and calculate the vertical elevation angle using the inverse tangent function ,Right now , the angle range is -90° to +90°, and negative values ​​indicate depression angles.

[0104] Step 305: Combining the horizontal azimuth angle and the vertical elevation angle into an azimuth elevation angle vector, and collecting the azimuth elevation angle vectors of all surrounding terrain obstacles to form a terrain reflection azimuth vector set;

[0105] In this step, the azimuth elevation vector refers to a binary data structure (horizontal azimuth, vertical elevation) that represents the spatial orientation of a single obstacle, and the units are all in degrees. The terrain reflection azimuth vector set refers to an ordered set of all obstacle azimuth elevation vectors, and the data structure is Used to store global terrain reflection path information.

[0106] In the embodiment of the present invention, the horizontal azimuth angle and vertical elevation angle Combined into two-tuples

[0107] Use it as the azimuth elevation vector; traverse all surrounding terrain obstacles to generate N vectors, and sort them according to the numbers of the surrounding terrain obstacles to form a terrain reflection azimuth vector set

[0108] .

[0109] The embodiments of the present invention break through the two-dimensional plane limitations of traditional digital elevation models through three-dimensional decomposition of coordinate difference component sets, accurately capturing the vertical undulation characteristics of the terrain; calculating the vertical elevation angle based on the proportional relationship between the horizontal projection length and the vertical difference, eliminating spherical projection errors, and reproducing the true geometric relationship of multipath reflection in mountainous environments; the azimuth-elevation vector set replaces spatial coordinates with angle values, improving storage efficiency and supporting real-time dynamic terrain loading.

[0110] The present invention provides a specific embodiment, step 104, extracting the main reflection direction from the terrain reflection azimuth vector set, combining the multipath scattering phase offset modulation to generate a directional reflection interference signal, specifically includes the following steps:

[0111] Step 401: Calculate obstacle distance parameters and size parameters corresponding to each azimuth elevation angle vector in the terrain reflection azimuth vector set;

[0112] In this step, the obstacle distance parameter refers to the three-dimensional straight-line distance between the terrain obstacle and the target navigation station. It is obtained through GPS ranging or lidar, and is measured in meters. It reflects the length of the radio wave propagation path. The obstacle size parameter refers to the maximum projected area of ​​the obstacle in the signal propagation direction. It is calculated by projecting the 3D model triangles, and is measured in square meters. It determines the reflection cross-sectional efficiency.

[0113] In an embodiment of the present invention, a straight-line distance between an obstacle and a navigation station is obtained by a ranging sensor as an obstacle distance parameter (unit: meter); and a maximum cross-sectional projection area of ​​the obstacle is calculated by a three-dimensional scanning model as a size parameter (unit: square meter).

[0114] Step 402: Calculate a distance attenuation factor based on the obstacle distance parameter, and simultaneously calculate a reflection gain factor based on the size parameter, and combine the distance attenuation factor to generate a directional strength weight value of each elevation angle vector.

[0115] In this step, the distance attenuation factor is a dimensionless attenuation coefficient generated based on the path loss principle, ranging from 0 to 1. The closer the distance, the larger the factor. The reflection gain factor is the amplification factor of the reflected signal due to the obstacle's size. It is calculated as the ratio of the obstacle's size parameter to the reference area (with an upper limit of 1). It is dimensionless and increases with the size. The directional strength weight is the product of the distance attenuation factor and the reflection gain factor. It comprehensively quantifies the obstacle's reflected energy intensity and is used to optimize the main reflection direction. A larger value indicates a more significant reflection contribution.

[0116] In an embodiment of the present invention, an inverse proportional operation is performed on the obstacle distance parameter to generate a distance attenuation factor; a size parameter is compared with a preset size threshold; when the size parameter is less than the preset size threshold, the size parameter and the unit gain coefficient are multiplied to obtain a reflection gain factor; when the size parameter is greater than or equal to the preset size threshold, a preset fixed gain value is used as the reflection gain factor; the distance attenuation factor and the reflection gain factor are multiplied to generate an original weight coefficient; the obstacle distribution density in the spatial sector corresponding to each azimuth elevation vector is obtained; all obstacle distribution densities and the original weight coefficients are weighted to generate a directional strength weight value for each azimuth elevation vector.

[0117] Step 403: Selecting the azimuth elevation angle vector with the largest directional strength weight value as the main reflection direction, and converting the multipath scattering phase offset into a phase rotation angle;

[0118] In this step, the dominant reflection direction refers to the azimuth-elevation angle vector with the largest directional strength weight. This contains binary data, horizontal azimuth and vertical elevation, representing the spatial orientation of the dominant interference path. The phase rotation angle is the physically operable quantity after unit conversion of the multipath scattering phase offset, expressed in radians, and is used to directly modify the signal phase component.

[0119] In the embodiment of the present invention, the directional intensity weight values ​​of all azimuth and elevation angle vectors are compared, and the vector corresponding to the maximum value is selected as the main reflection direction; at the same time, the multipath scattering phase offset is multiplied by an angle conversion coefficient, such as (π / 180), to generate a phase rotation angle in radians.

[0120] Step 404: Generate a basic reflected signal based on the transmit signal parameters, and perform a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal;

[0121] In this step, the basic reflected signal refers to the unmodulated carrier signal generated according to the transmission parameters, and the expression is , where A is the amplitude, f is the frequency, and t is the time variable. The phase-perturbed reflection signal is the modulated signal with the phase rotation angle superimposed on the basic reflection signal.

[0122] In an embodiment of the present invention, a carrier frequency value and an initial amplitude value are parsed from the transmission signal parameters; a reflection intensity correction coefficient is determined based on the obstacle distance parameter and size parameter corresponding to the main reflection direction; the initial amplitude value is multiplied by the reflection intensity correction coefficient to obtain the reflection signal amplitude value; based on the carrier frequency value and the reflection signal amplitude value, a sinusoidal carrier signal is generated and used as the basic reflection signal; and a phase rotation angle is superimposed on the phase component of the sinusoidal carrier signal to obtain a phase disturbance reflection signal.

[0123] Step 404: Associating and binding the phase disturbance reflection signal with the spatial orientation of the main reflection direction to generate a directional reflection interference signal;

[0124] In this step, the directional reflection interference signal refers to the phase disturbance signal bound to the spatial direction, and the data structure contains the signal waveform and spatial coordinates. , realizing the azimuth constraint of interference energy.

[0125] In the embodiment of the present invention, the horizontal azimuth angle θ and the vertical elevation angle φ of the main reflection direction are extracted to generate the spatial pointing coordinates ; The phase perturbation reflected signal Bind to the spatial pointing coordinates, output data structure {signal waveform: , spatial orientation: } as a directional reflection interference signal.

[0126] In the embodiment of the present invention, the path loss is quantified by the distance attenuation factor, and the reflection gain factor is used to characterize the reflection efficiency. The two are integrated to achieve physical-driven optimization of the dominant reflection path; the abstract phase offset is converted into an operational rotation angle, and phase perturbation is directly applied at the carrier signal level, breaking through the fidelity limitations of traditional statistical models; the binding mechanism of the interference signal and the spatial pointing method ensures the azimuth specificity of the reflection interference in the mountain fly-around scenario, and reduces the simulation error of the course deviation.

[0127] The present invention provides a specific embodiment, step 402, calculating a distance attenuation factor based on the obstacle distance parameter, and calculating a reflection gain factor based on the size parameter, and combining the distance attenuation factor to generate directional strength weight values ​​for each orientation elevation angle vector, specifically comprising the following steps:

[0128] Step 411: performing an inverse proportional operation on the obstacle distance parameter to generate a distance attenuation factor;

[0129] In the embodiment of the present invention, an inverse proportional operation is performed on the obstacle distance parameter D (unit: meter). Specifically, the preset distance reference value (typical value 1000 meters) is divided by the obstacle distance parameter D to generate a distance attenuation factor α = 1000 / D. The value range of this factor is , the closer the distance, the larger the factor.

[0130] Step 412: Compare the size parameter with a preset size threshold. When the size parameter is less than the preset size threshold, multiply the size parameter by a unity gain coefficient to obtain a reflection gain factor. When the size parameter is greater than or equal to the preset size threshold, use a preset fixed gain value as the reflection gain factor.

[0131] In this step, the preset size threshold refers to the critical value for determining the obstacle size level, typically ranging from 1000-2000m 2 , used to switch the gain calculation mode based on the radar cross-section characteristics. The unit gain coefficient is the linear conversion coefficient from the size parameter to the reflection gain factor. The preset fixed gain value is a uniform gain assignment for large obstacles, used to avoid gain overflow caused by oversized obstacles. It is obtained through microwave anechoic chamber calibration.

[0132] In the embodiment of the present invention, the size parameter A (unit: ) with a preset size threshold (typically 1500m 2 ) Comparison: If A < 1500, multiply the size parameter by the unity gain factor (typical value 0.001) to obtain the reflection gain factor =A×0.001; if A≥1500, the preset fixed gain value (typical value 1.5) is directly output as the reflection gain factor.

[0133] Step 413: multiplying the distance attenuation factor by the reflection gain factor to generate an original weight coefficient;

[0134] In this step, the original weight coefficient refers to the single obstacle reflection weight without considering the environment density, reflecting the reflection intensity of the isolated obstacle.

[0135] In the embodiment of the present invention, the distance attenuation factor and reflection gain factor Multiply to get the original weight coefficient .

[0136] Step 414: Obtain the obstacle distribution density in the space sector corresponding to each orientation and elevation angle vector;

[0137] In this step, the spatial sector is a three-dimensional conical area with an apex angle of 10° and an azimuth-elevation vector as its central axis. It defines the range of obstacle density statistics. The sector radius is equal to the maximum obstacle range. The obstacle distribution density refers to the concentration of other obstacles within the spatial sector. It is calculated by dividing the number of coordinate points within the sector by the maximum capacity and has a value range of 0 to 1.

[0138] In this embodiment of the present invention, the spatial sector (conical area with a vertex angle of 10°) pointed by the azimuth elevation vector is used as the statistical unit, and the number of three-dimensional coordinate points N of other obstacles in the sector is calculated. N is divided by the maximum capacity value of the sector (typical value 50) to generate the obstacle density parameter =min(N / 50,1).

[0139] Step 415: performing weighted calculation on all obstacle distribution densities and the original weight coefficients to generate directional strength weight values ​​of all orientation and elevation angle vectors;

[0140] In the embodiment of the present invention, the original weight coefficient γ is combined with the obstacle density parameter Add them together and multiply them by the weight normalization constant (typical value 0.5) to generate the final directional strength weight value , completing the energy weight distribution.

[0141] The embodiments of the present invention address path loss through a distance attenuation factor, while using a size gain factor that overcomes the limitations of a fixed reflection model. Preset size thresholds and fixed gain values ​​prevent weight distortion for large features (such as mountains). Spatial sector density parameters capture the clustering effect of obstacles, improving the accuracy of interference intensity restoration in terrain-obscured areas. Ultimately, the physically realistic distribution of interference energy in complex mountainous and urban scenarios is achieved, reducing the simulation error of the localizer signal offset.

[0142] The present invention provides a specific embodiment, in which step 404 generates a basic reflected signal based on the transmit signal parameters, and performs a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal, specifically comprising the following steps:

[0143] Step 421: parse the carrier frequency value and the initial amplitude value from the transmission signal parameters;

[0144] In this step, the carrier frequency value refers to the number of oscillations per second of the radio waves transmitted by the target navigation station. It is measured in megahertz (MHz) and determines the signal wavelength. It is obtained by parsing the frequency field in the transmission parameter package. The initial amplitude value refers to the raw signal strength value without environmental correction. It is measured in decibel milliwatts (dBm) and reflects the transmitter output power. It is obtained by parsing the amplitude field in the transmission parameter package.

[0145] In this embodiment of the present invention, the transmission signal parameters are read to extract the carrier frequency value (unit: megahertz) representing the number of radio wave oscillations and the initial amplitude value (unit: decibel milliwatt) representing the signal strength benchmark, thereby decoding the basic physical characteristics of the signal.

[0146] Step 422: determining a reflection intensity correction coefficient based on the obstacle distance parameter and size parameter corresponding to the main reflection direction;

[0147] In this step, the reflection intensity correction coefficient refers to the signal strength adjustment factor calculated based on the obstacle characteristics. It is obtained by multiplying the inverse of the distance parameter by the distance constant plus the size parameter divided by the size reference value. It is dimensionless and is used to compensate for terrain reflection loss.

[0148] In an embodiment of the present invention, based on the obstacle distance parameter (unit: meter) and size parameter (unit: square meter) associated with the main reflection direction, the reciprocal of the distance parameter is multiplied by a distance correction constant (typical value 0.8), and the size parameter is divided by a size reference value (typical value 1000) to obtain a size factor. The two are added together to generate a reflection intensity correction coefficient.

[0149] Step 423: performing a product operation on the initial amplitude value and the reflection intensity correction coefficient to generate a reflection signal amplitude value;

[0150] In this step, the reflected signal amplitude value refers to the actual signal strength value after terrain correction. The calculation formula is the initial amplitude value multiplied by the reflection strength correction coefficient, with the unit of decibel milliwatt, which determines the reflected signal energy intensity.

[0151] In the embodiment of the present invention, the reflected signal amplitude value = the initial amplitude value × the reflection intensity correction coefficient, the unit is decibel milliwatt, to complete the terrain adaptability adjustment of the signal intensity.

[0152] Step 424: Generate a sinusoidal carrier signal based on the carrier frequency value and the reflected signal amplitude value, and use the sinusoidal carrier signal as a basic reflected signal;

[0153] In this step, the sinusoidal carrier signal refers to a continuous waveform signal generated by a sine function, and its mathematical expression is A·sin(2πft), where A is the amplitude of the reflected signal, f is the carrier frequency, and t is the time variable. The phase component refers to the angle variable part of the sine function. , whose value change determines the waveform time offset, which is obtained by analyzing the phase term of the sine function.

[0154] In the embodiment of the present invention, a standard sinusoidal wave function is constructed with the carrier frequency value as the oscillation frequency and the reflected signal amplitude value as the peak amplitude: sinusoidal carrier signal (t) = reflected signal amplitude value × sin(2π × carrier frequency value × t). After discretization and sampling of this function, a sinusoidal carrier signal is generated as the basic reflected signal.

[0155] Step 425: superimposing the phase rotation angle onto the phase component of the sinusoidal carrier signal to obtain a phase disturbance reflection signal;

[0156] In the embodiment of the present invention, the phase component of the sinusoidal carrier signal (i.e., the angle variable in the sin function) is extracted, and the phase rotation angle is algebraically added to it to form a new phase component: new phase = original phase component + phase rotation angle; reconstruct the sinusoidal function: phase perturbation reflection signal = amplitude value × sin (2π × frequency value × t + new phase).

[0157] The embodiments of the present invention eliminate the mountain reflection intensity deviation of the fixed amplitude model in the traditional solution by dynamically correcting the signal amplitude based on the distance / size of the obstacle. The rotation angle is directly superimposed on the phase component of the carrier signal to avoid the phase nonlinear distortion introduced by the traditional filter. The sinusoidal carrier strictly follows the physical laws of radio signal propagation to ensure the electromagnetic authenticity of the course deviation simulation. Ultimately, high-fidelity reproduction of multipath interference signals in complex terrain is achieved, meeting the accuracy requirements of the instrument landing system anti-interference test.

[0158] The present invention provides a specific embodiment, step 105, generating a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal, and the transmission signal parameter, specifically includes the following steps:

[0159] Step 501: Convert the transmission signal parameters into a reference carrier signal, and perform amplitude modulation on the reference carrier signal to generate a reference navigation signal;

[0160] In this step, the reference carrier signal refers to an unmodulated sine wave generated based on the transmission parameters. Its expression is A·sin(2πft), where A is the initial amplitude and f is the carrier frequency. This signal is used for subsequent modulation processing. The reference navigation signal is a compliant amplitude-modulated signal. A 90Hz / 150Hz positioning audio modulation envelope is superimposed on the reference carrier signal to provide standard localizer guidance.

[0161] In this embodiment of the present invention, the carrier frequency and initial amplitude are extracted from the transmitted signal parameters to construct an unmodulated sinusoidal carrier signal as a reference carrier signal. According to the instrument landing system localizer signal specification, 90 Hz and 150 Hz positioning audio frequencies are modulated onto the carrier using amplitude modulation technology to output a reference navigation signal.

[0162] Step 502: performing waveform superposition on the directional reflected interference signal and the reference navigation signal to generate a synthetic interference signal;

[0163] In this step, the synthetic interference signal refers to the waveform superposition result of the reference navigation signal and the directional reflected interference signal, which includes the original signal and the reflected interference component.

[0164] In an embodiment of the present invention, the waveform data of the directional reflected interference signal is aligned with the waveform data of the reference navigation signal along the time axis; an algebraic addition operation is performed on each sampling point: synthetic interference signal (t) = reference navigation signal (t) + directional interference signal (t), thereby generating a synthetic interference signal containing the original signal and the reflected interference.

[0165] Step 503: converting the meteorological attenuation factor into a signal attenuation proportional coefficient, and performing a product operation on the signal attenuation proportional coefficient and the amplitude component of the synthetic interference signal to generate an attenuated interference signal;

[0166] In this step, the signal attenuation coefficient is the linearly converted value of the meteorological attenuation factor. It is dimensionless and ranges from 0 to 1. It is used for amplitude compression. The amplitude component refers to the parameter sequence in the signal waveform that represents the intensity change. It is the envelope data extracted through the Hilbert transform and is measured in decibel milliwatts. The attenuated interference signal is the output of the synthesized interference signal after amplitude compression, reflecting the signal weakening caused by weather.

[0167] In the embodiment of the present invention, the signal attenuation proportional coefficient = the meteorological attenuation factor × k (k is a conversion coefficient); an array of amplitude components of the synthetic interference signal is extracted, each array element is multiplied by the signal attenuation proportional coefficient, and an amplitude-compressed attenuated interference signal is output.

[0168] Step 504: obtaining a standard localizer reference azimuth, and calculating an azimuth deviation angle between the main reflection direction and the standard localizer reference azimuth;

[0169] In this step, the standard localizer reference bearing refers to the runway centerline magnetic heading angle published by the airport, obtained from the aviation database and expressed in degrees, typically ranging from 0° to 360°. The azimuth deviation angle, the angular difference between the horizontal bearing of the main reflection direction and the standard localizer, determines the signal phase offset.

[0170] In an embodiment of the present invention, a standard localizer reference azimuth (e.g., a magnetic heading of 120°) is read from an airport flight program database; the horizontal azimuth of the main reflection direction is extracted; and the difference between the two is calculated: azimuth deviation angle = |main reflection direction azimuth-standard localizer reference azimuth|.

[0171] Step 505: determining a signal phase offset according to the azimuth deviation angle, performing offset encoding processing on the attenuated interference signal, and generating a localizer deviation simulation signal;

[0172] In this step, the signal phase offset refers to the waveform time offset parameter calculated based on the azimuth deviation angle, expressed in radians. The localizer deviation simulation signal is the final output distorted signal, which adds a phase offset to the attenuated interference signal and is used to drive the flight simulator instrument display.

[0173] In this embodiment of the present invention, the azimuth deviation angle is multiplied by the angular phase conversion coefficient (typically 3 degrees / degree) to generate a signal phase offset. This offset is added to the entire phase component array of the attenuated interference signal, and the sine function waveform is reconstructed to output the final course deviation analog signal.

[0174] The embodiments of the present invention generate localizer signals by strictly adhering to amplitude modulation specifications, ensuring that the simulation basis complies with aviation standards. The time-domain waveform superposition technology retains the time-frequency characteristics of the reflected interference and the original signal, solving the phase distortion of traditional frequency-domain synthesis. Amplitude compression based on the attenuation proportional coefficient is used to reproduce the signal strength attenuation phenomenon caused by heavy rain. The phase offset is driven by the azimuth deviation angle to realize the essential connection between the physical azimuth deviation and signal distortion. Ultimately, it breaks through the simulation bottleneck of existing solutions in complex meteorological and terrain scenarios, and improves the fidelity of the localizer indication deviation.

[0175] Figure 2 The present invention provides a schematic structural diagram of a simulation system for an aviation radio signal jamming system. Figure 2 As shown, the system includes:

[0176] An acquisition module 21 is configured to acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of surrounding terrain obstacles, and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include rain intensity parameters and fog droplet concentration parameters;

[0177] A calculation module 22 is used to calculate the meteorological attenuation factor and the multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter;

[0178] An analysis module 23 is configured to analyze the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set;

[0179] An extraction module 24 is configured to extract a main reflection direction from the terrain reflection azimuth vector set and generate a directional reflection interference signal by combining the multipath scattering phase offset modulation;

[0180] The generating module 25 is configured to generate a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal, and the transmission signal parameter.

[0181] Figure 2 The simulation system of the aviation radio signal jamming system can perform Figure 1 The implementation principles and technical effects of the method for simulating an aviation radio signal jamming system described in the illustrated embodiment will not be elaborated upon. The specific manner in which the various modules and units perform operations in the simulation system for an aviation radio signal jamming system in the aforementioned embodiment have been described in detail in the relevant embodiments of the method and will not be further elaborated upon here.

[0182] In one possible design, Figure 2 The simulation system of an aviation radio signal jamming system of the embodiment shown can be implemented as a computing device, such as Figure 3 As shown, the computing device may include a storage component 31 and a processing component 32;

[0183] The storage component 31 stores one or more computer instructions, wherein the one or more computer instructions are called and executed by the processing component 32 .

[0184] The processing component 32 is used for the above Figure 1 The embodiment provides a method for simulating an aviation radio signal jamming system.

[0185] The processing component 32 may include one or more processors to execute computer instructions to complete all or part of the steps in the above method. Of course, the processing component may also be implemented as one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the above method.

[0186] The storage component 31 is configured to store various types of data to support operations at the terminal. The storage component can be implemented by any type of volatile or non-volatile memory device, or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk, or optical disk.

[0187] Of course, a computing device may also include other components, such as input / output interfaces, display components, communication components, etc.

[0188] The input / output interface provides an interface between the processing component and the peripheral interface module, which can be an output device, an input device, etc.

[0189] The communication component is configured to facilitate, among other things, wired or wireless communications between the computing device and other devices.

[0190] Among them, the computing device can be a physical device or an elastic computing host provided by a cloud computing platform, etc. In this case, the computing device can refer to a cloud server, and the above-mentioned processing components, storage components, etc. can be basic server resources rented or purchased from the cloud computing platform.

[0191] The embodiment of the present invention further provides a computer storage medium storing a computer program, which can achieve the above-mentioned Figure 1 A simulation method for an aviation radio signal jamming system according to the illustrated embodiment.

[0192] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0193] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they may be located in one location or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of the present embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0194] Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented using software plus a necessary general-purpose hardware platform, or of course, hardware. Based on this understanding, the essence of the above technical solution, or the portion that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, or an optical disk, and includes a number of instructions for causing a computer device (such as a personal computer, server, or network device) to execute the methods described in each embodiment or certain portions of the embodiments.

[0195] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for simulating an aviation radio signal jamming system, characterized in that: include: Acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of surrounding terrain obstacles, and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include rain intensity parameters and fog droplet concentration parameters; Calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter; Analyzing the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set; Extracting the main reflection direction from the terrain reflection azimuth vector set, and combining it with the multipath scattering phase offset modulation to generate a directional reflection interference signal; Based on the meteorological attenuation factor, the directional reflection interference signal and the transmission signal parameter, a localizer deviation simulation signal for an aviation radio signal jamming system is generated.

2. The method according to claim 1, characterized in that Calculating a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter includes: Based on the rainfall intensity parameter, generating an instantaneous rainfall intensity value sequence, and calculating the fluctuation amplitude of the instantaneous rainfall intensity value sequence; Based on the proportional relationship between the fluctuation amplitude and the preset rainfall intensity benchmark, a rainfall intensity persistence coefficient is determined, and combined with a preset atmospheric attenuation constant, a meteorological attenuation factor is generated; Analyzing the spatial distribution characteristics of the droplet concentration parameters to identify the three-dimensional position of the droplet concentration peak area; The spatial impact weight of the droplet is calculated according to the straight-line distance between the three-dimensional position and the target navigation station, and the multipath scattering phase offset is generated in combination with a preset phase disturbance base value.

3. The method according to claim 1, characterized in that According to the three-dimensional space coordinate set, analyzing the horizontal azimuth and vertical elevation angles of surrounding terrain obstacles relative to the target navigation station to generate a terrain reflection azimuth vector set, including: Calculating a coordinate difference component set between the coordinates of all obstacles and the coordinates of the target navigation station in the three-dimensional space coordinate set, wherein the coordinate difference component set includes a horizontal plane first direction difference, a horizontal plane second direction difference, and a vertical direction difference; determining the horizontal azimuth of the surrounding terrain obstacle according to the ratio of the first horizontal plane direction difference to the second horizontal plane direction difference; Calculating the square root of the sum of the squares of the first horizontal plane direction difference and the second horizontal plane direction difference, and using the square root of the sum of the squares as the horizontal projection length; determining the vertical elevation angle of the surrounding terrain obstacle based on a proportional relationship between the vertical direction difference and the horizontal projection length; The horizontal azimuth angle and the vertical elevation angle are combined into an azimuth elevation angle vector, and the azimuth elevation angle vectors of all surrounding terrain obstacles are collected to form a terrain reflection azimuth vector set.

4. The method according to claim 1, wherein Extracting a main reflection direction from the terrain reflection azimuth vector set, and combining it with the multipath scattering phase offset modulation to generate a directional reflection interference signal, including: Calculating obstacle distance parameters and size parameters corresponding to each azimuth elevation angle vector in the terrain reflection azimuth vector set; Calculating a distance attenuation factor based on the obstacle distance parameter and a reflection gain factor based on the size parameter, and generating a directional strength weight value of each elevation angle vector in combination with the distance attenuation factor; Selecting the azimuth elevation angle vector with the largest directional intensity weight value as the main reflection direction, and converting the multipath scattering phase offset into a phase rotation angle; generating a basic reflected signal based on the transmit signal parameter, and performing a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal; The phase disturbance reflection signal is associated and bound with the spatial orientation of the main reflection direction to generate a directional reflection interference signal.

5. The method according to claim 4, characterized in that Calculating a distance attenuation factor based on the obstacle distance parameter and a reflection gain factor based on the size parameter, and combining the distance attenuation factor to generate a directional strength weight value of each elevation angle vector, including: Performing an inverse proportional operation on the obstacle distance parameter to generate a distance attenuation factor; Comparing the size parameter with a preset size threshold, when the size parameter is less than the preset size threshold, multiplying the size parameter by a unity gain coefficient to obtain a reflection gain factor, and when the size parameter is greater than or equal to the preset size threshold, using a preset fixed gain value as the reflection gain factor; Multiplying the distance attenuation factor by the reflection gain factor to generate an original weight coefficient; Obtain the obstacle distribution density in the space sector corresponding to each azimuth elevation vector; All obstacle distribution densities and the original weight coefficients are weighted and calculated to generate directional strength weight values ​​of all orientation elevation vectors.

6. The method according to claim 4, characterized in that Generating a basic reflected signal based on the transmit signal parameter, and performing a phase rotation operation on the basic reflected signal according to the phase rotation angle to generate a phase-perturbed reflected signal, comprising: Resolving the carrier frequency value and the initial amplitude value from the transmission signal parameters; Determining a reflection intensity correction coefficient based on an obstacle distance parameter and a size parameter corresponding to the main reflection direction; Performing a product operation on the initial amplitude value and the reflection intensity correction coefficient to generate a reflection signal amplitude value; generating a sinusoidal carrier signal based on the carrier frequency value and the reflected signal amplitude value, and using the sinusoidal carrier signal as a basic reflected signal; The phase rotation angle is superimposed on the phase component of the sinusoidal carrier signal to obtain a phase disturbance reflection signal.

7. The method according to claim 1, characterized in that Generating a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal, and the transmission signal parameter, including: Converting the transmission signal parameters into a reference carrier signal, and performing amplitude modulation on the reference carrier signal to generate a reference navigation signal; Performing waveform superposition of the directional reflection interference signal and the reference navigation signal to generate a synthetic interference signal; Converting the meteorological attenuation factor into a signal attenuation proportional coefficient, and performing a product operation on the signal attenuation proportional coefficient and the amplitude component of the synthetic interference signal to generate an attenuated interference signal; Obtaining a standard localizer reference azimuth, and calculating an azimuth deviation angle between the main reflection direction and the standard localizer reference azimuth; A signal phase offset is determined according to the azimuth deviation angle, so as to perform offset coding processing on the attenuated interference signal and generate a course deviation simulation signal.

8. A simulation system for an aviation radio signal jamming system, characterized in that: include: An acquisition module is used to acquire the transmission signal parameters of the target navigation station, the three-dimensional spatial coordinate set of the surrounding terrain obstacles, and rain and fog meteorological parameters, wherein the rain and fog meteorological parameters include rain intensity parameters and fog droplet concentration parameters; A calculation module, configured to calculate a meteorological attenuation factor and a multipath scattering phase offset based on the rain intensity parameter and the fog droplet concentration parameter; an analysis module, configured to analyze the horizontal azimuth and vertical elevation of surrounding terrain obstacles relative to the target navigation station based on the three-dimensional spatial coordinate set to generate a terrain reflection azimuth vector set; An extraction module, configured to extract a main reflection direction from the terrain reflection azimuth vector set, and generate a directional reflection interference signal by combining the multipath scattering phase offset modulation; A generating module is used to generate a localizer deviation simulation signal for an aviation radio signal jamming system based on the meteorological attenuation factor, the directional reflection jamming signal and the transmission signal parameter.

9. A computing device, characterized in that It comprises a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component to implement a simulation method for an aviation radio signal interference system as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that A computer program is stored, and when the computer program is executed by a computer, a simulation method of an aviation radio signal jamming system as described in any one of claims 1 to 7 is implemented.

Images