A Pulse Signal Characteristic Recognition Algorithm Based on Multifractal Dimension

In low SNR condition, it is difficult to identify the radio transient characteristics of the signals. To solve this problem, a new recognition algorithm based on multifractal dimension characteristics is proposed. In fractal theory, multifractal dimension is the most sophisticated characterize that can describe the similar characteristics of the signals. Therefore, multifractal dimension is used in this paper to extract the subtle features of different impulse noise signals, in order to achieve the purpose of the classification and identification of the radiation source.


Introduction
In recent years, individual radiation source identification [1,2] has been developing rapidly as a nascent technology, which is an important topic in the field of noncooperative communication.Because of the continuous development of antireconnaissance technology and antijamming technology, we has found it difficult to get useful information through traditional intelligence information.The signal characteristics and individual subtle characteristics of individual radiation sources contain potential target information that is very meaningful.We can achieve the reconnaissance and surveillance of the goal by utilizing the subtle characteristics of the individual radiation sources.It is a new field with great research value.With the increasing complexity of the communication environment and the improvement of the information technology, the complexity of communication signals modulation types is becoming increasingly huge, while the differences between individuals radiation source are gradually decreasing.Therefore, the signal processing technology in recent years has become a hot issue in signal analysis areas and an integral part of signal recognition field.As the recognition results of modulation signals affect directly the restoration of the original signals, the automatic identification of communication signal plays a crucial role in the field of communication investigation [3], signal authentication [4], spectrum management [5], radio monitoring [6], electronic warfare [7], interference identification [8], and so forth.It is an important technology in radio spectrum management and intelligent signal processing.The basic task of communication signal recognition is to identify the modulation of the signals and estimate the modulation parameters and other basic information under the noisy environment or the multiple signal transmission conditions, finally, providing adequate basis for further processing and analyzing the signals.
A signal estimation and recognition method under white noise environment is put forward in literature [9], which obtains a good recognition performance.Therefore, in this paper, multifractal dimension is used to extract the different distribution of additional subtle impulse noise characteristics of radiation source, so as to achieve the purpose of classification and recognition.

Flowchart of Identification System
For the communication individual radiation source, even two models from the same manufacturer, the same individual categories, and the same working conditions with the same operating parameters, there are also subtle differences between individual radiation signals, due to the individual differences in their internal components.1.
First, the noises with different pulse widths and different pulse positions are added to different radiation sources which transmit the same signals.Then, the multifractal theory is used to extract the subtle characteristics of different signals and to establish the database of multiple fractal dimension features.Finally, multifractal dimension characteristics of the unknown signals which would be identified are obtained, and then gray relation value is calculated with the signal features in the database.The type of the radio in the database which has the maximum output relation degree is the type of the radiation source to be identified; thus the radiation source classification is realized.

Multifractal Dimension Theory.
Fractal dimension is an important indicator which is used to measure signals irregularity and explore the complexity of things.At present, it has been widely applied in image analysis [10], vibration signal fault diagnosis [11], and dynamics analysis [12].Onedimensional fractal dimension, box-counting dimension [13], can be sufficient for the description of the features of a simple signal.But, for complex signals, multifractal parameters are needed to describe the local characteristics of different levels.Therefore, the concept of multifractal dimension [14] is introduced in order to further understand the local characteristic of the signal and extract the signal features of different levels complexity.The proposed multifractal theory makes the subtle feature extraction of individual radio source possible.
Multifractal dimension is used to describe signal characteristics of different levels, and it is also used to discuss the probability distribution of parameters.It divides the research object (set its linear degree to 1) into  small areas.Suppose the line degree of areas  is   , and the density distribution function   of area  can be described by different scaling index   as Noninteger   can commonly be called a singular index, and its value is relevant to its region.
In order to get the distribution characteristics of a series of subsets, function   () is defined as follows: ( It is the weighting summation of probability in each area.Further, the generalized fractal dimension   is defined as where   () reflects the different effects of   .The formula shows that, when  ≫ 1, large probability area plays the leading role in the summation formula ∑  =1    and, at the same time,   () and   reflect the properties of big probability areas.Conversely, when  ≪ 1,   () and   reflect the properties of small probability areas (sparse areas).In this way, different  values reflect the properties of different probability characteristic areas.Through the weighting summation processing, it divides the signal into many areas with different singular degrees, and, therefore, the internal structure of the signal is hierarchical.Calculate   through phase space reconstruction and extract the different probability characteristics of the signal; the multilevel feature extraction results can be finally got in this paper.

Gray Relation Theory
where   is defined as the distinguishing coefficient.
Then the grey relation is defined as Here, ( 0 ,   ) is the grey relation of  0 and   .The gray relational ( 0 ,   ) is often written as  0 and the correlative coefficient ( 0 (),   ()) of  is written as  0 ().
Thus the calculation procedure of grey correlation can be got as follows.
(5) Calculate the gray correlation: Thus the gray relation values of different features are obtained, and the category of the signals to be recognized can be judged by comparing the value of the gray relation degree.

Simulation Results and Analysis
In order to verify the property of multifractal dimension, we randomly generate 1000 FSK signal samples with MATLAB.In Monte Carlo simulation experiment, two kinds of color noises are added onto FSK signals.Taking the two color noises as example, we suppose it to be the output of the internal device subtle characteristics, figure out the multifractal dimension curves, and verify the performance of the proposed algorithm.The  is varied from −10 to 10, which is used to get 21 multifractal dimensions and reflect the 21 levels of different subtle features of the signal.As we can see from Figure 2, the abscissa represents the dimension of reconstructed phase space, denoted by ln , and the ordinate represents the natural logarithm of ∑ 6 =1     , denoted by ln   .Simulation results are shown in Figure 2.
The simulation results show the following.
(1) The smaller the reconstructed phase space dimension is, the smaller the difference of different probability levels is.And it does not give further differences of the signals.When reconstructed phase space is larger, the signal is divided into different areas to study, and the probability of different characteristics areas is calculated, so more details of the signals are detected, and the differences are bigger.
(2) According to the fact that different  values reflect the characteristics of different probability characteristic areas, the  values of the curves take −10-10 from bottom to top, respectively.Thus it can be seen from the distribution of the curves that, when  is bigger, the differences between the signals are much more obvious.
As the distributions of different color noises are different, multifractal dimension can be used to detect the complexity changes of FSK signals, whose multifractal curves also show obvious differences.Then gray relation theory is used to identify the two kinds of color noise, and the results show the recognition rate of one hundred percent.
Taking FSK signal of communication individual radio source with the same distribution of impulse noise as an example, multifractal characteristics of FSK signal sequences with the same pulse distribution noise are discussed, but the subtle differences are that the pulse position and pulse width are not the same.Then the recognition performance of this algorithm is verified.The signal sequences consisted of 1024 points, where one signal produces an impulse noise at the place of 401-403 points, and the other signal sequences produce an impulse noise at the place of 701-703 points.Multifractal dimension curves are calculated and the simulation results are shown in Figure 3.The same signal sequences consist of 1024 points as described above, but, at the same place, the width of one impulse noise is 1 and the width of another impulse noise is 3.And the multifractal dimension curves are calculated and the simulation results are shown in Figure 4.
Conclusions can be drawn from simulation results in Figures 3 and 4 as follows.
(1) When the pulse positions of the noises are different, the increase of reconstructed phase space dimension will fade the characteristic of the noise.Only under lower phase space dimensions are the multifractal dimensions of signals different to some extent.(2) When  value is small, the differences of multifractal dimensions are also small, which makes it hard to classify the signals.Therefore, the differences of large probability area signals are more obvious.(3) The distribution curves show that the influence of the impulse noises with different widths is bigger than that of the impulse noises with different positions.Since the additive noises are randomly generated, the noises will be different at each time; however, the pulse position and pulse width of the noises are the same.Therefore, it has better universality and application under the condition of the same pulse position and pulse width.As described before, 1000 samples are randomly generated, and Monte Carlo simulation experiment is carried out to calculate the recognition effect.Due to the randomness of the noise, there are some differences on noise distribution.Therefore, there are some differences of the recognition rates at each time.In this experiment, we carry out the simulation for 6 times and randomly generate 2000 distribution sequences of the two kinds of noises every time.Then we use gray relation theory to classify and identify the signals.The entropy features of the noises are extracted under the same background, and then the signals could be recognized.The recognition rates of the two kinds of methods are shown in Tables 1 and 2, and the comparison of the simulation time is shown in Table 3.The comparison of Tables 1, 2, and 3 depicts that the results of using entropy feature method to measure the noises subtle characteristics are rougher than multifractal dimension.Therefore, the recognition rate is relatively lower compared with multifractal characteristics.Meanwhile, under the condition of randomly generated noise, the recognition rate of entropy feature is not stable, and the fluctuated variance is bigger.For the computational complexity, the algorithm of multifractal dimension is more complex than the algorithm of the entropy value.Furthermore, the computational complexity of multifractal dimensions is decided by the calculated dimension that we choose.We need to compute more dimensions if we want to achieve more accurate feature extraction.As a result, the amount of calculation will be greater.On the contrary, the amount of calculation will be less if we calculate fewer dimensions.Considering the better recognition rate and stability, multifractal characteristics have better application value in communication individual recognition.

Conclusion
It is difficult to extract the subtle features of the radio station under low SNR.In order to solve this problem, the paper proposes a subtle feature extraction algorithm based on multifractal dimension.

Figure 2 :
Figure 2: Multifractal curves of FSK signals mixed with different color noises.

FSK signal + pulse noise 1 FSK signal + pulse noise 2 − 20 Figure 3 :
Figure 3: Multifractal curves when the pulse positions of the noises are different.

Figure 4 :
Figure 4: Multifractal curves when the pulse widths of the noises are different.
First, for the same FSK signal of radio transmitter, different color noises are added to the signals to verify the feasibility of subtle features classification based on multifractal dimension.Then the FSK signals with different pulse widths and different pulse positions are used to imitate the different additional small internal radio noises, and multifractal dimension is used to classify the signals and verify the effect of recognition algorithm.Simulation results show that the multifractal dimension can characterize the small changes of signals in different scales to achieve the recognition of different radiation sources.It has great value in the project.
Therefore, it is possible to use the subtle characteristics to identify the signals of communication individual radiation source.Nature

Table 1 :
Recognition rates of two kinds of feature extraction methods under the condition of different position pulse.

Table 2 :
Recognition rates of two kinds of feature extraction methods under the condition of different pulse width.

Table 3 :
Recognition time of two feature extraction methods.