Abstract

We propose a modified transform domain communication system (TDCS) structure with weighted fractional Fourier transform (WFRFT) preprocessing over the channel with additive white Gaussian noise (AWGN) and narrowband interference (NBI). The majority of NBI can be removed by adding zeros to TDCS’s basic function (BF) amplitude vector at these subcarriers where strong NBI exists, and the corresponding residual NBI can be further greatly eliminated in time-frequency domain through WFRFT processing. To achieve the best bit error rate (BER) performance, we define as the normalised cross relationship coefficient between TDCS’s BF vector and residual NBI vector and prove that the corresponding WFRFT can be obtained by minimising the magnitude of . Numerical simulations are conducted to evaluate the proposed system’s performance with regard to NBI suppression.

1. Introduction

In current carrier schemes, multicarrier (MC) and single carrier (SC) communication systems have become the two major techniques. To make best use of these two techniques’ property, the long term evolution (LTE) system adopts orthogonal frequency division multiplexing (OFDM) for downlink transmission and single carrier frequency division multiple access (SC-FDMA) for uplink transmission [1]. Considering that OFDM and SC techniques have their individual advantages and disadvantages, researchers have done many attempts to converge these two carrier schemes [25]. Among them, the hybrid carrier (HC) system based on weighted-type fractional Fourier transform (WFRFT) is an inspiring attempt which can make a tradeoff between SC and MC. In addition, as a useful time-frequency analysing tool, WFRFT can make the original signal distributed more evenly in the time-frequency domain. The time-frequency distribution characteristics can greatly help mitigate intercarrier interference/intersymbol interference (ICI/ISI) [6] in additive white Gaussian noise (AWGN) or doubly selective channels [7] as long as the WFRFT order is set properly. Recently, our team has been extensively engaging in the related research on WFRFT. We explored WFRFT signals’ inherent characteristics, such as modulation recognition [8] and antiparameter scanning ability [9, 10]. Besides, we also demonstrated the feasibility and effectiveness of WFRFT’s applications in physical layer (PHY) secure communications [11] and covert satellite communications [1214]. In view of the above-mentioned facts, the WFRFT-based communication system will be a promising technology in the future communication system.

In the current broadband block-based transmission systems, most of them adopt the OFDM-based structure, then the WFRFT-based communication system also adopts the OFDM-based structure to be compatible with current communication structures. However, the OFDM-based communication systems are subject to the interference caused by narrowband interference (NBI), which has long attracted many researchers’ attention [15, 16]. As an important aspect of the anti-interference techniques, NBI suppression performance plays an important role in the system’s performance evaluation. Thus, it is necessary to carry out research on NBI suppression in WFRFT-based systems.

In [17], it has been pointed out that WFRFT precoding structure can make the residual NBI distributed more evenly in the time-frequency domain and then eliminate the residual NBI to most extent. Furthermore, in [18], it has also been testified that WFRFT processing can help suppress NBI in the traditional direct sequence spread system (DSSS), while how to obtain the optimal BER performance on has not been analysed.

To the best of our knowledge, according to the different ways to suppress the interference, these suppression methods can be classified into three categories: subtraction [1925], avoidance [2630], and transform process [3134]. In order to have a more direct understanding of the differences between all kinds of NBI suppression methods in OFDM-based systems, we list a table to describe the comparisons between different NBI suppression methods, which are shown in Table 1. In Table 1, as for the suppression effects, to show how the specific methods can help eliminate the NBI, we define three levels, namely, "A", "B", and "C". "A", "B", and "C" denote the suppression effects "almost all", "a large amount", and "moderate", respectively.

The basic principle of subtraction methods is to detect the NBI’s characteristic parameters and reconstruct the NBI based on these parameters. In particular, as for the compressive sense (CS) based subtraction methods, theoretically, NBI can be eliminated totally. However, as shown in Table 1, due to the reconstruction process, the prior information and iterative processes are often needed, which inevitably bring a great amount of computational burden. As for the avoidance methods, usually, by designing high performance filters to filter the NBI in time or frequency domain, the time/frequency domain filtering method can enhance the anti-interference property. In the meantime, it will also incur the performance reduction due to the useful signal’s energy loss in the filter’s stop bands and the high performance filters are also difficult to realise. Last, the transform process often does well in anti-NBI performance because of its capability to spread the interference into the corresponding transform domain. However, the universality of the transform process will also cause the insufficiency of NBI suppression.

Based on the above analyses, the traditional spread spectrum (SS) system has been proved an effective way to reduce the NBI or other interferences’ effects on the system’s bit error rate (BER) performance and has been widely used in the anti-interference communication system [35, 36]. Moreover, by utilising the basic function (BF) to spread the original transmitting signal in a wider bandwidth [37, 38], transform domain communication system (TDCS) can serve as a special spread spectrum system in some sense [39, 40]. Although TDCS can avoid the majority of the NBI by zeroing the amplitude on these carriers position where strong interference exists, we cannot guarantee the NBI have been fully mitigated, given that the residual NBI is still distributed in the time-frequency domain unevenly.

Based on the theory mentioned above, in this paper, the main idea of the proposed structure is to combine the two methods, namely, avoidance and transform process. Compared with the subtraction method, the proposed structure can suppress the NBI to a great extent with a proper WFRFT order , while it will not cause too much implementation complexity. A modified TDCS structure with the WFRFT preprocessing technique is provided to make full use of the techniques of TDCS and WFRFT. The proposed system is motivated twofold. First, by adopting the traditional TDCS technique, we can remove the most part of NBI. Second, through WFRFT, we can make the residual NBI dispersed in a wider time-frequency domain, which helps mitigate the residual NBI to the greatest extent.

This paper will be organised as follows: In Section 2, the overview of the modified TDCS structure with WFRFT preprocessing is described in detail. In Section 3, through theoretical deductions, the system’s BER performance has been derived, and the related principle has been established about how to select the optimal WFRFT order . In Section 4, numerical simulations are conducted to evaluate the performance of the proposed system. Conclusions are drawn in the last section.

Notation. , , , , and denote the real, imaginary, complex conjugate, transpose, and magnitude operations respectively. and are used to calculate the gradient and Hessian matrix of . , , and represent the identity matrix, normalised discrete Fourier transform (DFT), and inverse discrete Fourier transform (IDFT), respectively. means obeys the uniform distribution in , and means obeys the normal distribution with mean and variance .

2. System’s Description

Figure 1 illustrates the architecture of the proposed system combining TDCS spreading with WFRFT preprocessing. For description’s clarity, we will use Tx and Rx to denote the transmitting and receiving parts, respectively.

As for Tx, first, in the frequency domain expression, the basic function (BF) vector () in TDCS is composed of the spectrum sense ID () and pseudo phase vector (). Multiplied by the basic function vector , the original baseband signal can be transformed into , and . As a spread spectrum technique, TDCS makes every element in spread by the same BF vector , then we can omit the subscript for clarity’s sake and pay close attention to the spread vector , where is a variable generated from .

Then through the frequency mapping and WFRFT preprocessing, the original signal can be transmitted through the OFDM-based structure into the channel. As for Rx, we assume the corresponding spectrum sense ID and the pseudo phase vector are the same with the Tx; besides, considering that the system has the corresponding inverse processing modules, the related descriptions about Rx are not presented in detail. The main work of this paper is to construct the related frequency mapping matrix and analyse how different WFRFT order can influence the designed system’s capacity on NBI suppression. The channel is assumed to be only influenced by AWGN and NBI, so other adding cyclic prefix (CP) and frequency domain equalisation (FDE) modules are omitted. Detailed operating mechanism analyses on the designed system will be put as the following three subsections.

2.1. Spread Spectrum with TDCS

As for the BF vector in Figure 1, it consists of the spectrum sense ID and complex pseudo phase vector . By adding zeros to the amplitude vector at these carriers where strong NBI exists, we can remove most of the NBI in its frequency domain. Besides, through a phase mapping from a pseudo-number (PN) sequence, such as -sequence [41], can be regarded as random variables obeying uniform distribution in , . Then the corresponding time-domain BF vector can be given bywhere denotes the step-wise product.

According to the central limit theorem [42, 43], in a sense, in (1) can be regarded as a Gaussian-like sequence and then can serve as a spreading sequence, just like the PN sequence in the spread spectrum (SS) system [38]. However, unlike the traditional SS, the majority of interference will first be suppressed and further be spread in the frequency domain.

2.2. Time-Frequency Domain Dispersion with WFRFT

In the matrix form, the normalised 4-WFRFT can be described aswhere is the M-point normalised DFT matrix, whose -th element is ; are weighting coefficients, and , .

According to [17], the matrix can be implemented through the shift matrix , whose expression is given byThen, the WFRFT module can be implemented as Figure 2. As depicted in Figure 2, the S/P and P/S modules denote the serial-to-parallel and parallel-to-serial operations. And the DFT module and inverse module correspond matrix and , respectively. Furthermore, in Figure 1, through WFRFT, can be transformed into where is the frequency mapping matrix in Figure 1, details about will be shown in next subsection. In (4), are the frequency domain parts and can be regarded as the multicarrier components; are the time-domain parts and can be regarded as the single carrier components. That is the reason why WFRFT can make the original signal distributed in the time-frequency domain. Besides, different in (4) will produce different distribution characteristics.

