Enhanced interpolation-based interference rejection combining for black-space cognitive radio in time-varying channels

In this paper, we investigate multi-antenna interference rejection combing (IRC)-based black-space cognitive radio (BS-CR) operation in time-varying channels. The idea of BS-CR is to transmit secondary user (SU) signal in the same frequency band with the primary user (PU) such that SU’s power spectral density is clearly below that of the PU, and no significant interference is inflicted on the PU receivers. We explore the effects of interpolation and mobility on the novel blind IRC technique which allows such operation mode for effective reuse of the PU spectrum for relatively short-distance CR communication. We assume that both the PU system and the BS-CR use orthogonal frequency division multiplexing (OFDM) waveforms with common numerology. In this case, the PU interference on the BS-CR signal is strictly flat-fading at subcarrier level. Sample covariance matrix-based IRC adaptation is applied during silent gaps in CR operation. We evaluate the effect of the gap length on the link performance, using known PU channel-based interference covariance estimation as reference. We also propose an interpolation-based scheme for tracking the spatial covariance in time-varying channels, demonstrating significantly improved robustness compared to the earlier scheme. The performance of the proposed IRC scheme is tested considering terrestrial digital TV broadcasting (DVB-T) as the primary service. Also joint PU and co-channel CR interference cancellation is included in the study. The resulting interference suppression capability is evaluated with different PU interference power levels, silent gap durations, data block lengths, and CR device mobilities.

primary signal in the same time-frequency resources without causing objectionable interference, has received limited attention [5][6][7][8]. An underlay CR is ignorant about the existence of PUs in its frequency band. Commonly, it uses very low power spectral density and wide bandwidth, such that it does not cause interference to the PU transmission under any conditions [4]. BS-CR adapts its waveform and signal parameters depending on the ongoing PU transmissions and uses advanced signal processing techniques on the receiver side to facilitate low signal-to-interference ratio (SIR) at the receiver. BS-CR systems effectively reuse the spectrum for communication over short distances. It can operate with limited spectrum resources and can be used without any additional spectrum sensing.
As discussed in our previous papers [9,10], one of the major requirements for CR operation is to minimize the interference to the primary transmission system. In BS-CR, this is reached by setting the CR transmission power at a small-enough level. The most important factor that enables such a radio system is that stronger interference is easier to deal with as compared to weaker interference [11], if proper interference cancellation techniques are utilized. Previous studies from information theory provide theoretically achievable bounds for such cognitive radios [12]. Multi-antenna systems allow for spatiotemporal signal processing, which do not only improve the detection capability of the receiver but also improve performance in fading multipath channels with interference. Various methods of interference cancellation can be found in [13,[14][15][16][17][18] and the references therein. All other detection algorithms except the multi-user detector perform sub-optimally [13].
The interference rejection combining (IRC) receivers have the significant advantage in comparison with the other receivers in multi-user scenarios that they do not need detailed information about the interfering signals, such as modulation order and radio channel propagation characteristics. For CR scenarios, IRC receivers in general are simple and desirable compared to optimum detectors. IRC techniques are widely applied for mitigating co-channel interference, e.g., in cellular mobile radio systems like LTE-A [19]. The use of multiple antennas in CRs has been studied earlier, e.g., in [18]. Our initial study on this topic in highly simplified scenario with suboptimal algorithms was in [20], but to the best of our knowledge, IRC has not been applied to BS-CR (or underlay CR) elsewhere. In this current study, we develop the ideas that we presented in [9] under more practical situations and study the performance of our algorithms in more details. Notably, the scheme studied in [9] was found to be very sensitive to mobility, because its performance is critically affected by errors in spatial covariance estimation. In this work, we extent our previous studies on the effects of mobility and propose a scheme to improve the quality of the covariance estimation with time-varying channels using interpolation between sample covariance-based estimates.
In this paper, we consider BS-CR operation in the terrestrial TV frequency band, utilizing a channel with an ongoing relatively strong TV transmission. The PU is assumed to be active continuously. If the TV channel becomes inactive, this can be easily detected by each of the CR stations in the reception mode. Then the CR system may, for example, continue operation as a spectrum sensing-based CR system. In our case study, we focus on a scenario with multi-antenna CR receiver having IRC capability to mitigate co-channel interference generated by a single PU transmitter and multiple independent (uncoordinated) CR systems. Generally, interference covariance matrix estimation is a key element of IRC schemes. For this purpose, we propose to use sample covariance calculated from the received signal during silent periods (gaps) in the target CR transmission. More specifically, the contributions of this paper are listed below: • This paper extends and refines our earlier studies in papers [9,10] on IRC schemes operating under severe interference conditions, with applications in BS-CR. • Under such conditions, the direct estimation of the target CR channel is not possible.
After estimating the interference covariance matrix during silent gaps in the target CR transmission, our scheme extracts from the received signal a set of signals which span the interference-free subspace in spatial domain. The weight vector for the target CR signal is based on pilot-based target channel estimation in this subspace. • In BS-CR scenarios, the target CR link performance is very sensitive to the quality of the interference covariance estimate. Therefore, we propose to use linear interpolation of the interference covariance matrix to track the channel variations between consecutive interference covariance estimates, i.e., consecutive silent gaps. • Our case study focuses on a TV whitespace (TVWS) scenario. We consider a multiantenna CR receiver having IRC capability to mitigate co-channel interference generated by a single PU transmitter and multiple independent (uncoordinated) CR systems. The performance of such a system under different interference levels, mobilities, silent gap lengths, frame structures, and modulation orders is studied, and the improved robustness obtained through covariance interpolation is highlighted The rest of the paper is organized as follows: In Sect. 2, the BS-SC scenario and proposed IRC scheme are first explained. Then also the connections between the signal-tointerference ratios (SIRs) at the PU and CR receivers are analyzed based on the channel losses of the target and interference links. The system model and IRC process are formulated in Sect. 3. Section 4 presents the simulation setup and performance evaluation results. Finally, concluding remarks are presented in Sect. 5.

IRC-based black-space cognitive radio scenario-scenario and methods
In our basic scenario, illustrated in Fig. 1, we consider a CR receiver using multiple antennas to receive data from a single-antenna cognitive transmitter. The details of the receiver signal processing indicated in the figure are explained in Sect. 3. The CR operates within the frequency band of the PU, and the PU power spectral density (PSD) is very high in comparison with that of the CR. The primary transmission is assumed to be always present when the CR system is operating. The primary transmitter generates a lot of interference to the CR transmission, which operates closer to the noise floor of the primary receiver, and due to this, the primary communication link is protected. We consider frequency reuse over relatively small distances, such as an indoor CR system. The multi-antenna configuration studied here is that of single-input multiple output (SIMO). Other configurations, involving also transmit diversity in the CR link, are also possible, but they are left as a topic for future studies.
Here the PU is a cyclic prefix orthogonal frequency division multiplexing (CP-OFDM)-based DVB-T system [21]. The CR system is also an OFDM-based multi-carrier system using the same subcarrier spacing and CP length as the primary system. Thus, it has the same overall symbol duration. The CR receiver is assumed to be synchronized to the primary system in frequency and in quasi-synchronous manner also in time. The CP length is assumed to be sufficient to absorb the channel delay spread together with the residual timing offsets between the two systems observed at the CR receiver. Consequently, the subcarrier-level flat-fading circular convolution model for spatiotemporal channel effects applies to the target CR signal and to the PU interference signal as well. Then the IRC process can be applied individually for each subcarrier. Since the CR receiver observes the PU signal at very high SINR level, synchronization task is not particularly difficult and low-complexity algorithms can be utilized. Considering shortrange CR scenarios, the delay spread of the CR channel has a minor effect on the overall channel delay spread to be handled in the time alignment of the two systems. The same applies to the possible significant co-channel CR interferers, which are assumed to be at relatively short distance from the target CR system.
Basically, if all target and interfering CR stations are synchronized to their respective observed PU signals, they are also synchronized with each other in the quasi-synchronous manner, such that all multi-path components of all signals are within the CP length.
In addition to these, we assume that the secondary user experiences slow (e.g., pedestrian) mobility effects, due to either the mobility of devices or people moving around them. The effects of mobility on the CR system performance are evaluated in this study.
Both the primary and the CR systems use QAM subcarrier modulation, but usually with different modulation orders. The frame structures and pilot patterns of the PU and CR are independent. The received CR signal consists of contributions from the desired CR communication signal, co-channel CR interferers, and the primary transmission signal, the latter one constituting the dominating interference. Our proposed scheme includes two phases in the CR system operation as described in our previous work [9]. The spatial characteristics of the PU interference are modeled using multiantenna sample covariance matrix, which is estimated during silent gaps in the target CR transmission, independently for each active subcarrier. No explicit channel estimation of the PU channel is required. The CR channel is estimated from the partial IRC signals, from which the PU and other CR interferences have been effectively suppressed.

Methods
In this work, we study how channel fading with mobility affects the performance in the scenario described above. We consider adaptation of the IRC process using interpolation of the sample covariance matrices between the silent gaps. This is expected to improve the performance, allowing to increase the data block length between silent gaps, thus reducing the related overhead in throughput. The main elements of the scheme, covariance estimation and interpolation, IRC process, channel estimation [22], and MRCbased detection are formulated mathematically. Numerical performance results are obtained through Monte Carlo simulations based on a developed Matlab script, using a semi-analytical model based on PU channel knowledge as reference.

SIR analysis
Key parameters in BS-CR operation are the SIRs at the CR receiver and PU receiver. In the basic scenario with PU interference only, the SIRs can be expressed in dB units as follows: where P CR and P PU are the transmission powers of CR and PU, respectively, in dBm units, while the other parameters are transmission losses in dB units of relevant transmission links. The first part of the subscript indicates the transmitter, and second part the receiver of the corresponding link. From these equations, it is easy to derive the maximum SIR of the CR receiver in terms of the minimum SIR of a PU receiver, SIR PU,min : The maximum SIR of the CR link depends on the differences of the channel losses from CR TX to both receivers and from PU TX to both receivers. For example, if the losses of the main links (CR-CR and PU-PU) are equal to the losses of the corresponding interfering links (CR-PU and PU-CR), then the CR RX would be able to operate with the SIR of −SIR PU,min . If the main link losses are lower than those of the interference links, the SIR CR,max would be higher. The critical cases are: 1 If the losses of the main links (CR-CR and PU-PU) are much higher than the losses of the corresponding interfering links (CR-PU and PU-CR), then SIR of the CR receiver may be limited by the minimum acceptable PU SIR to a highly negative value. Then the SIR of the CR receiver may be limited to a highly negative value by the maximum transmit power of the CR device, which could be much lower than that of the PU transmitter. 2 If the PU transmitter is very close to the (short-range) CR system, then Then the SIR of the CR receiver may be limited to a highly negative value by the maximum transmit power of the CR device, which could be much lower than that of the PU transmitter. In both cases, the performance of the IRC scheme would be limited by the RF imperfections (e.g., receiver nonlinearity effects) in practice. The latter limitation appears especially in the TVWS application, but in rather limited geographical regions. Considering the first kind of limitation in the worst-case scenario, one should consider the hidden node margin (typically around 25-30 dB in the TVWS case [23] for the loss difference from PU (L PU−PU − L PU−CR ) , which would lead to very low SIR for the CR TX, unless the loss from the CR to the PU RX can be guaranteed to be much lower than that of the CR link. On the other hand, in short-range CR communication, especially if PU receivers use directive (possibly roof-top) antennas, the mentioned case where the SIR of CR RX is no less than −SIR PU,min should be possible.
Generally, the CR TX should have knowledge of the link losses in order to maximize its transmission power while not causing excessive interference to the PU link. This paper focuses on the physical layer capabilities of IRC-based BS-CR, and detailed procedures for controlling the CR operation remain as a topic for future studies. One possible way to enhance BS-CR operation would be cooperative spectrum sensing of CR devices to reliably estimate the PU power level in the CR operation region.

IRC for black-space cognitive radio
Based on the OFDM model mentioned above, subcarrier-wise detection is considered with flat-fading channel coefficients to get rid of the challenge of frequency selectivity in the IRC process.
In the SIMO configuration, the CR is assumed to have N receiver antennas and L different interference sources are assumed. Based on this model, the signal received by the CR can be formulated for each active subcarrier as follows: Here x T is a transmitted subcarrier symbol and h T = [h T,1 , h T,2 , . . . , h T,N ] T is the target channel vector with N receiver antennas in the CR, x I,l is the lth interfering signal, and h I,l = h I,l,1 , h I,l,2 , . . . , h I,l,N T is the channel vector for the lth interferer. The channel vectors consist of the complex channel gains from the corresponding transmit antenna to nth antenna of the CR receiver. Finally, η is the additive white Gaussian noise (AWGN) vector. In this generic system model, it is assumed that the PU is the dominant interferer, and the other interference sources are other CR systems introducing co-channel interference at relatively low power level.

Covariance matrix estimation
As it is illustrated in Fig. 1, the interference minimizing IRC weights are obtained during a silent period in the target CR transmission, while the co-channel interference sources are active. Due to that, Eq. (4) can be modified during silent gaps of CR operation as Here it is assumed that only interferences and noise are present during the silent period in the signal r observed by the CR. Linear combiner is used for the signals from different antennas with a weight process in detection as follows: where y is the detected signal, w is the weight vector with N elements, and superscript H denotes the Hermitian (complex-conjugate transpose).
Determining the optimum weight values is an optimization problem which can be solved with the linear minimum mean-squared error (LMMSE) criterion [21,22] that aims to minimize the mean-squared error with respect to the target signal x T , When knowledge of the covariance matrix is available, IRC can be applied. Two cases are considered below: (i) calculating interference covariance from known channel coefficients and (ii) sample covariance-based approaches.

Perfect channel information case
The deterministic known channel case (with perfect channel information) is considered as the theoretical bound in this study. Assuming that the channel vectors form the interferers are perfectly known, the noise plus interference covariance matrix can be calculated for each subcarrier as where P l is the variance of interferer l at the transmitter, P N is the noise variance that can be obtained from the SNR, and I is the identity matrix of size NxN. Assuming that also the channel vector for the target signal is known and the target and interfering signals are Gaussian and statistically independent from each other, the conventional LMMSE solution for the weight vector is [15] where P T is the target CR transmitters signal power and unit noise variance is assumed.
In the BS-CR scenario with a single dominant interferer, the estimation of the PU channel is relatively straightforward if the CR knows PUs pilot structure. However, in the case of multiple interferers, the channel vectors of all interferers should be estimated, which becomes quite challenging. Furthermore, the target channel cannot be estimated before the interference cancellation. Therefore, the perfect channel information case serves mainly as an ideal reference in performance comparisons, as will be seen in Sect. 4.

Sample covariance-based case without PU channel information
The IRC process starts from interference covariance matrix estimation during silent gaps in the receiver. It is difficult to have the perfect channel state information on the CR receiver side. Alternatively, the joint interference and noise covariance matrix can be estimated for each subcarrier by the sample covariance matrix of the corresponding received subcarrier signal (after CP removal and FFT) in the absence of the target transmission, i.e., during the silent gaps as.
Here m is the OFDM symbol index and M is the observation length in subcarrier samples, which is chosen equal to the length of the silent gap.

Linear interpolation for interference covariance tracking with mobility
Regarding the mobility aspects, there are significant differences in the effects of PU transmitter mobility, CR transmitter mobility, and CR receiver mobility. If the PU transmitter is stationery and CR receiver is stationary, the mobility of CR transmitter is easier to handle, because the dominating PU interference is stationary, and the variations in the noise and interference covariance matrix are only due to the co-channel CR interferes. However, even in this case, radio environment of the CR receiver may vary due to movement of people or vehicles nearby. Therefore, some tolerance to mobility is required also in such scenarios, at least with pedestrian mobilities. The mobility of PU transmitter or CR receiver makes the dominant interference time-varying, and in the BS-CR scenario, the CR link performance is very sensitive to quality of the PU interference covariance matrix estimate. Therefore, it is important to investigate these mobility effects and consider enhanced schemes for tracking the interference covariance with mobility.
While considering the sample covariance-based approach, increased observation (silent gap) length gives better PU interference covariance estimate in the stationary case or with low mobility. However, the channel variations during the silent gap affect critically the quality of the PU interference covariance matrix. Therefore, the optimum length of the silent gap (i.e., observation length) depends on the mobility, i.e., Doppler spread. Then we apply linear interpolation for the covariance matrix elements when calculating the weight vectors for the data symbols between two consecutive silent gaps.
There are two key parameters in this process, the silent gap length and the data block length between two consecutive gaps. Increasing the gap length improves the performance with low mobility but degrades the performance with higher mobility and (10) increases the overhead in throughput. Increasing the data block length increases the throughput but degrades the performance with mobility. These tradeoffs are investigated through simulations in Sect. 4 of this paper.

IRC process
As indicated in the model shown in Fig. 1, the CR channel cannot be estimated before the step of the interference cancellation, and the optimal steering vector for IRC cannot be directly calculated. Here we develop an IRC scheme which utilizes N orthogonal virtual steering vectors in the receiver's internal channel estimation process, which is based on pilot symbol structures typically used in OFDM systems. Actually, N-L vectors would be enough since IRC consumes L degrees of freedom. But this way the model is more straightforward and the IRC process is generic and robust since there is no need to estimate the number of interferers. Without loss of generality and for simplified computations, the following unit vectors are applied as the virtual steering vectors, Basically, the weight vectors obtained when applying the virtual steering vectors for w in Eq. (9) span the interference-free subspace in spatial domain. Assuming that the interference covariance is correctly estimated, using any linear combination h of these vectors instead of h T in (9) provides interference cancellation and maximizes signal power for the spatial CR channel h. The following process finds the combination of the virtual steering vectors for the actual target signal by first estimating the spatial channel vectors corresponding to the different virtual steering vectors. The obtained weight vectors are as follows: It should be noted that the denominator of Eq. (9) is a complex scaling coefficient, which will be included in the MRC weights.

Target channel estimation with linear interpolation and MRC combining
In the second stage, data symbols of the CR link are transmitted together with the training/pilot symbols. While receiving pilot symbols, the weighted output signals corresponding to each virtual steering vector are calculated as The IRC process cancels the interference from all of the weighted output signals, y n corresponding to different virtual steering vectors. For each subcarrier, the N channel coefficients for each of the weighted output signals can be estimated using the pilot symbols as follows:   where p is the transmitted pilot symbol value.

Linear interpolation for channel estimation
In the traditional pilot-based channel estimation process, it is required to use efficient interpolation techniques, such as Wiener interpolation, based on the channel information at pilot sub-carrier symbols. For simplicity and to avoid excessive received signal buffering over high number of OFDM symbols, we apply linear interpolation. The performance of linear interpolation technique is better than the piecewise-constant interpolation methods [24,25]. In the simulation studies of Sect. 4, a basic training symbol scheme is assumed where training symbols contain pilots in all active subcarriers, and the interpolation is done in time domain only, between two consecutive pilots in each subcarrier. The MRC weights for a data symbol are then calculated as where g V,n , n = 1, . . . , N , denote the corresponding interpolated channel estimates. Generally, the interpolated channel estimates are different for each subcarrier in each OFDM symbol. For simplicity, the subcarrier and OFDM symbol indexes are not included in the notation above. Then the interpolated channel estimates are given as where m is the OFDM symbol index, k is the subcarrier index, S is the pilot spacing, S p = m S · S S is the preceding training symbol index,S f = m S · S is the following training symbol index, and ⌊. ⌋ and ⌈. ⌉ stand for the floor and ceiling operations, respectively.

Combining for detection
In the final stage, while receiving data symbols, the equalized data symbols are calculated by maximum ratio combining the N samples obtained by applying the virtual steering vectors. The effective weight vectors for CR can be obtained as, where g = g V,1 , g V,2 , . . . , g V,N , T . It is enough to calculate and use this weight vector w CR , instead of separately applying the MRC weights on the samples obtained by the weight vectors [w V,n ] . The equalized data symbols are then calculated as follows: (15) g V,n = y n p = w H V,n · r p , n = 1, 2, . . . , N , In Eq. (19), g appears as the spatial channel estimate for the target CR transmitter. Overall, this is the zero-forcing IRC solution which maximizes the received signal power over normalized steering vectors.

Results and discussion
The simulations are carried out for the system setup explained in Sect. The number of spatial channel realizations simulated in these experiments is 300-1000. The ratio of CR and PU signal power levels at the CR receiver (referred to as the signal-to-interference ratio, SIR) is varied. In the case of co-channel CR interference, the average power levels of interfering and target CRs are the same at the target receiver and the channels are independent instances of the Vehicular A model with random timing offsets, while all multi-path delays remain with the CP. The lengths of the OFDM symbol frame and silent gap for interference covariance matrix estimation are also varied (expressed in terms of CP-OFDM symbol durations). A very basic training symbol scheme is assumed for the CR: training symbols contain pilots in all active subcarriers and the spacing of training symbols is 8 OFDM symbols. Frame length is selected in such a way that training symbols appear as the first and last symbol of each frame, along with other positions. Channel estimation uses linear interpolation between the training symbols. We have tested the BS-CR link performance with SIR values of {− 10, − 20, − 30} dB using silent gap durations of {8, 16, 32} OFDM symbols (known channel case also with 128 symbols), and data block lengths of {17, 25, 33, 41} OFDM symbols.

The effect of silent gap length
This subsection analyzes the BER performance of the proposed sample covariance-based IRC process using Eqs. (10), (12), (18), and (19), considering the known channel-based process as an ideal reference. As it was explained in Sect. 3, the known channel case assumes perfect knowledge of the interference channel and it provides a theoretical (19) performance bound for practical IRC schemes. It uses the LMMSE solution (8) with known channel-based covariance estimate (9). Otherwise, the receiver process is the same as in the sample covariance-based scheme, thus providing a theoretically achievable bound for the practical sample covariance-based approach. Below, the known channel model is applied with PU interference only, assuming no cochannel CR interferers. Figure 2 shows the BER performance considering both the known channel and sample covariance-based approaches in stationary case (no mobility). Here, the silent gap durations are 8, 16, 32, and 128 OFDM symbols and the data block length is 17. With zero-mobility, link performance is independent of the data block length and short block is used here mainly to reduce simulation time. It can be observed that the sample covariance-based simulation results converge to the corresponding known channel results with increasing gap length, since the covariance estimate is improved with increasing sequence length. This demonstrates that the known channel-based bound is theoretically achievable.
In Fig. 3, the BER performance of known channel and sample covariance approaches is shown considering SIR values of {− 10, − 20, − 30} dB. Similar to Fig. 2, the data block length is 17 and the gap duration is selected as 16 OFDM symbols. Zero-mobility case is considered also here. As seen in Fig. 3, the CR link performance of both known channel and sample covariance schemes is rather independent of the SIR.
A detailed comparison between the required SNR values of the known channel and sample covariance-based approaches for BER = 0.01 is provided in Table 1

Performance analysis of proposed scheme with mobility
This subsection reports the performance analysis of the proposed algorithms considering also the effects of mobility and different system parameters. It is noted that only sample covariance-based approach is considered in the following results. The results presented above still serve as ideal reference when evaluating the effects of mobility in different configurations. First we evaluate the performance with slow mobilities with and without pilot interpolation, still assuming that the interference is due to the PU only. Figure 4 shows the impact of covariance matrix interpolation on the BS-CR link performance. Here the data block length and gap duration are fixed to 17 and 16 OFDM symbols, respectively. This choice provides performance that is no more than 1 dB from the configuration reaching 1% or 10% BER with lowest SNR, among the tested configurations with even higher overhead. We can see that covariance interpolation provides significant improvement of robustness in time-varying channels. Focusing on the 1 -10% BER region, the performance with interpolation at 10 km/h mobility clearly exceeds the performance at 3 km/h without interpolation. However, for the 64QM case with -30 dB SIR, this is true only for BER of 10% or higher, due to the high error floor at very low  SIR and high mobility. We can also see that with stationary channel, the performance is practically independent of the SIR. (Compared with [9], the performance is slightly improved by fine-tuning the used algorithms.) In Fig. 5, the effect of silent gap duration is tested with 16 QAM modulation, SIR of − 20 dB, and data block length of 17. The overhead in data rate is about 58% and 44% for gap lengths of 16 and 8, respectively. The shorter gap length results in about 1.5 dB performance loss in the 1-10% range in stationary case and about 1.8-3.5 dB loss with 10 km/h mobility, compared to the gap length of 16. With 20 km/h mobility, the corresponding loss is about 2.2 dB at 10% BER, but longer gap leads to higher error floor, and the performance with shorter gap becomes better for BER below 3%.
Finally, we evaluate the link performance in the presence of co-channel CR interference, in addition to PU interference. While still assuming four antennas in the target CR receiver and − 20 dB SIR for the PU signal, also two interfering CR signals are included a b Fig. 4 Performance of QPSK and 64QAM BS-CR systems for SIR {− 10, − 30} dB with or without covariance matrix interpolation. a 0 and 3 km/h CR receiver mobilities. b 0 and 10 km/h CR receiver mobilities. Silent gap length of 16 symbols, and OFDM frame length of 17 symbols in the model. All CR signals are at the same average power level and their channels are independent instances of the Vehicular A channel model. The results are shown in Fig. 6, indicating that the impact of co-channel CRs on the target link performance is very minor.

Conclusion
The performance of black-space CR transmission links in the presence of strong interferences and mobility was investigated using spatial covariance interpolation between silent gaps. The interference rejection capability of IRC using multiple receive antennas  for various modulation orders under varying mobility and channel setups was studied. It was found that the IRC performs very well in the basic SIMO-type BS-CR scenario when stationary channel model is applicable, e.g., in fixed wireless broadband scenarios. However, the scheme is rather sensitive to the fading of the PU channel, e.g., due to people moving close to the CR receiver. Due to the strong interference level, the interference cancellation process is affected by relatively small errors in the covariance matrix estimate. For covariance estimation, the silent gap length in the order of 16-32 OFDM symbols provides the best performance with stationary channels, but even with 3 km/h mobility, the performance degrades greatly when considering SIR levels below -10 dB. The data block length should be of the same order or less, which leads to high overhead due to the silent gaps. Covariance interpolation was shown to greatly improve the robustness with time-varying channels, such that good link performance can be obtained with up to 20 km/h mobility at 700 MHz carrier frequency. This indicates that the proposed BS-CR scheme could be feasible at below 6 GHz frequencies with pedestrian mobilities. However, there is a significant tradeoff between link performance and overhead in data rate due to the silent gaps.
In the basic TV black-space scenario, there is only one strong TV signal present in the channel, in agreement with our assumption about the primary interference sources. DVB-T system allows also single-frequency network (SFN) operation and the use of repeaters to improve local coverage. In both cases, the primary transmissions can be seen as a single transmission, with a spatial channel that depends on the specific transmission scenario, and the proposed scheme is still applicable.
One important issue in the proposed scheme is its sensitivity to the nonlinearities of the CR receiver's analog front-end. Wide linear range is required in order to prevent nonlinear distortion from the high-power PU signal from degrading the CR link performance. This is a common issue with opportunistic CR operating in white spaces close to high-power PU channels and also with digital signal processing (DSP) intensive receiver architectures. An interesting technology in this context is advanced DSP algorithms for compensating the nonlinear effects of the receiver's analog front end [26]. On the other hand, sample covariance-based IRC may exhibit some capability to reject also the nonlinear distortion due to the strong PU signal. This remains as an important topic for future studies.
The scheme was also extended to scenarios where multiple CR systems are operating in the same region. If all CR systems are time-synchronized to the PU and they are at a relatively small distance from each other, they are also synchronized with each other, and can be handled by the IRC process as additional interference sources following the model of Eq. (3).
In future work, it is important to optimize the silent gap and data block lengths along with the modulation order to maximize throughput with given PU interference level and mobility. Lower-order modulations are more robust to errors in covariance estimation, allowing significantly lower gap and training overhead than higher order modulation. Complexity reduction of the covariance interpolation and IRC process is also an important topic for further studies. It is also worth to consider adaptation of the IRC process without silent gaps after the first one required for the initial solution. This would help to reduce the related overhead in throughput. One possible approach is to do this in a