Multiscale KF Algorithm for Strong Fractional Noise Interference Suppression in Discrete-Time UWB Systems

In order to suppress the interference of the strong fractional noise signal in discrete-time ultrawideband UWB systems, this paper presents a new UWB multi-scale Kalman filter KF algorithm for the interference suppression. This approach solves the problem of the narrowband interference NBI as nonstationary fractional signal in UWB communication, which does not need to estimate any channel parameter. In this paper, the received sampled signal is transformed throughmultiscale wavelet to obtain a state transition equation and an observation equation based on the stationarity theory of wavelet coefficients in time domain. Then through the Kalman filter method, fractional signal of arbitrary scale is easily figured out. Finally, fractional noise interference is subtracted from the received signal. Performance analysis and computer simulations reveal that this algorithm is effective to reduce the strong fractional noise when the sampling rate is low.


Introduction
Ultrawideband UWB radio is an emerging technology for future short-rang high-speed wireless communications.In recent years, the Federal Communication Commission FCC confined the working bandwidth of UWB systems in the 3.1-10.6GHz frequency band rang and the power spectral density PSD in less than −41.3 dBm/MHz 1, 2 .UWB signal is defined as any wireless signal that has a relative bandwidth more than 25% or an absolute bandwidth in excess of 500 MHz, satisfying PSD required by FCC.Compared with the existing radio communication technologies, the spread spectrum technique of UWB systems bears many advantages, for instance, low power consumption, low probability of intercept, strong ability of penetrating materials, and large space capacity.
Since a conventional UWB system transfers data by using ultrashort pulses, it possesses considerably large bandwidths.Usually, it causes mutual interference with the most existing electronic narrowband systems over the same frequency band inevitably, such as radio navigation or detection system, wireless local area network and global positioning system GPS 3 .Therefore, to guarantee the successful deployment of UWB, it is related to not only the deployment of efficient multiple access technologies but also narrow-band interference suppression technologies.
Recently a number of the narrowband interference NBI suppression technologies have been published thoroughly in many references.Some NBI cancellation solutions 4-11 have been utilized in UWB systems, for example, pulse waveform designing technology, minimum mean square error MMSE rake receiver technology, template waveform estimation technology, frequency domain technology, notch filter techniques, and so forth.These methods all selected simple and stable model as NBI model 12, 13 and used a large number of estimated channel parameters to eliminate interference.Moreover, they studied less about the nonstationary narrowband interference model and the reduction of channel parameter estimation.Hence, this paper proposes a new multi-scale Kalman filter algorithm for strong fractional noise interference suppression in time hopping TH UWB systems.As the nonstationary fractional signal is very strong and UWB signal mixed with white noise is very weak, the received signal in this paper is modeled as narrowband signal which concentrated mainly fractional signal energy.Then it performs the multi-scale wavelet analysis of the received signal and establishes the state space model of Kalman filter estimation 14-17 .The state space model includes a time domain state equation and observation equation.They are obtained by the stationarity of fractional wavelet coefficients.Lastly, we can estimate the fractional Brownian motion FBM , and a novel multiscale Kalman filter for UWB systems is proposed.Its performance analysis and computer simulations are also discussed in Section 5.
This paper is organized as follows.In Section 2, the discrete-time model of the received signal is defined.In Section 3, we introduce the FBM.A novel multi-scale Kalman filter algorithm is proposed in Section 4. The performance of this multi-scale Kalman filter is analyzed, and simulation results are given in Section 5.The conclusions are drawn in Section 6.

Discrete-Time UWB System Model
Consider a multiple-access TH-UWB system with single user.The transmitted UWB signal employing binary pulse amplitude modulation 2-PAM is provided by where ω t is a single link transmitted pulse waveform.T f is the pulse repetition period.T c is the TH unit time.N s is the number of repeated pulses over one symbol period.b j/N s ∈ { 1, −1} represents the user's information during the jth frame.If a symbol period duration is T b N s T f , the transmitted rate is R s 1/ N s T f .{C j } is the pseudorandom TH sequence which adds an additional time shift T c to each pulse, and 0 Assume d n stands for the nth symbol of single user.Then, 2.1 can be expressed as However, in practical situation, the received signal must be discrete.Suppose r t is sampled at the rate of 1/T s and an output r n r t | t nT s is received in discrete-time form.Then, the following received signal becomes where M is the number of sampling points in every symbol period.
Thus, in the nth symbol period, 2.4 can be changed into where S n X • η and X v 0 , v 1 , . . ., v M − 1 T .Here η αb n is a constant.In this paper, the estimated scheme is proposed to stress against the fractional narrowband interferer.The whole basic idea is to realize the prediction of FBM signal, so that it can be removed from the received signal to reduce the fractional interference.In theory, this paper preforms an analysis of FBM.However, in practice, the estimation of the pure FBM is difficult to achieve because the FBM is often unknown.And as the narrowband signal occupies most of the received signal energy, the prediction of r t is used to replace the estimation of FBM.On the other hand, the error signal r t − r t , where r t is an estimation of r t , is very important to eliminate the FBM.Kalman filter technique is also cited to estimate the FBM in this paper.It is shown in Figure 1.

Fractional Brownian Motion
Brownian motion is a kind of random, irregular motion.In the 1960s, a new theory of fractional Brownian motion, commonly referred to as the FBM, is advanced by Mandelbrot and Van Ness 18 .
Denote where Γ • is Gamma function and H 0 < H < 1 is a Hurst index.For the general Brownian motion, when H 1/2, it is a standard Brownian motion.Obviously, FBM, in a continuous Gaussian process with a zero mean, has the nonstationarity, self-similarity, and so forth 18, 19 .
Certainly, there are many methods of generating Brownian motion.For example, Weierstrass-Mandelbrot random function method, under certain conditions, emerges FBM through an almost everywhere nondifferentiable continuous function; in fractional Gaussian noise method, it is that the fractional Gaussian noise is represented as the incremental derivative of FBM; in covariance matrix transform method, Gaussian white noise is turned into a sequence with the covariance matrix through the Cholesky decomposition technology to get FBM; wavelet transform method, that is to say, in terms of adjacent times frequency conversion rate, makes Gaussian white noise into FBM by the inverse wavelet transform.However, this paper uses Weierstrass-Mandelbrot random function to obtain FBM.

Multiscale Kalman Filter Algorithm on Fractional Noise Model
With the purpose of NBI suppression to be successfully realized, we consider the discrete orthogonal wavelet decomposition of B H t .Then, 3.1 can be expressed in the following form: is a scale function of the wavelet transform, and the wavelet function is written by Ψ J,n t 2 −j/2 Ψ 2 −j t − n 20, 21 .If 2.4 is converted applying the wavelet decomposition theory, we have the following equation: where j 1, 2, . . ., J. u j n S j n N j n has the mean μ j S j n and variance Var u j n σ 2 u j .{d j n } is stable 22 .Then, due to the stationarity of wavelet coefficients {d j n }, this AR p model can be found in this study: where ε j m with a zero mean and variance σ 2 ε j , independent of S m and B H m , is AWGN.a i , i 1, 2, . . ., p are parameter coefficients.Define

4.4
For 4.3 , we have Gε m 1 .4.5 Equation 4.5 is one-step prediction from d M/2 j j,m to d M/2 j j,m 1 , so there M/2 j iterative steps of d j,n from d j,n 1 , namely, where

4.7
What is more, to facilitate the study below, we need to transform formula 2.5 into an observation equation.That is, r j,n d j,n S j,n N j,n . 4.8 Thus, if 4.6 is seen as a state transition equation in a dynamic space, combined with 4.8 , we can obtain a state space model for B H t : r j,n d j,n u j,n , 4.9 where u j,n S j,n N j,n .

4.10
To estimate the state space model 4.9 in the third part, the new Kalman filter algorithm is structured as follows. 1 Initialization: at n 0, initialize d j,0 and P 0 are given.

4.11
Therefore, this algorithm can predict all the high-frequency coefficients of the wavelet on any scale according to the KF iterative algorithm.It means that FBM can be calculated through the inverse wavelet transform to receive the interference suppression in UWB systems.

Simulation
In this section, the performance of our proposed NBI suppression method is surveyed.It is simulated by the Matlab7.1 software, and the pulse waveform ω t from the TH-UWB system is adopted as follows: where τ denotes the pulse forming factor.Suppose the parameters are T f 8 ns, N s 2, N h 5, M 160, and J 3, and the sampled time interval is 0.05 ns, which is low.Moreover, the Hausdorff-Besicovitch fractional dimension of the graph of Weierstrass-Mandelbrot random function is D 1.5.Therefore, there are some measured results for this simulation in the following portion.In Figure 2, the original fractional noise signal is shown.Figures 3  and 4 describe the estimated fractional signal and the noisy error signal, respectively.That is, the proposed new algorithm can effectively estimate the strong fractional noise from these results of Figures 1-3  Generally, the signal-to-interference ratio SINR can detect the performance of a algorithm very well.In order to preferably evaluate the performance of the new fractional noise suppression algorithm from another angle, SINR is improved in this paper: where e n B H n − B H n , which is the estimated error signal of B H n .σ 2 N is the variance of the white noise.Thus, 100 independent tests are done, where the lower sampling rate is 0.1 ns and D 1.2 is employed.From this simulation, it is indicated that the actual SINR has greatly improved in Figure 5. On the other hand, it is shown that this new method is effective in suppressing the nonstationary B H t for UWB systems, which is a strong NBI.

Conclusion
In this paper, the multi-scale KF solution is suggested to suppress the fractional noise interference in discrete-time UWB systems.Particularly, after a discrete-time received signal is modeled in the receiver and the fractional signal is transformed based on the wavelet analysis, the fractional signal can be estimated by Kalman filter, which is supposed as a strong nonstationary NBI process.From the simulation above, the algorithm estimates the fractional noise signal under the weak UWB signals by combining the theory of the multi-scale wavelet transform and the Kalman filter.So that it achieves the separation of the narrowband signal and wideband signals successfully.One of the advantages of this paper is not estimating any channel parameter.Simulation results show that this method is an effective narrowband interference suppression method of UWB systems.

Figure 1 :
Figure 1: Estimated scheme of fractional noise suppression in UWB systems.

Figure 3 :Figure 4 :
Figure 3: The estimated fractional signal based on multiscale wavelet transform.

Figure 5 :
Figure 5: Comparison of the performance of the ideal and proposed multi-scale KF method.
the number of transmitted symbols.Assume the nonstationary FBM model is used as NBI model.Then, while the transmitted UWB signal propagates through a flat channel and is corrupted by fractional noise signal B H t and additive white Gaussian noise AWGN N t , the observation signal is given by t t .α is the channel gain.N t ∼ N 0, σ 2 N t .S t N t denotes the wideband component of the received signal.