JP2006109279A - Underwater communication system - Google Patents

Underwater communication system Download PDF

Info

Publication number
JP2006109279A
JP2006109279A JP2004295516A JP2004295516A JP2006109279A JP 2006109279 A JP2006109279 A JP 2006109279A JP 2004295516 A JP2004295516 A JP 2004295516A JP 2004295516 A JP2004295516 A JP 2004295516A JP 2006109279 A JP2006109279 A JP 2006109279A
Authority
JP
Japan
Prior art keywords
communication system
synchronization
underwater communication
channel
synchronization point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2004295516A
Other languages
Japanese (ja)
Other versions
JP4341522B2 (en
JP2006109279A5 (en
Inventor
Kunio Hashiba
邦夫 橋場
Fumitaka Maeda
文孝 前田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2004295516A priority Critical patent/JP4341522B2/en
Publication of JP2006109279A publication Critical patent/JP2006109279A/en
Publication of JP2006109279A5 publication Critical patent/JP2006109279A5/ja
Application granted granted Critical
Publication of JP4341522B2 publication Critical patent/JP4341522B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide an underwater communication system in which stable communication can be performed even in a multipath phasing environment. <P>SOLUTION: In the underwater communication system employing QPSK modulation, a transmitting side adds the same synchronizing symbol to each of in-phase component and quadrature component of a QPSK signal in the head of a QPSK transmission data frame and transmits it, in a receiving side, a correlation operating unit 10 uses a quadrature detecting unit 12 to perform quadrature detection on a received signal to determine an auto-correlation function (f) between the in-phase component and the synchronizing symbol after quadrature detection, and an auto-correlation function (g) between the quadrature component and the synchronizing symbol after quadrature detection, and the square of the auto-correlation function (f) and the square of the auto-correlation function (g) are added to output a time function (17). In a synchronizing point detecting unit 20, the position of a peak value of the time function is detected as a synchronizing point, frame synchronism is established and in a QPSK demodulation unit 50, QPSK symbol data are outputted. Thus, stable synchronism can be realized. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は超音波を伝送媒体とした水中通信方式に関し、特に浅海域における多重反射等の影響に対して、安定した通信を行うことが可能な水中通信方式及び水中通信システムに関する。   The present invention relates to an underwater communication system using ultrasonic waves as a transmission medium, and more particularly to an underwater communication system and an underwater communication system capable of performing stable communication against the influence of multiple reflections in a shallow sea area.

従来、海水中でのデジタルデータ通信では、超音波を伝送媒体とした通信システムが用いられている。海中では、送信側から受信側に直接到来する直接波の他に、海面や海底等に反射して到来する遅延波が存在する。特に浅海域では、音波の距離減衰も小さいため、大きなパワーを有する多数の遅延波が生じ、データの伝送エラーを生じさせる大きな原因となる。以下、従来技術の水中通信システムについて、図9〜図11を用いて説明する。   Conventionally, in digital data communication in seawater, a communication system using ultrasonic waves as a transmission medium has been used. In the sea, in addition to the direct wave that arrives directly from the transmission side to the reception side, there is a delayed wave that is reflected and arrives at the sea surface or the sea floor. Particularly in the shallow water area, since the distance attenuation of the sound wave is small, a large number of delayed waves having a large power are generated, which causes a data transmission error. Hereinafter, the underwater communication system of a prior art is demonstrated using FIGS. 9-11.

図9は、従来技術における水中局と水上局からなる水中通信システムでの超音波の伝搬を説明する図、図10は、従来技術における水中通信システムでの実験用送波音圧の振幅スペクトル例を示す図、図11は、従来技術における水中通信システムでの受信音圧の振幅スペクトル例を示す図である。   FIG. 9 is a diagram for explaining the propagation of ultrasonic waves in an underwater communication system composed of an underwater station and a surface station in the prior art, and FIG. 10 is an example of an amplitude spectrum of an experimental transmission sound pressure in the underwater communication system in the prior art. FIG. 11 is a diagram illustrating an example of an amplitude spectrum of received sound pressure in an underwater communication system according to the related art.

図9において、水中局70で収集された海底地形データや物体探知データは、超音波の周波数帯域の変調信号に変換され、水中局70に備えられた送波器71から送信されて、水上局80の受波器81で受信される。受信された信号は水上局80でデジタル処理され、無線アンテナ82を介して、別の水上局や陸上の基地局に伝送される。受波器81で受信される超音波には、実線のパスRdに沿って、受波器81に直接到来する直接波の他、破線のパスRsのように海面1で反射して到来する遅延波、一点鎖線のパスRb1のように海底1で反射して到来する遅延波、二点鎖線のパスRb2のように海底2と海面1で反射して到来する遅延波等、多数の遅延波が存在する。これら遅延波は直接波に重畳して、受信信号の振幅や位相を変化させる、いわゆるマルチパス・フェージング環境を作り出すため、データの伝送エラーを生じさせる原因となる。   In FIG. 9, the seafloor topography data and object detection data collected by the underwater station 70 are converted into modulation signals in the ultrasonic frequency band, transmitted from the transmitter 71 provided in the underwater station 70, and It is received by 80 receivers 81. The received signal is digitally processed by the surface station 80 and transmitted to another surface station or a land base station via the radio antenna 82. The ultrasonic waves received by the receiver 81 are reflected by the sea surface 1 as shown by the broken line Rs in addition to the direct wave that directly arrives at the receiver 81 along the solid line Rd. A large number of delayed waves such as a delayed wave that is reflected by the seabed 1 such as a one-dot chain line path Rb1 and a delayed wave that is reflected by the seabed 2 and the sea surface 1 such as a two-dot chain line path Rb2 Exists. These delayed waves are superimposed on the direct wave to create a so-called multipath fading environment in which the amplitude and phase of the received signal are changed, which causes a data transmission error.

マルチパスの特性を表す例として、深さ5mの無響水槽内において実測したマルチパス・フェージング特性の計測例について説明する。実験では、深さ4mの位置に無指向性の送波用ハイドロホンを設置し、送波用ハイドロホンから水平距離6.6m隔てた位置に、無指向性の受波用ハイドロホンを設置した。図10は、送波用ハイドロホンから放射されたパルス音波の音圧振幅スペクトルであり、図11(a)〜図11(e)は、深さ66cmの位置から43mm間隔で計測した、受波用ハイドロホンにおける受信信号の音圧振幅スペクトルである。図10、図11(a)〜図11(e)から明らかなように、マルチパスの影響によって、受信位置では多数の不感領域(ヌル領域)が発生し、その発生状況も受信位置によって大きく変化することがわかる。このようなマルチパス・フェージング環境下においても安定した水中通信を可能とする技術の一例として、特許文献1に記載の「多重超音波による水中超音波伝送装置及び超音波伝送方法」が知られている。特許文献1に記載の技術では、海面反射に起因して発生する不感領域(ヌル領域)を、異なる周波数の搬送波を用いて補間することにより、データ受信の確実性を向上させる。   As an example representing multipath characteristics, a measurement example of multipath fading characteristics measured in an anechoic water tank having a depth of 5 m will be described. In the experiment, an omnidirectional transmission hydrophone was installed at a depth of 4 m, and an omnidirectional reception hydrophone was installed at a distance of 6.6 m from the transmission hydrophone. . FIG. 10 is a sound pressure amplitude spectrum of a pulsed sound wave radiated from a hydrophone for transmission. FIGS. 11A to 11E show received waves measured at intervals of 43 mm from a position of 66 cm in depth. It is a sound pressure amplitude spectrum of the received signal in the hydrophone. As is clear from FIGS. 10 and 11 (a) to 11 (e), a large number of insensitive areas (null areas) are generated at the reception position due to the influence of multipath, and the occurrence state varies greatly depending on the reception position. I understand that As an example of a technique that enables stable underwater communication even in such a multipath fading environment, “Underwater ultrasonic transmission apparatus and ultrasonic transmission method using multiple ultrasonic waves” described in Patent Document 1 is known. Yes. In the technique described in Patent Document 1, the insensitivity region (null region) generated due to sea surface reflection is interpolated using carrier waves having different frequencies, thereby improving the reliability of data reception.

また、別の例として、特許文献2に記載の「水中通信方式」が知られている。特許文献2に記載の技術は、スタートパルス、データパルス、エンドパルスで構成された受信信号のFSK(Frequency Shift Keying)変調波を周波数分析し、この変調波の周波数帯域内にスタートパルスが存在しているかどうかを検出して、スタートパルスの存在が確認されると、検出した信号の時間的に分離した周波数分析結果を用いて、コード判定を行う。   As another example, an “underwater communication system” described in Patent Document 2 is known. The technique described in Patent Document 2 analyzes the frequency of an FSK (Frequency Shift Keying) modulation wave of a received signal composed of a start pulse, a data pulse, and an end pulse, and the start pulse exists within the frequency band of this modulation wave. When the presence of the start pulse is confirmed, the code determination is performed using the frequency analysis result of the detected signal separated in time.

一方、装置構成の複雑さを避け、通信の安定性をより向上させるためのデータフレーム同期方式が非特許文献1に記載されている。非特許文献1に記載の技術では、変調方式としてQPSK(Quadrature Phase Shift Keying)変調に特化しており、送受信データを全てデジタルで扱い、送信データフレームのQPSK信号の先頭には、その同相成分と直交成分のそれぞれに、同期シンボルとしてバーカー符号が付加される。受信側では、直交検波後の複素受信信号の絶対値を求めることによって搬送波の位相誤差を除去し、その後バーカー符号の自己相関を利用して、同期点を検出している。これによって、位相同期ループやデジタル位相同期ループを必要とせず、装置構成の簡略化も図れるとしている。   On the other hand, Non-Patent Document 1 describes a data frame synchronization method for avoiding the complexity of the device configuration and further improving the stability of communication. The technique described in Non-Patent Document 1 specializes in QPSK (Quadrature Phase Shift Keying) modulation as a modulation method, handles all transmission / reception data digitally, and includes the in-phase component at the beginning of the QPSK signal of the transmission data frame. A Barker code is added as a synchronization symbol to each of the orthogonal components. On the reception side, the phase error of the carrier wave is removed by obtaining the absolute value of the complex reception signal after quadrature detection, and then the synchronization point is detected using the autocorrelation of the Barker code. This eliminates the need for a phase-locked loop or a digital phase-locked loop, and simplifies the device configuration.

特開平9−55707号公報JP-A-9-55707

特開2004−15762号公報JP 2004-15762 A Y.Choi,J.W.Park,S.M.Kim,Y.K.Lim, "Phase Coherent All-Digital Transmitter and Receiver for Underwater Acoustic Communication Systems", System Theory 2003, Proc. of the 35th Southeastern Symposium on, pp.79〜83(2003)Y. Choi, J. W. Park, S. M. Kim, Y. K. Lim, "Phase Coherent All-Digital Transmitter and Receiver for Underwater Acoustic Communication Systems", System Theory 2003, Proc. Of the 35th Southeastern Symposium on, pp. 79-83 (2003)

図9〜図11を用いて説明したように、特に浅海域では、音波の音速が遅いことと波長が短いことから、激しいマルチパス・フェージング環境が形成される。このフェージングによる音響伝搬路の周波数選択性は非常に細かく細分化され、多くの不感領域が形成される。また、不感領域の変化は水上局と水中局の位置関係に大きく依存する。このため、水中局や水上局が動揺したり移動したりすると、搬送波周波数が不感領域に含まれるような状態が頻繁に発生し、上述したような従来技術では、データを取得する以前に同期が取れず、通信の安定性を損なう可能性があった。   As described with reference to FIGS. 9 to 11, particularly in a shallow sea area, since the sound velocity of sound waves is slow and the wavelength is short, a severe multipath fading environment is formed. The frequency selectivity of the acoustic propagation path due to this fading is very finely divided, and many insensitive areas are formed. In addition, the change of the insensitive area largely depends on the positional relationship between the water station and the underwater station. For this reason, when the underwater station or the water station is shaken or moved, a state in which the carrier frequency is included in the insensitive area frequently occurs. In the conventional technology as described above, synchronization is performed before data is acquired. There was a possibility that the stability of communication could be lost.

また、非特許文献1に記載されているフレーム同期方式では、位相同期ループやデジタル位相同期ループを必要とせず、簡略な構成となっているものの、マッチドフィルタ処理の段階において、バーカー符号の2乗値による自己相関を利用しているため、バーカー符号本来の鋭い自己相関関数特性が活かされておらず、マルチパス・フェージング特性の変化に十分対応できない可能性があった。   The frame synchronization method described in Non-Patent Document 1 does not require a phase-locked loop or a digital phase-locked loop and has a simple configuration, but at the stage of matched filter processing, the square of the Barker code is used. Since the autocorrelation by value is used, the sharp autocorrelation function characteristic inherent in Barker code is not utilized, and there is a possibility that the change of the multipath fading characteristic cannot be sufficiently handled.

本発明の目的は、浅海域等におけるマルチパス・フェージング環境においても、安定した同期を実現し、安定して通信が可能な水中通信方式及び水中通信システムを提供することにある。   An object of the present invention is to provide an underwater communication system and an underwater communication system that realizes stable synchronization and enables stable communication even in a multipath fading environment in a shallow sea area or the like.

本発明の水中通信方式は、QPSK変調を用いた水中通信方式であって、送信側では、QPSK送信データフレームの先頭の、QPSK信号の同相成分と直交成分のそれぞれに同一の同期シンボルを付加して送信し、受信側では、受信信号を直交検波し、直交検波後の同相成分と同期シンボルとの自己相関関数fと、直交検波後の直交成分と同期シンボルとの自己相関関数gとを求め、自己相関関数fの2乗と自己相関関数gの2乗とを加算して得られる時間関数sのピーク値Pの位置を同期点として検出して、フレーム同期をとる。同期シンボルとしては、バーカー符号又はM系列符号を用いる。   The underwater communication system of the present invention is an underwater communication system using QPSK modulation, and the transmitting side adds the same synchronization symbol to each of the in-phase component and the quadrature component of the QPSK signal at the head of the QPSK transmission data frame. On the receiving side, the received signal is quadrature-detected, and an autocorrelation function f between the in-phase component and the synchronization symbol after quadrature detection and an autocorrelation function g between the quadrature component and the synchronization symbol after quadrature detection are obtained. The position of the peak value P of the time function s obtained by adding the square of the autocorrelation function f and the square of the autocorrelation function g is detected as a synchronization point to establish frame synchronization. A Barker code or an M-sequence code is used as the synchronization symbol.

本発明の水中通信方式は、受信側では複数の受信チャンネルを具備しており、各チャンネルにおいてピーク値Pを検出し、各チャンネルのピーク値Pが最も大きいチャンネルにおける時間関数sのピーク値Pの位置を同期点として検出して、又は、時間関数sのピーク値Pと、ピーク値Pの時間を含む一定時間範囲内における時間関数sの平均値Mとの比P/Mを各チャンネルについて求め、比P/Mの値が最も大きいチャンネルにおける時間関数sのピーク値Pの位置を同期点として検出して、同期点で全チャンネルのフレーム同期をとる。   The underwater communication system of the present invention has a plurality of receiving channels on the receiving side, detects the peak value P in each channel, and determines the peak value P of the time function s in the channel having the largest peak value P in each channel. The position P is detected as a synchronization point, or the ratio P / M between the peak value P of the time function s and the average value M of the time function s within a certain time range including the time of the peak value P is obtained for each channel. The position of the peak value P of the time function s in the channel with the largest ratio P / M value is detected as a synchronization point, and frame synchronization of all channels is established at the synchronization point.

本発明の水中通信方式は、連続する2つの送信データフレームにおける同期点の時間間隔と、送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、該当のデータフレームについて送信側に再送要求する、又は、一つ目の送信データフレームの同期点と送信データフレームの送信間隔から、二つ目の送信データフレームの同期点を決定してフレーム同期をとる。   In the underwater communication system of the present invention, when the absolute value of the difference between the synchronization point time interval between two consecutive transmission data frames and the transmission interval of the transmission data frame is larger than a predetermined value, the corresponding data The frame is synchronized by requesting the frame to retransmit the frame, or by determining the synchronization point of the second transmission data frame from the synchronization point of the first transmission data frame and the transmission interval of the transmission data frame.

本発明の水中通信方式は、適応等化器を備え、適応等化器は1チャンネル又はマルチチャンネルの判定帰還形適応等化器とする。   The underwater communication system of the present invention includes an adaptive equalizer, and the adaptive equalizer is a one-channel or multi-channel decision feedback type adaptive equalizer.

以上のような方式によって、浅海域等における激しいマルチパス環境化においても、確実な同期が可能となり、安定した通信が可能となる。   With the above method, even in a severe multipath environment such as in shallow water, reliable synchronization is possible and stable communication is possible.

位相同期ループやデジタル位相同期ループを用いることなく、マルチパス・フェージング環境でも、簡単な構成の同期方式により、安定した同期を実現し、安定した通信が可能な水中通信方式及び水中通信システムを提供できる。   Provides an underwater communication system and an underwater communication system that realizes stable synchronization and stable communication even in a multipath fading environment without using a phase locked loop or a digital phase locked loop. it can.

本発明の実施例では、浅海域における激しいマルチパス・フェージング環境下でも、確実なフレーム同期方式を提供するという目的を、位相同期ループやデジタル位相同期ループを用いない簡単な構成により実現する。また、受波器のアレイ化に伴い、さらなる通信の安定化、特により確実なフレーム同期が可能な水中通信方式を実現する。以下、本発明の実施例を図面を用いて詳細に説明する。
(実施例1)
図1は、本発明の水中通信方式の実施例における受信側での信号処理を説明するブロック図である。
In an embodiment of the present invention, the object of providing a reliable frame synchronization method even in a severe multipath fading environment in a shallow sea area is realized by a simple configuration that does not use a phase locked loop or a digital phase locked loop. In addition, with an array of receivers, an underwater communication system capable of further stabilization of communication, particularly more reliable frame synchronization, is realized. Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
Example 1
FIG. 1 is a block diagram for explaining signal processing on the receiving side in an embodiment of the underwater communication system of the present invention.

本発明の実施例による水中通信では、データ通信におけるデジタル変調方式としてQPSK変調を用いる。図1において、図示しない受波器で受信した変調信号は、図示しないプリアンプ、AGC回路、A/D変換器を経て、全てデジタル信号で処理される。デジタル信号に変換された受信信号は、1データフレーム長さ毎に相関演算部10に入力される。相関演算部10は、バンドパスフィルタ(BPF)11、直交検波部12の出力(同相成分)13Aが入力されるマッチドフィルタ14A、直交検波部12の出力13B(直交成分)が入力されるマッチドフィルタ14B、マッチドフィルタ14Aの出力(同相成分の自己相関結果)16Aが入力されその2乗値を求める2乗計算部15A、マッチドフィルタ14Bの出力(直交成分の自己相関結果)16Bが入力されその2乗値を求める2乗計算部15B、2乗計算部15Aの出力と2乗計算部15Bの出力との加算値を出力する加算器を含む。   In underwater communication according to an embodiment of the present invention, QPSK modulation is used as a digital modulation method in data communication. In FIG. 1, a modulated signal received by a receiver (not shown) is all processed as a digital signal through a preamplifier, an AGC circuit, and an A / D converter (not shown). The received signal converted into the digital signal is input to the correlation calculation unit 10 for each data frame length. The correlation calculation unit 10 includes a bandpass filter (BPF) 11, a matched filter 14A to which the output (in-phase component) 13A of the quadrature detection unit 12 is input, and a matched filter to which the output 13B (quadrature component) of the quadrature detection unit 12 is input. 14B, the output of the matched filter 14A (autocorrelation result of the in-phase component) 16A is input and the square calculation unit 15A for obtaining the square value thereof, and the output of the matched filter 14B (autocorrelation result of the quadrature component) 16B is input 2 An adder that outputs an addition value of the output of the square calculation unit 15B and the output of the square calculation unit 15B and the output of the square calculation unit 15B for obtaining a square value is included.

加算器の出力、即ち、相関演算部10の出力(相関性を示す時間関数)17は、同期点検出部20に入力される。同期点検出部20では、検出されたピーク値を参照して同期点を求め、この同期点をもとに制御信号30を介してスイッチ40をONの状態にし、予め定められたデータシンボル系列長の時間を経ると制御信号30を介してスイッチ40をOFFの状態にする。スイッチ40がONの状態になると、受信信号はQPSK復調部50へ伝送され、QPSKシンボルデータが出力される。次に、相関演算部10の動作について説明する。   The output of the adder, that is, the output (time function indicating the correlation) 17 of the correlation calculation unit 10 is input to the synchronization point detection unit 20. The synchronization point detection unit 20 obtains a synchronization point with reference to the detected peak value, turns on the switch 40 via the control signal 30 based on the synchronization point, and sets a predetermined data symbol sequence length. When the time elapses, the switch 40 is turned off via the control signal 30. When the switch 40 is turned on, the received signal is transmitted to the QPSK demodulator 50, and QPSK symbol data is output. Next, the operation of the correlation calculation unit 10 will be described.

図2は、本発明の実施例におけるデータフレームの構成を説明する図である。   FIG. 2 is a diagram for explaining the structure of a data frame in the embodiment of the present invention.

図2に示すように、送信側において形成されるデータフレーム100の構成は、データフレームの先頭に位置する同期シンボル系列101と、QPSKシンボルデータ系列102とからなる。同期シンボル系列101は、QPSK信号の同相成分Iと直交成分Qのそれぞれに、同じ同期シンボルを割り当てて作られている。このようなデータフレーム100が、送波器から送信され、音響伝搬路を経て受信され、更にデジタル化された後、相関演算部10に入力される。相関演算部10に入力された受信信号は、まずバンドパスフィルタ11で帯域制限され、直交検波部12に送られる。直交検波時における搬送波の位相誤差をφとすると、直交検波部12から復調される同相成分13AはIcosφとなり、直交検波部12から復調される直交成分13BはQsinφとなる。ここでデータフレーム100の同相成分Iと直交成分Qは等しいので、I=Qである。   As shown in FIG. 2, the configuration of the data frame 100 formed on the transmission side includes a synchronization symbol sequence 101 located at the head of the data frame and a QPSK symbol data sequence 102. The synchronization symbol sequence 101 is created by assigning the same synchronization symbol to each of the in-phase component I and the quadrature component Q of the QPSK signal. Such a data frame 100 is transmitted from the transmitter, received through the acoustic propagation path, further digitized, and input to the correlation calculation unit 10. The received signal input to the correlation calculation unit 10 is first band-limited by the bandpass filter 11 and sent to the quadrature detection unit 12. If the phase error of the carrier wave at the time of quadrature detection is φ, the in-phase component 13A demodulated from the quadrature detection unit 12 is Icosφ, and the quadrature component 13B demodulated from the quadrature detection unit 12 is Qsinφ. Here, since the in-phase component I and the quadrature component Q of the data frame 100 are equal, I = Q.

次に、同相成分13Aと直交成分13Bはマッチドフィルタ14A、14Bにそれぞれ伝送され、データフレーム100で付加された同期シンボル系列101と自己相関演算が行われる。位相誤差のない場合の自己相関演算の結果を同相成分についてI’、直交成分についてQ’と書くと、位相誤差φがある場合には、自己相関演算の結果にも搬送波との位相誤差φが含まれており、同相成分の自己相関結果16AはI’cosφ、直交成分の自己相関結果16BはQ’sinφとなる。このままでは同期点の検出が困難であるが、I’=Q’であることから、2乗計算部15A、15Bにおいて、自己相関結果16A、16Bをそれぞれ2乗しその結果を加算することで、(I’cosφ)+(Q’sinφ)=(I’cosφ)+(I’sinφ)=(I’)となり、同期点で高い相関性を示す時間関数17が得られる。 Next, the in-phase component 13A and the quadrature component 13B are transmitted to the matched filters 14A and 14B, respectively, and an autocorrelation calculation is performed with the synchronization symbol sequence 101 added in the data frame 100. When the autocorrelation calculation result when there is no phase error is written as I ′ for the in-phase component and Q ′ for the quadrature component, if there is a phase error φ, the phase error φ with the carrier wave is also included in the autocorrelation calculation result. The in-phase component autocorrelation result 16A is I'cosφ, and the quadrature component autocorrelation result 16B is Q'sinφ. Although it is difficult to detect the synchronization point as it is, since I ′ = Q ′, the square calculation units 15A and 15B respectively square the autocorrelation results 16A and 16B and add the results. (I'cosφ) 2 + (Q'sinφ) 2 = (I'cosφ) 2 + (I'sinφ) 2 = (I ') 2 , and the time function indicating a high correlation in the synchronous point 17 is obtained.

同期点検出部20では、時間関数17のピーク値Pを検出し、ピーク値Pの位置を同期点として、QPSK復調部50への入力がQPSKシンボルデータ系列102の先頭からになるように、スイッチ40をONの状態にする。なお、確実な同期には、符号サイドローブの小さな同期シンボル系列が望ましいので、同期シンボル系列101としては、バーカー符号やM系列符号を適用するのが好ましい。また、量子化誤差によって、同期点検出部20において求めた同期点においても搬送波の残留位相誤差が残る可能性があるため、QPSK復調部50の出力結果は図示しない適応等化器によって更に処理されることが望ましい。なお、適応等化器は、等化精度のよい判定帰還形適応等化器とするのが望ましい。以上のような構成によって、位相同期ループやデジタル位相同期ループを必要としない簡単な構成で、且つ確実なフレーム同期をとることができ、安定した通信を行うことができる。
(実施例2)
図3は、本発明の水中通信方式の他の実施例における受信側での信号処理を説明するブロック図である。
The synchronization point detector 20 detects the peak value P of the time function 17 and switches the position of the peak value P to the synchronization point so that the input to the QPSK demodulator 50 comes from the beginning of the QPSK symbol data sequence 102. 40 is turned on. Note that a synchronization symbol sequence having a small code side lobe is desirable for reliable synchronization, and therefore, it is preferable to apply a Barker code or an M sequence code as the synchronization symbol sequence 101. Further, since the residual error of the carrier wave may remain at the synchronization point obtained by the synchronization point detector 20 due to the quantization error, the output result of the QPSK demodulator 50 is further processed by an adaptive equalizer (not shown). It is desirable. Note that the adaptive equalizer is preferably a decision feedback adaptive equalizer with high equalization accuracy. With the above configuration, a simple configuration that does not require a phase-locked loop or a digital phase-locked loop, reliable frame synchronization can be achieved, and stable communication can be performed.
(Example 2)
FIG. 3 is a block diagram for explaining signal processing on the receiving side in another embodiment of the underwater communication system of the present invention.

図3において、図示しない受波器はアレイとなっており、異なるアレイ開口部から、Nチャンネル(Nは2以上の自然数)の受信信号が得られている。各チャンネルの受信信号は、1データフレーム長さ毎に各チャンネルに設けられた相関演算部10でそれぞれ処理され、各相関演算部10の出力は、同期点検出部21に入力される。同期点検出部21では、各相関演算部10からの情報に基づいて同期点が求められ、制御信号31を介して、各チャンネルに設けられたスイッチ40を全て同時にON、OFF制御する。スイッチ40によって通過した各チャンネルの受信信号は、各チャンネルに設けられたQPSK復調部50によって復調され、Nチャンネルの入力端を有するマルチチャンネル適応等化器60へ伝送され、マルチチャンネル適応等化器60によって、搬送波の残留位相誤差やマルチパス音響伝搬路のフェージング特性が補償され、誤りのない正しいQPSKシンボルデータが出力される。なお、マルチチャンネル適応等化器は、等化精度のよいマルチチャンネル判定帰還形適応等化器が望ましい。次に、同期点検出部21の同期点決定方法を図4、図5を用いて説明する。   In FIG. 3, the receiver not shown is an array, and N channel (N is a natural number of 2 or more) received signals are obtained from different array openings. The reception signal of each channel is processed by the correlation calculation unit 10 provided for each channel for each data frame length, and the output of each correlation calculation unit 10 is input to the synchronization point detection unit 21. The synchronization point detection unit 21 obtains a synchronization point based on information from each correlation calculation unit 10, and simultaneously controls ON / OFF of all the switches 40 provided in each channel via the control signal 31. The received signal of each channel passed by the switch 40 is demodulated by a QPSK demodulator 50 provided in each channel, and transmitted to a multi-channel adaptive equalizer 60 having an N-channel input terminal, and the multi-channel adaptive equalizer 60 compensates for the residual phase error of the carrier wave and the fading characteristics of the multipath acoustic propagation path, and outputs correct QPSK symbol data without error. The multi-channel adaptive equalizer is preferably a multi-channel decision feedback adaptive equalizer with good equalization accuracy. Next, the synchronization point determination method of the synchronization point detection unit 21 will be described with reference to FIGS.

図4は、本発明の他の実施例における同期点検出部の同期点決定方法を示すフローチャート図である。   FIG. 4 is a flowchart showing a synchronization point determination method of the synchronization point detection unit in another embodiment of the present invention.

図3に示す同期点検出部21には、各チャンネルに設けられた相関演算部10によって得られる時間関数s〜sが送られてくる。同期点検出部21では、まず時間関数s〜sのそれぞれについてピーク値P〜Pを求める。P=max{s}〜P=max{s}である。次に、ピーク値P〜Pの中から最も大きいピーク値PK1=max{P:k=1〜N}を求め、そのピーク値PK1をとるチャンネルK1の時間関数sK1を抽出する。次に、時間関数sK1のピーク値PK1の位置を同期点s K1=argmax{sK1}として求める。送信間隔Lの送信データフレームが連続して送られており、2番目以降のデータフレームを受信している場合には、直前のデータフレームの同期点sn−1 K1と、今求めた同期点s K1との時間間隔ΔT=s K1−sn−1 K1を求め、送信間隔Lと比較する。時間間隔ΔTと送信間隔Lとの差の絶対値│L−ΔT│が、予め定められた値X1以内にある時、即ち、│L−ΔT│≦X1である時は同期点s K1を採用し、│L−ΔT│≦X1でない時は、L+sn−1 K1を同期点として採用する。また、│L−ΔT│≦X1でない時は、s K1=L+sn−1 K1を同期点として採用するか、図示しない再送要求処理部へ再送処理開始の信号を送る。なお、予め定められた値X1は、シンボルレートTと、マルチチャンネル適応等化器60で等化の対象とされるシンボル数Jとの積T・J以下の値とすればよい。このような方式で同期点を求めることによって、時間的、空間的に激しく変化する浅海域マルチパス・フェージング環境下においても、正しくフレーム同期がとれる確率が向上し、安定した水中通信が可能となる。 The synchronization point detection unit 21 shown in FIG. 3 receives time functions s 1 to s N obtained by the correlation calculation unit 10 provided for each channel. The synchronization point detection unit 21 first obtains peak values P 1 to P N for each of the time functions s 1 to s N. P 1 = max {s 1 } to P N = max {s N }. Then, the peak value P 1 to P largest peak value among the N P K1 = max: seek {P K k = 1~N}, extracting the time function s K1 channel K1 taking the peak value P K1 To do. Next, the position of the peak value P K1 of the time function s K1 is obtained as a synchronization point s n K1 = argmax {s K1 }. When transmission data frames of the transmission interval L are continuously transmitted and the second and subsequent data frames are received, the synchronization point s n-1 K1 of the immediately preceding data frame and the synchronization point obtained now seeking s n time interval between K1 ΔT 1 = s n K1 -s n-1 K1, compared with the transmission interval L. When the absolute value | L−ΔT 1 | of the difference between the time interval ΔT 1 and the transmission interval L is within a predetermined value X1, that is, when | L−ΔT 1 | ≦ X1, the synchronization point s When n K1 is employed and | L−ΔT 1 | ≦ X1, L + s n−1 K1 is employed as the synchronization point. When | L−ΔT 1 | ≦ X1 is not satisfied, s n K1 = L + s n−1 K1 is adopted as a synchronization point, or a retransmission processing start signal is sent to a retransmission request processing unit (not shown). The predetermined value X1 may be a value equal to or less than the product T · J of the symbol rate T and the number J of symbols to be equalized by the multichannel adaptive equalizer 60. By finding the synchronization point in this way, the probability of correct frame synchronization is improved and stable underwater communication is possible even in shallow sea multipath fading environments that change drastically in time and space. .

図5は、本発明の他の実施例における同期点検出部の他の同期点決定方法を示すフローチャート図である。   FIG. 5 is a flowchart showing another synchronization point determination method according to another embodiment of the present invention.

図3に示す同期点検出部21には、各チャンネルに設けられた相関演算部10によって得られる時間関数s〜sが送られてくる。同期点検出部21では、まず、各チャンネルの時間関数s〜s毎のピーク値P〜Pをそれぞれ抽出する。P=max{s}〜P=max{s}である。次に、ピーク値P〜Pの位置を含む、予め定められた範囲の時間関数s〜sの平均値M〜Mの値を、各チャンネルについてそれぞれ計算する。M=mean{s}〜M=mean{s}である。求めたピーク値P〜P及び平均値M〜Mを用い、各チャンネルについて、これらの比P/M〜P/Mを求め、比P/M〜P/Mの中から最も大きい比PK2/MK2=max{P/M:k=1〜N}をとるチャンネルK2の時間関数sK2を抽出して、時間関数sK2のピーク値PK2の位置を同期点s K2=argmax{sK2}として求める。送信間隔Lの送信データフレームが連続して送られており、2番目以降のデータフレームを受信している場合には、直前のデータフレームの同期点sn−1 2Kと、今求めた同期点s K2との時間間隔ΔT=s K2−sn−1 K2を求め、送信間隔Lと比較する。時間間隔ΔTと送信間隔Lとの差の絶対値│L−ΔT│が、予め定められた値X2以内にある時、即ち、│L−ΔT│≦X2である時は同期点s K2を採用する。また、│L−ΔT│≦X2でない時は、s K2=L+sn−1 K2を同期点として採用するか、図示しない再送要求処理部へ再送処理開始の信号を送る。なお、同期シンボル系列101には、バーカー符号やM系列符号を適用する。また、予め定められた範囲X2は、同期シンボル系列101のシンボル長をLsとして、ピーク値P〜Pの各位値を中心に±Ls以上の範囲をとればよく、予め定められた値X2は、シンボルレートTと、マルチチャンネル適応等化器60で等化の対象とされるシンボル数Jとの積T・J以下の値とすればよい。以下、図5のフローチャートに基づく同期決定方式について、水深5mの無響水槽内で行った実験結果を、図6〜図8を用いて説明する。 The synchronization point detection unit 21 shown in FIG. 3 receives time functions s 1 to s N obtained by the correlation calculation unit 10 provided for each channel. In the synchronous point detecting unit 21 first extracts the peak value P 1 to P N for every time function s 1 ~s N for each channel, respectively. P 1 = max {s 1 } to P N = max {s N }. Then, including the location of the peak value P 1 to P N, the value of the average value M 1 ~M N time function s 1 ~s N of the predetermined range, calculating respectively for each channel. M 1 = mean {s 1 } to M N = mean {s N }. Using a peak value P 1 to P N and the average value M 1 ~M N obtained, for each channel, obtains these ratios P 1 / M 1 ~P N / M N, the ratio P 1 / M 1 ~P N / M N extracts the time function s K2 of the channel K2 taking the largest ratio P K2 / M K2 = max {P K / M K : k = 1 to N}, and the peak value of the time function s K2 The position of P K2 is obtained as a synchronization point s n K2 = argmax {s K2 }. When transmission data frames of the transmission interval L are continuously transmitted and the second and subsequent data frames are received, the synchronization point s n-1 2K of the immediately previous data frame and the synchronization point obtained now seeking s n time interval between K2 ΔT 2 = s n K2 -s n-1 K2, compared with the transmission interval L. When the absolute value | L−ΔT 2 | of the difference between the time interval ΔT 2 and the transmission interval L is within a predetermined value X2, that is, when | L−ΔT 2 | ≦ X2, the synchronization point s n K2 is adopted. When | L−ΔT 2 | ≦ X2, s n K2 = L + s n−1 K2 is adopted as a synchronization point, or a retransmission processing start signal is sent to a retransmission request processing unit (not shown). Note that Barker code or M-sequence code is applied to the synchronization symbol sequence 101. Further, the predetermined range X2 may be a range of ± Ls or more around each peak value P 1 to P N with the symbol length of the synchronization symbol sequence 101 as Ls, and a predetermined value X2 May be a value equal to or less than the product T · J of the symbol rate T and the number J of symbols to be equalized by the multi-channel adaptive equalizer 60. Hereinafter, the results of experiments performed in an anechoic water tank having a water depth of 5 m with respect to the synchronization determination method based on the flowchart of FIG. 5 will be described with reference to FIGS.

図6は、本発明の実施例における実験用データフレームの構成を示す図であり、無響水槽内の水深4mの位置に置かれた、無指向性の送波用ハイドロホンから送信した、実験用データフレームの構成図である。   FIG. 6 is a diagram showing the structure of an experimental data frame in the embodiment of the present invention, and an experiment transmitted from an omnidirectional wave transmitting hydrophone placed at a depth of 4 m in an anechoic water tank. It is a block diagram of a data frame for use.

図6に示すように、実験用データフレーム110は、2つの連続したデータフレーム111を含むように構成した。データフレーム111は、図2に示したデータフレーム100と同じ構成となっており、13シンボル長の13ビットバーカー符号系列112と、1154シンボル長のデータシンボル系列113とから構成されており、連続するデータフレーム111の間には、3シンボル長のポーズを設けた。なお、搬送波の周波数は70kHz、シンボルレートは30000ボーとした。   As shown in FIG. 6, the experimental data frame 110 was configured to include two consecutive data frames 111. The data frame 111 has the same configuration as that of the data frame 100 shown in FIG. 2, and is composed of a 13-symbol 13-bit Barker code sequence 112 and a 1154-symbol data symbol sequence 113, and is continuous. A pause of 3 symbol length was provided between the data frames 111. The frequency of the carrier wave was 70 kHz and the symbol rate was 30000 baud.

実験用データフレーム110を受信するための受波器は、5つの無指向性のハイドロホンからなる、CH1〜CH5までの5チャンネルの垂直アレイとし、各チャンネル間の間隔は搬送波の2波長分43mmとした。更に、垂直アレイをそのアレイ中心が水深75cm、送波用ハイドロホンとの水平距離5.8mの位置に垂直に配置した。   The receiver for receiving the experimental data frame 110 is a vertical array of 5 channels from CH1 to CH5 consisting of 5 omnidirectional hydrophones, and the interval between each channel is 43 mm for 2 wavelengths of the carrier wave. It was. Further, the vertical array was arranged vertically at the center of the array at a depth of 75 cm and a horizontal distance of 5.8 m from the hydrophone for transmission.

図7は、本発明の実施例の相関演算部により得られる時間関数の実験結果を示す特性図であり、垂直アレイの各チャンネル(水深の浅い方からCH1〜CH5とする)の受信信号波形を相関演算部10でそれぞれ処理することによって得た時間関数sCH1〜sCH5を示す。 FIG. 7 is a characteristic diagram showing the experimental results of the time function obtained by the correlation calculation unit of the embodiment of the present invention, and shows the received signal waveform of each channel of the vertical array (CH1 to CH5 from the shallower water depth). The time functions s CH1 to s CH5 obtained by processing in the correlation calculation unit 10 are shown.

図8は、本発明の実施例における同期点比較の実験結果を示す特性表を示す図であり、2番目に送信されたデータフレーム111に着目して、比P/M〜P/Mを時間関数sCH1〜sCH5から計算した同期点比較結果である。なお、予め定められた範囲X2は±100シンボル長とした。 FIG. 8 is a diagram showing a characteristic table showing the result of an experiment of synchronization point comparison in the embodiment of the present invention, and paying attention to the data frame 111 transmitted second, the ratios P 1 / M 1 to P 5 / M 5 is a synchronous point comparison result of calculating the time function s CH1 ~s CH5. The predetermined range X2 is set to ± 100 symbol length.

受信チャンネル数が1つの場合、特に水面からのマルチパス成分によって搬送波成分が干渉し、不感領域に入ってしまうと、時間関数sCH3のような演算結果が得られる可能性が出てくる。この場合、時間関数sCH3のピーク値における時間を同期点としてしまうと、ピーク値を取る位置は同期シンボル系列側ではなくデータシンボル系列側にあり、同期が取れなくなってしまう。一方、水中における70kHzの搬送波の波長は約21.5mmであり、非常に短いので、受信位置を少しずらすことによって不感領域を避けることができるので、複数の受信アレイの使用により、確実に同期がとれる。 When the number of reception channels is one, especially when the carrier component interferes with the multipath component from the water surface and enters the insensitive region, there is a possibility that an operation result such as the time function s CH3 is obtained. In this case, if the time at the peak value of the time function s CH3 is set as the synchronization point, the position where the peak value is taken is on the data symbol sequence side, not on the synchronization symbol sequence side, and synchronization is not achieved. On the other hand, the wavelength of a carrier wave of 70 kHz in water is about 21.5 mm, which is very short, so that a dead area can be avoided by slightly shifting the reception position, so that the use of multiple reception arrays ensures synchronization. I can take it.

また、時間関数sCH1〜sCH5のピーク値のみの比較では、CH1とCH3でほとんど変わらないものの、図8に示す、比P/M〜P/Mで比較することによって、符号サイドローブが小さく、鋭い自己相関関数を有するバーカー符号の特徴も利用でき、同期を取る基準とすべきチャンネルの決定が確実になされる。 Further, in the comparison of only the peak values of the time functions s CH1 to s CH5 , although there is almost no change between CH1 and CH3, by comparing with the ratios P 1 / M 1 to P 5 / M 5 shown in FIG. The characteristics of Barker codes with small side lobes and sharp autocorrelation functions can also be used to ensure the determination of the channel to be synchronized.

このように、マルチパス環境下においても、本同期点検出方式によって確実なフレーム同期が可能となることが実験的に確認された。   Thus, it was experimentally confirmed that reliable frame synchronization can be achieved by this synchronization point detection method even in a multipath environment.

以上のような構成により同期点を求めることによって、搬送波が不感領域に入ってしまう可能性を排除し、バーカー符号やM系列符号の自己相関関数の相関値だけでなく、低い符号サイドローブを有するという特徴も利用して同期点を求められる。従って、時間的、空間的に激しく変化する浅海域マルチパス・フェージング環境下においても、正しくフレーム同期がとれる確率が向上し、適応等化器によって搬送波の残留位相誤差も排除されるため、安定した水中通信が可能となる。   By obtaining the synchronization point with the above configuration, the possibility that the carrier wave enters the insensitive area is eliminated, and not only the correlation value of the autocorrelation function of the Barker code or the M-sequence code, but also a low code side lobe. Using this feature, the synchronization point can be obtained. Therefore, even in shallow sea multipath fading environments that change drastically in time and space, the probability of correct frame synchronization is improved, and the residual phase error of the carrier wave is eliminated by the adaptive equalizer. Underwater communication is possible.

QPSK変調を用いた本発明水中通信システムを纏めると次の通りである。   The underwater communication system of the present invention using QPSK modulation is summarized as follows.

本発明の水中通信システムでは、送信側では、QPSK送信データフレームの先頭の、QPSK信号の同相成分と直交成分のそれぞれに同一の同期シンボルを付加して送信し、受信側では、受信信号を直交検波し、直交検波後の同相成分と同期シンボルとの自己相関関数fと、直交検波後の直交成分と同期シンボルとの自己相関関数gとを求め、自己相関関数fの2乗と自己相関関数gの2乗とを加算して得られる時間関数sのピーク値Pの位置を同期点として検出して、フレーム同期をとる。同期シンボルには、バーカー符号又はM系列符号を用いる。   In the underwater communication system of the present invention, on the transmission side, the same synchronization symbol is added to each of the in-phase component and the quadrature component of the QPSK signal at the beginning of the QPSK transmission data frame, and transmission is performed on the reception side. The autocorrelation function f between the in-phase component and the synchronization symbol after detection and quadrature detection and the autocorrelation function g between the quadrature component and the synchronization symbol after quadrature detection are obtained, and the square of the autocorrelation function f and the autocorrelation function The position of the peak value P of the time function s obtained by adding the square of g is detected as a synchronization point, and frame synchronization is established. Barker code or M-sequence code is used for the synchronization symbol.

本発明の水中通信システムでは、受信側に複数の受信チャンネルを具備しており、各チャンネルにおいてピーク値Pを検出し、各チャンネルのピーク値Pが最も大きいチャンネルにおける時間関数sのピーク値Pの位置を同期点として検出して、又は、時間関数sのピーク値Pと、ピーク値Pの時間を含む一定時間範囲内における時間関数sの平均値Mとの比P/Mを各チャンネルについて求め、比P/Mの値が最も大きいチャンネルにおける時間関数sのピーク値Pの位置を同期点として検出して、同期点で全チャンネルのフレーム同期をとる。   In the underwater communication system of the present invention, the receiving side has a plurality of receiving channels, the peak value P is detected in each channel, and the peak value P of the time function s in the channel having the largest peak value P of each channel is detected. The position P is detected as a synchronization point, or the ratio P / M between the peak value P of the time function s and the average value M of the time function s within a certain time range including the time of the peak value P is obtained for each channel. The position of the peak value P of the time function s in the channel with the largest ratio P / M value is detected as a synchronization point, and frame synchronization of all channels is established at the synchronization point.

本発明の水中通信システムでは、連続する2つの送信データフレームにおける同期点の時間間隔と、送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、該当のデータフレームについて送信側に再送要求する、又は、一つ目の送信データフレームの同期点と送信データフレーム長から、二つ目の送信データフレームの同期点を決定してフレーム同期をとる。   In the underwater communication system of the present invention, when the absolute value of the difference between the synchronization point time interval between two consecutive transmission data frames and the transmission interval of the transmission data frame is larger than a predetermined value, the corresponding data The frame is synchronized by requesting the transmission side to retransmit the frame, or by determining the synchronization point of the second transmission data frame from the synchronization point of the first transmission data frame and the transmission data frame length.

本発明の水中通信システムでは、受信側に適応等化器を備え、適応等化器は1チャンネル又はマルチチャンネルの判定帰還形適応等化器としている。   In the underwater communication system of the present invention, an adaptive equalizer is provided on the receiving side, and the adaptive equalizer is a one-channel or multi-channel decision feedback type adaptive equalizer.

本発明の水中通信システムの代表的な効果は次の通りである。(1)位相同期ループやデジタル位相同期ループを必要としない簡単な構成で、同期シンボルの鋭い自己相関特性による確実な同期が可能となり、安定した通信を行なうことが可能となる。(2)激しく変化するマルチパス・フェージング環境化においても、正しい同期点が得られる確率が向上し、正確なデータの送受信が可能となる。(3)搬送波の残留位相誤差の影響を適応的に精度よく除去できるので、安定して正確な水中通信が可能となる。   Typical effects of the underwater communication system of the present invention are as follows. (1) With a simple configuration that does not require a phase-locked loop or a digital phase-locked loop, reliable synchronization can be achieved by the sharp autocorrelation characteristics of the synchronization symbols, and stable communication can be performed. (2) Even in a rapidly changing multipath fading environment, the probability that a correct synchronization point is obtained is improved, and accurate data transmission / reception becomes possible. (3) Since the influence of the residual phase error of the carrier wave can be adaptively and accurately removed, stable and accurate underwater communication is possible.

本発明の水中通信システムの構成の特徴は以下の通りである。   The characteristics of the configuration of the underwater communication system of the present invention are as follows.

(1)デジタル変調方式としてQPSK変調を用いる水中通信システムであって、受波器、前記受波器で受信された変調信号をデジタル信号に変換するA/D変換器、及び、前記デジタル信号が1データフレーム長さ毎に入力される相関演算部を含む複数の受信チャンネルと、前記各チャンネルの前記相関演算部の出力信号が入力される同期点検出部と、前記チャンネルのそれぞれに対応して設けられ前記チャンネルのそれぞれの前記A/D変換器の出力信号が入力されるスイッチ手段と、前記スイッチ手段のそれぞれに接続されるQPSK復調部と、前記各QPSK復調部の出力信号が入力されQPSKシンボルデータが出力されるマルチチャンネル適応等化器とを有し、前記同期点検出部において前記各相関演算部の出力信号に基づいて同期点が求められ、前記同期点検出部からの制御信号を介して、前記各チャンネルの前記スイッチ手段全てが同時にON、OFF制御され、前記各スイッチ手段を通過した前記各チャンネルの受信信号は、前記各チャンネルの前記QPSK復調部により復調され、前記マルチチャンネル適応等化器へ伝送され、前記マルチチャンネル適応等化器により、搬送波の残留位相誤差やマルチパス音響伝搬路のフェージング特性が補償され、QPSKシンボルデータが出力され、前記相関演算部のそれぞれは、前記相関演算部10に入力される受信信号を帯域制限するバンドパスフィルタと、前記バンドパスフィルタの出力信号を直交検波し、同相成分及び直交成分を出力する直交検波部と、前記同相成分が入力され前記同相成分の自己相関結果を出力する第1のマッチドフィルタと、前記直交成分が入力され前記直交成分の自己相関結果を出力する第2のマッチドフィルタと、前記第1のマッチドフィルタの出力の2乗値を求める第1の2乗計算部と、前記第2のマッチドフィルタの出力の2乗値を求める第2の2乗計算部と、前記第1の2乗計算部の出力信と前記第2の2乗計算部の出力信との加算値を出力する加算器を含み、前記加算器の出力を前記相関演算部のそれぞれ出力信号とすることを特徴とする水中通信システム。   (1) An underwater communication system using QPSK modulation as a digital modulation method, wherein a receiver, an A / D converter that converts a modulation signal received by the receiver into a digital signal, and the digital signal Corresponding to each of a plurality of reception channels including a correlation calculation unit input every data frame length, a synchronization point detection unit to which an output signal of the correlation calculation unit of each channel is input, Switch means provided for receiving an output signal of the A / D converter for each of the channels; a QPSK demodulator connected to each of the switch means; and an output signal of each QPSK demodulator for receiving QPSK A multi-channel adaptive equalizer that outputs symbol data, and based on the output signal of each correlation calculation unit in the synchronization point detection unit Through the control signal from the synchronization point detection unit, all the switch means of each channel are simultaneously ON and OFF controlled, and the received signal of each channel that has passed through each switch means, Demodulated by the QPSK demodulator for each channel and transmitted to the multi-channel adaptive equalizer, and the multi-channel adaptive equalizer compensates for residual phase error of the carrier wave and fading characteristics of the multipath acoustic propagation path, QPSK symbol data is output, and each of the correlation calculation units quadrature-detects the band-pass filter for band-limiting the received signal input to the correlation calculation unit 10 and the output signal of the band-pass filter, A quadrature detection unit that outputs a quadrature component and the in-phase component are input and an autocorrelation result of the in-phase component is output. A first matched filter, a second matched filter that receives the quadrature component and outputs an autocorrelation result of the quadrature component, and a first square that obtains a square value of the output of the first matched filter A calculation unit; a second square calculation unit for obtaining a square value of the output of the second matched filter; an output signal of the first square calculation unit; and an output signal of the second square calculation unit. An underwater communication system comprising: an adder that outputs an addition value of the correlation calculation unit, wherein the output of the adder is an output signal of each of the correlation calculation units.

QPSK変調を用いた水中通信システムの送信側では、QPSK送信データフレームの先頭の、QPSK信号の同相成分と直交成分のそれぞれに同一の同期シンボルを付加して送信する。QPSK変調を用いた水中通信システムの受信側では、受信信号を直交検波し、前記直交検波後の同相成分と前記同期シンボルとの自己相関関数fと、前記直交検波後の直交成分と前記同期シンボルとの自己相関関数gとを求め、前記自己相関関数fの2乗と前記自己相関関数gの2乗とを加算して得られる時間関数sのピーク値Pの位置を同期点として検出して、フレーム同期をとる。   On the transmission side of the underwater communication system using QPSK modulation, the same synchronization symbol is added to each of the in-phase component and the quadrature component of the QPSK signal at the beginning of the QPSK transmission data frame and transmitted. On the receiving side of the underwater communication system using QPSK modulation, the received signal is quadrature-detected, the autocorrelation function f between the in-phase component after the quadrature detection and the synchronization symbol, the quadrature component after the quadrature detection and the synchronization symbol And the position of the peak value P of the time function s obtained by adding the square of the autocorrelation function f and the square of the autocorrelation function g is detected as a synchronization point. Take frame synchronization.

前記同期点検出部には、前記各チャンネルの相関演算部10の出力信号であり、相関性を示す時間関数が送られ、同期点検出部では、前記各チャンネル毎の前記時間関数のピーク値が抽出され、前記ピーク値の位置を含む予め定められた範囲の時間関数の平均値の値が、前記各チャンネルについてそれぞれ計算され、求められた前記ピーク値及び前記平均値を用いて、前記各チャンネルについて、前記平均値に対する前記ピーク値の比が求められ、前記比が最大である前記チャンネルの前記時間関数を抽出して、前記比が最大である前記チャンネルの前記時間関数の前記ピーク値の位置を同期点s K2として求める。送信間隔Lの送信データフレームが連続して送られており、2番目以降のデータフレームを受信している場合には、直前のデータフレームの前記同期点sn−1 K2と、今回求めた前記同期点s K2との時間間隔ΔTが求められ、前記時間間隔ΔTが前記送信間隔Lと比較され、│L−ΔT│≦X2(X2:予め定められた値)である時は、同期点を今回求めた同期点s K2とし、│L−ΔT│≦X2でない時は、s K2=L+sn−1 K2を同期点とするか、又は、再送要求処理部へ再送処理開始の信号を送る。値X2は、同期シンボルのシンボル長をLsとして、前記ピーク値の各位値を中心に±L以上の範囲とし、シンボルレートTと、マルチチャンネル適応等化器で等化の対象とされるシンボル数Jとの積T・J以下の値とする。 The sync point detection unit is an output signal of the correlation calculation unit 10 of each channel, and a time function indicating the correlation is sent. In the sync point detection unit, the peak value of the time function for each channel is obtained. The average value of the time function in a predetermined range including the position of the peak value that is extracted is calculated for each channel, and the peak value and the average value are used to calculate each channel. A ratio of the peak value to the average value is determined, the time function of the channel having the maximum ratio is extracted, and the position of the peak value of the time function of the channel having the maximum ratio is extracted. Is obtained as a synchronization point s n K2 . When transmission data frames having a transmission interval L are continuously transmitted and the second and subsequent data frames are received, the synchronization point s n-1 K2 of the immediately preceding data frame and the current time obtained are obtained. time interval [Delta] T 2 between the synchronization point s n K2 is determined, the time interval [Delta] T 2 is compared with the transmission interval L, │L-ΔT 2 │ ≦ X2: when in (X2 predetermined value) is The synchronization point is the synchronization point s n K2 obtained this time, and when | L−ΔT 2 | ≦ X2, s n K2 = L + s n−1 K2 is set as the synchronization point or retransmitted to the retransmission request processing unit. Send a signal to start processing. The value X2 is a symbol to be equalized by a symbol rate T and a multi-channel adaptive equalizer, with the symbol length of the synchronization symbol being Ls, a range of ± L s or more around each peak value. It is set to a value equal to or less than the product T · J with the number J.

前記同期点検出部では、検出された前記ピーク値を参照して同期点を求め、この同期点をもとに制御信号を介して前記スイッチ手段をONの状態にし、予め定められたデータシンボルの系列長の時間を経過すると前記制御信号を介して前記スイッチ手段をOFFの状態にする。前記スイッチ手段がONの状態になると、受信信号はQPSK復調部へ伝送され、QPSKシンボルデータが出力される。   The synchronization point detection unit obtains a synchronization point with reference to the detected peak value, turns on the switch means via a control signal based on the synchronization point, and determines a predetermined data symbol. When the sequence length time elapses, the switch means is turned off via the control signal. When the switch means is turned on, the received signal is transmitted to the QPSK demodulator and QPSK symbol data is output.

(2)(1)の水中通信システムにおいて、前記同期シンボルは、バーカー符号又はM系列符号である。   (2) In the underwater communication system of (1), the synchronization symbol is a Barker code or an M-sequence code.

(3)(2)の水中通信システムにおいて、複数の前記受信チャンネルのそれぞれにおいて前記ピーク値を検出し、前記各チャンネルの前記ピーク値が最も大きいチャンネルにおける前記時間関数の、前記ピーク値の位置を同期点として検出して、前記同期点で全チャンネルのフレーム同期をとる。   (3) In the underwater communication system of (2), the peak value is detected in each of the plurality of reception channels, and the position of the peak value of the time function in the channel having the largest peak value of each channel is determined. Detection is performed as a synchronization point, and frame synchronization of all channels is performed at the synchronization point.

(4)(2)の水中通信システムにおいて、複数の前記受信チャンネルのそれぞれにおいて前記ピーク値を検出し、前記各チャンネルの前記時間関数のピーク値Pと、前記ピーク値の時間を含む一定時間範囲内における前記時間関数の平均値Mとの比P/Mを前記各チャンネルについて求め、前記比P/Mの値が最も大きい前記チャンネルにおける前記時間関数の、前記ピーク値Pの位置を同期点として検出して、前記同期点で全チャンネルのフレーム同期をとる。   (4) In the underwater communication system according to (2), the peak value is detected in each of the plurality of reception channels, and a certain time range including the peak value P of the time function of each channel and the time of the peak value A ratio P / M with respect to the average value M of the time function is obtained for each channel, and the position of the peak value P of the time function in the channel having the largest value of the ratio P / M is used as a synchronization point. Detection is performed, and frame synchronization of all channels is performed at the synchronization point.

(5)(2)の水中通信システムにおいて、連続する2つの前記送信データフレームにおける前記同期点の時間間隔と、前記送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、該当のデータフレームについて送信側に再送要求する。   (5) In the underwater communication system of (2), the absolute value of the difference between the time interval of the synchronization point in the two consecutive transmission data frames and the transmission interval of the transmission data frame is determined from a predetermined value. If it is larger, a retransmission request is made to the transmitting side for the corresponding data frame.

(6)(2)の水中通信システムにおいて、連続する2つの前記送信データフレームにおける前記同期点の時間間隔と、前記送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、一つ目の前記送信データフレームの同期点と前記送信間隔から、二つ目の前記送信データフレームの同期点を決定してフレーム同期をとる。   (6) In the underwater communication system of (2), the absolute value of the difference between the time interval of the synchronization point and the transmission interval of the transmission data frame in two consecutive transmission data frames is determined from a predetermined value. If larger, the synchronization point of the second transmission data frame is determined from the synchronization point of the first transmission data frame and the transmission interval, and frame synchronization is established.

マルチパス・フェージング環境でも、簡単な構成の同期方式により、安定した同期を実現し、安定した通信が可能な水中通信方式及び水中通信システムを提供できる。   Even in a multipath fading environment, it is possible to provide an underwater communication system and an underwater communication system that can achieve stable synchronization by using a synchronization method with a simple configuration.

本発明の水中通信システムの実施例における受信側での信号処理を説明するブロック図。The block diagram explaining the signal processing by the side of a receiver in the Example of the underwater communication system of this invention. 本発明の実施例におけるデータフレームの構成を説明する図。The figure explaining the structure of the data frame in the Example of this invention. 本発明の水中通信システムの他の実施例における受信側での信号処理を説明するブロック図。The block diagram explaining the signal processing by the receiving side in the other Example of the underwater communication system of this invention. 本発明の他の実施例における同期点検出部の同期点決定方法を示すフローチャート図。The flowchart figure which shows the synchronization point determination method of the synchronization point detection part in the other Example of this invention. 本発明の他の実施例における同期点検出部の他の同期点決定方法を示すフローチャート図。The flowchart figure which shows the other synchronization point determination method of the synchronization point detection part in the other Example of this invention. 本発明の実施例における実験用データフレームの構成を示す図。The figure which shows the structure of the data frame for experiment in the Example of this invention. 本発明の実施例の相関演算部により得られる時間関数の実験結果を示す特性図。The characteristic view which shows the experimental result of the time function obtained by the correlation calculating part of the Example of this invention. 本発明の実施例における同期点比較の実験結果を示す特性表を示す図。The figure which shows the characteristic table | surface which shows the experimental result of the synchronous point comparison in the Example of this invention. 従来技術における水中通信システムでの超音波の伝搬を説明する図。The figure explaining propagation of the ultrasonic wave in the underwater communication system in a prior art. 従来技術における水中通信システムでの実験用送波音圧の振幅スペクトル例を示す図。The figure which shows the amplitude spectrum example of the experimental transmission sound pressure in the underwater communication system in a prior art. 従来技術における水中通信システムでの受信音圧の振幅スペクトル例を示す図。The figure which shows the amplitude spectrum example of the received sound pressure in the underwater communication system in a prior art.

符号の説明Explanation of symbols

10…相関演算部、11…バンドパスフィルタ、12…直交検波部、13A…直交検波後の同相成分、13B…直交検波後の直交成分、14A、14B…マッチドフィルタ、15A、15B…2乗計算部、16A…同相成分の自己相関結果、16B…直交成分の自己相関結果、17…相関性を示す時間関数、20、21…同期点検出部、30、31…制御信号、40…スイッチ、50…QPSK復調部、60…マルチチャンネル適応等化器、70…水中局、71…送波器、80…水上局、81…受波器、82…無線アンテナ、100、111…データフレーム、101…同期シンボル系列、102…QPSKシンボルデータ系列、110…実験用データフレーム、112…バーカー符号系列、113…データシンボル系列。
DESCRIPTION OF SYMBOLS 10 ... Correlation calculation part, 11 ... Band pass filter, 12 ... Quadrature detection part, 13A ... In-phase component after quadrature detection, 13B ... Quadrature component after quadrature detection, 14A, 14B ... Matched filter, 15A, 15B ... Square calculation 16A: In-phase component autocorrelation result, 16B: Quadrature component autocorrelation result, 17 ... Time function indicating correlation, 20, 21 ... Synchronization point detection unit, 30, 31 ... Control signal, 40 ... Switch, 50 QPSK demodulator, 60 Multichannel adaptive equalizer, 70 Underwater station, 71 Transmitter, 80 Water station, 81 Receiver, 82 Radio antenna, 100, 111 Data frame, 101 Synchronization symbol sequence, 102... QPSK symbol data sequence, 110 .. Experimental data frame, 112 .. Barker code sequence, 113.

Claims (8)

QPSK変調を用いた水中通信方式であって、送信側では、QPSK送信データフレームの先頭の、QPSK信号の同相成分と直交成分のそれぞれに同一の同期シンボルを付加して送信し、受信側では、受信信号を直交検波し、前記直交検波後の同相成分と前記同期シンボルとの自己相関関数fと、前記直交検波後の直交成分と前記同期シンボルとの自己相関関数gとを求め、前記自己相関関数fの2乗と前記自己相関関数gの2乗とを加算して得られる時間関数sのピーク値Pの位置を同期点として検出して、フレーム同期をとることを特徴とする水中通信方式。   It is an underwater communication system using QPSK modulation, and on the transmission side, the same synchronization symbol is added to each of the in-phase component and the quadrature component of the QPSK signal at the beginning of the QPSK transmission data frame and transmitted, and on the reception side, The received signal is subjected to quadrature detection, an autocorrelation function f between the in-phase component after the quadrature detection and the synchronization symbol, an autocorrelation function g between the quadrature component after the quadrature detection and the synchronization symbol, and an autocorrelation function g are obtained. An underwater communication system characterized by detecting the position of the peak value P of the time function s obtained by adding the square of the function f and the square of the autocorrelation function g as a synchronization point to achieve frame synchronization . 請求項1に記載の水中通信方式において、前記同期シンボルは、バーカー符号又はM系列符号であることを特徴とする水中通信方式。   The underwater communication system according to claim 1, wherein the synchronization symbol is a Barker code or an M-sequence code. 請求項1又は請求項2に記載の水中通信方式において、受信側では複数の受信チャンネルを具備しており、各チャンネルにおいて前記ピーク値Pを検出し、各チャンネルの前記ピーク値Pが最も大きいチャンネルにおける前記時間関数sの、前記ピーク値Pの位置を同期点として検出して、前記同期点で全チャンネルのフレーム同期をとることを特徴とする水中通信方式。   3. The underwater communication system according to claim 1 or 2, wherein the receiving side has a plurality of receiving channels, the peak value P is detected in each channel, and the channel having the largest peak value P in each channel. The underwater communication system is characterized in that the position of the peak value P of the time function s is detected as a synchronization point, and frame synchronization of all channels is performed at the synchronization point. 請求項1又は請求項2に記載の水中通信方式において、受信側では複数の受信チャンネルを具備しており、前記時間関数sのピーク値Pと、前記ピーク値Pの時間を含む一定時間範囲内における前記時間関数sの平均値Mとの比P/Mを各チャンネルについて求め、前記比P/Mの値が最も大きいチャンネルにおける前記時間関数sの、前記ピーク値Pの位置を同期点として検出して、前記同期点で全チャンネルのフレーム同期をとることを特徴とする水中通信方式。   3. The underwater communication system according to claim 1 or 2, wherein the receiving side includes a plurality of receiving channels, and the peak value P of the time function s and the time of the peak value P are within a certain time range. A ratio P / M with respect to the average value M of the time function s in the channel is obtained for each channel, and the position of the peak value P of the time function s in the channel having the largest value of the ratio P / M is detected as a synchronization point. An underwater communication system characterized in that frame synchronization of all channels is performed at the synchronization point. 請求項1又は請求項2に記載の水中通信方式において、連続する2つの前記送信データフレームにおける前記同期点の時間間隔と、前記送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、該当のデータフレームについて送信側に再送要求することを特徴とする水中通信方式。   The underwater communication system according to claim 1 or 2, wherein an absolute value of a difference between a time interval of the synchronization point in two consecutive transmission data frames and a transmission interval of the transmission data frames is determined in advance. An underwater communication system, which requests retransmission of the corresponding data frame to the transmitting side when the value is larger than the specified value. 請求項1又は請求項2に記載の水中通信方式において、連続する2つの前記送信データフレームにおける前記同期点の時間間隔と、前記送信データフレームの送信間隔との差の絶対値が、予め定められた値より大きい場合に、一つ目の前記送信データフレームの同期点と前記送信間隔から、二つ目の前記送信データフレームの同期点を決定してフレーム同期をとることを特徴とする水中通信方式。   The underwater communication system according to claim 1 or 2, wherein an absolute value of a difference between a time interval of the synchronization point in two consecutive transmission data frames and a transmission interval of the transmission data frames is determined in advance. A frame synchronization is determined by determining the synchronization point of the second transmission data frame from the synchronization point of the first transmission data frame and the transmission interval. method. 請求項1又は請求項2に記載の水中通信方式において、適応等化器を備えることを特徴とする水中通信方式。   The underwater communication system according to claim 1 or 2, further comprising an adaptive equalizer. 請求項7に記載の水中通信方式において、前記適応等化器は1チャンネル又はマルチチャンネルの判定帰還形適応等化器であることを特徴とする水中通信方式。
8. The underwater communication system according to claim 7, wherein the adaptive equalizer is a one-channel or multi-channel decision feedback type adaptive equalizer.
JP2004295516A 2004-10-08 2004-10-08 Underwater communication method Expired - Fee Related JP4341522B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004295516A JP4341522B2 (en) 2004-10-08 2004-10-08 Underwater communication method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004295516A JP4341522B2 (en) 2004-10-08 2004-10-08 Underwater communication method

Publications (3)

Publication Number Publication Date
JP2006109279A true JP2006109279A (en) 2006-04-20
JP2006109279A5 JP2006109279A5 (en) 2006-12-28
JP4341522B2 JP4341522B2 (en) 2009-10-07

Family

ID=36378448

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004295516A Expired - Fee Related JP4341522B2 (en) 2004-10-08 2004-10-08 Underwater communication method

Country Status (1)

Country Link
JP (1) JP4341522B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120062587A (en) * 2010-12-06 2012-06-14 강릉원주대학교산학협력단 Apparatus and method for underwater wireless communications
CN103618575A (en) * 2013-12-11 2014-03-05 厦门大学 Automatic real-time frame synchronization method of frequency-modulated hydro-acoustic communication system
KR101917503B1 (en) * 2017-04-06 2018-11-09 국방과학연구소 Multiple underwater target distance estimation method and system using underwater acoustic signal included orthogonal code

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120062587A (en) * 2010-12-06 2012-06-14 강릉원주대학교산학협력단 Apparatus and method for underwater wireless communications
CN103618575A (en) * 2013-12-11 2014-03-05 厦门大学 Automatic real-time frame synchronization method of frequency-modulated hydro-acoustic communication system
CN103618575B (en) * 2013-12-11 2015-07-15 厦门大学 Automatic real-time frame synchronization method of frequency-modulated hydro-acoustic communication system
KR101917503B1 (en) * 2017-04-06 2018-11-09 국방과학연구소 Multiple underwater target distance estimation method and system using underwater acoustic signal included orthogonal code

Also Published As

Publication number Publication date
JP4341522B2 (en) 2009-10-07

Similar Documents

Publication Publication Date Title
Diamant et al. Low probability of detection for underwater acoustic communication: A review
JP7182806B2 (en) positioning system
US9261581B2 (en) Location system
US9503202B2 (en) Underwater acoustic array, communication and location system
US7310286B1 (en) System for undersea digital acoustic communications
KR101580427B1 (en) Doppler Frequency Estimation and receiving method for Time-varying Underwater Acoustic Communication Channel
JP2008306433A (en) Communication system, communication method, and base station apparatus
US20170026135A1 (en) Underwater acoustic array, communication and location system
JP2014106058A (en) Detection and ranging device, and interference signal identification method
Hursky et al. Point-to-point underwater acoustic communications using spread-spectrum passive phase conjugation
JP4341522B2 (en) Underwater communication method
US8699300B2 (en) Underwater communications
JP4333475B2 (en) Underwater communication system
US7218574B2 (en) High range rate signaling
JP2010060520A (en) Modulating/demodulating method for ultrasonic waves, distance detecting method, and communication method
RU2700005C1 (en) Method of estimating channel parameters in broadband hydroacoustic communication and a device for realizing said channel
CN114879143A (en) Cable underwater robot navigation positioning system and positioning method
Beaujean et al. HERMES—A high bit-rate underwater acoustic modem operating at high frequencies for ports and shallow water applications
Hursky et al. Passive phase-conjugate signaling using pulse-position modulation
JP4800886B2 (en) Receiving apparatus and receiving method
JP6051787B2 (en) Communication apparatus and communication method
Kuznetsov et al. Performance of the OFDM underwater acoustic communication system with high-rank subcarrier modulation and active beam pattern control
KR101790224B1 (en) Method for estimating synchronization position of training sequence and receiver using the same
US20240015477A1 (en) Method and Apparatus of An Automated Safety Response System in a Self-organizing, multi-networked cooperative NvisiLink Mesh with Echo Positioning
Pennec et al. On joint acoustic communication and positioning through MFSK-modulated signals and Costas arrays

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20060425

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061114

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061114

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20081209

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090106

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090303

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090324

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090525

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090616

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090629

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120717

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4341522

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120717

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130717

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees