Next Article in Journal
Simulation Modeling for Ships Entering and Leaving Port in Qiongzhou Strait Waters: A Multi-Agent Information Interaction Method
Next Article in Special Issue
Estimation of Source Range and Location Using Ship-Radiated Noise Measured by Two Vertical Line Arrays with a Feed-Forward Neural Network
Previous Article in Journal
An Anti-Occlusion Approach for Enhanced Unmanned Surface Vehicle Target Detection and Tracking with Multimodal Sensor Data
Previous Article in Special Issue
Underwater Mapping and Optimization Based on Multibeam Echo Sounders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Channel Estimation and Iterative Decoding for Underwater Acoustic OTFS Communication Systems

1
China Ship Scientific Research Center, Wuxi 214000, China
2
National Key Laboratory on Ship Vibration & Noise, Wuxi 214000, China
3
Taihu Laboratory of Deepsea Technological Science, Wuxi 214000, China
4
National Key Laboratory of Underwater Acoustic Technology, Harbin Engineering University, Harbin 150001, China
5
Key Laboratory of Marine Information Acquisition and Security, Harbin Engineering University, Ministry of Industry and Information Technology, Harbin 150001, China
6
College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, China
7
Sanya Nanhai Innovation and Development Base of Harbin Engineering University, Sanya 572024, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(9), 1559; https://doi.org/10.3390/jmse12091559
Submission received: 7 August 2024 / Revised: 2 September 2024 / Accepted: 3 September 2024 / Published: 5 September 2024
(This article belongs to the Special Issue Applications of Underwater Acoustics in Ocean Engineering)

Abstract

:
Orthogonal Time–Frequency Space (OTFS) is an innovative modulation method that ensures efficient and secure communication over a time-varying channel. This characteristic inspired us to integrate OTFS technology with underwater acoustic (UWA) communications to counteract the time-varying and overspread characteristics of UWA channels. However, implementing OTFS in UWA communications presents challenges related to overspread channels. To handle these challenges, we introduce a specialized OTFS system and offer frame design recommendations for UWA communications in this article. We propose a Doppler compensation method and a dual-domain joint channel estimation method to address the issues caused by severe Doppler effects in UWA communication. Additionally, we propose an OTFS system detection approach. This approach incorporates an iterative detection process which facilitates soft information exchange between a message passing (MP) detector and a low-density parity check (LDPC) decoder. By conducting simulations, we demonstrate that the proposed UWA OTFS system significantly outperforms Orthogonal Frequency-Division Multiplexing (OFDM), Initial Estimate Iterative Decoding Feedback (IE-IDF-MRC), and two-dimensional Passive Time Reversal Decision Feedback Equalization (2D-PTR-DFE) in UWA channels.

1. Introduction

Integrated sensing technology for space–air–ground–ocean (SAGO) contexts is a key research focus for next-generation mobile communications [1]. Underwater communication technology plays a crucial role in SAGO integrated sensing and information transmission systems. Acoustic waves are the sole effective medium for long-range underwater information transmission. However, the underwater acoustic (UWA) channel has numerous drawbacks, including limited bandwidths, significant delay spreads, frequency-selective fading, complex time variability, and susceptibility to the Doppler effect [2,3].
Orthogonal Frequency Division Multiplexing (OFDM) is a multi-carrier modulation technology in which the frequencies of multiplexed subcarriers are orthogonal to each other. This orthogonality ensures that the peak value of each subcarrier aligns with the zero value of other subcarriers. Consequently, even if subcarriers overlap in the spectrum, they do not interfere with each other at the peak of the desired subcarrier. Moreover, OFDM systems avoid bandwidth efficiency loss, which is common in systems using non-orthogonal carrier sets. Compared to earlier systems, OFDM significantly enhances the spectral efficiency of UWA communication systems. Its low implementation complexity in terms of modulation and demodulation has led to its widespread adoption in high-speed UWA communication and network systems [4,5,6,7,8,9,10,11,12].
However, in high-speed motion scenarios, rapid movement result in Doppler frequency shifts and expansions within the communication system [13]. Additionally, the carrier frequency offset (CFO) present in received signals can impair the orthogonality among subcarriers in OFDM systems, resulting in inter-carrier interference (ICI) which severely compromises system performance [14]. Existing approaches aim to mitigate this phenomenon through shortening the OFDM symbol duration to reduce channel variations within each symbol. However, this method’s primary drawback is reduced spectral efficiency due to symbol duration shortening [15]. Consequently, OFDM modulation is particularly vulnerable to Doppler spread when employed in UWA communication networks, resulting in severe system performance degradation. Addressing efficient data transmission in dynamic and mobile UWA communication scenarios is thus a critical challenge. Enhancing the spectral efficiency in UWA communication networks—particularly in high-Doppler and high-delay UWA channels—has emerged as a primary research focus for future advancements in UWA communication technology.
In 2017, Hadani et al. introduced a novel two-dimensional modulation technology termed Orthogonal Time–Frequency Space (OTFS) [16]. OTFS modulation presents a significant advancement in addressing the challenges posed by severe Doppler effects and channel impairments in high-speed mobile communication systems. Unlike traditional OFDM modulation in the time–frequency (TF) domain, OTFS operates in the delay-Doppler (DD) domain for information modulation. The DD domain exhibits robust tolerance to both time delays and Doppler shifts, offering the potential for full diversity—a key factor supporting reliable communication. Furthermore, OTFS modulation effectively transforms a time-varying fading channel into a stationary and non-fading channel in the DD domain, thereby maintaining nearly constant channel gains across all modulation symbols [16].
OTFS can be conceptualized either as an alternative multi-carrier modulation scheme or as a two-dimensional spread spectrum technique applied to QAM symbols, which can subsequently be used with OFDM or other multi-carrier modulation schemes such as Generalized Frequency Division Multiplexing (GFDM) for modulation. The integration of OTFS technology in UWA communication systems addresses the performance degradation observed in OFDM technology when operating in complex and rapidly changing UWA channels. However, as most existing UWA communication systems are designed for low mobility and low carrier frequency scenarios, OTFS technology introduces new challenges when designing the architectures of and algorithms for transmitters and receivers.
To fully leverage the potential of OTFS technology, several critical issues need to be addressed, including waveform design, channel estimation, detection techniques, and multi-antenna and multi-user designs. Solving these challenges is imperative for realizing the benefits of OTFS technology in UWA communication systems.
To date, research on OTFS modulation technology has primarily been concentrated in the domain of wireless communication, and it remains at an early stage of development. Moreover, there is a paucity of studies exploring the application of OTFS in the field of UWA communication [17,18,19,20,21,22,23]. Nevertheless, some scholars have begun investigating its applicability in UWA communication through leveraging insights from wireless communication research, yielding promising preliminary results. These efforts have primarily been centered on reducing receiver complexity, enhancing channel estimation, and improving signal detection; the research has largely been confined to simulation analyses and field experiments, with no significant progress towards engineering implementation or real-world application.
As an emerging modulation technique, OTFS technology has demonstrated significant potential, warranting further investigation to enhance the performance of UWA communication systems. Although some studies have explored the application of OTFS in UWA communication, the body of literature remains limited and lacks diversity in research directions. The challenge of effectively integrating the performance advantages of OTFS technology in the context of wireless communication into the UWA communication domain remains an area requiring more comprehensive research. This study aims to further combine the OTFS method with UWA communication in order to design OTFS systems based on the unique characteristics of UWA channels. The main contributions of this study are as follows:
For the OTFS system transmitter, frame structures tailored to the characteristics of UWA channels are designed by selecting a zero-padding (ZP) frame structure design method. Additionally, the size of the DD domain grid and the length of the guard interval are optimized to enhance the communication efficiency while maintaining performance. The proposed OTFS frame design can significantly reduce the bit error rate (BER) while maintaining the same data rate through careful frame optimization.
For the OTFS system receiver, we propose a Doppler compensation method and a dual-domain joint channel estimation method to tackle the issues caused by the strong Doppler effect in UWA communication. Additionally, we propose an OTFS system detection approach based on an iterative decoding feedback-message passing detector (IDF-MP). Based on careful consideration of the characteristics of UWA channels, the proposed OTFS system receiver is shown to provide improvements, when compared with the existing methods of OFDM, Initial Estimate Iterative Decoding Feedback (IE-IDF-MRC) and two-dimensional Passive Time Reversal Decision Feedback Equalization (2D-PTR-DFE), in UWA channels.
The remainder of this paper is organized as follows: Section 2 introduces the related work. Section 3 describes the UWA OTFS system, including OTFS modulation and the UWA OTFS channel model. In Section 4, we propose a comprehensive OTFS system design scheme based on the characteristics of UWA channels. In Section 5, we analyze the communication performance of the proposed UWA OTFS system using simulated UWA channels. Finally, Section 6 presents the concluding remarks and discusses directions for further research.

2. Related Work

To date, research on OTFS modulation technology has primarily focused on the realm of radio communication, and it remains in a nascent stage at present [15,16,17,24,25,26,27,28]. In recent years, numerous enhancements to the original OTFS modulation scheme have been proposed by researchers. For example, Raviteja et al. [15,16,24,25] established the mathematical framework of OTFS signals and derived the modulation and demodulation input–output equations in the DD domain. In [16], an OTFS modulation model for multi-channel dual dispersion channels was established, accounting for time delays within the maximum delay and limited Doppler shift range. However, the original OTFS input–output relationship assumes ideal pulse shaping waveforms to be biorthogonal in both time and frequency, which is unattainable due to the Heisenberg uncertainty principle [16]. To address this, the authors of [24] derived the input–output relationship for actual pulse-shaping waveforms based on the DD domain, extending it to any waveform. Moreover, in [15], it was demonstrated that using the message passing (MP) algorithm for symbol detection in OTFS results in a lower symbol error probability (SEP) compared to OFDM. Building on this, the authors of [24] proposed a scheme that adds a cyclic prefix (CP) only at the end of the OTFS frame, significantly reducing the overhead caused by the loss of biorthogonality of pulse shaping waveforms. Additionally, in [25], a channel estimation method utilizing embedded pilot symbols was introduced, showing that OTFS with non-ideal channel estimation outperforms OFDM with ideal channel estimation in terms of error performance. Further research has delved into the sparsity of the DD domain in OTFS systems [26,27,28]. For example, Tharaj Thaj et al. proposed a low-complexity iterative decision feedback equalizer (DFE) based on maximum ratio combining (MRC), demonstrating its favorable performance and complexity compared to state-of-the-art message passing detectors [17]. These advancements represent significant progress in enhancing OTFS technology, paving the way for its broader application in communication systems.
Compared with the field of terrestrial wireless communication, research on OTFS modulation technology in the field of UWA communication remains relatively limited [17,18,19,20,21,22,23]. Some scholars have explored its applicability in UWA communication through leveraging research findings from wireless communication and have achieved notable results, primarily focusing on receiver complexity reduction, channel estimation, and signal detection. Researchers at Xiamen University have proposed a single-input single-output (SISO) UWA communication scheme based on OTFS [19]. Simulation results comparing OTFS and OFDM in UWA channel environments in terms of bit error rate (BER), spectral efficiency, and peak-to-average power ratio (PAPR) indicated the superior performance of the OTFS-based UWA communication scheme in time-varying multipath UWA channels, compared to OFDM and discrete Fourier transform spread (DFTS) OFDM schemes. In 2019, Bocus et al. from the University of Bristol, U.K., introduced an OTFS system based on large-scale multiple-input multiple-output (MIMO) OFDM to investigate multi-user UWA communication schemes [20]. In the studies [19,20], OTFS-UWA communications in SISO and MIMO scenarios were explored, respectively, with simulation results demonstrating that OTFS outperforms OFDM. However, these studies assume that accurate channel information is known to the receiver. Researchers from Dalian University of Technology have proposed a time-reversal OTFS receiver with adaptive channel estimation [21] and evaluated its system performance in both analog and measured channels. Their results indicated that the proposed two-dimensional passive time reversal (PTR) receiver achieves improved system performance with low complexity. Su Hang et al. have designed a practical UWA OTFS system, incorporating Doppler compensation and channel estimation in DD domain [22]. Additionally, Wei Li et al. have discussed input–output relationships for OTFS under overspread channels while providing a frame design method tailored to the characteristics of UWA channels [23].
In summary, an increasing number of scholars have studied the performance of OTFS in UWA communication. However, most of these studies concentrate on theoretical applications in radio communications, often assuming that the channel is underspread, which is typically not the case for UWA channels that are frequently overspread. Moreover, most UWA OTFS system designs assume perfect channel estimation or only consider the influence of integer Doppler factors. While theoretical investigations into fractional Doppler shifts have been conducted, practical implementation remains challenging. Neglecting fractional Doppler shifts in UWA channels can lead to non-negligible inter-Doppler interference (IDI), severely compromising performance, particularly in UWA communication systems characterized by significant Doppler effects. Therefore, developing a practical OTFS design that is suitable for real-world UWA communication scenarios is imperative. Filling this gap was a key motivation of this work.

3. Underwater Acoustic OTFS System Model

3.1. OTFS Modulation

We assume that one UWA OTFS frame has M subcarriers and N symbols. Then, the bandwidth, signal duration, and symbol number are M Δ f , N T , and M N , respectively, where Δ f is the subcarrier spacing and T = 1 / Δ f is the symbol duration.
We place OTFS symbols x ( q , l ) into DD grids (size: M × N ) with delay resolution 1 / ( M Δ f ) and Doppler resolution 1 / ( N T ) , where q = 0 , 1 , , N 1 and l = 0 , 1 , , M 1 . After inverse symplectic finite Fourier transform (ISFFT), OTFS symbols x ( q , l ) in the DD domain are transformed into X ( n , m ) in the time–frequency domain:
X ( n , m ) = 1 M N q = 0 N 1 l = 0 M 1 x ( q , l ) e j 2 π ( n q N m l M ) .
Then, X ( n , m ) is transformed into the passband time domain signal s ˜ ( t ) after Heisenberg transformation and up-conversion:
s ˜ ( t ) = s ( t ) e j 2 π f c t = 1 M n = 0 N 1 m = 0 M 1 X ( n , m ) g ( t n T ) e j 2 π m Δ f ( t n T ) e j 2 π f c t ,
where g ( t ) is the rectangular window function and f c is the carrier frequency.

3.2. UWA OTFS Channel Model

In the DD domain, an UAW channel can be modeled as follows:
h ( τ , υ ) = i = 1 p h i δ ( τ τ i ) δ ( υ υ i ) ,
where h i , τ i , and υ i are the path gain, the path delay, and the Doppler shift of the ith path, respectively.
Then, we can obtain the delay shift and Doppler shift in the DD grid for the ith path:
τ i = l i M Δ f ,
υ i = q i N T = q υ i + q ˜ υ i N T ,
where l i is an integer delay tap when M is large enough and q υ i and q ˜ υ i are the integer part and the fractional part of Doppler shift, respectively, with q ˜ υ i ( 1 2 , 1 2 ] .
Based on the resolution in the DD domain, we define τ max < 1 Δ f and υ max < 1 T as the maximum delay and Doppler spread, respectively. Most OTFS system designs suppose that the channel is underspread [22,23]:
τ max υ max < < 1 .
Based on the underspread assumption, l i < M and N 2 < q i < N 2 must be met. The passband OTFS signal after passing through the UWA channel can be represented as follows:
r ˜ ( t ) = h ( τ , υ ) s ˜ ( t τ ) e j 2 π υ ( t τ ) d τ d υ + ω ˜ ( t ) = i p h i s ( t τ i ) e j 2 π f c ( t τ i ) e j 2 π υ i ( t τ i ) + ω ˜ ( t ) ,
where ω ˜ ( t ) denotes white Gaussian noise.
After down-conversion to baseband r ( t ) = r ˜ ( t ) e j 2 π f c t , we obtain the following:
r ( t ) = i p h i e j 2 π σ f c t s ( t τ i ) e j 2 π υ i ( t τ i ) + ω ( t ) ,
where h i = h i e j 2 π f c τ is the complex path gain and the exponential term e j 2 π σ f c t is caused by residual carrier frequency offset (CFO). After Wigner transform and SFFT, we can obtain the received symbols y ( q , l ) in the DD domain.

4. Underwater Acoustic OTFS System Design

This part introduces the design of the proposed OTFS system, as shown in Figure 1. It includes the design of the UWA OTFS frame (blue dashed line) and the UWA OTFS receiver (red and green dashed lines).

4.1. Frame Design

In order to avoid inter-block interference, we consider a zero-padded OTFS system design for UWA communications, as shown in Figure 1 [17]. The blue dashed line represents the design of ZP-OTFS frames. Setting the last 2 l g + 1 symbol vectors (rows) in the DD domain grid to zeros, we obtain an interleaved ZP guard band in the time domain. Pilot signals for channel estimation are inserted into the ZP section without additional data overhead. The ZP-OTFS frame in DD domain grid is defined as follows:
x ( q , l ) = x d ( q , l ) , 0 l < M ( 2 l g + 1 ) , 0 q N 1 0 , M ( 2 l g + 1 ) l M 1 , 0 q N 1 x p ( q , l ) , l = m p = M ( l g + 1 ) , q = n p [ 0 , N 1 ] ,
where x d ( q , l ) represents the OTFS data symbol at the lth row and qth column in the Doppler grid, x p ( q , l ) is the pilot symbol at the m p th row and n p th column, and 0 signifies the guard symbol outside of x d ( q , l ) and x p ( q , l ) in the Doppler grid. Assuming that l g is sufficiently long, channel effects on x d ( q , l ) and x p ( q , l ) will only extend into the ZP section, thereby effectively avoiding inter-block interference.

4.1.1. Design of M and N

In this study, the channel DD response is assumed to be constant, that is, the parameters h i , τ i , and υ i in Equation (3) are constant within one frame period [23]. Based on this assumption, the OTFS frame length T f cannot exceed the channel coherence time, which is typically short for a UWA channel. Consequently, N and M should not be excessively large simultaneously. Thus, we must carefully select N and M to ensure the robustness of UWA communication.
Typically, in conventional OTFS systems, a large M is chosen to meet the delay resolution requirements, resulting in a smaller N with bandwidth B = M Δ f and frame length T f = N T . This often leads to M > > N . However, in UWA communication, the system’s bandwidth is limited and the Doppler spread υ max is often significant. Consequently, a small N cannot satisfy the Doppler resolution requirement for an underspread condition [23].
Accordingly, we must set a suitable value of M value to avoid significant fractional delays, and appropriately increase N to enhance the system’s Doppler tolerance. Simultaneously, we must ensure that the frame duration T f remains within the channel coherence time to maintain reliable communication.

4.1.2. Design of l g

Zero-padding can completely avoid inter-block interference when l g > l m a x but comes at the cost of a very low communication rate. Thus, the means by which to balance the communication rate and robustness need to be carefully considered. As shown in Figure 2, which displays an example of the UWA channel power response, we can define a maximum additional delay τ β as the time at which the multi-path power decays from its initial value by a factor β dB, defined as follows:
10 lg P max P τ β d B = β ,
where l β dB is the delay tap based on τ β dB . Later lower-power paths can be ignored, due to their minimal impact.
In brief, we can set l g = l β dB to enhance the UWA communication rate and ensure communication robustness.

4.2. Synchronization and Doppler Compensation

The proposed UWA OTFS system employs a linear frequency modulation (LFM) signal for preamble detection and synchronization. However, relying solely on this preamble is inadequate for a UWA OTFS system designed for long delay spreads and significant Doppler effects, due to the requirement that the channel be underspread for OTFS to operate effectively. Therefore, a Doppler compensation method is proposed [29,30].
First, we resample the received signal as y ( t ) = r ( t 1 + a ^ ) , where a ^ is the general Doppler scaler:
y ( t ) = i p h i e j 2 π σ f c 1 + a ^ t s ( t τ i 1 + a ^ ) e j 2 π υ i ( t τ i ) + ω ( t 1 + a ^ ) ,
where υ i = υ i 1 + a ^ and τ i = τ i ( 1 + a ^ ) .
Second, we multiply y ( t ) by e j 2 π ε ^ t to remove the CFO, where ε ^ is the residual Doppler shift.
We estimate a ^ and ε ^ according to [31]. Through Doppler compensation, the influence of residual Doppler is significantly reduced, allowing the underspread channel assumption of Equation (6) to be satisfied.

4.3. Dual-Domain Joint Channel Estimation

Here, we propose a dual-domain joint channel estimation method, as illustrated in Figure 3. This method first estimates the path delay in the DD domain, then estimates the complete channel impulse response (CIR) in the DT domain based on the path delay estimation results.

4.3.1. Path Delay Estimation in DD Domain

As fractional Doppler shifts do not impact the estimation of delay taps, we can ignore the IDI interference caused by fractional Doppler shifts. The relationship between y ( q , l ) and x ( q , l ) with integer Doppler shifts was derived in [15], as follows:
y ( q , l ) = q = q max q l = 0 l g [ d ( q , l ) ξ q , l x ( ( q q ) mod N , ( l l ) mod M ) ] + w ( q , l ) ,
where q m a x is the range of fractional Doppler coefficients after Doppler compensation, ξ q , l = h ( q , l ) η ( q , l ) . d ( q , l ) { 0 , 1 } is the complex channel gain with the phase η ( q , l ) and the amplitude h ( q , l ) = h ( q , l ) e j 2 π q l N M , and d ( q , l ) { 0 , 1 } is the path indicator; for example, d ( q , l ) = 1 indicates that there is a path with Doppler tap q and delay tap l with corresponding path magnitude ξ q , l . Furthermore, w ( q , l ) denotes Gaussian noise and ‘(.)mod’ represents the modulus.
The relationship in Equation (12) can be transformed into vector form as follows:
Y p = X p h + w ,
in which:
h = [ ξ q max , 0 , ξ q max , 1 , , ξ q max , l g , ξ q max + 1 , 0 , ξ q max + 1 , 1 , , ξ q max , l g ] T ,
Y p = [ y ( 0 , m p ) , y ( 0 , m p + 1 ) , , y ( 0 , m p + l g ) , y ( 1 , m p ) , y ( 1 , m p + 1 ) , , y ( N 1 , m p + l g ) ] T ,
X p = [ x q max , 0 , x q max , 1 , , x q max , l g , x q max + 1 , 0 , x q max + 1 , 1 , , x q max , l g ] T ,
and
x q , l = x ( ( 0 q ) mod N , ( m p l ) mod M ) x ( ( 0 q ) mod N , ( m p + 1 l ) mod M ) x ( ( 0 q ) mod N , ( m p + l g l ) mod M ) x ( ( 1 q ) mod N , ( m p l ) mod M ) x ( ( 1 q ) mod N , ( m p + l g l ) mod M ) x ( ( N 1 q ) mod N , ( m p + l g l ) mod M ) ,
where h is the equivalent channel with sparsity, Y p is the received pilot signal, and X p is the sensing matrix, for l [ 0 , l g ] and q [ q max , q max ] .
Next, we can estimate the equivalent channel h using a compressive sensing algorithm, such as the orthogonal matching pursuit (OMP) method [29,32]. Then, the set of path delays { l ^ } is obtained according to the index of h .

4.3.2. CIR Estimation in DT Domain

According to [17], the channel impulse response in the DD domain in Equation (3) can be rewritten as follows:
h ( τ , υ ) = l L q Q l h ( l , q ) δ ( τ l T / M ) δ ( υ q Δ f / N ) ,
where L = { l i } is the set of distinct normalized delay shifts in the DD domain, Q l = { q i | l = l i } is the set of normalized Doppler shifts for each path with normalized delay shift l i , and the lth delay tap Doppler response is
h ( l , q ) = h i , if l = l i a n d q = q i 0 , otherwise
Then, the UWA CIR is as follows:
h ˜ ( τ , υ ) = υ h ( τ , υ ) e j 2 π υ ( t τ ) d υ = q Q l h ( l , q ) e j 2 π q Δ f N ( t τ ) .
Discretizing the channel after sampling at t = a T / M , 0 a N M 1 :
h ˜ ( l , a ) = q Q l h ( l , q ) e j 2 π q ( a l ) M N .
Setting a = m + n M , according to [17], the input–output relationship with the channel in the DT domain can be defined as follows:
y ( m , n ) = l L h ˜ ( l , m + n M ) x ( m l , n ) ,
where x ( m , n ) and y ( m , n ) are the transmitted and received symbols in the mth row, respectively, with size N × 1 , as shown in Figure 1.
Given an obtained path delay l ^ , let m = m p + l ^ , l ^ L , the input–output relationship of pilot signal can be expressed as follows:
y ( m p + l ^ , n ) = l L h ˜ ( l , m p + l ^ + n M ) x ( m p + l ^ l , n ) .
Based on Equation (9), x ( m p + l ^ l , n ) = 0 when 0 < | l l | l g (i.e., 0 < | l l ^ | l g ), due to the zero-padding. Thus, the CIR can be obtained as follows:
h ^ ( l , m p + l ^ + n M ) = y ( m p + l ^ , n ) x ( m p , n ) , f o r l ^ L ,
where x ( m p , n ) is the transmitted pilot signal. The UWA CIR in the DT domain can be represented using spline interpolation [22]:
h ^ ( l , a ) = s p l i n e _ i n t e r p o l a t e ( { h ^ ( l , m p + l ^ + n M ) | n [ 0 , N 1 ] } , [ m p + l ^ : M : N M ] , [ 0 : N M 1 ] ) .
Finally, based on h ^ ( l , a ) , the DD domain channel matrix H can be obtained according to [17].

4.4. Iterative Decoding Feedback-Based MP Detector

