# Research on Spectrum Feature Identification of Indoor Multimodal Communication Signal

Authors: Yunfei Chen; Yang Liu; Xintao Fan
Journal: Advances in Mathematical Physics
 (2021)

---

## Abstract

In order to solve the problem of large signal acquisition error caused by radio wave multipath effect in indoor environment, firstly, the signal source carried on the motion platform is collected for spectrum signal, and the signal processed by wavelet threshold denoising algorithms extracted and stored for spectrum feature extraction. Then, after data training and identification, the signal source is input into the system in random mode for identification. The experimental results show that the improved fuzzy clustering algorithm (FCA) is 12.7% higher than the spectrum envelope extraction method (SEEM) in the recognition rate of spectrum characteristics of different modes of signal source.


---

## Body

## 1. Introduction

In indoor positioning of motion platform, communication signal mode plays a key role in information exchange. Due to the characteristic of indoor radio wave transmission [1] of different modes of communication signal, the working mode can be automatically changed to adjust the channel when the communication [2] is not smooth. The positioning access point AP can identify the signal source on the sports platform to strengthen the positioning and tracking. The spectrum feature identification method is very important for the identification of wave signal is very important. The traditional envelope spectrum feature extraction method [3] can optimize the selection of signal source spectrum features [4] and effectively improve the signal detection accuracy of signal source. However, the feature extraction of the SEEM algorithm is easily limited by the experimental environment. When the indoor environment changes greatly, the error of feature data extracted by the SEEM algorithm is large.The FCA is a clustering algorithm based on the fuzzy system theory [5]. Bezdek proposed the FCA [6] as early as 1973, which originated from the improved early K-means clustering algorithm. The objective function constructed by the FCA can achieve the maximum similarity of samples belonging to the same cluster and the minimum correlation of samples between different clusters in the division process. However, the randomness of the FCA clustering center is large, which easily leads to local extremum [7]. At the same time, when the attributes of test samples are close to each other, it is difficult to classify them. All these need to guide the classification. The improved FCA based on pairwise regression fusion can be combined with the FCA in fuzzy classification, so as to improve the classification efficiency.


## 2. Related Work

The method of analyzing the signal characteristics by using the spectrum of the space wave signal received by the spectrum analyzer receiver has been successfully applied in many fields. At present, some scholars have studied the spectrum characteristics of outdoor GPS satellite signal [8], radar signal, and radio signal. However, there are few studies on the spectrum characteristics of indoor interference signal source. In Reference [9], the distortion of the emitter modulator and the nonlinear characteristics of the power amplifier are introduced into the classifier construction model, and the test results show that the spectrum distortion of the emitter signal exists. In Reference [10], the transient sparse feature of signal is used as the basis of evaluation mode to realize the identification of emitter signal feature. In Reference [11], the relevant features are extracted by analyzing the time-frequency characteristics of the emitter signal to realize the identification of transient characteristics. The above research results are significant, but none of them is involved in the research of indoor positioning signal source spectrum characteristics.


## 3. Methodology

### 3.1. Building Models

In order to solve the problem of large signal acquisition error caused by the multipath effect of indoor signal source, firstly, the spectrum of the radio signal is collected and processed on the indoor platform, and then, the random signal is input to the system to identify the large-scale signal. At last, the construction of signal acquisition and identification system is shown in Figure1.Figure 1
Structure diagram of signal acquisition and identification system.The known communication signal source is collected and input into the spectrum identification system by the SA44B spectrum analyzer. According to the four common communication modes of GSM, CDMA, DCS, and PHS, the corresponding spectrum features are extracted and stored in the database after signal detection feature extraction and identification. Then, after learning and training, a higher frequency spectrum feature recognition rate is achieved. In the online identification stage, the signal source with unknown working mode is input into the identification system, and the working mode of the output signal source is determined after identification.


### 3.2. Signal Acquisition and Processing

This section is to distinguish the source signal of the target for positioning; the four modes of spectrum signals of the signal source carried by the indoor sports platform are collected indoors. The spectrum of CDMA signal source has the characteristics of smooth waveform and obvious peak, as shown in Figure2. The spectrum of signal source in GSM mode is coarser than that in CDMA mode, but it has clear peak, as shown in Figure 3. The spectrum of DCS working mode is the most disordered among the four kinds of spectrum collected, as shown in Figure 4. In PHS mode, the peak to peak distance of communication signal spectrum is larger, so the characteristics are more obvious, as shown in Figure 5.Figure 2
Spectrum of CDMA mode.Figure 3
Spectrum of GSM mode.Figure 4
Spectrum of DCS mode.Figure 5
Spectrum of PHS mode.Obviously, the characteristics of the four modes of spectrum signals are complex, the optimal characteristic parameters of the spectrum are not easy to select, and there are many clutters in the spectrum, which affect the indoor positioning of the signal source.In order to prevent the edge information from blurring, a nonlinear filter is introduced, in which the median filter can keep the sharp edge well. The signal processed by median filter still contains noise, and the wavelet denoising method is better for reducing Gaussian white noise in radio wave signal. The wavelet denoising method [12] includes the modulus extremum denoising method, wavelet correlation denoising method, and wavelet threshold denoising method. The traditional wavelet hard thresholding function [13] and soft thresholding method [14] models are shown in(1)W⌢j,k=sgnWj,k·Wj,k−λ,Wj,k≥λ,0,Wj,k<λ,W⌢j,k=Wj,k,Wj,k≥λ,0,Wj,k<λ.The traditional wavelet hard thresholding function is easy to produce Gibbs oscillation [15], while the soft thresholding method is easy to produce “over smooth” distortion due to the constant difference of wavelet coefficients [16].(2)W⌢j,k=aWj,k+1−asgnWj,k·Wj,k−λexpWj,k−λδ,Wj,k≥λ,0,Wj,k<λ.As shown in formula (2), the threshold adjustment parameter a is a=1−e−αWj,k−λ2, W′ is the mean value of neighborhood wavelet coefficients, and δ is a positive number. When eWj,k⟶±λ is W⌢j,k⟶0, the threshold value Wj,k=±λ function is continuous, which can overcome the signal oscillation problem caused by the discontinuity of hard threshold; when Wj,k⟶±∞ is W⌢j,k/Wj,k⟶1, the new threshold function conforms to the characteristics of hard threshold function, which can overcome the hard threshold, so as to overcome the distortion caused by the constant difference of wavelet coefficients.When the threshold adjustment parametersα⟶0 and δ⟶0, the new wavelet threshold function is a soft threshold function; when the threshold adjustment parameters α⟶+∞ and δ⟶1, the new wavelet threshold function is a hard threshold function. By adjusting the parameters a and δ, we can obtain the signal processing results of different modes of the interference signal source. After denoising the spectrum of the four modes, we can normalize the Z-score [17], as shown in Figure 6.Figure 6
Spectrum normalization of four interference modes.


### 3.3. Spectrum Feature of Signal Extraction

The signal spectrum parameters are selected to extract the peak value, corresponding frequency, signal spectrum bandwidth, and kurtosis [18] of the signal transmitted by different mode signal sources as the characteristic parameters of the corresponding mode of the radio wave signal. The feature parameters of different modes of CDMA, GSM, DCS, and PHS of signal source are extracted step by step to construct the vector sequence required for the feature recognition model of radio wave signal. As shown in Table 1, the four modes of CDMA, GSM, DCS, and PHS are, respectively, marked as C=1~4.Table 1
Partial experimental characteristic data.
RSSIM (dBm)fM (103 MHz)fband (103 MHz)ψC-29.320.8670.02857.2311-30.510.8730.02756.5381-41.651.8060.06251.2672-40.371.8470.05850..3242-39.131.9300.04227.2253-38.421.9230.04526.8813-58.651.8410.00262.4474-60.311.8310.00163.6174…………


### 3.4. Analysis of Category Characteristics

Let the length of eigenvectors in feature setX=x1,x2,⋯,xn be l, and the number of eigenvectors be n; then, feature set X is divided into c fuzzy groups. In order to represent the degree to which each eigenvector belongs to c independent classes, the algorithm returns c clustering center matrix and membership matrix Uc×n=uij, where each element uij indicates that xi of X belongs to cj membership. The objective function is constructed as follows:(3)Ffcmu,c=∑j=1n∑i=1cuijmdij2xj,ci,s.t∑i=1cuij=1,∀j=1,⋯,n.Taking each row of feature set in Table1 as the corresponding feature vector, the feature vector matrix of M×N=4×8 can be constructed, and X is the input matrix of the fuzzy clustering algorithm. The number of clusters is set as c=4, the fuzzy weight parameter is m=2, the number of iterations is 500, and the iteration stop threshold parameter is ε=1.0×10−5 by four modes. The cluster center set C and membership matrix U of four kinds of eigenvectors can be obtained by operation.As shown in formula (3), C is the cluster center set, ci is the cluster center of class i, m is the fuzzy weight parameter, D is the initial number of clusters in fuzzy clustering, m=2 according to experience, and dlxi,ci=∑inxi−ci2 is the Euclidean distance from feature vector xi to class center ci. Through iterative calculation, the data is close to the optimal clustering center C=ci.(4)C=3.78530.49241.25821.36136.30720.56611.37131.84263.65190.63681.09111.19984.87300.66891.34611.3791,U=0.98610.99940.00650.02960.00010.01270.00740.00900.00000.00180.00000.97130.89360.00020.01150.00970.00070.00010.00820.00020.00320.01720.00000.00600.00250.98020.99990.00740.00030.02680.07260.99970.97070.98340.00200.0000.From the data in the membership matrixU, we can find the row number to which the eigenvector belongs, that is, the membership degree of the class label. The row code corresponding to the maximum value in each column of data is the category number to which the eigenvector belongs. According to matrix U, the maximum data values calculated in columns 1~2 appear in the first row, so the characteristic data in columns 1~2 in Table 1 are classified into the same category. Approximately, according to the maximum membership value of the characteristic data in columns 3~4, 5~6, and 7~8, the state class label can be determined.It can be seen that the above-extracted feature data has obvious category features, which can correctly distinguish the four mode state features of signal source, but the resolution is not clear enough, and more detailed classifiers need to be built to distinguish.


### 3.5. The Establishment of Classifier

The traditional FCA is influenced by kernel function, and fuzzy selection is not of probability significance. At the same time, it is hard to classify nonlinear samples. Therefore, for known feature setxii=1n, category tag sequence is lableii=1n, i∈1,2,⋯,n, and ωi is the weight parameter. The function of the soft interval classification optimization model of FCA is improved and constructed as follows:(5)minω,b12ω2+K∑i=1nℓlogyiωTxi+b−1.As shown in formula (5), where K>0 is a constant and ℓlogyiωTxi+b−1 is the logistic loss function, that is, ℓlogz=log1+exp−z. The introduction of the pair rate regression model has the advantage of output prediction marker and probability parameter, which is suitable for multiclassification problems in application.Therefore, it is equivalent to using the maximum likelihood method to solve the inner layer maximization in the Bayesian framework [19], as shown in(6)maxℓω,b=∑i=1nlogpyixi;ω,b.Letβ=ω,b, x⌢=x;1, formula (7) can be transformed into(7)minℓβ=∑i=1n−yiβTx⌢i+log1+eβTx⌢i.As shown in formula (7), it is a high-order continuous convex function of A. according to the Newton iterative method in convex optimization theory, the optimal solution β∗=argminβℓβ can be obtained, and then, the weights ω and b can be obtained. Furthermore, it is assumed that Pyiω follows Bernoulli distribution and ω follows Gauss prior distribution. In order to eliminate irrelevant noise component σ2, the likelihood estimation probability model is obtained, as shown in(8)Pyiω=∏i=1nσyxi;ωxi1−σyxi;ω1−xi.In order to reduce the learning degree of the model, a penalty parameterℝωi0,μi−1 is introduced for each weight to get the probability evaluation model.(9)Pωμ=∏i=0nℝωi0,μi−1.As shown in formula (9), μ=μ0,μ1,⋯,μnT is a sparse probability model constructed by n+1 dimensional hyperparameters. It can be seen that the improved FCA based on pairwise regression fusion has the advantage of estimating the output prediction marker and adding probability parameters at the same time.


### 3.6. Optimization of Classification Processing

Because the soft interval FCA classification model has high prediction cost, in order to reduce the overfitting risk of model training error, The SA algorithm is used to accelerate the convergence of FCA training process. The SA algorithm is easy to global optimization and to avoid local extremum, so it can achieve better global convergence.Step 1. Initialization program, in which the initial temperatureT0, initial solution vector ψ0, step size L0=1, solution space vector dimension 3, the number of iterations 200, and termination threshold ξ=1.0×10−3.Step 2. The FCA model is used to calculate the initial solutionωi, and the SA algorithm is called with ωi as the initial point to generate a new point ωi′=ωi+γL0 randomly, where γ∈0,1 is a random number and Pωiμ is calculated.Step 3. The fitnessφω,b is calculated by K-fold cross-validation, where the fold number K is set to 5, and the global optimal value is updated according to φω,b. If φω,b≤0, it has got the new solution is accepted. Otherwise, the new solution is accepted by probability p=exp−Δφω,b2i/T, where T is the current temperature control parameter. The annealing function is used to control the iteration speed Tn=αTn−1, and the annealing smoothing coefficient is set to α=0.8.Step 4. When the iteration satisfies criterionT<ξ or the number of iterations reaches, the program exits, and the fine classification of samples is obtained.


### 3.7. Train and Test

In order to verify the resolution and generalization ability of the above models, the “reserve method” is also used to identify and verify the samples collected from different mode signal sources. Each mode state is sampled in 30 groups, and a total of 120 groups of experimental data are collected, and their characteristic parameters are extracted. In order to prevent data overtraining, 100 groups of each group of signal feature set are randomly selected for model training. Using 20 groups of feature samples as verification set for classifier test, the iteration number is 500, then stop thresholdε=1.0×10−5, and model classification calculation to get confusion state matrix; according to the test results of confusion state matrix, the correct samples can be classified as table diagonal cumulative sum, a total of 18; the overall correct recognition rate can reach 90%, which can effectively identify the signals of four working modes.Among them, the mode signal with identification error is DCS. The reason is that the spectrum characteristics of the signal source in this frequency band are close to that of the adjacent signal frequency band. When the signal strength of the detection signal source is close to that of the adjacent frequency band, there will be misjudgment, but the misjudgment rate is low. The confusion state matrix calculated from the model classification is shown in Table2, which shows that the system can get enough accurate results for the specified signal source signal pattern identification.Table 2
Disturbance source signal characteristic identification results.
Test statusThe identification experiment resultsCDMADCSPHSMobileAccuracy rate (%)CDMA5000100DCS040190PHS0050100Mobile001490


## 3.1. Building Models

In order to solve the problem of large signal acquisition error caused by the multipath effect of indoor signal source, firstly, the spectrum of the radio signal is collected and processed on the indoor platform, and then, the random signal is input to the system to identify the large-scale signal. At last, the construction of signal acquisition and identification system is shown in Figure1.Figure 1
Structure diagram of signal acquisition and identification system.The known communication signal source is collected and input into the spectrum identification system by the SA44B spectrum analyzer. According to the four common communication modes of GSM, CDMA, DCS, and PHS, the corresponding spectrum features are extracted and stored in the database after signal detection feature extraction and identification. Then, after learning and training, a higher frequency spectrum feature recognition rate is achieved. In the online identification stage, the signal source with unknown working mode is input into the identification system, and the working mode of the output signal source is determined after identification.


## 3.2. Signal Acquisition and Processing

This section is to distinguish the source signal of the target for positioning; the four modes of spectrum signals of the signal source carried by the indoor sports platform are collected indoors. The spectrum of CDMA signal source has the characteristics of smooth waveform and obvious peak, as shown in Figure2. The spectrum of signal source in GSM mode is coarser than that in CDMA mode, but it has clear peak, as shown in Figure 3. The spectrum of DCS working mode is the most disordered among the four kinds of spectrum collected, as shown in Figure 4. In PHS mode, the peak to peak distance of communication signal spectrum is larger, so the characteristics are more obvious, as shown in Figure 5.Figure 2
Spectrum of CDMA mode.Figure 3
Spectrum of GSM mode.Figure 4
Spectrum of DCS mode.Figure 5
Spectrum of PHS mode.Obviously, the characteristics of the four modes of spectrum signals are complex, the optimal characteristic parameters of the spectrum are not easy to select, and there are many clutters in the spectrum, which affect the indoor positioning of the signal source.In order to prevent the edge information from blurring, a nonlinear filter is introduced, in which the median filter can keep the sharp edge well. The signal processed by median filter still contains noise, and the wavelet denoising method is better for reducing Gaussian white noise in radio wave signal. The wavelet denoising method [12] includes the modulus extremum denoising method, wavelet correlation denoising method, and wavelet threshold denoising method. The traditional wavelet hard thresholding function [13] and soft thresholding method [14] models are shown in(1)W⌢j,k=sgnWj,k·Wj,k−λ,Wj,k≥λ,0,Wj,k<λ,W⌢j,k=Wj,k,Wj,k≥λ,0,Wj,k<λ.The traditional wavelet hard thresholding function is easy to produce Gibbs oscillation [15], while the soft thresholding method is easy to produce “over smooth” distortion due to the constant difference of wavelet coefficients [16].(2)W⌢j,k=aWj,k+1−asgnWj,k·Wj,k−λexpWj,k−λδ,Wj,k≥λ,0,Wj,k<λ.As shown in formula (2), the threshold adjustment parameter a is a=1−e−αWj,k−λ2, W′ is the mean value of neighborhood wavelet coefficients, and δ is a positive number. When eWj,k⟶±λ is W⌢j,k⟶0, the threshold value Wj,k=±λ function is continuous, which can overcome the signal oscillation problem caused by the discontinuity of hard threshold; when Wj,k⟶±∞ is W⌢j,k/Wj,k⟶1, the new threshold function conforms to the characteristics of hard threshold function, which can overcome the hard threshold, so as to overcome the distortion caused by the constant difference of wavelet coefficients.When the threshold adjustment parametersα⟶0 and δ⟶0, the new wavelet threshold function is a soft threshold function; when the threshold adjustment parameters α⟶+∞ and δ⟶1, the new wavelet threshold function is a hard threshold function. By adjusting the parameters a and δ, we can obtain the signal processing results of different modes of the interference signal source. After denoising the spectrum of the four modes, we can normalize the Z-score [17], as shown in Figure 6.Figure 6
Spectrum normalization of four interference modes.


## 3.3. Spectrum Feature of Signal Extraction

The signal spectrum parameters are selected to extract the peak value, corresponding frequency, signal spectrum bandwidth, and kurtosis [18] of the signal transmitted by different mode signal sources as the characteristic parameters of the corresponding mode of the radio wave signal. The feature parameters of different modes of CDMA, GSM, DCS, and PHS of signal source are extracted step by step to construct the vector sequence required for the feature recognition model of radio wave signal. As shown in Table 1, the four modes of CDMA, GSM, DCS, and PHS are, respectively, marked as C=1~4.Table 1
Partial experimental characteristic data.
RSSIM (dBm)fM (103 MHz)fband (103 MHz)ψC-29.320.8670.02857.2311-30.510.8730.02756.5381-41.651.8060.06251.2672-40.371.8470.05850..3242-39.131.9300.04227.2253-38.421.9230.04526.8813-58.651.8410.00262.4474-60.311.8310.00163.6174…………


## 3.4. Analysis of Category Characteristics

Let the length of eigenvectors in feature setX=x1,x2,⋯,xn be l, and the number of eigenvectors be n; then, feature set X is divided into c fuzzy groups. In order to represent the degree to which each eigenvector belongs to c independent classes, the algorithm returns c clustering center matrix and membership matrix Uc×n=uij, where each element uij indicates that xi of X belongs to cj membership. The objective function is constructed as follows:(3)Ffcmu,c=∑j=1n∑i=1cuijmdij2xj,ci,s.t∑i=1cuij=1,∀j=1,⋯,n.Taking each row of feature set in Table1 as the corresponding feature vector, the feature vector matrix of M×N=4×8 can be constructed, and X is the input matrix of the fuzzy clustering algorithm. The number of clusters is set as c=4, the fuzzy weight parameter is m=2, the number of iterations is 500, and the iteration stop threshold parameter is ε=1.0×10−5 by four modes. The cluster center set C and membership matrix U of four kinds of eigenvectors can be obtained by operation.As shown in formula (3), C is the cluster center set, ci is the cluster center of class i, m is the fuzzy weight parameter, D is the initial number of clusters in fuzzy clustering, m=2 according to experience, and dlxi,ci=∑inxi−ci2 is the Euclidean distance from feature vector xi to class center ci. Through iterative calculation, the data is close to the optimal clustering center C=ci.(4)C=3.78530.49241.25821.36136.30720.56611.37131.84263.65190.63681.09111.19984.87300.66891.34611.3791,U=0.98610.99940.00650.02960.00010.01270.00740.00900.00000.00180.00000.97130.89360.00020.01150.00970.00070.00010.00820.00020.00320.01720.00000.00600.00250.98020.99990.00740.00030.02680.07260.99970.97070.98340.00200.0000.From the data in the membership matrixU, we can find the row number to which the eigenvector belongs, that is, the membership degree of the class label. The row code corresponding to the maximum value in each column of data is the category number to which the eigenvector belongs. According to matrix U, the maximum data values calculated in columns 1~2 appear in the first row, so the characteristic data in columns 1~2 in Table 1 are classified into the same category. Approximately, according to the maximum membership value of the characteristic data in columns 3~4, 5~6, and 7~8, the state class label can be determined.It can be seen that the above-extracted feature data has obvious category features, which can correctly distinguish the four mode state features of signal source, but the resolution is not clear enough, and more detailed classifiers need to be built to distinguish.


## 3.5. The Establishment of Classifier

The traditional FCA is influenced by kernel function, and fuzzy selection is not of probability significance. At the same time, it is hard to classify nonlinear samples. Therefore, for known feature setxii=1n, category tag sequence is lableii=1n, i∈1,2,⋯,n, and ωi is the weight parameter. The function of the soft interval classification optimization model of FCA is improved and constructed as follows:(5)minω,b12ω2+K∑i=1nℓlogyiωTxi+b−1.As shown in formula (5), where K>0 is a constant and ℓlogyiωTxi+b−1 is the logistic loss function, that is, ℓlogz=log1+exp−z. The introduction of the pair rate regression model has the advantage of output prediction marker and probability parameter, which is suitable for multiclassification problems in application.Therefore, it is equivalent to using the maximum likelihood method to solve the inner layer maximization in the Bayesian framework [19], as shown in(6)maxℓω,b=∑i=1nlogpyixi;ω,b.Letβ=ω,b, x⌢=x;1, formula (7) can be transformed into(7)minℓβ=∑i=1n−yiβTx⌢i+log1+eβTx⌢i.As shown in formula (7), it is a high-order continuous convex function of A. according to the Newton iterative method in convex optimization theory, the optimal solution β∗=argminβℓβ can be obtained, and then, the weights ω and b can be obtained. Furthermore, it is assumed that Pyiω follows Bernoulli distribution and ω follows Gauss prior distribution. In order to eliminate irrelevant noise component σ2, the likelihood estimation probability model is obtained, as shown in(8)Pyiω=∏i=1nσyxi;ωxi1−σyxi;ω1−xi.In order to reduce the learning degree of the model, a penalty parameterℝωi0,μi−1 is introduced for each weight to get the probability evaluation model.(9)Pωμ=∏i=0nℝωi0,μi−1.As shown in formula (9), μ=μ0,μ1,⋯,μnT is a sparse probability model constructed by n+1 dimensional hyperparameters. It can be seen that the improved FCA based on pairwise regression fusion has the advantage of estimating the output prediction marker and adding probability parameters at the same time.


## 3.6. Optimization of Classification Processing

Because the soft interval FCA classification model has high prediction cost, in order to reduce the overfitting risk of model training error, The SA algorithm is used to accelerate the convergence of FCA training process. The SA algorithm is easy to global optimization and to avoid local extremum, so it can achieve better global convergence.Step 1. Initialization program, in which the initial temperatureT0, initial solution vector ψ0, step size L0=1, solution space vector dimension 3, the number of iterations 200, and termination threshold ξ=1.0×10−3.Step 2. The FCA model is used to calculate the initial solutionωi, and the SA algorithm is called with ωi as the initial point to generate a new point ωi′=ωi+γL0 randomly, where γ∈0,1 is a random number and Pωiμ is calculated.Step 3. The fitnessφω,b is calculated by K-fold cross-validation, where the fold number K is set to 5, and the global optimal value is updated according to φω,b. If φω,b≤0, it has got the new solution is accepted. Otherwise, the new solution is accepted by probability p=exp−Δφω,b2i/T, where T is the current temperature control parameter. The annealing function is used to control the iteration speed Tn=αTn−1, and the annealing smoothing coefficient is set to α=0.8.Step 4. When the iteration satisfies criterionT<ξ or the number of iterations reaches, the program exits, and the fine classification of samples is obtained.


## 3.7. Train and Test

In order to verify the resolution and generalization ability of the above models, the “reserve method” is also used to identify and verify the samples collected from different mode signal sources. Each mode state is sampled in 30 groups, and a total of 120 groups of experimental data are collected, and their characteristic parameters are extracted. In order to prevent data overtraining, 100 groups of each group of signal feature set are randomly selected for model training. Using 20 groups of feature samples as verification set for classifier test, the iteration number is 500, then stop thresholdε=1.0×10−5, and model classification calculation to get confusion state matrix; according to the test results of confusion state matrix, the correct samples can be classified as table diagonal cumulative sum, a total of 18; the overall correct recognition rate can reach 90%, which can effectively identify the signals of four working modes.Among them, the mode signal with identification error is DCS. The reason is that the spectrum characteristics of the signal source in this frequency band are close to that of the adjacent signal frequency band. When the signal strength of the detection signal source is close to that of the adjacent frequency band, there will be misjudgment, but the misjudgment rate is low. The confusion state matrix calculated from the model classification is shown in Table2, which shows that the system can get enough accurate results for the specified signal source signal pattern identification.Table 2
Disturbance source signal characteristic identification results.
Test statusThe identification experiment resultsCDMADCSPHSMobileAccuracy rate (%)CDMA5000100DCS040190PHS0050100Mobile001490


## 4. Experimental Results and Analysis

### 4.1. The Comparison of Identification Rate

Through the SA44B spectrum receiver, the four modes of spectrum signals of the signal source carried by the indoor sports platform are collected indoors. Quantitative test experiments are carried out to verify the accuracy of feature selection of the SEEM algorithm and SA-FCA. The recognition rate is measured by extracting the percentage of correct samples in the unit number of features. The number of samples collected per unit is set to 200, and the proportion of correct eigenvalues and total eigenvalues is calculated, as shown in Figure7. When the number of spectrum increases to 300, the highest recognition rate can reach 93.1%. Then, the feature selection ability of the SEEM algorithm decreases with the increase of the number of features. When the number of features increases to 400, extracted by SA-FCA, the maximum feature recognition rate reaches 96%. From the comparison chart of recognition rate, it can be seen that the larger the number of selected features is not the better, and the number of features shows a downward trend when the number of features reaches 400.Figure 7
Comparison of SEEM and FCA identification rates.


### 4.2. Location Accuracy Analysis of Optimal Features

Furthermore, the optimized spectrum eigenvalues of the signal source are used for positioning experiments to verify the influence of the optimization algorithm on the positioning accuracy. Traditional NN algorithm, KNN algorithm [20], SEEM algorithm, and improved FCA algorithm are compared in indoor positioning experiments, as shown in Figure 8. The improved FCA feature selection method can achieve 78% in the positioning error of 1.3 m, which is 13% higher than the SEEM algorithm. The improved FCA feature selection method can achieve 83% in the positioning error of 1.5 m, which is 8% higher than the SEEM algorithm.Figure 8
Cumulative error distribution of SEEM algorithm and FCA.


## 4.1. The Comparison of Identification Rate

Through the SA44B spectrum receiver, the four modes of spectrum signals of the signal source carried by the indoor sports platform are collected indoors. Quantitative test experiments are carried out to verify the accuracy of feature selection of the SEEM algorithm and SA-FCA. The recognition rate is measured by extracting the percentage of correct samples in the unit number of features. The number of samples collected per unit is set to 200, and the proportion of correct eigenvalues and total eigenvalues is calculated, as shown in Figure7. When the number of spectrum increases to 300, the highest recognition rate can reach 93.1%. Then, the feature selection ability of the SEEM algorithm decreases with the increase of the number of features. When the number of features increases to 400, extracted by SA-FCA, the maximum feature recognition rate reaches 96%. From the comparison chart of recognition rate, it can be seen that the larger the number of selected features is not the better, and the number of features shows a downward trend when the number of features reaches 400.Figure 7
Comparison of SEEM and FCA identification rates.


## 4.2. Location Accuracy Analysis of Optimal Features

Furthermore, the optimized spectrum eigenvalues of the signal source are used for positioning experiments to verify the influence of the optimization algorithm on the positioning accuracy. Traditional NN algorithm, KNN algorithm [20], SEEM algorithm, and improved FCA algorithm are compared in indoor positioning experiments, as shown in Figure 8. The improved FCA feature selection method can achieve 78% in the positioning error of 1.3 m, which is 13% higher than the SEEM algorithm. The improved FCA feature selection method can achieve 83% in the positioning error of 1.5 m, which is 8% higher than the SEEM algorithm.Figure 8
Cumulative error distribution of SEEM algorithm and FCA.


## 5. Discussion

The classification analysis of eigenvalues extracted from radio wave signals is the basis of illegal signal identification. FCA model clustering is widely used as an unsupervised learning algorithm of the “hardening score” method. By adjusting the membership criterion, the individual feature vectors in the feature set are classified into subcategories of a cluster center one by one, and the similarity between different categories is the smallest. According to this, through the fuzzy clustering identification of the characteristics of CDMA, GSM, DCS, and PHS four different mode signals, a small amount of sample data is used to test the effectiveness of extracting the categories of the eigenvalues of radio wave signals.The FCA is a common classifier suitable for small sample sets and can also be extended to the field of multiclassification. Its core idea is to build a hyperplane model to classify the sample data in the feature space and remove the specific values. The difficulty is to separate the hyperplane close range samples with high confidence rate and maximum interval.


## 6. Conclusions and Future Work

Aiming at the problem of large positioning error caused by the working mode transformation of the signal source carried by the indoor sports platform, this paper first expounds the feature identification of the interference signal source from the new perspective of the spectrum detection and feature identification of the radio wave signal and studies the spectrum characteristics of the radio wave signal, after the denoising processing of the wavelet threshold method. Then, through the improved FCA model classification training and identification, the effect is better than the envelope method. Finally, the spectrum eigenvalues extracted from the identification are used for positioning. Indoor positioning experiments are carried out by traditional NN algorithm, KNN algorithm, seem algorithm, and improved FCA. The recognition rate of the collected samples is relatively ideal, and the positioning effect based on the feature data is better. In addition, there are other types of signal sources to be further studied, which require a large number of collected data to enrich the identification database, which need to be improved in the follow-up experiments.



---

*Source: 7913666-2021-11-16.xml*