# Performance Comparison of Doppler Scale Estimation Methods for Underwater Acoustic OFDM

Authors: Lei Wan; Zhaohui Wang; Shengli Zhou; T. C. Yang; Zhijie Shi
Journal: Journal of Electrical and Computer Engineering
 (2012)

---

## Abstract

Doppler scale estimation is one critical step needed by the resampling operation in acoustic communication
receivers. In this paper, we compare different Doppler scale estimation methods using
either cyclic-prefixed (CP) or zero-padded (ZP) orthogonal-frequency division-multiplexing (OFDM)
waveforms. For a CP-OFDM preamble, a self-correlation method allows for blind Doppler scale estimation
based on an embedded repetition structure while a cross-correlation method is available with
the knowledge of the waveform. For each received ZP-OFDM block, the existence of null subcarriers
allows for blind Doppler scale estimation. In addition, a pilot-aided method and a decision-aided method
are applicable based on cross-correlation with templates constructed from symbols on pilot subcarriers
only and from symbols on all subcarriers after data decoding, respectively. This paper carries out
extensive comparisons among these methods using both simulated and real experimental data. Further,
the applicabilities of these methods to distributed multiuser systems are investigated.


---

## Body

## 1.  Introduction

Underwater acoustic communications and networking have been under extensive investigation in recent years [1, 2]. Considerable progress on the physical layer communication techniques has been made for both single-carrier and multicarrier communications; see, for example, [3–19]. Relative to the radio channel which has relative short delay spread and slow time variation, underwater acoustic channels typically exhibit long delay spread and fast time variation. The latter brings significant Doppler effects to underwater acoustic communication systems, hence estimation of the Doppler scaling factor is one key receiver module [4, 20, 21].Typically, Doppler scale estimation is accomplished by inserting waveforms known to the receiver during the data transmission. Two popular approaches are described in the following.(i)
One approach is to use a pulse train which is formed by the repetition of aDoppler-insensitive waveform [22], such as linear-frequency-modulated (LFM) waveform [23] and hyperbolic-frequency-modulated (HFM) waveform [24]. A transmission format with one preamble and one postamble around the data burst is usually adopted [4, 20, 25], as shown in Figure 1. At the receiver side, by detecting the times-of-arrival of the preamble and postamble, thus the interval change in-between, an average Doppler scale estimate over the whole data burst can be obtained. Thanks to the Doppler-insensitive property of the waveforms, a single-branch-matched filtering operation is adequate even in the presence of Doppler distortion. However, this method is only suitable for offline processing due to the processing delay.(ii)
The other approach is to use aDoppler-sensitive waveform with a thumb-tack ambiguity function. A Doppler-sensitive waveform is usually transmitted as a preamble prior to the data burst, as shown in Figure 1. At the receiver side, a bank of correlators correlates the received signal with preambles prescaled by different Doppler scaling factors, and the branch with the largest correlation peak provides the estimated Doppler scale [25]. Typical Doppler-sensitive waveforms include Costa waveforms [26], m-sequence [27], and poly-phase sequence [28].Figure 1
Top: Doppler scale estimation with Doppler-insensitive waveforms. Bottom: Doppler scale estimation with Doppler-sensitive waveforms.In this paper, we focus on an underwater acoustic communication system using zero-padded orthogonal-frequency division-multiplexing modulation (ZP-OFDM), in which pilot subcarriers and null subcarriers are usually multiplexed with data subcarriers for channel estimation and residual Doppler shift mitigation, respectively [4]. A cyclic prefixed (CP) OFDM preamble is inserted prior to data transmission for detection, synchronization, and Doppler scale estimation [29]. This transmission format, as shown in Figure 2, has been implemented on DSP-based OFDM modem prototypes [30].Figure 2
The data burst structure considered in this paper, which consists of a special CP-OFDM preamble and multiple ZP-OFDM blocks.By exploiting the cyclic repetition structure of the CP-OFDM preamble, a blind estimation with a bank of self-correlators was proposed in [29]. However, it does not leverage the knowledge of the waveform itself which is known to the receiver. Taking this method as the first approach, one can easily construct the following Doppler scale estimators for the OFDM transmission in Figure 2.(i)
Cross-correlation with the CP-OFDM Preamble: Given the Doppler sensitivity of the OFDM waveform, a bank of cross-correlators can use the prescaled versions of the CP-OFDM waveform as local replicas.(ii)
Pilot-Aided Method for Each ZP-OFDM Block: By taking the waveform constituted by the pilot-subcarrier components as a replica of the transmitted signal, the Doppler estimation method using a bank of cross-correlators is directly applicable.(iii)
Null-Subcarrier Based Blind Estimation Method for Each ZP-OFDM Block: As an extension of the blind carrier frequency offset (CFO) estimation method [31], the receiver rescales the received waveform with different tentative Doppler scaling factors and uses the energy on the null subcarriers to find the best fit.(iv)
Decision-Aided Method for Each ZP-OFDM Block: Once a ZP-OFDM block is successfully decoded, the transmitted waveform corresponding to this block can be reconstructed at the receiver. Taking the reconstructed waveform as a local replica, the Doppler estimation method using a bank of correlators can be deployed to refine the Doppler scale estimation for this block. The refined Doppler scale estimate can be passed to the next block.The contributions of this paper are the following.(i)
We carry out extensive performance comparisons among the aforementioned Doppler estimation methods. Specifically, we focus on the OFDM transmission format in Figure2 in single-user transmissions. Both simulations and experimental results reveal that the correlation-based methods have a decent performance in the low SNR region, and the blind estimation methods can catch up or even outperform the correlation methods in the high SNR region. As a performance benchmark, the Cramer-Rao lower bound (CRLB) is also included for single-path channels.(ii)
We extend our investigation to a multiuser OFDM setting, where different users could have different Doppler scaling factors [32]. Simulation results show that the correlation-based methods are robust to the multiuser interference, while the blind method suffers severe performance degradation.The rest of this paper is as follows. Different Doppler scale estimation methods for CP-OFDM and ZP-OFDM waveforms are presented in Sections2 and 3, respectively. Simulation results of these methods are provided in Section 4, and experimental results are provided in Section 5. Extension to the multiuser scenario is described in Section 6. Conclusions are contained in Section 7.


## 2.  Doppler Scale Estimation with a CP-OFDM Preamble

Consider a CP-OFDM preamble structure in Figure2, which consists of two identical OFDM symbols of length T0 and a cyclic prefix of length Tcp in front, with the embedded structure(1)xcp(t)=xcp(t+T0),-Tcp≤t≤T0.LetB denote the system bandwidth, and define K0:=BT0 as the number of subcarriers. The baseband CP-OFDM signal is(2)xcp(t)=∑k=-K0/2K0/2-1d[k]ej2π(k/2T0)q(t),t∈[-Tcp,2T0],
where d[k] is the transmitted symbol on the kth subcarrier, and q(t) is a pulse shaping window,(3)q(t)={1,t∈[-Tcp,2T0].0,elsewhere.
The passband signal can be obtained as x̃cp(t)=2Re{xcp(t)ej2πfct}, where fc is the center frequency.Consider a multipath channel which consists ofNpa paths(4)h(t;τ)=∑p=1NpaAp(t)δ(t-τp(t)),
where Ap(t) and τp(t) denote the amplitude and delay of the pth path, respectively. Throughout this paper, we assume that the amplitude is constant within each OFDM block (about 200 ms for the system considered in this paper), that is, Ap(t)≈Ap, which leads to(5)h(t;τ)=∑p=1NpaApδ(t-τp(t)).After transmitting the passband signalx̃cp(t) through the multipath channel, the received passband signal ỹ(t) is converted to baseband as y(t)=LPF(ỹ(t)e-j2πfct), where LPF denotes the low pass filtering operation.


### 2.1.  Self-Correlation

If all the paths in the channel have the same Doppler scale factor(6)τp(t)=τp-at,
it is shown in [29] that the embedded structure in the received waveform becomes(7)y(t)=e-j2π(a/(1+a))fcT0y(t+T01+a),-Tcp-τmax⁡1+a≤t≤T01+a,
which has a repetition period T0/(1+a) regardless of the channel amplitudes.By exploiting the structure in (7), the time-of-arrival and the Doppler scale of the CP-OFDM symbol in the received signal can be jointly estimated via(8)(â,τ̂)=argmax⁡a,τ|∫0T0/(1+a)y(t+τ)y*(t+τ+T01+a)dt|,
which does not require the knowledge of the channel and the data symbols. This method can be implemented with a bank of self-correlators [29].


### 2.2.  Cross-Correlation

Rather than exploiting the structure of the CP-OFDM preamble, the cross-correlation-based method can be used, since the transmitted preamble is known at the receiver. Taking the basic unit of durationT0 as the template, the joint time-of-arrival and Doppler rate estimation can be achieved via(9)(â,τ̂)=argmax⁡a,τ|∫0T0y(t+τ)xcp*((1+a)t)e-j2πafctdt|.
This can be implemented via a bank of cross-correlators, where the branch yielding the largest peak provides the needed Doppler scale estimate.


## 2.1.  Self-Correlation

If all the paths in the channel have the same Doppler scale factor(6)τp(t)=τp-at,
it is shown in [29] that the embedded structure in the received waveform becomes(7)y(t)=e-j2π(a/(1+a))fcT0y(t+T01+a),-Tcp-τmax⁡1+a≤t≤T01+a,
which has a repetition period T0/(1+a) regardless of the channel amplitudes.By exploiting the structure in (7), the time-of-arrival and the Doppler scale of the CP-OFDM symbol in the received signal can be jointly estimated via(8)(â,τ̂)=argmax⁡a,τ|∫0T0/(1+a)y(t+τ)y*(t+τ+T01+a)dt|,
which does not require the knowledge of the channel and the data symbols. This method can be implemented with a bank of self-correlators [29].


## 2.2.  Cross-Correlation

Rather than exploiting the structure of the CP-OFDM preamble, the cross-correlation-based method can be used, since the transmitted preamble is known at the receiver. Taking the basic unit of durationT0 as the template, the joint time-of-arrival and Doppler rate estimation can be achieved via(9)(â,τ̂)=argmax⁡a,τ|∫0T0y(t+τ)xcp*((1+a)t)e-j2πafctdt|.
This can be implemented via a bank of cross-correlators, where the branch yielding the largest peak provides the needed Doppler scale estimate.


## 3.  Doppler Scale Estimation with Each ZP-OFDM Block

As described in [4], a ZP-OFDM signal design multiplexing pilot and null subcarriers with data subcarriers can effectively deal with fast channel variations. Assume that the ZP-OFDM system has K subcarriers. Let T denote the symbol duration and Tg the guard interval. The total OFDM block duration is thus Tbl:=T+Tg. Denote 𝒮D, 𝒮P, 𝒮N as the nonoverlapped sets formed by the data subcarriers, pilot subcarriers, and null subcarriers, respectively, which satisfy 𝒮D⋃𝒮P⋃𝒮N={-K/2,…,K/2-1}. The baseband transmitted ZP-OFDM signal can be expressed by(10)xzp(t)=∑k∈SD⋃SPd[k]ej2π(k/T)tg(t),t∈[0,Tbl],
where g(t) describes the zero-padding operation, that is,(11)g(t)={1,t∈[0,T],0,elsewhere.After transmitting the ZP-OFDM symbol through a multipath channel defined in (5), we denote ỹ(t) as the received passband signal, whose baseband version is y(t)=LPF(ỹ(t)e-j2πfct). The availability of null subcarriers, pilot subcarriers, and data subcarriers can be used for Doppler scale estimation.


### 3.1.  Null-Subcarrier-Based Blind Estimation

In [29], the null subcarriers in ZP-OFDM system are exploited to perform carrier frequency offset (CFO) estimation. Here in this paper, the same principle is used to estimate Doppler scale factor.Assume that coarse synchronization is available from the preamble. After truncating each ZP-OFDM block from the received signal, we resample one block with different tentative scaling factors. The total energy of frequency measurements at null subcarriers are used as a metric for the Doppler scale estimation(12)â=argmin⁡a∑k∈SN|∫0T+Tgy(t1+a)e-j2πafcte-j2π(k/T)tdt|2.
For each tentative a, a resampling operation is carried out followed by fast Fourier transform. A one-dimensional grid search leads to a Doppler scale estimate.


### 3.2.  Pilot-Aided Estimation

As introduced above, a set of subcarriers𝒮P is dedicated to transmit pilot symbols. Hence, the transmitted waveform xzp(t) is partially known, containing(13)xpilot(t)=∑k∈SPd[k]ej2π(k/T)tg(t),t∈[0,T].The joint time-of-arrival and Doppler scale estimation is achieved via(14)(â,τ̂)=argmax⁡a,τ|∫0T/(1+a)y(t+τ)xpilot*((1+a)t-τ)e-j2πafctdt|,
which can be implemented via a bank of cross-correlators.


### 3.3.  Decision-Aided Estimation

For an OFDM transmission with multiple blocks, the Doppler estimated in one block can be used for the resampling operation of the next block assuming small Doppler variation across blocks. After the decoding operation the receiver can reconstruct the transmitted time-domain waveform, by replacingd[k] by its estimate d̂[k], for all k∈𝒮D in. Denote the reconstructed waveform as x̂zp(t).Similar to the pilot-aided method, the decision-aided method performs the joint time-of-arrival and Doppler scale estimation via(15)(â,τ̂)=argmax⁡a,τ|∫0T/(1+a)y(t+τ)x̂zp*((1+a)t-τ)e-j2πafctdt|,
which again is implemented via a bank of cross-correlators. The estimated â can be used for the resampling operation of the next block.Remark 1.
Relative to the pilot-aided method, the decision-aided method leverages the estimated information symbols, thus is expected to achieve a better estimation performance. Assuming that all the information symbols have been successfully decoded, the decision-aided method has knowledge about both the data and pilot symbols. Let|𝒮P| and |𝒮D| denote the numbers of pilot and data symbols, respectively. Using the template x̂zp(t) constructed from (|𝒮P|+|𝒮D|) known symbols for cross correlation achieve a 10log10((|𝒮P|+|𝒮D|)/|𝒮P|) dB power gain in terms of noise reduction, relative to that using the template xpilot(t) constructed from |𝒮P| known symbols.


## 3.1.  Null-Subcarrier-Based Blind Estimation

In [29], the null subcarriers in ZP-OFDM system are exploited to perform carrier frequency offset (CFO) estimation. Here in this paper, the same principle is used to estimate Doppler scale factor.Assume that coarse synchronization is available from the preamble. After truncating each ZP-OFDM block from the received signal, we resample one block with different tentative scaling factors. The total energy of frequency measurements at null subcarriers are used as a metric for the Doppler scale estimation(12)â=argmin⁡a∑k∈SN|∫0T+Tgy(t1+a)e-j2πafcte-j2π(k/T)tdt|2.
For each tentative a, a resampling operation is carried out followed by fast Fourier transform. A one-dimensional grid search leads to a Doppler scale estimate.


## 3.2.  Pilot-Aided Estimation

As introduced above, a set of subcarriers𝒮P is dedicated to transmit pilot symbols. Hence, the transmitted waveform xzp(t) is partially known, containing(13)xpilot(t)=∑k∈SPd[k]ej2π(k/T)tg(t),t∈[0,T].The joint time-of-arrival and Doppler scale estimation is achieved via(14)(â,τ̂)=argmax⁡a,τ|∫0T/(1+a)y(t+τ)xpilot*((1+a)t-τ)e-j2πafctdt|,
which can be implemented via a bank of cross-correlators.


## 3.3.  Decision-Aided Estimation

For an OFDM transmission with multiple blocks, the Doppler estimated in one block can be used for the resampling operation of the next block assuming small Doppler variation across blocks. After the decoding operation the receiver can reconstruct the transmitted time-domain waveform, by replacingd[k] by its estimate d̂[k], for all k∈𝒮D in. Denote the reconstructed waveform as x̂zp(t).Similar to the pilot-aided method, the decision-aided method performs the joint time-of-arrival and Doppler scale estimation via(15)(â,τ̂)=argmax⁡a,τ|∫0T/(1+a)y(t+τ)x̂zp*((1+a)t-τ)e-j2πafctdt|,
which again is implemented via a bank of cross-correlators. The estimated â can be used for the resampling operation of the next block.Remark 1.
Relative to the pilot-aided method, the decision-aided method leverages the estimated information symbols, thus is expected to achieve a better estimation performance. Assuming that all the information symbols have been successfully decoded, the decision-aided method has knowledge about both the data and pilot symbols. Let|𝒮P| and |𝒮D| denote the numbers of pilot and data symbols, respectively. Using the template x̂zp(t) constructed from (|𝒮P|+|𝒮D|) known symbols for cross correlation achieve a 10log10((|𝒮P|+|𝒮D|)/|𝒮P|) dB power gain in terms of noise reduction, relative to that using the template xpilot(t) constructed from |𝒮P| known symbols.


## 4.  Simulation Results

The OFDM parameters are summarized in Table1. For CP-OFDM, the data symbols at all the 512 subcarriers are randomly drawn from a QPSK constellation. For ZP-OFDM, out of 1024 subcarriers, there are |𝒮N|=96 null subcarriers with 24 on each edge of the signal band for band protection and 48 evenly distributed in the middle for the carrier frequency offset estimation; |𝒮P|=256 are pilot subcarriers uniformly distributed among the 1024 subcarriers, and the remaining are |𝒮D|=672 data subcarriers for delivering information symbols. The pilot symbols are drawn randomly from a QPSK constellation. The data symbols are encoded with a rate-1/2 nonbinary LDPC code [33] and modulated by a QPSK constellation.Table 1
OFDM parameters in simulations.
System parametersCP-OFDMZP-OFDMCenter frequency:fc13 kHz13 kHzBandwidth:B4.88 kHz4.88 kHz# of subcarriers:K0=512K=1024Time duration:T0=104.86msT=209.72msGuard interval:Tcp=100msTg=40.3msThree UWA channel settings are tested.(i)
Channel Setting 1: A single-path channel:
(16)h(t,τ)=δ(t-[τ-at]).(ii)
Channel Setting 2: A multipath channel with Npa=15 paths, where all paths have one common Doppler scaling factor:
(17)h(t,τ)=∑p=1NpaApδ(t-[τp-at]).(iii)
Channel Setting 3:  A multipath channel with Npa=15 paths, where each path has an individual Doppler scaling factor:
(18)h(t,τ)=∑p=1NpaApδ(t-[τp-apt]).The interarrival time of paths follows an exponential distribution with a mean of 1 ms. The mean delay spread for the channels in and is thus 15 ms. The amplitudes of paths are Rayleigh distributed with the average power decreasing exponentially with the delay, where the difference between the beginning and the end of the guard time is20dB. For each path, the Doppler scale ap is generated from a Doppler speed vp (with unit m/s):(19)ap=vpc,
where c=1500m/s is the sound speed in water. In channel settings 1 and 2, the Doppler speed v is uniformly distributed within [-4.5,4.5] m/s. In channel setting 3, the Doppler speeds {vp} are randomly drawn from the interval [1.5-0.1,1.5+0.1] m/s.In channel settings 1 and 2, the ground truths ofv and a are known. We adopt the root-mean-squared-error (RMSE) of the estimated Doppler speed as the performance metric,(20)RMSE=E[|v̂-v|2]=E[|(â-a)c|2],
which has the unit m/s. In channel setting 3, different paths have different Doppler scales, while the Doppler scale estimator only provides one estimate. RMSE is hence not well motivated. With the estimated Doppler scale to perform the resampling operation, we will use the block-error-rate (BLER) of the ZP-OFDM decoding as the performance metric.


### 4.1.  RMSE Performance with CP-OFDM

For the single-path channel, Figure3 shows the RMSE performance of two estimation methods at different SNR levels. One can see a considerable gap between the self-correlation method and the cross-correlation method, while in the medium to high SNR region, both methods can provide a reasonable performance to facilitate receiver decoding.Figure 3
Performance of different estimators for the CP-OFDM preamble in single-path and multipath channels (channel settings 1 and 2).For the multipath channel with a single Doppler speed, Figure3 shows the RMSE performance of two estimation methods. One can see that the cross-correlation method outperforms the self-correlation method considerably in the low SNR region. However, the former suffers an error floor in the high SNR region, while the later does not.Relative to the RMSE performance in the single-path channel, a considerable performance degradation can be observed for the cross-correlation method in the multipath channel, whereas the performance of the self-correlation method is quite robust. The reason for the difference lies in the capability of the self-correlation method to collect the energy from all paths for Doppler scale estimation, while the cross-correlation method aims to get the Doppler scale estimate from only one path, the strongest path.


### 4.2.  RMSE Performance with ZP-OFDM

Figure4 shows the RMSE performance of three estimation methods for ZP-OFDM in single-path channels. In the low SNR region, one can see that the decision-aided method is the best, while the null-subcarrier-based blind method is the worst. As discussed in Remark 1, the decision-aided method achieves 10log10((|𝒮D|+|𝒮P|)/|𝒮P|)≈6dB power gain relative to the pilot-aided method. In the medium and high SNR region, the pilot-aided method suffers an error floor due to the interference from the data subcarriers, and the null-subcarrier-based blind method gets a good estimation performance. The Cramer-Rao lower bound (CRLB) with a known waveform is also included as the performance benchmark, whose derivation can be carried out similar to [34, 35].Figure 4
Performance of different estimators for ZP-OFDM in single-path channels (channel setting 1). The CRLB with all data known is included as a benchmark.Figure5 shows the RMSE performance of three methods in multipath channels with a common Doppler speed. For each realization, the Doppler scale, the path amplitudes, and delays are randomly generated. The RMSE corresponding to each method is calculated by averaging the estimation error over multiple realizations. Again, one can see that in the low SNR region, the decision-aided method has the best performance, while the null-subcarrier-based blind method is the worst. Different from the performance in the single-path channel, the decision-aided method has an error floor in the high SNR region, since it only picks up the maximum correlation peak of one path. On the other hand, the null-subcarrier method has robust performance in the presence of multiple paths.Figure 5
Performance of different estimators for ZP-OFDM in multipath channels with a common Doppler scale (channel setting 2).


### 4.3.  Comparison of Blind Methods of CP- and ZP-OFDM

The self-correlation method for the CP-OFDM preamble is closely related to the null-subcarrier-based blind method for ZP-OFDM. This can be easily verified by rewriting (2) as(21)xcp(t)=∑k=-K0K0-1d′[k]ej2π(k/(2T0))tq(t),t∈[-Tcp,2T0],
where d′[k]=0 when k is odd and d′[k]=d[k/2] when k is even. The cyclic repetition pattern in is generated by placing zeros on all odd subcarriers in a long OFDM symbol of duration 2T0. Hence, the self-correlation implementation could be replaced by the null-subcarrier-based implementation for the CP-OFDM preamble.Figure6 shows the performance comparison between the blind method for ZP-OFDM and that for CP-OFDM in the multipath channel with one Doppler scale factor, respectively. At low SNR, typically when it’s lower than 0 dB, the null-subcarrier-based method in CP-OFDM system has a better performance than that in the ZP-OFDM system, which is due to the fact that CP-OFDM system has 512 null subcarriers, more than 96 null subcarriers in the ZP-OFDM block. At high SNR, the null-subcarrier-based method in ZP-OFDM has better performance. The possible reason is that null subcarriers in ZP-OFDM are distributed with an irregular pattern, which could outperform the regular pattern in the CP-OFDM preamble.Figure 6
Null-subcarrier-based method in ZP-OFDM and CP-OFDM.


### 4.4.  BLER Performance with ZP-OFDM

With channels generated according to the channel setting 3, Figure7 shows the simulated BLER performance, where the received OFDM blocks are resampled with the Doppler scale estimates from different estimators and processed using the receiver from [4] and the LDPC decoder from [33]. At each SNR point, at least 20 block errors are collected.Figure 7
The BLER performance in multipath multi-Doppler channels (channel setting 3).It is expected that the OFDM system can only work when the useful signal power is above that of the ambient noise. Regarding the simulation results in Figure5, one can see that all the methods can reach a RMSE lower than 0.1 m/s. Hence, it is not surprising that these methods lead to quite similar BLER results as shown in Figure 7. This observation is consistent with the analysis in [29] that an estimation error of 0.1 m/s introduces a negligible error.


## 4.1.  RMSE Performance with CP-OFDM

For the single-path channel, Figure3 shows the RMSE performance of two estimation methods at different SNR levels. One can see a considerable gap between the self-correlation method and the cross-correlation method, while in the medium to high SNR region, both methods can provide a reasonable performance to facilitate receiver decoding.Figure 3
Performance of different estimators for the CP-OFDM preamble in single-path and multipath channels (channel settings 1 and 2).For the multipath channel with a single Doppler speed, Figure3 shows the RMSE performance of two estimation methods. One can see that the cross-correlation method outperforms the self-correlation method considerably in the low SNR region. However, the former suffers an error floor in the high SNR region, while the later does not.Relative to the RMSE performance in the single-path channel, a considerable performance degradation can be observed for the cross-correlation method in the multipath channel, whereas the performance of the self-correlation method is quite robust. The reason for the difference lies in the capability of the self-correlation method to collect the energy from all paths for Doppler scale estimation, while the cross-correlation method aims to get the Doppler scale estimate from only one path, the strongest path.


## 4.2.  RMSE Performance with ZP-OFDM

Figure4 shows the RMSE performance of three estimation methods for ZP-OFDM in single-path channels. In the low SNR region, one can see that the decision-aided method is the best, while the null-subcarrier-based blind method is the worst. As discussed in Remark 1, the decision-aided method achieves 10log10((|𝒮D|+|𝒮P|)/|𝒮P|)≈6dB power gain relative to the pilot-aided method. In the medium and high SNR region, the pilot-aided method suffers an error floor due to the interference from the data subcarriers, and the null-subcarrier-based blind method gets a good estimation performance. The Cramer-Rao lower bound (CRLB) with a known waveform is also included as the performance benchmark, whose derivation can be carried out similar to [34, 35].Figure 4
Performance of different estimators for ZP-OFDM in single-path channels (channel setting 1). The CRLB with all data known is included as a benchmark.Figure5 shows the RMSE performance of three methods in multipath channels with a common Doppler speed. For each realization, the Doppler scale, the path amplitudes, and delays are randomly generated. The RMSE corresponding to each method is calculated by averaging the estimation error over multiple realizations. Again, one can see that in the low SNR region, the decision-aided method has the best performance, while the null-subcarrier-based blind method is the worst. Different from the performance in the single-path channel, the decision-aided method has an error floor in the high SNR region, since it only picks up the maximum correlation peak of one path. On the other hand, the null-subcarrier method has robust performance in the presence of multiple paths.Figure 5
Performance of different estimators for ZP-OFDM in multipath channels with a common Doppler scale (channel setting 2).


## 4.3.  Comparison of Blind Methods of CP- and ZP-OFDM

The self-correlation method for the CP-OFDM preamble is closely related to the null-subcarrier-based blind method for ZP-OFDM. This can be easily verified by rewriting (2) as(21)xcp(t)=∑k=-K0K0-1d′[k]ej2π(k/(2T0))tq(t),t∈[-Tcp,2T0],
where d′[k]=0 when k is odd and d′[k]=d[k/2] when k is even. The cyclic repetition pattern in is generated by placing zeros on all odd subcarriers in a long OFDM symbol of duration 2T0. Hence, the self-correlation implementation could be replaced by the null-subcarrier-based implementation for the CP-OFDM preamble.Figure6 shows the performance comparison between the blind method for ZP-OFDM and that for CP-OFDM in the multipath channel with one Doppler scale factor, respectively. At low SNR, typically when it’s lower than 0 dB, the null-subcarrier-based method in CP-OFDM system has a better performance than that in the ZP-OFDM system, which is due to the fact that CP-OFDM system has 512 null subcarriers, more than 96 null subcarriers in the ZP-OFDM block. At high SNR, the null-subcarrier-based method in ZP-OFDM has better performance. The possible reason is that null subcarriers in ZP-OFDM are distributed with an irregular pattern, which could outperform the regular pattern in the CP-OFDM preamble.Figure 6
Null-subcarrier-based method in ZP-OFDM and CP-OFDM.


## 4.4.  BLER Performance with ZP-OFDM

With channels generated according to the channel setting 3, Figure7 shows the simulated BLER performance, where the received OFDM blocks are resampled with the Doppler scale estimates from different estimators and processed using the receiver from [4] and the LDPC decoder from [33]. At each SNR point, at least 20 block errors are collected.Figure 7
The BLER performance in multipath multi-Doppler channels (channel setting 3).It is expected that the OFDM system can only work when the useful signal power is above that of the ambient noise. Regarding the simulation results in Figure5, one can see that all the methods can reach a RMSE lower than 0.1 m/s. Hence, it is not surprising that these methods lead to quite similar BLER results as shown in Figure 7. This observation is consistent with the analysis in [29] that an estimation error of 0.1 m/s introduces a negligible error.


## 5.  Experimental Results

Thismobile acoustic communication experiment (MACE10) was carried out off the coast of Martha’s Vineyard, Massachusetts, June, 2010. The water depth was about 80 meters. The receiving array was stationary, while the source was towed slowly away from the receiver and then towed back, at a speed around 1 m/s. The relative distance of the transmitter and the receiver changed from 500 m to 4.5 km. Out of the two tows in this experiment, we only consider the data collected in the first tow. There are 31 transmissions in total, with a CP-OFDM preamble and 20 ZP-OFDM blocks in each transmission. We exclude one transmission file recorded during the turn of the source, where the SNR of the received signal is quite low.The CP-OFDM and ZP-OFDM parameters and signal structures are identical to that in the simulation, as listed in Table1.Figure8 shows the estimated Doppler speeds for ZP-OFDM blocks from different methods. Clearly, the Doppler speed fluctuates from block to block. Figure 9 shows the estimated channel impulse responses for two ZP-OFDM blocks from two data sets, where the time interval between these two data bursts is more than 1 hour. The channels have a delay spread about 20 ms but with different delay profiles.Figure 8
MACE10: Estimated Doppler speeds for 30 data bursts in MACE10, where each data  burst has 20 OFDM blocks. The time interval between two consecutive date bursts is around 4 mins.Estimated channel impulse responses for two different blocks at different bursts.
(a)
File ID: 1750155F1978_C0_S5(b)
File ID: 1750155F2070_C0_S5Based on the recorded files, we carried out two tests.(A) Test Case 1
In this test, we focus on one single file (file ID:1750155F1954_C0_S5) and compare the RMSE performance of different approaches by adding artificial noise to the recorded signal. The ground truth of the Doppler scale factor is not available. When computing the RMSE using (20) for each method, we use the estimated Doppler scale of the original file without adding the noise as the ground truth. Figure 10 shows the estimation performance of several approaches. Similar observations as the simulation results in Figures 3 and 5 are found.Figure 10
MACE10: Performance comparison of Doppler estimation approaches, file ID: 1750155F1954_C0_S5.(B) Test Case 2
In this test, we compare the BLER performance of an OFDM receiver where the resampling operation is carried out with different Doppler scale estimates from different methods.
Due to the relatively high SNR of the recorded signal, we create a semiexperimental data set by adding white Gaussian noise to the received signal. Defineσ̂2 as the estimated ambient noise power in the original recorded signal. Figure 11 shows the BLER performance with different Doppler estimation approaches by adding different amount of noises to the received files.
One can see that the methods for ZP-OFDM outperform the methods for CP-OFDM, as the Doppler scale itself is continuously changing from block to block, as illustrated in Figure8. Another interesting observation is that the null-subcarrier-based blind method has slight performance improvement relative to the pilot- and decision-aided methods. This agrees with the simulation results in Figure 5 that in the high SNR region, the blind estimation method does not suffer an error floor in the multipath channel, hence enjoys a better estimation performance.MACE10: BLER Performance using different Doppler estimation methods by adding artificial noise to the received signal,σ̂2  denoting the estimated ambient noise power.
(a)
Adding noise with powerσ̂2(b)
Added noise with power 2σ̂2


## 6.  Extension to Distributed MIMO-OFDM

If the transmitters in a multi-input multi-output (MIMO) system are co-located, the Doppler scales corresponding to all transmitters are similar, and hence a single-user blind Doppler scale estimation method would work well, as done in [10]. However, if the transmitters are distributed, for example in a system with multiple single-transmitter users, the Doppler scales for different users could be quite different, even with opposite signs [32]. We now investigate the performance of different Doppler scale estimation methods in the presence of multiuser interference. We will use the ZP-OFDM waveform as the reference design; similar conclusions can be applied to the CP-OFDM preamble. Only simulated data sets are used in the following tests.


### 6.1.  Pilot- and Decision-Aided Estimation

We simulate a two-user system. Each user generates a multipath channel according to channel setting 2 independently. The positions of pilot, null, and data subcarriers are the same for different users. The pilot and data symbols of different users are randomly generated and hence are different.Figure12 depicts the RMSE performance of the pilot- and decision-aided estimation methods. Compared with the performance in the single-user setting in Figure 5, there is performance degradation and the error floors are higher. However, both methods can achieve RMSE lower than 0.1 m/s at low SNR values. Hence, both methods have robust performance in the presence of multiuser interference.Figure 12
Pilot- and decision-aided Doppler scale estimation in a distributed two-user ZP-OFDM system.


### 6.2.  Null-Subcarrier-Based Blind Estimation

The null-subcarrier-based blind estimation method exploits the transmitted OFDM signal structure. Since all the users share the same positions of null subcarriers, there is a user-association problem even when multiple local minimums are found. We simulate a two-user system where the Doppler speeds of user 1 and user 2 are uniformly distributed within[-4.5,-0.5] m/s and [0.5,4.5] m/s, respectively. Without adding the ambient noise to the received signal, Figure 13 demonstrates both successful and failed cases using the objective function in (12). The objective functions in the single-user settings are also included for comparison. One can see that the multiuser interference degrades the estimation performance significantly. Hence, although the blind method developed for the single user case can be used to colocated MIMO-OFDM as in [10], it is not applicable to distributed MIMO-OFDM where different users have different Doppler scales.Illustration of the objective functions of the null-subcarrier-based method in a two-user system
(a)
Successful case(b)
Failed case


## 6.1.  Pilot- and Decision-Aided Estimation

We simulate a two-user system. Each user generates a multipath channel according to channel setting 2 independently. The positions of pilot, null, and data subcarriers are the same for different users. The pilot and data symbols of different users are randomly generated and hence are different.Figure12 depicts the RMSE performance of the pilot- and decision-aided estimation methods. Compared with the performance in the single-user setting in Figure 5, there is performance degradation and the error floors are higher. However, both methods can achieve RMSE lower than 0.1 m/s at low SNR values. Hence, both methods have robust performance in the presence of multiuser interference.Figure 12
Pilot- and decision-aided Doppler scale estimation in a distributed two-user ZP-OFDM system.


## 6.2.  Null-Subcarrier-Based Blind Estimation

The null-subcarrier-based blind estimation method exploits the transmitted OFDM signal structure. Since all the users share the same positions of null subcarriers, there is a user-association problem even when multiple local minimums are found. We simulate a two-user system where the Doppler speeds of user 1 and user 2 are uniformly distributed within[-4.5,-0.5] m/s and [0.5,4.5] m/s, respectively. Without adding the ambient noise to the received signal, Figure 13 demonstrates both successful and failed cases using the objective function in (12). The objective functions in the single-user settings are also included for comparison. One can see that the multiuser interference degrades the estimation performance significantly. Hence, although the blind method developed for the single user case can be used to colocated MIMO-OFDM as in [10], it is not applicable to distributed MIMO-OFDM where different users have different Doppler scales.Illustration of the objective functions of the null-subcarrier-based method in a two-user system
(a)
Successful case(b)
Failed case


## 7.  Conclusion

This paper compared different methods for Doppler scale estimation for a CP-OFDM preamble followed by ZP-OFDM data transmissions. Blind methods utilizing the underlying signalling structure work very well at medium to high SNR ranges, while cross-correlation-based methods can work at low SNR ranges based on the full or partial knowledge of the transmitted waveform. All of these methods are viable choices for practical OFDM receivers. In a distributed multiuser scenario, cross-correlation approaches are more robust against multiuser interference than blind methods.



---

*Source: 703243-2012-05-30.xml*