196 results about "Ionosphere" patented technology

The invention discloses a global ionosphere modeling method and device based on Beidou single-frequency ionosphere weighted PPP, and belongs to the technical field of GNSS navigation positioning and atmosphere modeling. The method comprises the following steps: acquiring real-time precise orbit and clock error information of a Beidou satellite and an initial ionospheric delay correction value; a single-frequency ionosphere weighted precision single point positioning (SF-IW-PPP) observation equation is constructed; establishing a space-time correlation constraint model of the ionized layer delay error and determining a prior variance; carrying out combined calculation by taking the previous variance as a constraint to obtain an ionosphere delay residual error, and further obtaining the total delay of the high-precision inclined ionosphere; and finally generating a global ionosphere grid model. According to the invention, high-precision and real-time global ionosphere modeling is realized by only needing a low-cost single-frequency receiver and combining Beidou real-time enhancement service, the hardware threshold and data processing delay are remarkably reduced, and the method is suitable for the fields of public navigation, intelligent transportation, space weather monitoring and the like.

The invention discloses an ionospheric scintillation monitoring method based on a geodetic measurement receiver, which belongs to the technical field of satellite navigation, and comprises the following steps of: accurately eliminating non-scintillation errors by constructing an inter-station single-difference function model and an inter-station single-difference random model, extracting carrier phase residual errors by using a single epoch ambiguity fixing technology, and finally obtaining an ionospheric scintillation monitoring result. Establishing a reliable ionospheric scintillation monitoring index; according to the ionospheric scintillation monitoring method based on the geodetic measurement type receiver, high-precision real-time monitoring of ionospheric scintillation can be realized through a single-epoch ambiguity fixing technology under the condition of not depending on an external precision product, and reliable technical support is provided for space weather monitoring.

The invention discloses a positioning method, device and equipment based on ionosphere error correction and a medium, and the method comprises the steps: obtaining a real-time observation data set of a target CORS system, calculating the real-time ionosphere delay of each observation satellite in each reference station, and obtaining a real-time ionosphere delay set; calculating longitude and latitude correction parameters of each reference station through the real-time ionosphere delay set and longitude and latitude change rate reference coefficients pre-acquired by each reference station, and updating a pre-acquired original ionosphere correction model through the longitude and latitude correction parameters to obtain a new ionosphere correction model; and generating a corrected ionosphere error through the new ionosphere correction model, and correcting positioning data acquired by the target CORS system through the corrected ionosphere error. According to the method, the problem that the change rule of the ionosphere on the longitude and latitude is ignored in the traditional positioning technology can be solved, so that the accuracy of the CORS system satellite positioning is improved.

The invention discloses a comprehensive quality evaluation method and system for Beidou deformation monitoring observation data in a complex environment, and relates to the technical field of satellite navigation positioning, and the method comprises the steps: carrying out the comprehensive evaluation of an elevating angle, a signal-to-noise ratio, an epoch difference and an environment model of the observation data of a Beidou single station, and determining the initial quality of the data; based on the double-difference observation values among the multiple stations, the data quality evaluation precision is further improved by combining a double-difference ambiguity resolving method; according to the method, ionosphere delay, three-difference errors and inter-epoch difference errors are used for data fusion, a more accurate quality evaluation result is obtained, and finally, the data are optimized through comprehensive analysis based on the quality evaluation results of various types of data so as to ensure accurate positioning and deformation monitoring. The technology aims at improving the accuracy and reliability of data quality evaluation by comprehensively analyzing influence factors of Beidou observation data in a complex environment, and provides support for geological disaster monitoring, smart city construction and other high-precision application scenes.

The invention relates to the technical field of engineering monitoring, in particular to an unmanned aerial vehicle dynamic monitoring and early warning system and method for the whole subway tunnel excavation process, and the system comprises a virtual base station differential positioning module which is used for carrying out the unmanned aerial vehicle dynamic monitoring and early warning according to a carrier phase observation value transmitted by a ground active control target and a satellite signal observation value obtained by an airborne receiving end, and calculating ionosphere residual values and troposphere residual values between the current position of the unmanned aerial vehicle and the three adjacent ground active control targets. According to the method, the pixel gradient vector of the thermal inertia matrix is further calculated, the edge candidate points exceeding the threshold value are extracted, and the edge candidate points are closed in the gradient vertical direction to form the geometric region, so that the fuzzy thermal field anomaly is converted into the quantized graph with the definite physical boundary, and the hidden danger area and the internal gradient mean value are combined for calculation; the risk early warning level can be dynamically generated according to the tunnel design tolerance, so that recognition and hierarchical management and control of hidden geological disasters are realized on the premise of not contacting the ground.

The invention provides a frequency point GNSS observation value determination method and device, and a storage medium. The method comprises the following steps: determining a satellite-ground distance difference value between a first frequency point of a known observation value and a second frequency point of an unknown observation value; according to the satellite-to-earth distance difference value between the first frequency point and the second frequency point, the observation value of the first frequency point, the state space representation correction number of the satellite and the center frequencies of the first frequency point and the second frequency point, the pseudo-range observation value of the second frequency point and / or the phase observation value of the second frequency point are / is determined, on one hand, the satellite-to-earth distance difference between different frequency points is fully considered; on the one hand, the precision of pseudo-range observation and / or phase observation values of the recovered frequency points is higher, the requirements for the first frequency point and the second frequency point are reduced, the universality of the method is improved, on the other hand, high-precision ionized layer delay in the SSR correction number and pseudo-range hardware delay and phase hardware delay of a satellite end are effectively utilized for error compensation, and the error compensation accuracy is improved. The implementation difficulty is reduced, and the feasibility of the method is improved.

Disclosed in the present invention are a method and apparatus for constructing an ionospheric mapping function model, a device, and a storage medium. The method comprises: acquiring historical slant total electron content data of an ionosphere; in the historical slant total electron content data, taking, as a corresponding historical vertical total electron content value, a historical slant total electron content value corresponding to a satellite elevation angle greater than a preset angle, and calculating a corresponding historical mapping function value; acquiring a corresponding Modified Julian Day, local time, satellite elevation angle and azimuth angle of an ionospheric pierce point; and on the basis of the historical mapping function value, the Modified Julian Day, the local time, the satellite elevation angle, and the azimuth angle of the ionospheric pierce point, training a preset LSTM model by taking the Modified Julian Day, the local time, the satellite elevation angle and the azimuth angle of the ionospheric pierce point as inputs and taking the historical mapping function value as an output, so as to obtain an ionospheric mapping function model. The present invention can improve the accuracy of ionospheric mapping functions.

The invention discloses a wide-lane ambiguity fixing method and device, a computer storage medium and a terminal, according to the embodiment of the invention, a high-precision double-difference ionized layer estimation value is obtained through joint calculation based on a multi-baseline ionized layer adjustment value and an ionized layer interpolation value, and the defect that a wide-lane observation equation analysis method is not easy to operate under the conditions of an ionized layer active area and a long baseline is effectively overcome. The delay of the ionized layer is obviously increased; when the difference between the ionosphere adjustment value and the ionosphere interpolation value is smaller than or equal to a preset threshold value, a pseudo observation equation established according to the ionosphere estimation value is fused into a Kalman filtering equation, measurement updating is carried out on a wide-lane ambiguity floating point solution, and under the condition that the efficiency advantage of direct resolving of a wide-lane observation equation analytical method is reserved, the wide-lane ambiguity floating point solution is obtained. The ambiguity floating point solution is constrained through the pseudo observation equation, the influence of the ionosphere on ambiguity fixation is effectively inhibited, and the reliability and resolving efficiency of ambiguity fixation are improved.

The invention relates to the technical field of satellite navigation meteorological remote sensing, and discloses a water vapor ionized layer snow depth integrated detection method and system. The method comprises the steps that direct and reflected signals of a Beidou satellite are received through fusion antenna array differentiation, carrier phase and pseudo-range observation values are obtained through carrier tracking loop synchronous demodulation, the total atmospheric rainfall amount and the total ionosphere electron content are extracted through a troposphere and ionosphere combined calculation algorithm, the snow depth is inverted through reflection path geometric difference, and the snow depth is calculated. And carrying out space-time synchronization constraint processing on the extended Kalman filter to realize multi-parameter integrated detection. The problems of low efficiency, insufficient precision and lack of cooperative processing caused by independent observation of water vapor, ionized layer and accumulated snow parameters in the prior art are solved, and the observation efficiency and the measurement precision of multi-parameter integrated detection are improved.

Described herein are solutions for real-time precise ionosphere corrections. An integrate and dump (I&D) server can receive a Networked Transport of Radio Technical Commission For Maritime Services (RTCM) via Internet Protocol (NTRIP) bit stream of a spherical harmonic model for precise point positioning (PPP). The model can be a vertical total electron content (VTEC) spherical harmonic expansion model and / or can include a set of ionosphere coefficients for correcting errors in satellite signals due to ionospheric conditions. The I&D server can generate a file of ionosphere coefficients, according to a refresh timer, and can send the file to a content delivery network (CDN) for distribution to user equipment (UEs) upon request. These and many other features and examples are described herein.

The invention discloses a method for constructing a regional polynomial ionosphere oblique delay model, which belongs to the field of ionosphere delay error correction of a global satellite navigation system, and comprises the following steps of: calculating a signal puncture point coordinate through a reference station latitude and longitude coordinate, and calculating a signal puncture point coordinate by taking the signal puncture point coordinate, observation time and a satellite elevation angle as joint independent variables; sTEC constructs a polynomial function model for dependent variables; the spatial-temporal dynamic characteristics and the geometric characteristics of the sight line of the satellite are synchronously incorporated into modeling, and the spatial-temporal change rule of the ionosphere region is accurately captured; by utilizing polynomial analysis characteristics, model parameters can be quickly output, a user can directly calculate the required STEC only by inputting own coordinates, time, a satellite elevation angle and a satellite azimuth angle, the process is remarkably simplified, and the real-time performance is improved. According to the invention, the defects of precision and efficiency in a complex regional scene in the prior art are overcome, and a more reliable ionosphere error correction scheme is provided for regional GNSS high-precision positioning.

The invention relates to a Beidou No.3 precise point positioning method based on double-frequency combination and analysis center product optimization, and belongs to the technical field of precise point positioning, and the method comprises the following steps: S1, constructing three ionosphere-free combinations B1C-B2a, B1C-B3I and B1I-B3I by using four frequency point signals B1C, B2a, B1I and B3I; s2, through comparing a root mean square value RMS of the positioning deviation, the final positioning precision and the convergence time, respectively comparing and evaluating precision point positioning PPP performance of BDS-3 in different frequency combinations under precision products issued by a plurality of analysis centers; and S3, finding out the ionosphere-free combination and analysis center product with the optimal PPP calculation performance.

The invention discloses a Beidou / GNSS differential service fault analysis method and system, and relates to satellite positioning. The method comprises the following steps: carrying out pseudo-range single-point positioning on Beidou / GNSS observation data to obtain a general position; reporting the approximate position to a Beidou / GNSS differential service platform to request Beidou / GNSS differential data of the current position; carrying out observation data quality analysis on the Beidou / GNSS observation data, wherein observation data quality analysis indexes comprise an integrity rate, a satellite number, a carrier-to-noise ratio, multiple paths, pseudo-range noise, carrier noise, a cycle slip ratio and a PDOP value; differential data quality analysis is carried out on the Beidou / GNSS differential data; rTK calculation is carried out on the Beidou / GNSS observation data and the Beidou / GNSS differential data; calculating a troposphere residual error and an ionosphere residual error; and performing fault detection according to the data. Aiming at the problem that Beidou / GNSS differential service fault positioning precision is low in the prior art, accurate identification and positioning of differential service faults are realized.

The invention discloses an ionosphere real-time monitoring method, device and system, and a storage medium, and the method comprises the steps: configuring a multi-frequency GNSS receiver on a low-orbit satellite, obtaining a navigation satellite signal in real time, rapidly calculating the total electron content of an oblique ionosphere through the dual-frequency carrier phase difference, and calculating the observation confidence based on a Kalman filtering model; further, the satellite end carries out anomaly monitoring through a dynamic threshold value, and respectively generates and downloads a compressed observation flow containing a quality index and original carrier phase observation data fragments in a time window before and after the anomaly is triggered; and the ground station performs weighted fusion on the compressed observation flow according to the confidence coefficient to generate high-precision TEC data. According to the technical scheme, high-precision TEC data products for ionosphere research and abnormal original data products for deep diagnosis are finally produced through on-satellite intelligent processing and selective downloading, communication bandwidth occupation can be remarkably reduced, and observation real-time performance and data reliability are improved.

