CN109884337A - A method of wind direction of ocean surface is detected using high-frequency ground wave radar - Google Patents
A method of wind direction of ocean surface is detected using high-frequency ground wave radar Download PDFInfo
- Publication number
- CN109884337A CN109884337A CN201910162632.0A CN201910162632A CN109884337A CN 109884337 A CN109884337 A CN 109884337A CN 201910162632 A CN201910162632 A CN 201910162632A CN 109884337 A CN109884337 A CN 109884337A
- Authority
- CN
- China
- Prior art keywords
- wind direction
- extra large
- state unit
- bragg peak
- large state
- 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.)
- Granted
Links
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
The present invention provides a kind of methods using high-frequency ground wave radar detection wind direction of ocean surface, belong to high-frequency ground wave radar detection technology field, include the following steps, according to the beam direction angle theta of the wind direction of extra large state unit and radar emissionwAbsolute value | θw|, calculate the absolute value wait seek extra large state unit and with the wind direction of the multiple extra large state units adjacent wait seek extra large same distance element of state unit or so and place beam direction angle | θw|i, i=1,2 ..., m, wherein m is the total number wait seek extra large state unit and with the multiple wave beam sea state units adjacent wait seek extra large same distance element of state unit or so；And constitute a set C；Seek the linear function approachedAndCalculate the wind direction wait seek extra large state unit According to the positive-negative relationship of the k value, the judgement wind direction wait seek extra large state unit, obtain the wind direction value wait seek extra large state unit, compared with the prior art, influence the invention avoids conventional method due to external interference to a certain extra large state unit and lead to ambiguity solution Problem-Error, method is more effective, anti-interference is stronger.
Description
Technical field
The present invention relates to high-frequency ground wave radar detection technology fields, in particular to a kind of to detect sea using high-frequency ground wave radar
The method of face wind direction.
Background technique
High-frequency ground wave radar can break through the limitation of earth curvature, sight target below be detected, in radar return
Often adulterate a large amount of interference and noise.Sea return is one of most important interference in higher-frequency radar target acquisition, it
Energy is often very high, seriously affects the detection of target, copes with it and is inhibited；But on the other hand, sea return includes
A large amount of extra large state information, be that extra large state inverting predominantly detects object.
A large amount of observed result shows that the mechanism of production of sea return can be explained by Bragg resonance principle, that is, works as height
When frequency electromagnetic waves are irradiated on coarse sea, electromagnetic wave can interact with sea occurs strong scattering, to generate one
Rank sea return.Theoretically, single order sea return can be formed a pair of symmetrical about zero-frequency on the frequency axis that radar return RD is composed
Energy spikes, i.e. Bragg (Prague) peak, by calculate left and right Bragg peak-to-peak amplitude ratio, available Ocean Wind-field
Wind direction information.
The inverting of wind direction includes two parts, extracts wind direction and radar by the peak positive and negative single order Bragg of sea return spectrum first
The absolute value of squint；Then the ambiguity of wind direction is eliminated.Extract at present the method for wind direction and radar beam angle relative at
It is ripe, but the method effect for eliminating wind direction ambiguity is all undesirable.In the existing method for eliminating wind direction ambiguity, including multi-beam
Minimal difference method, the difference for drafting the wind direction measured on adjacent several extra large state units should be equal to zero after summation or be approximately equal to zero,
But this method is very likely by the wind direction while ambiguity solution mistake of adjacent several extra large state units；It further include using wait seek extra large state list
The size relation of the angle of member and the extra large state unit in the adjacent beams of left and right and wind direction judges wait ask the wind direction of extra large state unit can
It can be worth, but this method anti-interference is not strong, when wait ask extra large state unit and wind direction angle or adjacent beams sea state unit and wind
Any one into angle, which is measured, staggers the time, it may appear that ambiguity solution Problem-Error.
Summary of the invention
In view of this, the present invention is directed to propose it is a kind of using high-frequency ground wave radar detection wind direction of ocean surface method, including with
Lower step:
Step 1, according to the beam direction angle theta of the wind direction of extra large state unit and radar emissionwAbsolute value | θw|, calculate to
Ask extra large state unit and with the wind direction of adjacent wait seek extra large same distance element of state unit or so multiple extra large state units and place wave beam side
To the absolute value of angle | θw|i, i=1,2 ..., m, wherein m be wait ask extra large state unit and with wait seek the extra large same distance element of state unit
The total number of the adjacent multiple wave beam sea state units in left and right；
Step 2, by step 1 obtain described in wait ask extra large state unit and with it is adjacent wait seek extra large same distance element of state unit or so
Multiple extra large state units wind direction and place beam direction angle absolute value | θw|i, i=1,2 ..., m constitutes a set C；
Step 3, the set C obtained according to step 2, seeks the linear function approachedWherein For
Wave beam number,For the wind direction of corresponding wave beam Shanghai state unit and the absolute value of beam direction angle, k is slope and k is not that 0, b is
Intercept；
Step 4, the wind direction wait seek extra large state unit is calculatedAndWherein,To be described wait seek extra large state
The wind direction of unit,For the direction wait wave beam where seeking extra large state unit, | θw|nFor the wind direction wait seek extra large state unit with
The absolute value of place beam direction angle, ± represent the uncertainty of wind direction；；
Step 5, the positive-negative relationship of the k value obtained according to step 3, wait seek extra large state unit described in judgment step 4
Wind direction obtains the wind direction value wait seek extra large state unit.
Preferably, judging the wind direction wait seek extra large state unit in step 5 and obtaining the wind direction wait seek extra large state unit
The method of value is as follows:
If k > 0, and when the angle of wave beam is sequentially increased, which isWork as wave beam
Angle when being sequentially reduced, which is
If k < 0, when the angle of wave beam is sequentially increased, which isWhen wave beam
When angle is sequentially reduced, judge that the sea state unit wind direction value is
Wherein the angle of the wave beam increases or reduces, with described wait ask the beam direction where extra large state unit for ginseng
It examines, is to increase in the angle of its clockwise wave beam, be to reduce in the angle of its anticlockwise wave beam.
Preferably, the absolute value of the beam direction angle of the wind direction and radar emission of the unit of sea state described in step 1 | θw|
Step obtains by the following method:
Step 11, according to the working frequency of radar, and theoretical Bragg frequency f is calculated by formula 2B, the formula 2 is
Wherein f0For the working frequency of radar；
Step 12, according to the wavelength of radar emission electromagnetic wave, and pass through the maximum value Δ f of the calculating frequency deviation of formula 3m, described
Formula 3 is
Wherein λ is the wavelength of radar emission electromagnetic wave, VmFor the maximum value of ocean current radial flow speed；
Step 13, according to the maximum value Δ f of Prague frequency deviation acquired in step 12mWith theory cloth obtained in step 11
Glug frequency fB, determine the search range of single order bragg peak；
Step 14, positive single order bragg peak and negative one are determined in the search range of the single order bragg peak described in step 13
Rank bragg peak；
Step 15, positive single order bragg peak and negative one rank bragg peak amplitude ratio R, the formula 4 are calculated according to formula 4
For
Wherein B+Be positive the range value of single order bragg peak, B-Be negative the range value of single order bragg peak.
Step 16, the positive single order bragg peak and negative one rank bragg peak Amplitude Ratio obtained according to formula 5 and step 15
Value R, obtains | θw|=2arctan [10R/(10s)], the formula 5 is R=10lg [tans(|θw|/2)], wherein s be extension because
Son, represents degree of scatter of the wave about wind direction, and the value range of s is [2,16].
Preferably, the maximum value V of ocean current radial flow speed described in step 12mIt is passed through by the extra large state data in radar detection sea area
Test acquisition.
Preferably, the search range of single order bragg peak described in step 13, the search model including positive single order bragg peak
The search range with negative one rank bragg peak is enclosed, and the search range of the positive single order bragg peak is [fB-Δfm,fB+Δ
fm], the search range of the negative one rank bragg peak is [- fB-Δfm,-fB+Δfm]。
Preferably, positive single order bragg peak and negative one rank bragg peak described in step 14 determine in the following manner,
Amplitude maxima is taken respectively in the search range of the positive single order bragg peak and the search range of the negative one rank bragg peak
Point is used as the positive single order bragg peak and the negative one rank bragg peak.
Preferably, calculating the linear function approached described in step 4 according to formula 7In k value and b value, institute
Stating formula 7 isWherein ajFor the coefficient of approximating function to be asked, and a1=k, a0=b.
Preferably, X in the formula 7jlAnd Yl, it is obtained by formula 8, the formula 8 is as follows:
WhereinWithFor selected linear independence function, xiFor the wave beam number where i-th of extra large state unit, yiFor xthiExtra large state list on wave beam
The wind direction of member and the absolute value of beam direction angle | θw|i, XjlFor the linear independence function on m pointWith's
The sum of products, YlFor on m pointAnd yiThe sum of products.
Compared with conventional method, the present invention solves wind direction value using the data information on the adjacent multiple extra large state units in left and right,
A kind of method using high-frequency ground wave radar detection wind direction of ocean surface is provided, avoids conventional method since external interference is to a certain
The influence of extra large state unit and lead to ambiguity solution Problem-Error, method is more effective, anti-interference is stronger.
Detailed description of the invention
The attached drawing for constituting a part of the invention is used to provide further understanding of the present invention, schematic reality of the invention
It applies example and its explanation is used to explain the present invention, do not constitute improper limitations of the present invention.In the accompanying drawings:
Fig. 1 is a kind of method flow diagram that wind direction of ocean surface is detected using high-frequency ground wave radar in the present invention；
Fig. 2 is the doppler echo spectrogram of the 47th range gate of the 4th wave beam in the embodiment of the present invention 1；
Fig. 3 is that wind direction of the 4th wave beam the 47th apart from gate cell solves schematic diagram in the embodiment of the present invention 1；
Fig. 4 is wind direction inversion result schematic diagram in the embodiment of the present invention 1；
Fig. 5 is the doppler echo spectrogram of the 30th range gate of the 4th wave beam in the embodiment of the present invention 2；
Fig. 6 is that wind direction of the 4th wave beam the 30th apart from gate cell solves schematic diagram in the embodiment of the present invention 2；
Fig. 7 is wind direction inversion result schematic diagram in the embodiment of the present invention 2.
Specific embodiment
The present invention will be described in detail below with reference to the accompanying drawings and embodiments.
The present invention is based in a fritter sea area wind direction be this identical or gradual rule, provide and a kind of utilize high frequency
The method that ground wave radar detects wind direction of ocean surface, will be wait seek the adjacent multiple wave beam sea state units of extra large state unit and same distance element or so
The absolute value of upwind and beam direction angle constitutes a set, approximating function is sought, according to the changing rule of approximating function
Wind direction is fuzzy to obtain true wind direction eliminating, specifically includes the following steps:
Step 1, according to the beam direction angle theta of the wind direction of extra large state unit and radar emissionwAbsolute value | θw|, calculate to
Ask extra large state unit and with the wind direction of adjacent wait seek extra large same distance element of state unit or so multiple extra large state units and place wave beam side
To the absolute value of angle | θw|i, i=1,2 ..., m, wherein m be wait ask extra large state unit and with wait seek the extra large same distance element of state unit
The total number of the adjacent multiple wave beam sea state units in left and right；
Step 2, by step 1 obtain described in wait ask extra large state unit and with it is adjacent wait seek extra large same distance element of state unit or so
Multiple extra large state units wind direction and place beam direction angle absolute value | θw|i, i=1,2 ..., m constitutes a set C；
Step 3, the set C obtained according to step 2, seeks the linear function approachedWherein For
Wave beam number,For the wind direction of corresponding wave beam Shanghai state unit and the absolute value of beam direction angle, k is slope and k is not that 0, b is
Intercept；
Step 4, the wind direction wait seek extra large state unit is calculatedAndWherein,To be described wait seek extra large state
The wind direction of unit,For the direction wait wave beam where seeking extra large state unit, | θw|nFor the wind direction wait seek extra large state unit with
The absolute value of place beam direction angle, ± represent the uncertainty of wind direction, i.e., when wind direction of ocean surface isWithWhen, radar observation to single order sea return be identical；
Step 5, the positive-negative relationship of the k value obtained according to step 3, wait seek extra large state unit described in judgment step 4
Wind direction obtains the wind direction value wait seek extra large state unit.
Wherein, the absolute value of the beam direction angle of the wind direction and radar emission of the unit of sea state described in step 1 | θw| it is logical
Cross the acquisition of following methods step:
Step 11, according to the working frequency of radar, and theoretical Bragg frequency f is calculated by formula 2B, the formula 2 is
Wherein f0For the working frequency of radar；
Step 12, according to the wavelength of radar emission electromagnetic wave, and pass through the maximum value Δ f of the calculating frequency deviation of formula 3m, described
Formula 3 is
Wherein λ is the wavelength of radar emission electromagnetic wave, VmFor the maximum value of ocean current radial flow speed, and the ocean current radial flow
The maximum value V of speedmIt is obtained by the extra large state data experience in radar detection sea area；And frequency deviation represents practical bragg peak relative to reason
Frequency deviation will not occur, when the radial velocity of ocean current is not 0 when the radial velocity of ocean current is 0 by the offset amplitude of bragg peak
When, then frequency deviation can occur.
Step 13, according to the maximum value Δ f of Prague frequency deviation acquired in step 12mWith theory cloth obtained in step 11
Glug frequency fB, determine the search range of single order bragg peak, and the search range of the single order bragg peak includes positive single order
The search range of bragg peak and the search range of negative one rank bragg peak, the search range of the positive single order bragg peak are
[fB-Δfm,fB+Δfm], the search range of the negative one rank bragg peak is [- fB-Δfm,-fB+Δfm]；
Step 14, positive single order bragg peak and negative one are determined in the search range of the single order bragg peak described in step 13
Rank bragg peak divides that is, in the search range of the positive single order bragg peak and the search range of the negative one rank bragg peak
Do not take Amplitude maxima point as the positive single order bragg peak and the negative one rank bragg peak；The wherein positive single order Bradley
Lattice peak correspond to towards radar direction wave, the negative one rank bragg peak correspond to far from radar direction wave, two
The size of the intensity of bragg peak reflects wind direction, when the intensity of positive single order bragg peak is greater than the intensity of negative one rank bragg peak
When, wind direction is shown as backwards to beam direction, then shows as wind direction on the contrary along beam direction, when the intensity at two peaks is identical, wind direction
Perpendicular to wave beam.
Step 15, positive single order bragg peak and negative one rank bragg peak amplitude ratio R are calculated according to formula 4, wherein R is dB
Format, the formula 4 are
Wherein B+Be positive the range value of single order bragg peak, B-Be negative the range value of single order bragg peak.
Step 16, the positive single order bragg peak and negative one rank bragg peak Amplitude Ratio obtained according to formula 5 and step 15
Value R, obtains | θw|=2arctan [10R/(10s)], the formula 5 is R=10lg [tans(|θw|/2)], wherein s be extension because
Son represents degree of scatter of the wave about wind direction, and the value range of s is [2,16], and preferred value is 4 in the present embodiment.
The linear function wherein approached described in step 3In k value and b value obtained by formula 7, the public affairs
Formula 7 isWherein ajFor the coefficient of approximating function to be asked, and a1=k, a0=b, in the formula 7
XjlAnd Yl, it is obtained by formula 8, the formula 8 is as follows:
WhereinWithFor selected linear independence function, xiFor the wave beam number where i-th of extra large state unit, yiFor xthiExtra large state list on wave beam
The wind direction of member and the absolute value of beam direction angle | θw|i, XjlFor the linear independence function on m pointWith's
The sum of products, YlFor on m pointAnd yiThe sum of products.
To sum up, if k > 0, angle absolute value between true wind direction and each wave beam | θw|i, i=1,2 ..., m present by
Cumulative big trend, when the angle of wave beam is sequentially increased, which isI.e.When the angle of wave beam is sequentially reduced, which is
I.e.
If k < 0, angle absolute value between true wind direction and each wave beam | θw|i, i=1,2 ..., m presentation gradually subtracts
Small trend, when the angle of wave beam is sequentially increased, which isI.e.When the angle of wave beam is sequentially reduced, which isI.e.
Wherein the angle of the wave beam increases or reduces, with described wait ask the beam direction where extra large state unit for ginseng
It examines, is to increase in the angle of its clockwise wave beam, be to reduce in the angle of its anticlockwise wave beam.
Therefore, the present invention solves wind direction value using the data information on the adjacent multiple extra large state units in left and right, provides one kind
Using the method for high-frequency ground wave radar detection wind direction of ocean surface, conventional method is avoided since external interference is to a certain extra large state unit
It influences and leads to ambiguity solution Problem-Error, method is more effective, anti-interference is stronger.
Embodiment 1
Certain batch data that the data source Mr. Yu's radar experiment centre used in the present embodiment was acquired on November 2nd, 2018,
Wherein the tranmitting frequency of radar is 8.4MHz, and wave number is 7.
In the doppler echo spectrogram of the 47th range gate of the 4th wave beam of Fig. 2, * point is labelled with the position at the peak Bragg, while can
Obtain the range value at the positive and negative peak Bragg.
Shadow region unit B represents the 4th wave beam the 47th apart from gate cell in Fig. 3, i.e., wait seek extra large state unit, shadow region is single
First A represents the 3rd wave beam the 47th apart from gate cell, and shadow region unit C represents the 5th wave beam the 47th apart from gate cell, wave beam a, b and
The angular dimension relationship of c isAnd acquire the positive and negative single order p-ratio R of A, B and C tri- extra large state unitsA、RBAnd RC
Respectively -6.4, -8.9, -10.6, the angle absolute value of the wind direction and wave beam a, b and c of A, B and C tri- extra large state units | θA|、|θB
|、|θC| it is respectively 69.4 degree, 61.6 degree, 57.1 degree, thus obtained approximating function isK < 0 at this time,
Illustrate that the angle absolute value between the extra large state unit wind directions of A, B and C tri- and wave beam a, b and c is presented the trend being gradually reduced, and wave
Beam angle degreeThe wind direction value for obtaining the sea state unit is
Usual wind direction refers to the direction of wind, as seen in Figure 4, the resulting wind direction of different sea state unit invertings it is identical or
Close, substantially southwester is also consistent with the weather forecast on the same day.
Embodiment 2
Certain batch data that the data source Mr. Yu's radar experiment centre used in the present embodiment was acquired on November 5th, 2018,
Wherein the tranmitting frequency of radar is 5.1MHz, and wave number is 7.
In the doppler echo spectrum of the 30th range gate of the 4th wave beam of Fig. 5, * point is labelled with the position at the peak Bragg, while can obtain
To the range value at the positive and negative peak Bragg.
The shadow region Fig. 6 unit B represents the 4th wave beam the 30th apart from gate cell, i.e., wait ask extra large state unit, shadow region unit
A represents the 3rd wave beam the 30th apart from gate cell, and shadow region unit C represents the 5th wave beam the 30th apart from gate cell, wave beam a, b and c
Angular dimension relationship beAnd acquire the positive and negative single order p-ratio R of A, B and C tri- extra large state unitsA、RBAnd RC
Respectively 14.5,16.6,20.7, A, the angle absolute value of the wind direction and wave beam a, b and c of B and C tri- extra large state units | θA|、|θB
|、|θC| it is respectively 133.1 degree, 137.8 degree, 146.2 degree, thus obtained approximating function isK > at this time
0, illustrate that the angle absolute value of the extra large state unit wind directions of A, B and C tri- and wave beam a, b and c are presented the trend being gradually increased, and wave beam
AngleThe wind direction value for obtaining the sea state unit is
Usual wind direction refers to the direction of wind, as seen in Figure 7, the resulting wind direction of different sea state unit invertings it is identical or
Close, substantially north wind is also consistent with the weather forecast on the same day.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Within mind and principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (8)
1. a kind of method using high-frequency ground wave radar detection wind direction of ocean surface, which comprises the following steps:
Step 1, according to the beam direction angle theta of the wind direction of extra large state unit and radar emissionwAbsolute value | θw|, it calculates wait ask sea
It state unit and is pressed from both sides with the wind direction of adjacent wait seek extra large same distance element of state unit or so multiple extra large state units and place beam direction
The absolute value at angle | θw|i, i=1,2 ..., m, wherein m be wait ask extra large state unit and with wait seek extra large same distance element of state unit or so
The total number of adjacent multiple wave beam sea state units；
Step 2, by step 1 obtain described in wait seek extra large state unit and with described wait ask extra large same distance element of state unit or so adjacent
Multiple extra large state units wind direction and place beam direction angle absolute value | θw|i, i=1,2 ..., m constitutes a set C；
Step 3, the set C obtained according to step 2, seeks the linear function approachedWherein For wave beam
Number,For the wind direction of corresponding wave beam Shanghai state unit and the absolute value of beam direction angle, k is slope and k be not 0, b is to cut
Away from；
Step 4, the wind direction wait seek extra large state unit is calculated Wherein,To be described wait seek extra large state unit
Wind direction,For the direction wait wave beam where seeking extra large state unit, | θw|nFor the wind direction wait seek extra large state unit and place wave
The absolute value of beam angular separation, ± represent the uncertainty of wind direction；
Step 5, the positive-negative relationship of the k value obtained according to step 3, wait seek the wind direction of extra large state unit described in judgment step 4,
Obtain the wind direction value wait seek extra large state unit.
2. the method according to claim 1 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that step 5
Wait seek the wind direction of extra large state unit and obtain described wait ask the method for the wind direction value of extra large state unit as follows described in middle judgement:
If k > 0, and when the angle of wave beam is sequentially increased, which isWhen the angle of wave beam
When degree is sequentially reduced, which is
If k < 0, when the angle of wave beam is sequentially increased, which isWhen the angle of wave beam
When being sequentially reduced, judge that the sea state unit wind direction value is
Wherein the angle of the wave beam increases or reduces, with described wait ask the beam direction where extra large state unit for reference,
The angle of its clockwise wave beam is to increase, and is to reduce in the angle of its anticlockwise wave beam.
3. the method according to claim 1 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that step 1
Described in sea state unit wind direction and radar emission beam direction angle absolute value | θw| step obtains by the following method:
Step 11, according to the working frequency of radar, and theoretical Bragg frequency f is calculated by formula 2B, the formula 2 is
Wherein f0For the working frequency of radar；
Step 12, according to the wavelength of radar emission electromagnetic wave, and pass through the maximum value Δ f of the calculating frequency deviation of formula 3m, the formula 3
For
Wherein λ is the wavelength of radar emission electromagnetic wave, VmFor the maximum value of ocean current radial flow speed；
Step 13, according to the maximum value Δ f of Prague frequency deviation acquired in step 12mWith theory Prague obtained in step 11
Frequency fB, determine the search range of single order bragg peak；
Step 14, positive single order bragg peak and negative one rank cloth are determined in the search range of the single order bragg peak described in step 13
Glug peak；
Step 15, positive single order bragg peak and negative one rank bragg peak amplitude ratio R are calculated according to formula 4, the formula 4 is
Wherein B+Be positive the range value of single order bragg peak, B-Be negative the range value of single order bragg peak.
Step 16, the positive single order bragg peak and negative one rank bragg peak amplitude ratio R obtained according to formula 5 and step 15,
Obtain | θw|=2arctan [10R/(10s)], the formula 5 is R=10lg [tans(|θw|/2)], wherein s is spreading factor, generation
Degree of scatter of the table wave about wind direction, the value range of s are [2,16].
4. the method according to claim 3 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that step 12
Described in ocean current radial flow speed maximum value VmIt is obtained by the extra large state data experience in radar detection sea area.
5. the method according to claim 4 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that step 13
Described in single order bragg peak search range, search range and negative one rank bragg peak including positive single order bragg peak are searched
Rope range, and the search range of the positive single order bragg peak is [fB-Δfm,fB+Δfm], the negative one rank bragg peak
Search range is [- fB-Δfm,-fB+Δfm]。
6. the method according to claim 5 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that step 14
Described in positive single order bragg peak and negative one rank bragg peak determine in the following manner, in searching for the positive single order bragg peak
Take Amplitude maxima point as described positive single order Prague in the search range of rope range and the negative one rank bragg peak respectively
Peak and the negative one rank bragg peak.
7. the method according to claim 1 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that according to public affairs
Formula 7 calculates the linear function approached described in step 3In k value and b value, the formula 7 isWherein ajFor the coefficient of approximating function to be asked, and a1=k, a0=b.
8. the method according to claim 7 using high-frequency ground wave radar detection wind direction of ocean surface, which is characterized in that the public affairs
X in formula 7jlAnd Yl, it is obtained by formula 8, the formula 8 is as follows:
WhereinWithFor selected linear independence function, xiFor the wave beam number where i-th of extra large state unit, yiFor xthiExtra large state list on wave beam
The wind direction of member and the absolute value of beam direction angle | θw|i,XjlFor the linear independence function on m pointWith's
The sum of products, YlFor on m pointAnd yiThe sum of products.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910162632.0A CN109884337B (en) | 2019-03-05 | 2019-03-05 | Method for detecting sea surface wind direction by using high-frequency ground wave radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910162632.0A CN109884337B (en) | 2019-03-05 | 2019-03-05 | Method for detecting sea surface wind direction by using high-frequency ground wave radar |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109884337A true CN109884337A (en) | 2019-06-14 |
| CN109884337B CN109884337B (en) | 2021-01-19 |
Family
ID=66930662
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910162632.0A Active CN109884337B (en) | 2019-03-05 | 2019-03-05 | Method for detecting sea surface wind direction by using high-frequency ground wave radar |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109884337B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113391287A (en) * | 2021-06-10 | 2021-09-14 | 哈尔滨工业大学 | High-frequency ground wave radar sea state data fusion method based on time sequence |
| CN113466854A (en) * | 2021-06-29 | 2021-10-01 | 哈尔滨工业大学 | High-frequency ground wave radar inversion vector flow velocity method based on ocean power model |
| CN115267717A (en) * | 2022-07-29 | 2022-11-01 | 哈尔滨工业大学 | Wind direction included angle inversion method, ambiguity resolution method, storage medium and equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2012133369A (en) * | 2012-08-03 | 2014-02-10 | Федеральное государственное бюджетное учреждение науки Институт прикладной астрономии Российской академии наук | DEVICE FOR REMOTE SPEED AND WIND DIRECTION |
| CN103616690A (en) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | Sea surface wind direction extraction method based on ship-borne high-frequency ground wave beyond visual range radar |
| CN108008392A (en) * | 2017-11-22 | 2018-05-08 | 哈尔滨工业大学 | A kind of ocean surface Wind field measurement method based on boat-carrying high-frequency ground wave radar |
| CN109143193A (en) * | 2018-08-30 | 2019-01-04 | 武汉大学 | A kind of high-frequency ground wave radar wave field inversion method based on second order peak |
-
2019
- 2019-03-05 CN CN201910162632.0A patent/CN109884337B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2012133369A (en) * | 2012-08-03 | 2014-02-10 | Федеральное государственное бюджетное учреждение науки Институт прикладной астрономии Российской академии наук | DEVICE FOR REMOTE SPEED AND WIND DIRECTION |
| CN103616690A (en) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | Sea surface wind direction extraction method based on ship-borne high-frequency ground wave beyond visual range radar |
| CN108008392A (en) * | 2017-11-22 | 2018-05-08 | 哈尔滨工业大学 | A kind of ocean surface Wind field measurement method based on boat-carrying high-frequency ground wave radar |
| CN109143193A (en) * | 2018-08-30 | 2019-01-04 | 武汉大学 | A kind of high-frequency ground wave radar wave field inversion method based on second order peak |
Non-Patent Citations (2)
| Title |
|---|
| WEIMIN HUANG等: "Measuring Surface Wind Direction by Monostatic HF Ground-Wave Radar at the Eastern China Sea", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 * |
| 王曙曜等: "OS081H 高频地波雷达系统海面风向反演实验研究", 《电子与信息学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113391287A (en) * | 2021-06-10 | 2021-09-14 | 哈尔滨工业大学 | High-frequency ground wave radar sea state data fusion method based on time sequence |
| CN113391287B (en) * | 2021-06-10 | 2023-09-01 | 哈尔滨工业大学 | High-frequency ground wave radar sea state data fusion method based on time sequence |
| CN113466854A (en) * | 2021-06-29 | 2021-10-01 | 哈尔滨工业大学 | High-frequency ground wave radar inversion vector flow velocity method based on ocean power model |
| CN115267717A (en) * | 2022-07-29 | 2022-11-01 | 哈尔滨工业大学 | Wind direction included angle inversion method, ambiguity resolution method, storage medium and equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109884337B (en) | 2021-01-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104777453B (en) | The Beam Domain Time-Frequency Analysis Method of naval vessel line spectrum noise source positioning | |
| CN109884337A (en) | A method of wind direction of ocean surface is detected using high-frequency ground wave radar | |
| Bai et al. | High-resolution three-dimensional imaging of spinning space debris | |
| Mayhan et al. | High resolution 3D" snapshot" ISAR imaging and feature extraction | |
| US20110223031A1 (en) | wind turbine, a blade therefor and a method of processing signals reflected therefrom | |
| CN102565784B (en) | Method of moving-target relocation and velocity ambiguity resolution based on velocity synthetic aperture radar (VSAR) system | |
| CN107656255B (en) | Ultra-wideband radar moving target two-dimensional positioning method based on multipath echoes | |
| CN109324322A (en) | A kind of direction finding and target identification method based on passive phased array antenna | |
| CN102736073B (en) | Method for computing range ambiguity of satellite-borne synthetic aperture radar (SAR) in universal mode | |
| CN113296073A (en) | Amplitude-phase joint blade clearance measurement method based on envelope constraint | |
| Gill et al. | On the development of a second-order bistatic radar cross section of the ocean surface: A high-frequency result for a finite scattering patch | |
| CN112904326A (en) | Satellite-borne passive positioning method based on virtual aperture | |
| CN109655827A (en) | A kind of concealed target localization method for single building corner | |
| Poulsen et al. | The scattering of VLF waves by localized ionospheric disturbances produced by lightning‐induced electron precipitation | |
| CN106959442B (en) | Method for extracting first-order echo spectrum of ground wave radar in strong interference environment based on multi-domain information | |
| CN106646466A (en) | Imaging method of weighted back projection algorithm based on principal component analysis | |
| Sun et al. | Parameter estimation method of walking human based on radar micro-Doppler | |
| CN105467373A (en) | Method for estimating physical dimension of conical target of broadband composite bistatic radar | |
| CN112198490B (en) | Ground clutter suppression method and device for airborne millimeter wave cloud detection radar and radar | |
| CN111273283A (en) | Trunk movement parameter extraction method based on human body three-dimensional micro Doppler signal | |
| Walsh et al. | The second-order high frequency radar ocean surface foot-scatter cross section for an antenna on a floating platform | |
| Reiten et al. | Properties of surface waves determined via bistatic terahertz impulse ranging | |
| Laurukevich et al. | Estimation of energy, spectral and polarimetric characteristics of meteorological echoes in DMRL-C | |
| Shen et al. | Simulation of HF radar cross sections for swell contaminated seas | |
| CN110109186A (en) | A kind of Coherent Noise in GPR Record three-dimensional Time-Frequency Analysis Method |
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 |