Next Article in Journal
The Need for Systems Awareness to Support Early-Phase Decision-Making—A Study from the Norwegian Energy Industry
Previous Article in Journal
A PROMETHEE MCDM Application in Social Inclusion: The Case of Foreign-Born Population in the EU
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Channel Identification Based on Cumulants, Binary Measurements, and Kernels

1
Department of Mathematics and Informatics, Sultan Moulay Slimane University, Po. Box 592, 23000 Beni Mellal, Morocco
2
Laboratoire d’Automatique de Caen, UNICAEN, ENSICAEN, Normandie University, 6, B. Marchal Juin, 14050 Caen, France
*
Author to whom correspondence should be addressed.
Systems 2021, 9(2), 46; https://doi.org/10.3390/systems9020046
Submission received: 24 May 2021 / Revised: 11 June 2021 / Accepted: 12 June 2021 / Published: 18 June 2021

Abstract

:
In this paper, we discuss the problem of channel identification by using eight algorithms. The first three algorithms are based on higher-order cumulants, the next three algorithms are based on binary output measurement, and the last two algorithms are based on reproducing kernels. The principal objective of this paper is to study the performance of the presented algorithms in different situations, such as with different sizes of the data input or different signal-to-noise ratios. The presented algorithms are applied to the estimation of the channel parameters of the broadband radio access network (BRAN). The simulation results confirm that the presented algorithms are able to estimate the channel parameters with different accuracies, and each algorithm has its advantages and disadvantages for a given situation, such as for a given S N R and data input. Finally, this study provides an idea of which algorithms can be selected in a given situation. The study presented in this paper demonstrates that the cumulant-based algorithms are more adequate if the data inputs are not available (blind identification), but the kernel- and binary-measurement-based methods are more adequate if the noise is not important ( S N R 16 dB).

1. Introduction

The identification of single-input single-output (SISO) systems in the blind case (without knowledge of the input system data) has been well studied by using higher-order statistics [1,2,3,4,5], without employing any restrictive hypotheses about channel zeros, channel-order overestimation errors, or additive noise color, as well as without increasing the data stream transmission rate. The shared characteristic of this category of techniques [4] is that they use a periodically time-varying precoder to provoke cyclo-stationary statistics at the transmitter (i.e., a simple periodic modulator or a filter bank) and take advantage of the received samples’ second-order cyclo-stationary statistics. Those algorithms have been called transmitter-induced multi-stationary-based techniques. Iterative strategies are used in some types of algorithms for blind channel identification. In these algorithms, an initial channel estimation is employed by using a symbol predictor to obtain provisional soft estimates of the symbol sequence as it is transmitted. A channel estimator uses these estimates to optimize the channel parameters. The symbol estimator then uses the improved channel estimates to optimize the symbol parameters, and so on. The second technique that will be discussed here deals with the identification of the system based on quantized measurement and is used in a wide range of applications, such as microfabricated devices [6,7,8].
The channels under investigation have a non-minimal phase and a selective frequency (i.e., normalized channels for MC-CDMA: BRAN). In the majority of wireless areas, there are numerous hitches in the channels, such as mountains, buildings, and barriers between the transmitter and the receiver. Many different propagation paths are created by reflections from these obstacles. Multi-path propagation, or a multi-path channel, is the term for this. Due to the reverberations and reflection between the transmitter and receiver, the frequency impulse response of this channel is not leveled (ideal case), but contains some cavities and bumps. The synchronization required between the transmitter and the receiver is also a problem in the communication process. To handle this phase estimation problem, some authors have employed higher-order cumulants (HOCs) in order to investigate the reliability of these techniques if the channel is influenced by a colored noise [2,9,10,11,12,13,14,15,16]. HOCs are a significant topic with multiple applications in the field of systems theory, such as in digital signal recognition [17], automatic modulation classification of real signals [18], and space-time block code identification [19].
HOCs are only used for non-linear and non-Gaussian processes because the cumulants of a Gaussian process are identical to zero [2,3,20]. There are many real-world applications that are really non-Gaussian [1,2,3,21]. In addition, the Fourier transforms of HOCs, which are called higher-order spectra (or polyspectra), provide an excellent technique for resolving the problem of equalization technology encountered in communication. The main attribute of HOCs in terms of equalization is that the channel phase information is available [4,5,22,23]. As a result, HOCs can be used to estimate the channel model parameters, even if the phase property (non-minimum phase (NMP) or minimum phase (MP)) of the channel or of the transmitted data is unknown (under the assumption of a non-Gaussian distribution) [1,2,3,24].
A technique founded on binary observation has been employed in a wide range of applications in industries and telecommunication, such as in chemical process sensors for a vacuum, the asynchronous transfer mode (ATM), pressure switches, liquid levels, and industrial sensors for brushless DC motors, shifting using a wire in automotive applications, and switching sensors for exhaust-gas oxygen ABS [7,25,26,27]. One of the most important problems in this wide range of applications is parameter estimation for systems with binary (or quantized) outputs [7,8,28,29,30,31,32]. Some examples of this field include sensor networks, telecommunications, and networked control systems. The third type of technique discussed here in terms of system identification is based on kernel adaptive filtering, which represents an emerging technique for nonlinear signal processing [33], classification, and machine learning (regularization networks [34], support vector machines (SVMs) [35]); the most emblematic attribute is the use of kernel functions, thus permitting the inexpensive operation of the algorithms in high-dimensional spaces.
The kernel methods are based on a central principle called the kernel trick, which gives a non-linear characteristic to many traditionally linear methods, as well as on the theory of reproducing kernel Hilbert spaces (RKHSs); such spaces provide an important framework for elegantly developing adaptive nonlinear filters. So far, several algorithms, including the kernel least mean square (KLMS) [36], kernel affine projection algorithm (KAPA) [37], kernel principal component analysis (KPCA) [38], and kernel recursive least square (KRLS) [39], have been proposed and effectively applied to nonlinear regression problems and classification. In addition, to increase the efficiency of kernel adaptive filtering techniques, some subtypes of these algorithms have also been introduced [40,41,42,43] for channel identification [44] and nonlinear system equalization.
In this paper, we focused on a comparison of the results of channel identification using some algorithms based on HOCs, binary measurement, or kernel-positive methods. We selected three algorithms that use the first two techniques (cumulants and binary measurement) and two algorithms that use the last technique (kernel methods). In order to study the efficiency of these algorithms, we used a practical—i.e., measurable—frequency-selective fading channel called the broadband radio access network (BRAN) (in particular, BRAN A, which is for transmission in indoor scenarios). This model radio channel was normalized by the European Telecommunications Standards Institute (ETSI) in [45,46]. Simulation results were obtained for different values of S N R and different data inputs.
This paper is arranged as follows. In Section 2, we discuss the identification problem. In Section 3, we present three algorithms for channel identification based on higher-order cumulants. In Section 4, we describe the presentation of the binary output measurement algorithms. In Section 5, we illustrate the concept of kernel methods for channel identification and the description of the KRLS and KLMS algorithms. Some simulations for evaluating the performance of the presented algorithms are shown in Section 6 and, finally, Section 7 concludes the paper.

2. Problem Statement

Let us consider the finite impulse response (FIR) system described in Figure 1.
Where, x k and s k are, respectively, the input and the output, y k is an internal signal (unknown for algorithms based on binary output measurement), and b k is an additive noise in this internal signal. We suppose that H has a finite impulse response of length n, where:
y k = i = 1 n h i x k i + b k
y k can be rewritten as:
y k = θ ϕ k + b k
or
y k = ϕ k θ + b k
where ϕ k = x k 1 x k 2 x k n R n and θ = h 1 h 2 h n R n
The binary-valued output, which has a fixed threshold C R , can be represented by:
s k = I ( y k ) = 1 , if y k C 1 , otherwise .

3. Channel Identification Based on Cumulants

In this section, three algorithms based on higher-order cumulants (HOCs) are presented. These algorithms can blindly identify the BRAN channel.

3.1. First Algorithm: Alg1

Safi and Zeroual (2004) demonstrated that they could identify the FIR parameters blindly based only on fourth-order cumulants.
In this paper, we briefly present this algorithm. The mth-order cumulants are [9]:
C m y ( t 1 , , t m 1 ) = γ m x i = 0 q h ( i ) h ( i + t 1 ) h ( i + t m 1 )
where q is the FIR system order and γ m x is the mth-order cumulant—at the origin—of the input FIR system. Based on Equation (5), they demonstrated that the cumulants of order m and n of the FIR output, y ( n ) , and the the impulse response coefficients, h ( i ) , with n > m , are linked by the following relationship:
j = 0 q h ( j ) C n y ( j + t 1 , , j + t m 1 , t m , , t n 1 ) = γ n x γ m x i = 0 q h ( i ) k = m n 1 h ( i + t k ) C m y ( i + t 1 , , i + t m 1 )
Based on Equations (5) and (6), Safi and Zeroual (2004) developed an algorithm (Alg1) based on fourth-order cumulants, which are given as follows:
h ( 0 ) h 2 ( q ) C 4 y ( 0 , 0 , i + α ) = j = 0 q h 3 ( j ) C 4 y ( q , q , j + α )
Based on the hypothesis that the system is causal of order q, j + α is in the interval 0 , q . For this, we have: 0 j + α q j α q j , with 0 j q . From this, two we obtain two inequalities:
q α q
From Equations (7) and (8), they obtained the following system of equations:
0 0 C 4 y ( q , q , 0 ) 0 C 4 y ( q , q , 0 ) C 4 y ( q , q , q ) 0 C 4 y ( q , q , q ) 0 0 × 1 h 2 ( q ) h 3 ( i ) h 2 ( q ) h 3 ( q ) h 2 ( q ) = C 4 y ( 0 , 0 , q ) C 4 y ( 0 , 0 , 0 ) C 4 y ( 0 , 0 , q )
In a compact form, the system of Equations (9) can be written as:
M b q = d
where M, b q , and d are defined in the system of Equations (9). The least squares ( L S ) solution, which is found by using the least squares of the system in Equation (10), is given by
b ^ q = ( M M ) 1 M d
This solution gives us an estimation of the quotient of the parameters h 3 ( i ) and h 3 ( q ) , i.e., b q ( i ) = h 3 ( i ) h 3 ( q ) ^ , i = 1 , , q . So, in order to estimate the parameters h ^ ( i ) , i = 1 , , q , they proceeded as follows:
  • The parameters h ( i ) for i = 1 , , q 1 are estimated from the values of b ^ q ( i ) using the following solution:
    h ^ ( i ) = s i g n b ^ q ( i ) ( b ^ q ( q ) ) 2 a b s ( b ^ q ( i ) ) ( b ^ q ( q ) ) 2 1 / 3
    with s i g n ( x ) = 1 , if x > 0 ; 0 , if x = 0 ; 1 , if x < 0 ;
    and a b s ( x ) = | x | indicates the absolute value of x.
  • The parameters h ^ ( q ) are estimated by:
    h ^ ( q ) = 1 2 s i g n b ^ q ( q ) a b s ( b ^ q ( q ) ) + ( 1 b ^ q ( 1 ) ) 1 / 2

3.2. Second Algorithm: Alg2

Based on Equation (5), the authors of [3] demonstrated the following equation:
j = 0 q h ( j ) C 4 y ( t 1 j , t 2 j , t 3 j ) = γ 4 x γ 2 x ϵ i = 0 q h ( i ) h ( i + t 2 t 1 ) h ( t 3 t 1 ) C 2 y ( t 1 i )
where ϵ = γ 4 x γ 2 x . If t 3 = t 1 , the Equation (14) becomes (with h ( 0 ) = 1 ):
j = 0 q h ( j ) C 4 y ( t 1 j , t 2 j , t 1 j ) = ϵ i = 0 q h ( i ) h ( i + t 2 t 1 ) C 2 y ( t 1 i )
By applying some algebra to Equations (13)–(15), they obtained the following equation:
C 4 y ( 2 q , 2 q , q ) C 4 y ( q , q , 0 ) C 4 y ( 2 q , 2 q , 2 q ) C 4 y ( q , q , q ) × h ( 0 ) h ( q ) = ϵ h ( 0 ) h ( q )
where ϵ = C 4 y ( q , q , q ) C 4 y ( q , 0 , 0 ) C 4 y ( q , q , 0 ) . Taking into account the properties of cumulants in FIR systems, i.e., C m y ( t 1 , , t k , , t m 1 ) = 0 , then t k > q for all k = 1 , , m 1 . Thus, the previous system of Equations (16) becomes:
C 4 y ( 2 q 1 , 2 q 1 , q 1 ) C 4 y ( q , q , 0 ) C 4 y ( 2 q 1 , 2 q 1 , q ) ϵ C 4 y ( q , q , 0 ) 0 0 0 C 4 y ( q , q , q ) ϵ × h ( 1 ) h ( q ) = ϵ C 4 y ( 2 q , 2 q , q ) 0 0
In more compact form, this can be written as:
M h q = d
The parameters h ( i ) , i = 1 , , q , can be estimated using the least squares ( L S ) as follows:
h q = ( M M ) 1 M d

3.3. Third Algorithm: Alg3

In this algorithm, the relationships between the third- and fourth-order cumulants and, eventually, the system parameters are established. Safi and Zeroual (2008) demonstrated the following relationship:
i = 0 q ϵ h ( i ) C 4 y ( t 1 , t 2 i , t 3 i ) = i = 0 q h ( i ) h ( i + t 3 t 2 ) C 3 y ( t 1 , t 2 i )
where ϵ = γ 3 x γ 4 x for different values of t 1 , t 2 , and t 3 Z . Based on this equation, they developed some algorithms.
In order to identify the channel parameters, h ( i ) , i = 1 , , q , they resolved the system of Equation (18) using the least squares (LS) method based on the following procedure.
In order to simplify the estimation of the channel parameters, they supposed that λ = t 1 and t = t 2 = t 3 . Thus, for each λ , based on the L S method, they could estimate the channel parameters from the relationship:
i = 0 q ϵ h ( i ) C 4 y ( λ , t i , t i ) = i = 0 q h 2 ( i ) C 3 y ( λ , t i )
For each λ , they could estimate the parameters h ( i ) using the following equation:
h ^ ( λ ) ( i ) = 1 2 ϵ h ^ ( i ) ϵ ^ + s i g n ( ϵ h ^ ( i ) ϵ ^ ) ( h ^ 2 ( i ) ) 1 / 2
where λ [ q , q ] .

4. Binary Output Measurement Algorithms

4.1. LIMBO Method

The LMS-based identification Method using Binary Observations (LIMBO) [30] is a recursive parameter estimation method that uses binary observations based on the LMS algorithm. The algorithm is based on the following assumption:
  • A.1: θ 1 θ > 0 .
  • A.2: ϕ k is a random process such that:
    The probability density function (pdf) of ϕ k ϕ k is non-zero on the unit sphere.
    ϕ k verifies the α -mixing condition [47].
  • A.3: b k represents the noise, which is uncorrelated with the input sequence.
The binary output defined in Equation (4) is used in the LIMBO algorithm with threshold C = 0 , and the unknown internal signal y k is defined in Equation (3).
The estimated output is:
s ^ k = I ( ϕ k θ ^ k )
where θ ^ k represents the estimation of θ at the time k. The LIMBO algorithm is based on the concept of minimizing the instantaneous error:
E k = 1 2 | s k s ^ k | y ^ k
θ ^ k + 1 is given by:
θ ^ k + 1 = θ ^ k μ α k | s k s ^ k | y ^ k ϕ k ω k
where
  • 0 < μ < 1
  • ω k = 1 2 μ ( 1 μ ) α k | s k s ^ k | y ^ k 2
  • α k = 1 ϕ k ϕ k

4.2. Recursive Least Squares (RLS) Method

The recursive identification algorithm presented in [28] is based on recursive least squares with a forgetting factor (RLS-FF), and it uses the binary output defined in Equation (4) and the threshold C = 0 . The algorithm is based on the following assumptions:
  • B.1: θ = 1 , where . is the l 2 norm.
  • B.2: At any time k, ϕ k θ 0 .
  • B.3: The noise sequence b k is an i.i.d. sequence of random variables with a mean of zero and finite covariance, and it is uncorrelated with the input sequence.
The identification is realized such that:
| s k ϕ k θ | < δ .
with δ < 1 .
The estimation of θ is given by:
θ ^ k = θ ^ k 1 + Γ k ε k / k 1
with
Γ k = P k 1 ϕ k σ k λ + ϕ k P k 1 ϕ k σ k P k = 1 λ ( I n + 1 Γ k ϕ k ) P k 1 ε k / k 1 = s k ϕ k θ ^ k 1 .
where
σ k = λ ϕ k P k 1 ϕ k ( | ε k / k 1 δ | 1 ) , if ( | ε k / k 1 | > δ ) and ( ϕ k P k 1 ϕ k > 0 ) 0 , otherwise .
Finally, 0 < λ < 1 is the forgetting factor.

4.3. Method Based on SVM

The method presented in [29] is based on the use of a supervised learning algorithm—a support vector machine (SVM)—through the formulation of the identification problem as a classification problem.
We assume that:
  • C.1: The static gain of the system is known.
  • C.2: y k is such that at any time k, y k C .
The algorithm has two steps:
Step 1: Find C and the solution θ of
min θ , C , e 1 2 θ 2 + γ k = 1 N e k , s . t s k ( θ ϕ C ) 1 e k , and e k 0 ; k [ 1 , N ] .
Step 2: Normalization of the static gain according to assumption C.1, where γ R 0 + is referred to as the regularization constant or penalty parameter.

5. Kernel-Based Channel Identification

The identification of nonlinear systems using functional estimates based on the notion of Hilbert spaces has led to very interesting results in the last decade. Note that linear problems with a finite order can be taken into consideration to address nonlinear adaptive filtering problems in which the measured input X has been mapped to a nonlinear Hilbert space H of infinite dimensions with a Mercer kernel [48]:
κ ( x i , x j ) = Ψ ( x i ) , Ψ ( x j ) H , ( x i , x j ) X 2
where Ψ maps X into a higher-dimensional space H with an inner product . , . H .
The sequence of samples is transformed by means of a feature map Ψ :
Ψ : X H x κ ( x , . )
To construct the reproducing kernel Hilbert space model, we will use the Gaussian kernel, which is typically a default selection due to its universal approximation property:
κ ( x i , x j ) = exp x i x j 2 2 σ 2 , x i , x j X .
where σ > 0 denotes the kernel bandwidth.
In order to represent the elements of H by using their coordinates, we should give a Hilbert space H an orthonormal base. The associated kernel must be a continuous, symmetric, normalized, and positive definite function κ : X × X R , where the subset X R N is compact.

5.1. Kernel Recursive Least Squares Algorithm

In this subsection, we present the kernel RLS algorithm [39,49,50,51]. It is an algorithm that belongs to the category of recursive least squares algorithms. The basic idea is to run the linear RLS algorithm in the kernel feature space, which is associated with a positive definite kernel κ , via the feature map Ψ ( . ) , which is defined in Equation (30), to generate the input data:
{ ( Ψ ( x 1 ) , s 1 ) , ( Ψ ( x 2 ) , s 2 ) , . . . ( Ψ ( x n ) , s n ) , . . . } .
The optimization problem of the kernel RLS algorithm can be expressed as:
min Ψ H k = 0 n 1 λ n k | s k Ψ ( x k ) | 2 + λ k Ψ H 2
where λ is a positive regularization parameter, s k is the desired response of the system, and Ψ ( x k ) designates the corresponding output of the system for x k .
According to the representation theorem [52], the solution to the optimization problem (33) can be written as:
Ψ ( . ) n = i = 1 m θ n , k κ ( . , x k )
where θ n = ( θ n , 1 , θ n , 2 , . . . , θ n , n ) is the model coefficient vector. By introducing Equation (34) into Equation (33), the optimal solution to the problem yields the following:
min θ Λ s k H k θ 2 + λ n θ K k θ
where ( H k ) k , i = κ ( x k , x i ) , for k , i = 1 , 2 , . . . , n , K k = ( κ ( x 1 , x n ) , κ ( x 2 , x n ) , . . . , κ ( x m , x n ) ) , and Λ is a diagonal matrix of size ( n , n ) , of which the ( k , k ) th element is λ n k . If we consider that H k is invertible, the solution of the problem (35) becomes:
θ n = P n H k Λ n s n
with P n = ( H k Λ n H k + λ n K k ) 1 . At the instant n + 1 , a new x n + 1 observation is detected at the input of the system, and we update H n + 1 and s n + 1 with Equations (37) and (38), respectively:
H n + 1 = H n h n + 1
s n + 1 = ( s 0 , s 2 , . . . , s n )
where h n + 1 = ( κ ( x n + 1 , x 1 ) , κ ( x n + 1 , x 2 ) , . . . , κ ( x n + 1 , x m ) ) . Then, the equations of the kernel RLS algorithm are updated according to:
G n + 1 = λ 1 P n + 1 h n + 1 1 + λ 1 h n + 1 P n + 1 h n + 1
θ n + 1 = θ n + G n + 1 ( s n + 1 h n + 1 θ n )
P n + 1 = λ 1 P n 1 G n + 1 h n + 1
where the term s n + 1 h n + 1 θ n represents the estimation error.

5.2. Kernel Least Mean Square Algorithm

In this paragraph, we describe the kernel LMS algorithm [36,53]. This algorithm belongs to the category of stochastic gradient algorithms. The fundamental concept is the execution of the linear least mean square algorithm in the kernel feature map, which is associated with the positively defined kernel κ via the feature map Ψ ( . ) to generate the input data. The update equations of the LMS algorithm are given by [54]:
e n = s n θ n 1 x n
θ n = θ n 1 + μ e n x n
where e n denotes the estimation error, θ n represents the weight vector estimate, and μ stands for the step parameter.
Let us now suppose that the Ψ function map is used to transform the sample sequence; then, we perform the least mean square procedure for the processed data sequence defined in Equation (32), so we obtain the kernel LMS:
e n = s n θ n 1 Ψ ( x n )
θ n = θ n 1 + μ e n Ψ ( x n )
The biggest distinction of LMS is when it is used in a space with potentially infinite-dimensional characteristics in Equation (45), which makes direct updating virtually impossible. Rather, we will use every θ n to relate to the initialization θ 0 :
θ n = θ 0 + μ i = 1 n e i Ψ ( x i ) .
θ n = μ i = 1 n e i Ψ ( x i ) . ( if we suppose that θ 0 = θ ( 0 ) = 0 )
The prediction solution is computed by using the kernel trick:
θ n , Ψ ( x n ) H = μ i = 1 n e i Ψ ( x i ) , Ψ ( x n ) = μ i = 1 n e i κ ( x i , x n ) ,
where κ ( x i , x n ) is a Mercer kernel and n is the number of training samples.

6. Simulation Results

In this section, we will evaluate the efficiency of the presented algorithms for various S N R values and for a fixed data input N.

6.1. Impulse Response Parameter Estimation

In Figure 2, it can be seen that Alg1, Alg2, and Alg3 gave good results for S N R = 0 dB and a data input of N = 2000 . These estimations of the channel parameters were closer to the measured values if the number of units of data was important and for S N R > 16 dB; this is due to an important property of the cumulant: The cumulant of a Gaussian process is zero. In the following (Figure 3 and Figure 4), we present the estimated parameters in the cases of S N R = 16 dB and 32 dB, respectively, with a fixed data input of N = 2000 .
Next, we tested the performance of the algorithms based on binary measurement for different values of S N R and for a fixed data input of N = 2000 in order to make a comparison with the algorithms based on cumulants.
In Figure 5, we present the results obtained for S N R = 0 dB and N = 2000 when using the methods based on binary measurement. We can observe in Figure 5 that the channel parameters estimated based on the SVM method were closer to the measured values than those obtained by using the RLS and LIMBO methods, which gave values that were very far from those measured.
In Figure 6, we show the estimations of the BRAN A parameters by using the RLS, LIMBO, and SVM methods in the case of S N R = 16 dB and with a data input of N = 2000 . In this case where the noise is not very important, we remark that the results obtained using the RLS and LIMBO methods were closer to those obtained in the case in which noise was important ( S N R = 0 dB); this result can be explained by the fact that the rounding errors due to the noise were minimal, which implies an improvement in the results obtained.
If the noise was lower (i.e., the S N R was high), such as with S N R = 32 dB, we observed that the estimated channel parameters were closer to the measured values (Figure 7), but if the noise was more important, the estimated channel parameters were significantly different from the measured values. Note that if the noise was equal to the power of the information signal (i.e., S N R = 0 dB), we could see that the estimated channel parameters were very different from the measured values; this is because of the lack of a filtering operation for the noise, which implies some rounding errors. However, the SVM method presented good results for different values of S N R ; this was due to the improved classification performance, as well as the ability to manipulate a large volume of data with a reasonable computational load. Furthermore, this strategy was one that was perfect for minimizing prediction error. However, due to its simpler and more efficient structure and learning procedure, the SVM has a second characteristic: It can learn efficiently with small amounts of data.
Figure 8, Figure 9 and Figure 10 show the estimation of the impulse response parameters of the BRAN A channel as a function of the path delays by using the kernel-based algorithms for different values of S N R and N = 2000 .
From Figure 8, Figure 9 and Figure 10, we can conclude that the estimated amplitude of the BRAN A channel was very close to the true value, especially for S N R 16 dB (see Figure 9 and Figure 10). Regarding the impulse response of the channel in a very noisy environment (Figure 8), the Gaussian noise influenced the estimated model parameters when using KLMS, and the noise slightly influenced the estimation of the impulse response parameters, especially when KRLS was used.

6.2. Magnitude and Phase Estimation

Next, we identified the amplitude and phase of the ETSI BRAN A channel for different values of S N R by setting the data input to N = 2000 .
Figure 11, Figure 12 and Figure 13 represent the results of the estimation of the amplitude and phase by comparing the measured parameters of the BRAN A channel and the parameters estimated using the algorithms based on cumulants (Alg1, Alg2, and Alg3). It should be noted that the algorithms based on cumulants provided remarkable performance. Indeed, if S N R < 16 dB, we obtained good results with the algorithms based on cumulants, but if S N R 16 dB, we observed that the algorithms based on binary measurement gave us better results than the algorithms based on cumulants. It must be emphasized that the algorithms based on cumulants are simple and have the advantage of identifying the parameters with a low calculation cost and without knowledge of the input (blind identification) in comparison with the methods based on binary measurement, which require the use of the input and output of the channel.
The estimations of the amplitude and phase of the BRAN A channel using the methods based on binary measurement with different S N R values are presented in Figure 14, Figure 15 and Figure 16. From these figures, we observe that the shapes of the amplitude and phase estimated with the SVM algorithm are the same as those of the measured data (Figure 16). We noticed a minor difference between the estimated and measured values for S N R < 16 dB (Figure 14 and Figure 15). These results are very important because we are using a very quickly fading channel.
The curves in Figure 17, Figure 18 and Figure 19 show the estimations of the amplitude and phase of the BRAN A channel by using kernel methods (KRLS and KLMS) for a data length of N = 2000 and for different values of S N R . In Figure 17, Figure 18 and Figure 19, we present the performance of the KRLS and KLMS algorithms. We observed that the KRLS algorithm performed better than the KLMS algorithm, which produced quite large errors. We can conclude from Figure 17 that the phase of the BRAN A channel’s impulse response estimated by the KRLS algorithm was closer to the true values. However, when using the KLMS algorithm, the difference between the measured and estimated phase was significantly larger.
The estimated parameters shown in Figure 12, Figure 13, Figure 15, Figure 16, Figure 18, and Figure 19 evidently show that both KLMS and KRLS achieved higher prediction accuracy than the binary-measurement-based algorithms and cumulant-based algorithms for all values of S N R . Meanwhile, the performance of KRLS was also marginally better than that of KLMS, although it had a greater computational complexity.

7. Conclusions

In this paper, we presented comparisons among three algorithms based on binary measurements, three algorithms based on higher-order cumulants, and two algorithms based on reproducing kernel Hilbert spaces. These algorithms were used to estimate the parameters of a broadband radio access network channel (BRAN A). Results were obtained for different S N R values and for different data inputs. From this study, we can conclude that both of the kernel algorithms are efficient and effective in identifying the amplitude and phase of the BRAN A channel compared to the other algorithms (those based on higher-order cumulants and binary measurements). In general, the algorithms based on binary measurement gave better results in the estimation of parameters if the S N R is higher than 16 dB (i.e., S N R 16 dB), but if the noise was more present (i.e., S N R < 16 dB), the algorithms based on cumulants gave better results than those obtained by algorithms based on binary measurement; this was because the cumulants of the Gaussian process (here, we supposed that the noise had a Gaussian distribution) were null. We can add that if the size of the data input is important, the estimation of the cumulants will be approximately without bias, which implies a good estimation of the channel parameters. Future work will focus on developing techniques for identifying communication channels by using kernel methods.

Author Contributions

Conceptualization, H.O., R.F. and A.D.; Data curation, R.F. and A.D.; Formal analysis, H.O. and R.F.; Investigation, R.F., A.D. and S.S.; Methodology, H.O., R.F., A.D. and S.S.; Resources, M.P. and M.F.; Software, R.F. and S.S.; Supervision, M.P. and M.F.; Validation, A.D. and S.S.; Visualization, H.O., R.F., M.P. and M.F.; Writing—original draft, R.F., H.O., A.D. and S.S.; Writing—review & editing, S.S., M.P. and M.F.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Safi, S.; Zeroual, A. MA system identification using higher order cumulants: Application to modelling solar radiation. Int. J. Stat. Comput. Simul. 2002, 72, 533–548. [Google Scholar] [CrossRef]
  2. Safi, S.; Zeroual, A. Blind parametric identification of non-Gaussian FIR systems using higher order cumulants. Int. J. Syst. Sci. 2004, 35, 855–867. [Google Scholar] [CrossRef]
  3. Safi, S.; Zeroual, A. Blind non-minimum phase channel identification using 3rd and 4th order cumulants. Int. J. Signal Process. 2008, 4, 158–168. [Google Scholar]
  4. Sadler, B.M.; Giannakis, G.B.; Lii, K.S. Estimation and detection in non-Gaussian noise using higher order statistics. IEEE Trans. Signal Process. 1994, 42, 2729–2741. [Google Scholar] [CrossRef]
  5. Wang, J.; He, Z. Criteria and algorithms for blind source separation based on cumulants. Int. J. Electron. 1996, 81, 1–14. [Google Scholar] [CrossRef]
  6. Yin, G.G.; Zhao, Y.; Zhang, J.F. Identification input design for consistent parameter estimation of linear systems with binary-valued output observations. IEEE Trans. Autom. Control 2008, 53, 867–880. [Google Scholar]
  7. Zhang, J.F.; Yin, G.G. System identification using binary sensors. IEEE Trans. Autom. Control 2003, 48, 1892–1907. [Google Scholar]
  8. Godoy, B.I.; Goodwin, G.C.; Agüero, J.C.; Marelli, D.; Wigren, T. On identification of FIR systems having quantized output data. Automatica 2011, 47, 1905–1915. [Google Scholar] [CrossRef]
  9. Safi, S.; Zeroual, A. Blind identification in noisy environment of nonminimum phase finite impulse response (FIR) system using higher order statistics. Syst. Anal. Model. Simul. 2003, 43, 671–681. [Google Scholar] [CrossRef]
  10. Antari, J.; Iqdour, R.; Safi, S.; Zeroual, A.; Lyhyaoui, A. Identification of quadratic non linear systems using higher order statistics and fuzzy models. In Proceedings of the IEEE International Conference On Acoustic, Speech and Signal Process (ICASSP), Toulouse, France, 14–19 May 2006; pp. 712–715. [Google Scholar]
  11. Zidane, M.; Safi, S.; Sabri, M. Extending HOC-based methods for identifying the diagonal parameters of quadratic systems. Signal Image Video Process. 2018, 12, 125–132. [Google Scholar] [CrossRef]
  12. Zidane, M.; Safi, S.; Sabri, M.; Boumezzough, A. Identification and equalization using higher order cumulants in MC-CDMA systems. In Proceedings of the 2014 5th Workshop on Codes, Cryptography and Communication Systems, WCCCS 2014, El Jadida, Morocco, 27–28 November 2014; pp. 81–85. [Google Scholar]
  13. Zidane, M.; Safi, S.; Sabri, M.; Boumezzough, A.; Frikel, M. Broadband radio access network channel identification and downlink MC-CDMA equalization. Int. J. Energy Inform. Commun. 2014, 5, 13–34. [Google Scholar] [CrossRef]
  14. Frikel, M.; Targui, B.; Hamon, F.; M’Saad, M. Adaptive equalization using controlled equal gain combining for uplink/downlink MC-CDMA systems. Int. J. Signal Process. 2008, 4, 230–237. [Google Scholar]
  15. Zidane, M.; Safi, S.; Sabri, M. Measured and estimated data of non-linear BRAN channels using HOS in 4G wireless communications. Data Brief 2018, 17, 1136–1148. [Google Scholar] [CrossRef]
  16. Nikias, C.L.; Mendel, J.M. Signal processing with higher-order spectra. IEEE Signal Process. Mag. 1993, 10, 10–37. [Google Scholar] [CrossRef]
  17. Wang, A.; Li, R. Research on digital signal recognition based on higher order cumulants. In Proceedings of the 2019 International Conference on Intelligent Transportation, Big Data & Smart City (ICITBS), Changsha, China, 12–13 January 2019; pp. 586–588. [Google Scholar]
  18. Simic, M.; Stankovic, M.; Orlic, V.D. Automatic Modulation Classification of Real Signals in AWGN Channel Based on Sixth-Order Cumulants. Radioengineering 2021, 30, 205. [Google Scholar] [CrossRef]
  19. Khosraviyani, M.; Kalbkhani, H.; Shayesteh, M.G. Higher order statistics for modulation and STBC recognition in MIMO systems. IET Commun. 2019, 13, 2436–2446. [Google Scholar] [CrossRef]
  20. Anderson, J.M.; Giannakis, G.B. Noisy input/output system identification using cumulants and the Steiglitz-McBride algorithm. IEEE Trans. Signal Process. 1996, 44, 1021–1024. [Google Scholar] [CrossRef]
  21. Kaiser, T. MA-model identification using modulated moment sequences. Signal Process. 1995, 47, 85–103. [Google Scholar] [CrossRef]
  22. Proakis, J.; Salehi, M. Digital Communications, 5th ed.; McGraw-Hill: New York, NY, USA, 2007. [Google Scholar]
  23. Ju, L.; Zhenya, H. Blind identification and equalization using higher-order cumulants and ICA algorithms. In Proceedings of the International Conference Neural Networks and Brain (ICNN B’98), Beijing, China, 13–15 October 1998. [Google Scholar]
  24. Zhang, X.D.; Zhang, Y.S. FIR system identification using higher order statistics alone. IEEE Trans. Signal Process. 1994, 42, 2854–2858. [Google Scholar] [CrossRef]
  25. Guo, J.; Zhao, Y.; Sun, C.Y.; Yu, Y. Recursive identification of FIR systems with binary-valued outputs and communication channels. Automatica 2015, 60, 165–172. [Google Scholar] [CrossRef]
  26. Guo, J.; Yin, G.; Zhao, Y.; Zhang, J.F. Identification of Wiener systems with quantized inputs and binary-valued output observations. Automatica 2017, 78, 280–286. [Google Scholar] [CrossRef]
  27. Yin, G.G.; Zhang, J.F. Joint identification of plant rational models and noise distribution functions using binary-valued observations. Automatica 2006, 42, 535–547. [Google Scholar]
  28. Pouliquen, M.; Menard, T.; Pigeon, E.; Gehan, O.; Goudjil, A. Recursive system identification algorithm using binary measurements. In Proceedings of the 2016 European Control Conference (ECC), Aalborg, Denmark, 29 June–1 July 2016; pp. 1353–1358. [Google Scholar]
  29. Goudjil, A.; Pouliquen, M.; Pigeon, E.; Gehan, O.; M’Saad, M. Identification of systems using binary sensors via support vector machines. In Proceedings of the 2015 54th IEEE Conference on Decision and Control (CDC), Osaka, Japan, 15–18 December 2015; pp. 3385–3390. [Google Scholar]
  30. Jafari, K.; Juillard, J.; Roger, M. Convergence analysis of an online approach to parameter estimation problems based on binary observations. Automatica 2012, 48, 2837–2842. [Google Scholar] [CrossRef] [Green Version]
  31. Guo, J.; Zhao, Y. Recursive projection algorithm on FIR system identification with binary-valued observations. Automatica 2013, 49, 3396–3401. [Google Scholar] [CrossRef]
  32. Oualla, H.; Pouliquen, M.; Frikel, M.; Safi, S. Comparison of algorithms for identification of IIR systems from binary measurements on the output. In Proceedings of the E3S Web of Conferences, EDP Sciences, Guizhou, China, 23–25 July 2021; p. 01049. [Google Scholar]
  33. Liu, W.; Principe, J.C.; Haykin, S. Kernel Adaptive Filtering: A Comprehensive Introduction; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  34. Girosi, F.; Jones, M.; Poggio, T. Regularization theory and neural networks architectures. Neural Comput. 1995, 7, 219–269. [Google Scholar] [CrossRef]
  35. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–6297. [Google Scholar] [CrossRef]
  36. Liu, W.; Pokharel, P.P.; Principe, J.C. The kernel least-mean-square algorithm. IEEE Trans. Signal Process. 2008, 56, 543–554. [Google Scholar] [CrossRef]
  37. Liu, W.; Principe, J.C. Kernel affine projection algorithms. EURASIP J. Adv. Signal Process. 2008, 2008, 1–12. [Google Scholar] [CrossRef] [Green Version]
  38. Schölkopf, B.; Smola, A.; Müller, K.R. Nonlinear component analysis as a kernel eigenvalue problem. Neural Comput. 1998, 10, 1299–1319. [Google Scholar] [CrossRef] [Green Version]
  39. Engel, Y.; Mannor, S.; Meir, R. The kernel recursive least-squares algorithm. IEEE Trans. Signal Process. 2004, 52, 2275–2285. [Google Scholar] [CrossRef]
  40. Chen, B.; Zhao, S.; Zhu, P.; Príncipe, J.C. Quantized kernel least mean square algorithm. IEEE Trans. Neural Netw. Learn. Syst. 2011, 23, 22–32. [Google Scholar] [CrossRef]
  41. Wang, S.; Zheng, Y.; Ling, C. Regularized kernel least mean square algorithm with multiple-delay feedback. IEEE Signal Process. Lett. 2015, 23, 98–101. [Google Scholar] [CrossRef]
  42. Chen, B.; Liang, J.; Zheng, N.; Príncipe, J.C. Kernel least mean square with adaptive kernel size. Neurocomputing 2016, 191, 95–106. [Google Scholar] [CrossRef] [Green Version]
  43. Slavakis, K.; Bouboulis, P.; Theodoridis, S. Online learning in reproducing kernel Hilbert spaces. In Signal Processing Theory and Machine Learning (Series Academic Press Library in Signal Processing); Chellappa, R., Theodoridis, S., Eds.; Academic: New York, NY, USA, 2014; pp. 883–987. [Google Scholar]
  44. Fateh, R.; Darif, A. Mean Square Convergence of Reproducing Kernel for Channel Identification: Application to Bran D Channel Impulse Response. In International Conference on Business Intelligence; Springer: Cham, Switzerland, 2021; pp. 284–293. [Google Scholar]
  45. Broadband Radio Acces Network (BRAN). Hyperlan type 2. In Requirements and Archtiectures for Wireless Broadband Access; European Telecommunications Standards Institute (ETSI): Sophia Antipolis, France, 1999. [Google Scholar]
  46. Broadband Radio Acces Network (BRAN). Hyperlan type 2. In Physical Layer; European Telecommunications Standards Institute (ETSI): Sophia Antipolis, France, 2001. [Google Scholar]
  47. Rosenblatt, M. A central limit theorem and a strong mixing condition. Proc. Natl. Acad. Sci. USA 1956, 42, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Aronszajn, N. Theory of reproducing kernels. Trans. Am. Math. Soc. 1950, 68, 337–404. [Google Scholar] [CrossRef]
  49. Tobar, F. Kernel-Based Adaptive Estimation: Multidimensional and State-Space Approaches. Ph.D. Thesis, Imperial College London, London, UK, 2014. [Google Scholar]
  50. Façanha, T.S.; Barreto, G.A.; Costa Filho, J.T. A Novel Kalman Filter Formulation for Improving Tracking Performance of the Extended Kernel RLS. Circuits Syst. Signal Process. 2021, 40, 1397–1419. [Google Scholar] [CrossRef]
  51. Saide, C.; Lengelle, R.; Honeine, P.; Richard, C.; Achkar, R. Nonlinear adaptive filtering using kernel-based algorithms with dictionary adaptation. Int. J. Adapt. Control. Signal Process. 2015, 29, 1391–1410. [Google Scholar] [CrossRef] [Green Version]
  52. Schölkopf, B.; Herbrich, R.; Smola, A.J. A generalized representer theorem. In International Conference on Computational Learning Theory; Springer: Berlin/Heidelberg, Germany, 2001; pp. 416–426. [Google Scholar]
  53. Li, K.; Principe, J.C. No-Trick (Treat) Kernel Adaptive Filtering Using Deterministic Features. 2019. Available online: https://arxiv.org/abs/1912.04530 (accessed on 16 June 2021).
  54. Martinek, R.; Rzidky, J.; Jaros, R.; Bilik, P.; Ladrova, M. Least Mean Squares and Recursive Least Squares Algorithms for Total Harmonic Distortion Reduction Using Shunt Active Power Filter Control. Energies 2019, 12, 1545. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Block diagram of the system model.
Figure 1. Block diagram of the system model.
Systems 09 00046 g001
Figure 2. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 0 dB and N = 2000 .
Figure 2. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 0 dB and N = 2000 .
Systems 09 00046 g002
Figure 3. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 16 dB and N = 2000 .
Figure 3. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 16 dB and N = 2000 .
Systems 09 00046 g003
Figure 4. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 32 dB and N = 2000 .
Figure 4. Parameters of BRAN A estimated using algorithms based on cumulants for S N R = 32 dB and N = 2000 .
Systems 09 00046 g004
Figure 5. Parameters of BRAN A estimated using RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) and SVM ( γ = 0.1 ) method for S N R = 0 dB and N = 2000 .
Figure 5. Parameters of BRAN A estimated using RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) and SVM ( γ = 0.1 ) method for S N R = 0 dB and N = 2000 .
Systems 09 00046 g005
Figure 6. Parameters of BRAN A estimated using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 16 dB and N = 2000 .
Figure 6. Parameters of BRAN A estimated using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 16 dB and N = 2000 .
Systems 09 00046 g006
Figure 7. Parameters of BRAN A estimated using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 32 dB and N = 2000 .
Figure 7. Parameters of BRAN A estimated using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 32 dB and N = 2000 .
Systems 09 00046 g007
Figure 8. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 0 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 8. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 0 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g008
Figure 9. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 16 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 9. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 16 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g009
Figure 10. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 32 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 10. Parameters of BRAN A estimated using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 32 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g010
Figure 11. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 0 dB and N = 2000 .
Figure 11. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 0 dB and N = 2000 .
Systems 09 00046 g011
Figure 12. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 16 dB and N = 2000 .
Figure 12. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 16 dB and N = 2000 .
Systems 09 00046 g012
Figure 13. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 32 dB and N = 2000 .
Figure 13. Estimations of the BRAN A magnitude and phase using algorithms based on cumulants for S N R = 32 dB and N = 2000 .
Systems 09 00046 g013
Figure 14. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 0 dB and N = 2000 .
Figure 14. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 0 dB and N = 2000 .
Systems 09 00046 g014
Figure 15. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 16 dB and N = 2000 .
Figure 15. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 16 dB and N = 2000 .
Systems 09 00046 g015
Figure 16. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 32 dB and N = 2000 .
Figure 16. Estimations of the BRAN A magnitude and phase using the RLS ( λ = 0.999 , δ = 0.9999 ) , LIMBO ( μ = 0.2 ) , and SVM ( γ = 0.1 ) methods for S N R = 32 dB and N = 2000 .
Systems 09 00046 g016
Figure 17. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 0 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 17. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 0 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g017
Figure 18. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 16 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 18. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 16 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g018
Figure 19. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 32 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Figure 19. Estimations of the BRAN A magnitude and phase using KRLS ( λ = 0.99 ) and KLMS ( μ = 0.05 ) for S N R = 32 dB and N = 2000 with the Gaussian kernel ( σ = 0.5 ) .
Systems 09 00046 g019
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oualla, H.; Fateh, R.; Darif, A.; Safi, S.; Pouliquen, M.; Frikel, M. Channel Identification Based on Cumulants, Binary Measurements, and Kernels. Systems 2021, 9, 46. https://doi.org/10.3390/systems9020046

AMA Style

Oualla H, Fateh R, Darif A, Safi S, Pouliquen M, Frikel M. Channel Identification Based on Cumulants, Binary Measurements, and Kernels. Systems. 2021; 9(2):46. https://doi.org/10.3390/systems9020046

Chicago/Turabian Style

Oualla, Hicham, Rachid Fateh, Anouar Darif, Said Safi, Mathieu Pouliquen, and Miloud Frikel. 2021. "Channel Identification Based on Cumulants, Binary Measurements, and Kernels" Systems 9, no. 2: 46. https://doi.org/10.3390/systems9020046

APA Style

Oualla, H., Fateh, R., Darif, A., Safi, S., Pouliquen, M., & Frikel, M. (2021). Channel Identification Based on Cumulants, Binary Measurements, and Kernels. Systems, 9(2), 46. https://doi.org/10.3390/systems9020046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop