Next Article in Journal
Electrochemical DNA-Sensor Based on Macrocyclic Dendrimers with Terminal Amino Groups and Carbon Nanomaterials
Next Article in Special Issue
Terahertz Meets AI: The State of the Art
Previous Article in Journal
Research and Application of Contactless Measurement of Transformer Winding Tilt Angle Based on Machine Vision
Previous Article in Special Issue
UAV-Aided Dual-User Wireless Power Transfer: 3D Trajectory Design and Energy Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Parameter Estimation Method and Closed-Form Solution Study for k-µ Channel Model

1
Institute of Electronic Engineering, China Academy of Engineering Physics, Mianyang 621999, China
2
High-Tech Institute, The First School, Rocket Force University of Engineering, Xi’an 710025, China
3
The School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2023, 23(10), 4760; https://doi.org/10.3390/s23104760
Submission received: 6 February 2023 / Revised: 27 April 2023 / Accepted: 8 May 2023 / Published: 15 May 2023
(This article belongs to the Special Issue System Design and Signal Processing for 6G Wireless Communications)

Abstract

:
This paper proposes a novel multi-parameter estimation algorithm for the k-µ fading channel model to analyze wireless transmission performance in complex time-varying and non-line-of-sight communication scenarios involving moving targets. The proposed estimator offers a mathematically tractable theoretical framework for the application of the k-µ fading channel model in realistic scenarios. Specifically, the algorithm obtains expressions for the moment-generating function of the k-µ fading distribution and eliminates the gamma function using the even-order moment value comparison method. It then obtains two sets of solution models for the moment-generating function at different orders, which enable the estimation of the k and µ parameters using three sets of closed-form solutions. The k and µ parameters are estimated based on received channel data samples generated using the Monte Carlo method to restore the distribution envelope of the received signal. Simulation results show strong agreement between theoretical and estimated values for the closed-form estimated solutions. Additionally, the differences in complexity, accuracy exhibited under different parameter settings, and robustness under decreasing SNR may make the estimators suitable for different practical application scenarios.

1. Introduction

In wireless communication systems, the quality of the channel plays a significant role in determining the overall communication quality. The channel quality indicator (CQI) is a critical parameter that determines the maximum achievable data rate of the system [1,2]. Advances in wireless communication technology have been focused on improving channel estimation performance and enhancing modulation schemes. Currently, 5G networks are still under active development and construction. To meet the high-performance requirements of 5G communication, such as achieving Gbps-level high-speed transmission, low end-to-end transmission latency of 1 ms, and high-bandwidth instant burst transmission, the communication system requires measuring CQI and analyzing it based on the channel state information. This analysis provides theoretical support for selecting system parameters such as modulation scheme, data block size, coding method, and data rate, to optimize the system’s throughput, latency, and bandwidth performance.
Accurately reflecting changes in the physical quantity of the channel in real time is a crucial aspect of the channel quality indicator (CQI). Commonly used metrics for assessing channel quality include the bit error rate (BER) and signal-to-interference-plus-noise ratio (SINR) [3], etc. These metrics are quantified factors calculated by the user through signal reception and processing, and their accuracy depends on the quality of the channel estimation. Since signal energy propagation in wireless channels is a stochastic process, channel characteristics are typically described using probability distribution models that can adapt to measured data in real-world environments. To model the channel, the usual approach is to estimate specific model parameters from the received measurement data and then evaluate the fit between the data and the model to determine the channel characteristics. Furthermore, the time-varying nature of channels in mobile communication necessitates that the channel estimation process controls computational delay to effectively reduce the real-time impact of measurement cycles on channel evaluation. Therefore, determining the probability distribution model depends not only on the sample size but also on the parameters to be estimated, the selected distribution, and the parameter estimation method.
However, in indoor communication, the highly dynamic and complex nature of channels makes it difficult to accurately analyze channel characteristics. This is mainly due to the rich transmission multipath and random shadowing effects caused by complex environments such as indoor and cabin spaces, which result in a large number of multipath propagation components in wireless signals between communication transceivers, and the LOS path cannot be guaranteed to exist. As a result, users have difficulty determining the typical channel types of signal fading and distribution changes in different positions while moving. This leads to a decrease in the performance of channel estimators and the problem of mismatch between the estimated and measured values of the channel estimation curve, which can seriously affect system performance prediction, adaptive modulation strategy generation, and transmission capacity optimization. Therefore, a more generalized representation of channel models is needed to address these challenges.
As research has progressed, new channel models have been proposed, including the k-µ distribution and k-µ fading distribution models, which are more generalized and versatile models for channel fading. The k-µ distribution model, in particular, has garnered attention due to its good agreement with measured channel data and its ability to degrade to different typical channel fading models (e.g., Rayleigh, Rice, and Nakagami-m channels) under certain parameter constraints [4], the development of typical fading channel models can be seen in Table 1. Since the estimation of channel models plays a crucial role in predicting the performance of 5G communication systems, it is important to study the k-µ fading channel model and develop reliable and effective algorithms for estimating its parameters, which can provide a new approach for modeling channels. This is essential to analyzing the theoretical performance limits of communication systems operating in k-µ fading channels and providing a solid theoretical basis for system design.
In recent years, there have been many studies on channel performance analysis and system performance optimization of typical fading channel models both domestically and internationally, as reported in [1,11,12,13]. However, there have been relatively few studies on parameter estimation of the k-µ distribution and k-µ shadowing distribution models [14,15,16,17]. In [14], the authors analyzed the error vector magnitude (EVM) performance of the k-µ fading channel, where the EVM parameter was used to evaluate the variation of channel quality, and the results could be used to predict the lower bound of channel quality based on k-µ fading channels. In [15], the authors assumed QAM modulation channels passing through k-µ shadow fading channels and proposed the average symbol error rate (ASER) expressions for orthogonal amplitude modulation (RQAM) and cross-orthogonal amplitude modulation (XQAM) in single-input single-output (SISO) and multiple-input multiple-output (MIMO) channel systems. In [16], the transmission characteristics of different receivers in k-µ shadowing fading channels were studied. The average bit error rate, outage probability, and channel capacity under three different receivers, namely, N-branch selection combining, dual-branch selection combining, and dual-branch maximal ratio combining, were derived. In [17], the authors used the k-µ fading channel model to describe complex data transmission scenarios in a bidirectional relay network, derived the outage probability under k-µ fading channels, and also derived the outage probability under various mixed fading channels.
In addition, the k-μ distribution and k-μ shadowing distribution have been widely used in the performance analysis of wireless communication and network systems. Through detailed analysis of the higher-order statistical characteristics of the frequency-agile system in the k-μ and k-μ shadow fading channels, reference [18] derived novel and accurate expressions for the high-order cumulants and moments of channel capacity under the generalized fading model, as well as simplified high-order statistical expressions for channel capacity in extremely low and high signal-to-noise ratio conditions. Reference [19] considered a master-slave and cognitive relay spectrum sharing network in which the channel follows a k-μ shadow fading distribution, derived the expression for the interruption probability of the cognitive relay network under the peak interference power constraint of the master user network, and analyzed the system performance of the cognitive relay network under different fading and shadowing parameters and power constraint conditions. Reference [20] studied the performance of energy-based spectrum sensing in k-μ and k-μ extreme fading channels under the Bayesian framework, and derived accurate expressions for the channel model error probability (as a convex combination of false alarm probability and signal detection probability). Reference [21] derived a novel closed-form analytical solution for the average quality metric (average receiver operating characteristic curve and area under the curve) of a generic k-μ shadow fading channel model by analyzing the average quality metric describing the energy detection process; moreover, to investigate the impact of human body occlusion on the channel in cellular device-to-device communication, reference [22] derived the probability density function and moment function of the channel model and applied it to link communication between devices in an outdoor urban environment at 868 MHz, and found that the shadow effect of the deterministic component varies significantly depending on the relative position of the user device to the body and link geometry. Reference [23] derived approximate expressions for the system interruption probability and channel capacity in a system where both the sensing transmission signal and the interference signal experience k-μ shadow fading under given conditions such as signal-to-noise ratio. These expressions are expressed in known hypergeometric functions and can be easily calculated numerically. However, while most of these studies evaluate channel transmission quality in specific scenarios through the analysis of system interruption probability, average bit error rate, and channel capacity based on k-μ distribution and k-μ shadow distribution, they do not provide effective channel model parameter estimation methods by assuming the range of critical channel parameters such as k and µ, making it difficult to provide adjustment strategies for communication system transmission scheme design and performance analysis optimization through accurate channel estimation and CQI quantization in practical scenarios, thus affecting system generality.
In the existing research on k-μ parameter estimators, Yacoub et al. [24] proposed closed-form expressions for k-μ parameters, and Ribeiro et al. further analyzed the asymptotic validity of this estimator [25]. However, our experiments have shown that this estimator is inefficient and prone to errors in estimation, which means the method in [24] has a higher probability of error when estimating channel parameters. [25] includes transcendental functions that incorporate special functions such as gamma and hypergeometric functions, which means that they cannot be solved analytically, but only through numerical methods. Additionally, a maximum likelihood-based α-η-μ distribution parameter estimator was proposed in [26] for inferring the characteristics of signal amplitude distribution. However, this study’s conclusion is not a closed-form expression, and the algorithm’s complexity is high. To address these issues, this paper proposes a new k-μ parameter estimator for the communication system’s k-µ fading channel model, providing theoretical support for distinguishing fading channel types. The main contributions of this paper are as follows: First, two mathematical models containing the unknown parameters k and µ are established using the moment-generating function. Second, three k-μ parameter moment estimators based on the k-µ fading model are derived, resulting in corresponding closed-form expressions for parameter estimation. The complexity and estimation accuracy of the estimator are analyzed. Finally, simulation analysis is performed on different parameters, ensuring high-precision and effective estimation of k and µ parameters under different k-value conditions.
The structure of this paper is organized as follows: In Section 2, we introduce the k-μ distribution-based fading channel model and the existing estimators. In Section 3, we provide the derivation process for three moment-based estimators. In Section 4, we utilize Monte Carlo simulation to generate test data and compare the computational complexity and estimation accuracy of the three moment-based estimators proposed in this paper with the estimator proposed by Yacoub. We also discuss the applicability of each estimator in different scenarios.