2.3. WFRFT Preprocessing Matrix

As illustrated in Figure 1, to combine the traditional TDCS with WFRFT, we add the WFRFT preprocessing module between the modules IDFT and BF generator. In addition, to implement the elimination of the major NBI in TDCS, we utilise the frequency mapping matrix and its transpose matrix to zero those carriers where strong NBI exists, and is expressed byIn (5), when the amplitude of NBI at the carrier frequency is larger than the preset threshold , the corresponding amplitude will be set to zero, and ; when the amplitude of NBI at the carrier frequency is less than the preset threshold , the corresponding amplitude will be set to one, and , which means the -th element in will be mapped into the -th element in .

As for Tx, the final transmitting signal can be derived asTo analyse the combinatorial effects caused by the related matrix and , for the sake of clarity, we defined the WFRFT preprocessing matrix as follows:After substituting (7) into (6), we can rewrite (6) asAs for Rx, due to the corresponding inverse prepressing modules for TDCS and WFRFT, the original transmitting signal can be recovered to the most extent. However, when it comes to NBI, the preprocessing matrix can eliminate the major part of NBI and further make the corresponding residual NBI dispersed in the time-frequency domain. Detailed analyses on the system’s performance will be conducted in the next section.

3. Theoretical Performance

3.1. Anti-Gaussian Noise

As for Gaussian noise () and considering that the normalised 4-WFRFT is a unitary linear transformation, then have the same distribution properties with [17]. Consider that the demapping matrix will first zero some carrier frequencies and then make the residual NBI distributed in time-frequency domain. According to [44], through demapping with , the corresponding spread gain is , and then the final signal-to-noise rate (SNR) can be expressed bywhere is the baseband PSK modulation order and is the original SNR without the spreading process with TDCS BF sequence.

3.2. Anti-NBI

In complex baseband form, the expression of NBI can be described as [45]where is the period of each data symbol, is bandwidth ratio of wideband system to NBI, , means to round the elements of to the nearest integers greater than or equal to , denotes the central frequency of NBI, and denotes the -th data symbol.

As for Rx, through demapping matrix , for a symbol period ( points) in the designed system, the residual NBI signal can be given byThen we can define NBI to signal ratio (ISR) within per bit with TDCS spreading as follows:Then, to analyse the effects of the matrix and TDCS spreader to the residual NBI, we also define the cross-correlation coefficient between TDCS spreading sequence and the residual NBI as follows:where and correspond the -th element of and , respectively. According to [46], when the residual NBI cannot be treated as Gaussian random variables, as for BPSK modulated system over channels with AWGN and other interference, the system’s BER can be approximately be expressed aswhere , , and denote original signal energy and the noise power spectral density and the equivalent interference energy for one bit duration ( points), respectively, and is the complementary error function, whose expression is .

In the complex baseband forms, without loss of generality, we will take Gray-coded QPSK as the baseband modulation to deduce the related BER expression. As the real part and imaginary part of the Gray-coded QPSK baseband signal have the same BER performance with equal probability, the final BER conditioned on the residual NBI can be given bywhere , , and are corresponding to (9), (12), and (13), respectively.

3.3. Optimal Order Selection

As for the optimisation program in (15), when NBI, ISR, and SNR are determined, BER will just be the function of . Different WFRFT order will produce different BER performances. Through strict mathematic deductions, we will prove that the optimal WFRFT order can be obtained by minimising the magnitude of in (15). The related deduction is put as follows.

Let , and . Using derivative method, , then we can get the following expressions:Usually, ISR is a positive number, then the equation is equivalent to the one . When , that is to say, , according to (16), we can easily find that , . Then according to [47], we can achieve the local minimum of in (15) when ().

Meanwhile, we can calculate the Hessian matrix for in (15) as follows:After substituting , into (17), we can get the following expression:Obviously, the Hessian matrix in (18) is positive definite according to the convex optimisation theory [48]; the local minimum at () is the global minimum in (15). Meanwhile, means that the residual NBI signal and the BF vector are totally unrelated to each other, which is almost impossible to make it in the actual communications. However, based on the analysis of (15)-(18), we can still draw a conclusion that there exists an obvious positive correlation between the designed system’s BER and the cross-correlation coefficient magnitude . So, to get the best resistance to the residual NBI in the designed system, we can keep as small as possible by selecting the WFRFT order properly.

Based on the theoretical deduction mentioned above, the optimal program of (15) can simply come down to the following optimisation:As shown in (19), to be compatible with the current SC and MC system, we only search the optimal in the range even though the classical WFRFT order have a period of 4. When , the TDCS with WFRFT preprocessing can degrade into the traditional TDCS system with the OFDM-based structure in [49], and when , the designed system operates with a DFT-preprocessing, similar to the mechanism in the SC-FDMA system. Certainly, we can utilise the exhaustive searching method or other fast iterative searching methods to get the optimal WFRFT order from 0 to 1.

4. Simulations and Discussion

To test the performance of the designed system in Figure 1 and validate the related theoretical analyses in Section 3, we conduct the following simulations with TDCS block length and Gray-coded QPSK baseband modulation. We assume that the channel is only influenced by AWGN and NBI. As for basic function (BF) in TDCS, the pseudo random phase is generated by a phase mapper from a -sequence generator. And as for NBI in (10), the spreading factor is set to 32, and the central frequency is set to 0.3 normalised block length. obeys Bernoulli distribution, and the initial phase is uniformly distributed in . The threshold for TDCS is set to , and is the maximum spectrum amplitude for NBI.

4.1. Time-Frequency Distribution

To analyse how WFRFT can influence the signal’s distribution in time-frequency domain, we illustrate different residual NBI’s amplitude distribution with different . As shown in Figure 3, for these samples in frequency domain where the corresponding amplitude is larger than the threshold (ZERO interval in each subfigures), the WFRFT preprocessing matrix in (7) can always zero the amplitudes at these carrier frequencies to remove the most strong interference of NBI. As depicted in Figures 3(a)3(d), different in will make the different residual NBI’s distribution in the other nonzero carriers position, which can further suppress the residual NBI to some extent.

4.2. BER Performance with JSR and SNR

To test the designed system’s BER performance of different selections against SNR and ISR, we set SNR in the range dB and ISR in the range dB. And the BER performances are averaged over 1000 Monte Carlo trails. Besides, exhausted searching method is adopted to get the optimal in (19), and optimal BER performance can be achieved at for the current simulation conditions. In Figure 4, the BER performance curves with three typical selections are exhibited to compare the suppression performances on the NBI. As depicted in Figure 4, the system at can greatly outperform the ones at both and with regard to NBI suppression. When ISR is larger enough, for example, when dB, the system with can gain 3-10 dB than the one at and 1 (SNR=8dB, and BER=). This is mainly because can be regarded as the major factor to affect the BER performance in (15), while can be reduced to the most extent at . When dB, there exists small difference between the systems with different selections, and the performance difference can be neglected especially when SNR is small enough (dB). This is mainly because the Gaussian noise has been the major factor to affect the final BER performance, and the NBI’s factor can be omitted compared to the Gaussian noise factor in (15).

4.3. Cross Relationship Coefficient

We also test the normalised cross-correlation coefficients in (13) between residual NBI () and BF () with different central frequency, . As shown in Figure 5, different central frequency will result in different optimisation of . This mainly because the final vectors for BF () in (1) and residual NBI () in (11) will also change when the amplitude vector change with different .

Besides, different vector for residual NBI will further result in quite different time-frequency domain distributions of residual NBI, just shown in Figure 4. The simulation results in Figure 5 show that the optimal WFRFT order should be decided and varies with the different NBI’s central frequency . As shown in Figure 5, especially when , the optimal WFRFT order is 0.3, and the current optimal result is just consistent with the BER performances in Figure 4, which further demonstrates the correctness of the related optimal theory in (19).

5. Conclusions

In this paper, we provide a modified TDCS structure with WFRFT preprocessing and analyse the system’s ability on NBI suppression. Through theoretical deduction, we also derive the system’s BER expression and prove that the best BER’s performance could be obtained by minimising the magnitude of the cross relationship coefficient between TDCS’s basic function vector and residual NBI vector. Compared with the traditional TDCS structure, final simulation results demonstrate that the modified structure with the optimal WFRFT order can exhibit best BER performance to suppress NBI to most extent. As for the optimal program in (19), the exhausted searching method is adopted in the simulation. Thus, in our future works, with regard to WFRFT order selection, the highly efficient searching algorithms need to be explored to achieve the optimal BER performance more effectively. Besides, to further evaluate the proposed system’s practicability, we will also proceed with analysis on the system’s anti-interception ability and antimultipath performances.

Data Availability

The data sets supporting the results of the article are included within the article and its additional files.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant 61271250) and the National Natural Science Foundation of China (Grant 61571460).