1. Introduction
Satellite broadcasting systems DVB-S2/S2X, reverse channels of interactive broadcasting data transmission systems, are characterized by a limited frequency band. Increasing the bit rate in such channels can be achieved by the application of FTN (Faster-Than-Nyquist) signals [
1,
2,
3,
4]. The energy spectra of these signals have a close to rectangular shape. The duration of FTN signals is longer than the duration of the message bit. The correlation properties of signals are determined by the level of intersymbol interference, which is caused by the superposition of neighboring signals on each other [
4,
5,
6]. Methods for constructing FTN signals are described in detail in [
2,
4,
6,
7,
8] and can be divided into two groups: filtering methods and optimization methods of the signal envelope shape [
9,
10,
11,
12,
13,
14]. Their common feature is that the duration of the signals is significantly (8–16 times) greater than the duration of message bits for the binary channel alphabet. The energy efficiency of using such signals is determined by their correlation properties and is characterized by energy costs for various processing algorithms.
Modern trends in increasing the data transmission rate in cases of limited frequency bands are implemented by increasing the bit rate and by increasing the volume of the channel alphabet [
4,
6,
7]. In this case, it is necessary to separate the concepts of the absolute data transmission rate and the transmission rate normalized to the frequency band. Restrictions on the computational complexity of the FTN signal demodulation algorithms are imposed for real-time data transmission systems. At the same time, when the data rate normalized to the band increases, there is no such limitation. An indicator of the effectiveness of the application of FTN signals in such conditions is the approach to the Shannon limit for continuous channels with additive Gaussian noise.
However, even in this case, the issue of computational complexity of FTN signal demodulation algorithms is debated. This is because additional intersymbol interference occurs as the bit rate of FTN signals increases. This interference leads to an increase in energy costs, the specific value of which is determined by the reception algorithms. Moreover, in this regard, most of the attention is paid to developing and obtaining the most effective algorithms, from the point of view of obtaining maximum reception reliability and ensuring minimal computational complexity in their implementation.
It is advisable to compare the reception efficiency under conditions of limited computational complexity of the proposed algorithms. Indeed, it is possible to obtain algorithms that provide low energy losses (for example, Viterbi-type algorithms [
15,
16]) but have high computational complexity. High complexity does not allow data processing in real time, which is often important, for example, when using portable receiving terminals of ground stations in DVB-S2/S2X satellite broadcasting systems. On the other hand, bit-by-bit algorithms are easy to implement, have low computational complexity and have significant energy losses. Moreover, energy losses increase with the increasing bit rate. In this regard, it is advisable to set the task of finding such processing algorithms that, with an increase in the data rate of FTN signals, allow an increase in energy losses of no more than specified for a fixed increase in computational complexity (for example, by more than 30%). Computational complexity can be estimated through the function
, which actually shows the dependence of the amount of calculations performed on the size of the input data/parameters of the algorithm. For example, in the case where the number of addition and multiplication operations (and, as a consequence, running time) depends linearly on the size
n of the input data, the algorithmic complexity can be defined as
. When using an algorithm that is optimal according to the maximum likelihood criterion for receiving “in general” sequences of length
N and signals with a volume of channel alphabet
M, in the general case it is necessary to make
comparisons of the observed implementation with all possible reference combinations of signals [
17]. Then, the computational complexity of such an algorithm can be estimated as exponential:
. In this way, it is possible to quantitatively compare different processing algorithms in terms of energy costs and computational complexity. In addition, it is possible to find the relationship between the increase in the number of addition and multiplication operations and the increase in the bit rate.
Energy costs are the most important in cases of increasing the bit rate for various processing algorithms. The energy loss minimization for FTN signals under such conditions has been the topic of numerous research works [
18,
19,
20,
21]. In these papers, the authors use coherent optimal and suboptimal processing algorithms: maximum likelihood sequence estimation (MLSE), the Viterbi algorithm, and so on. It has been shown that binary data can be transmitted at bit rates approximately 25% above the Nyquist limit without energy losses. The MLSE-based processing algorithm is implemented by enumerating all possible implementations of the FTN signal sequence. Solutions are obtained by finding the minimum Euclidean distance between the resulting implementation of the FTN signal sequence and all possible enumeration implementations. The computational complexity of the MLSE-based processing algorithm limits practical implementation because we need to use processing of the entire transmitted sequence of symbols, which is only possible when scale of the operation of the transmission system is not in real time. Nevertheless, it is possible to take the BER performance in this case as a reference for comparing other bit-by-bit processing algorithms with an assessment of energy losses relative to this algorithm.
The struggle to reduce computational complexity for FTN signal processing has led to the emergence of a large class of bit-by-bit processing algorithms that have been discussed in scientific papers over the past 15 to 20 years. This emphasizes the relevance of solving the problem of demodulating FTN signals using simple methods under conditions of significant intersymbol interference. Among such algorithms, iterative processing algorithms with decision feedback can be distinguished [
18,
19,
22,
23]. There are several modifications of these algorithms, both with a fixed observation bit interval and with an increased observation interval. In this case, the value of the received symbol is re-evaluated depending on the values of subsequent received symbols falling within the observation interval. In other words, the evaluations of previous and subsequent symbols in the message are used to make a decision about a given symbol received.
In [
18], the authors propose two methods for symbol sequence estimation in a received sequence: namely SSSSE (Sequential Character Sequence Estimation) and SSSgbKSE (Sequential Character Sequence Estimation with go-back-K Sequence Estimation). Simulation results show that the proposed methods are suitable for low intersymbol interference scenarios, allowing up to 11% increases in the data rate. However, the authors do not consider the possibility of further increasing the symbol rate beyond this value. The SSSSE and SSSgbKSE algorithms show a promising direction for low-computational-complexity FTN signal processing.
An interesting algorithm for processing FTN signals with a high level of intersymbol interference was proposed in [
24,
25]. This algorithm implements the rule of a decision-directed successive interference cancellation (DDSIC) scheme with Minimum Mean Square Error (MMSE) equalization for FTN signaling. FTN signals after equalization and decision making are regenerated and then subtracted from the original signals at the receiver input to estimate the level of intersymbol interference. This value is used for the next iteration of interference suppression. Based on the simulation results, the BER performance for the DDSIC scheme is better than those described in [
18]. The advantage of this algorithm is its low computational complexity.
In [
25], the authors propose a practical methodology for reducing intersymbol interference in an FTN signal sequence within linear equalization. The framework is based on a low-complexity linear equalizer and serial interference cancellation techniques in the frequency domain. This approach requires additional computational operations in the frequency domain and is characterized by high computational complexity.
Let us consider the energy costs for FTN signals with a binary alphabet as a quantitative illustration of the energy costs in case of increasing the data rate to 42% above the Nyquist barrier.
Table 1 shows the signal-to-noise ratios in dB at the input of the receiving device for various processing algorithms. From this table, we can see that the energy costs are in the range of signal-to-noise ratios of 7–12 dB for
for various algorithms even with a relatively small increase in the bit rate. These values will increase with increases in the data rate. The indicated signal-to-noise ratio values can be used as quantitative guidelines when obtaining new algorithms for processing FTN signals.
Thus, we can identify general trends in improving methods and algorithms for signal processing with an increase in data rates above the Nyquist barrier. The greatest attention is focused on algorithms for bit-by-bit iterative processing with decision feedback, which provide sufficiently high BER performance for receiving FTN signals along with significantly lower computational complexity than with MLSE. Taking these trends into account, we formulated research objectives:
Obtaining an iterative algorithm with decision feedback, using the generalized maximum likelihood criterion, that performs additional functions of estimating the values of subsequent symbols falling within the observation interval, providing a minimum number of computational procedures when determining the choice of maximum numbers in the sequence.
Determining the BER performance of FTN signals by simulation modeling at data rates higher than the Nyquist barrier by 2–3 times, with an assessment of energy costs for transmission channels with additive noise and channels with signal fading.
Determining the computational complexity of the algorithm to find the relationship between the increase in the number of addition and multiplication operations and the increase in the data transfer speed.
Experimentally assessing the level of additional energy losses during packet transmission of sequences of FTN signals due to inaccuracies in the phase and clock synchronization system.
To achieve this goal, it is first necessary to obtain algorithms for receiving FTN signals using an analytical method. In this case, it is necessary to use a criterion that provides the maximum likelihood and a minimum number of addition and multiplication operations as a criterion for the synthesis of the algorithm. As such, a generalized maximum likelihood ratio criterion is chosen, which is used in conjunction with an algorithm with decision feedback (reception of this and previous symbols) to process subsequent symbols falling within the observation interval.
This is discussed in
Section 2,
Section 3 and
Section 4 of this article and in
Appendix A.
Section 5 presents the results of simulation modeling of the BER performance of FTN signals for various channels. The results of experimental studies are presented in
Section 6. In the Conclusion, the main results of the research are presented.
2. Description of the FTN Signaling System
Signal
with the arbitrary form of amplitude pulse
with maximum value
and carrier frequency
is described by the following expression:
where
is the value of the modulation symbol at time interval
, which is equal to a transferring time of one information bit.
A random sequence of
N signals (
1) for an arbitrary value of symbol rate
can be written as
Random sequence (
2) provides the transfer of binary data at rate
. The amplitude pulse
determines the main spectral and correlation properties of random sequence (
2). For signal duration
, intersymbol interference appears even at a transmission rate below the Nyquist barrier. This interference increases significantly when transmitting at a rate that exceeds the Nyquist barrier. The appearance of significant intersymbol interference leads to a degradation in correlation properties and BER performance.
Let us consider a generalized block diagram of the data transmission system based on FTN signals (
Figure 1). A sequence of binary symbols arrives at the transmitter input from the message source at rate
. These binary symbols are converted into modulation symbols of the message in the “Mapper” block. For example, for binary phase modulation,
. A sequence of FTN signals on the basis of the amplitude pulse
is generated in the “FTN modulator”. The structure of this device is described in detail in [
8]. The symbol rate for forming FTN signals is determined by the “Symbol rate
R” block. Sequence of FTN signals (
2) goes to block “Channel” from the transmitter output. Two types of channels are used as a transmission channel for simulation and for experimental studies: a channel with constant parameters and additive Gaussian noise (AWGN)
with an average power spectral density
and a channel with multipath propagation. The random process
arrives at the input of the receiver from the output of the transmission channel. The realization of
is processed in accordance with one or another processing algorithm in the block “FTN demodulator”.
The choice of the reception algorithm is determined by the acceptable energy costs and the computational complexity of the implementation. We propose new bit-by-bit iterative processing algorithms with decision feedback, providing sufficiently high BER performance for FTN signals with limited computational complexity. A possible algorithm for processing the evaluation of the current symbol can be divided into two parts:
The first part is the reception of a given symbol, taking into account decisions about previous accepted symbols (decision feedback), which is implemented with minimal computational complexity;
The second part is the optimal processing of subsequent symbols falling within the observation interval.
Note that we cannot separate these parts of the algorithm. In this case, we make an analytical generation of the algorithm using the generalized maximum likelihood criterion. This general algorithm performs decision feedback processing and estimation of the values of subsequent symbols that fall within the observation interval. Also, this algorithm must provide a minimum number of computational procedures when determining the choice of maximum numbers in the sequence. The choice of this optimality criterion is due to the fact that the minimum possible number of addition and multiplication operations is required in the function of selecting the maximum numbers in the sequence to implement the proposed algorithm.
3. Proposed Algorithm with Decision Feedback with Estimation of Subsequent Symbols in Observation Interval by Maximum Likelihood Ratio
Let us present the process at the input of the receiver in the following form:
where
is a channel transfer coefficient;
and
are the sequences of previous and subsequent signals. Note that
and
depend on
. The indexes
i and
q denote the numbers of specific combinations of previous
I and subsequent
Q signals. For
and
,
and
.
The normalized symbol rate
is equal to
. In this case, the number of previous signals in the observation interval
is equal to
. The number of subsequent signals is equal to
. Then, for
in (
1),
where
and
are the channel symbols of previous (
i) and subsequent (
q) sequences of signals;
p is the ordinal number of the symbol in the sequence of signals. For example, values of
I and
Q are equal to 7 for
and have a symmetrical pulse shape
with duration
. Note that ISI increases significantly at such values of
R and affects the correlation properties.
Let us present
,
and
in terms of in-phase
and
and quadrature
and
baseband components:
Processing algorithms with decision feedback can be obtained provided that the functions
are known and the functions
are processed in accordance with the generalized maximum likelihood criterion. The derivation of this algorithm is given in
Appendix A.
The analytical notation of the algorithm with decision feedback with the evaluation of subsequent symbols in accordance with the maximum likelihood ratio criterion has the following form for
:
where
is the evaluation of the
p-th received symbol.
This algorithm is characterized, as will be shown later, by limited computational complexity with energy costs close to the MLSE algorithm (see
Table 1).
4. Realization of the Processing Algorithm
Analog processing of input process (
8) with transfer to the zero frequency of the spectrum of the sequence of FTN signals is carried out when implementing the preliminary reception algorithm. The radio modem is built on the SDR NI USRP 2920 platform and the signal spectrum is transferred from the carrier frequency to the zero frequency in the analog (hardware) part of the platform. Phase Locked Loop (PLL) and timing are used to achieve coherent reception. PLL and timing are implemented in the platform in the controller and processor in PCHost.
The software part of algorithm (
9) and its block diagram are shown
Figure 2 and
Figure 3. The software part determines the computational complexity of the receiving device. The input of the software part of algorithm (
9) receives digital samples from the function. These samples are distributed among the correlator blocks. The upper correlator (the left side of algorithm (
9)) performs processing on the observation interval
, equal to the duration of the FTN signal. The remaining correlators are used for preliminary assessment of subsequent symbols falling within the observation interval (the right side of inequality (
9)). The introduction of decision feedback when receiving previous symbols occurs in the block “Decision feedback with estimation of subsequent symbols in observation interval by maximum likelihood ratio” and the block “Correlation coefficient of previously received symbols”
in accordance with
in (
9). The correlator block process is performed with digital “sliding” integration by using memory integrators and continuously performing the integration process. This is due to the fact that the duration of the FTN signal exceeds the duration
T, so several consecutive signals will be included in integration processes during integration. Note that a digital matched filter can be used instead of a digital “sliding” integrator. The processing of subsequent signals is performed on the observation interval
in accordance with the maximum likelihood estimate in the digital block “Correlation coefficient of subsequent symbols”
.
Evaluations of the received symbols are formed at the output of the decision unit (DU). The threshold value of the DU is equal to . These values are formed by taking into account the correlation coefficients of the subsequent symbols and the correlation coefficients of the evaluations of the received symbols with the amplitude pulse shape . Blocks “” are needed to calculate the values of the correlation integral . Note that the value of depends on the coefficient , which is determined by the attenuation level of the channel.
The computational complexity of proposed algorithm (
9) is determined by the number of operations for the calculation of the maximum. The upper limit for the number of additions is equal to
and no multiplications, but there are many effective realizations. Note that the complexity of proposed algorithm (
9) does not depend on the length of the message packet. The complexity gain for the comparison with the processing algorithm based on the MLSE is more than two orders in terms of the memory consumption for the case of the length of the received message packet equal to
.
The computational complexity of algorithm (
9) is given in
Table 2. It is determined by the number of additions and multiplications and depends on the number
Q of subsequent symbols which fall within the observation interval. Also, it is determined by the number of samples per bit
. Compared to the MLSE, algorithm (
9) has a lower memory consumption since it does not depend on the packet length. For example, when packet length is
, the gain in memory consumption reaches more than two orders.
We can make a subconclusion that the proposed decision feedback with the estimation of subsequent symbols in the observation interval by the maximum likelihood ratio can be implemented on the basis of a universal digital platform using universal microprocessors (e.g., starting from family Intel 3000) as a PCHost. We obtained the values of the computational complexity of the proposed algorithm depending on signal parameters (signal duration, level of ISI). Using these values, it is possible to estimate the increase in the computational complexity under conditions of fixed energy losses and an increased transmission rate. In conditions of specific transmission rates up to 3 bps/Hz and signal duration (8–16)T, the computational complexity of the proposed algorithm is more than two orders lower compared to that of MLSE.
6. Simulation Results—Comparison of Energy Efficiency and Computational Complexity of Decision Feedback Detection Algorithm with Estimation of Subsequent Symbols in Observation Interval by Maximum Likelihood Ratio
Let us estimate the increase in the computational complexity of decision feedback algorithm (
9) caused by an increase in data transmission rate. We try to answer the following questions:
To what extent must the computational complexity of the digital demodulator of spectrally efficient signals be increased under conditions of a fixed BER and increased symbol rate?
What is the value of the energy loss caused by an increased symbol rate under conditions of a fixed BER and computational complexity?
By the energy loss, we mean the difference between the signal-to-noise ratios
in the case of the ISI-free BER performance and in the case of using proposed algorithm (
9) at
.
Then, it is possible to construct dependences between the energy losses of algorithm (
9) in relation to the potential BER depending on the number of addition and multiplication operations for different bit rates.
Figure 8 shows the values of the energy loss for the bit rates from
to
.
Figure 8 shows the energy loss values obtained as a result of simulation modeling for bit rates from
to
for FTN signals, the parameters of which are given in
Section 5. In this figure, the X-axis shows the number of addition and multiplication operations (computational complexity). Moreover, each point along this axis corresponds to the value of parameter
q in formula (
9). This parameter determines the number of subsequent signals falling within the observation interval, which are taken into account when implementing the scheme in
Figure 5. For example, for
and
(one previous symbol is used in the feedback), the energy loss is equal to 0.6 dB at the minimum values of the computational complexity. When
(four previous values of symbols are used in the feedback ), the energy loss practically becomes zero, that is, the BER coincides with the potential one. By increasing the transmission rate to
, approaching the potential BER is achieved by increasing the computational complexity, which reaches a value of
addition and multiplication operations.
With a further increase in the bit rate up to , the minimum energy loss will be equal to 1.2 dB with the number of operations . When the bit rate is reached, the number of operations increases up to with a minimum energy loss of 3.5 dB.
Thus, the following conclusions can be made for proposed processing algorithm (
9) with decision feedback. First, it is possible to estimate resulting energy losses from an increase in the bit rate at a fixed error probability (
) using the developed methodology. Second, the energy loss in relation to potential BER increases by approximately 2 dB with an increase in the bit rate speed by 1.5 times (above the Nyquist barrier). Third, an increase in the bit rate by 1.5 times leads to an increase in the energy loss from 1.2 dB to 3.7 dB for a fixed
with a fixed computational complexity (which is often the main criterion for the practical implementation of a digital demodulator). Finally, an increase in the bit rate by 1.5 times leads to an increase in computational complexity from
to
numbers of operations for acceptable increases in energy loss equal to 2 dB. Summarizing this conclusion, we can estimate how much additional increases in the number of computational operations will cost a gain in energy efficiency, for example, 1 dB. For the proposed processing algorithm with decision feedback, each decibel is achieved by increasing the number of computational operations by
.
8. Conclusions
We proposed a new bit-by-bit algorithm, (
9), for processing FTN signals designed for bit rates above the Nyquist barrier. The advantage of this algorithm is BER performance close to MLSE with minimal computational complexity. Algorithm (
9) is based on an iterative algorithm with decision feedback. It is complemented by an algorithm for receiving subsequent signals falling within the observation interval, which is optimal according to the generalized maximum likelihood criterion, providing a minimum number of computational procedures when determining the choice of maximum numbers in the sequence. It is shown that the computational complexity of the proposed algorithm is more than two orders less than the complexity of the MLSE algorithm.
The energy costs for receiving FTN signals at bit rates above the Nyquist barrier for transmission channels with additive noise (AWGN) and fading channels for proposed algorithm (
9) are given in
Table 7. This table contains the results of simulation and the results of experimental research.
In this paper, we give a numerical estimation of the increase in computational complexity for proposed algorithm (
9) at bit rates exceeding the Nyquist barrier. An increase in the bit rate by 1.5 times (from
to
) leads to an increase in energy losses from 1.2 dB to 3.7 dB in the area of
and under conditions of fixed computational complexity. If we fix the possible energy losses, for example, at a value of 2 dB, then an increase in the bit rate by 1.5 times (from
to
) leads to an increase in computational complexity by more than an order, from
to
operations. The same conclusion can be reformulated in another form: for proposed algorithm (
9), each decibel of energy gain is achieved by increasing the number of computational operations by
operations.
The application of proposed algorithm (
9) in the DVB standard of satellite broadcasting systems S2/S2X can provide increases in the bit rate without significant degradation in the BER with minimal computational complexity of the demodulator. This is especially important for mobile terminals of satellite broadcasting ground stations. Potential future directions for this research are related to the optimization of pulse shapes for different channels to achieve additional increases in spectral and energy efficiency without increasing computational complexity.