A Novel Parameter Estimation Method Based on a Tuneable Sigmoid in Alpha-Stable Distribution Noise Environments

In this paper, a novel method, that employs a fractional Fourier transform and a tuneable Sigmoid transform, is proposed, in order to estimate the Doppler stretch and time delay of wideband echoes for a linear frequency modulation (LFM) pulse radar in an alpha-stable distribution noise environment. Two novel functions, a tuneable Sigmoid fractional correlation function (TS-FC) and a tuneable Sigmoid fractional power spectrum density (TS-FPSD), are presented in this paper. The novel algorithm based on the TS-FPSD is then proposed to estimate the Doppler stretch and the time delay. Then, the derivation of unbiasedness and consistency is presented. Furthermore, the boundness of the TS-FPSD to the symmetric alpha stable (SαS) noise, the parameter selection of the TS-FPSD, and the feasibility analysis of the TS-FPSD, are presented to evaluate the performance of the proposed method. In addition, the Cramér–Rao bound for parameter estimation is derived and computed in closed form, which shows that better performance has been achieved. Simulation results and theoretical analysis are presented, to demonstrate the applicability of the forgoing method. It is shown that the proposed method can not only effectively suppress impulsive noise interference, but it also does not need a priori knowledge of the noise with higher estimation accuracy in alpha-stable distribution noise environments.


Introduction
The joint estimation of the Doppler and time delay of a noise-contaminated signal is a fundamental problem in radar and sonar systems and this has been extensively addressed for the case involving narrowband signals [1][2][3][4][5][6]. Weiss [7], Remley [8], and Qu [9] indicated that the narrowband model and the corresponding narrowband signal processing techniques are applicable when BT c/2v, where B is the bandwidth of the transmitted signal, T is the duration of the transmitted signal, v is the relative velocity between the target and the sensor, and c is the propagation speed of the signal. The wideband signals, such as linear frequency modulation (LFM) signals with a large time frequency-bandwidth product, are frequently used in sonar and radar systems because of their lower probability of interception. In many modern radar systems, however, a wideband signal is utilized, and the narrowband model is not appropriate. In applications where BT c/2v is invalid, the wideband model has to be employed. In wideband radar systems, the echo often contains a Doppler stretch (DS), and not a Doppler shift only, which results in difficulty in the parameters' estimation. For the determination of the range and the relative velocity of the target, the accurate estimation of these parameters is crucial.
Since the fractional Fourier transform (FRFT) spectrum of a LFM signal has a greater energy concentration characteristic, the FRFT, as a new time-frequency tool, has attracted more attention, and has been widely applied to the parameter estimation of LFM signal [10][11][12][13][14][15][16]. Zhao et al. propose a method for estimating the LFM signal by utilizing the pulse compression in both the time domain and the FRFT domain [10]. The parameters of the LFM signal can be estimated by a fractional Fourier transform [11,12], a fractional correlation [13], and a fractional power spectrum [14][15][16].
Up to now, in most parameter estimation methods for array signal processing, the additive noise is assumed to be Gaussian. Studies and experimental measurements have shown that a broad class of noise such as underwater acoustic noise, atmospheric noise, multiuser interference, and radar clutters in real world applications are non-Gaussian, primarily owing to impulsive phenomena [17,18]. Taking these scenarios into account, it is inappropriate to model the noise as Gaussian noise. Researchers have studied this impulsive nature, and shown that the symmetric alpha stable (SαS) processes are better models for impulsive noise than the Gaussian processes. The conventional algorithms based on second-order statistics degenerate severely in the SαS noise environment.
To reduce the SαS noise interference, many parameter estimation algorithms based on the fractional lower-order statistics (FLOS) have been proposed [16,[19][20][21]. However, these algorithms do need a priori knowledge of the SαS noise and have other limitations, where the characteristic exponent α and the fractional lower order of moments p must meet the need of 1 ≤ p < α or 0 < p < α/2, otherwise those algorithms performance can degrade seriously, or even become invalid, while the fractional lower-order moment value is not appropriate.
The Sigmoid function is widely used as a common nonlinear transform. The Sigmoid function can suppress impulsive noise interference, and this does not depend on a priori knowledge of the noise. Yu et al., propose a method based on generalized Sigmoid cyclic cross-ambiguity function to estimate the time delay and Doppler frequency shift in the impulsive noise and co-channel interference [22]. In [22], the signal is assumed to be the real signal. However, when the signal is a complex signal, the results do not always hold.
To handle this problem, a novel concept termed the tuneable Sigmoid transform fractional correlation (TS-FC) is proposed in this paper, and a relative method named the tuneable Sigmoid fractional power spectrum density (TS-FPSD) is presented, to fulfill the needs mentioned above. The novel algorithm based on the TS-FPSD is then proposed to estimate the Doppler stretch and time delay. In addition, we address unbiasedness and consistency by adding their corresponding derivation. Furthermore, the boundness of the TS-FPSD to the SαS noise, the parameter selection of the TS-FPSD, the feasibility analysis of the TS-FPSD, and the Cramér-Rao bound for parameter estimation are presented, to evaluate the performance of the proposed method. The proposed method does not need a priori knowledge of the alpha stable distribution noise. This paper is organized as follows. Section 2 presents a signal model of wideband echoes in an alpha-stable distribution noise environment. In Section 3, the analysis of the fractional power spectrum density is presented. In Section 4, a novel tuneable Sigmoid fractional correlation function (TS-FC) and a novel tuneable Sigmoid fractional power spectrum density function (TS-FPSD) are defined. In addition, unbiasedness and consistency are derived. In Section 5, a novel Doppler stretch and time delay estimation method based on TS-FPSD for the SαS noise is proposed. In addition, the boundness of the TS-FPSD to the SαS noise, parameter selection of the TS-FPSD, and a feasibility analysis of the TS-FPSD are analyzed, and the Cramér-Rao bound for parameter estimation is derived. In Section 6, the performance of the parameter estimation algorithm is studied through extensive numerical simulations. Finally, conclusions are drawn in Section 7.

The Signal Model
We consider x(t) as a LFM signal, defined as: A exp j2π f 0 t + 1 2 µ 0 t 2 , 0 < t < T 0, where f 0 and µ 0 are the initial frequency and the frequency modulation rate, respectively. A denotes the signal amplitude, T 0 denotes the modulation period, and T denotes the LFM pulse duration. In a wideband radar system, the echo from a wideband signal often contains a Doppler stretch (DS), and not a Doppler shift only. The echo y(t) can be expressed as [9]: where β l denotes the attenuation factor of the lth multipath,σ l is the Doppler stretch, τ l is the time delay (TD), and L denotes the number of the echo. The noise n(t) is a sequence of the independent and identically distributed (i.i.d) isotropic complex SαS random variable.

The SαS Distribution Noise Model
In order to evaluate the robustness of the proposed method, the α-stable distribution noise is modeled by a complex isotropic symmetric α-stable (SαS) noise distribution [16,17]. The characteristic function of the SαS distribution is defined as follows: where α(0 < α ≤ 2) is usually called the characteristic exponent. When α < 2, the distribution is algebraic-tailed with a tail constant α, implying infinite variance. The smaller it is, the heavier the tails of the density. When α = 2, the SαS distribution reduces to the Gaussian distribution. The parameter γ, usually called the dispersion, is a positive constant that is related to the scale of the distribution. The parameter γ plays a role that is analogous to that of the variance for a second-order process. The following proposition gives us a closed-form expression for the geometric power of the symmetric α-stable random variables [18]: where C ≈ 1.78. Since the α-stable distribution with α < 2 determines an infinite variance, we describe the signal-to-noise condition of SαS using the generalized signal-noise-ratio (GSNR) [17], which is defined as: where σ 2 x and γ are the variance of the underlying signal and dispersion of the SαS noise, respectively.

The Fractional Correlation Function and the Fractional Power Spectrum Density
The fractional Fourier transform (FRFT) is a generalization of the FT, and it can be interpreted as a rotation of the signal to any angles in the time-frequency plane [15,16]. The fractional correlation function (FC) of the signal x(t) is defined as [15]: where R xx (t + ξ, t) is the correlation function of the signal x(t), ξ denotes the delay, and ρ denotes the rotation angle in the FRFT domain. The fractional power spectrum density function (FPSD) of the signal x(t) is expressed as [15]: where A ρ = 1−j cot ρ 2π , m denotes the frequency in FRFT domain, and F ρ [·](m) denotes the fractional Fourier transform.
The continuous FRFT [11] of a signal x(t) with a rotation angle ρ, is defined as: where F b denotes the FRFT operator, b(0 < b ≤ 2) denotes the fractional order, ρ ≡ bπ/2, and K b (t, m) is the kernel function of the fractional Fourier transform. K b (t, m) can be expressed as:

The Parameter Estimation Based on the Fractional Power Spectrum Density Function
From Equations (2) and (6), the fractional correlationR ρ yy (ξ) of the echo signal y(t) can be expressed as (see Appendix A for details): whereR ρ yy,lq (ξ) denotes the fractional correlation function between y l (t) and y q (t),R ρ yn (ξ) denotes the fractional correlation functions between the echo signal y(t) and the noise n(t), andR ρ yn (ξ) is treated as a random interference.R ρ nn (ξ) denotes the fractional correlation function of the noise n(t). When 2πµ 0 σ 2 l = − cot ρ l , Equation (10) can be rewritten as: From Equation (11), we can find thatR ρ l yy (ξ) has the expression of a LFM signal, which has the characteristics of energy that is concentrated in the FRFT domain. However,R ρ l yy,lq (ξ),R ρ l yn (ξ), and R ρ l nn (ξ) do not have such characteristics in the FRFT domain; the amplitudes of TS-FPSD of different targets are very low at (ρ l , m l ), and these signals are not considered as randomly interfering. Therefore, the fractional power spectrum density P ρ yy (m) of the echo signal y(t) can be computed as: The P ρ yy (m) forms a pulse in the FRFT domain, and its peak value appears at (ρ l , m l ) as: Then, it follows directly from Equation (13) that the Doppler stretch σ l and the time delay τ l are estimated by [16]: When the alpha-stable distribution noise is added, this peak location algorithm may fail. The reason is that the alpha-stable distribution does not have a finite α-order moment and other higher-than-α-order moments, and the fractional correlation function and fractional power spectrum density function are based on a second-order moment. Accordingly, the fractional power spectrum density algorithm will become unbounded when the received signal contains an alpha-stable distribution noise. Therefore, we present a tuneable Sigmoid transform, to suppress the alpha-stable distribution noise interference.

The Tuneable Sigmoid Transform
The Sigmoid function is widely used as a common nonlinear transform [23][24][25]. Its definition is shown in Equation (15): where λ is the inclined coefficient to adjust x(t) at different scales. The parameter λ is used as a scale factor to fit various signals and noises. A proper λ will retain sufficient information regarding the time delay and Doppler stretch, and attenuate most of the outliers at the same time. Thus, the proper selection of λ will ensure the accuracy of the estimation. The selection of λ is analyzed in Section 5.3. Through the analysis, the inclined coefficient λ for the TS-FPSD is set as λ = 1 in all the later simulations of this paper. For a SαS process with a = 0, the Sigmoid transform has some properties as follows [22,26]: is a SαS process with β = 0 and a = 0, then Sigmoid[x(t)] has a symmetric distribution with zero mean in its probability density function.
Property 3. If x(t)is a SαS process with a = 0, then Sigmoid[x(t)] has a finite second-order moment with zero mean (referred to as second order moment process).
According to the first three properties, we can derive the following property as: , then X(t) has the same Doppler shift asx(t).
If we set F x (ω) is the Fourier transform of x(t) and F X (ω) is that of X(t), they satisfy the following The frequency shift between x(t) and x 2 (t) is the same as the frequency shift between X(t) and X 2 (t). The frequency shift could be arbitrary as the parameter c varies.
Therefore, X(t) has the same modulation characteristic as x(t).
, then X(t) has the same time delay as x(t).
Proof. Set x 1 (t) = x(t − D), then: That is to say, that the Sigmoid transform does not change the time delay contained in x(t).
Property 6. If x(t) is a periodic function, then X(t) has the same period as x(t).

Proof.
For a periodic signal x(t), satisfying x(t + T ) = x(t), where T is its period. Then: Accordingly, X(t) is periodic. Since [−x(t)] − 1 is a monotonic increasing function for x(t), the period of X(t) is the same as that of x(t).

Definition of the Tuneable Sigmoid-FC and the Tuneable Sigmoid-FPSD
Overcoming the limitations of the performance degradation of existing methods based on the fractional Fourier transform in the alpha-stable distribution noise-and the fractional lower-order statistics-based methods depends on the priori knowledge of noise. This paper presents two novel function definitions, the TS-FC and the TS-FPSD.
A novel fractional correlation,R TS xx,ρ (ξ), referred to as the TS-FC, is defined as: A novel power spectrum function P TS xx,ρ (m), referred to as the TS-FPSD, is defined as: The TS-FC in the discrete-time case can be defined as follows [15]: where T is the sampling period, t = nT. Let w = mT; the discrete-time TS-FPSD can be defined as: Figure 1 shows the spectrum of the FPSD and the TS-FPSD in the no-noise and impulsive noise environments, respectively. In the α-stable distribution noise environment with GSNR = 5 dB, α = 1.2, and L = 2, Figure 1a,b shows the time-frequency distribution of the fractional correlation (FC) and the FPSD of two LFM signals in the no-noise environment. From Figure 1a,b, we can find that the plots are smooth because there is no interference of noise. Figure 1c,d present that the time-frequency distribution of FC and FPSD of two echoes with the SαS noise. Figure 1e,f present the time-frequency distribution of the TS-FC and the TS-FPSD of two echoes with SαS noise. From Figure 1, it is clearly seen that identifying the correct peak location is not trivial, as the FPSD peak could not be distinguished from the SαS noise. Accordingly, the estimation performance of the FPSD method degraded severely in the SαS noise environment. After the application of the tuneable Sigmoid transformation, the SαS noise was suppressed effectively, and the TS-FPSD spectrum formed an obvious pulse in the FRFT domain. Thus, the method based on the TS-FPSD yielded better estimation performance. The Sigmoid transform can result in some information loss, same as other nonlinear transforms, but its capability for suppression impulsive noise is employed in this paper.
A novel power spectrum function

TS xx
Pm  , referred to as the TS-FPSD, is defined as: The TS-FC in the discrete-time case can be defined as follows [15]: where T is the sampling period, ; the discrete-time TS-FPSD can be defined as: Figure 1 shows the spectrum of the FPSD and the TS-FPSD in the no-noise and impulsive noise environments, respectively. In the α-stable distribution noise environment with GSNR = 5 dB, , and 2 L  , Figure 1a,b shows the time-frequency distribution of the fractional correlation (FC) and the FPSD of two LFM signals in the no-noise environment. From Figure 1a,b, we can find that the plots are smooth because there is no interference of noise. Figure 1c,d present that the timefrequency distribution of FC and FPSD of two echoes with the SS  noise. Figure 1e,f present the time-frequency distribution of the TS-FC and the TS-FPSD of two echoes with SS  noise. From Figure 1, it is clearly seen that identifying the correct peak location is not trivial, as the FPSD peak could not be distinguished from the SS  noise. Accordingly, the estimation performance of the FPSD method degraded severely in the SS  noise environment. After the application of the tuneable Sigmoid transformation, the SS  noise was suppressed effectively, and the TS-FPSD spectrum formed an obvious pulse in the FRFT domain. Thus, the method based on the TS-FPSD yielded better estimation performance. The Sigmoid transform can result in some information loss, same as other nonlinear transforms, but its capability for suppression impulsive noise is employed in this paper.

Unbiasedness and Consistency
In statistical theory, the bias (or bias function) of an estimator is the difference between the estimator's expected value and the true value of the parameter being estimated. In general, the bias is related to consistency, where consistent estimators are convergent and asymptotically unbiased (hence they converge to the correct value as the number of data points grows arbitrarily large) [27,28].
In this section, the unbiasedness and consistency of the TS-FC and TS-FPSD are analyzed, to evaluate the performance of the proposed method. The derivation of the unbiasedness and consistency are presented in details in Appendix B.

Unbiasedness and Consistency of the TS-FC
Let the estimator of the TS-FC According to the definition of the unbiasedness, we can obtain as follows: where:

Unbiasedness and Consistency
In statistical theory, the bias (or bias function) of an estimator is the difference between the estimator's expected value and the true value of the parameter being estimated. In general, the bias is related to consistency, where consistent estimators are convergent and asymptotically unbiased (hence they converge to the correct value as the number of data points grows arbitrarily large) [27,28].
In this section, the unbiasedness and consistency of the TS-FC and TS-FPSD are analyzed, to evaluate the performance of the proposed method. The derivation of the unbiasedness and consistency are presented in details in Appendix B.

Unbiasedness and Consistency of the TS-FC
Let the estimator of the TS-FCR TS xx,ρ (k) is: According to the definition of the unbiasedness, we can obtain as follows: where: For a given |k| value and as M → ∞ , we can obtain E R TS xx,ρ (k) =R TS xx,ρ (k), i.e., bia R TS xx,ρ (k) = 0. Therefore, the estimator of the TS-FC is the asymptotic unbiased estimation.
According to the definition of the consistency, we can obtain: When M → ∞ , we can obtain Var R TS xx,ρ (k) = 0.
In summary, according to the bia R TS xx,ρ (k) = 0 and Var R TS xx,ρ (k) = 0, we may draw a conclusion that the estimator of the TS-FC is the asymptotic consistent estimation.

Unbiasedness and Consistency of the TS-FPSD
Let the estimator of the TS-FPSD is: According to the definition of the unbiasedness, we can obtain: For a given |k| value and as M → ∞ ,E P TS xx,ρ (w) = P TS xx,ρ (w), i.e., bia P TS xx,ρ (w) = 0. Therefore, the estimator of the TS-FPSD is the asymptotic unbiased estimation.
According to the definition of the consistency, we can obtain: For a given |k| and as M → ∞ , we can get: In summary, according to bia P TS xx,ρ (w) = 0 and Var P TS xx,ρ (w) = 0, we may draw a conclusion that the estimator of the TS-FPSD is the asymptotic consistent estimation.

Joint Doppler Stretch and Time Delay Estimation
The echo signal y(t) with α-stable distribution noise can be expressed as: where the noise n(t) denotes the SαS noise.
According to the definition of the TS-FC, the TS-FCR TS yy,ρ (ξ) of the echo signal y(t) can be expressed as: (31) According to the definition of the TS-FPSD, the TS-FPSD P TS yy,ρ (m) of the signal y(t) can be expressed as The joint estimation method for the Doppler stretch σ l and the time delay τ l based on the TS-FPSD, is given by [16]:σ Accordingly, the estimation of the Doppler stretch and the time delay in wideband echoes for a LFM pulse radar under an alpha-stable distribution noise was achieved via the proposed tuneable Sigmoid fractional power spectrum density function. The steps involved in this process are as follows: Step 1 Obtain the echo signal y(t).
Step 4 Search for the peaks of P TS yy,ρ (m) and obtain the locations of these peaks (τ l ,σ l ), for l = 1, . . . , L.
Step 5 Estimate the Doppler stretch and time delay according to Equation (33).

The Boundness of the TS-FPSD to the SαS Noise
We consider z(t) as an observed signal, defined as: where s(t) denotes the signal, and the noise n(t) is a sequence of the i.i.d isotropic complex SαS random variable. According to the definition of the TS-FPSD, we can obtain: According to the properties of the tuneable Sigmoid transform, the SαS process with a = 0 can be transformed to a second-order moment process by the Sigmoid transform. Therefore, , which can guarantees the boundness of P TS zz,ρ (m) under the SαS noise. Furthermore, the transformation does not change the estimation result of the time delay and the Doppler frequency. Therefore, the TS-FPSD method can be used to estimate the parameters of the wideband echoes y(t) under the α-stable distribution noise.

Parameter Selection of the TS-FPSD
The inclined coefficient λ was used as a scale factor to fit various signals and noises. A proper λ will retain sufficient information that is associated with the time delay and Doppler stretch, and it will attenuate most of the impulsive noise at the same time. Thus, the proper selection of λ will ensure the accuracy of the estimation. According to Yu et al., the attenuation upon x(t) in the tuneable Sigmoid function changes as |x(t)| changes [26]. This concept is illustrated in Figure 2. When the signal x(t) is the real signal without noise, no matter what the value of the inclined coefficient λ is, |Sigmoid[x(t)]| < 1 is true, as illustrated in Figure 2a. Furthermore, when the real signal x(t) contains the real impulsive noise, no matter what the value of the inclined coefficient λ is, |Sigmoid[x(t)]| < 1 is also true, as illustrated is Figure 2b.
However, when the signal x(t) is a complex signal, the above results do not always hold. Figure 3 demonstrates the Sigmoid function curves of the complex signal with respect to λ, and the ratio between the real and imaginary component. Figure 3a shows that the tuneable Sigmoid function changes with respect to λ when the ratio between the real and imaginary components is 2. Figure 3b shows that the tuneable Sigmoid function changes with respect to x(t), with the same amplitude, but with a different ratio between the real and imaginary component when the inclined coefficient λ is 1. We found that the suppression capability may be increased by increasing the ratio between the real and imaginary components.
According to the properties of the tuneable Sigmoid transform, the SS  process with 0 a  can be transformed to a second-order moment process by the Sigmoid transform. Therefore, Furthermore, the transformation does not change the estimation result of the time delay and the Doppler frequency. Therefore, the TS-FPSD method can be used to estimate the parameters of the wideband echoes   yt under the  -stable distribution noise.

Parameter Selection of the TS-FPSD
The inclined coefficient  was used as a scale factor to fit various signals and noises. A proper  will retain sufficient information that is associated with the time delay and Doppler stretch, and it will attenuate most of the impulsive noise at the same time. Thus, the proper selection of  will ensure the accuracy of the estimation. According to Yu et al., the attenuation upon   xt in the tuneable Sigmoid function changes as   xt changes [26]. This concept is illustrated in Figure 2.
When the signal   xt is the real signal without noise, no matter what the value of the inclined coefficient  is, is true, as illustrated in Figure 2a. Furthermore, when the real signal   xt contains the real impulsive noise, no matter what the value of the inclined coefficient  is, is also true, as illustrated is Figure 2b.
However, when the signal   xt is a complex signal, the above results do not always hold. Figure 3 demonstrates the Sigmoid function curves of the complex signal with respect to  , and the ratio between the real and imaginary component. Figure 3a shows that the tuneable Sigmoid function changes with respect to  when the ratio between the real and imaginary components is 2. Figure   3b shows that the tuneable Sigmoid function changes with respect to   xt , with the same amplitude, but with a different ratio between the real and imaginary component when the inclined coefficient  is 1. We found that the suppression capability may be increased by increasing the ratio between the real and imaginary components.   We observed that   Sigmoid 1 xt    may not be true when   xt is a complex signal, as illustrated in Figure 3. Furthermore, we also found that the amplitude of the complex signal, the inclined coefficient  , and the ratio between the real and imaginary components had some effects in suppressing the impulsive noise capability of the tuneable Sigmoid function. Figure 4 shows that the suppression ability of the tuneable Sigmoid function for a SS  noise with GSNR 5dB  and 1.2   . In this simulation, the signal   xt is the complex signal with impulsive noise. The Sigmoid function employs the tuneable parameter  , which could be used to control the inclination of the curve. The outliers can be suppressed after the transformation. From Figure 4, we found that the suppression capability may be decreased by increasing the inclined coefficient  , and the tuneable Sigmoid function with 3   fails to suppress the interference of the  -stable distribution noise. A higher inclined coefficient  had negative impacts on the suppression capability of the impulsive noise. Therefore, the proper selection of  affected accuracy of the estimation.  We observed that |Sigmoid[x(t)]| < 1 may not be true when x(t) is a complex signal, as illustrated in Figure 3. Furthermore, we also found that the amplitude of the complex signal, the inclined coefficient λ, and the ratio between the real and imaginary components had some effects in suppressing the impulsive noise capability of the tuneable Sigmoid function. Figure 4 shows that the suppression ability of the tuneable Sigmoid function for a SαS noise with GSNR = 5dB and α = 1.2. In this simulation, the signal x(t) is the complex signal with impulsive noise. The Sigmoid function employs the tuneable parameter λ, which could be used to control the inclination of the curve. The outliers can be suppressed after the transformation. From Figure 4, we found that the suppression capability may be decreased by increasing the inclined coefficient λ, and the tuneable Sigmoid function with λ = 3 fails to suppress the interference of the α-stable distribution noise. A higher inclined coefficient λ had negative impacts on the suppression capability of the impulsive noise. Therefore, the proper selection of λ affected accuracy of the estimation. We observed that   Sigmoid 1 xt    may not be true when   xt is a complex signal, as illustrated in Figure 3. Furthermore, we also found that the amplitude of the complex signal, the inclined coefficient  , and the ratio between the real and imaginary components had some effects in suppressing the impulsive noise capability of the tuneable Sigmoid function. Figure 4 shows that the suppression ability of the tuneable Sigmoid function for a SS  noise with GSNR 5dB  and 1.2   . In this simulation, the signal   xt is the complex signal with impulsive noise. The Sigmoid function employs the tuneable parameter  , which could be used to control the inclination of the curve. The outliers can be suppressed after the transformation. From Figure 4, we found that the suppression capability may be decreased by increasing the inclined coefficient  , and the tuneable Sigmoid function with 3   fails to suppress the interference of the  -stable distribution noise. A higher inclined coefficient  had negative impacts on the suppression capability of the impulsive noise. Therefore, the proper selection of  affected accuracy of the estimation.

Feasibility Analysis of the TS-FPSD
According to Property 4 of the Sigmoid transform, the Sigmoid transform did not change the modulation characteristic of the signal, i.e., Sigmoid[x(t)] and x(t) had the same modulation characteristics. The simulation results are illustrated in Figure 5 below, to verify this property.

Feasibility Analysis of the TS-FPSD
According to Property 4 of the Sigmoid transform, the Sigmoid transform did not change the modulation characteristic of the signal, i.e.,     Sigmoid x t and   x t had the same modulation characteristics. The simulation results are illustrated in Figure 5 below, to verify this property. In summary, the parameters of the Doppler stretch and the time delay could be estimated by searching for the peak of the TS-FPSD. In summary, the parameters of the Doppler stretch and the time delay could be estimated by searching for the peak of the TS-FPSD.

Feasibility Analysis of the TS-FPSD
According to Property 4 of the Sigmoid transform, the Sigmoid transform did not change the modulation characteristic of the signal, i.e.,     Sigmoid x t and   x t had the same modulation characteristics. The simulation results are illustrated in Figure 5 below, to verify this property.

The Cramer-Rao Bound
In this section, we derived a novel explicit expression for the exact Cramer-Rao Bound (CRB) on the accuracy of estimating the signal model parameters.
The CRB expresses a lower bound for the variance of an unbiased estimate and is, in general, not too difficult to compute. By comparing the performance of an estimator to the CRB, we can often

The Cramer-Rao Bound
In this section, we derived a novel explicit expression for the exact Cramer-Rao Bound (CRB) on the accuracy of estimating the signal model parameters.
The CRB expresses a lower bound for the variance of an unbiased estimate and is, in general, not too difficult to compute. By comparing the performance of an estimator to the CRB, we can often have an indication on how close the estimator is to the optimum.
The echoes signal can be expressed as the following: First, we obtained closed-form expressions for all particular sub-blocks of the Fisher information matrix (FIM). The element i, j of the FIM for estimating the vector ξ = [σ, τ] T can be shown as [29][30][31]: We assumed that the noise was a sequence of the i.i.d isotropic complex SαS random variable. The geometric power S 0 is used to represent the power of symmetric α-stable random noise, i.e., Q n = S 0 I N . Since x(σ, τ, t)β and Q n depend on different elements of ξ, it is clear that FIM will be block diagonal with respect to the signal (ξ = [σ, τ] T ) and noise parameters. In particular, the first term of Equation (38) will give a nonzero result only for the noise block. Since we are concerned only with the CRB for the signal parameters, we need only consider the second term: Using Equations (38) and (39), the following explicit expressions for the blocks of the FIM are derived as follows: where: The expression for the CRB, shown in Equation (48), is obtained by substituting Equations (40)-(47) into Equation (39):

Simulation Results
In this section, we performed four types of simulation experiments to evaluate the relative performances of the FPSD [12], the FLOS-FPSD [16], and the TS-FPSD methods under the α-stable distribution noise, respectively.
The parameters of the transmitted LFM signal in the simulation are assumed as follows. The initial frequency f 0 = 0.2 f s and the modulation rate is set to µ 0 = 0.1 f 2 s /N. The sampling rate is set to f s = 1 MHz with a sampling length of N = 1000, T = 1 ms. The number of multipath is L = 2, and the Doppler stretch and time delay are set to σ 1 = 0.9,σ 2 = 1.1,τ 1 = 20/ f s and τ 2 = 60/ f s , respectively. The Root Mean Square Error (RMSE) is defined as: wherex 1 andx 2 are the estimation of x 1 and x 2 , and K is the Monte Carlo number. The RMSE of the time delay and Doppler stretch can be expressed as: and: The numbers of Monte Carlo runs was set to 200 in Simulations 2 and 3.

Simulation 1: Estimation Accuracy with Respect to λ
To evaluate the performance of the TD and DS with respect to λ in this simulation, the characteristic exponent α was set to α = 1.2 and GSNR = 5 dB. The RMSE was used to evaluate the performance of the TS-FPSD with respect to λ, as illustrated is Figure 7. The numbers of Monte Carlo runs was set to 200 in Simulations 2 and 3.

Simulation 1: Estimation Accuracy with Respect to 
To evaluate the performance of the TD and DS with respect to  in this simulation, the characteristic exponent  was set to and GSNR=5dB . The RMSE was used to evaluate the performance of the TS-FPSD with respect to  , as illustrated is Figure 7. From Figure 7, we observed that the performance estimation of the TD and DS using provided a better performance than that using other  values, for the case of the alpha-stable noise. Therefore, the inclined coefficient  for the TS-FPSD was set as 1   in all the later simulations of this paper.  Figure 7, we observed that the performance estimation of the TD and DS using λ ∈ [0.8, 2.7] provided a better performance than that using other λ values, for the case of the alpha-stable noise. Therefore, the inclined coefficient λ for the TS-FPSD was set as λ = 1 in all the later simulations of this paper.
Simulation 2: FPSD, FLOS-FPSD, and TS-FPSD for a Single Estimation   Figures 8 and 9 show the estimation results of the FPSD, FLOS-FPSD, and TS-FPSD for a single trial of data under the SαS noise with GSNR = 5 dB and α = 1.3 and α = 1.1. In order to show the peak location information and the performance of the algorithm more clearly, the 2D rotation angle plane and the 2D frequency plane are shown. The 2D rotation angle plane and 2D frequency plane could better demonstrate the peak location. In Figures 8 and 9, the red line denotes the true values of the rotation angle and frequency in the FRFT domain. Futhermore, Table 1 shows the comparison of three algorithms for impulsive noise suppression. The FLOS-FPSD failed to obtain the correct spectrum peak under the  -stable distribution noise GSNR 5dB  with 1.1   and 1.0 p  ; however, the TS-FPSD peak could be easily separated from the impulsive noise, as illustrated in Figure 9, mainly due to the fact that the Sigmoid transform could suppress the impulsive noise better than the FLOS. Therefore, the performance of the TS-FPSD outweighed those of the FLOS-FPSD. The proposed method based on the TS-FFPSD could effectively suppress impulsive noise interference, yielded an accurate peak estimation, and had a better estimation performance.     From Figure 8a, we found that the FPSD algorithm failed when the SαS noise occurred. The reason is that the FPSD method did not have the ability to suppress impulsive noise. Since the second-order moment of a SαS random variable with 0 < α < 2 does not exist, and the fractional correlation function was based on second-order moments, the performance of the FPSD degraded severely. The FLOS-FPSD algorithm, combining the fractional lower order statistics theory with the fractional power spectrum density function, effectively suppressed the α-stable distribution noise interference, so the FLOS-FPSD could obtain a clear peak under the SαS noise of GSNR = 5 dB, with α = 1.3 and p = 1.0. However, the FLOS-FPSD failed to obtain the correct spectrum peak under the SαS noise with α = 1.3 and p = 1.4, mainly due to the fact the fractional lower-order moment p value was not appropriate, as illustrated in Figure 8b,c. In fractional lower order statistics theory, the characteristic exponent of the noise must be estimated to ensure that 1 ≤ p < α or 0 < p < α/2, otherwise the algorithm performance can degrade seriously, and may even become invalid, while the fractional lower-order moment value is not appropriate.
The FLOS-FPSD failed to obtain the correct spectrum peak under the α-stable distribution noise GSNR = 5 dB with α = 1.1 and p = 1.0; however, the TS-FPSD peak could be easily separated from the impulsive noise, as illustrated in Figure 9, mainly due to the fact that the Sigmoid transform could suppress the impulsive noise better than the FLOS. Therefore, the performance of the TS-FPSD outweighed those of the FLOS-FPSD. The proposed method based on the TS-FFPSD could effectively suppress impulsive noise interference, yielded an accurate peak estimation, and had a better estimation performance.

Simulation 3: Estimation Accuracy with Respect to GSNR
To evaluate the performance of the time delay (TD) and the Doppler stretch (DS) in this simulation, the characteristic exponent α was set to α = 1.2 and the fractional lower order moment p was set to p = 1.1 and p = 1.4 for the FLOS-FPSD method, respectively. The resulting RMSE performance versus GSNR is illustrated in Figure 10.   Figure 10. From Figure 10, we can find that the FPSD method had a poor estimation performance with the SS  noise interference. On the other hand, combining the fractional lower-order statistics theory with the fractional power spectrum density, the FLOS-FPSD method with  From Figure 10, we can find that the FPSD method had a poor estimation performance with the SαS noise interference. On the other hand, combining the fractional lower-order statistics theory with the fractional power spectrum density, the FLOS-FPSD method with α = 1.2 and p = 1.1 could effectively suppress the SαS noise interference. Accordingly, the FLOS-FPSD method yielded a clear peak under the SαS noise. However, the performance was affected by the fractional lower-order moment p value. The FLOS-FPSD method, with α = 1.2 and p = 1.4 could not accurately estimate the parameters, because the fractional lower-order moment value was not appropriate. On the contrary, the TS-FPSD method could not suppress the SαS noise interference, employing the tuneable Sigmoid transform, but the estimation performance of the TS-FPSD also could not be affected by the fractional lower-order moment p value. Therefore, the performance of the TS-FPSD method outweighed those of the FLOS-based method.

Simulation 4: Estimation Accuracy with Respect to the Characteristic Exponent α
In this simulation, the GSNR was set to 5 dB, and the fractional lower-order moment p was set to p = 1.1 and p = 1.4 for the FLOS-FPSD method, respectively. Figure 11 shows the performance versus the characteristic exponent α. From Figure 11, we found that the FPSD algorithm had a better estimation performance when the characteristic exponent α was close to 2. The FLOS-FPSD method may suppress the α-stable distribution noise interference by employing the fractional lower-order statistic theory. The performance of the FLOS-FPSD method was shown to be better than that of the FPSD method. In this simulation, the GSNR was set to 5 dB, and the fractional lower-order moment p was set to 1.1 p  and 1.4 p  for the FLOS-FPSD method, respectively. Figure 11 shows the performance versus the characteristic exponent  . From Figure 11, we found that the FPSD algorithm had a better estimation performance when the characteristic exponent  was close to 2. The FLOS-FPSD method may suppress the  -stable distribution noise interference by employing the fractional lower-order statistic theory. The performance of the FLOS-FPSD method was shown to be better than that of the FPSD method. Since the FLOS and the Sigmoid transform methods could both suppress the impulsive noise, the suppression capacity of the FLOS method was not sufficient, and the Sigmoid function suppressed the outliers much harder than did the FLOS. Therefore, the estimation performance of the TS-FPSD algorithm was superior to that of FLOS-FPSD algorithm.

Conclusions
In this paper, two novel concepts, the tuneable Sigmoid transform fractional correlation function and the tuneable Sigmoid transform fractional power spectrum density function, are proposed to estimate the time delay and Doppler stretch of the wideband echoes of a LFM pulse radar signal under the presence of an impulsive noise environment. Then, the unbiased estimation and consistent estimation of the algorithm are derived. Furthermore, the boundness of the TS-FPSD to the SS  noise, the parameter selection of the TS-FPSD, the feasibility analysis of the TS-FPSD, and the Cramér-Rao bound for parameter estimation are presented, to evaluate the performance of the proposed method. The simulation results and theoretical analysis are presented to illustrate the validity of the foregoing method. It is clearly shown that the proposed method cannot only can effectively restrain impulsive noise interference, but it also does not depend on a priori knowledge of the noise. In addition, it yields a higher estimation accuracy and a lower computational complexity in the impulsive noise environment. Since the FLOS and the Sigmoid transform methods could both suppress the impulsive noise, the suppression capacity of the FLOS method was not sufficient, and the Sigmoid function suppressed the outliers much harder than did the FLOS. Therefore, the estimation performance of the TS-FPSD algorithm was superior to that of FLOS-FPSD algorithm.

Conclusions
In this paper, two novel concepts, the tuneable Sigmoid transform fractional correlation function and the tuneable Sigmoid transform fractional power spectrum density function, are proposed to estimate the time delay and Doppler stretch of the wideband echoes of a LFM pulse radar signal under the presence of an impulsive noise environment. Then, the unbiased estimation and consistent estimation of the algorithm are derived. Furthermore, the boundness of the TS-FPSD to the SαS noise, the parameter selection of the TS-FPSD, the feasibility analysis of the TS-FPSD, and the Cramér-Rao bound for parameter estimation are presented, to evaluate the performance of the proposed method. The simulation results and theoretical analysis are presented to illustrate the validity of the foregoing ρ yy (ξ) of the echo signal y(t) can be expressed as: whereR ρ yn (ξ) denotes the fractional correlation functions between the echo signal y(t) and the noise n(t), andR ρ yn (ξ) is treated as a random interference.R ρ nn (ξ) denotes the fractional correlation functions of the noise n(t).
Substituting Equations (A2)-(A5) into Equation (A1), we can rewrite Equation (A1) as follows: However,R ρ l yy,lq (ξ),R ρ l yn (ξ) andR ρ l nn (ξ) do not have the characteristics of energy that is concentrated in the FRFT domain, the amplitudes of the TS-FPSD of different targets are very low at (ρ l , m l ), and these signals are not considered as random interference. So we can obtain: In this section, the unbiasedness and consistency are derived. The definitions of the unbiasedness and consistency can be shown as: If the following holds: bia θ = E θ − θ = 0 (A13) Then,θ is an unbiased estimator of the parameter θ. Otherwise,θ is a biased estimator of θ. and if the following holds: Then,θ is a consistent estimator of the parameter θ. Otherwise,θ is not a consistent estimator of the parameter θ. According to the definition of the unbiasedness, we can obtain as follows: where: When |k| is given and M → ∞ , we can obtain E R TS xx,ρ (k) =R TS xx,ρ (k), i.e., bia R TS xx,ρ (k) = 0. Therefore the estimator of the TS-FC is the asymptotic unbiased estimation.
According to the definition of the consistency, we can obtain as follows: where:   When M → ∞ , we can obtain Var R TS xx,ρ (k) = 0.
In summary, according to the bia R TS xx,ρ (k) = 0 and Var R TS xx,ρ (k) = 0, we may draw the conclusion that the estimator of the TS-FC is the the asymptotic consistent estimation.