The method proposed in this study introduces an OTFS system detection approach based on an iterative decoding feedback message passing detector (IDF-MP). This method enhances the detection performance through an iterative process that enables soft information exchange between the MP detector and a low-density parity check (LDPC) decoder. It utilizes symbol-level soft information output by the detector and bit-level soft information output by the LDPC decoder to improve the information transmission efficiency.
As shown in Figure 4, the process commences with the utilization of the channel equivalent matrix obtained from the preceding channel estimation results. Subsequently, the OTFS DD domain signal and the channel equivalent matrix are fed into the MP detector. The MP algorithm utilized in this study, which has been previously detailed in [15], is specifically designed to leverage the sparse characteristics of the channel matrix H in the DD domain. We do not further elaborate on this algorithm in this paper.
Figure 4 illustrates the interconnection between the MP detector and the LDPC decoder within the iterative structure. The conversion of symbol-level soft information and bit-level soft information between the MP detector and LDPC decoder is facilitated by log-likelihood ratio (LLR) converters, as well as interleaving and de-interleaving sections.
In the iterative structure of the IDF-MP algorithm, there are three forms of loop counts:
The outer-loop iterations n o u t .
The inner-loop iterations of the MP detector n i n , M P .
The inner-loop iterations of the LDPC decoder n i n , L D P C .
The specific information transfer process between the MP detector and LDPC decoder is described in the following.

4.4.1. From the MP Detector to the LDPC Decoder

Within the n o u t outer loop, when the inner-loop iterations of the MP detector reach the maximum iterations N i n , M P , the symbol-level outer information output by the detector is as follows:
M e , M P n o u t ( x i = a j ) = p i N i n , M P ( x i = a j ) ,
where x i represents the ith symbol after OTFS modulation ( i = 1 , 2 , , N M ), A represents the set of modulation symbols, a j A ,   j = 1 , 2 , , Q , and Q is the size of the modulation symbol set A (i.e., Q = 16 when the 16 QAM is used for modulation). Furthermore, p i N i n , M P ( x i = a j ) denotes that when the inner-loop iterations of the MP detector reaches the maximum number of iterations N i n , M P , the probability of the symbol a j in the quadrature amplitude modulation (QAM) alphabet is calculated within the MP detector. This probability is the symbol-level output of the MP detector.
The symbol-level soft information output by the MP detector M e , M P n o u t ( x i = a i ) passes through the LLR converter to obtain the bit-level soft information.
L e , M P n o u t ( c i , k = a j , k ) = l n x i A k 1 M e , M P n o u t ( x i = a j ) x i A k 0 M e , M P n o u t ( x i = a j ) ,
where a j , k represents the bit value (0 or 1) of the kth bit in the jth element of the modulation symbol set.
Each modulation symbol a j A corresponds to a binary sequence with a length of l o g 2 | Q | . Therefore, c i , k represents the bit value (0 or 1) on the kth subcarrier of the ith symbol, where k = 0 , 1 , , l o g 2 | Q | . A k 1 ( A k 0 ) represents the modulation symbol set in which the kth subcarrier in the modulation symbol set A is 1 ( 0 ) .
The bit-level soft information L e , M P n o u t ( c i , k = a j , k ) output from the LLR in the n o u t outer loop is de-interleaved to obtain the input information L p r , L D P C n o u t ( c i , k = a j , k ) of the LDPC decoding part. L p r , L D P C n o u t ( c i , k = a j , k ) serves as the prior information of the LDPC decoder in the n o u t outer loop.

4.4.2. From the LDPC Decoder to the MP Detector

Within the n o u t outer loop, the bit-level soft information output from the LDPC decoder is transmitted when the outer-loop iterations of the LDPC decoder reach the maximum number of iterations N i n , L D P C .
L e , L D P C n o u t ( c i , k = a j , k ) = L p o , L D P C n o u t ( c i , k = a j , k ) L p r , L D P C n o u t ( c i , k = a j , k ) ,
where L p o , L D P C n o u t ( c i , k = a j , k ) is the cumulative probability value. Then, after passing through the interleaver, the bit-level soft information L e , L D P C n o u t ( c i , k = a j , k ) is transformed into the input information L p r , M P n o u t ( c i , k = a j , k ) of the detector. The prior probability value is calculated using Equations (29) and (30):
p ( c i , k = 0 ) = 1 1 + exp ( L p r , M P n o u t ( c i , k = a j , k ) ) ,
p ( c i , k = 1 ) = 1 1 + exp ( L p r , M P n o u t ( c i , k = a j , k ) ) .
We introduce the variable λ i , k and represent Equations (29) and (30) uniformly as shown in Equation (31). When c i , k = 0 , λ i , k = 1 ; meanwhile, when c i , k = 1 , λ i , k = 1 .
p ( c i , k = a j , k ) = 1 1 + exp ( λ i , k L p r , M P n o u t ( c i , k = a j , k ) ) = exp ( λ i , k L p r , M P n o u t ( c i , k = a j , k ) 2 ) exp ( L p r , M P n o u t ( c i , k = a j , k ) 2 ) + exp ( L p r , M P n o u t ( c i , k = a j , k ) 2 ) .
After passing through the LLR converter, the bit-level probability value L p r , M P n o u t ( c i , k = a j , k ) is converted into the probability value M p r , M P n o u t ( x i = a j ) for each modulation symbol. This probability value is then utilized as prior probability information input to the detector, thereby completing an iteration.
M p r , M P n o u t ( x i = a j ) = p ( x i = a j ) = Π k = 1 log 2 | Q | p ( c i , k = a j , k ) = Π k = 1 log 2 | Q | exp ( λ m , k L p r , M P n o u t ( c i , k = a j , k ) 2 ) exp ( L p r , M P n o u t ( c i , k = a j , k ) 2 ) + exp ( L p r , M P n o u t ( c i , k = a j , k ) 2 ) .
Finally, we repeat the iteration process until the maximum number of iterations N o u t is reached. Then, the LDPC decoder outputs bit-level LLR information L e , L D P C N o u t ( c i , k ) to make a decision and outputs bit information. The LDPC decoding algorithm in this article adopts the LLR-MP decoding algorithm. The computational complexity of the LDPC decoder is O log 2 ( | Q | ) N M , and the computational complexity of MP algorithm is O N M S Q , where S represents the number of non-zero elements present in each row or column of the DD domain channel matrix H . Thus, the overall computational complexity of the IDF-MP algorithm is N o u t ( N i n , M P O log 2 ( | Q | ) N M + N i n , L D P C O N M S Q ) . Algorithm 1 shows the process of IDF-MP detection algorithm based on LDPC.
Algorithm 1 Process of IDF-MP detection algorithm based on LDPC
Require:
  • Time DD domain channel matrix H , received signal y , maximum iteration times N o u t , N i n , M P , and N i n , L D P C ;
  • Set the outer-loop iteration times n o u t = 0 .
