CN115656946B - Incoherent scattering radar multi-beam calibration method and system for plasma line measurement - Google Patents

Incoherent scattering radar multi-beam calibration method and system for plasma line measurement Download PDF

Info

Publication number
CN115656946B
CN115656946B CN202211587456.3A CN202211587456A CN115656946B CN 115656946 B CN115656946 B CN 115656946B CN 202211587456 A CN202211587456 A CN 202211587456A CN 115656946 B CN115656946 B CN 115656946B
Authority
CN
China
Prior art keywords
calibration
radar
data
electron density
plasma line
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
Application number
CN202211587456.3A
Other languages
Chinese (zh)
Other versions
CN115656946A (en
Inventor
王俊逸
乐新安
王永辉
刘凡玉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202211587456.3A priority Critical patent/CN115656946B/en
Publication of CN115656946A publication Critical patent/CN115656946A/en
Application granted granted Critical
Publication of CN115656946B publication Critical patent/CN115656946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention belongs to the field of radar signal and data processing, particularly relates to a plasma line measurement-based incoherent scattering radar multi-beam calibration method and system, and aims to solve the problem that incoherent scattering radar multi-beam calibration cannot be independently realized conveniently and reliably without depending on external instruments and observation in the prior art. The invention includes: extracting the offset of the frequency of the plasma line at the peak value relative to the receiving frequency, and calculating the peak electron density of the ionized layer; low-pass filtering and down-sampling the original IQ data, calculating a power profile, denoising after accumulation, and calculating the electron density of an ionized layer before calibration by using a signal-to-noise ratio; and circularly calculating calibration coefficients for all the beams, substituting the obtained calibration coefficient sequences of different beams into a radar equation, and completing the calibration of the incoherent scattering radar multi-beam. The method does not depend on external instruments, has the advantages of accuracy, reliability, convenience and the like, and can be used for calibration, ionosphere detection, space weather research and the like of an incoherent scattering radar system.

Description

Incoherent scattering radar multi-beam calibration method and system for plasma line measurement
Technical Field
The invention belongs to the field of radar signal and data processing, and particularly relates to a incoherent scattering radar multi-beam calibration method and system based on plasma line measurement.
Background
The ionosphere is the ionized part of the earth's atmosphere, which serves as the inner boundary of the earth's magnetic layer and is distributed mainly in the region 60km to 1000km above the earth's surface. The ionosphere is coupled with other circle layers through various momentum and energy processes and is a key area for space weather research; meanwhile, the state of the ionized layer influences the propagation of radio waves, which has influence on satellite navigation and positioning, radio communication and the like, and has very important practical significance on human beings.
The incoherent scattering detection radar is the strongest ionosphere foundation detection means at present, can directly and simultaneously detect the plasma density, temperature, drift velocity (electric field) and components on almost the whole ionosphere height with high precision, and indirectly detect the temperature and wind field of the background neutral atmosphere. Incoherent scattering radar detection has the characteristic of high space-time resolution, and the calibration of equipment is critical to ensure the absolute measurement with high accuracy. According to the traditional method for calibrating by observing the altimeter, the altimeter data erected at the periphery of the radar needs to be used for comparing with the simultaneous observation of the incoherent scattering radar to finish calibration, and because the altimeter can only provide large-range observation with fixed direction, the incoherent scattering radar multi-beam is calibrated by using the observation of the altimeter, calibration errors are introduced to influence the data measurement precision by neglecting the space change.
In general, there is a need in the art for a method and a system for incoherent scattering radar multi-beam calibration that can directly perform device calibration through independent observation of incoherent scattering radar itself, so as to achieve convenient and reliable incoherent scattering radar multi-beam calibration.
Disclosure of Invention
In order to solve the above-mentioned problems in the prior art, that is, the prior art cannot independently implement convenient and reliable incoherent scattering radar multi-beam calibration without depending on external instruments and observation, the present invention provides an incoherent scattering radar multi-beam calibration method based on plasma line measurement, the calibration method comprising:
step S10, acquiring original IQ data detected by a plasma line of the incoherent scattering radar, calculating an accumulated signal power spectrum, extracting the offset of the frequency at the peak value of the plasma line relative to the signal receiving frequency, and calculating the peak electron density of an ionized layer based on the plasma line dispersion relation;
step S20, low-pass filtering and down-sampling are carried out on the original IQ data to obtain narrow-band ion line observation data, a power profile of the narrow-band ion line observation data is calculated, background noise is removed, a signal to noise ratio is calculated, and the original electron density of the ionized layer before calibration is calculated according to the signal to noise ratio, a radar equation of ionized layer detection and theoretical values of radar system parameters;
s30, comparing the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by the plasma line observation to obtain a calibration coefficient of radar system parameters;
and S40, circularly executing the steps S10 to S30, sequentially calculating calibration coefficients for observation of all beams to obtain calibration coefficient sequences of different beams, substituting the calibration coefficient sequences into a radar equation, and completing calibration of the incoherent scattering radar multi-beam.
In some preferred embodiments, step S10 includes:
s11, interpreting original IQ data detected by a plasma line of the incoherent scattering radar to obtain a data information file related to experimental parameters;
s12, acquiring an experimental waveform and coding information based on the data information file, decoding each frame of data based on a coding type, calculating an autocorrelation function, performing incoherent accumulation within a time resolution, and performing Fourier transform after adding a Hamming window to the accumulated autocorrelation function to obtain a power spectrum of a plasma line;
and S13, calculating the peak electron density of the ionized layer through the frequency offset at the peak on the plasma linear power spectrum based on the plasma linear dispersion relation under the influence of the magnetic field.
In some preferred embodiments, the plasma line dispersion relationship is expressed as:
Figure 118332DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 326459DEST_PATH_IMAGE002
is the frequency of the electron plasma and,
Figure 967656DEST_PATH_IMAGE003
is the frequency of electron cyclotron resonance,
Figure 653852DEST_PATH_IMAGE004
the offset of the plasma line frequency from the receive center frequency,
Figure 454318DEST_PATH_IMAGE005
is the angle between the beam and the earth's magnetic field.
In some preferred embodiments, the electron plasma frequency is expressed as:
Figure 771030DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 696261DEST_PATH_IMAGE007
in order to be of an electron mass,
Figure 186148DEST_PATH_IMAGE008
is the electron density of the ionosphere and,
Figure 841120DEST_PATH_IMAGE009
which is the dielectric constant in a vacuum, is,
Figure 391050DEST_PATH_IMAGE010
is the electron mass.
In some preferred embodiments, the electron cyclotron frequency is expressed as:
Figure 6839DEST_PATH_IMAGE011
wherein the content of the first and second substances,
Figure 34838DEST_PATH_IMAGE007
in order to be of an electron mass,
Figure 544317DEST_PATH_IMAGE012
in order to observe the magnetic induction intensity locally at the station,
Figure 265148DEST_PATH_IMAGE010
is the electron mass.
In some preferred embodiments, step S20 includes:
step S21, low-pass filtering is carried out on the original IQ data, and the filtered data is subjected to down-sampling according to the data bandwidth to obtain narrow-band ion line observation data;
step S22, calculating a power profile of the narrow-band ion line observation data, accumulating according to time resolution, removing noise from the accumulated power profile, and calculating a signal-to-noise ratio;
s23, generating a radar equation of ionosphere detection under the condition of neglecting the influence of the Debye length;
and S24, substituting the theoretical value of the radar system parameter into the radar equation, and calculating the original electron density of the ionized layer before calibration according to the signal-to-noise ratio.
In some preferred embodiments, the radar equation for ionospheric sounding is expressed as:
Figure 164971DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 934344DEST_PATH_IMAGE014
for the purpose of the signal-to-noise ratio,
Figure 236012DEST_PATH_IMAGE015
in order to be the radius of the electrons,
Figure 190062DEST_PATH_IMAGE016
in order to be the angle of polarization of the antenna,
Figure 842760DEST_PATH_IMAGE017
it is the temperature of the electrons that is,
Figure 150244DEST_PATH_IMAGE018
is the temperature of the ions, and is,
Figure 571999DEST_PATH_IMAGE019
in order to transmit the peak power of the signal,
Figure 696949DEST_PATH_IMAGE020
in order to obtain the gain of the antenna,
Figure 305785DEST_PATH_IMAGE021
in order to be the wavelength of the signal,
Figure 744857DEST_PATH_IMAGE022
for the slant distance between the radar antenna and the detection target,
Figure 958800DEST_PATH_IMAGE023
in order to be the speed of light,
Figure 192336DEST_PATH_IMAGE024
in order for the pulse width to be effective,
Figure 616364DEST_PATH_IMAGE025
is a constant of boltzmann's constant,
Figure 327968DEST_PATH_IMAGE026
in order to be the noise temperature of the radar system,
Figure 661997DEST_PATH_IMAGE027
in order to be a bandwidth of the communication channel,
Figure 66434DEST_PATH_IMAGE028
is the ionospheric raw electron density.
In another aspect of the present invention, a incoherent scattering radar multi-beam calibration system based on plasma line measurement is provided, the calibration system comprising:
the data acquisition module is configured to acquire original IQ data detected by a plasma line of the incoherent scattering radar;
the plasma line observed peak electron density calculation module is configured to extract the offset of the frequency at the plasma line peak relative to the signal receiving frequency and calculate the ionized layer peak electron density based on the plasma line dispersion relation;
the signal-to-noise ratio calculation module for ion line observation is configured to perform low-pass filtering and down-sampling on original IQ data to obtain narrow-band ion line observation data, calculate a power profile of the narrow-band ion line observation data, remove background noise and calculate a signal-to-noise ratio;
the ion line observation original electron density calculation module is configured to calculate the ionized layer original electron density before calibration according to the signal-to-noise ratio, the radar equation of ionized layer detection and the theoretical value of radar system parameters;
the calibration coefficient acquisition module is configured to compare the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by the plasma line observation to obtain the calibration coefficient of the radar system parameter;
the calibration coefficient sequence acquisition module is configured to iteratively calculate calibration coefficients for observation of all beams through the data acquisition module, the peak electron density calculation module for plasma line observation, the signal-to-noise ratio calculation module for ion line observation, the original electron density calculation module for ion line observation and the calibration coefficient acquisition module in sequence to obtain calibration coefficient sequences of different beams;
and the calibration module is configured to substitute the calibration coefficient sequences of the different beams into a radar equation to finish the calibration of the incoherent scattering radar multi-beam.
In a third aspect of the present invention, an electronic device is provided, including:
at least one processor;
and a memory communicatively coupled to at least one of the processors;
wherein the memory stores instructions executable by the processor for execution by the processor to implement the above-described non-coherent scatter radar multi-beam calibration method based on plasma line measurements.
In a fourth aspect of the present invention, a computer-readable storage medium is provided, which stores computer instructions for execution by the computer to implement the above-mentioned incoherent scattering radar multi-beam calibration method based on plasmonic line measurement.
The invention has the beneficial effects that:
(1) Compared with the traditional calibration method based on altimeter observation, the incoherent scattering radar multi-beam calibration method based on plasma line measurement does not need to rely on external instruments and observation, simplifies the calibration process and improves the convenience of incoherent scattering radar calibration.
(2) According to the incoherent scattering radar multi-beam calibration method based on the plasma ray measurement, the observation results of the plasma ray in different beam directions can truly reflect the change of radar system parameters in different observation angles, the actual measurement results of the plasma ray multi-beam are used for replacing the approximate processing of the single-beam measurement of the altimeter, and the accuracy of incoherent scattering radar multi-beam calibration is improved.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting embodiments thereof, made with reference to the accompanying drawings in which:
fig. 1 is a flow chart schematic diagram of a incoherent scattering radar multi-beam calibration method based on plasma line measurement according to the present invention.
Detailed Description
The present application will be described in further detail with reference to the following drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not to be construed as limiting the invention. It should be noted that, for convenience of description, only the portions related to the present invention are shown in the drawings.
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
The invention provides a plasma line measurement-based incoherent scattering radar multi-beam calibration method, which directly calibrates equipment through independent observation of an incoherent scattering radar, has the advantages of convenience and reliability, and can be used for calibrating an incoherent scattering radar system, such as ionosphere multi-parameter monitoring, space weather research and the like.
The invention relates to a multi-beam calibration method of incoherent scattering radar based on plasma line measurement, which comprises the following steps:
step S10, acquiring original IQ data detected by a plasma line of the incoherent scattering radar, calculating an accumulated signal power spectrum, extracting the offset of the frequency at the peak value of the plasma line relative to the signal receiving frequency, and calculating the peak electron density of an ionized layer based on the plasma line dispersion relation;
step S20, low-pass filtering and down-sampling are carried out on the original IQ data to obtain narrow-band ion line observation data, a power profile of the narrow-band ion line observation data is calculated, background noise is removed, a signal to noise ratio is calculated, and the original electron density of the ionized layer before calibration is calculated according to the signal to noise ratio, a radar equation of ionized layer detection and theoretical values of radar system parameters;
s30, comparing the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by the plasma line observation to obtain a calibration coefficient of radar system parameters;
and S40, circularly executing the step S10 to the step S30, sequentially calculating calibration coefficients for the observation of all the wave beams to obtain calibration coefficient sequences of different wave beams, and substituting the calibration coefficient sequences into a radar equation to finish the calibration of the incoherent scattering radar multi-wave beams.
In order to more clearly describe the incoherent scattering radar multi-beam calibration method based on plasma line measurement of the present invention, the following describes each step in the embodiment of the present invention in detail with reference to fig. 1.
The incoherent scattering radar multi-beam calibration method based on the plasma line measurement in the first embodiment of the invention comprises the following steps S10-S40, wherein the following steps are described in detail:
step S10, obtaining original IQ data detected by a plasma line of the incoherent scattering radar, calculating an accumulated signal power spectrum, extracting the offset of the frequency at the peak of the plasma line relative to the signal receiving frequency, and calculating the peak electron density of an ionized layer based on the plasma line dispersion relation, wherein the method specifically comprises the following steps:
and S11, interpreting the original IQ data detected by the incoherent scattering radar plasma line to obtain a data information file related to the experimental parameters.
According to the incoherent scattering detection principle, aiming at a transmitted signal, the system carries out data acquisition with large bandwidth and high sampling rate, and interprets an original data file to obtain a data information file and an original IQ data file which are related to experimental setting.
And S12, acquiring experimental waveforms and coding information based on the data information file, decoding each frame of data based on the coding type, calculating an autocorrelation function, performing incoherent accumulation within a time resolution, adding a Hamming window to the accumulated autocorrelation function, and performing Fourier transform to obtain a power spectrum of a plasma line.
The system acquires experimental waveforms and coding information of data according to a data information file, decodes each frame of data according to a coding type, calculates an autocorrelation function of each frame of data, performs incoherent accumulation within a time resolution, adds a Hamming window to the accumulated autocorrelation function, and finally obtains a corresponding plasma line power spectrum through Fourier transform.
And S13, calculating the ionized layer peak electron density through the frequency offset at the peak value on the plasma line power spectrum based on the plasma line dispersion relation under the magnetic field influence condition.
Wherein, the plasma line dispersion relation is shown as formula (1):
Figure 649862DEST_PATH_IMAGE001
(1)
wherein, the first and the second end of the pipe are connected with each other,
Figure 239192DEST_PATH_IMAGE002
is the frequency of the electron plasma and,
Figure 490045DEST_PATH_IMAGE003
is the frequency of electron cyclotron resonance,
Figure 3066DEST_PATH_IMAGE004
the offset of the plasma line frequency from the receive center frequency,
Figure 339370DEST_PATH_IMAGE005
is the angle between the beam and the earth's magnetic field.
The electron plasma frequency and the electron cyclotron frequency are expressed by the following equations (2) and (3), respectively:
Figure 720672DEST_PATH_IMAGE006
(2)
Figure 560452DEST_PATH_IMAGE011
(3)
wherein, the first and the second end of the pipe are connected with each other,
Figure 306692DEST_PATH_IMAGE007
in order to be of an electron mass,
Figure 802395DEST_PATH_IMAGE008
is the density of electrons in the ionized layer,
Figure 925072DEST_PATH_IMAGE009
is a dielectric constant in a vacuum, and,
Figure 947254DEST_PATH_IMAGE010
in order to be of an electron mass,
Figure 864395DEST_PATH_IMAGE012
is the magnetic induction local to the observation station.
Step S20, low-pass filtering and down-sampling the original IQ data to obtain narrow-band ion line observation data, calculating a power profile of the narrow-band ion line observation data, removing background noise, calculating a signal-to-noise ratio, and calculating the ionized layer original electron density before calibration according to the signal-to-noise ratio, a radar equation of ionized layer detection and theoretical values of radar system parameters, specifically comprising:
and S21, performing low-pass filtering on the original IQ data, and performing down-sampling on the filtered data according to the data bandwidth to obtain narrow-band ion line observation data.
Setting cut-off frequency, carrying out low-pass filtering on the original IQ data with large bandwidth, and then carrying out down-sampling on the filtered data according to the bandwidth to obtain narrow-band ion ray observation data.
And S22, calculating a power section of the narrow-band ion line observation data, accumulating according to the time resolution, removing noise from the accumulated power section, and calculating the signal-to-noise ratio.
And calculating the power profile of the power section, accumulating according to the time resolution, removing background noise from the accumulated power profile, and calculating the signal-to-noise ratio.
Step S23, generating a radar equation of ionosphere detection under the condition of neglecting the influence of Debye length, wherein the formula (4) is as follows:
Figure 909711DEST_PATH_IMAGE013
(4)
wherein the content of the first and second substances,
Figure 508183DEST_PATH_IMAGE014
in order to be the signal-to-noise ratio,
Figure 322555DEST_PATH_IMAGE015
is the radius of the electrons and the electron beam,
Figure 738493DEST_PATH_IMAGE016
in order to be the angle of polarization of the antenna,
Figure 271105DEST_PATH_IMAGE017
it is the temperature of the electrons that is,
Figure 407689DEST_PATH_IMAGE018
is the temperature of the ions, and is,
Figure 342147DEST_PATH_IMAGE019
in order to transmit the peak power of the signal,
Figure 601090DEST_PATH_IMAGE020
in order to obtain the gain of the antenna,
Figure 417736DEST_PATH_IMAGE021
in order to be the wavelength of the signal,
Figure 685906DEST_PATH_IMAGE022
for the slant distance between the radar antenna and the detection target,
Figure 412554DEST_PATH_IMAGE023
it is the speed of the light that is,
Figure 842398DEST_PATH_IMAGE024
for the effective pulse width to be useful,
Figure 411920DEST_PATH_IMAGE025
is a constant of boltzmann's constant,
Figure 952622DEST_PATH_IMAGE026
in order to be the noise temperature of the radar system,
Figure 799356DEST_PATH_IMAGE027
in order to be a bandwidth of the communication channel,
Figure 665681DEST_PATH_IMAGE028
is the ionospheric raw electron density.
And S24, substituting the theoretical value of the radar system parameter into the radar equation, and calculating the original electron density of the ionized layer before calibration according to the signal-to-noise ratio.
And S30, comparing the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by ion line observation to obtain a calibration coefficient of the radar system parameter.
And S40, circularly executing the steps S10 to S30, sequentially calculating calibration coefficients for observation of all beams to obtain calibration coefficient sequences of different beams, substituting the calibration coefficient sequences into a radar equation, and completing calibration of the incoherent scattering radar multi-beam.
Although the foregoing embodiments have described the steps in the foregoing sequence, those skilled in the art will understand that, in order to achieve the effect of the present embodiment, different steps are not necessarily performed in such a sequence, and may be performed simultaneously (in parallel) or in an inverse sequence, and these simple variations are within the scope of the present invention.
A incoherent scattering radar multi-beam calibration system based on plasma line measurement according to a second embodiment of the present invention, the calibration system comprising:
the data acquisition module is configured to acquire original IQ data detected by a non-coherent scattering radar plasma line;
the plasma line observed peak electron density calculation module is configured to extract the offset of the frequency at the plasma line peak relative to the signal receiving frequency and calculate the ionized layer peak electron density based on the plasma line dispersion relation;
the signal-to-noise ratio calculation module for ion line observation is configured to perform low-pass filtering and down-sampling on original IQ data to obtain narrow-band ion line observation data, calculate a power profile of the narrow-band ion line observation data, remove background noise and calculate a signal-to-noise ratio;
the ion line observation original electron density calculation module is configured to calculate the ionized layer original electron density before calibration according to the signal-to-noise ratio, the radar equation of ionized layer detection and the theoretical value of radar system parameters;
the calibration coefficient acquisition module is configured to compare the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by ion line observation to obtain the calibration coefficient of the radar system parameter;
the calibration coefficient sequence acquisition module is configured to iteratively calculate calibration coefficients for observation of all beams through the data acquisition module, the peak electron density calculation module for plasma line observation, the signal-to-noise ratio calculation module for ion line observation, the original electron density calculation module for ion line observation and the calibration coefficient acquisition module in sequence to obtain calibration coefficient sequences of different beams;
and the calibration module is configured to substitute the calibration coefficient sequences of the different beams into a radar equation to finish the calibration of the incoherent scattering radar multi-beam.
It can be clearly understood by those skilled in the art that, for convenience and simplicity of description, the specific working process and related description of the system described above may refer to the corresponding process in the foregoing method embodiment, and details are not described herein again.
It should be noted that, the incoherent scattering radar multi-beam calibration system based on plasma line measurement provided in the foregoing embodiment is only illustrated by dividing the functional modules, and in practical applications, the above functions may be allocated to different functional modules as needed, that is, the modules or steps in the embodiments of the present invention are further decomposed or combined, for example, the modules in the embodiments may be combined into one module, or may be further split into multiple sub-modules, so as to complete all or part of the functions described above. The names of the modules and steps involved in the embodiments of the present invention are only for distinguishing the modules or steps, and are not to be construed as unduly limiting the present invention.
An electronic device of a third embodiment of the present invention includes:
at least one processor;
and a memory communicatively coupled to at least one of the processors;
wherein the memory stores instructions executable by the processor for execution by the processor to implement the above-described non-coherent scatter radar multi-beam calibration method based on plasma line measurements.
A computer-readable storage medium of a fourth embodiment of the present invention stores computer instructions for execution by the computer to implement the incoherent scattering radar multi-beam calibration method based on plasmonic line measurements described above.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes and related descriptions of the storage device and the processing device described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
Those of skill in the art will appreciate that the various illustrative modules, method steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that programs corresponding to the software modules, method steps may be located in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. To clearly illustrate this interchangeability of electronic hardware and software, various illustrative components and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as electronic hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The terms "first," "second," and the like are used for distinguishing between similar elements and not necessarily for describing or implying a particular order or sequence.
The terms "comprises," "comprising," or any other similar term 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.
So far, the technical solutions of the present invention have been described in connection with the preferred embodiments shown in the drawings, but it is apparent to those skilled in the art that the scope of the present invention is not limited to these specific embodiments. Equivalent changes or substitutions of related technical features can be made by those skilled in the art without departing from the principle of the invention, and the technical scheme after the changes or substitutions can fall into the protection scope of the invention.

Claims (10)

1. A multi-beam calibration method for incoherent scattering radar based on plasma line measurement, the calibration method comprising:
s10, acquiring original IQ data detected by a plasma line of the incoherent scattering radar, extracting the offset of the frequency at the peak value of the plasma line relative to the signal receiving frequency, and calculating the electron density of the peak value of the ionized layer based on the plasma line dispersion relation;
step S20, low-pass filtering and down-sampling are carried out on the original IQ data to obtain narrow-band ion line observation data, a power profile of the narrow-band ion line observation data is calculated, background noise is removed, a signal to noise ratio is calculated, and the original electron density of the ionized layer before calibration is calculated according to the signal to noise ratio, a radar equation of ionized layer detection and theoretical values of radar system parameters;
s30, comparing the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by the plasma line observation to obtain a calibration coefficient of radar system parameters;
and S40, circularly executing the step S10 to the step S30, sequentially calculating calibration coefficients for the observation of all the wave beams to obtain calibration coefficient sequences of different wave beams, and substituting the calibration coefficient sequences into a radar equation to finish the calibration of the incoherent scattering radar multi-wave beams.
2. The incoherent scattering radar multi-beam calibration method based on plasmonic line measurements of claim 1, wherein step S10 comprises:
s11, interpreting original IQ data detected by a plasma line of the incoherent scattering radar to obtain a data information file related to experimental parameters;
s12, acquiring an experimental waveform and coding information based on the data information file, decoding each frame of data based on a coding type, calculating an autocorrelation function, performing incoherent accumulation within a time resolution, and performing Fourier transform after adding a Hamming window to the accumulated autocorrelation function to obtain a power spectrum of a plasma line;
and S13, calculating the ionized layer peak electron density through the frequency offset at the peak value on the plasma line power spectrum based on the plasma line dispersion relation under the magnetic field influence condition.
3. The method according to claim 2, characterized in that the plasmonic line dispersion relation, expressed as:
Figure 113201DEST_PATH_IMAGE001
wherein, the first and the second end of the pipe are connected with each other,
Figure 295921DEST_PATH_IMAGE002
is the frequency of the electron plasma and is,
Figure 403554DEST_PATH_IMAGE003
is the frequency of electron cyclotron resonance,
Figure 64343DEST_PATH_IMAGE004
the offset of the plasma line frequency from the receive center frequency,
Figure 878715DEST_PATH_IMAGE005
is the angle between the beam and the earth's magnetic field.
4. The method for incoherent scatter radar multi-beam calibration based on plasma line measurements according to claim 3, characterized in that the electron plasma frequency, expressed as:
Figure 904440DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 437052DEST_PATH_IMAGE007
in order to be of an electron mass,
Figure 963849DEST_PATH_IMAGE008
is the electron density of the ionosphere and,
Figure 898307DEST_PATH_IMAGE009
is a dielectric constant in a vacuum, and,
Figure 94933DEST_PATH_IMAGE010
is the electron mass.
5. The method according to claim 3, characterized in that the electron cyclotron frequency is expressed as:
Figure 114841DEST_PATH_IMAGE011
wherein the content of the first and second substances,
Figure 179749DEST_PATH_IMAGE007
in order to be of an electron mass,
Figure 703134DEST_PATH_IMAGE012
in order to observe the magnetic induction intensity locally at the station,
Figure 398558DEST_PATH_IMAGE010
is the electron mass.
6. The method according to claim 1, wherein step S20 comprises:
step S21, low-pass filtering is carried out on the original IQ data, and the filtered data is subjected to down-sampling according to the data bandwidth to obtain narrow-band ion line observation data;
step S22, calculating a power profile of the narrow-band ion line observation data, accumulating according to time resolution, removing noise from the accumulated power profile, and calculating a signal-to-noise ratio;
s23, generating a radar equation of ionosphere detection under the condition of neglecting the Debye length influence;
and S24, substituting the theoretical value of the radar system parameter into the radar equation, and calculating the original electron density of the ionized layer before calibration according to the signal-to-noise ratio.
7. The method according to claim 6, the radar equation for ionospheric sounding expressed as:
Figure 843446DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 384149DEST_PATH_IMAGE014
for the purpose of the signal-to-noise ratio,
Figure 355516DEST_PATH_IMAGE015
in order to be the radius of the electrons,
Figure 956261DEST_PATH_IMAGE016
in order to be the antenna polarization angle,
Figure 888445DEST_PATH_IMAGE017
it is the temperature of the electrons that is,
Figure 232839DEST_PATH_IMAGE018
is the temperature of the ions, and is,
Figure 58712DEST_PATH_IMAGE019
in order to transmit the peak power, the power,
Figure 95939DEST_PATH_IMAGE020
in order to obtain the gain of the antenna,
Figure 46577DEST_PATH_IMAGE021
as to the wavelength of the signal, is,
Figure 132345DEST_PATH_IMAGE022
for the slant distance between the radar antenna and the detection target,
Figure 750408DEST_PATH_IMAGE023
in order to be the speed of light,
Figure 20852DEST_PATH_IMAGE024
for the effective pulse width to be useful,
Figure 724366DEST_PATH_IMAGE025
is the boltzmann constant, and is,
Figure 348245DEST_PATH_IMAGE026
in order to be the noise temperature of the radar system,
Figure 86394DEST_PATH_IMAGE027
in order to be a bandwidth,
Figure 527740DEST_PATH_IMAGE028
is the ionospheric raw electron density.
8. An incoherent scatter radar multi-beam calibration system based on plasma line measurements, the calibration system comprising:
the data acquisition module is configured to acquire original IQ data detected by a non-coherent scattering radar plasma line;
the plasma line observation peak electron density calculation module is configured to extract the offset of the frequency at the plasma line peak relative to the signal receiving frequency and calculate the ionized layer peak electron density based on the plasma line dispersion relation;
the signal-to-noise ratio calculation module for ion line observation is configured to perform low-pass filtering and down-sampling on the original IQ data to obtain narrow-band ion line observation data, calculate a power profile of the narrow-band ion line observation data, remove background noise and calculate a signal-to-noise ratio;
the ionized layer original electron density calculation module is configured to calculate the ionized layer original electron density before calibration according to the signal to noise ratio, the radar equation of ionized layer detection and the theoretical value of radar system parameters;
the calibration coefficient acquisition module is configured to compare the ionized layer peak electron density calculated by plasma line observation with the ionized layer original electron density peak value before calibration calculated by the plasma line observation to obtain the calibration coefficient of the radar system parameter;
the calibration coefficient sequence acquisition module is configured to iterate through the data acquisition module, the peak electron density calculation module for plasma line observation, the signal-to-noise ratio calculation module for ion line observation, the original electron density calculation module for ion line observation and the calibration coefficient acquisition module to calculate calibration coefficients in sequence for all beam observations, and obtain calibration coefficient sequences of different beams;
and the calibration module is configured to substitute the calibration coefficient sequences of the different beams into a radar equation to finish calibration of the incoherent scattering radar multi-beam.
9. An electronic device, comprising:
at least one processor;
and a memory communicatively coupled to at least one of the processors;
wherein the memory stores instructions executable by the processor for execution by the processor to implement the plasmonic-measurement-based incoherent scatter radar multi-beam calibration method of any of claims 1-7.
10. A computer-readable storage medium having stored thereon computer instructions for execution by the computer to implement the plasmonic-measurement-based incoherent scatter radar multi-beam calibration method of any of claims 1-7.
CN202211587456.3A 2022-12-12 2022-12-12 Incoherent scattering radar multi-beam calibration method and system for plasma line measurement Active CN115656946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211587456.3A CN115656946B (en) 2022-12-12 2022-12-12 Incoherent scattering radar multi-beam calibration method and system for plasma line measurement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211587456.3A CN115656946B (en) 2022-12-12 2022-12-12 Incoherent scattering radar multi-beam calibration method and system for plasma line measurement

Publications (2)

Publication Number Publication Date
CN115656946A CN115656946A (en) 2023-01-31
CN115656946B true CN115656946B (en) 2023-03-10

Family

ID=85017430

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211587456.3A Active CN115656946B (en) 2022-12-12 2022-12-12 Incoherent scattering radar multi-beam calibration method and system for plasma line measurement

Country Status (1)

Country Link
CN (1) CN115656946B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117289234B (en) * 2023-11-24 2024-02-09 中国科学院地质与地球物理研究所 Method for extracting ionosphere electric field in low latitude region based on incoherent scattering radar

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2013157762A (en) * 2013-12-25 2015-06-27 Федеральное государственное бюджетное учреждение "Институт прикладной геофизики имени академика Е.К.Федорова" (ФГБУ "ИПГ") METHOD FOR DETERMINING ELECTRONIC CONCENTRATION IN A SPECIFIED IONOSPHERE AREA AND DEVICE FOR ITS IMPLEMENTATION
CN107255802A (en) * 2017-07-21 2017-10-17 西安电子科技大学 Method based on the peak structure inverting strong disturbance ionosphere parameter of ion line spectrum three
CN111257867A (en) * 2020-01-20 2020-06-09 中国科学院地质与地球物理研究所 Ionized layer detection system based on magnetized plasma incoherent scattering theory spectrum
CN111273267A (en) * 2020-02-19 2020-06-12 中国科学院地质与地球物理研究所 Signal processing method, system and device based on phased array incoherent scattering radar

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11557079B2 (en) * 2015-02-25 2023-01-17 Atmospheric and Space Technology Research Associates, LLC Iterative ray-tracing for autoscaling of oblique ionograms

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2013157762A (en) * 2013-12-25 2015-06-27 Федеральное государственное бюджетное учреждение "Институт прикладной геофизики имени академика Е.К.Федорова" (ФГБУ "ИПГ") METHOD FOR DETERMINING ELECTRONIC CONCENTRATION IN A SPECIFIED IONOSPHERE AREA AND DEVICE FOR ITS IMPLEMENTATION
CN107255802A (en) * 2017-07-21 2017-10-17 西安电子科技大学 Method based on the peak structure inverting strong disturbance ionosphere parameter of ion line spectrum three
CN111257867A (en) * 2020-01-20 2020-06-09 中国科学院地质与地球物理研究所 Ionized layer detection system based on magnetized plasma incoherent scattering theory spectrum
CN111273267A (en) * 2020-02-19 2020-06-12 中国科学院地质与地球物理研究所 Signal processing method, system and device based on phased array incoherent scattering radar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Artem Miroshnikov 等.Determination of the Incoherent Scatter Signal Statistical Characteristics.2019,全文. *
梁宇 等.磁化等离子体非相干散射理论谱中的共振线研究.2020,第63卷(第63期),全文. *

Also Published As

Publication number Publication date
CN115656946A (en) 2023-01-31

Similar Documents

Publication Publication Date Title
CN115656946B (en) Incoherent scattering radar multi-beam calibration method and system for plasma line measurement
CN104515909B (en) A kind of large antenna pattern measurement method based on correlation method
US8121222B2 (en) Systems and methods for construction of time-frequency surfaces and detection of signals
Perrenoud The computer-controlled solar radio spectrometer ‘IKARUS’
CN111273267B (en) Signal processing method, system and device based on phased array incoherent scattering radar
RU2436130C2 (en) Method and system for radar probing earth interior
CN114428225B (en) Multi-radiation source arrival angle measuring method and device based on quasi-matched filtering
EP3399334A1 (en) Object detecting device and sensor device
Yan et al. GPS-R L1 interference signal processing for soil moisture estimation: an experimental study
US20090009382A1 (en) Method for determining signal direction using artificial doppler shifts
Toups et al. Analysis of foliage-induced synthetic pattern distortions
CN111190135A (en) Calibration method suitable for any array
JP5562237B2 (en) Method and apparatus for monitoring storm activity on the surface of the earth in real time
CN113567916B (en) Correlation interferometer direction finding method, system, equipment and storage medium
RU2294546C2 (en) Method for identification of radio-radiation
Chen et al. Examination of EMC chamber qualification methodology for applications above 1 GHz using frequency domain mode filtering
RU2309425C2 (en) Method of forming calibration data for radio direction finder/ range finder (versions)
Chen et al. Extending site VSWR to millimeter wave using cylindrical mode filtering
CN114325700A (en) Satellite-borne multi-channel SAR moving target imaging method
Zhang et al. The influence of the terrain on height measurement using the GNSS interference signal
RU2470315C1 (en) Method for computer-interferometer detection-direction finding of radio signals with expanded spectrum
US10871552B1 (en) Interstitial filter system and method with improved performance
Voulgaris Guidelines for assessing HF radar capabilities and performance
GB2279456A (en) Sonar direction finding
CN115754968A (en) Data processing method of radar system and radar system

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