The invention provides a ground GNSS differential data correction method and device, and the method comprises the steps: obtaining a pseudo range and a carrier phase, calculating a whole cycle jump value of the carrier phase through employing an M-W combination function and an ionosphere delay change rate, and correcting a carrier phase value through employing the calculated whole cycle jump value of the carrier phase. The influence of station-satellite geometric distance is eliminated, the influence of ionosphere change is considered, and the method is suitable for cycle slip detection in a dynamic environment; according to the method, a self-adaptive threshold model is constructed on the basis of Turbo Edit cycle slip detection in a flight dynamic observation environment, cycle slip misjudgment at a low satellite elevation angle is reduced, and the influence of too low sampling rate on GF combined cycle slip detection performance is weakened to a certain extent; and a technical guarantee is provided for applying the GNSS differential data to the flight test.

The application discloses a link time delay calibration result uncertainty modeling evaluation method and device and a medium, and relates to the technical field of satellite navigation high-precision timing. The method comprises the following steps: selecting a time laboratory station and specifying a reference station, constructing a time transfer link, solving the receiver clock error based on an ionosphere-free combination model, and inversely deducing the hardware time delay of a station to be calibrated through clock error difference; the distribution is discriminated through time delay residual statistical characteristics, a model is matched to calculate a type B uncertainty, and an empirical type B uncertainty is obtained through weighted summation; a standard deviation of a calibration time sequence is calculated to obtain a type A standard uncertainty estimation value; a synthetic expanded uncertainty is calculated by fusing the empirical type B uncertainty and the type A standard uncertainty estimation value, the calibration result reliability is evaluated, and uncertainty quantitative analysis of time delay calibration is realized. The application provides a clear and reliable evaluation framework for link calibration error, and significantly improves the reliability and interpretability of the calibration result.

This invention discloses a projection function modeling method and system for geomagnetic soft partitioning and height perturbation correction. The method includes: acquiring puncture point observation samples; converting them to geomagnetic latitude and constructing periodic features; calculating the weights of low, middle, and high latitude sub-models using a soft partitioning mechanism with a fixed center and continuous weights; establishing the effective height function of each sub-model and weighted fusion to obtain the basic effective height; calculating the basic explicit projection function; inputting features into the data to drive the model, outputting the effective height correction, residual, and uncertainty; mapping the correction to an interpretable correction term based on sensitivity and combining it with the residual; and adaptively outputting the complete result or reverting to the basic model based on the uncertainty. This invention combines physical interpretability and data-driven capability, improving the accuracy, regional adaptability, and robustness of the projection function, and can be used for ionospheric delay correction in satellite navigation.

The invention discloses a satellite phase deviation resolving method considering signal distortion deviation, and relates to the technical field of satellite navigation and positioning. The method comprises the following steps: selecting a GNSS base station network to obtain double-frequency observation data; classifying receiver types and unifying same-kind signal distortion delay; signal distortion modeling is introduced based on an MW combined observation equation, and wide-lane satellite phase deviation is obtained through adjustment and ambiguity fixed solution in combination with a precise orbit and clock correction; and based on the ionosphere-free combination observation equation, combining the solved wide-lane deviation, and obtaining the narrow-lane satellite phase deviation and the ionosphere-free combination signal distortion deviation through adjustment and ambiguity fixation. According to whether the receiver type is calibrated or not, the signal distortion deviation is directly applied or estimated in real time, more comprehensive error compensation is achieved in an observation equation, and the positioning precision and the convergence speed are improved. According to the method, synchronous separation and accurate estimation of the wide-lane and narrow-lane phase deviation and the signal distortion deviation of the satellite are realized, and the estimation accuracy and physical consistency of the phase deviation are effectively improved.

The application discloses a method for realizing ionosphere aircraft stealth based on regulating and controlling electron distribution of tail area, first, evolution time scales under different non-extended parameters are calculated; and according to the superhot electron energy and particle flux per unit time required by different non-extended parameters, a radio frequency discharge plasma generator is controlled to generate plasma corresponding to the non-extended parameters; then, in the grid electric field, electrons enter the acceleration field area, and ions are directly discharged to the tail area of the aircraft without acceleration; the regulation and control duration under different non-extended parameters is determined by the evolution time scale corresponding to the non-extended parameter. The application changes the electron distribution of the ionosphere aircraft tail area by randomly injecting high-energy electrons into the tail area of the aircraft, regulates and controls the evolution law and characteristic scale of the non-linear local cavity in the tail area which cannot be eliminated, makes the anti-stealth method of monitoring the motion of the density cavity to detect the ionosphere stealth aircraft invalid, and thus realizes the stealth of the ionosphere stealth aircraft.

The present application relates to the field of geodetic surveying, and particularly relates to a method for determining altitude based on GNSS frequency signals of Chinese space station. Through obtaining data such as GNSS satellite precise orbit, clock error and earth rotation parameters from IGS analysis center, RINEX format data is obtained by using synchronous observation of space station and ground station receiver, and the gravity potential of space station observation moment is determined in combination with EIGEN-6C4 gravity field model. The frequency observation value is derived from the carrier phase observation value, the observation equation is established in the dual-frequency ionosphere elimination processing mode, the gravity potential difference between space station and ground station is solved, and finally the gravity potential and altitude of ground station are obtained. The present application can replace the traditional leveling and gravity measurement, significantly reduce the cost, overcome the environmental restriction, and support the cross-sea and cross-ocean height transfer and the global height datum unification.

This invention relates to a method and system for broadcasting PPP-RTK corrections. The method includes acquiring grid point information; determining correction information according to a dynamic partitioning method based on the satellite's orbital direction, satellite signal coverage area, and the grid point information; the correction information includes a header file, the grid point information, tropospheric corrections, and oblique ionospheric corrections; and uploading the correction information to an uplink station to control satellite broadcasting of the correction information. This invention solves the problems of satellite bandwidth allocation, correction area division, and correction data format encoding for PPP-RTK atmospheric delay corrections broadcast in low-Earth orbit satellite broadcasting areas.

This invention discloses an automatic method for extracting ionospheric TEC anomalies based on global seismic events. The method acquires a global seismic event catalog, global ionospheric TEC data, and space weather indices, performs data preprocessing, and for each seismic event, selects those occurring during quiescent periods based on space weather conditions. Ionospheric TEC data within fixed observation windows and space windows are extracted and differiated with contemporaneous background TEC data to obtain ionospheric TEC difference values. Normalized statistical analysis is then performed on the quiescent period TEC difference values ​​and their frequency corresponding to global seismic events to obtain statistical characteristics of ionospheric TEC anomaly rates and frequencies. Random seismic events are generated globally for comparison. This method is comprehensive, programmable, and suitable for automated processing of large global samples, possessing good scientific rigor, scalability, and engineering application value.

The invention discloses a broadcast ionosphere model adaptive selection method for a GNSS (Global Navigation Satellite System) terminal, which comprises the following steps: receiving navigation message information broadcasted by each GNSS system in real time, and obtaining various broadcast ionosphere model parameters; calculating the geographic latitude corresponding to the ionospheric puncture point according to the current position coordinate of the user, the satellite elevation angle and the azimuth angle; according to the geographic latitude and the time period of the current calculation time, selecting an optimal broadcast ionosphere model under the current space-time condition; and according to the optimal broadcast ionosphere model and the broadcast ionosphere model parameters, calculating to obtain a correction value of the ionosphere delay amount. According to the invention, by selecting the optimal broadcast ionosphere model in the current environment, the positioning precision is improved, the error of a certain model is prevented from being amplified, the advantages of each broadcast ionosphere model are exerted, and the stability and consistency of system positioning are improved; dynamic environment change is supported, and good space-time adaptive capacity is achieved; and model selection is realized only at a user side without modifying the broadcasting format of the existing satellite message.

The invention discloses a method for forecasting ionized layer delay by combining a ConvLSTM and ViT fusion model. The method comprises the following steps: acquiring single-moment TEC data provided by an ionized layer empirical model in a target area, and VTEC data and a mask layer acquired by a ground GNSS site in the same period; inputting the acquired data into the trained fusion model, and outputting a TEC forecast result in a future preset time period; wherein the fusion model is combined with a ConvLSTM model and a ViT model. According to the method, local space-time details are captured through ConvLSTM, and the structure of capturing global dependence and multi-input fusion auxiliary information in combination with ViT is adopted, so that the delay information of the ionosphere can be well predicted under the active condition of the ionosphere.

The invention provides a GNSS-R wide swath ionosphere vertical total electron content measurement method, device and equipment and a medium. The method comprises the following steps: measuring the propagation delay inequality of a direct signal and a reflection signal from the same navigation satellite through a satellite-borne GNSS-R dual-frequency receiver; calculating a first VTEC value below the satellite orbit at the specular reflection point according to the propagation delay inequality; according to the two-dimensional ionosphere model and the first VTEC value, calculating a second VTEC value and a third VTEC value below the satellite orbit at the incident puncture point and the reflection puncture point; for each target point in the specular reflection point, the incident puncture point and the reflection puncture point, calculating the total VTEC value of the vertical path based on the three-dimensional ionosphere model and the VTEC value below the satellite orbit at the target point, and taking the total VTEC value as the VTEC observation result of the single navigation satellite; and fusing VTEC observation results of a plurality of navigation satellites to generate wide-swath ionized layer VTEC observation data.

