New Wideband Beam-forming Method Used in Underwater Communication System

SU Wei(苏为)，HUANG Xiao-yan(黄晓燕)，CHENG En(程恩)，YUAN Fei(袁飞)，SUN Hai-xin(孙海信)

(Key Laboratory of Underwater Acoustic Communication and Marine Information Technology，Ministry of Education，Xiamen University，Xia’men 361005，Fujian，China)

Introduction

Some phase coherent modulation and demodulation techniques，such as PSK，can improve data rate greatly in underwater acoustic(UWA)communication，because the bandwidth of UWA channel is seriously limited.But compared with those incoherent techniques，such as FSK，they need higher SNR of the received signals for decoding the signals［1］.The use of beam-forming technique and the SNR gain obtained by spatial filters are demonstrated as an effective method to improve the communication performance if SNR is relatively low.And，switching-beam algorithms or adaptive beam-forming algorithms must be used in the receiver to track the direction of arrivals(DOA)of the incident signals.Compared with the adaptive beamforming algorithms，the switching-beam algorithms are more robust［2］.And，in order to receive undistorted PSK signals，a time domain wideband beam-forming structure with constant beam-width is proposed in Ref.［3］.It combined tapped-delay-lines with real coefficients FIR filters.And the real coefficients FIR filters are designed to realize amplitude and phase weighting at different frequency points.But，simulation results show that the real coefficients FIR filters proposed in Ref.［3］are unable to meet the requirements for complicated frequency response and the communication system with a large ratio of communication band to carrier frequency.In order to solve this problem，a time domain wideband beam-forming structure based on complex coefficient FIR digital filters is proposed in this paper.It is more suitable for PSK underwater acoustic communication than the structure proposed in Ref.［3］.

1 Design of Weights for Wideband and Constant Beam-width

An adaptive beam synthesizing algorithm is proposed in Ref.［4］，and an improved method is proposed in literature Ref.［5］.

Afterpth iteration，the designed weights of the al-gorithm proposed in Ref.［4］is

wherefp(θm)is the energy of pre-defined virtual interference source.

From Eq.(2)，it can be seen that the constraint of the mainlobe shape is not taken into account in the original algorithm.Thus，its stability is poor.An improved algorithm replacesV(θd)in Eq.(1)with(Vm+Wp－1)［5］.The formula of beam-forming weights afterpth iteration is

and，Vmcan be expressed as

whereθican be taken equally in the mainlobe，Aiis the constraint coefficient of the mainlobe.

The error between the array shape calculated after(p－1)th iteration and referenced array shapeD(θi)inθidirection is

According to Eq.(5)，the update constraint coefficientAiof the mainlobe is

个体对自身正当利益的追求是无可厚非的，毕竟“任何人如果不同时为了自己的某种需要和为了这种需要的器官而做事，他就什么也不能做”[3]286。 但在市场逐利本质的驱使下，受利益最大化原则的影响，个体极易因一己之私而产生短视、功利的思想和行为，沦为自身私欲的奴隶。 倘若各个利益主体都对一己之私锱铢必较，不惜牺牲他人或集体的利益，那么利益主体间就生成一种互抑状态。 在这种互抑状态下，各个利益主体看似享有争取自身利益的自由，实则经常性的由于多方牵制而难以真正实现自身的利益诉求。 这种多元利益主体之间自发形成的互抑状态若不加以引导和调试，必将激化本就错综复杂的利益冲突。

whereμis the iteration step.

2 Proposed Wideband Beam-forming Structure

2.1 Beam-forming Structure Based on Complex Coefficient FIR Digital Filter

The designed constant beam width weight of thekth sensor at the frequencyflcan be expressed as

whereaklandφklare the amplitude and phase of the weight respectively.

In literature Ref.［3］，the wideband signals received bykth sensor are weighted by real coefficient FIR digital filters with frequency responsesHk(ωl)=wkl .But，simulation results show that the real coefficient FIR digital filters can not meet the design requirements，because the frequency and phase responses of the required filters are complicated in PSK underwater acoustic communication.

Therefore，a time domain wideband beam-forming structure with constant beam width based on complex coefficient FIR digital filters is proposed in this paper，as shown in Fig.1.It compensates the time delay of integer sampling intervals with tapped-delay-lines first，calculates the complex envelopes of the received wideband signals by using spectrum shifting then，and weights the received signals with a set of complex coefficient FIR digital filters finally.

Fig.1 Time domain wideband beam-forming structure based on complex coefficient FIR digital filters

Obviously，the complex coefficient FIR digital filters can realize more complicate frequency responses than real coefficient FIR digital filters with the same order.The signals are processed as follows.

First，assuming that the received signal is the sum of a group of different frequency CW signals，thus，the compensated signal ofkth sensor can be written as

whereτkis the delay time introduced bykth tapped-de-lay-line.

Then，after spectrum shifting，the complex representation of the received signal envelop which contains high-frequency component can be expressed as

The complex envelopes of the received signals ofkth sensor can be achieved by low pass filter.

Finally，the signalsk(t)is weighted by complex coefficient FIR digital filters with responses shown in Eq.(7)，and the beam output signal is achieved by weighted sum of different sensor signals.

2.2 Design Algorithm of Complex Coefficient FIR Digital Filters

It can be proved that，if the input signal of filter is the complex envelop of received real signal，the frequency response of the complex coefficient FIR digital filter can be expressed as

wherehnis the coefficient of the designed filter，Nis the order of filter，andfsis the sampling rate of the system.

The adaptive［6］or optimization method［7］can be used to design the coefficients of the required FIR filters.For the frequencyfloutside the pass band，if using the adaptive method，the pre-defined frequency responses，including amplitude and phase，must be provided.It increases extra computational volume and has disadvantages of aberration and bad robustness.While using the optimization method，the amplitude frequency response at different frequencies will be constrained only.Therefore，it is more robust than the adaptive method.

The optimization method can be described as follows.

First，a group of frequencies in the range of［0，fs/2］is selected to make the frequency interval smaller in the communication band and larger outside the band.

Second，the optimization relational expression can be established by normL∞for the frequencyflin the pass band.

And the optimization relational expression established by normL2is

For the frequencyfloutside the pass band，it is necessary to add a constraint condition

whereais a constant used to constrain the side-lodes of FIR filters.

Finally，solving the above optimization problem can be solved.

In this paper，the software Sedumi and Yalmip are adopted to solve the above optimization problem.

3 Simulation Results

In the simulation，a volumetric array composed of 12 sensors，as shown in Fig.2(c)，is adopted.The communication band is 4 －10 kHz.The intersection angle between the direction of main beam and the base array normal is 15°.The weights are calculated by using the algorithm proposed in the Section 1 first.Then，the frequency band will be shifted from 4－10 kHz to 0－6 kHz to obtain the complex envelopes of the received signals.Finally，the real and complex coefficient FIR digital filters are used to fit the beam weights，respectively.The simulation results are shown in Fig.2.

Fig.2 Simulation results

Fig.2(a)and Fig.2(b)show the results of adopting 61-order real coefficient FIR digital filter and 31-order complex coefficient FIR digital filter to fit the weights of 0th sensor，respectively.the solid line in figure represents the expected frequency responses，i.e.the weights for wideband beam with constant width.The line labeled by‘* ’shows the fitting result of complex coefficient FIR digital filter.It essentially coincides with the expected value.And，the line labeled by‘△’shows the result of real coefficient FIR digital filter greatly.It differs from the expected value.It can be seen that the shape of the beam weighted by real coefficient filters have great distortion when the ratio of communication band to carrier frequency is larger than 0.85.

The beam shape weighted by the complex coefficient filters are shown in Fig.2(d).It can be seen clearly that the shapes of mainlobes meet the requirements well and the errors are less than 0.3 dB.

4 Conclusions

A novel wideband beam-forming structure with constant beam width based on complex coefficient FIR digital filters is proposed in this paper.The design algorithm of the complex coefficient FIR filters is presented also.The simulation results for a volumetric array composed of 12 sensors in the frequency range of 4－10 kHz show that the proposed structure is more suitable for the underwater acoustic communication situation where there is a larger ratio of communication band to carrier frequency.

［1］LIU Yun-tao，YANG Shen-yuan.Research of adaptive spatio-temporal DFE with embedded DPLL in high-speed underwater digital communication［J］.Journal of Harbin Engineering University，2005，26(5):658 －662.(in Chinese)

［2］Thompson D，Neasham J，Sharif B S，et al.Performance of coherent PSK receivers using adaptive combining，beamforming and equalizations in 50 km underwater acoustic channels［J］.IEEE Trans on Oceans，1996，25(2):845－850.

［3］YANG Yi-xin.Studies on beamforming and beamspace high resolution bearing estimation techniques in sonar systems［D］.Xi’an:Northwestern Polytechnical University，2002:26－35.(in Chinese)

［4］Olen C A.A numerical pattern synthesis algorithm for arrays［J］.IEEE Trans on Antennas and Propagation，1990，38:1666－1676.

［5］Anh P.A numerical pattern synthesis algorithm for arbitary arrays［J］.Antennas Propagation and EM Theory，2003，28(2):210－213.

［6］HE Bing.Designing oflinear phase FIR filters using modified NLMS algorithm［J］.Journal of Shanghai Jiaotong University，2000，(2):266 －268.(in Chinese)

［7］YAN She-feng，MA Yuan-liang.Optimal designing of space-time filters using Second order cone programming［J］.Science in China Series E:Information Sciences，2006，36(2):153－171.(in Chinese)