CN113985414A - Distributed target parameter estimation and signal processing method based on dual-polarization radar - Google Patents

Distributed target parameter estimation and signal processing method based on dual-polarization radar Download PDF

Info

Publication number
CN113985414A
CN113985414A CN202111274305.8A CN202111274305A CN113985414A CN 113985414 A CN113985414 A CN 113985414A CN 202111274305 A CN202111274305 A CN 202111274305A CN 113985414 A CN113985414 A CN 113985414A
Authority
CN
China
Prior art keywords
intensity
reflected wave
wave
polarized wave
polarized
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.)
Pending
Application number
CN202111274305.8A
Other languages
Chinese (zh)
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.)
UNIT 63892 OF PLA
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN202111274305.8A priority Critical patent/CN113985414A/en
Publication of CN113985414A publication Critical patent/CN113985414A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/024—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a distributed target parameter estimation and signal processing method based on a dual-polarization radar, which is implemented based on the dual-polarization radar equipment and comprises the following steps: when transmitting/receiving equipment of the dual-polarized radar equipment receives horizontal polarized waves or vertical polarized waves, an antenna transmits the horizontal polarized waves or the vertical polarized waves to the sky; the receiver receives the horizontal polarized wave or the vertical polarized wave received by the antenna, converts an analog receiving signal of the horizontal polarized wave or the vertical polarized wave into a digital receiving signal, and outputs the digital receiving signal to the signal processing device; the signal processing device calculates the intensity of the first reflected wave, the intensity of the second reflected wave and the intensity of the third reflected wave according to the received signals; calculating a linear depolarization ratio according to a radar reflectivity factor obtained by integrating the intensity of the reflected wave; the Doppler velocity of the target is calculated from the reflected wave intensity. The present invention estimates rainfall intensity and/or doppler velocity by using horizontally polarized waves and vertically polarized waves.

Description

Distributed target parameter estimation and signal processing method based on dual-polarization radar
Technical Field
The invention relates to the technical field of radar signal processing, in particular to a distributed target parameter estimation and signal processing method based on a dual-polarization radar.
Background
In a unipolar weather radar which is widely applied, the rainfall intensity can be estimated through the echo intensity by utilizing a relational expression between the echo intensity and the rainfall intensity. However, since the relationship is greatly related to the particle size distribution of raindrops, it is difficult for the monopole wave weather radar to accurately estimate the rainfall intensity.
Generally, the raindrops are approximately spherical in shape, which is like a sphere that is squeezed in a vertical direction. Raindrops with large particle sizes have large flatness, whereas raindrops with small particle sizes have smaller flatness and are almost spherical. Therefore, when dual polarization radar observation is performed using a horizontally polarized wave and a vertically polarized wave, information corresponding to the distribution of the droplet particle size can be obtained from the droplet shape information, and thus the rainfall intensity can be measured with high accuracy.
In the transmission and reception of a dual-polarized radar polarized wave, the Linear Depolarization Ratio (LDR) refers to the ratio of the backscatter power of the same polarized wave (co-polarized scatter power) to the backscatter power generated by another polarized wave perpendicular to the polarized wave (cross-polarized scatter power).
The linear depolarization ratio is used as an index for qualitatively determining the shape and/or direction of ice particles when a polarized wave perpendicular to the transmitted polarized wave is scattered (Z) in precipitation particleshv,Zvh) The linear depolarization ratio becomes an effective parameter when it is relatively strong in the bright band and the like.
When observing the linear depolarization ratio using the alternating dual-polarization radar, the following three types of polarized wave transmission/reception processing units need to be repeatedly operated:
(1) a "horizontally polarized transmission and horizontally polarized reception" unit that mainly performs transmission/reception processing for transmitting a horizontally polarized wave and receiving the horizontally polarized wave;
(2) a "vertical polarization transmitting and vertical polarization receiving" unit that mainly performs transmission/reception processing for transmitting and receiving vertical polarization waves;
(3) a "horizontally polarized transmission and vertically polarized reception" unit that mainly performs transmission/reception processing for transmitting a horizontally polarized wave and receiving a vertically polarized wave, or a "transmission of a vertically polarized wave and reception of a horizontally polarized wave" unit that mainly performs transmission/reception processing for transmitting a vertically polarized wave and receiving a horizontally polarized wave.
Since the radar operates in alternate and alternate transmission with different polarizations, in the case of repeatedly operating the above three types of polarized wave transmission/reception processing units, not only the transmission of a single polarized wave and the reception of a single polarized wave but also the transmission of a plurality of polarized waves and the reception of a plurality of polarized waves are simultaneously performed, and at least three pulse repetition cycles are required to complete the measurement, which results in a triple increase in the observation time of the radar.
Disclosure of Invention
To solve the above problems, it is an object of the present invention to provide a distributed target parameter estimation and signal processing method based on a dual polarization radar, which estimates rainfall intensity and/or doppler velocity by using a horizontally polarized wave and a vertically polarized wave.
In order to achieve the purpose, the invention adopts the following technical scheme:
a first dual polarization radar apparatus including a transmitting/receiving device, a signal processing device and a display device, the signal processing device being composed of a reflected wave intensity calculator, a linear depolarization ratio calculator and a velocity calculator; wherein the content of the first and second substances,
the transmitting/receiving device is used for repeatedly executing the following operations in sequence: a first-time transmission/reception processing unit to transmit a horizontally polarized wave and receive the horizontally polarized wave; a second transmission/reception processing unit to transmit the vertical polarized wave and receive the vertical polarized wave; and a third-time transmission/reception processing unit for transmitting the horizontally polarized wave and receiving the vertically polarized wave;
the reflected wave intensity calculator comprises a signal input processing circuit and a reflected wave intensity calculating circuit, wherein the signal input processing circuit is used for receiving a digital receiving signal output from the transmitting/receiving equipment;
a reflected wave intensity calculation circuit for sequentially performing the following operations: calculating the intensity V of the first reflected wave from the digital received signal of the horizontally polarized wavehhThe reflected wave intensity is a power value of the horizontal polarized wave received by the first time transmission/reception processing unit in the transmission/reception apparatus; calculating the intensity V of the second reflected wave from the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit in the transmission/reception apparatus; and calculating the intensity V of the third reflected wave from the digital received signal of the vertical polarized wavevhThe reflected wave intensity is a power value of the vertical polarized wave received by the third transmission/reception processing unit in the transmission/reception apparatus; by measuring the intensity V of the reflected wavehh、VvvAnd VvhIntegrating, calculating radar reflectivity factor Z in transmitting and receiving horizontally polarized waveshhRadar reflectivity coefficient Z in vertical polarized wave transmission and vertical polarized wave receptionvvAnd in transmission of horizontally polarized waves and reception of vertically polarized wavesRadar reflectivity factor Zvh(ii) a A radar reflectivity factor difference Zdr is also calculated, which is a radar reflectivity factor Z of ZdrhhAnd radar reflectivity factor ZvvThe ratio of (A) to (B); calculating radar reflectivity factor Zhh、ZvvAnd ZvhAnd transmitting the calculated radar reflectivity factor difference Zdr to a display device for display;
the linear depolarization ratio calculator comprises a linear depolarization ratio calculation circuit for calculating a linear depolarization ratio LDRvhThe ratio is the first reflected wave intensity V calculated by the reflected wave intensity calculatorhhAnd third reflected wave intensity VvhThe ratio between, i.e. the radar reflectivity factor ZhhAnd radar reflectivity factor ZvhThe ratio therebetween; calculating linear depolarization ratio LDRvhTransmitting the data to a display device for display;
the velocity calculator includes a velocity calculating circuit for calculating the intensity V of the first reflected wave from the first reflected wave using a pulse pair processing methodhhIntensity of the second reflected wave VvvAnd third reflected wave intensity VvhIn (1), calculating the Doppler velocity of the target
Figure BDA0003328888450000031
Will calculate the Doppler velocity
Figure BDA0003328888450000032
Transmitted to the display device 7 for display;
and the display device is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device.
Further, the above-mentioned velocity calculator calculates the autocorrelation lag of each of the transmitted polarized waves from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity at time intervals; and calculating a phase difference between the polarized waves repeatedly transmitted by the transmitting/receiving device based on the lag; and calculating the Doppler velocity, the phase difference between polarized waves and a Nyquist rate of the target according to the lag.
A second dual polarization radar apparatus including a transmitting/receiving device, a signal processing device and a display device, the signal processing device being composed of a reflected wave intensity calculator, a linear depolarization ratio calculator and a velocity calculator; wherein the content of the first and second substances,
the transmitting/receiving device is used for repeatedly executing the following operations in sequence: a first-time transmission/reception processing unit to transmit a horizontally polarized wave and receive the horizontally polarized wave; a second transmission/reception processing unit to transmit the vertical polarized wave and receive the vertical polarized wave; and a third-time transmission/reception processing unit for transmitting the vertically polarized wave and receiving the horizontally polarized wave;
the reflected wave intensity calculator comprises a signal input processing circuit and a reflected wave intensity calculating circuit, wherein the signal input processing circuit is used for receiving a digital receiving signal output from the transmitting/receiving equipment;
a reflected wave intensity calculation circuit for sequentially performing the following operations: calculating the intensity V of the first reflected wave from the digital received signal of the horizontally polarized wavehhThe reflected wave intensity is a power value of the horizontal polarized wave received by the first time transmission/reception processing unit in the transmission/reception apparatus; calculating the intensity V of the second reflected wave from the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit in the transmission/reception apparatus; and calculating the intensity V of the third reflected wave from the digital received signal of the horizontally polarized wavehvThe reflected wave intensity is a power value of the horizontal polarized wave received by the third transmission/reception processing unit in the transmission/reception apparatus; by measuring the intensity V of the reflected wavehh、VvvAnd VhvIntegrating, calculating radar reflectivity factor Z in transmitting and receiving horizontally polarized waveshhRadar reflectivity coefficient Z in vertical polarized wave transmission and vertical polarized wave receptionvvAnd radar reflectivity factor Z in transmission of vertically polarized waves and reception of horizontally polarized waveshv(ii) a A radar reflectivity factor difference Zdr is also calculated, the radar reflectivity factor difference ZdrI.e. radar reflectivity factor ZhhAnd radar reflectivity factor ZvvThe ratio of (A) to (B); calculating radar reflectivity factor Zhh、ZvvAnd ZhvAnd transmitting the calculated radar reflectivity factor difference Zdr to a display device for display;
the linear depolarization ratio calculator is used for calculating a linear depolarization ratio LDRhvA ratio between the second reflected wave intensity and the third reflected wave intensity; the ratio is the second reflected wave intensity V calculated by the reflected wave intensity calculator 11vvAnd third reflected wave intensity VhvThe ratio between, i.e. the radar reflectivity factor ZvvAnd radar reflectivity factor ZhvThe ratio therebetween; calculating linear depolarization ratio LDRhvTransmitting the data to a display device for display;
the velocity calculator calculates a Doppler velocity of the target from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity using a pulse pair processing method
Figure BDA0003328888450000051
The display device is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device
Figure BDA0003328888450000052
Further, the above-mentioned velocity calculator calculates the autocorrelation lag of each of the transmitted polarized waves from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity at time intervals; and calculating a phase difference between the polarized waves repeatedly transmitted by the transmitting/receiving device based on the lag; and calculating the Doppler velocity, the phase difference between polarized waves and a Nyquist rate of the target according to the lag.
Further, the above-mentioned transmission/reception device is composed of a transmitter, a transmission/reception switch, an antenna, and a receiver, wherein,
the transmitter is connected with the transmitting/receiving switch and is used for outputting a horizontal polarized wave or a vertical polarized wave to the transmitting/receiving switch;
the transmitting/receiving switch is connected with the transmitter, the antenna and the receiver and used for outputting the horizontal polarized wave or the vertical polarized wave input from the transmitter to the antenna and outputting the received horizontal polarized wave or the received vertical polarized wave to the receiver by the antenna;
and the antenna is connected with the transmitting/receiving switch and is used for transmitting the horizontal polarized wave or the vertical polarized wave output from the transmitting/receiving switch and simultaneously receiving the horizontal polarized wave or the vertical polarized wave reflected by the target and then returning the horizontal polarized wave or the vertical polarized wave to the receiver.
The receiver is connected with the transmitting/receiving switch and used for receiving the horizontal polarized wave or the vertical polarized wave received by the antenna, carrying out digital processing on the received signal, converting the analog received signal into a digital received signal and outputting the digital received signal to the signal processing device.
Further, the signal processing apparatus described above further includes a computer, and the operation programs executed by the reflected wave intensity calculator, the linear depolarization ratio calculator, and the velocity calculator are stored in a memory of the computer, and at least one processor of the computer is configured to execute the programs stored in the memory.
Further, the above-mentioned signal processing apparatus, which calculates the reflected wave intensity as the power value of the horizontal polarized wave and the vertical polarized wave by the reflected wave intensity calculator, calculates the linear depolarization ratio LDR from the reflected wave intensity by the linear depolarization ratio calculator, and calculates the doppler velocity of the target by the velocity calculator, based on the digital received signals of the horizontal polarized wave and the vertical polarized wave output from the transmitting/receiving device
Figure BDA0003328888450000061
Further, the display device is mainly composed of a liquid crystal display for displaying the linear depolarization ratio LDR and the doppler velocity calculated by the signal processing device
Figure BDA0003328888450000062
A distributed target parameter estimation and signal processing method based on dual-polarization radar is implemented based on the dual-polarization radar device and comprises the following steps:
ST1, when a transmitting/receiving switch of the transmitting/receiving equipment receives a horizontal polarized wave or a vertical polarized wave, the antenna transmits the horizontal polarized wave or the vertical polarized wave to the sky;
ST2, the receiver receives the horizontal polarized wave or the vertical polarized wave received by the antenna, converts an analog reception signal of the horizontal polarized wave or the vertical polarized wave into a digital reception signal, and outputs the digital reception signal to the signal processing device;
ST3, a reflected wave intensity calculator of the signal processing device calculates the intensity of the first reflected wave, the intensity of the second reflected wave, and the intensity of the third reflected wave based on the received signal;
ST4, calculating a linear depolarization ratio according to a radar reflectivity factor obtained by integrating the intensity of the reflected wave;
ST5, the Doppler velocity of the target is calculated from the intensity of the reflected wave.
Due to the adoption of the technical scheme, the invention has the following advantages:
according to the distributed target parameter estimation and signal processing method based on the dual-polarization radar, the linear depolarization ratio calculator is used for calculating the linear depolarization ratio through configuration, even under the condition that three types of polarized wave transmitting/receiving processing are repeatedly operated, the linear depolarization ratio can be calculated under the condition that the Nyquist rate is reduced, and the method is high in processing speed, short in time and high in accuracy.
Drawings
Fig. 1 is a block diagram of the construction of a dual polarized radar apparatus of the present invention;
fig. 2 is a block diagram of the signal processing apparatus in fig. 1;
FIG. 3 is a schematic block circuit diagram of the signal processing apparatus of FIG. 2;
FIG. 4 is a hardware block diagram of a computer included in the signal processing apparatus of FIG. 1;
FIG. 5 is a flow chart of a distributed target parameter estimation and signal processing method based on dual-polarization radar according to the present invention;
fig. 6 is a schematic diagram of an embodiment of a transmitting/receiving apparatus of the present invention that transmits and receives polarized waves;
FIG. 7 is
Figure BDA0003328888450000071
And
Figure BDA0003328888450000072
a phase relationship diagram therebetween;
FIG. 8 is a diagram in which each time interval for transmitting and receiving polarized waves by the transmitting/receiving apparatus of FIG. 1 has a length TsExamples of (2) are shown.
Detailed Description
The technical solution of the present invention will be further described in detail with reference to the accompanying drawings and examples.
As shown in fig. 1, 2, 3, a dual polarized radar apparatus comprises a transmitting/receiving apparatus 1, a signal processing means 6 and a display means 7, wherein,
the transmitting/receiving equipment 1 consists of a transmitter 2, a transmitting/receiving switch 3, an antenna 4 and a receiver 5, wherein the transmitter is connected with the transmitting/receiving switch and used for outputting a horizontal polarized wave or a vertical polarized wave to the transmitting/receiving switch; the transmitting/receiving switch is connected with the transmitter, the antenna and the receiver and used for outputting the horizontal polarized wave or the vertical polarized wave input from the transmitter to the antenna and outputting the received horizontal polarized wave or the received vertical polarized wave to the receiver by the antenna; the antenna is connected with the transmitting/receiving switch and is used for transmitting the horizontal polarized wave or the vertical polarized wave output by the transmitting/receiving switch and simultaneously receiving the horizontal polarized wave or the vertical polarized wave reflected by a target (such as raindrops existing in the sky) and then returning the horizontal polarized wave or the vertical polarized wave to the receiver; the receiver is connected with the transmitting/receiving switch and used for receiving the horizontal polarized wave or the vertical polarized wave received by the antenna, carrying out digital processing on the received signal, converting the analog received signal into a digital received signal and outputting the digital received signal to the signal processing device;
the signal processing device 6 consists of a reflected wave intensity calculator 11, a linear depolarization ratio calculator 12 and a speed calculator 13; wherein the content of the first and second substances,
the transmitting/receiving device is used for repeatedly executing the following operations in sequence: a first-time transmission/reception processing unit to transmit a horizontally polarized wave and receive the horizontally polarized wave; a second transmission/reception processing unit to transmit the vertical polarized wave and receive the vertical polarized wave; and a third-time transmission/reception processing unit for transmitting the horizontally polarized wave and receiving the vertically polarized wave;
the reflected wave intensity calculator includes a signal input processing circuit 21 and a reflected wave intensity calculating circuit 22, the signal input processing circuit 21 receiving the digital received signal output from the transmission/reception apparatus 1;
a reflected wave intensity calculation circuit for sequentially performing the following operations: calculating the intensity V of the first reflected wave from the digital received signal of the horizontally polarized wavehhThe reflected wave intensity is a power value of the horizontal polarized wave received by the first time transmission/reception processing unit in the transmission/reception apparatus; calculating the intensity V of the second reflected wave from the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit in the transmission/reception apparatus; and calculating the intensity V of the third reflected wave from the digital received signal of the vertical polarized wavevhThe reflected wave intensity is a power value of the vertical polarized wave received by the third transmission/reception processing unit in the transmission/reception apparatus; by measuring the intensity V of the reflected wavehh、VvvAnd VvhIntegrating, calculating radar reflectivity factor Z in transmitting and receiving horizontally polarized waveshhRadar reflectivity coefficient Z in vertical polarized wave transmission and vertical polarized wave receptionvvAnd radar reflectivity factor Z in transmission of horizontally polarized waves and reception of vertically polarized wavesvh(ii) a A radar reflectivity factor difference Zdr is also calculated, the radar reflectivity factor difference Zdr being the radarUp to the reflectance factor ZhhAnd radar reflectivity factor ZvvThe ratio of (A) to (B); calculating radar reflectivity factor Zhh、ZvvAnd ZvhAnd transmitting the calculated radar reflectivity factor difference Zdr to a display device for display;
the linear depolarization ratio calculator includes a linear depolarization ratio calculation circuit 23 for calculating a linear depolarization ratio LDRvhThe ratio is the first reflected wave intensity V calculated by the reflected wave intensity calculator 11hhAnd third reflected wave intensity VvhThe ratio between, i.e. the radar reflectivity factor ZhhAnd radar reflectivity factor ZvhThe ratio therebetween; calculating linear depolarization ratio LDRvhTransmitted to the display device 7 for display;
the velocity calculator includes a velocity calculating circuit 24 for reflecting the wave intensity V from the first time using a Pulse Pair Processing (PPP) methodhhIntensity of the second reflected wave VvvAnd third reflected wave intensity VvhIn (1), calculating the Doppler velocity of the target
Figure BDA0003328888450000091
Will calculate the Doppler velocity
Figure BDA0003328888450000092
Transmitted to the display device 7 for display;
and the display device is connected with the signal processing device and is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device.
As shown in fig. 4, the signal processing apparatus 6 further includes a computer, and the operation programs executed by the reflected wave intensity calculator 11, the linear depolarization ratio calculator 12, and the velocity calculator 13 are stored in a memory 31 of the computer, and at least one processor 32 of the computer is configured to execute the programs stored in the memory 31.
As shown in FIG. 6, each time interval of the transmission and reception of the polarized wave by the transmission/reception device 1 is given a length TsOfExamples are shown.
According to (3i-2) -TsPerforms transmission of a horizontally polarized wave and reception of the horizontally polarized wave as a first transmission/reception processing unit in accordance with (3i-1) -TsPerforms transmission of a vertically polarized wave and reception of the vertically polarized wave at intervals of time as a second transmission/reception processing unit, and in accordance with 3i-TsAs a third-time transmission/reception processing unit, the transmission of the horizontally polarized wave and the reception of the vertically polarized wave are performed at intervals of time; in this embodiment, I ═ l, 2, 3.
I denotes the number of sample pairs, i.e. the number of samples in each of the following processes: transmission of a horizontally polarized wave and reception of a horizontally polarized wave, transmission of a vertically polarized wave and reception of a vertically polarized wave, and transmission of a horizontally polarized wave and reception of a vertically polarized wave.
Setting the transmission/reception device 1, the following operations are repeatedly performed in order: transmitting a horizontally polarized wave and receiving the horizontally polarized wave as a first transmission/reception processing unit; transmitting a vertically polarized wave and receiving the vertically polarized wave as a second transmission/reception processing unit; the transmission of the horizontally polarized wave and the reception of the vertically polarized wave serve as a third transmission/reception processing unit.
An example in which the dual polarization radar apparatus is targeted to raindrops will be described in detail with reference to fig. 6.
Since the transmitter 2 in the transmission/reception apparatus 1 executes the transmission/reception processing units in the following order: transmitting horizontally polarized wave and receiving horizontally polarized wave → transmitting vertically polarized wave and receiving vertically polarized wave, transmitting transmitter with TsTransmitting a horizontally polarized wave → a vertically polarized wave → a horizontally polarized wave → … → a horizontally polarized wave → a vertically polarized wave → a horizontally polarized wave for a transmission time interval.
When receiving a horizontally polarized wave or a vertically polarized wave from the transmitter 2, the transmission/reception switch 3 outputs the horizontally polarized wave or the vertically polarized wave to the antenna 4.
When receiving a horizontally polarized wave or a vertically polarized wave from the transmission/reception switch 3, the antenna 4 transmits the horizontally polarized wave or the vertically polarized wave to the sky, see step ST1 in fig. 5.
The antenna 4 transmits a horizontally polarized wave or a vertically polarized wave to the sky, receives the horizontally polarized wave or the vertically polarized wave reflected by the raindrops as a target, and then returns the horizontally polarized wave or the vertically polarized wave.
Antenna with TsThe reception time interval of (1) receives a horizontally polarized wave → a vertically polarized wave → a horizontally polarized wave → a vertically polarized wave.
The transmission/reception switch 3 outputs the horizontally polarized wave or the vertically polarized wave received by the antenna 4 to the receiver 5.
The receiver 5 receives the horizontally polarized wave or the vertically polarized wave received by the antenna 4, converts an analog reception signal of the horizontally polarized wave or the vertically polarized wave into a digital reception signal, and outputs the digital reception signal to the signal processing device 6, see step ST2 in fig. 5.
A reflected wave intensity calculator 11 of the signal processing device 6 that calculates a reflected wave intensity Vrin, which is a power value of the horizontally polarized wave received by the first transmission/reception processing unit in the transmission/reception apparatus 1, from the digital received signal of the horizontally polarized wave; and simultaneously calculating the reflected wave intensity V according to the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit; calculating reflected wave intensity V from digital received signal of vertical polarized wavevhThe reflected wave intensity is the power value of the vertical polarized wave received by the third transmission/reception processing unit, see step ST3 in fig. 5.
The calculation of the reflected wave intensity V performed by the reflected wave intensity calculator 11 will be specifically described belowhh、VvvAnd VvhAnd (4) processing.
Since raindrops as an observation object generally exist randomly in the sky, the digital reception signal output from the transmission/reception device 1 also has randomness.
Therefore, the reflected wave intensity V is calculatedhh、VvvAnd VvhAt this time, the reflected wave intensity calculator 11 performs averaging (digital integration) on digital received signals (statistically independent received signals) obtained by transmission of polarized waves (multiple transmission of pulses), thereby suppressing fluctuation of the estimated value of each reflected wave intensity and improving the accuracy of each reflected wave intensity.
Digital received signals as the above-mentioned averaging processing targets, which include: a digital received signal of a horizontally polarized wave received a plurality of times by the first transmission/reception process, a digital received signal of a vertically polarized wave received a plurality of times by the second transmission/reception process, and a digital received signal of a vertically polarized wave received a plurality of times by the third transmission/reception process.
In this regard, the relationship among the voltage matrix of the backscattered reception signal (i.e., the reception signal of the antenna 4), the transmission matrix of the polarized wave, and the backscatter matrix of the polarized wave is generally given by the following equations (1) to (3) as relational expressions:
Figure BDA0003328888450000111
Figure BDA0003328888450000112
Figure BDA0003328888450000113
in the formula: vhA voltage that is a horizontally polarized wave backscatter receive signal;
Vva voltage that is a vertical polarized wave backscatter receive signal;
Shha scattering parameter in "transmission of a horizontally polarized wave and reception of a horizontally polarized wave";
Svvfor transmission of "vertically polarized wavesAnd scattering parameters in reception of vertically polarized waves ";
Svhscattering parameters in "transmission of horizontally polarized wave and reception of vertically polarized wave";
Shvscattering parameters in "transmission of vertically polarized wave and reception of horizontally polarized wave";
λ is the wavelength of the polarized wave;
g is the gain of the antenna 4;
r is a distance from the dual-polarization radar device to a raindrop as an observation object;
Mhis a parameter which is set to 1 when a horizontally polarized wave is transmitted, while being set to 0 when a vertically polarized wave is transmitted;
Mvis a parameter which is set to 0 when a horizontally polarized wave is transmitted and to 1 when a vertically polarized wave is transmitted;
K0is a propagation constant in free space;
Khis a space propagation constant which depends on the horizontally polarized wave;
kvis a spatial propagation constant that depends on the vertically polarized wave.
The reflected wave intensity calculator 11 calculates the sum of the voltages of the first backscattered received signal, i.e., the digital received signal of the horizontally polarized wave received a plurality of times by the first transmission/reception process, i.e., the reflected wave intensity V in the nth samplehh(n) is as shown in the following formula (4).
The reflected wave intensity calculator 11 also calculates the sum of the voltages of the second backscattered received signal, i.e., the digital received signal of the vertically polarized wave received a plurality of times by the second transmission/reception process, i.e., the reflected wave intensity V in the (n + l) -th samplevv(n +1) as shown in the following formula (5).
The reflected wave intensity calculator 11 further calculates the sum of voltages of a third backscattered received signal, i.e., a digital received signal of a vertically polarized wave received a plurality of times by the third transmission/reception process, i.e., the fourthReflected wave intensity V in (n +2) samplesvh(n +2) as shown in the following formula (6). These reflected wave intensities are calculated for each segment (grid point) in azimuth for a certain observation region and range.
Figure BDA0003328888450000121
Figure BDA0003328888450000122
Figure BDA0003328888450000123
Wherein i is an index of the scattered wave at each resolution unit (i-th);
Iiis a matrix containing the initial phase.
Due to the distance r1(n + l) is obtained by reacting TsviA distance r from the previous onei(n) are added to move a point at the distance by one distance, and thus, the distance r is expressed using the following formula (7)i(n+l)。
ri(n+1)＝ri(n)+Tsvi (7)
After the reflected wave intensity V is calculatedhh(n)、Vvv(n +1) and Vvh(n +2) thereafter, the reflected wave intensity calculator 11 measures the reflected wave intensity Vhh(n)、Vvv(n +1) and Vvh(n +2) performing integration; thereby determining radar reflectivity factors Z corresponding to the transmission and reception of the horizontally polarized waveshhThe radar reflectivity factor is the intensity of the reflected wave Vhh(n) the intensity of the reflected wave after integration; radar reflectivity factor Z corresponding to transmission and reception of vertically polarized wavesvvThe radar reflectivity factor is the intensity of the reflected wave Vvv(n +1) the intensity of the reflected wave after integration; and transmission of horizontally polarized waves and vertically polarized wavesCorresponding radar reflectivity factor Z in receptionvhThe radar reflectivity factor is the intensity of the reflected wave Vvh(n +2) the integrated reflected wave intensity.
The user can arbitrarily determine the range of the integration in consideration of the accuracy of the radar apparatus and the like. For example, it is assumed that the integration range is the length T of each time interval for transmitting and receiving polarized waves by the transmission/reception device 1sAnd the known beam range (range of radar rotation) provided by the transmitting/receiving device 1.
The reflected wave intensity calculator 11 also calculates a radar reflectivity factor difference Zdr, which is a radar reflectivity factor ZhhAnd radar reflectivity factor ZvvThe ratio of (c) is shown in the following equation (8).
Figure BDA0003328888450000131
The reflected wave intensity calculator 11 further calculates a radar reflectivity factor ZhhAnd ZvhOutputs to the linear depolarization ratio calculator 12, and also outputs the radar reflectivity factor Zhh、ZvvAnd ZvhAnd the radar reflectance factor difference Zdr are displayed on the display device 7.
When receiving the radar reflectivity factor Z from the reflected wave intensity calculator 11hhAnd ZvhThereafter, the linear depolarization ratio calculator 12 calculates the radar reflectivity factor ZhhAnd ZvhSubstituted into the following equation (9) to calculate the linear depolarization ratio LDRvhSee step ST4 in fig. 5.
Figure BDA0003328888450000141
In calculating linear depolarization ratio LDRvhThereafter, the linear depolarization ratio calculator 12 displays the linear depolarization rate LDR on the display device 7vh。
When the reflected wave intensity calculator 11 calculates the reflected wave intensity Vhh、VvvAnd VvhThe velocity calculator 13 uses a pulse pair method based on the reflected wave intensity Vhh、VvvAnd VvhCalculating the Doppler velocity of a target
Figure BDA0003328888450000142
See step ST5 in fig. 5.
Specifically, the velocity calculator 13 calculates the reflected wave intensity V from the reflected wave intensity calculator 11hh、VvvAnd VvhRepeatedly transmitting T of polarized wave according to the transmitting/receiving apparatus 1sTime intervals, calculating the autocorrelation lag of each transmitted polarized wave
Figure BDA0003328888450000143
According to the lag
Figure BDA0003328888450000144
Calculating the phase difference between polarized waves repeatedly transmitted by the transmitting/receiving apparatus 1
Figure BDA0003328888450000145
And according to
Figure BDA0003328888450000146
Calculating the Doppler velocity of a target
Figure BDA0003328888450000147
Phase difference between polarized waves
Figure BDA0003328888450000148
And nyquist rate vaThe velocity represents the Doppler velocity allowed to be observed
Figure BDA0003328888450000149
Range).
Hereinafter, the doppler velocity performed by the velocity calculator 13 will be specifically described
Figure BDA00033288884500001410
The calculation process of (2).
In an autocorrelation estimation method (pulse pair processing method) in polarized wave alternate transmission, hysteresis is used
Figure BDA00033288884500001411
Represented by the following formula (10). In general terms, the amount of the solvent to be used,
Figure BDA00033288884500001412
referred to as a "lag 1",
Figure BDA00033288884500001413
known as "lag 2"
Figure BDA00033288884500001414
In equation (10), M represents the number of sample pairs, i.e., the number of samplings in each of the following transmit/receive processing units: transmission of a horizontally polarized wave and reception of a horizontally polarized wave, transmission of a vertically polarized wave and reception of a vertically polarized wave, and transmission of a horizontally polarized wave and reception of a vertically polarized wave.
Shown in formula (10)
Figure BDA00033288884500001415
Is decomposed into three formulas as shown in the following formulas (11) to (13). Specifically, as shown in the formula (10)
Figure BDA00033288884500001416
Can be made of
Figure BDA00033288884500001417
And
Figure BDA00033288884500001418
to be specified.
Figure BDA0003328888450000151
Figure BDA0003328888450000152
Figure BDA0003328888450000153
Figure BDA0003328888450000154
Represents the covariance of the horizontally polarized wave when the horizontally polarized wave is received, i.e., the covariance between the horizontally polarized wave received by the first transmission/reception process and the vertically polarized wave received by the second transmission/reception process.
Figure BDA0003328888450000155
The covariance of the vertically polarized wave, that is, the covariance between the vertically polarized wave received through the second transmission/reception process and the vertically polarized wave received through the third transmission/reception process is expressed.
Figure BDA0003328888450000156
Represents the covariance of the horizontally polarized wave when the vertically polarized wave is received, i.e., the covariance between the vertically polarized wave received by the third transmission/reception process and the horizontally polarized wave received by the first transmission/reception process.
Thus, as shown in equation (10)
Figure BDA0003328888450000157
Represented by the following formula (14).
Figure BDA0003328888450000158
By using the reflected wave intensity V in the formula (4)hh(n) andreflected wave intensity V in equation (5)vv(n +1) as shown in the formula (11)
Figure BDA0003328888450000159
Represented by the following formula (15).
Figure BDA00033288884500001510
Due to 2 (k) of phase0+kh)ri(n) is uniformly distributed in the range of 0 to 2 pi, so that the expectation value of the exponent part in the formula (15) is zero in the case other than i ═ l.
Thus, as shown in the formula (15)
Figure BDA00033288884500001511
Represented by the following formula (16).
Figure BDA0003328888450000161
Because r in the pair formula (16)i(n) when averaging, the average is the midpoint of all distances, ri(n) may be set to r0(intermediate distance). Furthermore, because in the pair viWhen averaging, the average is one of all speeds, so viMay be set to v (average velocity).
Figure BDA0003328888450000162
Because 2 (k)h-kv)r0Equal to the phase difference between the horizontally polarized wave and the vertically polarized wave in equation (17)
Figure BDA00033288884500001636
Therefore, it is not only easy to use
Figure BDA00033288884500001637
Furthermore, because of-2 (k)0+kv)Tsv is equal to the doppler shift amount ψ d of the vertically polarized wave, so-2 (k0+ kv) Tsv ═ ψ d.
FIG. 7 shows
Figure BDA0003328888450000163
And
Figure BDA0003328888450000164
the phase relationship between them, in the example given in fig. 7, the pair (T) is omitted for simplification of the drawings) And will be described
Figure BDA0003328888450000165
And
Figure BDA0003328888450000166
are respectively represented as
Figure BDA0003328888450000167
And
Figure BDA0003328888450000168
as shown in figure 7 of the drawings,
Figure BDA0003328888450000169
and
Figure BDA00033288884500001610
or
Figure BDA00033288884500001611
A phase difference therebetween of
Figure BDA00033288884500001638
And is
Figure BDA00033288884500001612
Or
Figure BDA00033288884500001613
And
Figure BDA00033288884500001614
or
Figure BDA00033288884500001615
A phase difference therebetween of
Figure BDA00033288884500001616
In addition to this, the present invention is,
Figure BDA00033288884500001617
or
Figure BDA00033288884500001618
Has a phase of psid。
Therefore, the phase difference between the polarized waves repeatedly transmitted by the transmission/reception device 1
Figure BDA00033288884500001619
By passing
Figure BDA00033288884500001620
And
Figure BDA00033288884500001621
two thirds of the parameter in between, or
Figure BDA00033288884500001622
And
Figure BDA00033288884500001623
two thirds of the parameter in between.
Specifically, the phase difference between polarized waves repeatedly transmitted by the transmission/reception device 1
Figure BDA00033288884500001624
Can be expressed by the following formula.
Figure BDA00033288884500001625
For this reason, as a phase difference between polarized waves repeatedly transmitted by the transmission/reception device 1
Figure BDA00033288884500001626
An example of the calculation by the velocity calculator 13 is given, in which the phase difference is
Figure BDA00033288884500001627
And
Figure BDA00033288884500001628
two thirds of the parameter in between, or
Figure BDA00033288884500001629
And
Figure BDA00033288884500001630
two thirds of the parameter in between. On the other hand, the velocity calculator calculates
Figure BDA00033288884500001631
And
Figure BDA00033288884500001632
an average of two-thirds of the parameters in between, and
Figure BDA00033288884500001633
and
Figure BDA00033288884500001634
an average value of two-thirds of the parameters as a phase difference between polarized waves repeatedly transmitted by the transmitting/receiving apparatus 1
Figure BDA00033288884500001635
Specifically as shown in the following formula
Figure BDA0003328888450000171
After calculating the phase difference between the polarized waves
Figure BDA0003328888450000172
Thereafter, the speed calculator 13 calculates the speed according to the hysteresis
Figure BDA0003328888450000173
Is contained in
Figure BDA0003328888450000174
Calculating the Doppler velocity of a target
Figure BDA0003328888450000175
Phase difference between polarized waves
Figure BDA0003328888450000176
And the Nyquist rate vaThe velocity indicating the Doppler velocity allowed to be observed
Figure BDA0003328888450000177
The range of (c) is shown in the following equation (20).
Figure BDA0003328888450000178
In the present embodiment, equation (20) and calculating Doppler velocity
Figure BDA0003328888450000179
Is the same (in the case of performing a process of transmitting and receiving a single polarized wave or a process of simultaneously transmitting and receiving a polarized wave), and the nyquist rate v is set to be equal toaThe velocity indicating the Doppler velocity allowed to be observed
Figure BDA00033288884500001710
The range of (a) is restored to the nyquist rate which is the same as that obtained in the process of transmitting a single polarized wave and receiving a single polarized wave or simultaneously transmitting a polarized wave and simultaneously receiving a polarized wave.
Therefore, even when three types of polarized wave transmission/reception processing units are repeatedly executed, the Doppler velocity of the target
Figure BDA00033288884500001711
It is still possible to perform the calculation using the nyquist rate which is the same as that obtained in the process of transmitting a single polarized wave and receiving a single polarized wave or simultaneously transmitting a polarized wave and simultaneously receiving a polarized wave.
As can be seen from the above description, the configuration of the embodiment includes the linear depolarization ratio calculator 12 for determining the radar reflectivity factor Z in the horizontally polarized wave transmission and the horizontally polarized receptionhh(the radar reflectance factor is a factor for the reflected wave intensity V calculated by the reflected wave intensity calculator 11hh(n) reflected wave intensity after integration), and radar reflectivity factor Z in transmission of horizontally polarized waves and reception of vertically polarized wavesvh(the radar reflectivity factor is the intensity of the reflected wave Vvh(n +2) intensity of reflected wave after integration), and calculating linear depolarization ratio LDRvhThe linear depolarization ratio is a radar reflectivity factor ZhhAnd radar reflectivity factor ZvhThe ratio of. Therefore, there is an advantage in that even when three types of polarized wave transmission/reception processing elements are repeatedly executed, it is possible to calculate the linear depolarization ratio LDR while preventing the nyquist rate va from decreasingvh。
In the present embodiment, the transmission/reception apparatus 1 gives "horizontally polarized transmission and horizontally polarized reception" as the first transmission/reception processing, "vertically polarized transmission and vertically polarized reception" as the second transmission/reception processing, and "horizontally polarized transmission and vertically polarized reception" as the third transmission/reception processing, respectively, in the example.
In this case, the reflected wave intensity calculator 11 does not calculate the sum of the voltages of the backscattered received signals corresponding to the vertically polarized wave digital received signals received a plurality of times by the third transmission/reception process, that is, the reflected wave intensity V in the (n +2) th samplevh(n + 2); but calculateThe voltage sum of the backscatter received signals corresponding to the horizontally polarized wave digital received signals received a plurality of times by the fourth transmission/reception process, i.e., the reflected wave intensity V in the (n +2) th samplehv(n+2)。
The reflected wave intensity calculator 11 also calculates the reflected wave intensity Vhv(n +2) is integrated to determine the radar reflectivity factor Z in the vertical polarization transmission and horizontal polarization receptionhvThe radar reflectivity factor being the intensity of the reflected wave Vhv(n +2) reflected wave intensity after integration.
Although the range of the integration can be arbitrarily determined by the user in consideration of the accuracy of the radar apparatus or the like, it is assumed here that the integration range is transmitted and received by the transmission/reception apparatus 1 for each time interval length T of the polarized wavesAnd a given known beam range (range of radar rotation) of the transmitting/receiving device 1.
In addition, the reflected wave intensity calculator 11 calculates a radar reflectance factor ZvvAnd ZhvOutputs to the linear depolarization ratio calculator 12, and also outputs the radar reflectivity factor Zhh，ZvvAnd ZhvAnd the radar reflectance factor difference Zdr are displayed on the display device 7.
When receiving the radar reflectivity factor Z from the reflected wave intensity calculator 11vvAnd ZhvThe linear depolarization ratio calculator 12 calculates the radar reflectivity factor ZvvAnd ZhvSubstituting into the following equation (21) to calculate the linear depolarization ratio LDRhv。
Figure BDA0003328888450000181
The linear depolarization ratio calculator 12 calculates a linear depolarization ratio LDRhvThereafter, the linear depolarization ratio LDR is displayed on a display device 7hv。
When the reflected wave intensity calculator 11 calculates the reflected wave intensity Vhh、VvvAnd VhhvThe velocity calculator 13 uses the reflected wave intensity V by using a pulse pair methodhh、VvvAnd VhhvCalculating the Doppler velocity of a target
Figure BDA0003328888450000197
Specifically, the velocity calculator 13 calculates the reflected wave intensity V from the reflected wave intensity calculator 11hh、VvvAnd VhhvCalculating the autocorrelation lag in each of the transmitted polarized waves at the time interval Ts during which the transmitting/receiving apparatus 1 repeatedly transmits the polarized waves
Figure BDA0003328888450000191
Then according to the lag
Figure BDA0003328888450000192
Calculating the phase difference between polarized waves repeatedly transmitted by the transmitting/receiving apparatus 1
Figure BDA0003328888450000193
And according to hysteresis
Figure BDA0003328888450000194
Calculating the Doppler velocity of a target
Figure BDA0003328888450000195
Phase difference between polarized waves
Figure BDA0003328888450000196
And nyquist rate va。
The above description is only a preferred embodiment of the present invention, and not intended to limit the present invention, and all equivalent changes and modifications made within the scope of the claims of the present invention should fall within the protection scope of the present invention.

Claims (9)

1. A dual-polarized radar apparatus, characterized by: the device comprises a transmitting/receiving device, a signal processing device and a display device, wherein the signal processing device consists of a reflected wave intensity calculator, a linear depolarization ratio calculator and a speed calculator; wherein the content of the first and second substances,
the transmitting/receiving device is used for repeatedly executing the following operations in sequence: a first-time transmission/reception processing unit to transmit a horizontally polarized wave and receive the horizontally polarized wave; a second transmission/reception processing unit to transmit the vertical polarized wave and receive the vertical polarized wave; and a third-time transmission/reception processing unit for transmitting the horizontally polarized wave and receiving the vertically polarized wave;
the reflected wave intensity calculator comprises a signal input processing circuit and a reflected wave intensity calculating circuit, wherein the signal input processing circuit is used for receiving a digital receiving signal output from the transmitting/receiving equipment;
a reflected wave intensity calculation circuit for sequentially performing the following operations: calculating the intensity V of the first reflected wave from the digital received signal of the horizontally polarized wavehhThe reflected wave intensity is a power value of the horizontal polarized wave received by the first time transmission/reception processing unit in the transmission/reception apparatus; calculating the intensity V of the second reflected wave from the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit in the transmission/reception apparatus; and calculating the intensity V of the third reflected wave from the digital received signal of the vertical polarized wavevhThe reflected wave intensity is a power value of the vertical polarized wave received by the third transmission/reception processing unit in the transmission/reception apparatus; by measuring the intensity V of the reflected wavehh、VvvAnd VvhIntegrating, calculating radar reflectivity factor Z in transmitting and receiving horizontally polarized waveshhRadar reflectivity coefficient Z in vertical polarized wave transmission and vertical polarized wave receptionvvAnd radar reflectivity factor Z in transmission of horizontally polarized waves and reception of vertically polarized wavesvh(ii) a A radar reflectivity factor difference Zdr is also calculated, which is a radar reflectivity factor Z of ZdrhhAnd radar reflectivity factor ZvvThe ratio of (A) to (B); calculating radar reflectivity factor Zhh、ZvvAnd ZvhAnd transmitting the calculated radar reflectivity factor difference Zdr to a display deviceCentering and displaying;
the linear depolarization ratio calculator comprises a linear depolarization ratio calculation circuit for calculating a linear depolarization ratio LDRvhThe ratio is the first reflected wave intensity V calculated by the reflected wave intensity calculatorhhAnd third reflected wave intensity VvhThe ratio between, i.e. the radar reflectivity factor ZhhAnd radar reflectivity factor ZvhThe ratio therebetween; calculating linear depolarization ratio LDRvhTransmitting the data to a display device for display;
the velocity calculator includes a velocity calculating circuit for calculating the intensity V of the first reflected wave from the first reflected wave using a pulse pair processing methodhhIntensity of the second reflected wave VvvAnd third reflected wave intensity VvhIn (1), calculating the Doppler velocity of the target
Figure FDA0003328888440000021
Will calculate the Doppler velocity
Figure FDA0003328888440000022
Transmitted to the display device 7 for display;
and the display device is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device.
2. The dual polarized radar apparatus of claim 1, wherein: a velocity calculator that calculates an autocorrelation lag of each of the transmitted polarized waves from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity at time intervals; and calculating a phase difference between the polarized waves repeatedly transmitted by the transmitting/receiving device based on the lag; and calculating the Doppler velocity, the phase difference between polarized waves and a Nyquist rate of the target according to the lag.
3. A dual-polarized radar apparatus, characterized by: the device comprises a transmitting/receiving device, a signal processing device and a display device, wherein the signal processing device consists of a reflected wave intensity calculator, a linear depolarization ratio calculator and a speed calculator; wherein the content of the first and second substances,
the transmitting/receiving device is used for repeatedly executing the following operations in sequence: a first-time transmission/reception processing unit to transmit a horizontally polarized wave and receive the horizontally polarized wave; a second transmission/reception processing unit to transmit the vertical polarized wave and receive the vertical polarized wave; and a third-time transmission/reception processing unit for transmitting the vertically polarized wave and receiving the horizontally polarized wave;
the reflected wave intensity calculator comprises a signal input processing circuit and a reflected wave intensity calculating circuit, wherein the signal input processing circuit is used for receiving a digital receiving signal output from the transmitting/receiving equipment;
a reflected wave intensity calculation circuit for sequentially performing the following operations: calculating the intensity V of the first reflected wave from the digital received signal of the horizontally polarized wavehhThe reflected wave intensity is a power value of the horizontal polarized wave received by the first time transmission/reception processing unit in the transmission/reception apparatus; calculating the intensity V of the second reflected wave from the digital received signal of the vertical polarized wavevvThe reflected wave intensity is a power value of the vertical polarized wave received by the second transmission/reception processing unit in the transmission/reception apparatus; and calculating the intensity V of the third reflected wave from the digital received signal of the horizontally polarized wavehvThe reflected wave intensity is a power value of the horizontal polarized wave received by the third transmission/reception processing unit in the transmission/reception apparatus; by measuring the intensity V of the reflected wavehh、VvvAnd VhvIntegrating, calculating radar reflectivity factor Z in transmitting and receiving horizontally polarized waveshhRadar reflectivity coefficient Z in vertical polarized wave transmission and vertical polarized wave receptionvvAnd radar reflectivity factor Z in transmission of vertically polarized waves and reception of horizontally polarized waveshv(ii) a A radar reflectivity factor difference Zdr is also calculated, which is a radar reflectivity factor Z of ZdrhhAnd radar reflectivity factor ZvvThe ratio of (A) to (B); calculating radar reflectivity factor Zhh、ZvvAnd ZhvAnd transmitting the calculated radar reflectivity factor difference Zdr to a display device for display;
the linear depolarization ratio calculator is used for calculating a linear depolarization ratio LDRhvA ratio between the second reflected wave intensity and the third reflected wave intensity; the ratio is the second reflected wave intensity V calculated by the reflected wave intensity calculator 11vvAnd third reflected wave intensity VhvThe ratio between, i.e. the radar reflectivity factor ZvvAnd radar reflectivity factor ZhvThe ratio therebetween; calculating linear depolarization ratio LDRhvTransmitting the data to a display device for display;
the velocity calculator calculates a Doppler velocity of the target from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity using a pulse pair processing method
Figure FDA0003328888440000031
The display device is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device
Figure FDA0003328888440000032
4. The dual polarized radar apparatus of claim 3, wherein: a velocity calculator that calculates an autocorrelation lag of each of the transmitted polarized waves from the first reflected wave intensity, the second reflected wave intensity, and the third reflected wave intensity at time intervals; and calculating a phase difference between the polarized waves repeatedly transmitted by the transmitting/receiving device based on the lag; and calculating the Doppler velocity, the phase difference between polarized waves and a Nyquist rate of the target according to the lag.
5. A dual polarized radar apparatus according to claim 1 or claim 3, wherein: the transmitting/receiving device is composed of a transmitter, a transmitting/receiving switch, an antenna and a receiver, wherein,
the transmitter is connected with the transmitting/receiving switch and is used for outputting a horizontal polarized wave or a vertical polarized wave to the transmitting/receiving switch;
the transmitting/receiving switch is connected with the transmitter, the antenna and the receiver and used for outputting the horizontal polarized wave or the vertical polarized wave input from the transmitter to the antenna and outputting the received horizontal polarized wave or the received vertical polarized wave to the receiver by the antenna;
and the antenna is connected with the transmitting/receiving switch and is used for transmitting the horizontal polarized wave or the vertical polarized wave output from the transmitting/receiving switch and simultaneously receiving the horizontal polarized wave or the vertical polarized wave reflected by the target and then returning the horizontal polarized wave or the vertical polarized wave to the receiver.
The receiver is connected with the transmitting/receiving switch and used for receiving the horizontal polarized wave or the vertical polarized wave received by the antenna, carrying out digital processing on the received signal, converting the analog received signal into a digital received signal and outputting the digital received signal to the signal processing device.
6. A dual polarized radar apparatus according to claim 1 or claim 3, wherein: the signal processing apparatus further includes a computer, and operation programs executed by the reflected wave intensity calculator, the linear depolarization ratio calculator, and the velocity calculator are stored in a memory of the computer, and at least one processor of the computer is configured to execute the programs stored in the memory.
7. The dual polarized radar apparatus of claim 1, 3 or 6, wherein: the signal processing device calculates the intensity of the reflected wave as the power value of the horizontal polarized wave and the vertical polarized wave by a reflected wave intensity calculator based on the digital received signals of the horizontal polarized wave and the vertical polarized wave output from the transmitting/receiving device, calculates a linear depolarization ratio LDR by a linear depolarization ratio calculator based on the intensity of the reflected wave, and calculates the Doppler velocity of the target by a velocity calculator
Figure FDA0003328888440000051
8. A dual polarized radar apparatus according to claim 1 or claim 3, wherein: the display device is mainly composed of a liquid crystal display and is used for displaying the linear depolarization ratio LDR and the Doppler velocity calculated by the signal processing device
Figure FDA0003328888440000052
9. A distributed target parameter estimation and signal processing method based on dual-polarized radar, which is implemented based on the dual-polarized radar device of any one of claims 1 to 8, and is characterized in that: which comprises the following steps:
ST1, when a transmitting/receiving switch of the transmitting/receiving equipment receives a horizontal polarized wave or a vertical polarized wave, the antenna transmits the horizontal polarized wave or the vertical polarized wave to the sky;
ST2, the receiver receives the horizontal polarized wave or the vertical polarized wave received by the antenna, converts an analog reception signal of the horizontal polarized wave or the vertical polarized wave into a digital reception signal, and outputs the digital reception signal to the signal processing device;
ST3, a reflected wave intensity calculator of the signal processing device calculates the intensity of the first reflected wave, the intensity of the second reflected wave, and the intensity of the third reflected wave based on the received signal;
ST4, calculating a linear depolarization ratio according to a radar reflectivity factor obtained by integrating the intensity of the reflected wave;
ST5, the Doppler velocity of the target is calculated from the intensity of the reflected wave.
CN202111274305.8A 2021-10-29 2021-10-29 Distributed target parameter estimation and signal processing method based on dual-polarization radar Pending CN113985414A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111274305.8A CN113985414A (en) 2021-10-29 2021-10-29 Distributed target parameter estimation and signal processing method based on dual-polarization radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111274305.8A CN113985414A (en) 2021-10-29 2021-10-29 Distributed target parameter estimation and signal processing method based on dual-polarization radar

Publications (1)

Publication Number Publication Date
CN113985414A true CN113985414A (en) 2022-01-28

Family

ID=79744601

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111274305.8A Pending CN113985414A (en) 2021-10-29 2021-10-29 Distributed target parameter estimation and signal processing method based on dual-polarization radar

Country Status (1)

Country Link
CN (1) CN113985414A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114637001A (en) * 2022-05-17 2022-06-17 艾索信息股份有限公司 Target detection method, device and equipment of radar and storage medium
CN117991198A (en) * 2024-04-07 2024-05-07 成都远望科技有限责任公司 Single-shot double-receiving top-sweeping cloud radar same-frequency interference identification method and device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114637001A (en) * 2022-05-17 2022-06-17 艾索信息股份有限公司 Target detection method, device and equipment of radar and storage medium
CN114637001B (en) * 2022-05-17 2022-08-30 艾索信息股份有限公司 Target detection method, device and equipment of radar and storage medium
CN117991198A (en) * 2024-04-07 2024-05-07 成都远望科技有限责任公司 Single-shot double-receiving top-sweeping cloud radar same-frequency interference identification method and device
CN117991198B (en) * 2024-04-07 2024-06-11 成都远望科技有限责任公司 Single-shot double-receiving top-sweeping cloud radar same-frequency interference identification method and device

Similar Documents

Publication Publication Date Title
KR102041767B1 (en) Dual polarized radar device and radar signal processing method
Zhang et al. Phased array radar polarimetry for weather sensing: A theoretical formulation for bias corrections
CN113985414A (en) Distributed target parameter estimation and signal processing method based on dual-polarization radar
US7898459B2 (en) Adaptive specific differential phase in dual-polarization radar
US9261583B2 (en) Systems and methods for calibrating dual polarization radar systems
CN107976660B (en) Missile-borne multi-channel radar ultra-low-altitude target analysis and multi-path echo modeling method
Vivekanandan et al. Radar reflectivity calibration using differential propagation phase measurement
WO1992022829A1 (en) Turbulence radar system
Mega et al. A low-power high-resolution broad-band radar using a pulse compression technique for meteorological application
Schvartzman et al. Signal processing and radar characteristics (SPARC) simulator: A flexible dual-polarization weather-radar signal simulation framework based on preexisting radar-variable data
WO2013181246A1 (en) Wide band clear air scatter doppler radar
Jameson et al. Estimation of propagation-differential phase shift from sequential orthogonal linear polarization radar measurements
US11614539B2 (en) Apparatus and method for removing noise for observation information of weather radar
Seliga et al. Differential reflectivity measurements in rain: First experiments
CN113281709A (en) Radar performance evaluation method based on area coupling forecasting system
CN103645475B (en) The spaceborne TOPSAR of complete polarization improves the method for cross polarization RASR
Falconi et al. Weather Radar Data Processing and Atmospheric Applications: An overview of tools for monitoring clouds and detecting wind shear
NAKAGAWA et al. P5. 8 DEVELOPMENT OF A NEW C-BAND BISTATIC POLARIMETRIC RADAR AND OBSERVATION OF TYPHOON EVENTS
RU2449312C1 (en) Panoramic radar method of determining parameters of state of ocean surface layer from satellite
CN108594235B (en) Method and system for improving reflectivity coefficient resolution of azimuth radar
Jacobi et al. Meteor radar wind and temperature measurements over Collm (51.3 N, 13 E) and comparison with co-located LF drift measurements during autumn 2004
Goldhirsh et al. Rain Measurement Results Derved from a Two-Polarzation Frequency-Diversity S-band Radar at Wallops Island, Virginia
Averyanova et al. Polarimetric method to discriminate the drops on sizes within the resolution volume
Honggang et al. The requirement of range sidelobe for dual-frequency precipitation radar on FengYun-3 Precipitation Measurement Satellite
Nakamura et al. Dual-polarization rain observation by X-and Ka-band dual-wavelength radar

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230103

Address after: No.17, Zhoushan Road, Jianxi District, Luoyang City, Henan Province

Applicant after: UNIT 63892 OF PLA

Address before: No.104 Tianjin Road, Jianxi District, Luoyang City, Henan Province

Applicant before: Dai Huanyao