Advanced Coding For Underwater Communication
Advanced Coding For Underwater Communication
1. INTRODUCTION
In this section we present some of the theoretical concepts of the Pattern Time
Delay Shift Coding (PDS) technique.
The major challenge of hydro acoustic communication in underwater
horizontal acoustic channels is the multipath interference. There are two
aspects of the multipath interference presented in the underwater acoustic
channel: code's multipath propagation time delay spread and received wave
pattern distortion by multipath interference.
The model of the impulse response function of the coherent multipath
channel is given by:
I
N·!
h( T) = Aob( T- To)+ Aib( T- T i ) (1)
i=!
Where Ai and Ti are sound ray parameters through the receiver
corresponding to the amplitude and arrival delay time respectively.
The sound ray determining the impulse response function is named acoustic
Eigenray.
If signal transmitted from the sound source is denoted by z (t), the received
signal in the multipath channel is given by:
N-!
S( t ) = AoZ( t - To ) + I AiZ( t - Ti ) + n( t ) (2)
i=!
Where the first term on the right is the direct wave, the second are refracted
and reflected with boundaries. When the signals of multipath arrivals are
overlapped with the direct arrival at the same time, it will derive interference,
so the waveform of composite signal is different from the transmitted signal.
The third term is noise; including environment noise and ship self noise.
When the delay time difference between the multipath ray and the direct
ray is larger than the code's duration of a symbol, the adjacent codes (symbols)
will be overlapped and interfered to cause "intersymbol multipath
interference". When the delay time difference is smaller than the code's delay
time resolution, "inner symbol multi path interference" will be introduced.
The simplest method to overcome the intersymbol multipath interference is to
introduce enough rest-periods between transmitting of adjacent symbols, that
is, the symbol interval To is larger than the multipath delay spread Td .
However, such method results in a low communication rate.
The other approach is to use Frequency Shift Keying (FSK) method. In this
approach the relationship between the numbers of used frequency points and
229
the code's duration is typical given by:
M =11 Td 11+1 (3)
To
Where II. II is round-off operator, To is the code's duration, Td is maximum
delay spread.
At the receiver system, M filters are used to separate the symbols. In the
time interval of To, every filter only outputs the code's direct wave and
multipath arrivals of same corresponding frequency, so the intersymbol
interference is overcome. Although, the FSK is a well known reliable
communication system used for multipath channel conditions, but it still
requires broad bandwidth especially for the high bit rate digital
communication conditions.
The PDS separates symbol with some pattern codes. At the receiver end,
copy-correlators are used to separate the codes. To every pattern code, the
copy-correlator's reference signals are different. The copy-correlator outputs
a large correlation peak only when the received wave pattern is the same with
reference wave pattern. If the cross correlation coefficients between every
two codes are low enough, the copy-correlator can choose the corresponding
code. In the time interval of To, the copy-correlator outputs the peak of one
pattern code's direct wave and multipath waves, so the intersymbol
interference is overcome.
In PDS system design, the key to suppress intersymbol multi path
interference is to design a set of pattern codes, the cross correlation
coefficients of which are low enough. If there are M kind of pattern codes, the
delay time spread can be overcome is Td , when Td =M To. The PDS system
uses pattern to separate codes, so it uses narrower frequency bandwidth, and it
is more reliable than the FSK system. The PDS system codes with the delay
time difference between the information code and its correcting code.
The multipath interference problem can be mitigated by well-known
channel compensation technique such adaptive channel equalization
technique. Such techniques are effective to overcome the intersymbol
interference and inner symbol interference to a certain extent. Their recovery
of the transmitted signal from the received signal is via estimates of the
channel's impulse response. However, this will decrease the influence of
multipath interference and decode more correctly. Hence, neither of the above
methods is effective to cancel the effects of inner symbol interference with
extremely small delay time difference (multipath opposite-phase interference,
i.e., the direct wave is 1t phase shift to the reflected wave).
230
code 2 n
~ -.0 ~ Block N
The first code is named frame synchronization code, which provides time
reference for decoding. The first block codes "Block 0" are correcting codes
that provide time delay reference and detection threshold reference for
corresponding pattern information recoding. The correcting codes may
suppress the measuring error of time delay difference due to the inner symbol
multipath interference. It's noted that the correcting codes have the same
pattern with their corresponding information codes, so time delay deviation
due to inner symbol multipath interference in the slowly varying channel is
the same. If the correcting codes are used as time delay reference, the inner
symbol multipath interference has little influence to measuring the time delay
difference. The block codes after 'Block 0" are the information codes used for
transmitting data or messages. The code pulse width is Tc ' the time width used
for coding information is 1'" and the code interval is To = ~ + T" the time
delay 'difference between information code and its reference code represents
the messages. The number of pattern of frequency modulated (PM) codes is
M in every block. The cross correlation coefficients between every
two-pattern codes must be low enough. Their normalization cross correlation
coefficients must be less than 0.35 at least.
The general block diagram of the PDS digital communication system is
shown in Fig.2.
231
"~~i~~"~····*"·
source + Channel I
shift
The source coder produces data sequence. The channel coder modulates
them into PDS time delay shift codes relatively to their correcting codes. The
received signal is distorted by multipath interference and noise. Before the
beginning of received message, copy-correlator searches the frame
synchronization signal. Once the PM signal is detected, the copy-correlator
opens a time window immediately and produces the corresponding pattern
reference signal in every time window. The copy-correlation detects the
code's correlation peak from the noise background. The detected code has the
same pattern with the reference signal, so the time delay difference between
the detected code's correlation peak and its correcting code's correlation peak
can be precisely measured. "Block 0" provides the time delay difference
reference for every pattern code. The information codes' time delay shift
sequence is translated into data sequence or messages. For every pattern code
of "block 0", the amplitudes of output peak of the copy-correlator are
different. They are time-varying and influenced by inner symbol multipath
interference and the channel's propagation loss. That the pattern codes'
detection threshold value can be determined by the corresponding correcting
codes' correlation peak multiplying a constant, which is smaller than one, it is
named adaptive threshold. The adaptive threshold can suppress intersymbol
interference and noise. Because every two pattern codes' cross correlation
coefficient is low enough, the adaptive threshold can suppress the
intersymbol multipath interference.
232
Since the PDS communication system uses time delay shift as information
coding, hence the time delay difference estimation represent one of the key
techniques in the decoding process of the system. The higher is the estimation
precision required, the higher would be the communication rate.
This process may be derived from the following formula:
The PDS system communication rate is given by [22]:
1 -1 R( k - 1) + R( k + 1)
OJ=-cos (8)
rs 2R( k)
U7 tillS}
0l)4 . - , - - -
Q1l:!
0lr1
.aD1
.(Ul2
.o..tXl
.0.,:)< t i n
D s:.o ~ '1IXl %lXl
fig.(3) Time delay estimation results for :WOO times in dependent experiments
(cosine interpolation time delay estimation simulation result)
Hence, one block of correcting codes is needed to transmit during the channel
coherent time length. In generally, it is only needed to transmit the correcting
codes one times in a few ten seconds if the relative speed between the
transmitter and receiver is not high. The other factor to influence the time
delay estimation is the relative motion between the transmitter and the
receiver. Motion can fasten the channel's varying, so correcting code
transmitting must increase correspondingly. Doppler effect due to relative
motion will introduce time delay measuring deviation. LFM signal's
ambiguity function is coupled to Doppler and time delay. For different value
of Doppler, the time delay corresponding to the correlation peak position is
different. Because the correcting code and the information code have the
same pattern, their time delay estimation deviations caused by Doppler are the
same and the time delay difference deviation can be eliminated.
,,I
1\
1 I
I
'\J
\ -I
1 '\
\
I
1
1 \
\,t.,p"iifill1ii61ijWffi..d'filOW"[Link].dUtlnC6ou"[Link],il:,,".J 1 ,\1
.0
BOTTO""" REF1...E~
---r---.--':'!
.0
I, I I 11 _ _ 20.0""
_ ::0.0,....
:.1 0","_
I
",7,....
I' DI- DD.C
1I't_
1'\0-
::tOO..,...
2:JO.O~
In the Fig. above the upper window on the left displays a ray trace diagram
which is a set of ray trace transmitted from the source with equalized angle
interval to understand the acoustic field's basic characteristic. The eigenray
cluster is shown in the middle window on the left. The eigenrays are a series
of special sound ray passing through the receiver which the parameters (Ai> 'ti»
determine the channel's impulse response function are shown in the lower
window on the left of the Fig. The calculation error of time delay of all
eigenray is less than 3/ls. The number of egienray is related to the following
effects: transducer's directivity, geometric position of the transmitter and
receiver, sea depth and profile of sound velocity etc. Once the channel
impulse response function is calculated, the simulation program can predict
the received signals form for the transmitted source. The sound velocity
profile is shown in the upper window on the right. The simulation program
enables to input actual measuring sound velocity profile, so communication
quality in any condition could be observed. The middle window on the right
displays bottom reflection coefficient curve. The simulation program enables
to input arbitrary bottom reflection characteristics. The simulation parameters
are listed in the lower window on the right. All parameters could be changed
to analyze the communication quality. Some of the simulation results of this
simulated model are shown in Fig.5 and Fig.6 respectively.
.•
237
~ ~~
~L""""''''''''LE
. ,.eo 1~10
'-~~"~~=====~~:: !
eo "
1.00 ...--.::.r1····I---n'
:::: :::: =-->····I'
[Link] ···I-=wr-···'
::! ::::..I~··I
...... ·~I
::::::: ::..I·'·-----=I·
.... ·..·.. ·· =
:~: :::::!:::::::: ..::!::::
·· .. ·..::: :::!::::: ::::ru··! \0' ._ ........... ~ ............ _\
=..-·,·
-......
-·· 1 "OTTO'" REFLECf
....- ~O.D",
.... - 2C.O",
,~& Illrl':~j:II:I:li~~;;·~~~wrn
...riL...... ,.............. .. I
......-- , '.0",
_1.De
........ 4 .. .:::::: . :::::: .:::::
j ..•
:::;:. :::::.
•
::::::::. :: :~: •. :::::;::::::::::;::::::::::::::::::::
R1:.--~oo.O!"'l"\
[~I__llilN~W~~·~m~'.~
o 18~ ........ s~- 11.0dB
~,,
Ct: 1.81
Cr-: ,.el
coo.: CI
I~I~"--------"~---"-----~---~J~o~
~~~~~~ec l~ "
!
'-OO ~
.:::1::::1:::::::1::1::1::::::::::::::::::::::::;::;::;::::::I!' ~.9C!fT"OU
./'::::I~ITI:::
•••• _ - . . . '" •••• _ ......._ •• _ ••. _ .......................... _ ••••••••• _ .... [Link]
•• } • • • • • • • • • • • • • • u
....... .. ~
~~:~ ~:~: ~:~~ ::~~ .~:: :::~:~: :~::::y:~:~:::I::::::::L,~~~:~ :::~'~~1~~~::
.Jl
_ \ .00 ( ::
I '" 11!;I.I~I;lllli~~~~.1
r.- 20..otn
Rv-D-OO..om
-...- ,0
-2.13 C 'I:I~""" ~- [Link]
~""'"
11iIRmiwill~:j,~~i
.. ~ .. __ 1
Ct: ,eel
Cr: l.e1
1"::: o
I
fig. (6) PSD acoustic communication result of simulation
(inner symbol constructive multipath interference)
In these Figs. (from top to bottom), the left windows show eigenray
diagram, transmitted PDS code signal, received wave disturbed by noise and
238
multipath interference and output signal form the copy-correlator respectively.
while the right window represent the sound velocity profile, bottom reflection
coefficient and main parameters of the simulation system.
Fig.5 and Fig.6 show the simulation results of inner symbol multi path
interference. The former considers the inner symbol destructive multipath
interference, while the later considers the inner symbol constructive multipath
interference. On the account of above, the former received signal is not more
regular than the later, and its amplitude is smaller. The copy-correlator output
peak amplitude of the former is smaller, and its correlation peak waveform is
distorted as shown in Fig.7. The simulation results of above channel
conditions are listed in the Table 2 and Table 3 respectively. It can be
concluded from these results that in the inner symbol destructive multipath
interference case, the time delay estimation deviation is increased
considerably with a corresponding error rate increase. In order to reach high
reliable communication effect in such channel, higher SNR is needed.
J-_.
H(Y)
W
~f o
=c ,)'!-\cvy>o
w ~
A? =-'v\rvv. . ~~'""--
~ ~ ~
c /v,
~
--=~
~
Y(ms)
Table 2 Time delay estimation deviation and error rate in inner symbol
destructive mul --. ---- ------------- -------"
,-----~-------
Table 3 Time delay estimation deviation and error rate in inner symbol
Itioath interference (simul .
292 bitls 383 bitls 453 bitls
~
SR a,(ms)
17dB -.. . . 0.003248 I 0
I Pe o,(ms) I Pe
0.003113 I 0
o,(ms) I Pe
0.003182 I 0
239
In this section the experimental set up and the underwater test results of the
system are presented. The aim of these underwater tests is to evaluate the
performance of the PDS communication system and to study the influence of
the tested underwater channel on the performance of the system. The block
diagram of the experimental test configuration is shown in Fig.8.
Wide
Compu Pattern Signal Compu
-band Hydro
-ter genera -ter
trans phone
(1 ) -tor (2)
-ducer
R(t)
D.B
~. ~
o·
o
-~-
$
)~~10 15 20 2$ 1I5 4D
t(ms)
, J\~, ,
10 15
--~
_rv-f
20
! -
2$
~([Link]
30 1I5 4D
t(ms)
DB J
Q.8
0-4
D.2
DE crJ"V\o=c ocyJVVV"V\!V'..,,,,,,, . ."c",r- it(ms)
o 10 15 20 2$ 30 35 4D
The Figs above show the correlation peak of the frame synchronization code
varying with the distances tested. The distances (270m, 460m) with high
correlation peak as shown in Fig.12 correspond to the case of inner symbol
constructive multipath interference. For the test distance (300m) with low
correlation peak as shown also in Fig.12 correspond to the case of the inner
symbol destructive multipath interference. In the condition of destructive
multipath interference the maximum error rate were at the level of 10"4.
242
Rmax(t)
,eo
I'Gr......."'~"
150 ,
120 . . . ~~~~
., I
90
so
ro ~-
t(s;
o ~I------~----~------~----~------~------~----~----~
o 100 200 XX) 400 500 600 700 BOO
~~~~
2 ~t----------------------------------------------------~
'5~~vJ\ I
0.5 .-
\
o ~I----~--~------------------------~------------~
o 200 4O.J EM) 800 10Xl 1200 1400 1600 1800 range(m)
The results shown in Fig. 13 confirm that that the simulation results of the
system's model correlate well with the corresponding results of the actual
underwater trial conditions. Fig.I4 and Fig.I5 give the trial's correlation
output waveform every second during a testing period of lOs. This illustrates
that the channel is stable enough even if the receiver ship drifts at speed of 3
knot.
243
! ~
<J)
-.
6
§'
0
bh
'CXI
~ 4)
~
~ -1 a
}
~
110
0
I-<
-«iI::
0
4)
0
U
--
2D
fig (14)
"<:t
.......
." . ." ."
'-'
E
•. .. OIl
R(tJ
0
0
d 0":1 d•
:~
a:: ::I " tH
D
0
D..I
0.-
,
0._
0
"I
0.1
::I
0
D.B
D._
<>
::I
D
d
D-
O
D..I
o
::I
D-
D
D.B
D
0
d
D-
<>
D.B
D
d
D.-
D.-
D.B
~ d 0
D.I
~
D.B
0._
~
0
0
~ " 0 Q 0
244
R(t)
. I\
f,
o.~ l-
0.' I-
~ rJ \r\~
:
o
_A~
o.e
0.' I-
o
~J\~~
0.1 I-
0.' I- 0" •
o
0.1
0.'
o
-
_~J\~,_
-
~J\A.~
0.8
0.' l-
o
~J\~~
O.B I-
0.' -
,
o
o.e I-
0.4
a
- _J~~
-
_r\!\~r-
0.8
0.4 I-
o
_.J\,\~ _
0.8
-
0.4
o
-
_J\~~
0.8
04
o t(ms)
o 10 \2 14
l' II 31
6. CONCLUSIONS
7. REFERENCES