The invention provides an ionosphere parameter inversion method and device, and belongs to the technical field of short-wave communication, and the method comprises the steps: enabling the actually-measured amplitude and phase difference to be mapped into a low-order wave mode excitation amplitude and a complex feature value through the synchronous receiving of the same very-low-frequency transmitting station signal at different transmitting and receiving distances on the ground; and then fast convergence solution is carried out by using a Newton-Raphson method, and ionosphere electron concentration parameters are circularly corrected by using a gradient descent method, so that the residual error of a theoretical calculation result and a measured value is stabilized within a preset threshold, thereby getting rid of excessive dependence on a traditional ionosphere mode, remarkably improving the inversion precision of low ionosphere parameters, and improving the reliability of the low ionosphere. The method has the advantages of simple operation, low cost and long-term operation.

The application discloses a mixed order regional ionospheric model construction method, which comprises the following steps: establishing a mixed order regional ionospheric model based on observation data; according to whether the ionosphere is in a quiet period or an active period, using different orders and model parameter update intervals to establish corresponding regional ionospheric models; designing a model order self-adaptive adjustment strategy to realize switching of the mixed order regional ionospheric model to a first-order polynomial model or a second-order polynomial model; and unifying generation approaches of mixed order regional ionospheric model parameters to realize estimation of first-order / second-order polynomial model parameters. The mixed order regional ionospheric model has the advantages that, compared with the first-order polynomial method, the problem of low modeling accuracy in the ionosphere active period is overcome; compared with the second-order polynomial method, the limitations of model overfitting and modeling instability in the ionosphere quiet period are avoided; and in terms of parameter generation and modeling accuracy, the mixed order regional ionospheric model combines the advantages of the second-order model and the first-order model, thereby ensuring the accuracy and applicability of regional ionospheric modeling.

A method and system for high-quality ionosphere spatial interpolation is designed to improve outlier detection and accuracy estimation in regional networks to accelerate convergence of PPP-RTK positioning services. STEC estimations from a regional receiver network undergo a sequence of anomaly detection steps including jump detection by Alpha-Beta filtering, outlier detection by a posteriori residuals from differenced STEC spatial approximation and by cross-validation of differenced STEC values. The method further includes generating an estimation of ionosphere activity indicator which is further used in outlier detection. On-the-fly estimation of the indicator provides more robust STEC outliers detection because it accounts for systematic periodic changes in ionosphere activity and for spatial correlation of the ionosphere.

Various techniques for wireless communication are disclosed. In an aspect, a user equipment (UE) can receive a signal of a first frequency band from a satellite vehicle (SV), estimate a first ionospheric delay residual error based on the signal of the first frequency band, calculate a first pseudorange measurement and a first carrier phase measurement based on the first ionospheric delay residual error, and estimate a position using the first pseudorange measurement and the first carrier phase measurement. In some aspects, the ionospheric delay residual error is estimated via a Klobuchar equation. In some aspects, the position is estimated using ultra-long baseline real-time kinematic (RTK) positioning.

This invention discloses a method, apparatus, device, and medium for calculating electromagnetic waves in a geo-ionospheric waveguide. It obtains the pulse operating frequency set and medium parameter set of electromagnetic waves in a layered geo-ionospheric waveguide model. The electromagnetic waves consist of multiple electromagnetic waves to be calculated in different frequency bands. Based on the pulse operating frequency and medium parameter set of each electromagnetic wave to be calculated, a target calculation method for the receiving electric field is determined from multiple preset electric field calculation methods. The target transmission coefficient is determined based on the transmitting electric field at the emission point of each electromagnetic wave to be calculated, the layered geo-ionospheric waveguide model, and the target calculation method for the receiving electric field. The time-domain waveform of the electromagnetic waves is calculated based on the target transmission coefficients of multiple electromagnetic waves to be calculated and a preset electromagnetic pulse frequency domain calculation formula. In this way, the most suitable target calculation method for the receiving electric field is determined based on the pulse operating frequency set and medium parameter set of the electromagnetic waves to be calculated in different frequency bands, enabling accurate and rapid calculation of electromagnetic waves in a layered geo-ionospheric waveguide.