2. k-µ Fading Channel Models

In wireless communication systems, the received signal at the receiver can be represented as [27]
y ( t ) = x ( t ) h ( t ) + n ( t )
where x ( t ) is the transmit signal, h ( t ) is the channel impulse response (CIR) signal, and n ( t ) is the noise. The CIR signal describes the energy attenuation and time delay of countless electromagnetic waves that undergo optical phenomena such as reflection and diffraction in the channel, and then arrive at the receiver. Energy attenuation and time delay are two uncorrelated random variables that follow different probability distributions. Various channel fading models (such as Rice, Nakagami-M, etc.) use statistical methods to describe the energy attenuation random variable. Assuming that a known frequency-domain signal X ( f ) is transmitted from the transmitter and undergoes multipath fading, the received signal filtered at the receiver can be represented in the frequency domain as Y ( f ) = X ( f ) H ( f ) + N ( f ) . The estimated channel state information (CSI) obtained from it is H ^ ( f ) = X 1 ( f ) [ Y ( f ) N ( f ) ] , which is the channel frequency response (CFR). By performing the inverse Fourier transform on H ^ ( f ) , the channel impulse response (CIR) can be obtained and denoted as h ^ ( t ) , and r = | h ( t ) | .
The classical Rician channel model assumes that the obstacles in the channel follow a uniform distribution, and it defines the signal that arrives at the receiver as a cluster; there is no time correlation between each electromagnetic wave, and the energy fading of the electromagnetic wave obeys a normal distribution and the phase follows a uniform distribution. However, under practical scenario conditions, the electromagnetic waves that arrive at the receiver will be time-dependent because of the non-uniform distribution of the obstacles in the channel, and the CIR will show multiple energy peaks in the experimental data. According to the multiple-cluster effects of the transmitted electromagnetic waves arriving at the receiving end, the electromagnetic waves can be divided into a set of clusters; the number of clusters is denoted by μ, and each cluster is individually subject to the Rice distribution. This new method for modeling electromagnetic waves was first proposed by Yacoub [10] and is called the k-µ fading distribution model. The cumulative electrical frequency amplitude (envelope) of the µ clusters of electromagnetic waves that arrive at the receiver satisfies the following physical model [24].
r 2 = i = 1 μ ( p i + x i ) 2 + i = 1 μ ( q i + y i ) 2
where r 2 is the instantaneous power of the carrier signal at the receiving end, μ is the number of clusters of electromagnetic waves. pi, qi are the mean values of the homogeneous and quadrature components of the i-th cluster of multipath waves, respectively, and p i 2 , q i 2 are the real and imaginary parts of the i-th cluster of the deterministic component signal. xi, yi are the real and imaginary parts of the scattered wave signal in the i-th cluster signal, respectively, which are the mutually independent Gaussian processes and satisfy the relationship of E [ x i ] = E [ y i ] = 0 and E [ x i 2 ] = E [ y i 2 ] = σ 2 . Therefore, the probability distribution function (PDF) of the k-µ distribution can be expressed as [24]
f ( α ) = 2 μ ( 1 + k ) ( μ + 1 ) / 2 α μ k ( μ 1 ) / 2 e ( μ k + μ ( 1 + k ) α 2 ) I μ 1 ( 2 μ k ( 1 + k ) α )
where α is the normalized envelope and satisfies the relationship of α = r / E [ r 2 ] . According to the probability distribution conversion relationship in probability theory, the PDF of the k-µ fading distribution with respect to the envelope r can be expressed as:
f ( r ) = f ( α ) Ω = 2 μ ( 1 + k ) ( μ + 1 ) / 2 r μ Ω μ + 1 k ( μ 1 ) / 2 e ( μ k + μ ( 1 + k ) r Ω 2 ) I μ 1 ( 2 μ k ( 1 + k ) r Ω )
the parameter Ω represents the average signal power Ω = E [ r 2 ] . When the number of clusters µ is a natural number, k represents the ratio of the deterministic component to the non-line-of-sight path component (or scatter component). When k = 0 and µ = M, the k-µ fading distribution degenerates to the Nakagami-M distribution. When k = K and µ = 1, the k-µ fading distribution becomes the Rician distribution. When k = 0 and µ = 1, the k-µ fading distribution becomes the Rayleigh distribution. Applying E [ r n ] = r ¯ n E [ α n ] = E [ r 2 ] n E [ α n ] from ([19], Equation (6)) to the function from ([20], Equation (5)), the moment-generating function (MGF) for the k-µ distribution can be expressed as:
E [ r n ] = E ( r 2 ) n Γ ( μ + n / 2 ) exp ( k μ ) Γ ( μ ) [ ( 1 + k ) μ ] n / 2 F 1 1 ( μ + n / 2 , μ ; k μ )
Here, n denotes the order of the moment function, F 1 1 (   ·   ,   ·   ;   · ) represents the Kummer confluent hypergeometric function, and Γ ( ) represents the gamma function. In previous research, Yacoub has proposed a closed-form estimator for the parameter k in ([24], Equation (12)), and demonstrated the mapping relationship between k and μ in ([24], Equation (9)), as shown below:
k ^ y a c o u b = 2 ( E [ α 4 ] 1 ) 2 E [ α 4 ] 2 E [ α 4 ] E [ α 6 ] 2 1
μ ^ y a c o u b = 1 V a r [ α 2 ] 1 + 2 k ^ y a c o u b ( 1 + k ^ y a c o u b ) 2
In the above equation, Var[·] denotes the variance operator, and V a r [ α 2 ] =   V a r [ r 2 ] / E 2 [ r 2 ] . Based on the equation given in reference E [ r n ] = r ¯ n E [ α n ]   = E [ r 2 ] n E [ α n ] ([19], Equation (6)), when the total signal power is E [ r 2 ] = 1 , the variable α in Equations (6) and (7) can be exchanged with the envelope r .

3. The Derivation and Design of the k-μ Estimators

The moment estimation method has always been one of the most commonly used methods for extracting channel parameters, owing to its advantages of closed-form simplicity, low algorithm complexity, and small estimation error. It has been widely used in the research of fading channel model parameter estimators, such as the K-factor estimator based on the Rice distribution [28], the M-parameter estimator based on the Nakagami distribution [29], and the k-parameter estimator based on the k-µ distribution [25], among others. Moreover, the various-order moment values of the signal envelope can be easily obtained through experimental data. As the sample size approaches infinity, the theoretical n-th moment value and the expected value of the n-th power of the measured data have unbiased asymptotic properties. Therefore, in this paper, we choose the moment estimation method to estimate the k and µ parameters of the k-µ distribution.
From the moment functions of the k-μ distribution in Equation (5), it is evident that each order of moment functions contains gamma functions and Kummer confluent hypergeometric functions. The expansion of the right-hand side of the moment function equation is a complex and irreducible combination polynomial, making it difficult to solve for k and μ using inverse function methods or substitution-based elimination techniques. At this point, we observe that even-order moments share the same functional structure. By utilizing two even-order moment functions as the numerator and denominator, we can eliminate the gamma function terms and exp(-) terms in Equation (5), thus obtaining the simplified second-, fourth-, and sixth-order moments of the kappa-mu distribution (for the detailed derivation process, please refer to Appendix A).
E ( r 2 ) = Ω
E ( r 4 ) = Ω 2 1 + 1 + 2 k μ ( 1 + k ) 2
E ( r 6 ) = Ω 3 μ ( μ + 1 ) ( μ + 2 ) + 3 k μ ( μ + 1 ) ( μ + 2 ) + 3 ( k μ ) 2 ( μ + 2 ) + ( k μ ) 3 ( 1 + k ) 3 μ 3
Equations (8)–(10) form a channel parameter estimation model with respect to variables k and μ. In this model, Ω represents the total signal power, while E ( r 2 ) , E ( r 4 ) , E ( r 6 ) can be estimated from either empirical measurements or simulation data. Transforming Equation (9) into a function with respect to μ, we obtain:
μ = E 2 ( r 2 ) V a r [ r 2 ] ( 1 + k ) 2 ( 1 + 2 k ) = η ( 1 + k ) 2 ( 1 + 2 k )
where η = E 2 ( r 2 ) / V a r [ r 2 ] , Var[·] denotes the variance operator. By substituting Equation (11) into the sixth-order moment function given by Equation (10), we can simplify it to obtain Equation (12) as:
E ( r 6 ) E 3 ( r 2 ) = ( 2 + 6 k ) ( 1 + k ) η 2 ( 1 + 2 k ) 2 + 3 η + 1
Expanding Equation (12) allows us to obtain that
[ 4 ( E ( r 6 ) E 3 ( r 2 ) η 2 3 η η 2 ) 6 ] k 2 + [ 4 ( E ( r 6 ) E 3 ( r 2 ) η 2 3 η η 2 ) 8 ] k + ( E ( r 6 ) E 3 ( r 2 ) η 2 3 η η 2 ) 2 = 0 ,
Considering the initial conditions that k and μ must both be greater than zero, we can eliminate an invalid solution concerning k. Consequently, our first closed-form estimate for k and μ derived in this paper is as follows, denoted by a delta subscript ( Δ ):
k ^ Δ = Δ + Δ 2 2 Δ μ ^ Δ = η ( 1 + 2 k ^ Δ ) ( 1 + k ^ Δ ) 2 , η = E 2 ( r 2 ) V a r ( r 2 ) Δ = 4 2 ( E ( r 6 ) [ E ( r 2 ) ] 3 η 2 3 η η 2 ) ,
Simultaneously, we define the received power E ( r 2 ) = P 1 + P 2 in the k-µ distribution, where P1 represents the energy of the line of sight (LOS) path and P2 corresponds to the energy of the non-line-of-sight (NLOS) paths. Parameter k = 1 n ( p i 2 + q i 2 ) / 1 n ( x i 2 + y i 2 ) = P 1 / P 2 , the second-, fourth-, and sixth-order moments can be further transformed as follows (for a detailed derivation process, please refer to Appendix B):
E ( r 2 ) = U 2 = P 1 + P 2
E ( r 4 ) = U 4 = μ + 1 μ U 2 2 1 μ P 1 2
E ( r 6 ) = U 6 = ( μ + 1 ) ( μ + 2 ) μ 2 U 2 3 3 ( μ + 2 ) μ 2 P 1 2 U 2 + 4 μ 2 P 1 3
Equations (15)–(17) constitute a solution model with respect to the LOS path energy P1 and μ. And U2, U4, and U6 can be obtained from either measurement data or simulation data.
Let E [ r 2 ] = U 2 , E [ r 4 ] = U 4 , E [ r 6 ] = U 6 , and transform Equation (16) into the following form of function:
μ = U 2 2 P 1 2 U 4 U 2 2
Substituting this function into Equation (17) allows us to eliminate μ, resulting in an equation concerning P1:
0 = [ U 2 3 3 U 2 ( U 2 2 U 4 ) U 6 ] P 1 4 + 4 [ U 4 2 2 U 4 U 2 2 + U 2 4 ] P 1 3 + [ U 2 3 ( U 2 2 3 U 4 ) 3 U 2 ( 2 U 4 2 + U 2 4 3 U 4 U 2 2 ) + 2 U 6 U 2 2 ] P 1 2 + U 2 3 ( 2 U 4 2 U 4 U 2 2 ) U 6 U 2 4 ,
By solving the aforementioned equation, we can obtain the second closed-form estimate for k and µ derived in this paper, as presented below:
k ^ A = P ^ 1 U 2 P ^ 1 μ ^ A = U 2 2 P ^ 1 2 U 4 U 2 2 , P ^ 1 = b + 2 A E B 2 B A 4 a
Subsequently, we can derive the third closed-form estimate for k and µ in this paper, as shown below:
k ^ B = P ^ 1 U 2 P ^ 1 μ ^ B = U 2 2 P ^ 1 2 U 4 U 2 2 , P ^ 1 = b + D 4 a  
The parameters in the model satisfy the following relationship:
a = [ U 2 3 3 U 2 ( U 2 2 U 4 ) U 6 ] , b = 4 [ U 4 2 2 U 4 U 2 2 + U 2 4 ] , c = [ U 2 3 ( U 2 2 3 U 4 ) 3 U 2 ( 2 U 4 2 + U 2 4 3 U 4 U 2 2 ) + 2 U 6 U 2 2 ] , d = [ U 2 3 ( 2 U 4 2 U 4 U 2 2 ) U 6 U 2 4 ] , A = D 2 3 F , B = D F 9 E 2 , D = 3 b 2 8 a c , E = b 3 + 4 a b c   , F = 3 b 4 + 16 a 2 c 2 16 a b 2 c 64 a 3 d ,
and we summarize the estimation process of each estimator as below:
In Figure 1, the input to the estimator is the channel envelope data r, which consists of N sample points. At a certain moment, the instantaneous envelope r is the superposition of the useful signal and noise, which are temporarily indistinguishable. Consequently, in practical applications, it is necessary to filter the received signals. However, our estimator is also applicable to the study of noisy environments. In cases where the input signal is represented by instantaneous signal-to-noise ratios (SNRs) γ , the estimator proposed in this study remains equally valid and applicable.
In signal measurement experiments conducted in practical scenarios, the transmitter sends a pseudorandom sequence, and the channel information is obtained by correlation operating the received data with the known transmitted pseudorandom sequence. The calculation process for the CFR is denoted as H ^ ( f ) = Y ( f ) c o n j ( X ( f ) ) / | X ( f ) | 2 , where X ( f ) represents the transmitted frequency domain signal (PN sequence), Y ( f ) denotes the down-converted IQ signal of the received frequency domain signal, and c o n j ( ) signifies the complex conjugate. At this point, H ^ ( f ) encompasses the transmitter system response H T ( f ) , free-space channel response H C ( f ) , and receiver system response H R ( f ) . The noise of the transmitter system is contained within H T ( f ) , the noise in the propagation medium is contained within H C ( f ) , and the noise of the receiver system is contained within H R ( f ) . The channel system can be represented as Y ( f ) = H T ( f ) H C ( f ) H R ( f ) X ( f ) . To further calibrate the system, the transmitter and receiver are directly connected via a high-frequency cable to obtain the system response, where the received frequency domain signal is Y s y s ( f ) = H T ( f ) H R ( f ) X ( f ) . Consequently, the calibrated transmitter and receiver channel frequency response is denoted as H ^ C ( f ) = Y ( f ) / Y y s y ( f ) . Substituting the calculation process of the CFR, the calibrated transmitter and receiver channel frequency response is represented as H ^ C ( f ) = H ^ ( f ) X ( f ) / Y y s y ( f ) . Performing the inverse Fourier transform on H ^ C ( f ) results in the calibrated channel impulse response h ^ ( t ) , which is a discrete variable. By taking the magnitude of h ^ ( t ) , the discrete envelope signal r = [ r 1 , r 2 , r 3 , , r n ] can be obtained. By calculating the second-, fourth-, and sixth-order moments of the envelope signal and incorporating them into the estimation process depicted in Figure 1, the k and μ parameters can be estimated.

4. Performance Analysis of Closed-Form Estimator

4.1. Computational Complexity Analysis of Closed-Form Estimator

