收藏切换
The multifractal spectrum of a sea clutter using a random walk model
收藏切换
PDF
Jingbo HE1, *, Jianghu XU1
Acta Oceanologica Sinica | 2017, 36(9) : 23 - 26
Less
收藏切换
Acta Oceanologica Sinica | 2017, 36(9): 23-26
The multifractal spectrum of a sea clutter using a random walk model
Full
Jingbo HE1, *, Jianghu XU1
Affiliations
  • 1 Electronics Engineering College, Naval University of Engineering, Wuhan 430033, China
Published: 2017-09-01 doi: 10.1007/s13131-017-1107-y
Outline
收藏切换

The radar echo signal of a sea clutter is nonlinear, nonstationary and time varying. A multifractal measure analysis can describe the local singularity of a physics system. The random walk model of a sea clutter scattering is analysed to disclose the intrinsic physical characteristics and laws of the sea clutter. Stochastic differential equations are given for the physical quality of the sea clutter. A diffusion process model is established using $ I t {\mathop o \limits^ \frown } $ formula. The singularity of the random walk model is tested by a multifractal spectroscopy, and the accuracy of this model is proven by the multifractal spectroscopy of a real-life IPIX radar data set. Thus, the random walk model is effective for describing the dynamics mechanism of the sea clutter.

random walk  /  sea clutter  /  multifractal
Jingbo HE, Jianghu XU. The multifractal spectrum of a sea clutter using a random walk model[J]. Acta Oceanologica Sinica, 2017 , 36 (9) : 23 -26 . DOI: 10.1007/s13131-017-1107-y
The sea clutter refers to the surface of a radar-scattering echo and is the most complicated form of a radar clutter. After years of testing, researchers have found that the probability distribution of the sea clutter is a Gaussian model. They have proposed a logarithmic normal distribution, Weibull distribution, K-distribution and other distribution models (Rosenberg and Bocquet, 2015; Velotto et al., 2014; Güntürkün, 2015; Jakeman and Pusey, 1976; Li et al., 2014). These models are close to the sea clutter probability distribution form to some extent, but they do not clearly describe the generating mechanism of the sea clutter. Only the statistical models for the first-order and second-order characteristics of the sea clutter are described. The time-varying characteristic of the sea clutter is not depicted (Granström et al., 2015; Suresh et al., 2015; Xing et al., 2014; Xiong et al., 2014). Consequently, precise modelling, analysis and processing of the nonlinear, nonstationary and time-varying characteristics of the sea clutter are difficult (Wu et al., 2014; Zhang et al., 2014; Xu et al., 2014).
The multifractal measure analysis of different physical systems can completely describe the local singularity (Guan et al., 2010; Liu et al., 2012); the statistical properties of such singular measure in the characterisation can reveal complex heterogeneous structure system by determining the singularity spectrum. Jakeman and Tough (1987) proposed electromagnetic scattering of the random walk model to reveal the dynamics mechanism of the sea clutter. Considering the different parameters of the model, the evolution of the sea clutter can be all types of a statistical distribution model (e.g., Rayleigh distribution, Weibull distribution and K-distribution). The random walk model must accurately depict the time-varying characteristic of the sea clutter. On this basis, it utilises the random walk model to describe the time-varying electromagnetic scattering characteristics of the sea clutter using the stochastic differential equation model.
The electromagnetic scattering signal of the sea clutter is derived by the stochastic differential equation based on a random walk. According to Jakeman random walk theory, the electromagnetic scattering model of the sea clutter (Jakeman and Tough, 1987; Feng et al., 2007) can be expressed as
$\begin{array}{*{20}{l}} {{\text{Ψ} _t}} & { = \mathop {\lim }\limits_{{N_t} \to \infty } \left\{ {\frac{1}{{{{(\overline N )}^{1/2}}}}\sum\limits_{j = 1}^{{N_t}} {{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} } \right\}}\\& { = \mathop {\lim }\limits_{{N_t} \to \infty } \left\{ {{{\left( {\frac{{{N_t}}}{{\overline N }}} \right)}^{1/2}}\frac{1}{{{{({N_t})}^{1/2}}}}\sum\limits_{j = 1}^{{N_t}} {{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} } \right\}}\\ & { = {{({x_t})}^{1/2}}{\gamma _t}}\,\,\,\,\, ,\end{array}$
with $\phi _t^{(j)} = {\Delta ^{(j)}} + {B^{1/2}}W_t^{(j)}$. ${\Delta ^{(j)}}$ obeys a uniform random distribution, and $W_t^{(j)}$ denotes the Wiener process. ${x_t} = \mathop {\lim }\limits_{{N_t} \to \infty } \left[ {{N_t}/\overline N} \right]$; and ${\gamma _t} = \mathop {\lim }\limits_{{N_t} \to \infty } \left[ {\sum\limits_{j = 1}^{{N_t}} {{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} /{{({N_t})}^{1/2}}} \right]$. In studying the stochastic differential equation of γt, the following definition is given:
$\varepsilon _t^{({N_t})} = \sum\limits_{j = 1}^{{N_t}} {{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} .$
According to the $ I t {\mathop o \limits^ \frown } $ formula,
$\begin{array}{*{20}{l}}{{\rm{d}}\varepsilon _t^{(N)}} & { = \sum\limits_{j = 1}^N {[{\rm{id}}\phi _t^{(j)} - \frac{1}{2}{{({\rm{d}}\phi _t^{(j)})}^2}]{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} }\\{} & { = \sum\limits_{j = 1}^N {[{\rm{i}}{B^{1/2}}{\rm{d}}W_t^{(j)} - \frac{1}{2}B{\rm{d}}t]{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} }\\{} & { = \sum\limits_{j = 1}^N {{\rm{i}}{B^{1/2}}{\rm{d}}W_t^{(j)}{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} - \frac{1}{2}B{\rm{d}}t\sum\limits_{j = 1}^N {{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} }\\{} & { = V - \frac{1}{2}B\varepsilon _t^{(N)}{\rm{d}}t},\end{array}$
where
$\begin{array}{*{20}{l}}V & { = \sum\limits_{j = 1}^N {{\rm{i}}{B^{1/2}}{\rm{d}}W_t^{(j)}{{\rm{e}}^{{\rm{i}}\phi _t^{(j)}}}} }\\{} & { = {B^{1/2}}[{\rm{i}}\sum\limits_{j = 1}^N {\cos \phi _t^{(j)}{\rm{d}}W_t^{(j)}} - \sum\limits_{j = 1}^N {\sin \phi _t^{(j)}{\rm{d}}W_t^{(j)}} ]}\\{} & { = {B^{1/2}}[{\rm{i}}{\sigma _c}{\rm{d}}W_t^{(c)} - {\sigma _s}{\rm{d}}W_t^{(s)}]}\\{} & { = {B^{1/2}}\sigma [{\rm{i}}\frac{{{\sigma _c}}}{\sigma }{\rm{d}}W_t^{(c)} - \frac{{{\sigma _s}}}{\sigma }{\rm{d}}W_t^{(s)}]}.\end{array}$
When $N \to \infty $, ${\sigma _c}/\sigma = {\sigma _s}/\sigma = 1/\sqrt 2 $, the Wiener process is defined as
${\rm{d}}{\xi _t} = \frac{1}{{\sqrt 2 }}({\rm{id}}W_t^{(c)} - {\rm{d}}W_t^{(s)}),$
which satisfies
$\left\{ {\begin{array}{*{20}{l}}{|{\rm{d}}\xi {|^2} = {\rm{d}}t}\\{{{({\rm{d}}\xi )}^2} = 0}\end{array}} \right..$
The resulting equation is
${\rm{d}}{\gamma _t} = - \frac{1}{2}B{\gamma _t}{\rm{d}}t + {B^{\frac{1}{2}}}{\rm{d}}{\xi _t}.$
According to a birth-death-immigration model (Wu et al., 2014), satisfies
${\rm{d}}{x_t} = {\rm A}(\text{α} - {x_t}){\rm{d}}t + {(2{\rm A}{x_t})^{1/2}}{\rm{d}}W_t^{(x)},$
where α and A are constants.
Let ${r_t} = {({x_t})^{1/2}}$, we can obtain
${\rm{d}}{r_t} = {\rm A}\left[ {\frac{{2(\text{α} - {x_t}) - 1}}{{4{r_t}}}} \right]{\rm{d}}t + {\left(\frac{\rm A}{2}\right)^{1/2}}{\rm{d}}W_t^{(x)},$
where rt is satisfied using the generated stochastic differential equations.
In polar coordinates ${\varPsi _t} = {R_t}{{\rm{e}}^{{\rm{i}}{\theta _t}}}$, according to the theory of stochastic differential, ${\rm{d}}({X _t}{Y _t}) = {X _t}{\rm{d}}{Y _t} + {Y _t}{\rm{d}}{X _t} + {\rm{d}}{X _t}{\rm{d}}{Y _t}$; thus,
${\rm{d}}{z_t} = {\varPsi _t}^*{\rm{d}}{\varPsi _t} + {\varPsi _t}{\rm{d}}{\varPsi _t}^* + {\left| {{\rm{d}}{\varPsi _t}} \right|^2}.$
According to Eq. (1),
$\begin{array}{*{20}{l}}{{\varPsi _t}^*{\rm{d}}{\varPsi _t} + {\varPsi _t}{\rm{d}}{\varPsi _t}^*}{ = {x_t}({\gamma _t}^*{\rm{d}}{\gamma _t} + {\gamma _t}{\rm{d}}{\gamma _t}^*) + 2|{\gamma _t}{|^2}{r_t}{\rm{d}}{r_t}}\\\,\,\, \quad \quad { = - B{z_t}{\rm{d}}t + {B^{1/2}}{x_t}({\gamma _t}^*{\rm{d}}{\xi _t} + {\gamma _t}{\rm{d}}{\xi _t}^*) + \frac{{2{z_t}}}{{{r_t}}}{\rm{d}}{r_t}}.\end{array}$
Combining Eqs (9)–(11), yields
$\begin{array}{l}{\rm{d}}{z_t} = \left[ {B({x_t} - {z_t}) + \frac{{A{z_t}(\alpha - {x_t})}}{{{x_t}}}} \right]{\rm{d}}t\\ \qquad + {\left( {2B{x_t}{z_t} + \frac{{2Az_t^2}}{{{x_t}}}} \right)^{1/2}}{\rm{d}}W_t^{(z)}.\end{array}$
We can then obtain the sea clutter scattering signal energy part (square amplitude) of the stochastic differential equation.
IPIX radar sea clutter data are provided by the McMaster University and have become a global study of the theory of a sea clutter standard database. The database records all types of a heavy sea condition, all types of a radar frequency and different polarisation modes of the sea clutter data. The database is thus suitable to be used for the in-depth study of the sea clutter. A number of studies have cited the database depicted in Fig. 1.
Using the public IPIX on the coherent radar sea clutter data by the Canadian McMaster University, this study selects 100 sea clutter data records for the experiments; these data records are shown in Fig. 2.
The multifractal spectrum is analysed, and the result is shown in Fig. 3. From the diagram, the IPIX coherent radar-measured sea clutter data of the minimum fractal index fall in between (0.6, 1.0), and the maximum fractal index falls in between (1.0, 1.4).
The multifractal spectrum analysis must also be conducted to validate the sea clutter stochastic differential equation model. According to Eq. (12) types of the sea clutter simulation data, the parameter settings are as follows: B=1, dt=0.001 and σ=1. The 100 data records are shown in Fig. 4.
The multifractal spectrum analysis results are shown in Fig. 5. The diagram indicates that the sea clutter random walk model data yield a minimum fractal index between (0.6, 1.0) and a maximum fractal index mainly between (1.0, 1.4), with consistent IPIX radar-measured sea clutter data of the multifractal spectrum. Thus, the stochastic differential equation based on the random walk can effectively describe the dynamics mechanism of the sea clutter.
The sea clutter is nonlinear, nonstationary and time varying. The multifractal measure analysis is the local singularity of the complete description of a physical system. On the basis of the classic model of sea clutter electromagnetic scattering, that is, the random walk model, the characteristics of the sea clutter are systematically analysed using the theory of stochastic differential physical in this study. The sea clutter time-varying electromagnetic scattering characteristics of stochastic differential equations are described. The sea clutter scattering signal amplitude and phase of the diffusion process model are obtained using the generated formulas. The multifractal spectrum is analysed. The accuracy of this model is verified by the Canadian McMaster IPIX radar data set. The experimental results show that the stochastic differential model based on the random walk is an effective algorithm to describe the mechanism of the sea clutter.
  • The National Natural Science Foundation of China under contract No. 61401493; the Naval University of Engineering Natural Science Foundation of China under contract No. HGDQNSQJJ15003.
Feng Tao, Field T R, Haykin S. 2007. Stochastic differential equation theory applied to wireless channels. IEEE Transactions on Communications, 55(8): 1478–1483
Granström K, Natale A, Braca P, et al. 2015. Gamma Gaussian inverse wishart probability hypothesis density for extended target tracking using X-band marine radar data. IEEE Transactions on Geoscience and Remote Sensing, 53(12): 6617–6631
Guan Jian, Liu Ningbo, Zhang Jian, et al. 2010. Multifractal correlation characteristic of real sea clutter and low-observable targets detection. Journal of Electronics & Information Technology (in Chinese), 32(1): 54–61
Güntürkün U. 2015. Bivariate empirical mode decomposition for cognitive radar scene analysis. IEEE Signal Processing Letters, 22(5): 603–607
Jakeman E, Pusey P N. 1976. A model for non-Rayleigh sea echo. IEEE Transactions on Antennas and Propagation, 24(6): 806–814
Jakeman E, Tough R J A. 1987. Generalized K distribution: a statistical model for weak scattering. Journal of the Optical Society of America A, 4(9): 1764–1772
Li Yang, Lv Xiaowen, Liu Kuisheng, et al. 2014. Fractal-based weak target detection within sea clutter. Acta Oceanologica Sinica, 33(9): 68–72
Liu Ningbo, Huang Yong, Guan Jian, et al. 2012. Fractal analysis of real sea clutter in frequency domain. Journal of Electronics & Information Technology (in Chinese), 34(4): 929–935
Rosenberg L, Bocquet S. 2015. Application of the Pareto plus noise distribution to medium grazing angle sea-clutter. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 8(1): 255–261
Suresh P, Thayaparan T, Venkataramaniah K. 2015. Fourier-bessel transform and time-frequency-based approach for detecting manoeuvring air target in sea-clutter. IET Radar, Sonar & Navigation, 9(5): 481–491
Velotto D, Soccorsi M, Lehner S. 2014. Azimuth ambiguities removal for ship detection using full polarimetric X-band SAR data. IEEE Transactions on Geoscience and Remote Sensing, 52(1): 76–88
Wu Gengkun, Ji Guangrong, Ji Tingting, et al. 2014. Study of electromagnetic scattering from two-dimensional rough sea surface based on improved Wen’s spectrum. Acta Physica Sinica (in Chinese), 63(13): 134203
Xing Hongyan, Zhu Qingqing, Xu Wei. 2014. A method of weak target detection based on the sea clutter. Acta Physica Sinica (in Chinese), 63(10): 100505
Xiong Gang, Zhang Shuning, Zhao Huichang. 2014. Time-varying singularity spectrum distribution of sea clutter based on wavelet leaders. Acta Physica Sinica (in Chinese), 63(15): 150503
Xu Xinyu, Zhang Yushi, Li Xin, et al. 2014. KK distribution modeling with L band low grazing sea clutter. Systems Engineering and Electronics (in Chinese), 36(7): 1304–1308
Zhang Yushi, Xu Xinyu, Yin Yalei, et al. 2014. Research on amplitude statistics of l-band low grazing angle sea clutter. Journal of Electronics & Information Technology (in Chinese), 36(5): 1044–1048
Year 2017 volume 36 Issue 9
PDF
50
27
Cite this Article
BibTeX
Article Info
doi: 10.1007/s13131-017-1107-y
  • Receive Date:2016-08-07
  • Online Date:2026-04-16
  • Published:2017-09-01
Article Data
Affiliations
History
  • Received:2016-08-07
  • Accepted:2016-11-07
Funding
The National Natural Science Foundation of China under contract No. 61401493; the Naval University of Engineering Natural Science Foundation of China under contract No. HGDQNSQJJ15003.
Affiliations
    1 Electronics Engineering College, Naval University of Engineering, Wuhan 430033, China

Corresponding:

References
Share
https://castjournals.cast.org.cn/joweb/aos/EN/10.1007/s13131-017-1107-y
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数

Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)

Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT