Joint Blind Identification of the Number of Transmit Antennas and MIMO Schemes Using Gerschgorin Radii and FNN
Abstract
Blind enumeration of the number of transmit antennas and blind identification of multipleinput multipleoutput (MIMO) schemes are two pivotal steps in MIMO signal identification for both military and commercial applications. Conventional approaches treat them as two independent problems, namely the source number enumeration and the presence detection of spacetime redundancy, respectively. In this paper, we develop a joint blind identification algorithm to determine the number of transmit antennas and MIMO scheme simultaneously. By restructuring the received signals, we derive three quantified features based on the signal subspace to determine the number of transmit antennas and identify spacetime redundancy. Then, a Gerschgorin radiibased method and a feedforward neural network are employed to calculate these three features, and a minimal weighted norm1 distance metric is utilized for decision making. In particular, our approach can identify additional MIMO schemes, which most previous works have not considered, and is compatible with both singlecarrier and orthogonal frequency division multiplexing (OFDM) systems. Simulation results verify the viability of our proposed approach for singlecarrier and OFDM systems and demonstrate its favorable identification performance for a short observation period with acceptable complexity.
I Introduction
Blind identification of information signals’ parameters of a transmitter from received signals has important applications both in military and civilian communication systems. In the context of military applications, this parametric knowledge can help an attacker carry out electronic warfare operations, such as surveillance and jamming signal selection. Moreover, blind identification has found wide applications in civilian reconfigurable systems including softwaredefined and cognitive radios[1]. Multipleinput multipleoutput (MIMO) and orthogonal frequency division multiplexing (OFDM) technologies are adopted in cellular and WiFi standards because of their excellent performance. Different from the identification of singleantenna systems, the blind identification of MIMO or MIMOOFDM signals requires the enumeration of the number of transmit antennas[2, 3, 4, 5, 6] and identification of MIMO schemes[7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23] for singlecarrier or OFDM systems.
The identification of MIMO schemes is the process of classifying the spatial multiplexing (SM) or transmit diversity (TD) codes, namely spacetime block codes (STBC), which utilize spacetime redundancy to reduce the error rate. Previous works on the identification of MIMO schemes include [7, 8, 9, 10, 11, 12, 13, 14, 15, 16] for singlecarrier systems and [17, 18, 19, 20, 21, 22, 23] for OFDM systems. Regarding the identification of MIMO schemes for singlecarrier systems, previous works follow either likelihoodbased [7] or featurebased [8, 9, 10, 11, 12, 13, 14, 15, 16] methods. The former relies on the likelihood function of the received signals to quantify the spacetime redundancy and classify different STBCs. The latter detects the presence of the spacetime redundancy at some specific timelag locations based on the features of signal statistics or cyclic statistics; however, it can only identify a small number of STBC types due to identical features for several SFBCs. Although [16] can identify 11 types of STBCs utilizing the feature of secondorder cyclostationary statistics, it requires the number of transmit antennas and channel coefficients as apriori information.
As for OFDM systems, there are two main approaches to combine the TD codes with OFDM signals. The first approach is STBCOFDM, where the diversity STBC is implemented over consecutive OFDM symbol intervals. STBCOFDM has been adopted in several indoor MIMO wireless standards, such as WiFi[24, 25], owing to its excellent performance. Another TD coding scheme used in OFDM systems is spacefrequency block coding (SFBC)OFDM, where the diversity SFBC is employed over consecutive subcarriers of an OFDM symbol. Several cellular wireless standards supporting high mobility, such as LTE[26] and WiMAX [27], favor SFBCOFDM over STBCOFDM because of its superior performance in a high mobility environment [28]. The identification approaches for singlecarrier systems fail to identify MIMO schemes of OFDM systems under frequencyselective fading channels due to multipath effects. Previous works on the identification of STBCOFDM systems include [17, 18, 19], which detect the presence of spacetime redundancy based on the peaks of the crosscorrelation functions between two receive antennas in the timedomain. Specifically, [17, 18] use different crosscorrelation functions, while [19] employs a cyclic crosscorrelation function with a specific timelag during adjacent OFDM symbols to detect the spacetime redundancy. However, the approaches for identifying STBCOFDM signals cannot be directly applied to SFBCOFDM signals since the peaks of the crosscorrelation functions between two adjacent OFDM symbols are difficult to detect. The previous works on the identification of SFBCOFDM systems include [20, 21, 22, 23]. In [20], the idea of detecting the peak of the crosscorrelation function between two receive antennas is extended to identify SFBCOFDM signals, which uses a specific timelag during the same OFDM symbol. In [21], we utilize a crosscorrelation function between two receive antennas at adjacent OFDM subcarriers to improve the identification performance by detecting both the space and frequency redundancy. In [22, 23], we use the random matrix theory to identify 5 types of SFBCs by detecting the spacefrequency redundancy at adjacent OFDM subcarriers. However, most previous works can only discriminate between a few MIMO schemes since they only consider detecting the presence of the redundancy. Specifically, a widelyused TD code, namely the frequency switched transmit diversity (FSTD)[26, 27], and several nonorthogonal STBCs/SFBCs can not be discriminated by the previous nonlikelihoodbased methods. Although the likelihoodbased method [7] can identify more STBCs, it does not work in a frequencyselective fading environment.
In the existing literature, the enumeration of the number of transmit antennas and identification of MIMO schemes are handled as two independent problems. The enumeration problem of the number of transmit antennas is formulated as the enumeration of independent channel pathways between transmit antennas and receive antennas in general. Previous works on the identification of the number of transmit antennas mainly fall into two classes, namely second [2, 3, 4, 5] and higherorder statisticsbased methods [6]. Basically, the secondorder statisticsbased methods analyze eigenvalues or eigenvectors of the covariance matrix of the received signals to determine the number of transmit antennas by distinguishing between the signal and noise subspaces. These methods also fall under two categories, theoreticcriteriabased algorithms[2] and hypothesistestingbased algorithms [3, 4, 5]. Methods in the first category determine the number of transmit antennas by minimizing the KullbackLeibler distance metric. Reference [2] introduces two classical calculations of the KullbackLeibler metric, the Akaike information criterion (AIC) and the minimum description length (MDL), for the enumeration of the number of transmit antennas. Methods in the second category transform the problem into a detection problem, which compares an elaborately constructed statistic with a threshold. Furthermore, for the class of higherorder statisticsbased methods, the sole existing algorithm [6] constructs a fourthorder decision statistic of the received signals with only one receive antenna by using the feature of timevarying block fading channels. To the best of our knowledge, no method exists in the literature for the joint blind identification of the number of transmit antennas and MIMO schemes. It is our main goal in this paper to fill this research gap.
As mentioned earlier, the quantification of the spacetime/frequency redundancy can classify more STBCs/SFBCs since the redundancy of some STBCs/SFBCs is in the same location. On the other hand, this quantification is also needed in the enumeration of the number of transmit antennas. In this paper, a quantified featurebased joint blind identification algorithm of the number of transmit antennas and MIMO schemes is proposed. Three different quantified features for the number of transmit antennas and redundancy are derived from the eigenvalue analysis of the covariance matrix of the received signals at adjacent symbols or OFDM symbols/subcarriers with multiple receive antennas. A Gerschgorin radiibased method and a feedforward neural network (FNN) are applied to calculate these features, and a minimal weighted norm1 distance metric is proposed to determine the number of transmit antennas and MIMO schemes. The proposed algorithm does not require apriori knowledge of the signal parameters, such as channel coefficients, modulation type or noise power.
The main contributions of this paper are the following:

The proposed algorithm jointly identifies the number of transmit antennas and MIMO schemes, which has not been considered in the previous works.

The scenarios of singlecarrier and OFDM, including STBCOFDM and SFBCOFDM are all investigated in this paper, unlike previous works.

Unlike the existing algorithms, more STBC/SFBC types, such as the orthogonal STBCs/SFBCs (OSBC) with the same rate, FSTD, quasiorthogonal STBC/SFBC (QOSBC) and nonorthogonal STBCs/SFBCs (SBC), are identified by the proposed algorithm thanks to the analysis of quantified features.

A Gerschgorin radiibased method and an FNN are efficiently combined to calculate the quantified features. Furthermore, we extend the investigation to OFDM systems.

Simulation results are presented to demonstrate the viability of the proposed algorithm both in singlecarrier and OFDM systems, with different system parameters.
This paper is organized as follows. In Section II, the system model is introduced. Then, Section III derives the three quantified features. The proposed algorithm is described in Section IV. The simulations results are presented in Section V. Finally, conclusions are drawn in Section VI.
Notation: The following notation is used throughout the paper. The superscripts , and denote the complex conjugate, transposition and conjugate transposition, respectively. represents the probability of the event . indicates the statistical expectation. and denote the real and imaginary parts, respectively. , and denote the identity matrix, zero vector and zero matrix, respectively. , and are the set of natural numbers, positive integers and complex numbers, respectively. The notation represents the symbol at the th transmit or receive antenna. The notation denotes the cardinality of the set . denotes the trace of a matrix. Conventionally, and denote the Euler constant and natural logarithm, respectively. Finally, denotes the complexity order.
Ii System Model
Iia Signal Model of MIMO SingleCarrier System
We consider a MIMO singlecarrier wireless communication system employing TD or SM with transmit antennas and receive antennas, as shown in Figure 1 (a). As a special case, singleantenna systems are also considered. The transmitted data symbols are drawn from an PSK (PhaseShiftKeying) or QAM (Quadrature Amplitude Modulation), , signal constellation. Subsequently, the modulated symbol stream is parsed into a data block of symbols, denoted by the vector (). A TD/SM encoder takes the row of an codeword matrix, denoted by , to span consecutive time slots and maps every column of the matrix into different transmit antennas. In this paper, the codewords include the singleantenna, Alamouti (AL), SM, 7 types of OSBC[29, 30, 31], one type of QOSBC[32], FSTD in LTE[26] and 3 types of SBC in WiMAX[27] (see Appendix A). Then, mapped signals are transmitted after the pulse shaping and carrier modulation operations.
The receiver is assumed to perfectly estimate the synchronization parameters.^{1}^{1}1Several algorithms can be applied to blindly estimate synchronization parameters, such as those in [33, 34] for the singlecarrier system and in [35, 36] for the OFDM system. The sensitivity analysis to model mismatches is provided in Section V. We construct an transmit vector and an receive vector, denoted by and , which represent the transmitted and received signals at the th () time slot, respectively. The channel is assumed to be flatfading and characterized by an matrix of Rayleigh fading coefficients, denoted by
(1) 
where represents the channel coefficient between the th transmit antenna and th receive antenna. The channel matrix is assumed to be of fullcolumn rank and the channel gains remain constant over the observation interval. Then, the th received signal is described by the following model
(2) 
where the vector represents a white Gaussian noise vector with zeromean and covariance . The first processed sample is assumed to be the start of a TD code block, which allows simplifications of the following mathematical expressions. However, extensions of the proposed methods can be easily obtained when this assumption does not hold.
IiB Signal Model of MIMOOFDM System
IiB1 STBCOFDM System
Consider a MIMOOFDM wireless communication system with subcarriers and a cyclic prefix (CP) of length , as shown in Figure 1 (b). Different from the singlecarrier system, the TD/SM encoder puts data blocks, denoted by , on consecutive subcarriers with the same operation as the singlecarrier system. At the receiver side, the received OFDM symbol is converted into a frequencydomain block via an point fast Fourier transform (FFT) after removing the CP. We construct an transmit vector and an receive vector, denoted by , and , which represent the transmitted and received signals at the th time slot and th () subcarrier, respectively. The channel is assumed to be frequencyselective fading and the th subchannel is characterized by an matrix of Rayleigh fading coefficients, denoted by . Then, the th received signal at the th subcarrier is described by the following model
(3) 
where the vector represents a frequencydomain white Gaussian noise vector at the th subcarrier.
IiB2 SFBCOFDM System
The system model of the SFBCOFDM system is similar to the STBCOFDM system, with the difference that the SFBC encoder takes the row of the codeword matrix to span consecutive subcarriers directly. The three mappings are shown in Figure 1 (c).
Iii Quantified Redundancy Features of Different Numbers of Transmit Antennas and MIMO Schemes
In this section, three quantified features are defined as discriminating features of different MIMO schemes, by considering the dimension of the subspace of restructured received signals. In addition, we summarize these features for different numbers of transmit antennas and MIMO schemes in Table I. Without loss of generality, we consider the singlecarrier system first and then extend the analysis to the OFDM system.
Iiia Number of Transmit Antennas Feature
Let us construct a timedomain receive window to observe the received signals at adjacent time slots. The window length is set to two since it is the finest granularity to quantify the features of different MIMO schemes, while a larger window length results in failure to distinguish many MIMO schemes. By using (2), the th received signal block inside the window is expressed as
(4) 
where the transmitted signal block of adjacent time slots is denoted by an matrix , and the noise vector is denoted by . The covariance matrix of the received block is
(5) 
where is the covariance matrix of the th transmitted block. If we denote the eigenvalues of by () in descending order, then we have
(6) 
Proposition 1: The smallest ordered eigenvalues of are all equal to , i.e., if the rank of is .
Proof: See Appendix B.
The corresponding eigenvectors of the eigenvalues form a basis for the signal subspace. Define the following subset of signal subspace eigenvalues of
(7) 
By sliding the window, the cardinality of the set with even subscript (, ) can be used as the first type of the quantified feature since the cardinality of for FSTD is equal to two. Therefore, we define this cardinality as the number of transmitantenna feature (NTAF), denoted by
(8) 
Note that the NTAF, , is the discriminating feature for different numbers of transmit antennas.
IiiB Number of Independent Complex Symbols Feature
Let us vectorize the th signal and noise block inside the window in (4) as follows
(9) 
Then, the th received vectorized block is rewritten as
(10) 
where the matrix is
(11) 
Like (5), the covariance matrix of the th received vectorized block is
(12) 
where is the covariance matrix of the th transmitted vectorized block. The eigenvalues of are
(13) 
Next, we define the notion of linearly independent random symbols. Let be the vector observations of the random variables . Then, we define as linearly independent random symbols if the equation
(14) 
can only be satisfied by for [37]. For example, assuming whose elements are complex random variables, then has five linearly independent complex random symbols, i.e., , , , , . Assuming , then has four linearly independent real random symbols, i.e., , , , .
Proposition 2: Assume that is the number of linearly independent complex random symbols of a transmitted vectorized block. The smallest ordered eigenvalues of are all equal to , i.e., .
Proof: See Appendix B.
Similarly, let us define a subset of the signal subspace eigenvalues of as
(15) 
By sliding the window, we can select the cardinalities of three sets, namely with subscripts , , and () as the second type of quantified features to identify more MIMO schemes in the MIMO scheme pool. Then, we define three linearly independent complex randomsymbol features (ICSF), which represent the cardinalities of the three sets, denoted by
(16a)  
(16b)  
(16c) 
IiiC Number of Independent Real Symbols Feature
By stacking the real and imaginary parts of the block inside the window in (9), we obtain a transmitted, a received stacked block and the noise block as
(17) 
Similarly, the th received stacked block is rewritten as
(18) 
where the matrix is given by
(19) 
and corresponds to the Kronecker product. Then, the covariance matrix of is
(20) 
where is the covariance matrix of the th transmitted stacked block. The eigenvalues of are
(21) 
Proposition 3: Assume that is the number of linearly independent real random symbols of a transmitted stacked block. The smallest ordered eigenvalues of are all equal to , i.e., .
Proof: See Appendix B.
Let us define the following subset of the signal subspace eigenvalues of as
(22) 
By sliding the window, the cardinalities of the sets with subscripts and () are chosen to be the third type of quantified features since the cardinalities of the other sets can not be used to identify several TD codes with different transmission rates. Then, we define two linearly independent real randomsymbol feature (IRSF) representing these cardinalities given by
(23a)  
(23b) 
Actually, the ICSF and IRSF can quantify the spacetime redundancy. For different MIMO schemes with the same number of transmit antennas, the transmitted signal block inside the window has more spacetime redundancies, or in other words, the block transmits some identical symbols due to the structure of the codeword matrix. Therefore, smaller ICSFs and IRSFs values are calculated at the receiver since these symbols are linearly dependent.
Table I shows all signal features using different numbers of transmit antennas and common MIMO schemes.
Number of  MIMO  Code  

Antennas  schemes  rate  
1  Single  1  1  2  2  2  4  4 
antenna  
2  1  2  4  4  4  4  4  
2  2  2  4  4  4  8  8  
3  1/2  3  4  4  4  8  8  
3  3/4  3  3  5  5  6  6  
3  3/4  3  5  3  3  6  6  
3  3/4  3  5  5  3  6  6  
3  2  3  6  6  6  8  8  
3  3  3  6  6  6  12  12  
4  1/2  4  4  4  4  8  8  
4  3/4  4  5  5  5  6  6  
4  1  4  8  6  8  8  8  
4  1  4  8  4  8  8  8  
4  1  4  4  4  4  4  4  
4  2  4  8  8  8  8  8  
4  3  4  8  8  8  12  12  
4  4  4  8  8  8  16  16 
IiiD Quantified Features in the OFDM System
IiiD1 STBCOFDM Case
Since STBCOFDM is a direct extension of STBC singlecarrier, we can use the th received signal block at the th subcarrier, denoted by , to derive the NTAF, ICSFs and IRSFs as we described previously.
IiiD2 SFBCOFDM Case
Regarding the SFBCOFDM system, we construct and slide a frequencydomain receive window to observe the received signals at adjacent subcarriers and calculate the quantified features. Thus, the th received block, vectorized block and stacked block in (4), (9) and (17) inside the window are, respectively, modified to
(24) 
Iv Proposed Blind Identification Algorithm
In this section, we use a Gerschgorin radiibased method and an FNN to calculate the quantified features and employ a minimal weighted norm1 distance metric to discriminate between these features. Different from the original Gerschgorin radiibased method in [38], the radii of the circles are compressed after a similarity transformation, and then, an FNN is used to calculate the quantified features. Additionally, extensions to the STBCOFDM and SFBCOFDM systems are proposed by combining the calculations from different subcarriers after an analysis of the channel response. Our motivation for modifying the original algorithm is that the algorithm can not identify ICSFs and IRSFs well since the Gerschgorin circles difference between the signal and noise subspace gets smaller, which causes more identification errors.
Iva Proposed Algorithm for SingleCarrier System
The covariance matrices estimators are given by
(25a)  
(25b)  
(25c)  
(25d)  
(25e)  
(25f) 
where is the number of symbols. For convenience, we employ a common notation, , to represent the estimated covariance matrices in (25). Assume that is a matrix. First, we partition the estimated covariance matrix as
(26) 
where the reduced covariance matrix is a leading principal submatrix of obtained by removing the last row and column of . Then, the reduced covariance matrix can be decomposed by its eigenstructure as follows
(27) 
where is a diagonal matrix constructed from the eigenvalues of as
(28) 
and is a unitary matrix formed by the corresponding eigenvectors as follows
(29) 
Then, we construct the following unitary transformation matrix
(30) 
where . The transformation of is
(31) 
where for . Thus, the radius of the th Gerschgorin circle is
(32) 
In order to scale the radii of Gerschgorin circles in proportion, we construct the diagonal matrix
(33) 
where is the mean of the eigenvalues given by
(34) 
The matrix can be similarly transformed into
(35) 
Practically, the centers of the Gerschgorin circles are fixed while their radii are relatively squeezed by the factor . The Gerschgorin circles of the noise subspace are squeezed more than those of the signals since the noise circles radii are generally smaller than . Then, the radii of the compressed Gerschgorin circles for , are denoted by
(36) 
After extracting the radii of the compressed Gerschgorin circles, the identification problem can be considered as a classification problem of the Gerschgorin circles which determines how many Gerschgorin circles are those of the signal subspace. However, the radii of the Gerschgorin circles of the signal and noise subspaces have a wide range under different conditions, including signaltonoise ratio (SNR), the number of processed samples and number of receive antennas, which results in the nonlinearity between the inputs and outputs of the learning system. An artificial neural network (ANN) is ideally suited for nonlinear and complex relationships, and does not impose any restrictions on the input variables, unlike other prediction techniques. Given that the FNN is a popular family of ANN owing to its simple structure, it can be used to develop high performance signal identification solutions [39, 40].
In this paper, we use an FNN, as shown in Figure 3, to determine the number of the signalsubspace Gerschgorin circles. The FNN includes an input layer, a hidden layer and an output layer. After the SNR is normalized, the LevenbergMarquardt backpropagation algorithm [41] is utilized to train the FNN by using test data. Then, the feature value is a fitting function of the radii of compressed Gerschgorin circles given by
(37) 
The quantity represents the calculated feature, , , , , or depending on the corresponding covariance matrix in Equations (25a)(25f). Since the sizes of the covariance matrices and the eigenvalues after decomposition have large differences for the three quantified features, the numbers and values (the distributions of values of the radii) of the FNN inputs are significantly different for different features. To enhance performance, we use three trained FNNs to determine , and , respectively.
Finally, the number of transmit antennas, , and the MIMO scheme, , are determined by a minimal weighted norm1 distance metric given by
(38) 
where the notation {CODE} refers to the set of all schemes listed in Table I. The reason for employing a norm1 distance metric is that it is more robust against outliers than other distance metrics [42]. The weighting coefficients are chosen to eliminate the bias caused by the features with larger values and equal to the least common multiple of the NTAF, sum of ICSFs and sum of IRSFs for singleantenna, which is equal to 24, divided by these values resulting in the weighting coefficients of 1, 6, 8, respectively.^{2}^{2}2 It is unfair for the coefficient of the NTAF if the MIMO scheme with the spacetime redundancy is chosen here, since the NTAF does not quantify the spacetime redundancy and has the same value for different MIMO schemes with the same number of transmit antennas. For clarity, the main steps of the proposed algorithm are summarized subsequently.
IvB Extension to OFDM Systems
For an OFDM system, each frequency subchannel can be reasonably assumed to be a quasistatic flatfading channel since the subchannel width is designed to be much less than the channel’s coherence bandwidth constraint. The frequency responses of adjacent subchannels are close to each other, especially when increasing the number of subchannels under a given total bandwidth constraint. Therefore, we can rewrite the subchannel frequency response at the th subcarrier as
(39) 
where the tiny increment, , has the form
(40) 
Figure 3 shows the frequency response for a frequencyselective fading channel which consists of 3 independent taps with an exponential power delay profile[18], . In this figure, the number of subchannels is set to 64, and the two yaxes represent the amplitude and phase responses, respectively. From the figure, we can reasonably assume that four consecutive subcarriers have similar responses, expressed as
(41) 
IvB1 STBCOFDM Case
For convenience, let us use a new variable to rearrange the subscript indices of subcarriers, denoted by , where represents the ceiling function. Assume that OFDM symbols are observed at the receiver. Signals at four consecutive subcarriers are serially incorporated into a data block and the th data block is denoted by
(42) 
According to the assumption of (41), the block can be restructured into a received block, vectorized block and stacked block via (4), (9) and (17), respectively, and be processed by the operations from (26) to (36) to compute the radii of the compressed Gerschgorin circles, denoted by . Different data blocks transmit over independent channels and have independent data so that we can combine results from different detectors of the data blocks with different , expressed as
(43) 
where is the number of detectors. Finally, the minimal weighted norm1 distance metric in (38) decides on the number of transmit antennas and the MIMO scheme using the trained FNNs.
IvB2 SFBCOFDM Case
The estimators of different covariance matrices need to be rewritten as follows
(44a)  
(44b)  
(44c)  
(44d)  
(44e)  
(44f) 
where is constrained by the conditions and . We use detectors to calculate the radii of the compressed Gerschgorin circles, for each adjacent subcarriers and combine detectors as
(45) 
Then, the radii are input into the trained FNNs and the features are computed. Finally, the number of transmit antennas and MIMO scheme are determined by the minimal weighted norm1 distance metric in (38).
V Simulation Results
Va Simulation Setup
System  /  CP  Modulation  MIMO scheme  Channel  
pool  
singlecarrier  8  256        4PSK  17 types  Flatfading with Rayleigh fading coefficients 
OFDM  8  100  128  32  10  4PSK  17 types  Frequencyselective fading consisting of 4 independent Rayleigh taps with power delay profile 
Monte Carlo simulations are conducted to evaluate the performance of the proposed algorithm. We consider three systems, namely, the singlecarrier, STBCOFDM, and SFBCOFDM systems. Unless otherwise mentioned, the default values of the number of receive antennas, observed samples and OFDM symbols, subcarriers, detectors, CP length, modulation type, MIMO scheme pool and channel parameters are as listed in Table II. The noise is assumed zeromean additive white Gaussian with variance . The total power of transmitted signals is constrained to regardless of the number of transmit antennas and the SNR is defined as [7]. The average probabilities of correct identification of the number of transmit antennas and MIMO scheme, respectively, are calculated as follows
(46) 
respectively, and used as performance measures. The MIMO scheme pool is defined to contain the 17 types of MIMO schemes listed in Table I. The simulation of each MIMO scheme was run for 1000 trials.
VB Simulation Results for SingleCarrier System
We used 2048 received signal samples with different SNRs (at point B in Figure 1 (a)) to train FNNs before the simulation of the singlecarrier system.
VB1 Performance Evaluation
System  Algorithm  MIMO scheme pool 

singlecarrier  [8, 9, 10, 11, 12, 13, 14, 15]  types 
[16]  types  
[7]  types  
Proposed algorithm  types  
STBCOFDM  [17, 18, 19]  types 
Proposed algorithm  types  
SFBCOFDM  [20, 21, 22, 23]  types 
Proposed algorithm  types 
First, we compare the size of the MIMO scheme pool between existing algorithms and the proposed algorithm, as shown in Table III. The data on the existing algorithms are based on the survey [1]. Table III shows that the featurebased algorithms[8, 9, 10, 11, 12, 13, 14, 15] are only able to identify less than 5 types of MIMO schemes since many MIMO schemes have the same redundancy locations. The algorithm in [16] can identify 11 types of MIMO schemes owing to the preprocessing operation which leads to a finer discriminatory capability. However, this capability depends on apriori information including the number of transmit antennas and channel coefficients. The algorithm in [7] utilizes the code rate to construct a likelihood function which quantifies MIMO schemes with different code rates. There are three approaches introduced in [7]. The first two approaches require apriori information including the number of transmit antennas, channel coefficients and noise power, while the third approach referred to as codesparameter (COP) based, only requires the number of transmit antennas. The proposed algorithm can identify 17 types of MIMO schemes, and thus has a wider applicability.
For a fair comparison, we compare the combination of AIC/MDL[2] and COP[7] with the proposed algorithm for the identification of both the number of transmit antennas and MIMO scheme using the same parameters of the singlecarrier system described in Table II. Figure 4 shows that the algorithms in [2, 7] have a better performance in the lowSNR regime. The reason is that the channel in these algorithms is assumed to be a flatfading channel which leads to an accurate eigenvalue analysis, while the proposed algorithm employs the heuristic method so that it is compatible with the OFDM system. In addition, the probabilities of correct identification of AIC and COP do not converge to one due to the inconsistency of AIC[2] and COP has a smaller pool size (only identifies 13 types).
Algorithm  Computational complexity  Matlab runtime (Average) 

AIC+COP  
MDL+COP  
Proposed  
Proposed (OFDM) 
From a practical point of view, it is important to analyze the computational complexity of the proposed algorithm, which is , where floating point operations are needed for the eigenvalue decomposition. It is worth noting that this complexity has the same order as those of AIC/MDL or COP, since they require similar operations including the covariance matrix construction and eigenvalue decomposition. We also verify the practical runtime of these algorithms using a computer with a Core i77700T CPU, 16 GB RAM (the simulation software is Matlab 2017b). The average run time of the proposed algorithm is about 1.55 , while the combination of AIC/MDL and COP takes about 2.97 or 3.12 . The complexity results are shown in Table IV.
VB2 Effect of the Number of Processed Samples
VB3 Effect of the Number of Receive Antennas
Figure 6 demonstrates that the probability of correct identification increases with the number of receive antennas for the assumed default simulation parameters. In fact, a large number of receive antennas enhances the estimation performance of the signal subspace dimension, since signals are mapped into a higher dimensional space where it is easy to discriminate between the features. It should be noted that for , the probability of correct identification decreases in the high SNR regime. In this case, the number of the input nodes of the NTAF FNN is 4, which is equal to the theoretical when . In the low SNR regime, the noise effect dominates the signal eigenvalues and the radii are . The NTAF FNN classifies these 4 circles as the Gerschgorin circles of the signal subspace and determines the result to be 4. In the high SNR regime, the radii are and is classified as the Gerschgorin circle of the noise subspace by the NTAF FNN.
VB4 Effect of Modulation Type
Figure 8 shows the effect of the modulation type on the identification performance. We evaluate the performance with four modulation schemes: 4PSK, 8PSK, 16QAM, 64QAM. These modulations are mandatory for most of the wireless standards. The results demonstrate that the performance does not depend on the modulation type.
VB5 Effect of Timing Offset
Results presented so far assumed perfect timing synchronization. Now, we evaluate the performance of the proposed algorithm under timing offset. The effect of the timing offset effect can be modeled as a twopath channel [43], where is the normalized timing offset when the whole sampling period is one. Figure 8 shows that the timing offset has a small effect on the performance of the identification of the number of transmit antennas, while the effect can be significant on the performance of the identification of MIMO scheme. The reason is that the timing offset destroys the orthogonality of the STBC, which introduces extra terms for the ICSF and IRSF and leads to the overestimation of the features.
VB6 Effect of Frequency Offset
Figure 10 illustrates the identification performance of the proposed algorithm for different frequency offsets. The frequency offset, , is normalized to the signal bandwidth. The identification of the MIMO scheme is impacted by the frequency offset, while the enumeration of the number of transmit antennas is robust with respect to this model mismatch. This is because the frequency offset rotates complex signals and destroys the orthogonality of STBCs, while it does not impact the independence of channels between transmit and receive antennas.
VC Simulation Results for STBCOFDM System
For robust OFDM signals identification, we retrain the FNNs using OFDM signals instead of singlecarrier signals. The number of subcarriers, OFDM symbols and detectors are set to 256, 500 and 64, respectively. Our proposed algorithm can identify a larger MIMO scheme pool than existing algorithms, as shown in Table III. In addition, the complexity of the proposed algorithm is due to the use of detectors, as shown in Table IV.
VC1 Effect of the Number of Detectors
Figure 10 demonstrates the viability of the proposed algorithm for STBCOFDM systems and presents the identification performance for different numbers of detectors, denoted by . The performance improves significantly as the number of detectors increases from 1 to 16 with diminishing performance gains beyond 16 detectors. This result indicates that the combination in (40) converges fast with increasing . It should also be mentioned that employing one detector makes the proposed algorithm degenerate into the singlecarrier system. It can be observed that the performance curve for in Figure 10 is between the curves of and in Figure 6 since one detector handles 400 samples for the OFDM systems.
VC2 Effect of the Number of OFDM Symbols
Figure 12 shows that increasing the number of OFDM symbols, , improves the performance of the proposed algorithm, which performs well even under a limited observation period.
VD Simulation Results for SFBCOFDM System
We use the same FNNs as in the STBCOFDM system to identify the number of transmit antennas and MIMO scheme in the SFBCOFDM system. The proposed algorithm can identify more MIMO schemes than previous works, as shown in Table III, and has a similar computational complexity to the STBCOFDM system.
Figure 12 verifies the viability of the proposed algorithm for the SFBCOFDM system and illustrates the identification performance for different numbers of detectors. It can be observed that the proposed algorithm for the SFBCOFDM system has a close performance to that for the STBCOFDM system. This is because the detector combinations (see Equations (40) and (43)) are the same for these two systems with the same parameters.
Vi Conclusions
We proposed a novel joint blind identification algorithm for the number of transmit antennas and MIMO schemes. After restructuring the received signals, three quantified features based on the dimension of the signal subspace, namely, NTAF, ICSF and IRSF, are derived to discriminate between different numbers of transmit antennas and MIMO schemes. Then, we proposed a neuralnetwork Gerschgorin radiibased method to calculate the three features and used a minimal weighted norm1 distance metric for decision making. Taking advantage of the quantified features and the neuralnetwork Gerschgorin radiibased method, the proposed algorithm can identify a large number of MIMO schemes and applies to both singlecarrier and OFDM systems. In addition, the proposed algorithm has an acceptable computational complexity and does not require apriori information on the channel coefficients, modulation type or noise power. The simulation results demonstrated the viability of the proposed algorithm for a short observation period both in the singlecarrier and OFDM systems, and showed an acceptable performance in the presence of small timing and frequency offsets.
Appendix A SFBC Examples
The code matrix of the SM and AL are, respectively, defined as