Algorithmic complexity is a criterion for evaluating the execution efficiency of an effective algorithm. When the sample size is N, the total number of basic operations performed by different algorithms is denoted as time complexity, while the required memory space is denoted as space complexity. Necessary procedures that both types of complexities need to address are not considered in the comparison of complexities. We present the complexity analysis of four closed-form solutions in Table 2:
In Table 2, O 1 ( ) represents the number of basic operations required for computing other variables (excluding k and μ) before the closed-form calculation. For example, the O 1 ( ) value for Estimation Method k ^ - μ ^ Δ indicates the complexity of solving for variable V a r , η , Δ ; similarly, the O 1 ( ) value for Estimation Method k ^ - μ ^ B also represents the complexity of solving for variable a, b, c, D, P1. O 2 ( ) denotes the complexity of the k and μ estimation functions, such as the closed-form expressions for Method k ^ - μ ^ Y a c o u b given by Equations (6) and (7). O 3 ( ) refers to the memory space required to store temporary variables during the estimation process. For instance, Estimation Method k ^ - μ ^ Y a c o u b only needs to store the variance operator Var[·], thus O 3 ( ) = 1 . Since each moment estimation requires the calculation of the second-order, fourth-order, and sixth-order moments of the envelope, the computational cost of E ( r 2 ) , E ( r 4 ) , E ( r 6 ) is not included in the complexity analysis.
To further analyze time complexity, we employ the Monte Carlo method to generate channel data following the kappa-mu distribution and compare estimation efficiency based on the time taken for the estimation. As shown in Figure 1, we use the channel impulse response h(t) as the data sample. The channel data generated using the Monte Carlo method in this section does not include noise.
We set the sample size to 10,000, with each round of estimation experiments using the same data. The efficiency of the estimator is evaluated based on the total runtime of multiple estimations, and the runtime is recorded in Table 3.
As can be seen from Table 2 and Table 3, Estimation Method k ^ - μ ^ Y a c o u b has the optimal complexity, followed by Method k ^ - μ ^ Δ , while Method k ^ - μ ^ A is the least efficient. However, in the authors’ numerous simulation studies, it was found that when the estimation error of Method E ( r 2 ) , E ( r 4 ) , E ( r 6 ) is relatively large, the estimator k ^ - μ ^ Y a c o u b may produce k and μ values that are less than zero or complex numbers, which can be considered erroneous estimates. In practical applications, repeated estimations may generally consume more time.

4.2. The Simulation for Estimator Accuracy Analysis

Due to the probability density function (PDF) of the κ-μ fading channel model being controlled by the model parameters κ and μ, existing research has not yet been able to explain the conversion relationship between these two parameters. Consequently, it is not possible to evaluate the accuracy of the estimation algorithm using the deviation distance ( k , μ ) ( k ^ , μ ^ ) 2 between the theoretical and estimated values. Therefore, this paper employs the error f P D F ( k , μ ) f P D F ( k ^ , μ ^ ) 2 between the theoretical PDF curve and the estimated PDF curve, as well as the root mean square error (RMSE), to evaluate the estimation accuracy of the algorithm. The RMSE reflects the degree to which a set of data deviates from the true data. As the RMSE approaches 0, it indicates a smaller deviation between the theoretical PDF curve and the estimated PDF curve, resulting in higher estimation accuracy. The mathematical expression for RMSE is given as 1 n f P D F ( k , μ , r i ) f P D F ( k ^ , μ ^ , r i ) 2 / n , where r i r n represents the minimum and maximum values of the envelope. In the experiments of this paper, the gradient Δ r = r i + 1 r i is 0.01. In the equation, f P D F ( k , μ , r i ) is the theoretical probability (assumed to be P i ) when the envelope is equal to r i , and f P D F ( k ^ , μ ^ , r i ) is the estimated probability (assumed to be P ^ i ) when the envelope is equal to r i . The calculation method of RMSE can be expanded as ( P 1 P ^ 1 ) 2 + ( P 2 P ^ 2 ) 2 + + ( P i P ^ i ) 2 + + ( P n P ^ n ) 2 / n   0.5 .
Simulation experiment samples were generated using the Monte Carlo method, producing channel data that follows the κ-μ distribution. In each estimation experiment, four estimators employed the same set of generated channel data, with a sample size of 10,000. Four simulation experiments were conducted with the following parameter settings: κ = 2 and μ = 1; κ = 4 and μ = 1; κ = 4 and μ = 2; and κ = 6 and μ = 2. The estimated PDF curves were plotted based on the estimation results and compared with the theoretical PDF curves. The comparison and error plots are illustrated in Figure 2. The estimated values and root mean square errors were rounded to the 12th decimal place and recorded in Table 4.
In Figure 2, the left column of the plots demonstrates the comparison between the PDF curves calculated from the estimated k and μ parameters and the theoretical PDF curves. The right column of plots shows the errors between the four estimated PDF curves and the theoretical PDF curves. In Table 3, RMSEpdf refers to the root mean square error between the PDF curve and the theoretical PDF curve, while the last column displays the ranking of estimation accuracy for each estimator in each experiment.
Throughout numerous experiments conducted during the research process, we discovered that the k, μ, and RMSEpdf values for estimators k ^ - μ ^ Δ and k ^ - μ ^ A were nearly identical, with differences only appearing beyond the seventh decimal place. Consequently, we consider them to have the same estimation accuracy.
Meanwhile, in multiple experiments, Estimator k ^ - μ ^ Y a c o u b frequently exhibits instances of negative or complex values for k and µ, which we consider to be incorrect estimations. The cause of this phenomenon is the lack of stability in Estimator k ^ - μ ^ Y a c o u b , and there are three factors contributing to its instability:
Firstly, in Equation (6) k ^ y a c o u b , the denominator of the term 2 E [ α 4 ] 2 E [ α 4 ] E [ α 6 ] often becomes negative, resulting in complex-valued estimates.
Secondly, Estimator k ^ y a c o u b violates two principles of numerical computation that should be adhered to. Ideally, the subtraction of two numbers with close values should be avoided. However, in the denominator of term k ^ y a c o u b in Equation (6), the values of E [ α 6 ] are very close to that of 2 E [ α 4 ] 2 E [ α 4 ] .
Last, it is preferable to avoid dividing a large number by a small one. Nonetheless, in Equation (6) k ^ y a c o u b , the numerator 2 ( E [ α 4 ] 1 ) is several orders of magnitude larger than the denominator 2 E [ α 4 ] 2 E [ α 4 ] E [ α 6 ] .
Therefore, under the experimental scenarios presented in Table 4, Estimator k ^ - μ ^ Δ and k ^ - μ ^ A outperforms other estimators. Moreover, the computation time of Estimator k ^ - μ ^ Δ is smaller than that of Estimator k ^ - μ ^ A , making Estimator k ^ - μ ^ Δ the most optimal choice.
Additionally, we conducted another set of comparative experiments. Assuming k = 0 and µ taking values of 1, 1.5, and 2, we generated channel data using the Monte Carlo method. In each estimation experiment, the four estimators used the same data samples, with a sample size of 10,000. The estimated PDF curves, theoretical PDF curves, and their errors f P D F ( k , μ ) f P D F ( k ^ , μ ^ ) 2 are depicted in Figure 3. Furthermore, the estimation results and root mean square errors (RMSE) are recorded in Table 5.
In Figure 3, the experiments for different µ values under the condition of k = 0 are presented. The left image displays a comparison between the estimated PDF curves and the theoretical PDF curves, while the right image illustrates the errors between the estimated PDF curves and the theoretical PDF curves. In Table 5, the RMSEpdf column shows the root mean square errors between the PDF curves and the theoretical PDF curves; values closer to 0 indicate higher estimation accuracy. The last column demonstrates the ranking of estimation accuracy for each estimator in every experiment.
In Table 5, we still observe that the k, µ, and RMSEpdf values of Estimator k ^ - μ ^ Y a c o u b and k ^ - μ ^ Y a c o u b are almost identical, with differences only occurring after the seventh decimal place. Consequently, we can assert that these two estimators have the same estimation accuracy, regardless of the variations in k and µ.
At the same time, Estimator k ^ - μ ^ Y a c o u b continues to exhibit incorrect estimates with negative or complex values for k and µ. In cases of correct estimation, as shown in Table 5, the RMSEpdf of the estimator is sufficiently small.
From the last column in Table 5, it can be observed that under the condition of k = 0, Estimator k ^ - μ ^ B has the best estimation accuracy, while Estimator k ^ - μ ^ Y a c o u b has the worst. Therefore, Estimator k ^ - μ ^ Y a c o u b is more suitable for situations with extremely poor channel quality, where there are numerous obstacles between the transmitting and receiving antennas and no line-of-sight component (k approaching 0).
In summary, through complexity analysis and simulation, we find that Estimator k ^ - μ ^ Y a c o u b has the lowest estimation efficiency. When k approaches 0, using Estimator k ^ - μ ^ B is the best choice; in other cases, Estimator k ^ - μ ^ Δ is the most optimal option.

4.3. Impact of SNR on Estimator Performance Analysis

Noise can cause an increase in the bias and variance of the estimators, thereby reducing the accuracy and efficiency of the estimation. In order to analyze the impact of noise on the accuracy of the estimators, we assume signal-to-noise ratios (SNRs) of 30 dB, 25 dB, 20 dB, and 15 dB and conduct simulation analysis under the assumption of channel parameters k = 3.5 and μ = 1.5. In this case, the envelope signal satisfies the relationship r a d d N o i s e = r k μ S i g n a l + r n o i s e , r k μ S i g n a l is generated by the Monte Carlo method, and r n o i s e represents additive Gaussian white noise. For ease of comparison, the analysis results here are presented with a precision of four decimal places. In each experiment, the four estimators use the same sample data, and the sample size is 10,000.
In multiple simulation experiments, we find that when the SNR is 30 dB, the estimation method k ^ - μ ^ Y a c o u b has about a 40% probability of not producing valid results; when the SNR is 25 dB, the estimation method k ^ - μ ^ Y a c o u b has about a 60% probability of not producing valid results; when the SNR is 20 dB, the estimation method k ^ - μ ^ Y a c o u b has about a 90% probability of not producing valid results; and when the SNR is 15 dB, the estimation method k ^ - μ ^ Y a c o u b can no longer produce valid results. In all the above experiments, the channel parameter estimator we proposed demonstrated significantly lower error probabilities under the aforementioned signal-to-noise ratio levels, indicating its robustness in various noise conditions. In the figures below, we present the valid estimation results.
From Table 6, we can see that the parameters k and µ have already become less than 0, and the k ^ - μ ^ Y a c o u b estimation method can no longer produce valid results. Moreover, in Figure 4d, it can be observed that the curve obtained from the k ^ - μ ^ Y a c o u b estimation method has been completely deformed. This indicates that the k ^ - μ ^ Y a c o u b method will fail when the noise interference is significant.
From the RMSEpdf column in Table 6, it can be seen that as the SNR decreases, the estimation errors of the estimators gradually increase. When the SNR drops to 15 dB, the estimation errors of all estimators are essentially unacceptable. Therefore, we can draw the conclusion that our proposed estimators have limitations in their performance when the SNR is low. When using the estimators proposed in this paper as well as Yacoub’s estimator, it is necessary to first filter out the noise or increase the signal gain. Furthermore, as the SNR decreases, the RMSE values of all estimators gradually increase. At the 15 dB level, k ^ - μ ^ Y a c o u b estimator cannot be computed, and with its higher estimation error probability, it becomes unsuitable for practical applications. Moreover, although the RMSE of the k ^ - μ ^ B estimator is higher than that of the k ^ - μ ^ Δ / k ^ - μ ^ A estimator under higher SNR conditions, as the SNR decreases, the performance of the k ^ - μ ^ B estimator declines more slowly than that of the k ^ - μ ^ Δ / k ^ - μ ^ A estimator. Therefore, at 15 dB, the RMSE performance of the k ^ - μ ^ B estimator is better than that of the k ^ - μ ^ Δ / k ^ - μ ^ A estimator, which intuitively reflects the correct trend of estimation accuracy changes for different channel parameter estimators under different SNR conditions. In the methods proposed in this paper, the k ^ - μ ^ Δ / k ^ - μ ^ A estimator is used for communication scenarios with better channel conditions, while the k ^ - μ ^ B estimator is suitable for non-line-of-sight (NLOS) scenarios with poorer channel conditions and lower received SNR. This is consistent with the conclusions reached.

4.4. Impact of Pilot Signal on Estimator Performance Analysis

In general, the amount of information required by the pilot signal depends on the specific estimation methods and the channel environment. The impact of the pilot signal on the estimation performance also varies with different estimation algorithms and channel conditions. Basically, a higher pilot signal power or density can lead to better channel estimation accuracy, but also at the cost of higher overhead and decreased spectral efficiency. On the other hand, a lower pilot signal power or density can result in lower accuracy and higher estimation error, especially in fast-varying channels.
Consequently, it is indeed true that with a longer pilot signal, more samples can be obtained, and as the sample size increases, the estimations of k and μ parameters will approach unbiasedness and become more accurate. In practical applications, extending the pilot signal duration and increasing the sampling rate can enhance estimation accuracy. However, this comes at the cost of consuming more communication resources, which in turn results in a reduced data rate for wireless communication. It is essential to weigh these trade-offs carefully.
To address this question, we conducted related parameter measurements and fitting simulation experiments in an indoor environment. The wireless communication scenario parameters are set as follows: The operating frequency band is the C-band (5 GHz); assuming a Rician channel with a channel length of 20, the transmitter employs a pilot PN sequence of 2048 symbols (25% of the total amount of transmitted data), BPSK modulation, and a code rate of 20 Mbps transmission scheme. On the receiver. At the receiving end, we used a receiving demodulation and processing device equipped with high-speed sampling, fiber-optic broadband transmission, and high-speed storage, achieving a transmission rate of 5 Gbps. The sampling frequency can reach up to 500 Msps.
The left column in Figure 5 shows the time on the x-axis and signal power on the y-axis (within a one-second acquisition period, we detected 10 instances of robust signal waves; however, only a single representative example is displayed here). On the right column, we estimate the received signal, where the gray area represents the frequency histogram of envelope r, and the two curves are estimation curves, with the x-axis representing the normalized envelope and the y-axis representing the probability. In the right column, using the estimator proposed in this paper and defining µ = 1, we can obtain the estimated curve—“Rice distribution” (by taking advantage of the k-µ distribution transforming into the Rice distribution when µ = 1). Using the estimator without restricting the values of k and μ, we can obtain the curve—“k-μ distribution”. As can be seen from the graphs, when the sample size is sufficient, the k-µ distribution approximates the Rice distribution. When the sample size is not enough, the k-µ distribution estimator exhibits excellent adaptability.
As Figure 5 indicates, as the number of sampling points increases, the k-µ distribution becomes more closely aligned with the actual probability distribution, and the greater the number of samples, the more evident the alignment. In this experimental scenario, when the number of samples N obtained from pilot signal sampling is greater than 8000, the PDF estimation curve generated by the sample points essentially matches the histogram and theoretical curve. Consequently, a consistent correlation between the theoretical analysis and experimental measurement results concerning the impact of the pilot signal on estimator performance is observed.

5. Conclusions

The paper introduces a novel and efficient methodology for estimating multiple parameters in the k-µ fading channel model. Two sets of solutions are derived using moment-generating functions based on second-, fourth-, and sixth-order moments. Three sets of closed-form solutions are designed for estimating the k and µ parameters based on the previously derived solutions. Performance analysis of the proposed parameter estimators using simulation results demonstrates strong consistency between theoretical and estimated values. The proposed parameter estimators address some of the shortcomings of existing estimators through simulation and comparison, including large estimation errors, susceptibility to errors, and fragility to noise interference. The estimators can calculate k and µ parameters based on the envelope and SNR signals. The work presents a mathematically tractable framework for estimating characteristic parameters of wireless channels in small-scale environments. It is expected to allow for real-time and accurate analysis of complex scenarios, such as moving targets, time-varying conditions, or non-line-of-sight (NLOS) extreme communication, based on k-µ fading channel model parameters. Simulation results demonstrate that the proposed parameter estimator k ^ - μ ^ Δ is appropriate for low latency and limited resource overhead communication scenarios, and the other estimator k ^ - μ ^ B is suitable for scenarios with poor channel quality, numerous obstacles between antennas, or no LOS communication path available (k approaches zero, NLOS communication scenarios).

Author Contributions

Conceptualization, Z.F. and J.T.; formal analysis, Z.J. and D.H.; investigation, Z.J., X.L., P.F. and D.H.; methodology, Z.F. and J.T.; project administration, J.T.; supervision, X.L. and P.F.; validation, Z.F., Z.J. and J.T.; visualization, Z.J., P.F. and D.H.; writing—original draft, Z.F., Z.J. and J.T.; writing—review and editing, Z.F., X.L., P.F., D.H. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the CAEP Foundation under Grant CX20200010, National Natural Science Foundation of China under grant number 61871084, and National Natural Science Foundation of China under grant number 62271109.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The transformation form of confluent hypergeometric function given in [30]:
F 1 1 ( α , ϒ ; z ) = e z F 1 1 ( ϒ α , ϒ ; z )
The moment function of the k-µ distribution can be further expressed as:
E [ r n ] = E ( r 2 ) n Γ ( μ + n / 2 ) exp ( k μ ) Γ ( μ ) [ ( 1 + k ) μ ] n / 2 F 1 1 ( μ + n / 2 , μ ; k μ )
where n is the order of the moment function, Γ ( ) is the gamma function, and Γ ( n + 1 ) = n ! . F 1 1 ( ,   ;   ) is the Kummer confluent hypergeometric function, with the expressions ([30], Section 9.212, Equation (1)):
F 1 1 ( α , β ; x ) = n = 0 ( α ) n x n ( β ) n n ! , ( α ) n = 1 ( α ) ( α + 1 ) ( α + n 1 )
Substituting the gamma function and the Kummer confluent hypergeometric function into Equation (A2), setting n to 2 and 4, and eliminating the gamma function term using the second-order moments compared to the fourth-order distance:
E ( r 2 ) E ( r 4 ) = E ( r 2 ) Γ ( μ + 1 ) Γ ( μ ) [ ( 1 + k ) μ ] 1 F 1 1 ( 1 , μ ; k μ ) E ( r 2 ) 2 Γ ( μ + 2 ) Γ ( μ ) [ ( 1 + k ) μ ] 2 F 1 1 ( 2 , μ ; k μ ) = μ ! [ ( 1 + k ) μ ] E ( r 2 ) ( μ + 1 ) ! F 1 1 ( 1 , μ ; k μ ) F 1 1 ( 2 , μ ; k μ ) = ( 1 + k ) μ E ( r 2 ) ( μ + 1 ) [ 1 + ( 1 ) μ ( k μ ) 1 ! + ( 1 ) ( 0 ) μ ( μ + 1 ) ( k μ ) 2 2 ! + ( 1 ) ( 0 ) ( 1 ) μ ( μ + 1 ) ( μ + 2 ) ( k μ ) 3 3 ! + ( 1 ) ( 0 ) ( 1 ) ( 2 ) μ ( μ + 1 ) ( μ + 2 ) ( μ + 3 ) ( k μ ) 4 4 ! + ] [ 1 + ( 2 ) μ ( k μ ) 1 ! + ( 2 ) ( 1 ) μ ( μ + 1 ) ( k μ ) 2 2 ! + ( 2 ) ( 1 ) ( 0 ) μ ( μ + 1 ) ( μ + 2 ) ( k μ ) 3 3 ! + ( 2 ) ( 1 ) ( 0 ) ( 1 ) μ ( μ + 1 ) ( μ + 2 ) ( μ + 3 ) ( k μ ) 4 4 ! + ] = ( 1 + k ) μ E ( r 2 ) ( μ + 1 ) ( 1 + k ) ( μ + 1 ) + 2 k ( μ + 1 ) + k 2 μ ( μ + 1 ) = 1 E ( r 2 ) μ ( 1 + k ) 2 μ ( 1 + k ) 2 + 1 + 2 k
After organizing the above equations, the expression for the fourth-order moment is given by
E ( r 4 ) = E ( r 2 ) 2 1 + 1 + 2 k μ ( 1 + k ) 2
Equivalently, we then set n to be 2 and 6, and using the second-order moments compared with the sixth-order distance to eliminate the gamma function term, we can obtain the ratio relationship between the second-order moments and the sixth-order moments as follows:
E ( r 2 ) E ( r 6 ) = E ( r 2 ) Γ ( μ + 1 ) Γ ( μ ) [ ( 1 + k ) μ ] 1 F 1 1 ( 1 , μ ; k μ ) E ( r 2 ) 3 Γ ( μ + 3 ) Γ ( μ ) [ ( 1 + k ) μ ] 3 F 1 1 ( 3 , μ ; k μ ) = μ ! [ ( 1 + k ) μ ] 2 E ( r 2 ) 2 ( μ + 2 ) ! F 1 1 ( 1 , μ ; k μ ) F 1 1 ( 3 , μ ; k μ ) = [ ( 1 + k ) μ ] 2 E ( r 2 ) 2 ( μ + 1 ) ( μ + 2 ) [ 1 + ( 1 ) μ ( k μ ) 1 ! ] [ 1 + ( 3 ) μ ( k μ ) 1 ! + ( 3 ) ( 2 ) μ ( μ + 1 ) ( k μ ) 2 2 ! + ( 3 ) ( 2 ) ( 1 ) μ ( μ + 1 ) ( μ + 2 ) ( k μ ) 3 3 ! ] = 1 E ( r 2 ) 2 μ 3 ( 1 + k ) 3 μ ( μ + 1 ) ( μ + 2 ) + 3 k μ ( μ + 1 ) ( μ + 2 ) + 3 ( k μ ) 2 ( μ + 2 ) + ( k μ ) 3
After the organization of Equation (A6), the expression for the sixth-order moments can be equivalent to
E ( r 6 ) = E ( r 2 ) 3 μ ( μ + 1 ) ( μ + 2 ) + 3 k μ ( μ + 1 ) ( μ + 2 ) + 3 ( k μ ) 2 ( μ + 2 ) + ( k μ ) 3 ( 1 + k ) 3 μ 3

Appendix B

The total LOS path energy of the k-µ distribution is now defined as P1 and the total NLOS path energy as P2. In the physical model of the k-µ distribution, we have the relationship 1 n ( p i 2 + q i 2 ) / 1 n ( x i 2 + y i 2 ) = P 1 / P 2 , the total received power E ( r 2 ) = P 1 + P 2 . As a result, the fourth-order moment function of the k-µ distribution (Equation (A5)) can be reformulated as:
E ( r 4 ) = E ( r 2 ) 2 1 + 1 + 2 k μ ( 1 + k ) 2 = ( P 1 + P 2 ) 2 μ ( P 1 + P 2 P 2 ) 2 + 1 + 2 P 1 P 2 μ ( P 1 + P 2 P 2 ) 2 = P 2 2 μ ( P 1 + P 2 P 2 ) 2 + P 2 2 1 + P 2 2 2 P 1 P 2 μ = ( P 1 + P 2 ) 2 + 1 μ ( P 2 2 + 2 P 1 P 2 + P 1 2 ) 1 μ P 1 2 = μ + 1 μ U 2 2 1 μ P 1 2
Similarly, Equation (A7) in Appendix A, that is, the fourth-order moment function of the k-µ distribution, can be reformulated as:
E ( r 6 ) = E ( r 2 ) 3 μ ( μ + 1 ) ( μ + 2 ) + 3 k μ ( μ + 1 ) ( μ + 2 ) + 3 ( k μ ) 2 ( μ + 2 ) + ( k μ ) 3 ( 1 + k ) 3 μ 3 = ( P 1 + P 2 ) 3 μ ( μ + 1 ) ( μ + 2 ) + 3 μ ( μ + 1 ) ( μ + 2 ) P 1 P 2 + 3 μ 2 ( μ + 2 ) P 1 P 2 2 + μ 3 P 1 P 2 3 μ 3 ( P 1 + P 2 ) 3 P 2 3 = P 2 3 μ ( μ + 1 ) ( μ + 2 ) + 3 μ ( μ + 1 ) ( μ + 2 ) P 1 P 2 + 3 μ 2 ( μ + 2 ) P 1 P 2 2 + μ 3 P 1 P 2 3 μ 3 = 3 μ + 2 μ 2 P 2 3 + 3 3 μ + 2 μ 2 P 1 P 2 2 + 6 μ P 1 2 P 2 + ( P 1 + P 2 ) 3 = ( μ + 1 ) ( μ + 2 ) μ 2 U 2 3 3 ( μ + 2 ) μ 2 P 1 2 U 2 + 4 μ 2 P 1 3

References

  1. Domanovitz, E.; Erez, U. Simple Bounds for the Symmetric Capacity of the Rayleigh Fading Multiple Access Channel. IEEE Trans. Wirel. Commun. 2020, 19, 4594–4603. [Google Scholar] [CrossRef]
  2. Yin, H.; Guo, X.; Liu, P.; Hei, X.; Gao, Y.J.A. Predicting Channel Quality Indicators for 5G Downlink Scheduling in a Deep Learning Approach. arXiv 2020, arXiv:2008.01000. [Google Scholar]
  3. Hanzo, L. Quadrature Amplitude Modulation: From Basics to Adaptive Trellis-Coded, Turbo-Equalised and Space-Time Coded OFDM, CDMA and MC-CDMA Systems, 2nd ed.; John Wiley & Sons: Chichester, UK, 2015. [Google Scholar]
  4. Moreno-Pozas, L.; Lopez-Martinez, F.J.; Paris, J.F.; Martos-Naya, E. The κµ Shadowed Fading Model: Unifying the κ-µ and η-µ Distributions. IEEE Trans. Veh. Technol. 2016, 65, 9630–9641. [Google Scholar] [CrossRef]
  5. Wang, X.-B.; Yu, Z.-W.; Hu, J.-Z. Entanglement-distribution maximization over one-sided Gaussian noisy channels. Phys. Rev. A 2010, 82, 022316. [Google Scholar] [CrossRef]
  6. Liao, R.-F.; Wen, H.; Wu, J.; Song, H.; Pan, F.; Dong, L. The Rayleigh Fading Channel Prediction via Deep Learning. Wirel. Commun. Mob. Comput. 2018, 2018, 6497340. [Google Scholar] [CrossRef]
  7. Jin, X.Y.; Zhang, Y.; Pan, C. A New Simulation Model of Rician Fading Channel. DEStech Trans. Eng. Technol. Res. 2016. [Google Scholar] [CrossRef]
  8. Shi, L.; Guo, Z.; Liu, Y.; Li, X. A novel simulation method for Nakagami-m fading channel. In Proceedings of the 2014 XXXIth URSI General Assembly and Scientific Symposium (URSI GASS), Beijing, China, 16–23 August 2014; pp. 1–4. [Google Scholar]
  9. Mazid-Ul-Haque, M.; Sohidul, S. Data Rate Limit in Low and High SNR Regime for Nakagami-q Fading Wireless Channel. Int. J. Adv. Comput. Sci. Appl. 2020, 11, 2020. [Google Scholar] [CrossRef]
  10. Yacoub, M.D. The κ-μ Distribution: A General Fading Distribution. In Proceedings of the Vehicular Technology Conference, Atlantic City, NJ, USA, 7–11 October 2001; Volume 3, pp. 1427–1431. [Google Scholar]
  11. Sarkar, S.; Adhikary, S.; Rahman, M.M. Optimizing channel capacity in a MIMO–PSK wireless rayleigh fading channel. Microsyst. Technol. 2019, 25, 1957–1963. [Google Scholar] [CrossRef]
  12. Kim, B.; Kim, J.; Chae, H.; Yoon, D.; Choi, J. Deep Neural Network-Based Automatic Modulation Classification Technique. In Proceedings of the 2016 International Conference on Information and Communication Technology Convergence (ICTC), Jeju, Republic of Korea, 19–21 October 2016; pp. 579–582. [Google Scholar] [CrossRef]
  13. Qian, B.; Xie, S.; Huang, J.; Feng, Y. Performance Analysis of OFDM System under Rice Channel. J. Phys. Conf. Ser. 2019, 1237, 032010. [Google Scholar] [CrossRef]
  14. Yang, F.; Zeng, X.; Mao, H.; Jian, X.; Du, D.; Tan, X.; Gao, Y. Nondata-aided error vector magnitude performance analysis over κ−μ shadowed fading channe. J. Commun. 2018, 39, 177–188. [Google Scholar] [CrossRef]
  15. Bilim, M. QAM signaling over κ-µ shadowed fading channels. Phys. Commun. 2019, 34, 261–271. [Google Scholar] [CrossRef]
  16. Bilim, M. A Performance Study on Diversity Receivers over κ-µ Shadowed Fading Channels. AEU Int. J. Electron. Commun. 2019, 112, 152934. [Google Scholar] [CrossRef]
  17. Kong, J.; Zeng, H.; Zhang, L.; She, Q.; Luo, Y.; Qin, H. Outage Performance Analysis of Cognitive Two-way Relay Network in κ-µ Fading Channel. Radio Commun. Technol. 2021, 47, 243–248. [Google Scholar] [CrossRef]
  18. Zhang, J.; Chen, X.; Peppas, K.; Li, X.; Liu, Y. On High-Order Capacity Statistics of Spectrum Aggregation Systems over κ-µ and κ-µ shadowed Fading Channels. IEEE Trans. Commun. 2016, 65, 935–944. [Google Scholar] [CrossRef]
  19. Kumar, S.; Chouhan, S. Outage probability analysis of cognitive decode-and-forward relay networks over κ-µ Shadowed Channels. In Proceedings of the 2015 21st Asia-Pacific Conference on Communications (APCC), Kyoto, Japan, 14–16 October 2015; pp. 433–437. [Google Scholar]
  20. Sofotasios, P.C.; Rebeiz, E.; Zhang, L.; Tsiftsis, T.A.; Cabric, D.; Freear, S. Energy Detection Based Spectrum Sensing Over κ-µ and κ-µ Extreme Fading Channels. IEEE Trans. Veh. Technol. 2013, 62, 1031–1040. [Google Scholar] [CrossRef]
  21. Gvozdarev, A.S. Average quality description of energy detection for generalised multipath shadowed κ-µ fading channels. J. Electromagn. Waves Appl. 2021, 35, 2163–2175. [Google Scholar] [CrossRef]
  22. Cotton, S. Human Body Shadowing in Cellular Device-to-Device Communications: Channel Modeling Using the Shadowed Fading Model. Sel. Areas Commun. IEEE J. 2015, 33, 111–119. [Google Scholar] [CrossRef]
  23. Kumar, S. Approximate Outage Probability and Capacity for κ-μ Shadowed Fading. IEEE Wirel. Commun. Lett. 2015, 4, 301–304. [Google Scholar] [CrossRef]
  24. Yacoub, M.D. The κ-μ distribution and the η-μ distribution. IEEE Antennas Propag. Mag. 2007, 49, 68–81. [Google Scholar] [CrossRef]
  25. Ribeiro, A.; Conforti, E. Asymptotically Efficient Moment-Based Estimator of the Parameter for the Distribution. Antennas Wirel. Propag. Lett. IEEE 2015, 14, 598–601. [Google Scholar] [CrossRef]
  26. Batista, F.P.; Souza, R.A.A.d.; Ribeiro, A.M.O. Maximum likelihood estimator for the α-η-μ fading environment. In Proceedings of the 2016 IEEE Radio and Wireless Symposium (RWS), Austin, TX, USA, 24–27 January 2016; pp. 133–136. [Google Scholar]
  27. Wang, J.; Cui, Y.; Sun, H.; Feng, X.; Han, G.; Wen, M.; Li, J.; Esmaiel, H. K-Factor Estimation for Wireless Communications Over Rician Frequency-Flat Fading Channels. IEEE Wirel. Commun. Lett. 2021, 10, 2037–2040. [Google Scholar] [CrossRef]
  28. Tang, P.; Zhang, J.; Molisch, A.F.; Smith, P.J.; Shafi, M.; Tian, L. Estimation of the K-Factor for Temporal Fading From Single-Snapshot Wideband Measurements. IEEE Trans. Veh. Technol. 2019, 68, 49–63. [Google Scholar] [CrossRef]
  29. Giunta, G.; Hao, C.; Orlando, D. Estimation of Rician K-Factor in the Presence of Nakagami-m Shadowing for the LoS Component. IEEE Wirel. Commun. Lett. 2018, 7, 550–553. [Google Scholar] [CrossRef]
  30. Zwillinger, D.; Jeffrey, A. Table of Integrals, Series, and Products; Elsevier Science: Amsterdam, The Netherlands, 2007. [Google Scholar]
Figure 1. Flow chart of the parameter estimation [24].
Figure 1. Flow chart of the parameter estimation [24].
Sensors 23 04760 g001
Figure 2. Estimation accuracy analysis of k, µ parameter estimation methods based on Monte Carlo method generated k-µ channel data; (a) Comparison and error of theoretical PDF and estimated PDF curves at k = 2, µ = 1; (b) Comparison and error of theoretical PDF and estimated PDF curves at k = 4, µ = 1; (c) Comparison and error of theoretical PDF and estimated PDF curves at k = 4, µ = 2; (d) Comparison and error of theoretical PDF and estimated PDF curves at k = 6, µ = 2.
Figure 2. Estimation accuracy analysis of k, µ parameter estimation methods based on Monte Carlo method generated k-µ channel data; (a) Comparison and error of theoretical PDF and estimated PDF curves at k = 2, µ = 1; (b) Comparison and error of theoretical PDF and estimated PDF curves at k = 4, µ = 1; (c) Comparison and error of theoretical PDF and estimated PDF curves at k = 4, µ = 2; (d) Comparison and error of theoretical PDF and estimated PDF curves at k = 6, µ = 2.
Sensors 23 04760 g002aSensors 23 04760 g002b
Figure 3. Estimation accuracy analysis of k, µ parameter estimation algorithm based on Monte Carlo method generated k-µ channel data; (a) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 1; (b) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 1.5; (c) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 2.
Figure 3. Estimation accuracy analysis of k, µ parameter estimation algorithm based on Monte Carlo method generated k-µ channel data; (a) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 1; (b) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 1.5; (c) Comparison the error of theoretical PDF and estimated PDF curves at k = 0, µ = 2.
Sensors 23 04760 g003
Figure 4. Impact analysis of estimator performance based on different SNR conditions; (a) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 30 dB; (b) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 25 dB; (c) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 20 dB; (d) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 15 dB.
Figure 4. Impact analysis of estimator performance based on different SNR conditions; (a) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 30 dB; (b) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 25 dB; (c) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 20 dB; (d) Error comparison of theoretical and estimated PDF curves at k = 3.5, μ = 1.5, SNR = 15 dB.
Sensors 23 04760 g004aSensors 23 04760 g004b
Figure 5. Impact analysis of estimator performance based on different SNR conditions: (a) The power delay spectrum and PDF values of parameters k and µ when N = 8500; (b) The power delay spectrum and PDF values of parameters k and µ when N = 8000; (c) The power delay spectrum and PDF values of parameters k and µ when N = 5000; (d) The power delay spectrum and PDF values of parameters k and µ when N = 2600; (e) The power delay spectrum and PDF values of parameters k and µ when N = 1000.
Figure 5. Impact analysis of estimator performance based on different SNR conditions: (a) The power delay spectrum and PDF values of parameters k and µ when N = 8500; (b) The power delay spectrum and PDF values of parameters k and µ when N = 8000; (c) The power delay spectrum and PDF values of parameters k and µ when N = 5000; (d) The power delay spectrum and PDF values of parameters k and µ when N = 2600; (e) The power delay spectrum and PDF values of parameters k and µ when N = 1000.
Sensors 23 04760 g005aSensors 23 04760 g005b
Table 1. The development process of the fading channel model.
Table 1. The development process of the fading channel model.
Model NumberModel TypeDegradation to Other ModelsReference
1One-sided Gaussian [5]
2Rayleigh distribution [6]
3Rician distribution2[7]
4Nakagami-m distribution2[8]
5Hoyt (Nakagami-q) distribution [9]
6k-µ distribution1, 2, 3, 4[10]
7k-µ shadowed fading1, 2, 3, 4, 5, 6, 7[4]
Table 2. Algorithm complexity characteristics table.
Table 2. Algorithm complexity characteristics table.
Estimator NameTime Complexity T ( n ) Space Complexity S ( n ) Estimation Efficiency
k ^ - μ ^ Y a c o u b [24] O 1 ( 2 ) + O 2 ( 16 ) = O ( n 0 ) O 3 ( 1 ) = O ( n 0 ) 1
k ^ - μ ^ Δ O 1 ( 15 ) + O 2 ( 10 ) = O ( n 0 ) O 3 ( 3 ) = O ( n 0 ) 2
k ^ - μ ^ A O 1 ( 101 ) + O 2 ( 8 ) = O ( n 0 ) O 3 ( 10 ) = O ( n 0 ) 4
k ^ - μ ^ B O 1 ( 50 ) + O 2 ( 8 ) = O ( n 0 ) O 3 ( 5 ) = O ( n 0 ) 3
Table 3. Total runtime of the estimator characteristics table.
Table 3. Total runtime of the estimator characteristics table.
Estimator NameNumber of Samples 10,000Ranking
Estimated 1 TimeEstimated 50 Time
k ^ - μ ^ Y a c o u b [24]0.0013720.0178951st
k ^ - μ ^ Δ 0.0014810.0181352nd
k ^ - μ ^ A 0.0044150.0484014th
k ^ - μ ^ B 0.0025820.0284653rd
Table 4. Performance comparison table of estimators for k-µ fading channels when µ = 1.
Table 4. Performance comparison table of estimators for k-µ fading channels when µ = 1.
NumberFactorkµRMSEpdfTopFigure
1Truth value21//Figure 2a
Estimation Method k ^ - μ ^ Y a c o u b [24]7.3458032589210.4038662537170.0286172136973rd
k ^ - μ ^ Δ 1.8369355212341.0566791558910.0047322421111st
k ^ - μ ^ A 1.8369355212351.0566791558910.0047322421111st
k ^ - μ ^ B 0.0717070627881.8114097991050.0535106589914th
2Truth value41//Figure 2b
Estimation Method k ^ - μ ^ Y a c o u b [24]1.5576134149531.7836437776020.0168638522103rd
k ^ - μ ^ Δ 13.6385642063930.3704512359680.0097151646561st
k ^ - μ ^ A 13.6385642061030.3704512359750.0097151646561st
k ^ - μ ^ B 0.0950994354422.7861540898650.0487169357614th
3Truth value42//Figure 2c
Estimation Method k ^ - μ ^ Y a c o u b [24]0.5817018674974.9509406937800.0266592320233rd
k ^ - μ ^ Δ 2.6232050119292.7151577142360.0108153506471st
k ^ - μ ^ A 2.6232050119092.7151577142490.0108153506471st
k ^ - μ ^ B 0.0796604086935.6751804091160.0376237138204th
4Truth value62//Figure 2d
Estimation Method k ^ - μ ^ Y a c o u b [24]0.8813371397795.7879757183360.0197748897123rd
k ^ - μ ^ Δ 3.2551055533713.0704589365370.0084965336641st
k ^ - μ ^ A 3.2551055534893.0704589364630.0084965336641st
k ^ - μ ^ B 0.0834463138627.3584768456200.0349702875994th
Table 5. Performance comparison table of estimators for k-µ fading channels when k = 0.
Table 5. Performance comparison table of estimators for k-µ fading channels when k = 0.
NumberFactorkµRMSEpdfTopFigure
1Truth value01//Figure 3a
Estimation Method k ^ - μ ^ Y a c o u b [24]2.6975299727100.4911486587900.0946538479934th
k ^ - μ ^ Δ 1.6588164678010.6501381792590.0658286635722nd
k ^ - μ ^ A 1.6588164678010.6501381792590.0658286635722nd
k ^ - μ ^ B 0.0690358137891.0600387945780.0193380292541st
2Truth value01.5//Figure 3b
Estimation Method k ^ - μ ^ Y a c o u b [24]1.1712626820651.0629441795150.0439703662334th
k ^ - μ ^ Δ 0.6069085897221.3020133079910.0214775543242nd
k ^ - μ ^ A 0.6069085897231.3020133079910.0214775543242nd
k ^ - μ ^ B 0.0370250640651.5167082867520.0036928163351st
3Truth value02//Figure 3c
Estimation Method k ^ - μ ^ Y a c o u b [24]2.4428537844560.9714941411330.0566003624964th
k ^ - μ ^ Δ 0.3955583805751.8422341812660.0121513411002nd
k ^ - μ ^ A 0.3955583805761.8422341812660.0121513411002nd
k ^ - μ ^ B 0.0242424783562.0020438862810.0004346398101st
Table 6. The comparison of estimation results under different SNR conditions.
Table 6. The comparison of estimation results under different SNR conditions.
SNR (dB)FactorkµRMSEpdfTopFigure
Truth Value3.51.5
30 dBEstimation Method k ^ - μ ^ Y a c o u b 8.26500.75670.0109/Figure 4a
k ^ - μ ^ Δ / k ^ - μ ^ A 2.38381.87360.00811
k ^ - μ ^ B 0.07773.70020.03872
25 dBEstimation Method k ^ - μ ^ Y a c o u b 1.64242.20410.0205/Figure 4b
k ^ - μ ^ Δ / k ^ - μ ^ A 2.79581.63920.01971
k ^ - μ ^ B 0.08093.56290.04022
20 dBEstimation Method k ^ - μ ^ Y a c o u b 2.82941.46640.0554/Figure 4c
k ^ - μ ^ Δ / k ^ - μ ^ A 5.73610.88520.05571
k ^ - μ ^ B 0.09033.19840.06052
15 dBEstimation Method k ^ - μ ^ Y a c o u b −1.4877−20.7926//Figure 4d
k ^ - μ ^ Δ / k ^ - μ ^ A 2.71611.24830.11532
k ^ - μ ^ B 0.08032.66520.11291
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, J.; Fan, Z.; Ji, Z.; Li, X.; Fei, P.; Hou, D. Multi-Parameter Estimation Method and Closed-Form Solution Study for k-µ Channel Model. Sensors 2023, 23, 4760. https://doi.org/10.3390/s23104760

AMA Style

Tian J, Fan Z, Ji Z, Li X, Fei P, Hou D. Multi-Parameter Estimation Method and Closed-Form Solution Study for k-µ Channel Model. Sensors. 2023; 23(10):4760. https://doi.org/10.3390/s23104760

Chicago/Turabian Style

Tian, Jie, Zhongqing Fan, Zhengyu Ji, Xianglu Li, Peng Fei, and Dong Hou. 2023. "Multi-Parameter Estimation Method and Closed-Form Solution Study for k-µ Channel Model" Sensors 23, no. 10: 4760. https://doi.org/10.3390/s23104760

APA Style

Tian, J., Fan, Z., Ji, Z., Li, X., Fei, P., & Hou, D. (2023). Multi-Parameter Estimation Method and Closed-Form Solution Study for k-µ Channel Model. Sensors, 23(10), 4760. https://doi.org/10.3390/s23104760

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