CN117349575B - Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field - Google Patents
Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field Download PDFInfo
- Publication number
- CN117349575B CN117349575B CN202311642421.XA CN202311642421A CN117349575B CN 117349575 B CN117349575 B CN 117349575B CN 202311642421 A CN202311642421 A CN 202311642421A CN 117349575 B CN117349575 B CN 117349575B
- Authority
- CN
- China
- Prior art keywords
- heating
- ionosphere
- frequency
- low frequency
- current distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000010438 heat treatment Methods 0.000 title claims abstract description 166
- 239000005433 ionosphere Substances 0.000 title claims abstract description 108
- 230000005855 radiation Effects 0.000 title claims abstract description 73
- 238000004364 calculation method Methods 0.000 title claims abstract description 33
- 230000005284 excitation Effects 0.000 title claims description 13
- 238000009826 distribution Methods 0.000 claims abstract description 63
- 238000000034 method Methods 0.000 claims description 42
- 230000006870 function Effects 0.000 claims description 37
- 238000003860 storage Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 230000005358 geomagnetic field Effects 0.000 claims description 10
- 230000005684 electric field Effects 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 13
- 230000006872 improvement Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Abstract
The invention discloses a calculation method and a device for exciting a very low frequency radiation field by heating a difference frequency ionosphere, comprising the following steps: calculation of current distribution: calculating current distribution generated by difference frequency ionosphere modulation heating; calculation of electric dipole radiation field: calculating the very low frequency radiation field excited by the electric dipole in the anisotropic ionized layer, and taking the very low frequency radiation field as a kernel function of current distribution; calculation of the radiation field excited by difference frequency ionosphere heating: and calculating the relative coordinates between each current unit and the observation point in the current distribution, multiplying the current distribution by a kernel function, and integrating the heating area to obtain the very low frequency radiation field excited by the differential frequency ionosphere heating. The invention provides an effective calculation method of the very low frequency radiation field excited by the difference frequency ionosphere heating, which has the advantages of less calculation time consumption and high calculation precision, and can provide theoretical support for practical engineering application of the difference frequency ionosphere heating.
Description
Technical Field
The invention relates to the technical field of antennas and electromagnetic fields, in particular to a calculation method and a calculation device for a difference frequency ionosphere heating excitation very low frequency radiation field.
Background
Very low frequency (Very Low Frequency, VLF) is an electromagnetic wave frequency band between 3kHz and 30kHz, and very low frequency wave has important value for practical applications such as remote communication navigation, ionosphere detection, geological exploration and the like, however, most of traditional ground-based very low frequency radiation antennas are required to be built with large ground emission stations, so that the cost is huge and the radiation efficiency is very limited. High Frequency (HF) electromagnetic waves are radio waves with frequencies in the range of 3MHz to 30MHz, which frequency band is also called the short wave band, and are mainly used in the aeronautics industry, near field communication, radio operators, and weather broadcasting stations, etc. Artificial amplitude modulation methodThe ionosphere is heated by using frequency electromagnetic waves, which is an emerging technology for generating very low frequency waves, and natural current in the ionosphere is modulated and oscillated by manual modulation and heating, so that a virtual very low frequency radiation antenna is formed in the ionosphere. The very low frequency signal excited by the mode has the advantages of small volume, flexible maneuvering, difficult attack and the like compared with a ground-based very low frequency transmitting system if the mode is used for communication. The differential frequency ionosphere heating is a typical heating method for heating the ionosphere by manual modulation, and the method is based on a high ionosphere prime power nonlinear heating theory, and adopts two high-frequency continuous waves with different frequencies to perform modulation heating on the ionosphere (namely, a heating device is divided into two groups, and one group has the emission frequency ofIs transmitted by another group +.>Is a high frequency signal of>For the operating frequency of the excited very low frequency signal), the very low frequency radiation signal can be effectively excited under the condition that the natural current of the background of the ionized layer is weak.
The computational problem of the very low frequency radiation field generated by differential frequency ionosphere heating has not been fully studied. Since the ionosphere in which the modulated heating region is located exhibits very strong anisotropic properties in the very low frequency band, the ionosphere should be considered as an anisotropic plasma medium, but this important factor has not been considered in the past. The anisotropic nature of the ionosphere in the very low frequency band makes modeling characterization and mathematical formula derivation of the very low frequency radiation field in the plasma medium more complex and difficult, and a specific and effective method for calculating the very low frequency radiation field excited by heating the difference frequency ionosphere under the anisotropic condition is not yet available.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a method and a device for calculating a very low frequency radiation field excited by heating a difference frequency ionosphere, which can effectively solve the problem of calculating the very low frequency radiation field excited by heating the difference frequency ionosphere under the anisotropic condition.
In order to achieve the above purpose, the present invention adopts the following technical scheme: a calculation method for a difference frequency ionosphere heating excitation very low frequency radiation field comprises the following steps:
(1) According to parameters of the ionosphere and the heating device, calculating current distribution generated by difference frequency ionosphere modulation heating in the heating area;
(2) Calculating a very low frequency (VLF: 3-30 kHz) radiation field excited by an electric dipole in any direction in the anisotropic ionosphere, and taking the radiation field as a kernel function of current distribution generated by differential frequency ionosphere modulation heating;
(3) According to the current distribution obtained in the step (1), calculating the relative coordinates between each current unit and the observation point in the current distribution, multiplying the current distribution by the kernel function obtained in the step (2), and integrating the whole heating area to obtain the very low frequency radiation field excited by the difference frequency ionosphere heating.
Further, the heating device in the step (1) is two high-frequency (HF: 3-30 MHz) heating waves, namely a first high-frequency heating wave and a second high-frequency heating wave, and the working frequencies are respectivelyAnd->Wherein->Is the operating frequency of the excited very low frequency electromagnetic wave.
Further, in the step (1), the parameters of the ionosphere include an ionosphere electron density and a geomagnetic field strength; the parameters of the heating device include the operating frequency of the high-frequency heating wave, the operating frequency of the excited very low frequency electromagnetic wave, the electric field strength of the high-frequency heating wave, and the phase difference between the first high-frequency heating wave and the second high-frequency heating wave.
Further, the heating area in the step (1) is located in the F layer of the ionosphere, an xyz three-dimensional coordinate system is established, a normal direction of the heating area perpendicular to the ground is defined as a z-axis direction, and current distribution generated by differential frequency ionosphere modulation heating is expressed as:
wherein,representing the position vector of each current cell in the current distribution, respectively>Representing the current distribution resulting from differential ionosphere heating, which at the same time has +.>Component sum->A component; />，/>Respectively indicate->Direction and->Unit vector of direction,/>Representing the imaginary part of the complex number,/, and>indicating the charge of the electron, ">Indicating electron quality->Indicating the operating angular frequency of the high-frequency heating wave,indicating electron cyclotron frequency, ">Representing the intensity of the geomagnetic field +.>Represents the electron density of the ionosphere, < >>Represents the electric field strength of the high-frequency heating wave, < + >>Representing the phase difference between the first high-frequency heating wave and the second high-frequency heating wave, +.>And->The electric field phases of the first high-frequency heating wave and the second high-frequency heating wave are respectively.
Further, in the step (2), a three-dimensional fourier transform is performed on maxwell's equations satisfied by electric dipoles located at any angle with respect to the geomagnetic field in the uniform anisotropic ionosphere to obtain a middle edge of the ionosphereDirection and edge->Very low frequency radiation fields excited by directional electric dipoles, which are used as current distributions +.>Component sum->Kernel functions of the components are combined +.>Andand (3) representing.
Further, in the step (3), the current is distributedComponent sum->The components are multiplied and summed with their corresponding kernel functions respectively, and the components are integrated as a integrand function over the entire heating area to obtain a very low frequency radiation field excited by differential frequency ionosphere modulation heating, the expression of which is:
wherein,position vector representing the viewpoint, +.>Representing the very low frequency radiation field excited by the difference frequency ionosphere heating, the corresponding field component and kernel function>，/>The components of the representation are consistent; />Representing the volume of the heating zone>Representing the volumetric differentiation of the heating zone,/->，/>Respectively representing the current distribution generated by differential ionosphere modulated heating>Component sum->Component (F)>Representing wave impedance in vacuum, +.>Representing wavenumber in vacuum, < >>A diagonal component representing an ionospheric dielectric constant matrix; />For the relative coordinates between the current cells and the observation point, wherein +.>Representing the relative lateral propagation distance between the current cell and the observation point,/->Indicating the relative azimuth angle between the current cell and the observation point, < >>Representing the relative longitudinal propagation distance between the current cell and the observation point.
Further, the relative coordinates between each current unit and the observation pointThe relative propagation distance between each current unit in the current distribution generated by the difference frequency ionosphere modulation heating and the observation point is calculated and converted into a cylindrical coordinate system.
The invention also provides a computing device for exciting the very low frequency radiation field by heating the difference frequency ionosphere, which comprises:
the current distribution calculation module is used for calculating current distribution generated by difference frequency ionosphere modulation heating in the heating area according to parameters of the ionosphere and the heating device;
the kernel function calculation module is used for calculating a very low frequency radiation field excited by an electric dipole in any direction in the anisotropic ionized layer and taking the very low frequency radiation field as a kernel function of current distribution generated by modulating and heating the difference frequency ionized layer;
and the very low frequency radiation field calculation module is used for calculating the relative coordinates between each current unit and the observation point in the current distribution according to the current distribution, multiplying the current distribution by the kernel function and integrating the whole heating area to obtain the very low frequency radiation field excited by the differential frequency ionosphere heating.
The invention also provides a computer readable storage medium storing a computer program which when executed by a processor implements a method of calculating a difference frequency ionosphere heating excitation very low frequency radiation field as described above.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the calculation method of the difference frequency ionosphere heating excitation very low frequency radiation field when executing the program.
The beneficial effects of the invention are as follows:
(1) The effective calculation method of the very low frequency radiation field excited by the difference frequency ionosphere heating is provided, and the current distribution generated by the difference frequency modulation heating and the kernel function are integrated to directly estimate the radiation field, so that the calculation complexity is reduced, the calculation time is less, the accuracy is high, and theoretical guidance and basis can be provided for practical engineering application;
(2) The method takes the key factor of anisotropy of the ionosphere in the very low frequency range into consideration, fills the blank of the method for calculating the very low frequency radiation field generated by ionosphere modulation heating under the anisotropic condition to a certain extent, has the advantages of clear physical meaning, less calculation time consumption, high calculation precision and the like, and can be used for analysis and calculation in practical difference frequency ionosphere heating application;
(3) The method has good ductility, and can be easily extended to calculate radiation fields excited by other ionosphere heating modes by determining current distribution generated by ionosphere modulation heating and relative propagation distance between current elements and observation points.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort to a person skilled in the art.
FIG. 1 is a schematic flow chart of an embodiment of the method of the present invention;
FIG. 2 is a schematic diagram of the geometry of a difference frequency ionosphere heating excitation very low frequency radiation field;
FIG. 3 is a schematic diagram showing the relative propagation distances between the current cell and the observation point;
FIG. 4 is a schematic block diagram of an embodiment of the apparatus of the present invention;
fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the invention. Rather, they are merely examples of apparatus and methods consistent with aspects of the invention as detailed in the accompanying claims.
The present invention will be described in detail with reference to the accompanying drawings. The features of the examples and embodiments described below may be combined with each other without conflict.
The invention relates to a calculation method for exciting a very low frequency radiation field by heating a difference frequency ionosphere, which is shown in figure 1 and comprises the following steps:
(1) Calculation of current distribution: according to the parameters of the ionized layer and the heating device, the current distribution generated by the difference frequency ionized layer modulation heating is calculated.
Wherein the parameters of the ionosphere include ionosphere electron density and geomagnetic field intensity, which are generally related to the height at which the heating region is located; the parameters of the heating device include the operating frequency of the high-frequency heating wave, the operating frequency of the excited very low frequency electromagnetic wave, the electric field strength of the high-frequency heating wave, and the phase difference between the first high-frequency heating wave and the second high-frequency heating wave.
FIG. 2 is a schematic diagram showing the geometry of the excited very low frequency radiation field by the difference frequency ionosphere heating in accordance with the embodiment of the present invention; wherein the heating area is positioned on the F layer of the ionized layer, and the normal direction of the heating area perpendicular to the ground is set as follows when the xyz three-dimensional coordinate system is establishedAxial direction, geomagnetic field and +.>The included angle of the shaft is set as->. The operating frequencies of the two sets of heating devices are respectively set to +.>And->Wherein->For the operating frequency of the second high-frequency heating wave, < >>For the operating frequency of the second high-frequency heating wave, < >>Is the operating frequency of the excited very low frequency electromagnetic wave. Thus, based on the theory of high ionospheric mass dynamic nonlinear heating, the current distribution resulting from differential ionospheric modulation heating can be expressed as:
wherein,representing the position vector of each current cell in the current distribution, respectively>Representing the current distribution resulting from differential ionosphere heating, which at the same time has +.>Component sum->A component; />，/>Respectively indicate->Direction and->Unit vector of direction,/>Representing the imaginary part of the complex number,/, and>indicating the charge of the electron, ">Indicating electron quality->Indicating the operating angular frequency of the high-frequency heating wave,indicating electron cyclotron frequency, ">Representing the intensity of the geomagnetic field +.>Represents the electron density of the ionosphere, < >>Represents the electric field strength of the high-frequency heating wave, < + >>Representing the phase difference between the first high-frequency heating wave and the second high-frequency heating wave, +.>And->The electric field phases of the first high-frequency heating wave and the second high-frequency heating wave are respectively. The parameters can be determined according to the actual parameters of the ionosphere and the heating device, and then the current distribution generated by the differential frequency ionosphere modulation heating can be obtained according to the expression of the current distribution.
(2) Calculation of electric dipole radiation field: calculating a very low frequency radiation field excited by an electric dipole in any direction in the anisotropic ionized layer, and taking the very low frequency radiation field as a kernel function of current distribution generated by heating;
specifically, the Maxwell equation satisfied by the electric dipole at any angle with the geomagnetic field in the uniform anisotropic ionosphere is subjected to three-dimensional Fourier transform to obtain the middle edge of the ionosphereDirection and edge->Very low frequency radiation fields excited by directional electric dipoles, which are used as current distributions generated by differential ionosphere modulated heating>Component sum->Kernel functions of the components are combined +.>And->And (3) representing.
(3) Computing a very low frequency radiation field excited by difference frequency ionosphere heating: according to the current distribution obtained in the step (1), calculating the relative coordinates between each current unit and the observation point in the current distribution, multiplying the current distribution by the kernel function obtained in the step (2), and integrating the whole heating area to obtain the very low frequency radiation field excited by the difference frequency ionosphere heating.
FIG. 3 is a schematic diagram showing the relative propagation distances between a current cell and an observation point according to an embodiment of the present invention; the relative coordinates between each current unit and the observation point can be obtained by calculating the relative propagation distance between each current unit and the observation point in the current distribution generated by the difference frequency ionosphere modulation heating and transforming the relative propagation distance to the observation point into a cylindrical coordinate systemAs shown in fig. 3; wherein (1)>Representing the relative lateral propagation distance between the current cell and the observation point,/->Indicating the relative azimuth angle between the current cell and the observation point, < >>Representing the relative longitudinal propagation distance between the current cell and the observation point.Then, multiplying and summing the y component and the x component of the current distribution obtained in the step (1) with the corresponding kernel functions obtained in the step (2), respectively, and integrating the y component and the x component as the integrated functions over the whole heating area to obtain a very low frequency radiation field excited by the differential ionosphere modulation heating, wherein the very low frequency radiation field is expressed as:
wherein,position vector representing the viewpoint, +.>Representing the very low frequency radiation field excited by the difference frequency ionosphere heating, the corresponding field component and kernel function>，/>The components of the representation are consistent; as shown in fig. 2->Representing the volume of the heating zone>，/>，/>Indicating the length, width, height dimensions of the heating zone, respectively,/->Representing the volumetric differentiation of the heating zone,/->Representing ionosphere lower boundary to groundHeight of face, < >>， />Respectively representing the y-component and the x-component of the current distribution generated by differential ionosphere modulated heating,/->Representing wave impedance in vacuum, +.>Representing wavenumber in vacuum, < >>Representing the diagonal components of the ionospheric dielectric constant matrix.
The invention can be applied to the great application of generating very low frequency signals by adopting a manual modulation and heating ionosphere means to realize the very low frequency underwater communication, navigation and the like. The difference frequency ionosphere heating is a typical heating method for heating the ionosphere by manual modulation, and can realize effective radiation of very low frequency signals under the condition that the natural current of the ionosphere is weak. Therefore, the accurate calculation of the very low frequency radiation field excited by the differential frequency ionosphere heating has important guiding significance for device optimization design, radiation efficiency improvement, optimal parameter determination and the like in practical engineering application.
The invention also provides a computing device for exciting a very low frequency radiation field by difference frequency ionosphere heating, as shown in fig. 4, comprising:
the current distribution calculation module is used for calculating current distribution generated by difference frequency ionosphere modulation heating in the heating area according to parameters of the ionosphere and the heating device;
the kernel function calculation module is used for calculating a very low frequency radiation field excited by an electric dipole in any direction in the anisotropic ionized layer and taking the very low frequency radiation field as a kernel function of current distribution generated by modulating and heating the difference frequency ionized layer;
and the very low frequency radiation field calculation module is used for calculating the relative coordinates between each current unit and the observation point in the current distribution according to the current distribution, multiplying the current distribution by the kernel function and integrating the whole heating area to obtain the very low frequency radiation field excited by the differential frequency ionosphere heating.
It should be noted that, the embodiment of the apparatus shown in this embodiment is matched with the content of the embodiment of the method, and reference may be made to the content of the embodiment of the method, which is not described herein again.
The present invention also provides a computer readable storage medium storing a computer program operable to perform a method of calculating a difference frequency ionosphere heating excitation very low frequency radiation field as provided in fig. 1 above.
The invention also provides a schematic block diagram of the electronic device shown in fig. 5, which corresponds to fig. 1. At the hardware level, as shown in fig. 5, the electronic device includes a processor, an internal bus, a network interface, a memory, and a nonvolatile storage, and may of course include hardware required by other services. The processor reads the corresponding computer program from the nonvolatile memory into the memory and then runs the computer program to realize the calculation method of the difference frequency ionosphere heating excitation very low frequency radiation field shown in the figure 1. Of course, other implementations, such as logic devices or combinations of hardware and software, are not excluded from the present invention, that is, the execution subject of the following processing flows is not limited to each logic unit, but may be hardware or logic devices.
Improvements to one technology can clearly distinguish between improvements in hardware (e.g., improvements to circuit structures such as diodes, transistors, switches, etc.) and software (improvements to the process flow). However, with the development of technology, many improvements of the current method flows can be regarded as direct improvements of hardware circuit structures. Designers almost always obtain corresponding hardware circuit structures by programming improved method flows into hardware circuits. Therefore, an improvement of a method flow cannot be said to be realized by a hardware entity module. For example, a programmable logic device (ProgrammableLogic Device, PLD) (e.g., field programmable gate array (Field Programmable Gate Array, FPGA)) is an integrated circuit whose logic function is determined by the programming of the device by a user. A designer programs to "integrate" a digital system onto a PLD without requiring the chip manufacturer to design and fabricate application-specific integrated circuit chips. Moreover, nowadays, instead of manually manufacturing integrated circuit chips, such programming is mostly implemented by using "logic compiler" software, which is similar to the software compiler used in program development and writing, and the original code before being compiled is also written in a specific programming language, which is called Hardware description language (Hardware DescriptionLanguage, HDL), and HDL is not only one but a plurality of kinds, such as ABEL (Advanced BooleanExpression Language), AHDL (Altera Hardware Description Language), confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java HardwareDescription Language), lava, lola, myHDL, PALASM, RHDL (Ruby Hardware DescriptionLanguage), etc., VHDL (Very-High-Speed Integrated CircuitHardware Description Language) and Verilog are currently most commonly used. It will also be apparent to those skilled in the art that a hardware circuit implementing the logic method flow can be readily obtained by merely slightly programming the method flow into an integrated circuit using several of the hardware description languages described above.
The controller may be implemented in any suitable manner, for example, the controller may take the form of, for example, a microprocessor or processor and a computer readable medium storing computer readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, application specific integrated circuits (Application Specific Integrated Circuit, ASIC), programmable logic controllers, and embedded microcontrollers, examples of which include, but are not limited to, the following microcontrollers: ARC 625D, atmel AT91SAM, microchip PIC18F26K20, and Silicone Labs C8051F320, the memory controller may also be implemented as part of the control logic of the memory. Those skilled in the art will also appreciate that, in addition to implementing the controller in a pure computer readable program code, it is well possible to implement the same functionality by logically programming the method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc. Such a controller may thus be regarded as a kind of hardware component, and means for performing various functions included therein may also be regarded as structures within the hardware component. Or even means for achieving the various functions may be regarded as either software modules implementing the methods or structures within hardware components.
The system, apparatus, module or unit set forth in the above embodiments may be implemented in particular by a computer chip or entity, or by a product having a certain function. One typical implementation is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
For convenience of description, the above devices are described as being functionally divided into various units, respectively. Of course, the functions of each element may be implemented in the same piece or pieces of software and/or hardware when implementing the present invention.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of computer-readable media.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises the element.
It will be appreciated by those skilled in the art that embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, the present specification may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present description can take the form of a computer program product on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
The description may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for the device embodiments, since they are substantially similar to the method embodiments, the description is relatively simple, and reference is made to the description of the method embodiments in part.
The above embodiments are merely for illustrating the design concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement the same, the scope of the present invention is not limited to the above embodiments. Therefore, all equivalent changes or modifications according to the principles and design ideas of the present invention are within the scope of the present invention.
Claims (8)
1. A method for calculating a difference frequency ionosphere heating excitation very low frequency radiation field, comprising the steps of:
(1) According to parameters of the ionosphere and the heating device, calculating current distribution generated by difference frequency ionosphere modulation heating in the heating area;
the heating area in the step (1) is positioned on an F layer of the ionosphere, an xyz three-dimensional coordinate system is established, a normal direction of the heating area perpendicular to the ground is defined as a z-axis direction, and current distribution generated by differential frequency ionosphere modulation heating is expressed as:
wherein r 'represents a position vector of each current unit in the current distribution, and J (r') represents the current distribution generated by differential frequency ionosphere heating, which has both a y component and an x component; y, x represents a unit vector in the y direction and the x direction, i represents an imaginary part of a complex number, e represents a charge amount of electrons, m e Representing electron mass, omega h Indicating the operating angular frequency, ω, of the high-frequency heating wave e ＝eB 0 /m e Indicating electron cyclotron frequency, B 0 Represents the intensity of the geomagnetic field, N e Representing the electron density of the ionosphere, E 0 Represents the electric field strength of the high-frequency heating wave, phi=phi 1 -φ 2 Represents a phase difference phi between the first high-frequency heating wave and the second high-frequency heating wave 1 And phi 2 The electric field phases of the first high-frequency heating wave and the second high-frequency heating wave are respectively;
(2) Calculating a very low frequency radiation field excited by an electric dipole in any direction in the anisotropic ionosphere, and taking the very low frequency radiation field as a kernel function of current distribution generated by modulating and heating the difference frequency ionosphere;
(3) According to the current distribution obtained in the step (1), calculating the relative coordinates between each current unit and the observation point in the current distribution, multiplying the current distribution by the kernel function obtained in the step (2), and integrating the whole heating area to obtain a very low frequency radiation field excited by the differential frequency ionosphere heating;
in the step (3), the y component and the x component of the current distribution are multiplied by the corresponding kernel functions respectively and summed, so that the y component and the x component are used as the integrated functions and the integral is carried out on the whole heating area, and the very low frequency radiation field excited by the differential frequency ionosphere modulation heating can be obtained, wherein the expression is as follows:
wherein r represents the position vector of the observation point, E (r) represents the very low frequency radiation field excited by the difference frequency ionosphere heating, and the corresponding field component and the kernel function K (y) ，K (x) The components of the representation are consistent; v 'represents the volume of the heating region, dV' represents the volumetric differential of the heating region, J y (r′)，J x (r') represents the y-component and the x-component, η, respectively, of the current distribution produced by differential ionosphere modulated heating 0 Representing wave impedance, k, in vacuum 0 Representing wavenumber in vacuum, β represents the pair of ionospheric dielectric constant matricesAngular line components;is the relative coordinate between each current unit and the observation point, wherein ρ is as follows * Representing the relative lateral propagation distance between the current cell and the observation point,/->Indicating the relative azimuth angle, z, between the current cell and the observation point * Representing the relative longitudinal propagation distance between the current cell and the observation point.
2. The method of claim 1, wherein the heating device in the step (1) is two high-frequency heating waves, namely a first high-frequency heating wave and a second high-frequency heating wave, and the working frequencies are f h +f and f h Where f is the operating frequency of the excited very low frequency electromagnetic wave.
3. The method of claim 2, wherein in the step (1), the parameters of the ionosphere include ionosphere electron density and geomagnetic field strength; the parameters of the heating device include the operating frequency of the high-frequency heating wave, the operating frequency of the excited very low frequency electromagnetic wave, the electric field strength of the high-frequency heating wave, and the phase difference between the first high-frequency heating wave and the second high-frequency heating wave.
4. The method of claim 1, wherein in the step (2), the maxwell's equations satisfied by electric dipoles located at arbitrary angles to the geomagnetic field in the uniform anisotropic ionosphere are three-dimensionally fourier transformed to obtain the very low frequency radiation fields excited by the electric dipoles in the ionosphere in the y-direction and in the x-direction as the kernel functions of the y-component and the x-component of the current distribution, respectively, and K is used (y) And K (x) And (3) representing.
5. The method of claim 1, wherein the relative coordinates between the current units and the observation pointThe relative propagation distance between each current unit in the current distribution generated by the difference frequency ionosphere modulation heating and the observation point is calculated and converted into a cylindrical coordinate system.
6. A computing device for difference frequency ionosphere heating excitation of very low frequency radiation fields, comprising:
the current distribution calculation module is used for calculating current distribution generated by difference frequency ionosphere modulation heating in the heating area according to parameters of the ionosphere and the heating device;
the kernel function calculation module is used for calculating a very low frequency radiation field excited by an electric dipole in any direction in the anisotropic ionized layer and taking the very low frequency radiation field as a kernel function of current distribution generated by modulating and heating the difference frequency ionized layer;
and the very low frequency radiation field calculation module is used for calculating the relative coordinates between each current unit and the observation point in the current distribution according to the current distribution, multiplying the current distribution by the kernel function and integrating the whole heating area to obtain the very low frequency radiation field excited by the differential frequency ionosphere heating.
7. A computer readable storage medium, characterized in that the storage medium stores a computer program which, when executed by a processor, implements a method for calculating a difference frequency ionosphere heating excitation very low frequency radiation field according to any of claims 1-5.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements a method of calculating a difference frequency ionosphere heating excitation very low frequency radiation field as claimed in any one of claims 1 to 5 when the program is executed by the processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311642421.XA CN117349575B (en) | 2023-12-04 | 2023-12-04 | Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311642421.XA CN117349575B (en) | 2023-12-04 | 2023-12-04 | Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117349575A CN117349575A (en) | 2024-01-05 |
| CN117349575B true CN117349575B (en) | 2024-03-22 |
Family
ID=89357790
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311642421.XA Active CN117349575B (en) | 2023-12-04 | 2023-12-04 | Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117349575B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117574097B (en) * | 2024-01-08 | 2024-04-09 | 之江实验室 | Method and device for calculating radiation field of multipoint feed antenna |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5777476A (en) * | 1995-12-08 | 1998-07-07 | Papadopoulos; Konstantinos | Ground global tomography (CGT) using modulation of the ionospheric electrojets |
| CN113642208A (en) * | 2021-07-11 | 2021-11-12 | 西北工业大学 | Calculation method for radiation field distribution of underwater very-low-frequency symmetrical oscillator antenna array |
| CN113673070A (en) * | 2020-05-14 | 2021-11-19 | 北京机械设备研究所 | Rapid calculation method for radiation electromagnetic field of electrical antenna in any posture |
| CN114397705A (en) * | 2021-12-10 | 2022-04-26 | 西安理工大学 | Method for predicting time-dependent change of very-low-frequency electric wave field intensity with high precision |
| CN114741839A (en) * | 2022-03-02 | 2022-07-12 | 西北工业大学 | FDTD method for analyzing transmission of very-low frequency electromagnetic wave in earth-ionized layer |
| CN115902344A (en) * | 2022-09-28 | 2023-04-04 | 西安电子科技大学杭州研究院 | Current distribution acquisition method for satellite-borne square-ring antenna, computer device and computer-readable storage medium |
| CN116720389A (en) * | 2023-08-10 | 2023-09-08 | 之江实验室 | Calculation method, device and medium for radiation field of spaceborne array antenna |
| CN116776626A (en) * | 2023-07-04 | 2023-09-19 | 西安理工大学 | Method for calculating field intensity of very low frequency electric wave in earth-anisotropic ionized layer waveguide |
| CN116961725A (en) * | 2023-07-04 | 2023-10-27 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Method for calculating propagation delay of very low frequency electromagnetic wave in anisotropic ionized layer |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060044176A1 (en) * | 2004-08-27 | 2006-03-02 | Bae Systems Advanced Technologies, Inc. | ELF/VLF wave generator using a virtual vertical electric dipole |
| JP2007263789A (en) * | 2006-03-29 | 2007-10-11 | Fujitsu Ltd | Electromagnetic wave analysis program, electromagnetic wave analyzer, and electromagnetic wave analysis method |
| US8299936B2 (en) * | 2008-12-18 | 2012-10-30 | Bae Systems Information And Electronic Systems Integration Inc. | Method and apparatus for establishing low frequency/ultra low frequency and very low frequency communications |
-
2023
- 2023-12-04 CN CN202311642421.XA patent/CN117349575B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5777476A (en) * | 1995-12-08 | 1998-07-07 | Papadopoulos; Konstantinos | Ground global tomography (CGT) using modulation of the ionospheric electrojets |
| CN113673070A (en) * | 2020-05-14 | 2021-11-19 | 北京机械设备研究所 | Rapid calculation method for radiation electromagnetic field of electrical antenna in any posture |
| CN113642208A (en) * | 2021-07-11 | 2021-11-12 | 西北工业大学 | Calculation method for radiation field distribution of underwater very-low-frequency symmetrical oscillator antenna array |
| CN114397705A (en) * | 2021-12-10 | 2022-04-26 | 西安理工大学 | Method for predicting time-dependent change of very-low-frequency electric wave field intensity with high precision |
| CN114741839A (en) * | 2022-03-02 | 2022-07-12 | 西北工业大学 | FDTD method for analyzing transmission of very-low frequency electromagnetic wave in earth-ionized layer |
| CN115902344A (en) * | 2022-09-28 | 2023-04-04 | 西安电子科技大学杭州研究院 | Current distribution acquisition method for satellite-borne square-ring antenna, computer device and computer-readable storage medium |
| CN116776626A (en) * | 2023-07-04 | 2023-09-19 | 西安理工大学 | Method for calculating field intensity of very low frequency electric wave in earth-anisotropic ionized layer waveguide |
| CN116961725A (en) * | 2023-07-04 | 2023-10-27 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Method for calculating propagation delay of very low frequency electromagnetic wave in anisotropic ionized layer |
| CN116720389A (en) * | 2023-08-10 | 2023-09-08 | 之江实验室 | Calculation method, device and medium for radiation field of spaceborne array antenna |
Non-Patent Citations (4)
| Title |
|---|
| Radiation Pattern of a VLF Linear Antenna in an Anisotropic Magnetoplasma;Tong He et al.;《IEEE Transaction on Antennas and Propagation》;20220530;全文 * |
| 双频双波束加热电离层激发甚低频/极低频辐射理论分析;杨巨涛;李清亮;王建国;郝书吉;潘威炎;;物理学报;20170131;第66卷(第01期);全文 * |
| 各向异性电离层中甚低频波的近场散射理论和分层介质中的辐射场研究;林朱红;《万方数据库》;20221213;全文 * |
| 电离层中甚低频天线/天线阵及电磁环境建模;张雪薇;《万方数据库》;20221213;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117349575A (en) | 2024-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117349575B (en) | Calculation method and device for difference frequency ionosphere heating excitation very low frequency radiation field | |
| Ergul et al. | Comparison of integral-equation formulations for the fast and accurate solution of scattering problems involving dielectric objects with the multilevel fast multipole algorithm | |
| Rahmat-Samii et al. | A plane-polar approach for far-field construction from near-field measurements | |
| CN106096267B (en) | A kind of EM scattering of cavities characteristic quick calculation method | |
| Ahmed et al. | Numerical dispersion analysis of the unconditionally stable three-dimensional LOD-FDTD method | |
| CN116720389B (en) | Calculation method, device and medium for radiation field of spaceborne array antenna | |
| Mao et al. | A hybrid implicit-explicit spectral FDTD scheme for oblique incidence problems on periodic structures | |
| Oetting et al. | Near-to-far-field transformation by a time-domain spherical-multipole analysis | |
| Ju et al. | Investigation on the characteristics of the envelope FDTD based on the alternating direction implicit scheme | |
| CN103235193A (en) | Numerical method of satellite electromagnetic scattering characteristics within millimeter wave band | |
| Pham-Xuan et al. | On comparison of integral equation approaches for indoor wave propagation | |
| Malik et al. | Alternative mathematical design of vector potential and radiated fields for parabolic reflector surface | |
| Ergul et al. | Iterative solutions of hybrid integral equations for coexisting open and closed surfaces | |
| Wait | Propagation effects for electromagnetic pulse transmission | |
| Liu et al. | Ambiguity analysis and resolution for phase-based 3D source localization under given UCA | |
| Wu et al. | Computationally efficient complex envelope approximate Crank–Nicolson scheme and its open region problem for anisotropic gyrotropic plasma | |
| Conway | New exact solution procedure for the near fields of the general thin circular loop antenna | |
| Zeyde | Generalized approach for the rotating axisymmetric bodies electromagnetic scattering problems | |
| Västberg et al. | Signal intensity in the geometrical optics approximation for the magnetized ionosphere | |
| Heh et al. | Complex-envelope LOD-FDTD method for ionospheric propagation | |
| El Adraoui et al. | Shift operator-TLM method for modeling gyroelectric media | |
| CN116842601B (en) | Method for determining field distribution on arbitrary three-dimensional curved section in circular waveguide | |
| Smith et al. | Transient fields of offset reflector | |
| Yang et al. | A new method to analyze the EM wave propagation characteristics in the hypersonic sheath | |
| Seliga | Numerical full wave solution techniques for the calculation of low frequency plane wave fields in the ionosphere |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |