Mathematical Problems in Engineering

Volume 2015, Article ID 719025, 7 pages

http://dx.doi.org/10.1155/2015/719025

## Symbol Estimation Algorithm for MIMO Underwater Acoustic Communication System Based on Multiplicative Noise Model

^{1}College of Engineering, Ocean University of China, Qingdao 266100, China^{2}Department of Education, Ocean University of China, Qingdao 266100, China

Received 25 September 2014; Accepted 22 January 2015

Academic Editor: Hsuan-Ling Kao

Copyright © 2015 Ling Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The stochastic and time-varying underwater acoustic (UWA) channels are usually affected by serious multipath delays, energy loss and distortion factors, thus making the modeling and estimation of the UWA channel challenging problems in the research community. Based on the analysis of the UWA channel, the system with multiplicative noise (SMN) model is established to characterize the complicated factors such as random time-variation, nonlinearity, and energy attenuation. As to the multiple-input multiple-output (MIMO) UWA communication, the complicated SMN model is established for MIMO UWA channels; based on which, the transmitted symbols are estimated according to the optimal recursive filtering algorithm. The algorithm is derived based on the projection theorem, which is optimal in the sense of linear minimum variance, and can overcome the intersymbol interference and noise pollution efficiently. The optimal algorithm is computed recursively, which has the advantage of computation-efficiency and can track the random variation of the fast time-varying channel gain dynamically. Simulation results have validated the effectiveness of the algorithm. The model and the algorithm can be extended flexibly to certain practical problems, such as the joint channel and symbol estimation in underwater acoustic communication systems.

#### 1. Introduction

Underwater acoustic (UWA) communication usually works in very complex transmission environment with seriously noisy background and very limited bandwidth. In addition, the complex multipath propagation and time-frequency selective fading are prone to causing serious distortion of UWA communication signals, which set obstacles to the high-speed and reliable transmission of UWA communication. Hence, how to characterize the underwater channel and estimate the transmitted symbol have become the focus of current research [1–4].

In recent years, there is an increasing demand for high data rate underwater acoustic communication. Thus, some advanced and sophisticated technologies from land wireless communication are introduced into the field of underwater coherent communication, such as multiple input multiple output (MIMO) technology [5–8] and orthogonal frequency division multiplexing (OFDM) technology [9–11]. MIMO technology especially fully exploits space resources and transmits signals from various transmitting and receiving apparatus occupying the same frequency band, without additional receiver transmitting power and spectrum resource; therefore it can greatly increase channel capacity and reduce the bit error rate (BER) of system.

Underwater acoustic channel is an integral part of UWA communication system. Considering the various constraints and the time-space-frequency variability of the channel, UWA channel is much more complex than the radio channel. Usually, there exist reflection and refraction in the propagation of acoustic wave in seawater as marine environment is inhomogeneous, which results in a lot of acoustic signals superimposed in the receiving end and the so-called multipath effect leads to intersymbol interference (ISI) and signal amplitude fading. Currently, the adaptive equalization technology is commonly used to overcome this interference [12–15]. But the equalization algorithm cannot perform well under complex and random time-varying communication environment.

Recently, the research on the multiplicative noise in UWA communication has become a hot spot [16–19] because the multiplicative noise can describe many factors in practical systems such as nonlinearity, packet dropout, and frequency offset. Considering the various kinds of interference in UWA channel, it can be characterized by a system with multiplicative noise (SMN) model. In [20, 21], the UWA channel is seen as SMN with additive and multiplicative noises, and the underwater target is regarded as a special scatterer. Therefore, the problem of underwater target detection can be converted into the estimation of the scattering coefficient sequence constituted by the target and the scatterer nearby. Conventionally, the SMN model is applied in single input multiple output (SIMO) UWA system and seldom considers the MIMO system. This paper focuses on the modeling of MIMO UWA channel and derives the optimal estimation algorithm based on the new model, which can be applied for more complex UWA communication system.

In this paper, based on the characteristics of the UWA channel, the multiplicative noise is introduced to model the MIMO UWA communication system in order to track the fast time-varying UWA channel dynamically; that is, a novel state-space model with multiplicative noise is established to describe the channel variation and achieve effective estimations of the transmitted symbols. Based on the system model, the optimal recursive filtering algorithm is derived for complicated multichannel system with multiplicative noise in the form of general stochastic matrix, which can solve the optimal estimation problem of systems under more complicated environment. The filtering algorithm is optimal under the linear minimum variance criterion and is computed recursively, which is suitable for high speed UWA communication. In the modeling process, the transmitted symbols are included in the state vector. Therefore, according to this correspondence between the state vector in the system model and the transmitted signals from various transmitters, the symbol detection problem can be converted into the filtering estimation of the system state.

The paper is arranged as follows. The SMN model is established to characterize the MIMO UWA channel in Section 2. In Section 3, the symbol estimation algorithm is developed for MIMO UWA system. The performance of the established model and the present algorithm is analyzed in Section 4. Simulation results are shown to validate the algorithm in Section 5. Finally, conclusions are drawn in Section 6.

#### 2. MIMO UWA Channel Modeling Based on SMN

For an MIMO UWA communication system employing transducer sources and hydrophone receivers, the discrete-time baseband equivalent signal received at the th receiver is expressed bywhere is the transmitted symbol from the th transducer at time instant ; is the th fading channel coefficient between the th transducer and the th receiver at time ; is the maximum length of channels; and is the zero-mean additive white Gaussian noise on the th receiver.

If the UWA channel has ideal transmission characteristics, the slowly varying and stable compositions will be the dominant channel characteristics [22]; thereby the channel parameters remain unchanged during a packet or even few packets of data; that is, . Thus (1) can be rewritten as

Otherwise, for underwater remote communication systems, UWA channel usually has more complex transmission characteristics [16]. Signals from transmitters along different paths arrive at the receiving end at different times. Therefore, the received signal is the superposition of all these signals, causing waveform distortion. This effect of UWA channel includes not only marine environment noise and target radiation noise as additive interferences, but also a variety of distortions and information energy attenuation as multiplicative interferences due to the random time-varying characteristics of UWA channel. Broadly speaking, those factors that bring the received signal multiplicative interference can be characterized by the multiplicative noise. Therefore, in this paper, based on the random characteristics of the UWA channel, the multiplicative noise is applied to characterize the complicated and time-varying UWA channel, and the convolution model with multiplicative noise in the received signal can be expressed aswhere the multiplicative noise is of Gaussian distribution with the mathematic expectation and known covariance.

From (3), let

For underwater acoustic communication system, define

From (3)–(6), we have the following SMN model:where

#### 3. Symbol Estimation Algorithm of MIMO UWA Communication System

During the process of SMN modeling, we can see that , . For binary phase shift keying (BPSK) modulation system, the transmitted symbol is taken independently and equally from ; hence, it is clear that the mean of is 0 and the variance is 1. Based on the SMN model (7), we can obtain the state estimation according to the optimal filtering algorithm. Because of the correspondence between the transmitted signals of transducer sources and the state vector , the estimation of the transmitted signals can be acquired easily. For convenience of understanding, the algorithm is listed as follows.

*The Optimal Recursive Filtering Algorithm for SMN*. Considerwhere denotes Hermitian transpose of a vector or matrix, is the th element of matrix , and ;

#### 4. Performance Analysis

##### 4.1. Optimality

The recursive filtering algorithm proposed in this paper is based on the well-known projection theorem [23], which is optimal in the sense of linear minimum variance. Therefore, the theoretical performance of the estimation algorithm is guaranteed such as the optimality and the convergence rate, which will be illustrated by the simulation examples in the next section. Furthermore, the model and the corresponding algorithm can be more flexibly extended to solve some practical systems under complex environment interfered by multiplicative noises, such as the symbol estimation problem of the MIMO UWA communication system with packet dropout and the estimation of Doppler frequency offset in MIMO OFDM systems.

##### 4.2. Recursive Computation

The algorithm is computed recursively with computational complexity , where is the number of the transducers, that is, the dimension of the state vector. Since the algorithm has low computation cost, it can track the fast time-varying UWA channel dynamically, which allows the UWA channel gain changing even over one symbol interval, while the equalization algorithm cannot. The convergence rate of the proposed algorithm will be verified by the simulation in the next section. The above advantages of the new algorithm will better adapt to the demand of the future high-speed underwater acoustic communication system.

#### 5. Simulation and Comparison

It is known that UWA channels change randomly over time, so we have some simulation results for the optimal recursive filtering algorithm when the multiplicative noise is a general stochastic matrix with every component related to each other at the same time index. In addition, based on the characteristics of UWA channel, we assume that it obeys Rayleigh distribution. Since the multiplicative noise can describe the fading, nonlinearity, and distortion of the channel, it is reasonable to model the Rayleigh fading channel. The received signal obtained from the multipath Rayleigh fading channel is estimated by the proposed optimal filtering algorithm in this paper and the adaptive equalization algorithm from [12], respectively, and the BER performance of these two algorithms is compared under the same signal to noise ratio (SNR) in Figures 4, 5, and 6.

##### 5.1. Simulation of the Proposed Optimal Filtering Algorithm

For UWA communication system with 2 transmitters and 2 receivers, the least square (LS) estimation of the channel impulse response can be obtained by the method in [24]. Set . Based on the characteristics of UWA channel, we presume that the mathematic expectation of the multiplicative noise is 0.9; the variance and the covariance of each component in are 0.02.

Figure 1 depicts the plot of BER versus SNR by BPSK modulation using the proposed algorithm to estimate transmitted signals from the two transducer sources. Figure 2 shows BER versus SNR by quadrature phase shift keying (QPSK) modulation from the two transducer sources using the present algorithm.