Analysis method for multi-user random access signals

Abstract

An analysis method for multi-user random access signals is disclosed, which solves the connection problem between a base station and a user equipment when the narrow-band IoT uplink signal transmission is implemented. The analysis method of the present invention utilizes the detection threshold value to effectively determine multiple user equipment to be signaled; then, for each detected user equipment, an effective method based on the phase difference is used to estimate their synchronization parameters, i.e. time-of-arrival (ToA) and residual carrier frequency offset (RCFO). Therefore, according to the present invention, random access signals can be received correctly and efficiently and the user equipment related information can be obtained at the same time to facilitate subsequent communications.

Description

This application claims priority for Taiwan patent application no. 107116680 filed on May 16, 2018, the content of which is incorporated by reference in its entirely.

BACKGROUND OF THE INVENTION

Field of the Invention

The present invention relates to an analysis method for random access signals, particularly to an analysis method for random access signals of multiple user equipment in an NB-IoT uplink system.

Description of the Related Art

In the NB-IoT standard of the 3G Partnership Project (3GPP), a brand new channel, i.e. the Narrowband Physical Random Access Channel (NPRACH), is defined in the uplink signals of user equipment to increase the transmission distance between user equipment and base stations. The NPRACH adopts a preamble signal in the format of single-tone frequency hopping (STFH). The preamble signal is the first signal user equipment transmits to a base station, also called Message-1, expressing that the user equipment intends to link with the base station. Therefore, it is an important problem for the base station side: how to correctly and effectively receive the NPRACH preamble signal and simultaneously acquire the related information of the user equipment to facilitate subsequent communication.

In the NPRACH structure, a preamble signal is normally composed of 4 symbol groups. While the preamble signal is used as the unit of channel hopping, each symbol group is composed of 5 symbols and a cyclic prefix. While there is a plurality of user equipment, the NB-IoT system having a bandwidth of 180 kHz can support the wireless access of at most 48 user equipment and the transmissions of the NPRACH preamble signals thereof to the base station. The kth user equipment (UE_k) is defined to carry an n_init (k) parameter, wherein k ranges from 0 to 47. Thereby, the channel hopping can be undertaken according to the preamble signals generated by the single-tone symbol groups. Thus, the NPRACH of the ith symbol of the nth symbol group of UE_k is a single-tone signal, and the subcarrier frequency thereof is f_k(n)×Δf(n), wherein f_k(n)∈[0,47], expressing the index of the channel hopping carrier of the nth symbol group of the kth user equipment.

FIG. 1 shows 48 possible channel hopping groups generated according to the 3 GPP R13 NPRACH standard, wherein {f_k(n), k=0˜47, n=0˜31}. As shown in FIG. 1, anyone of the channel hopping groups would not collide with another of the channel hopping groups. In other words, each of the hopping groups is different from the others of the hopping groups. In the single-tone fundamental frequency vector signal, which is transmitted by UE_k and has a length of L×1, the ith symbol of the nth symbol group may be expressed bys_k,n,i(I)=exp(j2πf_k(n)1/L)  (1)

The time-domain vector signal of the nth symbol group may be expressed byg_k,n=CP{REP₅(s_k,n,i)}  (2)wherein REP_r(a) expresses the repeated concatenation of the same column vector a for r times; CP {●} expresses the operation mode of CP, resulting in the length of 5L+L_cp. Four symbol groups constitute a preamble signal, and the NPRACH signal may consist of a number of preamble signals. Similar to the conventional LTE uplink system, the whole NPRACH fundamental frequency signal has an additional Δf/2 frequency shift, which is converted by up-conversion to the central frequency that can be transmitted.

In the signal-receiving mode of a base station eNodeB (eNB) involving a plurality of user equipment, as different user equipment (UEs) may have different distances to the base station eNB, there are different delays for the NPRACH signals transmitted therebetween and different time-of-arrivals (ToAs). D_k ∈[0, L_cp−1]. The channel between each user equipment and the base station eNB can be expressed by a single-tap flat Rayleigh fading channel. Because the channel varies very slowly, the signal fading of UE_k can be simplified to be a complex fading coefficient h_k˜CN(0,σ_k²), wherein σ_k² expresses the average received signal power of UE_k at the eNB. Each user equipment also has a residual carrier frequency offset (RCFO) parameter η_k, which has only a very small value. The signal received by eNB is formed by superimposing the fading NPRACH preamble signals of multiple user equipment and has AWGN (Additive White Gaussian Noise) with an power of σ_n². Therefore, the input time-domain signal-to-noise power ratio SNR for UE_k can be expressed bySNR_i(k)=σ_k²/σ_n²  (2)

At eNB, the received RF signal is first down converted and frequency shifted by −Δf/2. Then, the baseband receiver performs the CP removal and taking 512-point FFT for each symbol in the symbol groups according to the timing of eNB. As the residual carrier frequency offset (RCFO) is quite small in practice, the inter-carrier interference (ICI) and the multiple access interference (MAI) can be neglected. After some derivation, it can be shown that the post-FFT signal vector of the ith symbol of the nth symbol group has its f_k(n)^−th FFT bin as follow:

$\begin{array}{cc}\begin{array}{c}R\left(n,i,f_k\left(n\right)\right)=Lh_kS_L\left({\eta }_k\right)\times \left[\begin{array}{c}{e}^{-j2\pi D_k{f}_k\left(n\right)/L}\times e^{j2\pi {\eta }_{k}t\left(n,i\right)/L}\times \\ e^{j\left({\theta }_k+{\varphi }_k+{\phi }_k\right)}\times B\left({\eta }_k\le 0\right)e^{-j\pi }\end{array}\right]+W\left(n,i\right)\\ =Lh_kS_L\left({\eta }_k\right)e^{j{\beta }_{k,n,i}}+W\left(n,i\right)\end{array}& \left(4\right)\end{array}$ wherein

$S_L\left(\eta \right)=\frac{\sin \left(\pi \eta \right)}{L\sin \left(\pi \eta /L\right)};$ t(n,i)=n(5L+L_cp)+iL; ϕ_k=−πη_k(L−1)/L; ϕ_k=∠h_k is the channel phase coefficient; θ_k˜∪(0,2π) is the carrier phase offset; if η_k≤0, B(η_k≤0)=1; alternatively, b(η_k≤0)=0; β_k,n,i=2π(η_kt(n,i)−D_kf_k(n))/L+θ_k+ϕ_k×φ_k−B(η_k≤0)π; W(n,i) is the post-FFT AWGN noise samples with power L×σ_n².

In the case that UE_k is absent, R(n,i,ϕ_k (n))=W(n,i) and is expressed by CN(0,Lσ_n²). In the case that UE_k exists, the residual carrier frequency offset is very small (S_L(η_k)≈1). Thus is obtained Equation (5):R(n,i,ϕ_k(n))=L|h_k|e^jβk,n,i+W(n,i)  (5)which may be expressed by CN(0, L²σ_k²+Lσ_n²). The absolute value L|h_k| of the signal component is not fluctuating over the symbols. Thus, the post-FFT SNR of UE_k becomes SNR_o=(L²σ_k²)/(Lσ_n²)=L×SNR_i, which is L times of SNR_i. Therefore, the single-tone signal can be coherently summed up.

Based on the above discussion, the problems of detection and estimation in the multi-UE NPRACH are summarized as follows: in order to catch the presence and identity the parameter {n_init(k)} of the UEs, the performance of detection should achieve an false alarm probability P_F≤0.1% and a detection probability P_D>99%; it is also very important for detecting all user equipment to precisely and efficiently estimate ToA and RCFO.

In 2016, Lin et al. proposed a paper “Random Access Preamble Design and Detection for 3GPP Narrowband IoT systems” published in IEEE Wireless Communications Letters, vol. 5, no. 6, pp. 640-643, December 2016, which is to integrally deal with the problems of detection and estimation in NPRACH, and intends to use the peak values of 2D-FFT to detect the existence of user equipment and estimate ToA/RCFO. However, the conventional technology involves several problems: it is hard to perform the critical 2D-FFT of the matrix of the fundamental frequency data; the practical computation is much more complicated than the proposed method; the paper does not described the detection threshold and performance analysis in detail.

Accordingly, the present invention proposes an analysis method for multi-user random access signals to solve the problems of the conventional technologies.

SUMMARY OF THE INVENTION

The primary objective of the present invention is to provide an analysis method for multi-user random access signals, which is to solve the linking problems of user equipment and base stations while an NB-IoT uplink signal is transmitted, including the problem of detecting signal accessing and the problem of estimating the time synchronization parameter and the frequency synchronization parameter, whereby can be achieved an overall detection performance with an false alarm probability P_F≤0.1% and a detection probability P_D>99%, and whereby all user equipment can be efficiently detected, and whereby ToA and RCFO can be precisely estimated, wherefore the present invention can achieve better performance with less computation.

Another objective of the present invention is to provide a low-complexity analysis method for multi-user random access signals, which has advantages of higher precision, less computation and definite threshold values.

In order to achieve the abovementioned objectives, the present invention proposes an analysis method for multi-user random access signals, which applies to an NB-IoT uplink system, and which comprises steps: receiving the preamble signals of the random access signals from a plurality of user equipment, detecting a plurality of symbol groups of each preamble signal and acquiring the corresponding average power, and comparing each average power with a detection threshold value to determine whether the user equipment intends to access the base station; after determining the user equipment intending to access the base station, acquiring the phase trace of the preamble signal of each detected user equipment, and calculating parameters of ToA and RCFO according to the phase difference of adjacent symbol groups of the phase trace.

While the average power is greater than the detection threshold value, it means that the corresponding user equipment intends to access the base station. While the average power is lower than the detection threshold value, it means that the corresponding user equipment does not intend to access the base station.

In the present invention, the abovementioned detection threshold value is a Neyman Pearson threshold value. The detection threshold value is determined using a false alarm level and a decision delay of the random access signal in cooperation with a detected noise power. The decision delay is the number of all the symbol groups of the random access signal. For example, each preamble signal has four symbol groups; the number of the preamble signal multiplied by four is the number of all the symbol groups.

Below, embodiments will be described in detail in cooperation with the attached drawings to make easily understood the objectives, technical contents and accomplishments of the present invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram schematically showing 48 possible channel hopping groups generated according to the 3GPP R13 NPRACH standard;

FIG. 2 is a diagram schematically showing a multi-user system applied to NB-IoT according to one embodiment of the present invention;

FIG. 3 is a flowchart of an analysis method according to one embodiment of the present invention;

FIG. 4 is a diagram showing the relationship between the scale factor A(N; α=0.1%) and the decision delay N;

FIG. 5 is a diagram showing theoretical curves of the relationship of P_D,AWGN to SNR_i for different decision delays at α=0.1% according to one embodiment of the present invention; and

FIG. 6 is a diagram showing theoretical curves of the relationship of P_D,Fading to SNR_i for different decision delays at α=0.1% according to one embodiment of the present invention.

DETAILED DESCRIPTION OF THE INVENTION

Refer to FIG. 2, in the random access process of the NB-IoT uplink system, the first piece of uplink signal of the user equipment 12 is called NPRACH, which uses a single-tone frequency hopping preamble sequence. The time unit of frequency hopping is a symbol group. According to the cellular coverage and the distance between the base station and the user equipment, a base station (eNB) 10 specifies the length of the preamble sequence, which is to be sent out by the user equipment 12, in the information block of the downlink system thereof. A detection threshold value is set in the NPRACH receiver algorithm to achieve the false alarm probability and the detection probability, which are regulated by the NB-IoT standard, and to estimate the parameters of ToA and RCFO.

The present invention proposes an analysis method for multi-user random access signals, which is extensively applicable to the NB-IoT uplink system. Refer to FIG. 2 and FIG. 3, the analysis method of the present invention comprises Steps 10-30. In Step S10, the base station 10 receives the preamble signals of the random access signals from a plurality of user equipment 12, detects a plurality of symbol groups of each preamble signal, and acquires the corresponding average energies. In other words, the base station 10 collects the bits generated by the post-FFT of all the symbol groups and detects the sufficient statistics of the user equipment to determine the average power.

In Step S20, determine whether there is any user equipment 12 intending to access the base station 10 according to the corresponding average power and a detection threshold value. In Step S22, the average power is lower than the threshold value, and the corresponding user equipment 12 does not intend to access the base station 10; then, the process would not proceed to the next step but ends herein. In Step S24, the average power is higher than the threshold value, and the corresponding user equipment 12 intends to access the base station 10; then the process proceeds to Step S30. The abovementioned detection threshold value is a Neyman Pearson threshold value. The Neyman Pearson threshold value is determined using a false alarm level and a decision delay, which are obtained beforehand from the random access signal, in cooperation with a detected noise power. The decision delay is the number of all the symbol groups of the random access signal.

In Step S30, while having detected the user equipment 12 intending to access the base station 10, the base station 10 acquires the phase trace of the preamble signal of the corresponding user equipment 12 and calculates the ToA parameter and the RCFO parameter of the user equipment 12 according to the phase difference of the adjacent symbol groups of the phase trace. In Step S30, it is according to the phase trace that the ToA parameters and the RCFO parameters are sequentially estimated. In other words, the phase difference induced by channel hopping of the symbol groups is calculated firstly, and all phase differences corresponding to all preamble signals are averaged, whereby to acquire the RCFO parameter corresponding to the user equipment; next, the average phase of each symbol group is calculated to acquire the average phase difference corresponding to each symbol group; the average phase differences are summed up to obtain the related ToA parameter.

After the technical characteristics of the present invention have been described above, the principles of the present invention will be described thereinafter to prove that the analysis method for multi-user random access signals is practicable and easy to practice.

As there are 48 distinct and possible hopping groups, the test data Q_k of each user equipment is used to compensate for the bits of the corresponding post-FFT. The test data Q_k of UE_k is defined asQ_k={R(n,i,f_k(n)),n=0, . . . ,N−1,i=0,1, . . . 4}wherein N=4P and means the number of all the symbol groups; P is the number of the preamble signals for detection; N is also regarded as the decision delay in the user equipment detection process. For different coverage class, the base station eNB may designate P to be 1, 2, 4, 8, 16, 32, 64, and 128, which are respectively corresponding to the decision delays N of 4, 8, 16, 32, 64, 128, 256, and 512. Therefore, if the user equipment UE_k has a lower receiving SNR, it is necessary to collect more data, i.e. obtain larger N, leading to a longer decision delay for the farther or weaker user equipment.

As the hopping groups are orthogonal to each other, the superimposed NPRACH detection problems can be decoupled into a parallel of single-UE detection problem. In order to test the presence of UE_k, the energy of the bits of the post-FFT of N symbol groups is collected; then, the sufficient statistics of detecting the user equipment is used to determine the average power:

$\begin{array}{cc}{P}_k\left(N\right)=\frac{1}{5N}\sum _{n=0}^{N-1}\sum _{i=0}^4{R\left(n,i,f_k\left(n\right)\right)}^2& \left(7\right)\end{array}$

Next, the decision rule used to determine the presence of UE_k can be simplified comparing P_k(N) with a threshold value λ, which may be expressed by

$\begin{array}{cc}{P}_k\left(N\right)\stackrel{{\mathrm{UE}}_k\mathrm{present}}{\underset{{\mathrm{UE}}_k\mathrm{abscent}}{\gtrless }}\lambda & \left(8\right)\end{array}$

The receiver operation characteristic (ROC) or the performance test can be expressed by the false alarm probability P_F P_D=Pr(P_k(N)>λ/UE_k is present)  (9)

The detection probability P_D is expressed byP_D=Pr(P_k(N)>λ/UE_k is present)  (10)

Therefore, the main problem of NPRACH detection is to specify the threshold value λ. The NB-IoT standard demands that the detection probability P_D must exceed 99% and the false alarm probability P_F should not exceed 0.1%. In order to achieve the two conditions, the present invention uses the Neyman Pearson rule to resolve the detection problem, which is expressed bymax{P_D}, such that P_F≤α  (11)

It means that the decision rule is most powerful at the significant level α for the threshold value λ. In the case of NB-IoT, it may be selected that α=0.1%. According to the Neyman Pearson rule, the threshold value λ is the function of the decision delay N and has two parameters of a specified false alarm level α and a noise power P_n=Lσ_n². Suppose that the receiver can detect the noise power P_n. Thus, the threshold value λ must satisfy the equation:∫_λ^∞g_n(x)dx=α  (12)wherein g_n(x) is a probability density function (PDF) of test statistics P_k(N) under the noise-only case. Since P_k(N) is the average power of 10×N independent real Gaussian random variables (RVs) with identical distribution N(0, P_n/2). Therefore, PDF g_n(x) is determined by a scaled central Chi-square distribution with 10×N degrees of freedom as follows:g_n(x)=κf_c(κx;10N)  (13)wherein

$\begin{array}{cc}{f}_c\left(x;m\right)=\frac{1}{2^{k/2}\Gamma \left(k/2\right)}{x}^{\frac{k}{2}-1}{e}^{-\frac{x}{2}}u\left(x\right)& \left(14\right)\end{array}$ which is a standard χ_m² PDF with in degrees of freedom, wherein κ=(10N)/(Lσ_η²)=10N/P_n is a scale factor. The cumulative distribution function (CDF) of the standard χ_m² PDF is determined by

$\begin{array}{cc}{F}_c\left(x;m\right)={\int }_0^x{f}_c\left(\lambda ;m\right)d\lambda =\frac{1}{\Gamma \left(\frac{k}{2}\right)}\gamma \left(\frac{k}{2},\frac{x}{2}\right)u\left(x\right)& \left(15\right)\end{array}$ wherein γ(s,t)=∫₀^xt^s−1e^=tdt is a lower incomplete Gamma function.

It is learned from Equation (11): the false alarm level α and the decision delay N can be used to definitely determine the optimized Neyman Pearson threshold value λ_o with the equation:

$\begin{array}{cc}\begin{array}{c}{\lambda }_o=\frac{1}{\kappa }{F}_c^{-1}\left(1-\alpha ;10N\right)\\ =P_n\times \left(\frac{1}{10N}{F}_c^{-1}\left(1-\alpha ;10N\right)\right)\\ =P_n\times A\left(N;\alpha \right)\end{array}& \left(16\right)\end{array}$ wherein F_x⁻¹(x; m) is the inverse function of χ_m² CDF, i.e. F_c⁻¹(F_c(x; m); m)=x.

From Equation (16), it is noted to determine the threshold for a given α and N. The average noise power P_n of post-FFT can be precisely measured using those noise-only resource grids. The threshold value is obtained via multiplying P_n by the scale factor A(N; α). FIG. 4 shows the relationship of the scale factor A(N; α) and the decision delay N. The scale factor is always greater than one, but it decreases with the increase of N.

The 48 average energies of the channel hopping mode along different user equipment can be worked out using only 48 values of 512-point FFT. The scale factor A(N; α=0.1%) can be calculated beforehand and picked up from the memory for application. Therefore, the present invention can realize NPRACH detection in very low complexity and very low power consumption.

Next is deduced the detection performance in the AWGN channel and the Rayleigh fading channel. For the AWGN channel, suppose that the composite fading coefficient h_k=1 and that R(n,i,ϕ_k (n)) has a non-zero average value Le^jβk,n,i. In the AWGN channel, the signal power of post-FFT is P₀=L², which is nonrandom. With the signal and noise, a scaled non-central Chi-square distribution PDF may be used to express the statistical data P_k(N) as follows:g_s(x)=κf_nc(κx;10N,γ₀)  (17)

• wherein k is the same as the scale factor in Equation (13).

$\begin{array}{cc}{f}_{\mathrm{nc}}\left(x;m,{\gamma }_0\right)=\frac{1}{2}{e^{-\left(x+\gamma \right)/2}\left(\frac{x}{{\gamma }_0}\right)}^{k/4-1/2}{I}_{k/2-1}\left(\sqrt{{\gamma }_{0}x}\right)u\left(x\right)& \left(18\right)\end{array}$ Equation (18) is a standard non-central Chi-square distribution PDF with m degrees of freedom, wherein the non-central parameter cis given by the following equation:γ₀=10N×P₀/P_n=10NL/σ_n²  (19)

• l_v(x) is the modified Bessel function of the first kind with degree ν. CDF of f_nc(x; m, γ₀) is expressed by

$\begin{array}{cc}{F}_{\mathrm{nc}}\left(x;m\right)={\int }_0^x{f}_{\mathrm{nc}}\left(\lambda ;m,{\gamma }_0\right)d\lambda =1-Q_{\frac{k}{2}}\left(\sqrt{{\gamma }_0},\sqrt{x}\right)& \left(20\right)\end{array}$ wherein

$\begin{array}{cc}{Q}_M\left(a,b\right)={\int }_b^{\infty }{x\left(\frac{x}{a}\right)}^{M-1}\exp \left(-\frac{x^2+a^2}{2}\right)I_{M-1}\left(\mathrm{ax}\right)\mathrm{dx}& \left(21\right)\end{array}$ Equation (21) is a generalized Marcum Q-function; a and b are the parameters substituted into the equation for computation.

After the Neyman Pearson threshold value λ₀ is acquired, the theoretical detection probability P_D under the AWGN channel may be derived asP_D,AWGN(γ₀)=Q_5N(√{square root over (γ₀)},√{square root over (10N×λ_o/P_n)})  (22)

P_D,AWGN(γ₀) may be used to calculate the detection probability under the Rayleigh fading channel Letting σ_k²=1 without loss of generality, then ρ=|h_k|² is an exponential redundancy version (RV) with unity mean. The average signal power of post-FFT P_s=E[ρP₀]=P₀ is unchanged and the non-central parameter becomes γ=ργ₀. Therefore, the detection probability under the Rayleigh fading channel may be obtained by averaging P_D,AWGN(γ₀) over the PDF of ρ and may be expressed by the following implicit integral expression:P_D,Fading(γ₀)=∫₀^∞e^−ρQ_5N(√{square root over (ργ₀)},√{square root over (10N×λ_o/P_n)})dρ  (23)Equation (23) may have a complicated closed-form formula. However, the present invention can use numerical integration to calculate the theorectical detection performance under the Rayleigh fading channel.

Refer to FIG. 5 and FIG. 6 respectively showing the theorectical curves of the relationships of P_D,AWGN and P_D,Fading to SNR_i for different decision delays at α=0.1%. In order to verify the analyses, the performance curves is simulated with the Monte-Carlo method. The results show that the simulation curves almost coincide with the theorectical curves. It is observed in these curves: a longer preamble signal, i.e. a larger N, can significantly improve the UE detection performance of eNB, especially while SNR is low. In other words, eNB can use a longer NPRACH configuration to extend the coverage thereof.

After NPRACH detection has been discussed above, the present invention proposes an algorithm that can effectively estimate the synchronization parameters of ToA and RCFO. The algorithm integrates the abovementioned detection methods, addressing only the detected user equipment and using the decoupling method described below to estimate the parameters of ToA and RCFO.

For the detected user equipment UE_k, the algorithm starts from its unwrapped phase trace of R(n,i, f_k(n)) expressed as:q_k,n,i=unwrap{arg{R(n,i,f_k(n))}}  (24)In the noise-free case, the phase trace can be obtained using q_k,n,i=β_k,n,i=−2πD_kf_k(n)−2πη_k(n,i)+C in Equation (4), wherein C is a constant phase. As the first two phase terms are directly related to RCFO and ToA parameters, the phase trace, together with suitable phase differences, can be used to estimate the two synchronization parameters. Further, average operation can be used to decrease the estimation variance caused by AWGN.

For each symbol group of 5 symbols, there are four adjacent phase differences that can be calculated:ε_k,n,i=q_k,n,i+1−q_k,n,i for i=1,2,3,4,  (25)Next, by averaging ε_k,n,i, over all n and i in the preamble, it can obtain the RCFO estimate for UE_k, as shown by the equation:

$\begin{array}{cc}{\hat{\eta }}_k=\frac{1}{2\pi }\times \frac{1}{4N}\sum _{n=0}^{N-1}\sum _{i=1}^4{\varepsilon }_{k,n,i}& \left(26\right)\end{array}$ Next, in order to estimate ToA parameter, the derivative terms of RCFO is removed from q_k,n,i to obtain Equation (27):q_k,n,i=q_k,n,i−2π{circumflex over (η)}_kΔft_n,i  (27)wherein t_n,i=[(5n+i)L+(n+1)L_cp]T_s is the starting time instant of the ith symbol of the nth symbol group. The average phase of each symbol group q_k,n,i may be expressed by

$\begin{array}{cc}{z}_{k,n}=\frac{1}{5}\sum _{i=0}^4{\stackrel{\_}{q}}_{k,n,i}& \left(28\right)\end{array}$ Then, the sum of the differences of z_k,n is divided by the sum of the differences of channel hopping to obtain ToA, as shown by the equation:

$\begin{array}{cc}{\hat{D}}_k=\frac{{\alpha }_k}{2\pi T_s}\times \frac{\sum _{n=1}^{N-1}z_{k,n+1}-z_{k,n}}{\Delta f\times \sum _{n=1}^{N-1}f_k\left(n+1\right)-f_k\left(n\right)}& \left(29\right)\end{array}$ wherein the sign of ToA may be expressed by

$\begin{array}{cc}{\alpha }_k=\mathrm{sign}\left(\frac{-\sum _{n=1}^{N-1}\left(z_{k,n+1}-z_{k,n}\right)}{\sum _{n=1}^{N-1}\left(f_k\left(n+1\right)-f_k\left(n\right)\right)}\right)=+1\mathrm{or}-1& \left(30\right)\end{array}$

Hence, for each detected UE, the above joint RCFO/ToA synchronization algorithm is straight forward, easily implemented, and computationally efficient.

In conclusion, the present invention proposes an analysis method for multi-user random access signals, which can solve the linking problem between the base stations and the user equipment, and which can make the overall detection performance achieve a false alarm probability P_F≤0.1% and a detection probability P_D>99%, and which can accurately and efficiently detect each user equipment and estimate the synchronization parameters of ToA and RCFO to facilitate subsequent communications. Therefore, the present invention is a high-precision and low-computation burden analysis method with a definite threshold value.

The embodiments are described above to demonstrate the technical contents and characteristics of the present invention to enable the persons skilled in the art to understand, make, and use the present invention. However, these embodiments are only to exemplify the present invention but not to limit the scope of the present invention. Any equivalent modification or variation according to the spirit of the present invention is to be also included by the scope of the present invention.

Claims

1. An analysis method for multi-user random access signals, which applies to an uplink system of a narrowband Internet of thing (NB-IoT), and which comprises steps: receiving preamble signals of random access signals from a plurality of user equipment, detecting a plurality of symbol groups of each said preamble signal and acquiring corresponding average power, and comparing said average power with a detection threshold value to determine said user equipment intending to access signals; andacquiring a phase trace of said preamble signal corresponding to each of said user equipment intending to access signals, and calculating parameters of Time of Arrival (ToA) and Residual Carrier Frequency Offset (RCFO) according to phase differences of adjacent said symbol groups of said phase trace, wherein said detection threshold value is a Neyman Pearson threshold value, at least one base station collects bits generated by post-FFT (Fast Fourier Transform) of said symbol groups and detects sufficient statistics of said user equipment to determine said average power, and while said average power is greater than said detection threshold value, said user equipment corresponding to said average power intends to access said base station; while said average power is lower than said detection threshold value, said user equipment corresponding to said average power does not intend to access said base station.

2. The analysis method for multi-user random access signals according to claim 1, wherein said threshold value is determined using a false alarm level and a decision delay of said random access signals, in cooperation with a detected noise power; said decision delay is a count of all said symbol groups of said random access signal.

3. The analysis method for multi-user random access signals according to claim 1, wherein said preamble signals includes 4 symbol groups.

4. The analysis method for multi-user random access signals according to claim 3, wherein said step of calculating parameters of Time of Arrival (ToA) and Residual Carrier Frequency Offset (RCFO) according to said phase differences of adjacent said symbol groups of said phase trace further includes steps: calculating phase differences caused by channel hopping of said symbol groups;averaging all said phase differences corresponding to all said preamble signals to obtain said parameter of Residual Carrier Frequency Offset (RCFO);calculating average phase of each said symbol group to obtain an average phase difference corresponding to each said symbol group; andsumming up said average phase differences to obtain said parameter of Time of Arrival (ToA).