METHOD FOR ESTIMATING NONDIRECTIONAL WAVE SPECTRUM FROM SEA ECHOES OF MULTIPLE HIGH RADAR FREQUENCIES

Abstract

The disclosure provides a method for estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies. The method includes: dividing the radar detection area into a plurality of fan-shaped units at an equal range interval and angle interval according to the distance resolution and the angular resolution of an HF radar; obtaining the Doppler spectrum from the sea echo of a single radar frequency at a fan-shaped unit by performing the first fast Fourier transform (FFT) in distance dimension, the second FFT in Doppler frequency dimension and the digital beamforming; extracting the positive first-order peak and the negative first-order peak from the aforementioned Doppler spectrum by the peak-searching method; and selecting the stronger first-order peak σ.sub.R.sup.(1)(ω) ; dividing the second-order spectrum on the stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum.

Claims

1. A method for directly estimating a nondirectional wave spectrum from sea echoes of multiple HF radar frequencies, the method comprising: (a) dividing a radar detection area into a plurality of fan-shaped units at an equal range interval and angle interval according to a distance resolution and an angular resolution of a multi-frequency HF radar, wherein the multi-frequency HF radar is capable of simultaneously operating at more than one frequency in the HF band; (b) obtaining a Doppler spectrum from a sea echo of a single radar frequency at a fan-shaped unit by performing a first fast Fourier transform (FFT) in distance dimension, a second FFT in Doppler frequency dimension and a digital beamforming in a signal processor unit, extracting a positive first-order peak and a negative first-order peak from the Doppler spectrum by a peak-searching method and selecting a stronger first-order peak σ.sub.R.sup.(1)(ω) from the positive and negative first-order peaks in a computer processor; (c) dividing a second-order spectrum on a stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum and separating the outer second-order spectrum on the stronger first-order peak side from the Doppler spectrum in the computer processor according to a Doppler frequency range of the outer second-order spectrum; (d) calculating R.sub.f(ω) which is a function of the outer second-order spectrum divided by a first-order peak energy for the stronger side of the Doppler spectrum of the single radar frequency at the fan-shaped unit in the computer processor; (e) calculating A.sub.f which is defined as a single-frequency coefficient matrix of the nondirectional wave spectrum by linearizing the outer second-order spectrum on the side of the stronger first-order peak in the computer processor; (f) repeating the above steps (b) to (e), combining the sea echoes of different radar frequencies to construct a matrix R which is a function of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum of multiple radar frequencies and merging coefficient matrices of the nondirectional wave spectrum from multiple radar frequencies into a matrix A in the computer processor; and (g) calculating a pseudo-inverse A.sup.+ of the matrix A by a singular value decomposition and estimating the nondirectional wave spectrum at the fan-shaped unit from the matrix R and the pseudo-inverse A.sup.+ in the computer processor.

2. The method of claim 1, wherein in (b), the Doppler spectrum at the fan-shaped unit is defined as: σ(ω), where ω represents a Doppler frequency generated by the motion of ocean waves to the multi-frequency HF radar; and σ(ω) represents a wave energy distribution at a different value of ω; wherein in (b), the positive first-order peak and negative first-order peak are extracted from the Doppler spectrum at the fan-shaped unit using the peak searching method; the first-order peaks are defined as two peaks in the Doppler spectrum which are roughly symmetrically distributed on both sides of zero frequency; the first-order peaks are generated by the Bragg scattering of the waves of half the radar wavelength that are either advancing directly towards the multi-frequency HF radar or receding directly from the multi-frequency HF radar; wherein (b) comprises: searching for a point of a maximum amplitude in the Doppler frequency range [0.6ω.sub.B,1.4ω.sub.B] of the Doppler spectrum at the fan-shaped unit as the peak of the positive first-order peak, and recording the Doppler frequency of the peak of the positive first-order peak as ω.sub.P+, wherein ω.sub.B=√{square root over (2gk.sub.0)}, is a Bragg frequency, wherein $k_0=\frac{2\pi f}{c}$ is a radar wavenumber; c is a speed of light; and f is a radar frequency; searching for a local minimum point inside the peak of the positive first-order peak where the Doppler frequency meets the inequation ω.sub.P+−0.2ω.sub.B≤ω<ω.sub.P+, and denoting the Doppler frequency of the local minimum point inside the peak of the positive first-order peak as ω.sub.L+; searching for the local minimum point outside the peak of the positive first-order peak where the Doppler frequency satisfies the inequation ω.sub.P+<ω≤0.2ω.sub.B+ω.sub.P+, and recording the Doppler frequency of the local minimum point outside the peak of the positive first-order peak as ω.sub.R+; intercepting the Doppler spectrum at the fan-shaped unit with Doppler frequency [ω.sub.L+, ω.sub.R+] as the positive first-order peak; searching for the point of the maximum amplitude in the Doppler frequency range [−1.4ω.sub.B, −0.6ω.sub.B] of the Doppler spectrum at the fan-shaped unit as the peak of the negative first-order peak, and recording the Doppler frequency of the peak of the negative first-order peak as ω.sub.P−; searching for the local minimum point inside the peak of the negative first-order peak where the Doppler frequency meets the inequation ω.sub.P−<ω≤ω.sub.P−+0.2ω.sub.B, and denoting the Doppler frequency of the local minimum point inside the peak of the negative first-order peak as ω.sub.L−; searching for the local minimum point outside the peak of the negative first-order peak where the Doppler frequency satisfies the inequation ω.sub.P−−0.2φ.sub.B≤ω<ω.sub.P−, and recording the Doppler frequency of the local minimum point outside the peak of the negative first-order peak as ω.sub.R−; intercepting the Doppler spectrum at the fan-shaped unit with Doppler frequency [ω.sub.R−, ω.sub.L−] as the negative first-order peak; and comparing the amplitude of the peak of the positive first-order peak with that of the negative first-order peak, and selecting the first-order peak with a larger amplitude of the peak point as the stronger first-order peak σ.sub.R.sup.(1)(ω).

3. The method of claim 1, wherein in (c), the second-order spectrum originated from the second-order scattering of ocean waves and radar waves is a continuum with lower amplitude than the first-order peaks and distributed around the first-order peaks; wherein (c) comprises the steps of: (1) dividing the second-order spectrum on the stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum; the Doppler frequency range of the outer second-order spectrum on the stronger first-order peak side is given as $\{\begin{array}{c}{\omega }_{c+}<\omega \le 1.4{\omega }_B,\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ -1.4{\omega }_B\le \omega <{\omega }_{c-},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array};$ the Doppler frequency range of the inner second-order spectrum on the stronger first-order peak side is given as $\{\begin{array}{c}0.6{\omega }_B\le \omega <{\omega }_{c+},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ {\omega }_{c-}<\omega \le -0.6{\omega }_B,\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array};$ where ω is a Doppler frequency generated by the motion of ocean waves to a radar; ω.sub.B=√{square root over (2gk.sub.0)} is a Bragg frequency; k.sub.0 is a radar wavenumber; ${\omega }_{c+}=\frac{\stackrel{{\omega }_{R+}}{\underset{{\omega }_{L+}}{.Math.}}\omega .Math.{\sigma }_R^{\left(1\right)}\left(\omega \right)}{\stackrel{{\omega }_{R+}}{\underset{{\omega }_{L+}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)}\mathrm{is}a\mathrm{centroid}\mathrm{frequency}\mathrm{of}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak};\mathrm{and}$ ${\omega }_{c-}=\frac{\stackrel{{\omega }_{L-}}{\underset{{\omega }_{R-}}{.Math.}}\omega .Math.{\sigma }_R^{\left(1\right)}\left(\omega \right)}{\stackrel{{\omega }_{L-}}{\underset{{\omega }_{R-}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)}\mathrm{is}a\mathrm{centroid}\mathrm{frequency}\mathrm{of}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak};$ (2) if the stronger first-order peak obtained in (b) of claim 1 is the positive first-order peak, taking the centroid frequency ω.sub.c+, as a starting point of Doppler frequency, taking a cutoff frequency of 1.4ω.sub.B as an endpoint, extracting the Doppler spectrum at the fan-shaped unit σ(ω) in the frequency range of (ω.sub.c+,1.4ω.sub.B], and setting the Doppler spectrum at the fan-shaped unit σ(ω) in the frequency range of (ω.sub.c+,ω.sub.R+] to zero to obtain the outer second-order spectrum on the stronger first-order peak side σ.sub.R.sup.(2)(ω); if the stronger first-order peak obtained in (b) of claim 1 is the negative first-order peak, taking the centroid frequency ω.sub.c− as the starting point of Doppler frequency, taking a cutoff frequency of −1.4ω.sub.B as the endpoint, extracting the Doppler spectrum at the fan-shaped unit σ(ω) in the frequency range of [−1.4ω.sub.B,ω.sub.c−), and setting the Doppler spectrum at the fan-shaped unit σ(ω) in the frequency range of [ω.sub.R−,ω.sub.c−) to zero to obtain the outer second-order spectrum on the stronger first-order peak side σ.sub.R.sup.(2)(ω).

4. The method of claim 1, wherein in (d), the aforementioned function R.sub.f(ω) is given by the equation: $\{\begin{array}{c}{R}_f\left(\omega \right)=\frac{{\sigma }_R^{\left(2\right)}\left(\omega \right)}{\stackrel{{\omega }_{R+}}{\underset{{\omega }_{L+}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{\Delta \omega }},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ R_f\left(\omega \right)=\frac{{\sigma }_R^{\left(2\right)}\left(\omega \right)}{\stackrel{{\omega }_{L-}}{\underset{{\omega }_{R-}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{\Delta \omega }},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array};$ where f in a subscript indicates a radar frequency; σ.sub.R.sup.(1)(ω) is the stronger first-order peak; σ.sub.R.sup.(2) (ω) is the outer second-order spectrum on the stronger first-order peak side; ω.sub.L+ is the Doppler frequency of the local minimum point inside the peak of the positive first-order peak; ω.sub.R+ is the Doppler frequency of the local minimum point outside the peak of the positive first-order peak; ω.sub.L− is the Doppler frequency of the local minimum point inside the peak of the negative first-order peak; ω.sub.R− is the Doppler frequency of the local minimum point outside the negative first-order peak; and Δω is a frequency resolution of the Doppler spectrum σ(ω).

5. The method of claim 1, wherein (e) comprises the steps of: (1) linearizing the outer second-order spectrum on the side of the stronger first-order peak; at the Doppler frequencies near the first-order peak that satisfy the condition ω.sub.B<|ω|≤1.4ω.sub.B, the direction of one of the two sets of ocean waves that generate the second-order scattering with the radar vector is approximately equal to that of the Bragg wave vector, so the outer second-order spectrum σ.sub.R.sup.(2) (ω) is linearized as σ.sub.RL.sup.(2)(ω): ${\sigma }_{RL}^{\left(2\right)}\left(\omega \right)=\underset{\theta =0}{\overset{2\pi }{.Math.}}{2}^8\pi k_0^4{.Math.\Gamma .Math.}^{2}S\left(m\stackrel{.fwdarw.}{k}\right)S\left(-2m^{\prime }\stackrel{.fwdarw.}{k_0}\right)\frac{{\left(2k_0\right)}^4}{k^{\mathrm{\prime 4}}}{y}^{*3}{.Math.\frac{\partial y}{\partial h}.Math.}_{\theta ,y=y^{*}}.Math.\mathrm{\Delta \theta },m=m^{\prime }=1\mathrm{or}-1$ where k.sub.0 is a radar wavenumber; k and k′ are the wavenumber of two sets of ocean waves which generate the second-order scattering with the radar wave vector {right arrow over (k.sub.0)}; θ is an angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; Δθ is a discrete interval of θ; S{right arrow over ((⋅))} is a directional wave spectrum; if m=m′=1, σ.sub.R.sup.(2)(ω) indicates the outer second-order spectrum on the positive first-order peak side; if m=m′=−1, σ.sub.R.sup.(2)(ω) indicates the outer second-order spectrum on the negative first-order peak side; ω is a Doppler frequency; Γ is a coupling coefficient; y, y* and h are intermediate variables defined for the convenience of calculation, wherein y=√{square root over (k)} and h+mg.sup.1/2y+m′g.sup.1/2√{square root over (k′)}, g is the acceleration of gravity; and y* is the solution of constant Doppler frequency contours denoted by ω−h=0; (2) obtaining A.sub.f (θ) which is defined as the single-frequency coefficient matrix of the directional wave spectrum by calculating the theoretical value of a function of the linearized outer second-order spectrum σ.sub.RL.sup.(2)(ω) divided by the stronger first-order peak energy σ.sub.R.sup.(1)(ω); $\frac{{\sigma }_{RL}^{\left(2\right)}\left(\omega \right)}{{\sigma }_R^{\left(1\right)}\left({\omega }_B\right)}=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)S\left(m\stackrel{.fwdarw.}{k}\right)\mathrm{\Delta \theta },\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}$ $\frac{{\sigma }_{RL}^{\left(2\right)}\left(\omega \right)}{{\sigma }_R^{\left(1\right)}\left(-{\omega }_B\right)}=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)S\left(m\stackrel{.fwdarw.}{k}\right)\mathrm{\Delta \theta },\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}$ where A.sub.f (θ) is the single-frequency coefficient matrix of the directional wave spectrum; f in the subscript is the radar frequency; S{right arrow over ((⋅))} is the directional wave spectrum; k is the wavenumber of a set of waves that generate second-order scattering with the radar wave vector; θ is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; Δθ is the discrete interval of θ; σ.sub.R.sup.(1)(ω.sub.B) and σ.sub.R.sup.(1)(−ω.sub.B) are the theoretical values of the first-order peaks; the theoretical value of the stronger first-order peak is given as:
σ.sub.R.sup.(1)(ω.sub.B)=2.sup.6πk.sub.0.sup.4S(−2{right arrow over (k.sub.0)}), when σ.sub.R.sup.(1)(ω) is the positive first-order peak
σ.sub.R.sup.(1)(−ω.sub.B)=2.sup.6πk.sub.0.sup.4S(2{right arrow over (k.sub.0)}), when σ.sub.R.sup.(1)(ω) is the negative first-order peak whereby the single-frequency coefficient matrix of the directional wave spectrum A.sub.j(θ) is given as: $A_f\left(\theta \right)={.Math.\Gamma .Math.}^2\frac{4{\left(2k_0\right)}^4{y}^{*3}}{{\left(k^{\prime }\right)}^4}{.Math.\frac{\partial y}{\partial h}.Math.}_{\theta ,y=y^{*}}$ where Γ is the coupling coefficient; k.sub.0 is the radar wavenumber; k′ is the wavenumber of a set of ocean waves which generate the second-order scattering with the radar wave vector; θ is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; y, y* and h are intermediate variables defined for the convenience of calculation, wherein y=√{square root over (k)} and h=mg.sup.1/2y+m′g.sup.1/2{right arrow over (k′)}, wherein g is the acceleration of gravity; and y* is the solution of constant Doppler frequency contours denoted by ω−h=0; (3) calculating A.sub.f which is defined as the single-frequency coefficient matrix of the nondirectional wave spectrum by discretizing θ into n pieces at equal intervals and summing all terms of θ in the following equation: $A_f=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)G\left(\theta \right)\mathrm{\Delta \theta }$ where A.sub.f(θ) is the single-frequency coefficient matrix of the directional wave spectrum; f in the subscript is the radar frequency; θ is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector ${\stackrel{.fwdarw.}{k}}_0;\mathrm{\Delta \theta }=\frac{2\pi }{n}$ is the discrete interval of θ; G(θ) is the direction distribution function of the directional wave spectrum.

6. The method of claim 1, wherein (f) comprises the step of: merging R.sub.f(ω) which is defined as the function of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum of each frequency and A.sub.f which is defined as the single-frequency coefficient matrix of the nondirectional wave spectrum into a matrix R and a matrix A respectively according to the following equation:
R=[R.sub.f.sub.1(ω)R.sub.f.sub.2(ω)R.sub.f.sub.3(ω)R.sub.f.sub.4(ω)].sup.T
A=[A.sub.f.sub.1A.sub.f.sub.2A.sub.f.sub.3A.sub.f.sub.4].sup.T where R.sub.f.sub.1(ω), R.sub.f.sub.2(ω), R.sub.f.sub.2(ω) and R.sub.f.sub.2(ω) are the functions of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum extracted from the sea echoes of radar frequencies f.sub.1, f.sub.2, f.sub.3 and f.sub.4, respectively; A.sub.f.sub.1, A.sub.f.sub.2, A.sub.f.sub.3, and A.sub.f.sub.4 are the coefficient matrices of the nondirectional wave spectrum from the sea echoes of radar frequencies f.sub.1, f.sub.2, f.sub.3 and f.sub.4, respectively.

7. The method of claim 1, wherein (g) comprises the steps of: calculating the pseudo-inverse A.sup.+ of the matrix A by the singular value decomposition; and estimating the nondirectional wave spectrum S(k) from the matrix R and the pseudo-inverse A.sup.+ by the equation S(k)=A.sup.+R.

8. The method of claim 1, being carried out in an apparatus, the apparatus comprising a multi-frequency HF radar and a computer processor: the multi-frequency HF radar comprising: a transmitter for transmitting electromagnetic wave signals with multiple frequencies between 3-30 MHz in the HF band to the sea surface; a receiver for receiving sea echoes of multiple radar frequencies originated from the electromagnetic wave signals modulated by sea surface movements; and a signal processor unit for obtaining each Doppler spectrum from the sea echo of each radar frequency at a fan-shaped unit by performing the first FFT in distance dimension, the second FFT in Doppler frequency dimension and the digital beamforming; the computer processor providing code segments for: (1) extracting a positive first-order peak and a negative first-order peak from a Doppler spectrum of a single radar frequency by a peak-searching method and selecting a stronger first-order peak σ.sub.R.sup.(1)(ω) from the positive and negative first-order peaks; (2) dividing a second-order spectrum on a stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum and separating the outer second-order spectrum on the stronger first-order peak side from the Doppler spectrum according to a Doppler frequency range of the outer second-order spectrum; (3) calculating R.sub.f(ω) which is a function of the outer second-order spectrum divided by a first-order peak energy for the stronger side of the Doppler spectrum of the single radar frequency at the fan-shaped unit; (4) calculating A.sub.f which is defined as a single-frequency coefficient matrix of the nondirectional wave spectrum by linearizing the outer second-order spectrum on the side of the stronger first-order peak; (5) repeating the above steps (1) to (4), combining the sea echoes of different radar frequencies to construct a matrix R which is a function of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum of multiple radar frequencies and merging coefficient matrices of the nondirectional wave spectrum from multiple radar frequencies into a matrix A; and (6) calculating a pseudo-inverse A.sup.+ of the matrix A by a singular value decomposition and estimating the nondirectional wave spectrum at the fan-shaped unit from the matrix R and the pseudo-inverse A.sup.+.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0054] The disclosure will be more clearly understood from the following description when read in conjunction with the accompanying drawings of which:

[0055] FIG. 1 is a flow chart showing the method for directly estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies according to an embodiment of the disclosure.

[0056] FIGS. 2A-2D are graphs showing the comparison between the nondirectional wave spectrum estimated from the simulated sea echoes of a single radar frequency and the theoretical wave spectrum according to an embodiment of the disclosure.

[0057] FIG. 3 is a graph showing the comparison between the nondirectional wave spectrum estimated from the simulated sea echoes of multiple radar frequencies and the theoretical wave spectrum according to an embodiment of the disclosure.

[0058] FIGS. 4A-4C are graphs showing the comparison among the nondirectional wave spectra estimated from the sea echoes collected by a multi-frequency HF radar system according to an embodiment of the disclosure.

[0059] FIG. 5 shows a block diagram of an apparatus for carrying out the method for directly estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies according to an embodiment of the disclosure.

DETAILED DESCRIPTION

[0060] FIG. 5 shows a block diagram of an apparatus for carrying out a method for estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies. The apparatus comprises a multi-frequency HF radar and a computer processor. The multi-frequency HF radar comprises a transmitter, a receiver and a signal processor unit. The transmitter transmits electromagnetic wave signals of multiple frequencies between 3-30 MHz in the HF band to the sea surface. The receiver receives sea echoes from multiple radar frequencies originated from the electromagnetic wave signals modulated by ocean surface movements. The signal processor unit obtains each Doppler spectrum from the sea echo of each radar frequency at a fan-shaped unit by performing a first fast Fourier transform (FFT) in distance dimension, a second FFT in Doppler frequency dimension and a digital beamforming. The computer processor provides code segments for carrying out the method for estimating the nondirectional wave spectrum. The transmitter and receiver are respectively connected to the signal processor unit through communication links A and B to realize the full-duplex communication. The signal processor unit is connected to the computer processor through the communication link C to realize the full-duplex communication. Communication links A and B include but are not limited to a wireless connection and a hardwired connection. Communication link C is a network connection, such as Transmission Control Protocol/Internet Protocol (TCP/IP). The multi-frequency HF radar and the computer processor may be installed far away from each other.

[0061] In operation, the signal processor unit may control the transmitter to transmit frequency modulated interrupted continuous wave (FMICW) chirps to the ocean surface through the link A. In a Doppler sampling period, the transmitter may transmit FMICW chirps with different frequencies in chronological order, and then the receiver receives sea echoes modulated by ocean surface movements with corresponding frequencies. The receiver communicates with the signal processor unit through the link B and sends the received sea echo to the signal processor unit. The signal processor unit performs the first FFT in distance dimension and the second FFT in Doppler frequency dimension and the digital beamforming to obtain the Doppler spectrum in a fan-shaped unit. The signal processor unit sends the Doppler spectrum in a fan-shaped unit to the computer processor through the link C, and then the code segments for estimating the nondirectional wave spectrum provided by the computer processor are executed.

[0062] FIG. 1 show the outlined steps of the method for estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies, comprising the following steps (a)-(g).

[0063] The step (a) is dividing a radar detection area into a plurality of fan-shaped units at an equal range interval and angle interval according to a distance resolution and an angular resolution of a multi-frequency HF radar, wherein the multi-frequency HF radar is capable of simultaneously operating at more than one frequency in the HF band.

[0064] The step (b) is obtaining a Doppler spectrum from a sea echo of a single radar frequency at a fan-shaped unit by performing a first fast FFT in distance dimension, a second FFT in Doppler frequency dimension and a digital beamforming in a signal processor unit, extracting a positive first-order peak and a negative first-order peak from the Doppler spectrum by a peak-searching method and selecting a stronger first-order peak σ.sub.R.sup.(1)(ω) from the positive and negative first-order peaks in a computer processor.

[0065] In the step (b), the Doppler spectrum at the fan-shaped unit is defined as: σ(ω), where ω represents a Doppler frequency generated by the motion of ocean waves to the multi-frequency HF radar; and σ(ω) represents a wave energy distribution at the different value of ω.

[0066] The positive first-order peak and negative first-order peak are extracted from the Doppler spectrum at the fan-shaped unit using the peak searching method. The first-order peaks are defined as two peaks in the Doppler spectrum which are roughly symmetrically distributed on both sides of zero frequency. The first-order peaks are generated by the Bragg scattering of the waves of half the radar wavelength that are either advancing directly towards the multi-frequency HF radar or receding directly from the multi-frequency HF radar The detailed steps of the step (b) are as follows:

[0067] searching for a point of a maximum amplitude in the Doppler frequency range [0.6ω.sub.B,1.4ω.sub.B] of the Doppler spectrum at the fan-shaped unit as the peak of the positive first-order peak, and recording the Doppler frequency of the peak of the positive first-order peak as ω.sub.P+, wherein ω.sub.B=√{square root over (2gk.sub.0)} is a Bragg frequency, wherein

[00011]$k_0=\frac{2\pi f}{c}$

is a radar wavenumber; c is a speed of light; and f is a radar frequency;

[0068] searching for a local minimum point inside the peak of the positive first-order peak where the Doppler frequency meets the inequation ω.sub.P+−0.2ω.sub.B≤ω<ω.sub.P+, and denoting the Doppler frequency of the local minimum point inside the peak of the positive first-order peak as ω.sub.L+;

[0069] searching for the local minimum point outside the peak of the positive first-order peak where the Doppler frequency satisfies the inequation ω.sub.P+<(ω) 0.2ω.sub.B+ω.sub.P+, and recording the Doppler frequency of the local minimum point outside the peak of the positive first-order peak as ω.sub.R+;

[0070] intercepting the Doppler spectrum at the fan-shaped unit with Doppler frequency [ω.sub.L+,ω.sub.R+] as the positive first-order peak;

[0071] searching for the point of the maximum amplitude in the Doppler frequency range [−1.4ω.sub.B,−0.6ω.sub.B] of the Doppler spectrum at the fan-shaped unit as the peak of the negative first-order peak, and recording the Doppler frequency of the peak of the negative first-order peak as ω.sub.P−;

[0072] searching for the local minimum point inside the peak of the negative first-order peak where the Doppler frequency meets the inequation ω.sub.P−<ω≤ω.sub.P−+0.2ω.sub.B, and denoting the Doppler frequency of the local minimum point inside the peak of the negative first-order peak as ω.sub.L−;

[0073] searching for the local minimum point outside the peak of the negative first-order peak where the Doppler frequency satisfies the inequation ω.sub.P− −0.2ω.sub.B≤ω<ω.sub.P−, and recording the Doppler frequency of the local minimum point outside the peak of the negative first-order peak as ω.sub.R−;

[0074] intercepting the Doppler spectrum at the fan-shaped unit with Doppler frequency [ω.sub.R−, ω.sub.L−] as the negative first-order peak; and

[0075] comparing the amplitude of the peak of the positive first-order peak with that of the negative first-order peak, and selecting the first-order peak with a larger amplitude of the peak point as the stronger first-order peak σ.sub.R.sup.(1)(ω).

[0076] The step (c) is dividing a second-order spectrum on the stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum, and separating the outer second-order spectrum on the stronger first-order peak side from the Doppler spectrum in the computer processor according to a Doppler frequency range of the outer second-order spectrum;

[0077] In the step (c), the second-order spectrum originated from the second-order scattering of ocean waves and radar waves is a continuum with lower amplitude than the first-order peaks and distributed around the first-order peaks.

[0078] The step (c) comprises the following sub-steps:

[0079] dividing the second-order spectrum on the stronger first-order peak side into an inner second-order spectrum and an outer second-order spectrum;

[0080] The Doppler frequency range of the outer second-order spectrum on the stronger first-order peak side is given as

[00012]$\{\begin{array}{c}{\omega }_{c+}<\omega \le 14{\omega }_B,\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ -1.4{\omega }_B\le \omega <{\omega }_{c-},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array}$

[0081] The Doppler frequency range of the inner second-order spectrum on the stronger first-order peak side is given as

[00013]$\{\begin{array}{c}0.6{\omega }_B\le \omega <{\omega }_{c+},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ {\omega }_{c-}<\omega \le -0.6{\omega }_B,\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array}$

[0082] where ω is a Doppler frequency generated by the motion of ocean waves to the multi-frequency HF radar; ω.sub.B=√{square root over (2gk.sub.0)} is a Bragg frequency; k.sub.0 is the radar wavenumber;

[00014]${\omega }_{c+}=\frac{\underset{{\omega }_{L+}}{\stackrel{{\omega }_{R+}}{.Math.}}\omega .Math.{\sigma }_R^{\left(1\right)}\left(\omega \right)}{\underset{{\omega }_{L+}}{\stackrel{{\omega }_{R+}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)}$

is a centroid frequency of the positive first-order peak; and

[00015]${\omega }_{c-}=\frac{\underset{{\omega }_{R-}}{\stackrel{{\omega }_{L-}}{.Math.}}\omega .Math.{\sigma }_R^{\left(1\right)}\left(\omega \right)}{\underset{{\omega }_{R-}}{\stackrel{{\omega }_{L-}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)}$

is a centroid frequency of the negative first-order peak.

[0083] if the stronger first-order peak obtained in the step (b) is the positive first-order peak, taking the centroid frequency ω.sub.c+ as a starting point of Doppler frequency, taking a cutoff frequency of 1.4ω.sub.B as an endpoint, extracting the Doppler spectrum σ(ω) at the fan-shaped unit in the frequency range of (ω.sub.c+,1.4ω.sub.B], and setting the Doppler spectrum σ(ω) at the fan-shaped unit in the frequency range of (ω.sub.c+,ω.sub.R+] to zero to obtain the outer second-order spectrum on the stronger first-order peak side σ.sub.R.sup.(2)(ω);

[0084] if the stronger first-order peak obtained in the step (b) is the negative first-order peak, taking the centroid frequency ω.sub.c− as the starting point of Doppler frequency, taking a cutoff frequency of −1.4ω.sub.B as the endpoint, extracting the Doppler spectrum σ(ω) at the fan-shaped unit in the frequency range of [−1.4ω.sub.B,ω.sub.c−), and setting the Doppler spectrum σ(ω) at the fan-shaped unit in the frequency range of [ω.sub.R−,ω.sub.c−) to zero to obtain the outer second-order spectrum on the stronger first-order peak side σ.sub.R.sup.(2)(ω).

[0085] The step (d) is calculating R.sub.f(ω) which is a function of the outer second-order spectrum divided by a first-order peak energy for the stronger side of the Doppler spectrum of the single radar frequency at the fan-shaped unit in the computer processor. The function R.sub.f(ω) can be given as

[00016]$\{\begin{array}{c}{R}_f\left(\omega \right)=\frac{{\sigma }_R^{\left(2\right)}\left(\omega \right)}{\underset{{\omega }_{L+}}{\stackrel{{\omega }_{R+}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{\Delta \omega }},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}\\ R_f\left(\omega \right)=\frac{{\sigma }_R^{\left(2\right)}\left(\omega \right)}{\underset{{\omega }_{R-}}{\stackrel{{\omega }_{L-}}{.Math.}}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{\Delta \omega }},\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}\end{array}$

[0086] where f in a subscript indicates a radar frequency; σ.sub.R.sup.(1)(ω) is the stronger first-order peak; σ.sub.R.sup.(2)(ω) is the outer second-order spectrum on the stronger first-order peak side; ω.sub.L+ is the Doppler frequency of the local minimum point inside the peak point of the positive first-order peak; ω.sub.R+ is the Doppler frequency of the local minimum point outside the peak point of the positive first-order peak; ω.sub.L− is the Doppler frequency of the local minimum point inside the peak point of the negative first-order peak; φ.sub.R− is the Doppler frequency of the local minimum point outside the peak point of the negative first-order peak; and Δω is a frequency resolution of the Doppler spectrum σ(ω).

[0087] The step (e) is calculating the single-frequency coefficient matrix of the nondirectional wave spectrum by linearizing the outer second-order spectrum on the side of the stronger first-order peak in the computer processor.

[0088] The step of calculating the single-frequency coefficient matrix of the nondirectional wave spectrum comprises:

[0089] linearizing the outer second-order spectrum on the side of the stronger first-order peak;

[0090] at the Doppler frequencies near the first-order peak that satisfy the condition ω.sub.B<|ω|≤1.4ω.sub.B, the direction of one of the two sets of ocean waves that generate the second-order scattering with the radar vector is approximately equal to that of the Bragg wave vector, so the outer second-order spectrum σ.sub.R.sup.(2)(ω) can be linearized as σ.sub.RL.sup.(2)(ω)

[00017]${\sigma }_{RL}^{\left(2\right)}\left(\omega \right)=\underset{\theta =0}{\overset{2\pi }{.Math.}}{2}^8\pi k_0^4{.Math.\Gamma .Math.}^{2}S\left(m\stackrel{.fwdarw.}{k}\right)S\left(-2m^{\prime }\stackrel{.fwdarw.}{k_0}\right)\frac{{\left(2k_0\right)}^4}{k^{\mathrm{\prime 4}}}{y}^{*3}{.Math.\frac{\partial y}{\partial h}.Math.}_{\theta ,y=y^{*}}.Math.\mathrm{\Delta \theta },m=m^{\prime }=1\mathrm{or}-1$

[0091] where

[00018]$k_0=\frac{2\pi f}{c}$

is a radar wavenumber; k=(0,0.3] and k′=(k.sup.2+4kk.sub.0 cosθ+(2k.sub.0).sup.2).sup.1/2 are the wavenumber of two sets of ocean waves which generate the second-order scattering with the radar wave vector {right arrow over (k.sub.0)}; θ=[0,2π] is an angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; Δθ is a discrete interval of θ; S{right arrow over ((⋅))} is a directional wave spectrum; if m=m′=1, σ.sub.R.sup.(2)(ω) indicates the outer second-order spectrum on the positive first-order peak side; if m=m′=−1, σ.sub.R.sup.(2)(ω) indicates the outer second-order spectrum on the negative first-order peak side;

[00019]$\Gamma =\frac{1}{2}\left[\frac{\left(\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}_0\right)\left({\stackrel{.fwdarw.}{k}}^{\prime }.Math.{\stackrel{.fwdarw.}{k}}_0\right)/k_0^2-2\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }}{\sqrt{\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }}-k_0\Delta }\right]-\frac{i}{2}\left[k+k^{\prime }-\frac{\left({\mathrm{kk}}^{\prime }-\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }\right)\left({\omega }^2+{\omega }_B^2\right)}{{\mathrm{mm}}^{\prime }\sqrt{{\mathrm{kk}}^{\prime }}\left({\omega }^2-{\omega }_B^2\right)}\right]$

is a coupling coefficient, wherein m, m′, k.sub.0, k, k′, {right arrow over (k)} and {right arrow over (k.sub.0)} are the same as described above, {right arrow over (k)} is one of two sets of ocean waves which generate the second-order scattering with the radar wave vector {right arrow over (k.sub.0)}, ω is the Doppler frequency, ω.sub.B=√{square root over (2gk.sub.0)} is the Bragg frequency, i is an imaginary unit and Δ is the normalized impedance of sea surface; y, y* and h are intermediate variables defined for the convenience of calculation, wherein y=√{square root over (k)} and h=mg.sub.1/2y m′g.sup.1/2√{square root over (k′)}, wherein g is the acceleration of gravity; and y* is the solution of constant Doppler frequency contours denoted by ω−h=0.

[0092] (2) obtaining A.sub.f(θ) which is defined as the single-frequency coefficient matrix of the directional wave spectrum by calculating the theoretical value of a function of the linearized outer second-order spectrum σ.sub.R.sup.(2)(ω) divided by the stronger first-order peak energy σ.sub.R.sup.(1)(ω).

[00020]$\frac{{\sigma }_{RL}^{\left(2\right)}\left(\omega \right)}{{\sigma }_R^{\left(1\right)}\left({\omega }_B\right)}=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)S\left(m\stackrel{.fwdarw.}{k}\right)\mathrm{\Delta \theta },\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{positive}\mathrm{first}-\mathrm{order}\mathrm{peak}$$\frac{{\sigma }_{RL}^{\left(2\right)}\left(\omega \right)}{{\sigma }_R^{\left(1\right)}\left(-{\omega }_B\right)}=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)S\left(m\stackrel{.fwdarw.}{k}\right)\mathrm{\Delta \theta },\mathrm{when}{\sigma }_R^{\left(1\right)}\left(\omega \right)\mathrm{is}\mathrm{the}\mathrm{negative}\mathrm{first}-\mathrm{order}\mathrm{peak}$

[0093] where A.sub.f(θ) is the single-frequency coefficient matrix of the directional wave spectrum; f in the subscript is the radar frequency; S{right arrow over ((⋅))} is the directional wave spectrum; k is the wavenumber of a set of waves that generate second-order scattering with the radar wave vector; θ=[0,2π] is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; Δθ is the discrete interval of θ; σ.sub.R.sup.(1)(ω.sub.B) or σ.sub.R.sup.(1)(−ω.sub.B) is the theoretical value of the first-order peaks. The theoretical value of the stronger first-order peak is given as:

σ.sub.R.sup.(1)(ω.sub.B)=2.sup.6πk.sub.0.sup.4S(−2{right arrow over (k.sub.0)}), when σ.sub.R.sup.(1)(ω) is the positive first-order peak

σ.sub.R.sup.(1)(−ω.sub.B)=2.sup.6πk.sub.0.sup.4S(2{right arrow over (k.sub.0)}), when σ.sub.R.sup.(1)(ω) is the negative first-order peak

[0094] whereby the single-frequency coefficient matrix of the directional wave spectrum A.sub.j(θ) is given as:

[00021]$A_f\left(\theta \right)={.Math.\Gamma .Math.}^2\frac{4{\left(2k_0\right)}^4{y}^{*3}}{{\left(k^{\prime }\right)}^4}{.Math.\frac{\partial y}{\partial h}.Math.}_{\theta ,y=y^{*}}$

[0095] where

[00022]$\Gamma =\frac{1}{2}\left[\frac{\left(\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}_0\right)\left({\stackrel{.fwdarw.}{k}}^{\prime }.Math.{\stackrel{.fwdarw.}{k}}_0\right)/k_0^2-2\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }}{\sqrt{\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }}-k_0\Delta }\right]-\frac{i}{2}\left[k+k^{\prime }-\frac{\left({\mathrm{kk}}^{\prime }-\stackrel{.fwdarw.}{k}.Math.{\stackrel{.fwdarw.}{k}}^{\prime }\right)\left({\omega }^2+{\omega }_B^2\right)}{{\mathrm{mm}}^{\prime }\sqrt{{\mathrm{kk}}^{\prime }}\left({\omega }^2-{\omega }_B^2\right)}\right]$

is the coupling coefficient;

[00023]$k_0=\frac{2\pi f}{c}$

is the radar wavenumber; k′=(k.sup.2+4kk.sub.0 cos θ+(2k.sub.0).sup.2).sup.1/2 is the wavenumber of a set of ocean waves which generate the second-order scattering with the radar wave vector; θ−[0,2π] is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector {right arrow over (k.sub.0)}; y, y* and h are intermediate variables defined for the convenience of calculation, wherein y=√{square root over (k)} and h=mg.sup.1/2y+m′g.sup.1/2√{square root over (k′)}, wherein g is the acceleration of gravity; and y* is the solution of constant Doppler frequency contours denoted by ω−h=0.

[0096] (3) calculating A.sub.f which is defined as the single-frequency coefficient matrix of the nondirectional wave spectrum by discretizing θ into n pieces at equal intervals and summing all terms of θ in the following equation:

[00024]$A_f=\underset{\theta =0}{\overset{2\pi }{.Math.}}{A}_f\left(\theta \right)G\left(\theta \right)\mathrm{\Delta \theta }$

[0097] where A.sub.f(θ) is the single-frequency coefficient matrix of the directional wave spectrum; f in the subscript is the radar frequency; θ−[0,2π] is the angle between the ocean wave vector {right arrow over (k)} and the radar wave vector

[00025]${\stackrel{.fwdarw.}{k}}_0;\mathrm{\Delta \theta }=\frac{2\pi }{n}$

is the discrete interval of θ; G(θ) is the direction distribution function of the directional wave spectrum. The cardioid distribution function

[00026]$G\left(\theta \right)={\cos }^{2s}\left(\frac{\theta -{\theta }^{*}}{2}\right)/{\int }_{-\pi }^{\pi }{\cos }^{2s}\left(\frac{\theta }{2}\right)d\theta$

can be selected, where θ* is the wind direction and s=2 is the spread factor.

[0098] The step (f) is repeating the above steps (b) to (e); combining the sea echoes of different radar frequencies to construct a matrix R which is a function of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum of multiple radar frequencies and merging coefficient matrices of the nondirectional wave spectrum from multiple radar frequencies into a matrix A in the computer processor.

[0099] The detailed process of the step (f) is as follows:

[0100] R.sub.f(ω) which is defined as the function of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum of each frequency and A.sub.f which is defined as the single-frequency coefficient matrix of the nondirectional wave spectrum are merged into a matrix R and a matrix A respectively according to the following equation:

R=[R.sub.f.sub.1(ω)R.sub.f.sub.2(ω)R.sub.f.sub.3(ω)R.sub.f.sub.4(ω)].sup.T

A=[A.sub.j.sub.1A.sub.f.sub.2A.sub.f.sub.3A.sub.f.sub.4].sup.T

[0101] where R.sub.f.sub.1(ω), R.sub.f.sub.2(ω), R.sub.f.sub.2(ω) and R.sub.f.sub.2(ω) are the functions of the outer second-order spectrum divided by the first-order peak energy for the stronger side of the Doppler spectrum extracted from the sea echoes of radar frequencies f.sub.1, f.sub.2, f.sub.3 and f.sub.4, respectively; and A.sub.f.sub.1, A.sub.f.sub.2, A.sub.f.sub.3, and A.sub.f.sub.4 are the coefficient matrices of the nondirectional wave spectrum from the sea echoes of radar frequencies fj, f.sub.3 and f.sub.4 respectively.

[0102] The step (g) is calculating the pseudo-inverse A.sup.+ of the matrix A by a singular value decomposition; and estimating the nondirectional wave spectrum at the fan-shaped unit from the matrix R and the pseudo-inverse A.sup.+ in the computer processor.

[0103] The aforementioned nondirectional wave spectrum in the step (g) is given by the equation:

S(k)=A.sup.+R

[0104] Because the matrix A is not square, only the pseudo-inverse of the matrix A can be calculated. The singular value decomposition can be carried out on the matrix A to obtain its pseudo-inverse, and then the nondirectional wave spectrum S (k) can be estimated from the matrix R and the matrix A.sup.+.

[0105] As an example, the Pierson-Moskowitz (PM) spectrum and the cardioid distribution function are selected to simulate the Doppler spectrum of the sea echoes of radar frequency 8 MHz, 13 MHz, 19 MHz and 25 MHz respectively. The wind speed is set to 10 m/s, and the wind direction is set to 135°. The comparison between the nondirectional wave spectrum estimated from the simulated sea echoes of a single radar frequency and the theoretical wave spectrum is shown in FIGS. 2A-2D. As shown in FIG. 3 the nondirectional wave spectrum estimated from the sea echoes of multiple radar frequencies is in good agreement with that of the theoretical values. FIGS. 2A-2D and FIG. 3 demonstrate that the estimation result of the nondirectional wave spectrum from the multiple frequencies provides reasonable agreement with the theoretical values compared with that from the sea echoes at 8 MHz, 13 MHz, 19 MHz, and 25 MHz.

[0106] The method provided by the disclosure is applied to a dataset collected by a multi-frequency HF radar system operating at 8.267 MHz and 19.2 MHz. The wave height estimated by the method from the radar data is compared with that measured by a buoy to evaluate the performance of the method. The comparison among the nondirectional wave spectra estimated from the sea echoes collected with a multi-frequency HF radar system at 08:00 AM on Jul. 10, 2015 is shown in FIGS. 4A-4C. The wave height extracted by the nondirectional wave spectrum directly estimated from the multifrequency data is more consistent with that measured by the buoy than that estimated from the single-frequency data.

[0107] The method provided by the disclosure for directly estimating the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies breaks through the limitation of the sea states that can be detected at a fixed radar frequency by combining echo signals from multiple radar frequencies, whereby the robustness and accuracy of estimating the nondirectional wave spectrum are greatly enhanced and the potential of the method for monitoring the complicated and changeable sea surface is recognized. Therefore, it is feasible to estimate the nondirectional wave spectrum from the sea echoes of multiple HF radar frequencies, and it has higher accuracy and wide application prospect compared with the conventional method for inverting the nondirectional wave spectrum from the single-frequency sea echo.

[0108] It will be obvious to those skilled in the art that changes and modifications may be made, and therefore, the aim in the appended claims is to cover all such changes and modifications.