Ensure:
  • Bit information;
1:
When n o u t is less than N o u t , start MP detection until the number of internal loops in the detector reaches N i n , M P ;
2:
Calculate the symbol-level outer information M e , M P n o u t ( x i = a j ) of the detector output using Equation (26);
3:
Convert M e , M P n o u t ( x i = a j ) to bit-level information L e , M P n o u t ( c i , k = a j , k ) using Equation (27);
4:
After de-interleaving, the L p r , L D P C n o u t ( c i , k = a j , k ) is obtained as the prior information for the LDPC decode;
5:
Start LDPC decoding until the number of internal loops in the decoder reaches N i n , L D P C ;
6:
Calculate the bit-level LLR information L e , L D P C n o u t ( c i , k = a j , k ) output from the LDPC decoder using formula Equation (28);
7:
Obtain the bit-level input information L p r , M P n o u t ( c i , k = a j , k ) of the MP detector through an interleaver;
8:
Using Equation (32), convert the bit level information L p r , M P n o u t ( c i , k = a j , k ) into the probability values of each modulation symbol, as the prior probability information M p r , M P n o u t ( x i = a j ) input to the MP detector;
9:
n o u t = n o u t + 1 , return to step 2 if n o u t is less than N o u t , until n o u t reaches N o u t ;
10:
return Make a decision to output bit information based on L e , L D P C N o u t ( c i , k = a j , k ) ;

5. Performance Evaluation

In this section, the performance of our approach is evaluated. Table 1 shows the parameters used in the simulation experiments. The UWA channels were modeled following the methodology described in [33]. The water depth was set to 100 m, with the transmitter positioned at a depth of 30 m and the receiver at 50 m. The horizontal distance between the transmitter and receiver was 2 km. The sound speed profile utilized for channel simulation is illustrated in Figure 5.
Figure 6a illustrates the CIR without Doppler compensation with relative movement speed v = 10 kn. It is evident that the maximum channel delay spread exceeds 15 ms. The relative motion between the receiver and transmitter results in strong time variation in the UWA channel. Figure 6b depicts the Doppler dispersion in the DD domain distributed from 2 Hz to 3 Hz, indicating an overspread UWA channel. After Doppler compensation, Figure 6c shows that the CIR in the DD domain becomes sparser compared to the DT domain, confirming that the modeled channel becomes underspread.
Figure 7 compares the BER performance with different values of N and M under varying v. The frame lengths were identical, and the communication rates were set to the same value of 2.333 kbps. The relative speed between the transmitter and receiver was specified as v = 5 kn and v = 10 kn. In the v = 5 kn or v = 10 kn cases, OTFS frames with larger values of N exhibited higher performance due to their higher Doppler tolerance, provided M is large enough to avoid significant fractional delay. As the Doppler spread increases, the system with larger N shows substantial superiority.
Based on these analyses, setting M = 384 and N = 24 proved to be effective, allowing the system to better handle large Doppler spreads and obtain enhanced performance while maintaining a reasonable frame size.
Figure 8 compares the obtained performance based on different ZP lengths when setting the value of L g to 3 dB, 7 dB, or 10 dB.
The communication rate with maximum ZP length L m a x , maximum additional delay L 10 d B , maximum additional delay L 7 d B , and maximum additional delay L 3 d B were 2.3333 kbps, 2.5197 kbps, 2.8329 kbps, and 2.9579 kbps, respectively. We can see that setting L g to the maximum additional delay resulted in a higher communication rate than when using the maximum delay. As β in Equation (10) decreases, the communication rate improves. However, when β dB decreases, the ZP length may not sufficiently prevent inter-block interference, thereby diminishing its effectiveness and increasing the BER.
In this simulation, setting L g = L 10 d B improved the communication rate while ensuring communication reliability. Therefore, subsequent simulations used the parameters M = 384 , N = 24 , and L g = L 10 d B .
Figure 9 shows the BER performance of the IDF-MP algorithm with different numbers of outer-loop iterations. The outer-loop iteration number n o u t was set to 0, 1, 2, or 3, and the relative speed v was set to 5 kn or 10 kn. It can be observed that, regardless of whether v = 5 kn or 10 kn, the BER decreased as the number of outer-loop iterations increased. There was a minimal difference in BER performance between three iterations and two iterations. Therefore, subsequent simulations employed two outer-loop iterations.
To comprehensively analyze the IDF-MP algorithm proposed in this study, we compared its BER performance with OFDM, an initial estimate iterative decoding feedback (IE-IDF-MRC) algorithm [17], and a two-dimensional passive time reversal decision feedback equalization (2D-PTR-DFE) algorithm [21]. The IDF-MP algorithm is built upon the proposed dual-domain joint channel estimation technique, distinguishing it from other algorithms that rely on the OMP method for channel estimation. Both the IE-IDF-MRC and IDF-MP algorithms were configured with an outer-loop iteration count of 2, and LDPC codes were employed for channel coding across all algorithms, achieving a communication rate of 2.5 kbps.
As shown in Figure 10, the results demonstrated that the OTFS communication system outperformed the OFDM system, regardless of whether v = 5 kn or 10 kn, with a more pronounced performance difference observed under severe channel time variations. Specifically, the IDF-MP algorithm proposed in this study exhibited the best performance. This superiority can be attributed to the IDF-MP algorithm’s utilization of the soft information output from LDPC decoding as prior information for the MP detector, facilitating iterative information exchange that minimizes information loss and enhances the system’s detection and decoding performance.

6. Conclusions

In this study, we explored the application of OTFS modulation in UWA communications, particularly for overspread UWA channels. At the transmitter of the proposed UWA OTFS system, we introduced a frame structure design which is tailored with respect to the characteristics of UWA channels, in order to avoid inter-block interference and mitigate the impact of overspread channels. At the receiver, we introduced a dual-domain joint channel estimation method and an iterative decoding feedback-based detector to enhance the robustness of UWA OTFS communication systems.
Our simulation results demonstrated the feasibility of the proposed OTFS system in overspread UWA channels. Through meticulous frame design, the OTFS system exhibits significant performance improvements in the case of strong overspread. Moreover, leveraging the proposed channel estimation method and detector, our OTFS system achieved enhanced communication performance when compared to conventional systems such as the OFDM, IE-IDF-MRC, and 2D-PTR-DEF methods.
Our future work will involve deploying the proposed OTFS communication system in real-world UWA environments in order to further analyze its operational performance and effectiveness.

Author Contributions

Conceptualization, L.L.; methodology, L.L.; software, L.L.; validation, C.M.; formal analysis, L.L.; writing—original draft preparation, L.L.; writing—review and editing, C.M. and Y.D.; visualization, C.M. and X.Q.; supervision, Y.D. and X.L.; project administration, X.L. and X.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Stable Supporting Fund of Acoustic Science and Technology Laboratory (Grant No. JCKYS2024604SSJS012).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to thank the anonymous reviewers for their careful assessment of our work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhao, H.; Ji, F.; Wang, Y.; Yao, K.; Chen, F. Space–Air–Ground–Sea Integrated Network with Federated Learning. Remote Sens. 2024, 16, 1640. [Google Scholar] [CrossRef]
  2. Liu, H.; Ma, L.; Wang, Z.; Qiao, G. Channel Prediction for Underwater Acoustic Communication: A Review and Performance Evaluation of Algorithms. Remote Sens. 2024, 16, 1546. [Google Scholar] [CrossRef]
  3. Stojanovic, M.; Preisig, J. Underwater acoustic communication channels: Propagation models and statistical characterization. IEEE Commun. Mag. 2009, 47, 84–99. [Google Scholar] [CrossRef]
  4. Wang, C.; Yin, J.; Du, P.; Guo, L. Application oforthogonal frequency division multiplexing incognitive underwater communication. J. Acoust. Soc. Am. 2012, 132, 2015. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Wang, C.; Yin, J.; Sheng, X. Research on multilevel differential amplitude and phase-shift key in gin convolution-coded orthogonal frequency division multiplexing underwater communication system. J. Acoust. Soc. Am. 2012, 132, 2015. [Google Scholar] [CrossRef]
  6. Qiao, G.; Liu, L.; Ma, L. Adaptive downlink OFDMA system with low-overhead and limited feedback in time-varying underwater acoustic channel. IEEE Access 2019, 7, 12729–12741. [Google Scholar] [CrossRef]
  7. Liu, L.; Cai, L.; Ma, L.; Qiao, G. Channel state information prediction for adaptive underwater acoustic downlink OFDMA system: Deep neural networks based approach. IEEE Trans. Veh. Technol. 2021, 70, 9063–9076. [Google Scholar] [CrossRef]
  8. Liu, L.; Ma, C.; Duan, Y. A channel temporal correlation-based optimization method for imperfect underwater acoustic channel state information. Phys. Commun. 2023, 58, 102021. [Google Scholar] [CrossRef]
  9. Ma, C.; Liu, L.; Wang, Q.; Zhang, W.; Liu, X. High-performance deep-sea long-range underwater acoustic communication: Deconvolved conventional beamforming based approach. Phys. Commun. 2024, 64, 102339. [Google Scholar] [CrossRef]
  10. Wu, F.; Tian, T.; Su, B.; Song, Y. Hadamard–Viterbi Joint Soft Decoding for MFSK Underwater Acoustic Communications. Remote Sens. 2022, 14, 6038. [Google Scholar] [CrossRef]
  11. Fang, T.; Wang, Q.; Zhang, L.; Liu, S. Modulation Mode Recognition Method of Non-Cooperative Underwater Acoustic Communication Signal Based on Spectral Peak Feature Extraction and Random Forest. Remote Sens. 2022, 14, 1603. [Google Scholar] [CrossRef]
  12. Zhou, M.; Wang, J.; Feng, X.; Sun, H.; Qi, J.; Lin, R. Neural-Network-Based Equalization and Detection for Underwater Acoustic Orthogonal Frequency Division Multiplexing Communications: A Low-Complexity Approach. Remote Sens. 2023, 15, 3796. [Google Scholar] [CrossRef]
  13. Wang, Y.; Wang, G.; He, R.; Ai, B.; Tellambura, C. Doppler Shift and Channel Estimation for Intelligent Transparent Surface Assisted Communication Systems on High-Speed Railways. IEEE Trans. Commun. 2023, 71, 4204–4215. [Google Scholar] [CrossRef]
  14. Wang, T.; Proakis, J.G.; Masry, E.; Zeidler, J.R. Performance degradation of OFDM systems due to Doppler spreading. IEEE Trans. Wirel. Commun. 2006, 5, 1422–1432. [Google Scholar] [CrossRef]
  15. Raviteja, P.; Phan, K.T.; Hong, Y.; Viterbo, E. Interference cancellation and iterative detection for orthogonal time frequency space modulation. IEEE Trans. Wirel. Commun. 2018, 17, 6501–6515. [Google Scholar] [CrossRef]
  16. Hadani, R.; Rakib, S.; Tsatsanis, M.; Monk, A.; Goldsmith, A.J.; Molisch, A.F.; Calderbank, R. Orthogonal time frequency space modulation. In Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC), San Francisco, CA, USA, 19–22 March 2017. [Google Scholar]
  17. Thaj, T.; Viterbo, E. Low Complexity Iterative Rake Decision Feedback Equalizer for Zero-Padded OTFS systems. IEEE Trans. Veh. Technol. 2020, 69, 15606–15622. [Google Scholar] [CrossRef]
  18. Francis, J.; Chivurala, P.; Koilpillai, R.D. Performance of OTFS and OCDM schemes in underwater acoustic communication channels. In Proceedings of the OCEANS Conference, Chennai, India, 21–24 February 2022. [Google Scholar]
  19. Feng, X.; Esmaiel, H.; Wang, J.; Qi, J.; Zhou, M.; Qasem, Z.A.; Sun, H.; Gu, Y. Underwater acoustic communications based on OTFS. In Proceedings of the 15th IEEE International Conference on Signal Processing (ICSP), Beijing, China, 6–9 December 2020. [Google Scholar]
  20. Bocus, M.; Doufexi, A.; Agrafiotis, D. Performance of OFDM-based massive MIMO OTFS systems for underwater acoustic communication. IET Commun. 2020, 14, 588–593. [Google Scholar] [CrossRef]
  21. Jing, L.; Zhang, N.; He, C.; Shang, J.; Liu, X.; Yin, H. OTFS underwater acoustic communications based on passive time reversal. Appl. Acoust. 2022, 185, 108386. [Google Scholar] [CrossRef]
  22. Hang, S.; Li, W. OTFS for Underwater Acoustic Communications: Practical System Design and Channel Estimation. In Proceedings of the OCEANS 2022, Hampton Roads, VA, USA, 17–20 October 2022. [Google Scholar]
  23. Li, W.; Lin, B.; Guo, R.; Hao, Z. OTFS for Underwater Acoustic Communications: Frame Design and Channel Estimation. In Proceedings of the OCEANS 2023—MTS/IEEE U.S. Gulf Coast, Biloxi, MS, USA, 25–28 September 2023. [Google Scholar]
  24. Raviteja, P.; Hong, Y.; Viterbo, E.; Biglieri, E. Practical pulse-shaping waveforms for reduced-cyclic-prefix OTFS. IEEE Trans. Veh. Technol. 2018, 68, 957–961. [Google Scholar] [CrossRef]
  25. Raviteja, P.; Phan, K.; Hong, Y.; Viterbo, E. Embedded delay-Doppler channel estimation for orthogonal time frequency space modulation. In Proceedings of the 88th IEEE Vehicular Technology Conference (VTC-Fall), Chicago, IL, USA, 27–30 August 2018. [Google Scholar]
  26. Srivastava, S.; Singh, R.K.; Jagannatham, A.K.; Hanzo, L. Bayesian learning aided sparse channel estimation for orthogonal time frequency space modulated systems. IEEE Trans. Veh. Technol. 2021, 70, 8343–8348. [Google Scholar] [CrossRef]
  27. Wei, Z.; Yuan, W.; Lit, S.; Yuant, J.; Ngt, D.W. A new off-grid channel estimation method with sparse Bayesian learning for OTFS systems. In Proceedings of the IEEE Global Communications Conference (GLOBECOM), Madrid, Spain, 7–11 December 2021. [Google Scholar]
  28. Shen, W.; Dai, L.; An, J.; Fan, P.; Heath, R.W. Channel estimationfor orthogonal time frequency space (OTFS) massive MIMO. IEEE Trans. Signal Process. 2019, 67, 4204–4217. [Google Scholar] [CrossRef]
  29. Berger, R.; Zhou, S.; Preisig, J.; Willett, P. Sparse channel estimation for multicarrier underwater acoustic communication: From subspace methods to compressed sensing. IEEE Trans. Signal Process. 2010, 58, 1708–1721. [Google Scholar] [CrossRef]
  30. Ma, L.; Qiao, G.; Liu, S. A combined Doppler scale estimation scheme for underwater acoustic OFDM system. J. Comput. Acoust. 2015, 23, 1540004. [Google Scholar] [CrossRef]
  31. Li, W.; Zhou, S.; Willett, P.; Zhang, Q. Preamble Detection for Underwater Acoustic Communications Based on Sparse Channel Identification. IEEE J. Ocean. Eng. 2019, 44, 256–268. [Google Scholar] [CrossRef]
  32. Ma, L.; Zhou, S.; Qiao, G.; Liu, S.; Zhou, F. Superposition coding for downlink underwater acoustic OFDM. IEEE J. Ocean. Eng. 2017, 42, 175–187. [Google Scholar] [CrossRef]
  33. Qarabaqi, P.; Stojanovic, M. Statistical Characterization and Computationally Efficient Modeling of a Class of Underwater Acoustic Communication Channels. IEEE J. Ocean. Eng. 2013, 38, 701–717. [Google Scholar] [CrossRef]
