JP3032776B2 - Passage position detection method and passage position detection device - Google Patents

Passage position detection method and passage position detection device

Info

Publication number
JP3032776B2
JP3032776B2 JP35059297A JP35059297A JP3032776B2 JP 3032776 B2 JP3032776 B2 JP 3032776B2 JP 35059297 A JP35059297 A JP 35059297A JP 35059297 A JP35059297 A JP 35059297A JP 3032776 B2 JP3032776 B2 JP 3032776B2
Authority
JP
Japan
Prior art keywords
magnetic
value
axis
signal
detector
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.)
Expired - Lifetime
Application number
JP35059297A
Other languages
Japanese (ja)
Other versions
JPH11183109A (en
Inventor
重光 倉野
好之 野本
邦彦 吉見
晴之 森下
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.)
Kyosan Electric Manufacturing Co Ltd
Original Assignee
Kyosan Electric Manufacturing Co 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 Kyosan Electric Manufacturing Co Ltd filed Critical Kyosan Electric Manufacturing Co Ltd
Priority to JP35059297A priority Critical patent/JP3032776B2/en
Publication of JPH11183109A publication Critical patent/JPH11183109A/en
Application granted granted Critical
Publication of JP3032776B2 publication Critical patent/JP3032776B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、磁気的条件が未知
の磁性体目標物が、その位置や状態が静止し定点である
磁気検出器から一定範囲内を通過したかどうかをその磁
性体の大小にかかわらず確実に検出できるようにした通
過位置検出方法及び通過位置検出装置に関するものであ
る。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for determining whether a magnetic target whose magnetic condition is unknown has passed a fixed range from a magnetic detector whose position or state is stationary and fixed. The present invention relates to a passing position detecting method and a passing position detecting device which can surely detect a passing position regardless of the size.

【0002】[0002]

【従来の技術】地磁気中の磁性体は磁性体自身の磁化或
いは地磁気によって磁化されている。従って磁性体から
の磁気を検出して当該磁性体の存否を検出することはで
きるが、磁性体の磁化の程度が様々な場合には磁気検出
器の側方、すなわち直横に対しどの程度の距離を通過し
たかを検出することは不可能であった。
2. Description of the Related Art A magnetic body in the terrestrial magnetism is magnetized by its own magnetization or terrestrial magnetism. Therefore, the presence or absence of the magnetic material can be detected by detecting the magnetism from the magnetic material.However, when the degree of magnetization of the magnetic material is various, the degree of the side of the magnetic detector, that is, It was impossible to detect whether the vehicle had passed the distance.

【0003】これを解決するものとして特開昭55−1
62078号が提案されている。すなわち3軸磁気検出
器を用いたこの磁性体の検出方法は、Z軸を垂直軸とす
る直交座標系(X,Y,Z)の3軸磁気検出器に対しZ
軸方向が一定のXY平面を移動する磁性体目標物のZ軸
方向最大磁気成分(絶対値)Hv とその時のXY平面の
磁気成分Hh との比C=Hv /Hh をとり、当該比Cを
しきい値と比べることにより磁気検出器の側方に対しど
の程度の距離を通過したかを検出するようにしている。
[0003] To solve this problem, Japanese Patent Application Laid-Open No. Sho 55-1
No. 62078 has been proposed. That is, this magnetic body detection method using a three-axis magnetic detector is based on a three-axis magnetic detector of a rectangular coordinate system (X, Y, Z) having the Z axis as a vertical axis.
The ratio C = Hv / Hh of the maximum magnetic component (absolute value) Hv in the Z-axis direction of the magnetic target moving on the XY plane having a constant axial direction and the magnetic component Hh on the XY plane at that time is calculated. By comparing with a threshold value, it is detected how far the side of the magnetic detector has passed.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、磁性体
のZ軸方向最大磁気成分(絶対値)Hv とその時のXY
平面の磁気成分Hh との比C=Hv /Hh をとり、当該
比Cをしきい値と比べる従来の磁性体の検出方法では、
直上通過時の磁気測定、例えば事前の測定による基準値
の決定を必要とすると共に、XY平面の磁気成分Hh が
0のとき当該比C=Hv /Hh が無限大となる欠点があ
り、また磁気的条件、すなわち磁性体目標物の磁気能率
の伏角や偏角が検出結果の精度に影響を及ぼす欠点があ
った。
However, the maximum magnetic component (absolute value) Hv in the Z-axis direction of the magnetic material and XY at that time
In a conventional method for detecting a magnetic substance, a ratio C = Hv / Hh with respect to a plane magnetic component Hh is taken, and the ratio C is compared with a threshold value.
There is a drawback that the magnetic measurement at the time of passing right above, for example, the determination of a reference value by prior measurement is required, and when the magnetic component Hh on the XY plane is 0, the ratio C = Hv / Hh becomes infinite. However, there is a disadvantage that the magnetic conditions, that is, the inclination and the declination of the magnetic efficiency of the magnetic target affect the accuracy of the detection result.

【0005】本発明は、上記の欠点を解決することを目
的としており、事前の磁気測定を必要とせず、演算処理
の過程で無限大が生じる処理がなく、かつ船体の大きさ
や地磁気などの磁性体の磁気的条件の影響を受けにくい
磁性体目標物の通過位置検出方法及び通過位置検出装置
を提供することを目的としている。
SUMMARY OF THE INVENTION An object of the present invention is to solve the above-mentioned drawbacks, does not require a prior magnetic measurement, has no processing in which infinity occurs in the course of arithmetic processing, and has a magnetic property such as the size of a hull and geomagnetism. It is an object of the present invention to provide a method and a device for detecting a passing position of a magnetic target which is hardly affected by magnetic conditions of a body.

【0006】[0006]

【課題を解決するための手段】上記の目的を解決するた
めに、本発明の通過位置検出方法は、Z軸を垂直軸とす
る直交座標系(X,Y,Z)の各軸上に指向性を有する
磁気センサがそれぞれ配置されてなるその位置が固定の
3軸磁気検出器を備え、当該3軸磁気検出器に対しZ軸
方向が一定のXY平面を移動する磁性体目標物の予め定
められた範囲内の通過を検出する通過位置検出方法にお
いて、3軸磁気検出器のサンプリングされた各検出信号
から磁性体目標物の合成磁気信号値を求めると共に、そ
の最大値となるタイミングを求め、サンプリング信号毎
に、Z軸磁気センサの検出信号の絶対値|Hz|から、
X軸磁気センサ及びY軸磁気センサの検出信号から得ら
れる水平成分値√(Hx2+Hy2)を減じた値Hn2(i)
=|Hz|−√(Hx2+Hy2)(i=1,2,……,
n)を演算してその演算結果を記憶手段に記憶し、磁性
体目標物の合成磁気信号値が最大となったとき、記憶手
段に最も新しく記憶されたものから順に前に逆上り記憶
手段に記憶されている予め定められた数の各値Hn2を合
成磁気信号の最大値Ht0でそれぞれ割った値S(i)=
Hn2/Ht0を求め、この各値S(i)=Hn2/Ht0と、
予め定められた重み付け値とを基にニューラルネットワ
ーク演算処理を行って波形識別を行い、波形識別の演算
値から、磁性体目標物の予め定められた範囲内の通過を
検出するようにしたことを特徴としている。
In order to solve the above-mentioned object, a passing position detecting method according to the present invention is directed to a method of detecting a passing position on each axis of a rectangular coordinate system (X, Y, Z) having a Z axis as a vertical axis. A magnetic target having a fixed position and a three-axis magnetic detector having a fixed position, and a magnetic target moving in an XY plane having a constant Z-axis direction with respect to the three-axis magnetic detector. In a passing position detection method for detecting passing within a given range, a combined magnetic signal value of a magnetic target is obtained from each sampled detection signal of the three-axis magnetic detector, and a timing at which the maximum value is obtained, From the absolute value | Hz | of the detection signal of the Z-axis magnetic sensor for each sampling signal,
A value Hn 2 (i) obtained by subtracting the horizontal component value √ (Hx 2 + Hy 2 ) obtained from the detection signals of the X-axis magnetic sensor and the Y-axis magnetic sensor.
= | Hz | −√ (Hx 2 + Hy 2 ) (i = 1, 2,...,
n) is calculated and the result of the calculation is stored in the storage means. When the composite magnetic signal value of the magnetic target becomes maximum, the reversely stored data is stored in the storage means in order from the most recently stored in the storage means. A value S (i) = a predetermined number of stored values Hn 2 divided by the maximum value H t0 of the composite magnetic signal, respectively.
Seeking Hn 2 / H t0, this values S (i) = Hn 2 / H t0,
Performing neural network arithmetic processing based on a predetermined weight value and performing waveform identification, and detecting the passage of a magnetic target within a predetermined range from the arithmetic value of the waveform identification. Features.

【0007】そして本発明の通過位置検出装置は、Z軸
を垂直軸とする直交座標系(X,Y,Z)の各軸上に指
向性を有する磁気センサがそれぞれ配置されてなるその
位置が固定の3軸磁気検出器を備え、当該3軸磁気検出
器に対しZ軸方向が一定のXY平面を移動する磁性体目
標物の予め定められた範囲内の通過を検出する通過位置
検出装置において、予め定められたレベル以上の磁気信
号を検出したとき出力を発生する磁気レベル検出器と、
磁気レベル検出器の出力で3軸磁気検出器の各検出信号
をサンプリングさせるサンプリング回路と、サンプリン
グ回路が出力するサンプリング信号毎に、3軸磁気検出
器の各検出信号から磁性体目標物の合成磁気信号を求め
る磁気信号合成手段と、磁気信号合成手段で求められた
合成磁気信号を基に、当該合成磁気信号の最大値Ht0
検出する最大値検出手段と、サンプリング回路が出力す
るサンプリング信号毎に、Z軸磁気センサの検出信号の
絶対値|Hz|から、X軸磁気センサ及びY軸磁気セン
サの検出信号から得られる水平成分値√(Hx2+Hy2
を減じた値Hn2(i)=|Hz|−√(Hx2+Hy2
(i=1,2,……,n)を演算する演算手段と、演算
手段の演算結果を記憶する記憶手段と、最大値検出手段
が上記合成磁気信号の最大値Ht0を検出したとき、当該
記憶手段に最も新しく記憶されたものから順に前に逆上
り記憶されている予め定められた数の各値Hn2を読み出
し、この値Hn2を最大値検出手段で検出された合成磁気
信号の最大値Ht0で割る計算Hn2/Ht0を実行する計算
手段と、計算手段で計算された各計算値S(i)=Hn2
/Ht0と予め定められた重み付け値とを基にニューラル
ネットワーク演算処理を行い、波形識別を行う波形識別
処理手段と、波形識別処理手段が処理した演算値と予め
定められた基準値とを比較する比較手段とを備えたこと
を特徴としている。
In the passing position detecting device of the present invention, the position at which the magnetic sensors having directivity are arranged on each axis of the orthogonal coordinate system (X, Y, Z) having the Z axis as the vertical axis is determined. A passing position detecting device including a fixed three-axis magnetic detector and detecting the passage of a magnetic target that moves on an XY plane in which the Z-axis direction is constant with respect to the three-axis magnetic detector within a predetermined range. A magnetic level detector that generates an output when detecting a magnetic signal of a predetermined level or more,
A sampling circuit for sampling each detection signal of the three-axis magnetic detector with an output of the magnetic level detector; and for each sampling signal output by the sampling circuit, a composite magnetism of a magnetic target object from each detection signal of the three-axis magnetic detector. A magnetic signal synthesizing means for obtaining a signal, a maximum value detecting means for detecting a maximum value H t0 of the synthetic magnetic signal based on the synthetic magnetic signal obtained by the magnetic signal synthesizing means, and a sampling signal output by the sampling circuit. From the absolute value | Hz | of the detection signal of the Z-axis magnetic sensor, the horizontal component value √ (Hx 2 + Hy 2 ) obtained from the detection signals of the X-axis magnetic sensor and the Y-axis magnetic sensor
Hn 2 (i) = | Hz | −√ (Hx 2 + Hy 2 )
(I = 1, 2,..., N), storage means for storing the calculation result of the calculation means, and when the maximum value detection means detects the maximum value H t0 of the composite magnetic signal, A predetermined number of each value Hn 2 stored in the storage unit in reverse order is read out in order from the most recently stored in the storage unit, and this value Hn 2 is read out of the composite magnetic signal detected by the maximum value detection unit. Calculating means for executing the calculation Hn 2 / H t0 divided by the maximum value H t0 , and each calculated value S (i) = Hn 2 calculated by the calculating means
/ H t0 and a waveform identification processing means for performing a neural network arithmetic processing based on a predetermined weighting value to identify a waveform, and comparing the arithmetic value processed by the waveform identification processing means with a predetermined reference value And comparing means.

【0008】ここで図17を用いて、磁性体目標物15
をその例として船体としたときの進行方向とその船体磁
気との関係について説明すると、図17(I)は磁性体
目標物15の針路が北の方向、(II)は磁性体目標物
15の針路が南の方向に向かっている場合をそれぞれ示
しており、(I),(II)の各上側の図は磁性体目標
物15を上らか見たときのもの、その各下側の図は磁性
体目標物15を横から見たときのものをそれぞれ表して
いる。
[0008] Here, referring to FIG.
17 (I), the course of the magnetic target 15 is in the north direction, and (II) is the direction of the magnetic target 15. The upper part of each of (I) and (II) shows a case where the course is heading in the south direction, and the lower part of FIG. Represents the magnetic substance target 15 when viewed from the side.

【0009】磁性体目標物15の船体磁気の永久成分H
P (図17(I)の(イ))、地磁気の方向およびその
大きさを点線で表したときの当該磁性体目標物15に誘
導されるその誘導成分をHI (図17(I)の(ロ))
としたとき、磁性体目標物15の船体磁気はこれらの合
成成分として求められ、図17(I)の(ハ)のHP
I (HP ,HI はいずれもベクトル)で表わすことが
できる。
The permanent component H of the hull magnetism of the magnetic target 15
P ((a) in FIG. 17 (I)), and the induction component guided to the magnetic target 15 when the direction and the magnitude of the geomagnetism are represented by dotted lines are H I (FIG. 17 (I)). (B))
When the hull magnetic magnetic body target 15 is determined as the synthetic components, FIG. 17 (I) (iii) H P +
H I (H P and H I are both vectors).

【0010】一方、磁性体目標物15の針路が図17
(I)と逆方向の南の方向に向かっている場合、磁性体
目標物15の船体磁気の永久成分HP (図17(II)
の(イ))、その位置での地磁気の方向およびその大き
さは点線図示の如く変わらないので、このときの磁性体
目標物15の船体磁気は、上記船体磁気の永久成分HP
(図17(II)の(イ))と当該磁性体目標物15に
誘導されるその誘導成分HI (図17(II)の
(ロ))との合成成分となり、図17(II)の(ハ)
P +HI となる。
On the other hand, the course of the magnetic target 15 is shown in FIG.
If toward the opposite direction of the south direction (I), permanently component hull magnetic magnetic body target 15 H P (FIG. 17 (II)
Of (b)), since the geomagnetic direction and its magnitude at that position is not changed as dotted lines shown, the hull magnetic magnetic body target 15 at this time, the permanent components of the hull magnetic H P
((A) of FIG. 17 (II)) and its induction component H I ((b) of FIG. 17 (II)) induced to the magnetic substance target 15 and become a composite component. (C)
The H P + H I.

【0011】一般に磁性体目標物15の船体磁気の大き
さは、磁性体目標物15の大きさ、つまり船体の大きさ
に比例する。しかしながら、上述の図17の(I),
(II)から明らかなようにその進行方向がどちらを向
いていてもZ方向成分の大きさは同一となり、Z成分に
なおしてしまえば磁性体目標物15の進行方向について
は変わりはなくなり、Hn2=|Hz|−√(Hx2+H
y2)は磁性体目標物15の針路方向に関係がなくなる。
In general, the magnitude of the hull magnetism of the magnetic target 15 is proportional to the size of the magnetic target 15, that is, the size of the hull. However, FIG.
As is apparent from (II), the magnitude of the Z-direction component becomes the same regardless of the direction of travel, and the travel direction of the magnetic target 15 remains unchanged if the Z component is replaced. 2 = | Hz | -√ (Hx 2 + H
y 2 ) has no relation to the course direction of the magnetic target 15.

【0012】そこで、本発明では、船体磁気を磁気双極
子モーメント(磁気能率)で表示すると共に、3軸磁気
検出器で検出された各検出信号を演算処理する際、船体
の大きさに無関係とすべく値Hn2(i)=|Hz|−√
(Hx2+Hy2)(i=1,2,……,n)を合成磁気信
号の最大値Ht0で割って無次元化を行い、その上でニュ
ーラルネットワーク演算を用いた波形識別処理により波
形識別を行い、3軸磁気検出器から一定範囲内であると
する基準値と比較するようにしているので、磁気的条件
が未知の磁性体目標物が3軸磁気検出器から一定範囲内
を通過したかどうかをその磁性体の大小にかかわらず、
また磁気的条件が未知の場合でも確実に検出できる。
Therefore, in the present invention, the hull magnetism is displayed as a magnetic dipole moment (magnetic moment), and when each detection signal detected by the three-axis magnetic detector is processed, it is independent of the size of the hull. The value Hn 2 (i) = | Hz | −√
(Hx 2 + Hy 2 ) (i = 1, 2,..., N) is divided by the maximum value H t0 of the composite magnetic signal to make the dimensionless, and then the waveform is determined by the waveform identification processing using the neural network operation. Since the target is identified and compared with a reference value that is within a certain range from the three-axis magnetic detector, a magnetic target whose magnetic condition is unknown passes within a certain range from the three-axis magnetic detector. Regardless of the size of the magnetic material,
Further, even if the magnetic condition is unknown, it can be reliably detected.

【0013】[0013]

【発明の実施の形態】図1は本発明に係る通過位置検出
装置の一実施例構成を示している。同図において、3軸
磁気検出器1はZ軸を垂直軸とする直交座標系(X,
Y,Z)の各軸上に指向性を有するX軸磁気センサ2,
Y軸磁気センサ3,Z軸磁気センサ4を備えている。こ
れらの磁気センサは、例えばフラックスゲート型やリン
グコア型のものが用いられる。
FIG. 1 shows the configuration of an embodiment of a passing position detecting device according to the present invention. In FIG. 1, a three-axis magnetic detector 1 has a rectangular coordinate system (X,
X-axis magnetic sensor having directivity on each of Y, Z) axes 2,
A Y-axis magnetic sensor 3 and a Z-axis magnetic sensor 4 are provided. As these magnetic sensors, for example, a flux gate type or a ring core type is used.

【0014】3軸磁気検出器1が予め定められたレベル
以上の磁気信号を検出したとき、出力を発生する磁気レ
ベル検出器5が設けられており、当該磁気レベル検出器
5の出力発生でサンプリング回路6が駆動され、サンプ
リング回路6からサンプリング信号が出力するようにな
っている。
When the three-axis magnetic detector 1 detects a magnetic signal of a predetermined level or more, a magnetic level detector 5 for generating an output is provided, and sampling is performed when the magnetic level detector 5 generates an output. The circuit 6 is driven, and the sampling circuit 6 outputs a sampling signal.

【0015】X軸磁気センサ2,Y軸磁気センサ3及び
Z軸磁気センサ4には、検出されたそれぞれの信号をそ
れぞれディジタル化するA/D変換器7が対応して設け
られている。A/D変換器7でディジタル化されたX軸
磁気センサ2,Y軸磁気センサ3及びZ軸磁気センサ4
の検出信号は、磁気信号合成手段8と演算手段9とに入
力されるようになっている。
The X-axis magnetic sensor 2, the Y-axis magnetic sensor 3, and the Z-axis magnetic sensor 4 are provided with corresponding A / D converters 7 for digitizing the respective detected signals. X-axis magnetic sensor 2, Y-axis magnetic sensor 3, and Z-axis magnetic sensor 4 digitized by A / D converter 7
Is input to the magnetic signal synthesizing means 8 and the calculating means 9.

【0016】いま、図3に示されている様に、磁性体目
標物15の磁気能率をM、磁気能率Mの伏角,偏角をそ
れぞれβ,γとし、図4図示の如く、3軸磁気検出器1
に対しZ軸方向が一定のXY平面を定速度Vで移動する
磁性体目標物15について、X軸磁気センサ2で検出さ
れた検出信号をHx、Y軸磁気センサ3で検出された検
出信号をHy、Z軸磁気センサ4で検出された検出信号
をHzとしたとき、磁気信号合成手段8は各磁気センサ
の検出信号を基に、当該磁性体目標物15の合成磁気信
号値√(Hx2+Hy2+Hz2)を演算するようになってい
る。
Now, as shown in FIG. 3, the magnetic moment of the magnetic target 15 is represented by M, and the inclination and declination of the magnetic moment M are represented by β and γ, respectively, as shown in FIG. Detector 1
On the other hand, the detection signal detected by the X-axis magnetic sensor 2 is Hx, and the detection signal detected by the Y-axis magnetic sensor 3 is Hx. Assuming that the detection signal detected by the Hy and Z axis magnetic sensors 4 is Hz, the magnetic signal synthesizing means 8 generates a composite magnetic signal value √ (Hx 2) of the magnetic target 15 based on the detection signals of the respective magnetic sensors. + Hy 2 + Hz 2 ).

【0017】磁気レベル検出器5の出力信号をサンプリ
ング回路6が受けると、当該サンプリング回路6は磁気
信号合成手段8と演算手段9とにサンプリング信号を送
出し、磁気信号合成手段8及び演算手段9に上記A/D
変換器7でディジタル化された各検出信号Hx、Hy、
Hzを取り込ませる。
When the output signal of the magnetic level detector 5 is received by the sampling circuit 6, the sampling circuit 6 sends the sampling signal to the magnetic signal synthesizing means 8 and the arithmetic means 9, and the magnetic signal synthesizing means 8 and the arithmetic means 9 The above A / D
Each detection signal Hx, Hy, digitized by the converter 7
Hz.

【0018】磁気信号合成手段8は、サンプリング回路
6が出力するサンプリング信号毎に、取り込まれる各検
出信号Hx、Hy、Hzから上記の合成磁気信号値√
(Hx2+Hy2+Hz2)を演算し、次段に設けられた最大
値検出手段10に当該合成磁気信号値√(Hx2+Hy2
Hz2)を送出する。
The magnetic signal synthesizing means 8 calculates, for each sampling signal output from the sampling circuit 6, the above-mentioned synthesized magnetic signal value か ら from each of the captured detection signals Hx, Hy and Hz.
(Hx 2 + Hy 2 + Hz 2 ) is calculated, and the maximum value detection means 10 provided at the next stage is used to calculate the composite magnetic signal value √ (Hx 2 + Hy 2 +
Hz 2 ).

【0019】演算手段9は、サンプリング回路6が出力
するサンプリング信号毎に、取り込まれる各検出信号H
x、Hy、Hzを基に、Z軸磁気センサ4の検出信号の
絶対値|Hz|から、X軸磁気センサ2及びY軸磁気セ
ンサ3の検出信号から得られる水平成分値√(Hx2+H
y2)を減じた値Hn2、すなわちHn2(i)=|Hz|−
√(Hx2+Hy2)(i=1,2,……,n)の演算を行
い、その演算結果Hn2(i)をメモリ11に順に記憶さ
せる。
The calculating means 9 is provided for each sampling signal output from the sampling circuit 6 and for each detection signal H taken in.
The horizontal component value √ (Hx 2 + H) obtained from the detection signals of the X-axis magnetic sensor 2 and the Y-axis magnetic sensor 3 from the absolute value | Hz | of the detection signal of the Z-axis magnetic sensor 4 based on x, Hy, and Hz.
The value Hn 2 obtained by subtracting the y 2), i.e. Hn 2 (i) = | Hz | -
The operation of √ (Hx 2 + Hy 2 ) (i = 1, 2,..., N) is performed, and the operation result Hn 2 (i) is sequentially stored in the memory 11.

【0020】最大値検出手段10は、磁気信号合成手段
8から送られてくる合成磁気信号値Ht(k)=√(H
x2+Hy2+Hz2)と一つ前に取り込まれた合成磁気信号
値Ht(k−1)=√(Hx2+Hy2+Hz2)とを順次比
較し、合成磁気信号値√(Hx2+Hy2+Hz2)の最大値
t0を求めるようになっている。
The maximum value detecting means 10 generates a composite magnetic signal value Ht (k) = √ (H
x 2 + Hy 2 + Hz 2 ) and the previously captured combined magnetic signal value Ht (k−1) = √ (Hx 2 + Hy 2 + Hz 2 ) are sequentially compared, and the combined magnetic signal value √ (Hx 2 + Hy) is compared. and it obtains the maximum value H t0 of 2 + Hz 2).

【0021】そして最大値検出手段10が合成磁気信号
値√(Hx2+Hy2+Hz2)の最大値Ht0を求めたとき、
すなわち最大値検出手段10が上記合成磁気信号の最大
値H t0を検出したとき、最大値検出手段10は上記合成
磁気信号の最大値Ht0を検出した旨の信号とその合成磁
気信号の最大値Ht0とを計算手段12に送出する。
The maximum value detecting means 10 outputs the composite magnetic signal
Value √ (HxTwo+ HyTwo+ HzTwo) Maximum value Ht0When asked for
That is, the maximum value detecting means 10 determines the maximum value of the composite magnetic signal.
Value H t0Is detected, the maximum value detecting means 10
Maximum value of magnetic signal Ht0Signal of the detection of
Maximum value of air signal Ht0Are sent to the calculating means 12.

【0022】計算手段12は、メモリ11に記憶された
上記演算手段9による演算結果の各値Hn2を、最も新し
く記憶されたものから順に前に逆上り予め定められた数
Nだけメモリ11から順に読み出すと共に、この各値H
n2を最大値検出手段10で検出された合成磁気信号の最
大値Ht0で割る計算のHn2/Ht0を順に行い、その各計
算値S(i)=Hn2/Ht0(i=1,2,……,N)を
波形識別処理手段13に入力させるようになっている。
The calculating means 12 reverses the values Hn 2 of the result of the calculation by the calculating means 9 stored in the memory 11 by a predetermined number N in order from the most recently stored one to the previous one. The values are read out in order and
n 2 the turn performs Hn 2 / H t0 calculation dividing the maximum value H t0 of the detected resultant magnetic signal with a maximum value detecting means 10, the respective calculated values S (i) = Hn 2 / H t0 (i = ,..., N) are input to the waveform identification processing means 13.

【0023】この各値Hn2を合成磁気信号の最大値Ht0
で割ることにより、無次元化された磁性体モデルとする
ことができ、磁性体目標物15の磁気能率Mの大きさに
関係なく取り扱えるようになる。
The respective values Hn 2 are converted to the maximum value H t0 of the composite magnetic signal.
By dividing by, the dimensionless magnetic material model can be obtained, and the magnetic material target 15 can be handled regardless of the magnitude of the magnetic efficiency M.

【0024】波形識別処理手段13では、計算手段12
から入力された所定数N個の値Hn2と予め定められた重
み付け値とを基にニューラルネットワーク演算処理を行
い、積と和からなる演算の波形識別のための波形識別処
理を行う。波形識別処理手段13で処理された演算結果
値は比較手段14に出力される。この波形識別処理手段
13については次の図2で詳しく説明する。
In the waveform discriminating means 13, the calculating means 12
Based on a predetermined number N values Hn 2 entered the predetermined weighting value performs neural network processing from the waveform discrimination processing for waveform identification operation consisting of product and sum. The calculation result value processed by the waveform identification processing means 13 is output to the comparison means 14. This waveform identification processing means 13 will be described in detail with reference to FIG.

【0025】比較手段14には、3軸磁気検出器1の側
方に対しどの程度の距離までの通過を、磁性体目標物1
5の予め定められた範囲内の通過とするかを意味する値
が予め基準値として設定されており、当該基準値と波形
識別処理手段13で処理された演算値とが比較手段14
で比較される。波形識別処理手段13で処理された演算
値が基準値より大きいとき、磁性体目標物15が予め定
められた範囲内の通過を意味する「1」が比較手段14
から出力され、波形識別処理手段13で処理された演算
値が基準値より小さいとき、磁性体目標物15が予め定
められた範囲外の通過を意味する「0」が比較手段14
から出力されるようになっている。
The comparison means 14 determines how far the side of the three-axis magnetic detector 1 should pass to reach the magnetic target 1
5 is set in advance as a reference value, and the calculated value processed by the waveform identification processing unit 13 is compared with the comparison value by the comparison unit 14.
Are compared. When the calculated value processed by the waveform identification processing means 13 is larger than the reference value, “1” meaning that the magnetic substance target 15 passes within a predetermined range is set to the comparison means 14.
When the calculated value output from the waveform identification processing means 13 is smaller than the reference value, “0” indicating that the magnetic target 15 passes outside a predetermined range is set to the comparison means 14.
Output from

【0026】図2は波形識別処理手段の一実施例概念説
明図を示しており、波形識別処理手段13としてニュー
ラルネットワークが用いられたものである。同図におい
て、波形識別処理手段13のニューラルネットワーク
は、N個の入力ユニット21−1ないし21−N、N個
の隠れユニット22−1ないし22−N及び1個の出力
ユニット23で構成されている。そして入力ユニット2
1−1ないし21−Nには、図1で説明した様に最大値
検出手段10が合成磁気信号の最大値Ht0を検出したと
き、メモリ11に格納されている最も新しく記憶された
ものから順に前に逆上り、予め定められた数N個の各値
Hn2を合成磁気信号の最大値Ht0で割った値S(i)=
Hn2/Ht0(i=1,2,……,N)が、それぞれ入力
される。
FIG. 2 is a conceptual explanatory diagram of one embodiment of the waveform discrimination processing means, in which a neural network is used as the waveform discrimination processing means 13. In the figure, the neural network of the waveform identification processing means 13 is composed of N input units 21-1 to 21-N, N hidden units 22-1 to 22-N, and one output unit 23. I have. And input unit 2
1-1 to 21-N, when the maximum value detecting means 10 detects the maximum value Ht0 of the composite magnetic signal as described with reference to FIG. Going back up in order, a value S (i) = a predetermined number N of values Hn 2 divided by the maximum value H t0 of the composite magnetic signal =
Hn 2 / H t0 (i = 1, 2,..., N) are input, respectively.

【0027】最大値検出手段10が合成磁気信号の最大
値Ht0を検出したとき、メモリ11の中で最も新しく記
憶された演算結果の値をHn2(t1 )とすると、計算手
段12によって求められた計算結果の値S(t1 )=H
n2(t1 )/Ht0が入力ユニット21−1に入力され、
メモリ11の中で次に新しい演算結果の値S(t2 )=
Hn2(t2 )/Ht0が入力ユニット21−2に入力さ
れ、メモリ11の中でその次に新しい演算結果の値S
(t3 )=Hn2(t3 )/Ht0が入力ユニット21−3
に入力され、以下同様にして演算結果の値S(tN )=
Hn2(tN )/Ht0が入力ユニット21−Nに入力され
る。
When the maximum value detecting means 10 detects the maximum value H t0 of the composite magnetic signal, and the value of the operation result most recently stored in the memory 11 is Hn 2 (t 1 ), the calculating means 12 The calculated value S (t 1 ) = H
n 2 (t 1 ) / H t0 is input to the input unit 21-1;
The next newest operation result value S (t 2 ) in the memory 11 =
Hn 2 (t 2 ) / H t0 is inputted to the input unit 21-2, and the next new operation result value S in the memory 11 is obtained.
(T 3 ) = Hn 2 (t 3 ) / H t0 is the input unit 21-3
Is input, the following Similarly the operation result value S (t N) =
Hn 2 (t N ) / H t0 is input to the input unit 21-N.

【0028】入力ユニット21−1ないし21−Nと隠
れユニット22−1ないし22−Nとの間は、図示の如
くA1 ,A2 ,A3 ……,B1 ,B2 ,B3 ……,
1 ,C 2 ,C3 ……等の重み付け値で重み付けされて
それぞれリンクされており、また隠れユニット22−1
ないし22−Nと出力ユニット23との間も、図示の如
くM1 ,M2 ,M3 ……の重み付け値で重み付けされて
リンクされている。
The input units 21-1 to 21-N are hidden.
Between the units 22-1 to 22-N as shown in the figure.
A1, ATwo, AThree............ B1, BTwo, BThree……,
C1, C Two, CThreeWeighted by weighting values such as
Each is linked, and the hidden unit 22-1
, And between the output unit 23 and the output unit 23 as shown in FIG.
K1, MTwo, MThreeWeighted by the weight value of ......
Linked.

【0029】このような値S(t1 )ないし値S
(tN )を入力とするニューラルネットワークで重み付
け処理、すなわち積と和とからなるニューラルネットワ
ーク演算を行うことにより、出力ユニット23から或る
値の演算値が出力される。
The values S (t 1 ) to S
By performing weighting processing using a neural network having (t N ) as an input, that is, performing a neural network operation consisting of a product and a sum, an output value of a certain value is output from the output unit 23.

【0030】ニューラルネットワークで用いられる上記
の各重みA1 ,A2 ,A3 ……,B 1 ,B2 ,B3
…,C1 ,C2 ,C3 ……等とM1 ,M2 ,M3 ……
は、図5ないし図14に示された如き代表的な磁気モデ
ルでのHn2/Ht0波形パターンから求められたものであ
る。
The above used in the neural network
Weight A of1, ATwo, AThree............ B 1, BTwo, BThree
…, C1, CTwo, CThree…… etc and M1, MTwo, MThree......
Are typical magnetic models as shown in FIGS.
Hn in leTwo/ Ht0Is determined from the waveform pattern.
You.

【0031】図5ないし図14は条件を変え計算によっ
てシミュレートしたHn2/Ht0波形図をそれぞれ示しい
る。図5ないし図14において、横軸は磁性体目標物の
合成磁気信号値Ht=√(Hx2+Hy2+Hz2)が最大と
なる位置からの距離、すなわち図16の最近点Pの位置
から磁性体目標物15の進行方向と逆に向う距離(−
L)を表し、縦軸はHn2のHt0に対する割合を表してい
る。そして図5ないし図9は、磁性体目標物15の磁気
能率の伏角(β)、偏角(γ)の磁気的条件及び深度が
同一で進入角α=0°のときの横距離Wを0、10、2
0、30、40mに変えたものであり、図10ないし図
14は、磁性体目標物15の磁気能率の伏角(β)、偏
角(γ)の磁気的条件及び深度が同一で進入角α=0°
(αは図15参照)のときの横距離Wを0、10、2
0、30、40mに変えたものである。図5ないし図9
のそれぞれは、伏角(β)を変えた状態で図10ないし
図14のそれぞれに対応している。
FIGS. 5 to 14 show Hn 2 / Ht0 waveform diagrams simulated by calculation with changing conditions. In FIGS. 5 to 14, the horizontal axis represents the distance from the position at which the composite magnetic signal value Ht = √ (Hx 2 + Hy 2 + Hz 2 ) of the magnetic substance target becomes the maximum, that is, from the position of the closest point P in FIG. Distance opposite to the direction of travel of the body target 15 (−
L) represents the vertical axis represents the ratio of H t0 of Hn 2. 5 to 9 show that the lateral distance W when the approach angle α = 0 ° and the magnetic conditions of the inclination (β) and the declination (γ) of the magnetic moment of the magnetic substance target 15 are the same and the approach angle α = 0 ° is 0. , 10, 2
FIGS. 10 to 14 show that the magnetic condition of the magnetic efficiency of the magnetic target 15 is equal to the inclination (β) and the angle of deviation (γ), and the approach angle α = 0 °
(Α is shown in FIG. 15), the horizontal distance W is 0, 10, 2
It was changed to 0, 30, and 40 m. 5 to 9
Correspond to FIGS. 10 to 14 in a state where the inclination (β) is changed.

【0032】その内の図5は3軸磁気検出器1の側方
(横距離W)0mのときのもの、つまり磁性体目標物1
5が3軸磁気検出器1の真上を通過するときのもので、
図6は3軸磁気検出器1の側方(横距離W)10mのと
きのもの、図7は3軸磁気検出器1の側方(横距離W)
20mのときのもの、図8は3軸磁気検出器1の側方
(横距離W)30mのときのもの、図9は3軸磁気検出
器1の側方(横距離W)40mのときのものである。
FIG. 5 shows the case where the side (lateral distance W) of the three-axis magnetic detector 1 is 0 m, that is, the magnetic target 1.
5 when passing just above the three-axis magnetic detector 1,
FIG. 6 shows the case where the side (lateral distance W) of the three-axis magnetic detector 1 is 10 m, and FIG. 7 shows the side (lateral distance W) of the three-axis magnetic detector 1.
FIG. 8 shows the case when the side (lateral distance W) of the three-axis magnetic detector 1 is 30 m, and FIG. 9 shows the case when the side (lateral distance W) is 40 m of the three-axis magnetic detector 1. Things.

【0033】またその内の図10は3軸磁気検出器1の
側方(横距離W)0mのときのもの、つまり磁性体目標
物15が3軸磁気検出器1の真上を通過するときのもの
で、図11は3軸磁気検出器1の側方(横距離W)10
mのときのもの、図12は3軸磁気検出器1の側方(横
距離W)20mのときのもの、図13は3軸磁気検出器
1の側方(横距離W)30mのときのもの、図14は3
軸磁気検出器1の側方(横距離W)40mのときのもの
である。
FIG. 10 shows a case where the side (lateral distance W) of the three-axis magnetic detector 1 is 0 m, that is, when the magnetic target 15 passes right above the three-axis magnetic detector 1. FIG. 11 shows the side (lateral distance W) 10 of the three-axis magnetic detector 1.
m, FIG. 12 shows a case where the three-axis magnetic detector 1 is at a side (lateral distance W) of 20 m, and FIG. 13 shows a case where the three-axis magnetic detector 1 is at a side (lateral distance W) of 30 m. 14 in FIG.
This is the case when the side (lateral distance W) of the axial magnetic detector 1 is 40 m.

【0034】図16は本発明の一実施例フローチャート
を示しており、当該図16を参照しながら図1の動作を
次に説明する。波形識別処理手段13には、上記説明の
代表的な磁気モデルでのHn2/Ht0波形パターンから求
めた各重みA1 ,A2 ,A3 ……,B1 ,B2 ,B3
…,C 1 ,C2 ,C3 ……等とM1 ,M2 ,M3 ……が
予め設定されている。比較手段14には、3軸磁気検出
器1の側方に対しどの程度の距離までの通過を、磁性体
目標物15の予め定められた範囲内の通過とするか定め
る基準値が予め設定されている。
FIG. 16 is a flowchart of one embodiment of the present invention.
The operation of FIG. 1 is described with reference to FIG.
This will be described next. The waveform identification processing means 13 includes
Hn in typical magnetic modelTwo/ Ht0From the waveform pattern
Each weight A1, ATwo, AThree............ B1, BTwo, BThree
…, C 1, CTwo, CThree…… etc and M1, MTwo, MThree……But
It is set in advance. The comparison means 14 has a three-axis magnetic detection
The distance up to the side of the vessel 1
Determines whether to pass the target 15 within a predetermined range
The reference value is set in advance.

【0035】このような状態の下で、3軸磁気検出器1
に対しZ軸方向が一定のXY平面を定速度で移動する磁
性体目標物15が3軸磁気検出器1に接近し(ステップ
1)、磁気レベル検出器5が3軸磁気検出器1のX軸磁
気センサ2,Y軸磁気センサ3,Z軸磁気センサ4のい
ずれかの検出信号の中で、予め定められたレベル以上の
磁気信号を検出したとき(ステップ2)、磁気レベル検
出器5は出力を発生する。磁気レベル検出器5の出力発
生でサンプリング回路6が駆動され、サンプリング回路
6からサンプリング信号が出力する。
Under these conditions, the three-axis magnetic detector 1
In contrast, the magnetic target 15 moving at a constant speed on the XY plane in which the Z-axis direction is constant approaches the three-axis magnetic detector 1 (step 1), and the magnetic level detector 5 moves the X of the three-axis magnetic detector 1 to X. When a magnetic signal of a predetermined level or more is detected in any of the detection signals of the axis magnetic sensor 2, the Y axis magnetic sensor 3, and the Z axis magnetic sensor 4 (Step 2), the magnetic level detector 5 Generate output. The sampling circuit 6 is driven by the output of the magnetic level detector 5, and the sampling circuit 6 outputs a sampling signal.

【0036】X軸磁気センサ2,Y軸磁気センサ3及び
Z軸磁気センサ4で検出されたそれぞれの検出信号は、
対応して設けられているA/D変換器7でそれぞれディ
ジタル化され、このX軸磁気センサ2,Y軸磁気センサ
3及びZ軸磁気センサ4の検出信号は、サンプリング回
路6からのサンプリング信号発生毎に、磁気信号合成手
段8と演算手段9とに取り込まれる(ステップ3)。
Each detection signal detected by the X-axis magnetic sensor 2, the Y-axis magnetic sensor 3, and the Z-axis magnetic sensor 4 is
The detection signals of the X-axis magnetic sensor 2, the Y-axis magnetic sensor 3, and the Z-axis magnetic sensor 4 are digitized by corresponding A / D converters 7, respectively. Each time, it is taken into the magnetic signal synthesizing means 8 and the arithmetic means 9 (step 3).

【0037】磁気信号合成手段8では、各磁気センサの
検出信号Hx、Hy、Hzを基に当該磁性体目標物15
の合成磁気信号値√(Hx2+Hy2+Hz2)が演算され、
演算手段9では、取り込まれた各検出信号Hx、Hy、
Hzを基に、Z軸磁気センサ4の検出信号の絶対値|H
z|から、X軸磁気センサ2及びY軸磁気センサ3の検
出信号から得られる水平成分値√(Hx2+Hy2)を減じ
た値Hn2、すなわち、Hn2=|Hz|−√(Hx2+H
y2)の演算を行い、その演算結果の値Hn2をメモリ11
に記憶させる(ステップ4)。
In the magnetic signal synthesizing means 8, based on the detection signals Hx, Hy, and Hz of the respective magnetic sensors, the magnetic target 15
The composite magnetic signal value √ (Hx 2 + Hy 2 + Hz 2 ) of
In the arithmetic means 9, each detection signal Hx, Hy,
Hz, the absolute value of the detection signal of the Z-axis magnetic sensor 4 | H
A value Hn 2 obtained by subtracting the horizontal component value √ (Hx 2 + Hy 2 ) obtained from the detection signals of the X-axis magnetic sensor 2 and the Y-axis magnetic sensor 3 from z |, that is, Hn 2 = | Hz | −√ (Hx 2 + H
y 2 ), and stores the result Hn 2 of the result in the memory 11
(Step 4).

【0038】最大値検出手段10は、磁気信号合成手段
8から送られてくる合成磁気信号値Ht(k)=√(H
x2+Hy2+Hz2)と一つ前に取り込まれた合成磁気信号
値Ht(k−1)=√(Hx2+Hy2+Hz2)とを比較
し、合成磁気信号値√(Hx2+Hy2+Hz2)の最大値H
t0を求める処理を行う。最大値検出手段10が上記合成
磁気信号の最大値Ht0を検出したとき(ステップ5)、
最大値検出手段10は合成磁気信号の最大値Ht0を検出
した旨の信号とその合成磁気信号の最大値Ht0とを計算
手段12に送出する。計算手段12はメモリ11に記憶
された上記演算手段9による演算結果の値Hn2を、最も
新しく記憶されたものから順に前に逆上り予め定められ
た数Nだけメモリ11から順に読み出し(ステップ
6)、そのN個の値Hn2(t1 )ないしHn2(tN )を
合成磁気信号の最大値Ht0で割ってS(i)=Hn2/H
t0(i=1,2,……,N)の値を求め(ステップ
7)、そのN個の計算値S(i)を波形識別処理手段1
3に入力させる。
The maximum value detecting means 10 generates a composite magnetic signal value Ht (k) = √ (H
x 2 + Hy 2 + Hz 2 ) is compared with the previously acquired combined magnetic signal value Ht (k−1) = H (Hx 2 + Hy 2 + Hz 2 ), and the combined magnetic signal value √ (Hx 2 + Hy 2). + Hz 2 ) maximum value H
Perform processing to find t0 . When the maximum value detecting means 10 detects the maximum value Ht0 of the composite magnetic signal (step 5),
The maximum value detection means 10 sends to the calculation means 12 a signal indicating that the maximum value H t0 of the composite magnetic signal has been detected and the maximum value H t0 of the composite magnetic signal. The calculating means 12 reads the value Hn 2 of the result of the calculation by the calculating means 9 stored in the memory 11 in reverse order from the most recently stored one by a predetermined number N from the memory 11 in order (step 6). ), And dividing the N values Hn 2 (t 1 ) to Hn 2 (t N ) by the maximum value H t0 of the composite magnetic signal to obtain S (i) = Hn 2 / H
The value of t0 (i = 1, 2,..., N) is determined (step 7), and the N calculated values S (i) are obtained by the waveform identification processing means 1.
3 is input.

【0039】波形識別処理手段13では、計算手段12
で計算された上記N個の値S(t1)ないしS(tN
と予め設定されている上記重み付け値A1 ,A2 ,A3
……,B1 ,B2 ,B3 ……,C1 ,C2 ,C3 ……等
とM1 ,M2 ,M3 ……とを基に、積と和とからなるニ
ューラルネットワーク演算の波形識別処理を行う(ステ
ップ8)。
In the waveform discriminating means 13, the calculating means 12
N values S (t 1 ) to S (t N ) calculated by
And the weighting values A 1 , A 2 , A 3
.., B 1 , B 2 , B 3, ..., C 1 , C 2 , C 3, etc. and M 1 , M 2 , M 3 ,. (Step 8).

【0040】ここで、計算手段12から波形識別処理手
段13に入力される計測値S(t1)ないし値S
(tN )を基に、波形識別処理手段13がニューラルネ
ットワーク演算の波形識別処理を行うとき、3軸磁気検
出器1が海底の所定の位置に配置される前に予め学習さ
せた学習練度に基づいた学習により、どの様な波形を形
成するかの予測認識が行われ、その結果としての3軸磁
気検出器1の側方の横距離Wが無次元化された形で出力
される。
Here, the measured value S (t 1 ) to the value S inputted from the calculating means 12 to the waveform identification processing means 13 are shown.
When the waveform discriminating means 13 performs the waveform discriminating process of the neural network operation based on (t N ), the learning skill previously learned before the three-axis magnetic detector 1 is arranged at a predetermined position on the sea floor. Is performed to predict what kind of waveform is formed, and the resulting lateral distance W of the three-axis magnetic detector 1 is output in a dimensionless form.

【0041】波形識別処理手段13で処理された演算値
は、無次元化されたものであるので元の大きさに戻す演
算が実行される。すなわち波形識別処理手段13で処理
された演算値を基に、分散率η=演算値/センサの計測
レンジからHn2=η×Ht0を演算して元の大きさに戻さ
れる。
The operation value processed by the waveform discrimination processing means 13 is a dimensionless one, so that an operation for returning to the original size is executed. That is, based on the calculated value processed by the waveform identification processing means 13, Hn 2 = η × Ht0 is calculated from the dispersion ratio η = calculated value / sensor measurement range, and is returned to the original size.

【0042】一方、別の測定手段で得られる磁性体目標
物15の速度をV、計測時間をtとしたとき、図15に
示されているP点からの距離Lは、L=V×tで表され
るので、当該計測時間tから、磁性体目標物15が3軸
磁気検出器1の側方を通過する時点(タイミイング)が
得られ、その時点での3軸磁気検出器1の側方の横距離
Wが求められる。
On the other hand, when the speed of the magnetic target 15 obtained by another measuring means is V and the measuring time is t, the distance L from the point P shown in FIG. 15 is L = V × t From the measurement time t, the time (timing) at which the magnetic substance target 15 passes the side of the three-axis magnetic detector 1 is obtained, and the side of the three-axis magnetic detector 1 at that time is obtained. Is determined.

【0043】そして比較手段14で上記3軸磁気検出器
1の側方に対しどの程度の距離までの通過を、磁性体目
標物15の予め定められた範囲内の通過とするかを定め
る基準値と比較される(ステップ9)。波形識別処理手
段13で処理された上記の演算値が基準値より大きいと
き、磁性体目標物15が予め定められた範囲内の通過を
意味する「1」が比較手段14から出力され(ステップ
10)、波形識別処理手段13で処理された演算値が基
準値より小さいとき、磁性体目標物15が予め定められ
た範囲外の通過を意味する「0」が比較手段14から出
力される(ステップ11)。
Then, the comparison means 14 sets a reference value for determining how far the side of the three-axis magnetic detector 1 should pass through within a predetermined range of the magnetic target 15. Is compared with (step 9). When the calculated value processed by the waveform discrimination processing means 13 is larger than the reference value, "1" indicating that the magnetic target 15 passes within a predetermined range is output from the comparison means 14 (step 10). If the calculated value processed by the waveform identification processing means 13 is smaller than the reference value, "0" indicating that the magnetic target 15 passes outside a predetermined range is output from the comparison means 14 (step). 11).

【0044】なおステップ5で合成磁気信号値√(Hx2
+Hy2+Hz2)が最大値Ht0でないときには、当該合成
磁気信号の最大値Ht0が検出されるまでステップ3ない
しステップ5が繰り返される。
In step 5, the composite magnetic signal value √ (Hx 2
+ Hy 2 + Hz 2) is at not the maximum value H t0, the step 3 to step 5 to a maximum value H t0 of the resultant magnetic signal is detected is repeated.

【0045】また磁気レベル検出器5が3軸磁気検出器
1のX軸磁気センサ2,Y軸磁気センサ3,Z軸磁気セ
ンサ4のいずれかの検出信号でも、予め定められたレベ
ル以上の磁気信号を検出できないときには(ステップ
2)、磁性体目標物が15定点である3軸磁気検出器1
から一定範囲内を通過しないものとして処理され、終了
となる。
In addition, the magnetic level detector 5 detects a magnetic signal of a predetermined level or more even when the detection signal of any of the X-axis magnetic sensor 2, the Y-axis magnetic sensor 3, and the Z-axis magnetic sensor 4 of the three-axis magnetic detector 1 is detected. When the signal cannot be detected (step 2), the three-axis magnetic detector 1 in which the magnetic target is 15 fixed points
Is processed as if it does not pass within a certain range from, and the process ends.

【0046】上記説明では、波形識別処理手段13で処
理される上記計算値S(t1 )ないし値S(tN )の各
値が、無次元化されたものであるので元の大きさに戻す
演算を実行しているが、無次元化されたままの無次元の
下で次のようにして処理することもできる。
In the above description, since each of the calculated values S (t 1 ) to S (t N ) processed by the waveform discriminating means 13 is made dimensionless, it is reduced to the original size. Although the return operation is performed, the following operation can be performed under the dimensionless state as it is.

【0047】すなわち、計算手段12から波形識別処理
手段13に入力される計測値S(t 1 )ないし値S(t
N )を基に、波形識別処理手段13がニューラルネット
ワーク演算の波形識別処理を行うとき、3軸磁気検出器
1が海底の所定の位置に配置される前に予め学習させた
学習練度に基づいた学習により、どの様な波形を形成す
るかの予測認識が行われ、その結果としての3軸磁気検
出器1の側方の横距離Wが無次元化された形で出力され
る。
That is, the waveform discriminating process is performed by the calculating means 12.
The measured value S (t 1) Or value S (t
N), The waveform identification processing means 13
3-axis magnetic detector when performing waveform discrimination processing of work calculation
Learned before 1 was placed in place on the sea floor
What kind of waveform is formed by learning based on learning skill
Or three-axis magnetic detection as a result.
The lateral distance W on the side of the output unit 1 is output in a dimensionless form.
You.

【0048】この出力と、上記無次元化された比率で
の、上記3軸磁気検出器1の側方に対しどの程度の距離
までの通過を磁性体目標物15の予め定められた範囲内
の通過とするかを定める無次元化対応の基準値(このと
きの当該無次元化対応の基準値は無次元化された比率の
基準値を最大値Ht0で割った値)とが、比較手段14で
比較される(ステップ9)。以下同様にして、波形識別
処理手段13で処理された上記の演算値が無次元化対応
の基準値より大きいとき、磁性体目標物15が予め定め
られた範囲内の通過を意味する「1」が比較手段14か
ら出力され(ステップ10)、波形識別処理手段13で
処理された演算値が基準値より小さいとき、磁性体目標
物15が予め定められた範囲外の通過を意味する「0」
が比較手段14から出力される(ステップ11)。
The output and the distance to the side of the three-axis magnetic detector 1 at the dimensionless ratio within a predetermined range of the magnetic target 15 are determined. The reference value corresponding to the non-dimensionalization that determines whether or not to pass (the reference value corresponding to the non-dimensionalization at this time is a value obtained by dividing the reference value of the non-dimensionalization ratio by the maximum value Ht0 ) is compared with the comparison means. A comparison is made at 14 (step 9). Similarly, when the above calculated value processed by the waveform identification processing means 13 is larger than the reference value corresponding to the non-dimensionalization, "1" meaning that the magnetic substance target 15 passes through a predetermined range. Is output from the comparison means 14 (step 10), and when the operation value processed by the waveform identification processing means 13 is smaller than the reference value, "0" meaning that the magnetic target 15 passes outside a predetermined range.
Is output from the comparison means 14 (step 11).

【0049】[0049]

【発明の効果】以上説明した如く、本発明によれば、磁
性体目標物を磁気双極子モーメントで表示した上で、検
出された磁気信号を演算処理するに当たり、無次元化を
行い、磁性体目標物の大きさに無関係となし、かつ磁性
体目標物の磁気センサへの接近方向の違いの通過位置検
出への影響を消去したHn2/Ht0波形を波形識別のため
の要素となし、様々な磁性体の磁気的条件(磁気能率の
伏角及び偏角)及び深度を想定した磁気モデルでパター
ン、すなわちHn2/Ht0波形を作成し、これから求めた
重みを用い、積と和とからなるニューラルネットワーク
演算によって演算値を求め、そして基準値と比較するよ
うにしたので、磁気的条件が未知の磁性体目標物が、定
点である磁気検出器から一定範囲内を通過したかどうか
をその磁気双極子モーメントの大小にかかわらず、また
磁気的条件が未知の場合でも確実に検出できる。そして
事前の磁気測定を必要とせず、かつ磁性体目標物の磁気
的条件の影響を受けにくい通過位置検出が可能となる。
As described above, according to the present invention, a magnetic target is displayed as a magnetic dipole moment, and then the detected magnetic signal is subjected to arithmetic processing to perform dimensionless processing. An Hn 2 / Ht0 waveform, which is irrespective of the size of the target and eliminates the influence of the difference in the direction of approach of the magnetic target to the magnetic sensor on the detection of the passing position, as an element for waveform identification; create different magnetic conditions of magnetic pattern in a magnetic model assuming the (dip and declination of magnetic moment) and the depth, i.e. Hn 2 / H t0 waveform, using a weight determined from this, and a product and sum The calculated value is obtained by a neural network calculation and compared with a reference value, so that it is determined whether a magnetic target whose magnetic condition is unknown has passed within a certain range from a fixed point magnetic detector. Magnetic dipole Regardless of the moment, also the magnetic condition can be detected reliably even if unknown. Then, it is possible to detect a passing position without requiring a prior magnetic measurement and hardly affected by magnetic conditions of the magnetic target.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明に係る通過位置検出装置の一実施例構成
である。
FIG. 1 is a configuration of an embodiment of a passage position detecting device according to the present invention.

【図2】波形識別処理手段の一実施例概念説明図であ
る。
FIG. 2 is a conceptual explanatory diagram of one embodiment of a waveform identification processing means.

【図3】磁性体目標物の磁気双極子モーメントMの伏角
と偏角説明図である。
FIG. 3 is an explanatory view of a dip angle and a declination angle of a magnetic dipole moment M of a magnetic target.

【図4】磁性体目標物と3軸磁気検出器との直交座標系
(X,Y,Z)における配置関係説明図である。
FIG. 4 is an explanatory diagram of a positional relationship between a magnetic target and a three-axis magnetic detector in a rectangular coordinate system (X, Y, Z).

【図5】計算によってシミュレートしたHn2/Ht0波形
図である。
FIG. 5 is a Hn 2 / Ht0 waveform diagram simulated by calculation.

【図6】条件を変え計算によってシミュレートしたHn2
/Ht0波形図である。
FIG. 6: Hn 2 simulated by calculation with changing conditions
FIG. 14 is a / Ht0 waveform diagram.

【図7】条件を変え計算によってシミュレートしたHn2
/Ht0波形図である。
FIG. 7: Hn 2 simulated by calculation with changing conditions
FIG. 14 is a / Ht0 waveform diagram.

【図8】条件を変え計算によってシミュレートしたHn2
/Ht0波形図である。
FIG. 8: Hn 2 simulated by calculation with changing conditions
FIG. 14 is a / Ht0 waveform diagram.

【図9】条件を変え計算によってシミュレートしたHn2
/Ht0波形図である。
FIG. 9: Hn 2 simulated by calculation with changing conditions
FIG. 14 is a / Ht0 waveform diagram.

【図10】条件を変え計算によってシミュレートしたH
n2/Ht0波形図である。
FIG. 10 shows H simulated by calculation with changing conditions.
It is n 2 / H t0 waveform.

【図11】条件を変え計算によってシミュレートしたH
n2/Ht0波形図である。
FIG. 11 is a graph simulating H by changing conditions.
It is n 2 / H t0 waveform.

【図12】条件を変え計算によってシミュレートしたH
n2/Ht0波形図である。
FIG. 12 is a graph simulating H by changing conditions.
It is n 2 / H t0 waveform.

【図13】条件を変え計算によってシミュレートしたH
n2/Ht0波形図である。
FIG. 13 is a graph simulating H by changing conditions and calculating.
It is n 2 / H t0 waveform.

【図14】条件を変え計算によってシミュレートしたH
n2/Ht0波形図である。
FIG. 14 is a graph simulating H by changing conditions.
It is n 2 / H t0 waveform.

【図15】X−Y平面での磁性体目標物の進行方向と座
標を説明しているX−Y平面説明図である。
FIG. 15 is an XY plane explanatory view illustrating a traveling direction and coordinates of a magnetic target on the XY plane.

【図16】本発明の一実施例フローチャートである。FIG. 16 is a flowchart of one embodiment of the present invention.

【図17】磁性体目標物の進行方向と船体磁気との磁気
能率についてのZ軸方向成分説明図である。
FIG. 17 is a diagram illustrating a Z-axis direction component of the magnetic efficiency between the traveling direction of the magnetic target and the hull magnetism.

【符号の説明】[Explanation of symbols]

1 3軸磁気検出器 2 X軸磁気センサ 3 Y軸磁気センサ 4 Z軸磁気センサ 5 磁気レベル検出器 6 サンプリング回路 7 A/D変換器 8 磁気信号合成手段 9 演算手段 10 最大値検出手段 11 メモリ 12 計算手段 13 波形識別処理手段 14 比較手段 15 磁性体目標物 Reference Signs List 1 3-axis magnetic detector 2 X-axis magnetic sensor 3 Y-axis magnetic sensor 4 Z-axis magnetic sensor 5 Magnetic level detector 6 Sampling circuit 7 A / D converter 8 Magnetic signal synthesizing means 9 Computing means 10 Maximum value detecting means 11 Memory 12 Calculation means 13 Waveform identification processing means 14 Comparison means 15 Magnetic target

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI // G01R 33/02 G01R 33/02 Q (72)発明者 吉見 邦彦 神奈川県横浜市鶴見区平安町2丁目29番 地の1 株式会社 京三製作所内 (72)発明者 森下 晴之 神奈川県横浜市鶴見区平安町2丁目29番 地の1 株式会社 京三製作所内 (56)参考文献 特開 平2−74892(JP,A) 特開 昭55−162078(JP,A) 特開 昭60−71910(JP,A) (58)調査した分野(Int.Cl.7,DB名) G01B 7/00 G01C 15/00 G01V 3/08 G01V 3/10 G01R 33/02 G01C 21/00 ────────────────────────────────────────────────── ─── Continuing on the front page (51) Int.Cl. 7 Identification symbol FI // G01R 33/02 G01R 33/02 Q (72) Inventor Kunihiko Yoshimi 2-29 Heiancho, Tsurumi-ku, Yokohama-shi, Kanagawa 1 Kyosan Seisakusho Co., Ltd. (72) Inventor Haruyuki Morishita 2-29-2 Heian-cho, Tsurumi-ku, Yokohama-shi, Kanagawa 1 Kyosan Seisakusho Co., Ltd. (56) References JP-A-2-74892 (JP, A JP-A-56-162078 (JP, A) JP-A-60-71910 (JP, A) (58) Fields investigated (Int. Cl. 7 , DB name) G01B 7/00 G01C 15/00 G01V 3 / 08 G01V 3/10 G01R 33/02 G01C 21/00

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 Z軸を垂直軸とする直交座標系(X,
Y,Z)の各軸上に指向性を有する磁気センサがそれぞ
れ配置されてなるその位置が固定の3軸磁気検出器を備
え、当該3軸磁気検出器に対しZ軸方向が一定のXY平
面を移動する磁性体目標物の予め定められた範囲内の通
過を検出する通過位置検出方法において、 3軸磁気検出器のサンプリングされた各検出信号から磁
性体目標物の合成磁気信号値を求めると共に、その最大
値となるタイミングを求め、 サンプリング信号毎に、Z軸磁気センサの検出信号の絶
対値|Hz|から、X軸磁気センサ及びY軸磁気センサ
の検出信号から得られる水平成分値√(Hx2+Hy2)を
減じた値Hn2(i)=|Hz|−√(Hx2+Hy2)(i
=1,2,……,n)を演算してその演算結果を記憶手
段に記憶し、 磁性体目標物の合成磁気信号値が最大となったとき、記
憶手段に最も新しく記憶されたものから順に前に逆上り
記憶手段に記憶されている予め定められた数の各値Hn2
を合成磁気信号の最大値Ht0でそれぞれ割った値S
(i)=Hn2/Ht0を求め、 この各値S(i)=Hn2/Ht0と、予め定められた重み
付け値とを基にニューラルネットワーク演算処理を行っ
て波形識別を行い、 波形識別の演算値から、磁性体目標物の予め定められた
範囲内の通過を検出するようにしたことを特徴とする通
過位置検出方法。
1. An orthogonal coordinate system (X,
A magnetic sensor having directivity on each axis of (Y, Z) is provided with a three-axis magnetic detector whose position is fixed, and an XY plane having a fixed Z-axis direction with respect to the three-axis magnetic detector; A passage position detecting method for detecting passage of a moving magnetic target within a predetermined range, wherein a combined magnetic signal value of the magnetic target is obtained from each sampled detection signal of the three-axis magnetic detector. , The timing at which the maximum value is obtained, and the horizontal component value √ () obtained from the detection signals of the X-axis magnetic sensor and the Y-axis magnetic sensor from the absolute value | Hz | of the detection signal of the Z-axis magnetic sensor for each sampling signal. Hx 2 + Hy 2 ) Hn 2 (i) = | Hz | − (Hx 2 + Hy 2 ) (i
= 1, 2,..., N) and stores the result of the calculation in the storage means. When the composite magnetic signal value of the magnetic target reaches the maximum value, the storage means starts with the most recently stored value. A predetermined number of values Hn 2 previously stored in the reverse storage means in order.
Are divided by the maximum value H t0 of the composite magnetic signal, respectively.
(I) = seeking Hn 2 / H t0, the each value S (i) = Hn 2 / H t0, performs waveform identification by performing a neural network processing based on the weighting value set in advance, the waveform A passing position detecting method, wherein a passing of a magnetic target within a predetermined range is detected from a calculated value of the identification.
【請求項2】 Z軸を垂直軸とする直交座標系(X,
Y,Z)の各軸上に指向性を有する磁気センサがそれぞ
れ配置されてなるその位置が固定の3軸磁気検出器を備
え、当該3軸磁気検出器に対しZ軸方向が一定のXY平
面を移動する磁性体目標物の予め定められた範囲内の通
過を検出する通過位置検出装置において、 予め定められたレベル以上の磁気信号を検出したとき出
力を発生する磁気レベル検出器と、 磁気レベル検出器の出力で3軸磁気検出器の各検出信号
をサンプリングさせるサンプリング回路と、 サンプリング回路が出力するサンプリング信号毎に、3
軸磁気検出器の各検出信号から磁性体目標物の合成磁気
信号を求める磁気信号合成手段と、 磁気信号合成手段で求められた合成磁気信号を基に、当
該合成磁気信号の最大値Ht0を検出する最大値検出手段
と、 サンプリング回路が出力するサンプリング信号毎に、Z
軸磁気センサの検出信号の絶対値|Hz|から、X軸磁
気センサ及びY軸磁気センサの検出信号から得られる水
平成分値√(Hx2+Hy2)を減じた値Hn2(i)=|H
z|−√(Hx2+Hy2)(i=1,2,……,n)を演
算する演算手段と、 演算手段の演算結果を記憶する記憶手段と、 最大値検出手段が上記合成磁気信号の最大値Ht0を検出
したとき、当該記憶手段に最も新しく記憶されたものか
ら順に前に逆上り記憶されている予め定められた数の各
値Hn2を読み出し、この値Hn2を最大値検出手段で検出
された合成磁気信号の最大値Ht0で割る計算Hn2/Ht0
を実行する計算手段と、 計算手段で計算された各計算値S(i)=Hn2/Ht0
予め定められた重み付け値とを基にニューラルネットワ
ーク演算処理を行い、波形識別を行う波形識別処理手段
と、 波形識別処理手段が処理した演算値と予め定められた基
準値とを比較する比較手段とを備えたことを特徴とする
通過位置検出装置。
2. An orthogonal coordinate system (X,
A magnetic sensor having directivity on each axis of (Y, Z) is provided with a fixed three-axis magnetic detector at a fixed position, and an XY plane having a fixed Z-axis direction with respect to the three-axis magnetic detector; A passage level detection device that detects passage of a moving magnetic target within a predetermined range, a magnetic level detector that generates an output when a magnetic signal of a predetermined level or more is detected, and a magnetic level. A sampling circuit for sampling each detection signal of the three-axis magnetic detector at the output of the detector;
Magnetic signal synthesizing means for obtaining a synthetic magnetic signal of a magnetic target from each detection signal of the axial magnetism detector; and a maximum value Ht0 of the synthetic magnetic signal based on the synthetic magnetic signal obtained by the magnetic signal synthesizing means. A maximum value detecting means for detecting, and a Z signal for each sampling signal output from the sampling circuit.
Hn 2 (i) = | obtained by subtracting the horizontal component value √ (Hx 2 + Hy 2 ) obtained from the detection signals of the X-axis magnetic sensor and the Y-axis magnetic sensor from the absolute value | Hz | H
z |-を (Hx 2 + Hy 2 ) (i = 1, 2,..., n); a storage unit for storing the calculation result of the calculation unit; upon detection of the maximum value H t0, reads out the storage unit most recently stored values Hn number of predetermined which are sequentially before the reverse uplink storage from those in 2, the maximum value the value Hn 2 Calculation Hn 2 / H t0 divided by the maximum value H t0 of the composite magnetic signal detected by the detection means
Means for performing a neural network operation process based on each of the calculated values S (i) = Hn 2 / Ht0 and a predetermined weight value, and performing waveform identification for performing waveform identification. A passing position detecting device comprising: a processing unit; and a comparing unit that compares a calculated value processed by the waveform identification processing unit with a predetermined reference value.
JP35059297A 1997-12-19 1997-12-19 Passage position detection method and passage position detection device Expired - Lifetime JP3032776B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP35059297A JP3032776B2 (en) 1997-12-19 1997-12-19 Passage position detection method and passage position detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP35059297A JP3032776B2 (en) 1997-12-19 1997-12-19 Passage position detection method and passage position detection device

Publications (2)

Publication Number Publication Date
JPH11183109A JPH11183109A (en) 1999-07-09
JP3032776B2 true JP3032776B2 (en) 2000-04-17

Family

ID=18411534

Family Applications (1)

Application Number Title Priority Date Filing Date
JP35059297A Expired - Lifetime JP3032776B2 (en) 1997-12-19 1997-12-19 Passage position detection method and passage position detection device

Country Status (1)

Country Link
JP (1) JP3032776B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100555668B1 (en) * 2003-11-11 2006-03-03 삼성전자주식회사 Fluxgate sensor for calibrating azimuth at slope and calibration method thereof
JP4631622B2 (en) * 2005-09-05 2011-02-16 株式会社島津製作所 Magnetic measurement system
JP4958424B2 (en) * 2005-11-07 2012-06-20 ユニバーサル特機株式会社 Electric field detection method, apparatus, electric field detection method program, moving object position estimation detection method, apparatus, and moving object position estimation detection method program
JP5109344B2 (en) * 2005-11-17 2012-12-26 株式会社島津製作所 Magnetic measurement system
JP5151540B2 (en) * 2008-02-21 2013-02-27 株式会社島津製作所 Target probe
JP5228555B2 (en) * 2008-03-19 2013-07-03 株式会社島津製作所 Target body exploration device
JP6766333B2 (en) * 2015-10-06 2020-10-14 愛知製鋼株式会社 Micromagnetic material detection sensor and foreign matter detection device
JP7435036B2 (en) 2020-03-03 2024-02-21 株式会社島津製作所 magnetic detection system

Also Published As

Publication number Publication date
JPH11183109A (en) 1999-07-09

Similar Documents

Publication Publication Date Title
US3644825A (en) Magnetic detection system for detecting movement of an object utilizing signals derived from two orthogonal pickup coils
US6476610B1 (en) Magnetic anomaly sensing system and methods for maneuverable sensing platforms
CA2311323C (en) Magnetic object tracking based on direct observation of magnetic sensor measurements
JP2012528304A (en) Method and apparatus for determining position by magnetism
JP3032776B2 (en) Passage position detection method and passage position detection device
CN110221302A (en) Environmental detection device and its modification method, system, portable equipment and storage medium
US7577544B2 (en) Method and system for improving target localization and characterization
JP2520952B2 (en) Vehicle direction detector
JP2009053056A (en) Magnetic moving body speed detector
JPS60135814A (en) Azimuth detecting apparatus
CA1326708C (en) Magnetic system for locating avalanche victims and other lost objects
JP3007935B2 (en) Passage position detection method and passage position detection device
JPS5830557B2 (en) How to detect magnetic substances
Wiegert et al. Magnetic anomaly sensing system for mine countermeasures using high mobility autonomous sensing platforms
JPH07198407A (en) Magnetism measuring instrument
Gutnik et al. Data-Driven Underwater Navigation workshop: AUV Close-Range Localization and Guidance Employing an Electro-Magnetic Beacon
JPH0778540B2 (en) Magnetic detection device capable of position localization
Janošek et al. Compact magnetic gradiometer and its astatization
JPH05274589A (en) Vehicle detector and traffic volume measuring instrument
JPH06325288A (en) Vehicle detecting device and traffic quantity mesuring equipment
KR20170092356A (en) System for calibrating azimuth of 3-axis magnetic sensor
JP2961945B2 (en) Magnetometer
JPH07174831A (en) Magnetism measuring apparatus
JP2542673Y2 (en) Simulated magnetic field signal generator for adjustment and inspection of object detector
JPH0650248B2 (en) Mobile body orientation detector

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

Free format text: PAYMENT UNTIL: 20080218

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20090218

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20090218

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20100218

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20100218

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20110218

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20120218

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20120218

Year of fee payment: 12

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

Free format text: PAYMENT UNTIL: 20130218

Year of fee payment: 13

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

Free format text: PAYMENT UNTIL: 20140218

Year of fee payment: 14

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

Free format text: PAYMENT UNTIL: 20140218

Year of fee payment: 14

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

EXPY Cancellation because of completion of term