Figure 1. Zero-padded OTFS system.
Figure 1. Zero-padded OTFS system.
Jmse 12 01559 g001
Figure 2. The UWA channel power response.
Figure 2. The UWA channel power response.
Jmse 12 01559 g002
Figure 3. The dual-domain joint UWA OTFS channel estimation method.
Figure 3. The dual-domain joint UWA OTFS channel estimation method.
Jmse 12 01559 g003
Figure 4. The structure of IDF-MP.
Figure 4. The structure of IDF-MP.
Jmse 12 01559 g004
Figure 5. Sound speed profile.
Figure 5. Sound speed profile.
Jmse 12 01559 g005
Figure 6. The same UWA Channel in (a) DT domain, (b) DD domain, and (c) DD domain after Doppler compensation.
Figure 6. The same UWA Channel in (a) DT domain, (b) DD domain, and (c) DD domain after Doppler compensation.
Jmse 12 01559 g006
Figure 7. BER comparison of OTFS frames with different M and N under different relative speeds.
Figure 7. BER comparison of OTFS frames with different M and N under different relative speeds.
Jmse 12 01559 g007
Figure 8. BER comparison of OTFS frames with different ZP lengths.
Figure 8. BER comparison of OTFS frames with different ZP lengths.
Jmse 12 01559 g008
Figure 9. BER comparison of the IDF-MP method under different iterations.
Figure 9. BER comparison of the IDF-MP method under different iterations.
Jmse 12 01559 g009
Figure 10. BER comparison of the IDF-MP, OFDM, IE-IDF-MRC, and 2D-PTR-DEF methods.
Figure 10. BER comparison of the IDF-MP, OFDM, IE-IDF-MRC, and 2D-PTR-DEF methods.
Jmse 12 01559 g010
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersValue
Modulation modeQPSK
Channel coding1/2 LDPC
Carrier frequency10 kHz
Bandwidth4 kHz
Max path delay τ max 35 ms
Relative speed of transmitter and receiver v5 kn/10 kn
Number of subcarriers M384/512
Number of OTFS symbols N24/16
Sampling frequency48 kHz
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, L.; Ma, C.; Duan, Y.; Liu, X.; Qing, X. Channel Estimation and Iterative Decoding for Underwater Acoustic OTFS Communication Systems. J. Mar. Sci. Eng. 2024, 12, 1559. https://doi.org/10.3390/jmse12091559

AMA Style

Liu L, Ma C, Duan Y, Liu X, Qing X. Channel Estimation and Iterative Decoding for Underwater Acoustic OTFS Communication Systems. Journal of Marine Science and Engineering. 2024; 12(9):1559. https://doi.org/10.3390/jmse12091559

Chicago/Turabian Style

Liu, Lei, Chao Ma, Yong Duan, Xinyu Liu, and Xin Qing. 2024. "Channel Estimation and Iterative Decoding for Underwater Acoustic OTFS Communication Systems" Journal of Marine Science and Engineering 12, no. 9: 1559. https://doi.org/10.3390/jmse12091559

APA Style

Liu, L., Ma, C., Duan, Y., Liu, X., & Qing, X. (2024). Channel Estimation and Iterative Decoding for Underwater Acoustic OTFS Communication Systems. Journal of Marine Science and Engineering, 12(9), 1559. https://doi.org/10.3390/jmse12091559

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop