JP4815859B2 - Magnetic model calculation method - Google Patents
Magnetic model calculation method Download PDFInfo
- Publication number
- JP4815859B2 JP4815859B2 JP2005130776A JP2005130776A JP4815859B2 JP 4815859 B2 JP4815859 B2 JP 4815859B2 JP 2005130776 A JP2005130776 A JP 2005130776A JP 2005130776 A JP2005130776 A JP 2005130776A JP 4815859 B2 JP4815859 B2 JP 4815859B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic
- data
- actual measurement
- virtual
- measured
- 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 - Fee Related
Links
Images
Landscapes
- Measuring Magnetic Variables (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Description
この発明は，磁性構造体の外部磁界測定を行い、その実測値をもとに、磁性構造体の磁気モデルを計算する磁気モデル計算方法に関する。 The present invention relates to a magnetic model calculation method for measuring an external magnetic field of a magnetic structure and calculating a magnetic model of the magnetic structure based on the actually measured value.
船舶の磁性体からなる船体のＸ，Ｙ，Ｚ方向の外部磁界を打ち消すために、船体内に複数個の消磁コイルを設けると共にその船体内に複数個の磁気検知器からなる磁気監視部を設置し、各磁気検知器から測定された船内磁界に基づいて算出した船外磁気モーメントと、予め測定算出した各消滅コイル効果による船外磁気モーメントとから、外部磁界を最小にする消磁電流を決定して、各消磁コイルに通電して、消磁する方法が開示されている（例えば特許文献１参照）。 In order to cancel the external magnetic field in the X, Y, and Z directions of the ship's magnetic body, a plurality of degaussing coils are installed inside the ship and a magnetic monitoring unit consisting of a plurality of magnetic detectors is installed in the ship's body. The demagnetizing current that minimizes the external magnetic field is determined from the outboard magnetic moment calculated based on the inboard magnetic field measured from each magnetic detector and the outboard magnetic moment due to each extinction coil effect measured and calculated in advance. Thus, a method of demagnetizing each demagnetizing coil by energizing each demagnetizing coil is disclosed (for example, see Patent Document 1).
従来、この種の艦艇などの消磁を行うために、実測した外部磁界から、任意の位置の磁界を推定する方法として、（１）磁界が距離比率のべき乗に比例して減衰・増大するという理論を用いた「べき乗計算方式」、（２）外部磁界から逆計算で磁気源を求め任意の点の磁界を再計算する方式（「マルチダイポール方式」、「長球調和関数方式」）などがある。
一例として、マルチダイポール方式について説明する。今磁気構造体である艦艇が、図６に示すように位置ｎ＝１〜６に配置される６個のダイポールモーメントＭ１、Ｍ２、・・・、Ｍ６で表せる磁気源であるとする。このうちの１つのダイポールモーメントＭにつき図７に示すようにｒ離れた点Ｐ（ｘ，ｙ，ｚ）の磁界は次式で表せる。
Conventionally, in order to demagnetize this type of ship, etc., as a method of estimating the magnetic field at an arbitrary position from the actually measured external magnetic field, (1) The theory that the magnetic field attenuates and increases in proportion to the power of the distance ratio (2) A method of calculating a magnetic source by inverse calculation from an external magnetic field and recalculating a magnetic field at an arbitrary point (“multi-dipole method”, “long spherical harmonic function method”), etc. .
As an example, a multi-dipole system will be described. Assume that the ship, which is a magnetic structure, is a magnetic source represented by six dipole moments M1, M2,..., M6 arranged at positions n = 1 to 6 as shown in FIG. As shown in FIG. 7, the magnetic field at a point P (x, y, z) separated by r per one dipole moment M can be expressed by the following equation.
Ｈｘ＝１００×〔（３ｘ２‐ｒ２）・Ｍｘ＋３ｘｙ・Ｍｙ＋３ｘｚ・Ｍｚ〕／ｒ５〔ｎＴ〕
Ｈｙ＝１００×〔３ｘｙ・Ｍｘ＋（３ｙ２‐ｒ２）・Ｍｙ＋３ｙｚ・Ｍｚ〕／ｒ５〔ｎＴ〕
Ｈｚ＝１００×〔３ｘｚ・Ｍｘ＋３ｙｚ・Ｍｙ＋（３ｚ２‐ｒ２）・Ｍｚ〕／ｒ５〔ｎＴ〕
Ｍｘ：ダイポールモーメントのＸ成分
Ｍｙ：ダイポールモーメントのＹ成分
Ｍｚ：ダイポールモーメントのＺ成分
ｒ：ダイポールモーメントと観測点Ｐとの距離 ｒ＝（ｘ２＋ｙ２＋ｚ２）０．５
ダイポールモーメントＭ及び、その各成文Ｍｘ、Ｍｙ、Ｍｚと観測点Ｐとの関係を図７に示す。
Hx = 100 × [(3x 2 −r 2 ) · Mx + 3xy · My + 3xz · Mz] / r 5 [nT]
Hy = 100 × [3xy · Mx + (3y 2 −r 2 ) · My + 3yz · Mz] / r 5 [nT]
Hz = 100 × [3xz · Mx + 3yz · My + (3z 2 −r 2 ) · Mz] / r 5 [nT]
Mx: X component of dipole moment My: Y component of dipole moment Mz: Z component of dipole moment r: Distance between dipole moment and observation point P r = (x 2 + y 2 + z 2 ) 0.5
FIG. 7 shows the relationship between the dipole moment M and its respective sentences Mx, My, Mz and the observation point P.
マルチダイポール方式とは、ダイポールモーメントが艦艇内に複数個分布していると仮定し、それぞれのダイポールモーメントからの磁界の積算で磁界を近似するという考え方である。例えば、図６のようにダイポールの個数がｎ＝６個の場合、その磁界は次の式で与えられる。
Ｈｘ（全体）＝１００×〔（３ｘ１ ２‐ｒ１ ２）・Ｍｘ１＋３ｘ１ｙ１・Ｍｙ１＋３ｘ１ｚ１・Ｍｚ１〕／ｒ１ ５
＋１００×〔（３ｘ２ ２‐ｒ２ ２）・Ｍｘ２＋３ｘ２ｙ２・Ｍｙ２＋３ｘ２ｚ２・Ｍｚ２〕／ｒ２ ５
＋・・・・・・・
＋１００×〔（３ｘ６ ２‐ｒ６ ２）・Ｍｘ６＋３ｘ６ｙ６・Ｍｙ６＋３ｘ６ｚ６・Ｍｚ６〕／ｒ６ ５〔ｎＴ〕
Ｈｙ（全体）＝１００×〔３ｘ１ｙ１・Ｍｘ１＋（３ｙ１ ２‐ｒ１ ２）・Ｍｙ１＋３ｙ１ｚ１・Ｍｚ１〕／ｒ１ ５
＋１００×〔３ｘ２ｙ２・Ｍｘ２＋（３ｙ２ ２‐ｒ２ ２）・Ｍｙ２＋３ｙ２ｚ２・Ｍｚ２〕／ｒ２ ５
＋・・・・・・・
＋１００×〔３ｘ６ｙ６・Ｍｘ６＋（３ｙ６ ２‐ｒ６ ２）・Ｍｙ６＋３ｙ６ｚ６・Ｍｚ６〕／ｒ６ ５〔ｎＴ〕
Ｈｚ（全体）＝１００×〔３ｘ１ｚ１・Ｍｘ１＋３ｙ１ｚ１・Ｍｙ１＋（３ｚ１ ２‐ｒ１ ２）・Ｍｚ１〕／ｒ１ ５
＋１００×〔３ｘ２ｚ２・Ｍｘ２＋３ｙ２ｚ２・Ｍｙ２＋（３ｚ２ ２‐ｒ２ ２）・Ｍｚ２〕／ｒ２ ５
＋・・・・・・・
＋１００×〔３ｘ６ｚ６・Ｍｘ６＋３ｙ６ｚ６・Ｍｙ６＋（３ｚ６ ２‐ｒ６ ２）・Ｍｚ６〕／ｒ６ ５〔ｎＴ〕
各ダイポールモーメントの３軸成分を各配置した複数の磁気センサの実測磁界から最小自乗法等で逆算して求め、求められた各ダイポールモーメントから、任意の水深の磁界を逆計算する。
Hx (whole) = 100 × [(3x 1 2 −r 1 2 ) · Mx 1 + 3x 1 y 1 · My 1 + 3x 1 z 1 · Mz 1 ] / r 1 5
+ 100 × [(3x 2 2 -r 2 2 ) · Mx 2 + 3x 2 y 2 · My 2 + 3x 2 z 2 · Mz 2 ] / r 2 5
+ ...
+ 100 × [(3x 6 2 -r 6 2 ) · Mx 6 + 3x 6 y 6 · My 6 + 3x 6 z 6 · Mz 6 ] / r 6 5 [nT]
Hy (whole) = 100 × [3x 1 y 1 · Mx 1 + (3y 1 2 −r 1 2 ) · My 1 + 3y 1 z 1 · Mz 1 ] / r 1 5
+ 100 × [3x 2 y 2 · Mx 2 + (3y 2 2 -r 2 2 ) · My 2 + 3y 2 z 2 · Mz 2 ] / r 2 5
+ ...
+ 100 × [3x 6 y 6 · Mx 6 + (3y 6 2 -r 6 2 ) · My 6 + 3y 6 z 6 · Mz 6 ] / r 6 5 [nT]
Hz (overall) = 100 × [3x 1 z 1 · Mx 1 + 3y 1 z 1 · My 1 + (3z 1 2 −r 1 2 ) · Mz 1 ] / r 1 5
+ 100 × [3x 2 z 2 · Mx 2 + 3y 2 z 2 · My 2 + (3z 2 2 -r 2 2 ) · Mz 2 ] / r 2 5
+ ...
+ 100 × [3x 6 z 6 · Mx 6 + 3y 6 z 6 · My 6 + (3z 6 2 -r 6 2 ) · Mz 6 ] / r 6 5 [nT]
A triaxial component of each dipole moment is obtained by back-calculating from the measured magnetic fields of a plurality of magnetic sensors arranged by the least square method or the like, and a magnetic field at an arbitrary water depth is back-calculated from each obtained dipole moment.
上記した従来の実測した外部磁界から、任意の位置の磁界を推定する方法において、磁気源を求めた後に実測データ位置（磁気センサの配置位置）の磁界を再計算（あてはめ計算という）を行う際、実測データの存在する位置については、実測データとほぼ同等の計算値が得られるが、実測データのない位置（磁気センサの配置されない位置）についての計算値は異常な値を計算結果として得ることがある。これは最小自乗法において実測データだけとの誤差を最小としてしまうような磁気源の計算をしてしまうため磁気モーメント自体が角度、大きさにおいて異常な値が計算されることによる。この現象は「発散」あるいは「食いつぶし現象」と呼ばれるものであり、特に位置測定方式では、船底下の磁気センサ数がコスト面などから限られるため、実測データのない位置では、等磁界曲線で実際にはありえない磁界が計算されることがある。 In the above-described conventional method for estimating the magnetic field at an arbitrary position from the actually measured external magnetic field, the magnetic field at the measured data position (positioning position of the magnetic sensor) is recalculated (referred to as fitting calculation) after obtaining the magnetic source. For the position where the actual measurement data exists, a calculated value almost equivalent to the actual measurement data can be obtained, but the calculated value for the position without the actual measurement data (position where the magnetic sensor is not arranged) should be obtained as an abnormal result. There is. This is because, in the least square method, the magnetic source is calculated so as to minimize the error with only the actually measured data, and therefore the magnetic moment itself is calculated as an abnormal value in angle and size. This phenomenon is called “divergence” or “crushing phenomenon”. In particular, in the position measurement method, the number of magnetic sensors under the ship bottom is limited in terms of cost and so on. A magnetic field that is not possible is calculated.
この発明は上記問題点に着目してなされたものであって、最小自乗法で得られた正規
方程式の解、すなわち磁気源の発散を抑えることができ、より精度の良い近傍磁界計算値や等磁界曲線を得ることができる磁気モデル計算方法を提供することを目的とする。
The present invention has been made paying attention to the above problems, and can solve the normal equation obtained by the method of least squares, that is, suppress the divergence of the magnetic source, and can calculate a near magnetic field calculation value with higher accuracy, etc. It is an object to provide a magnetic model calculation method capable of obtaining a magnetic field curve.
この発明は、磁性構造体に対して複数の外部磁界を実測し、この複数の実測データに基づいて磁性構造体における仮想的な磁気源を計算し、実測データを測定した実測位置以外の任意の位置の磁界を計算する磁気測定システムにおける磁気モデル計算方法において、前記実測位置に対して実測データのない仮想位置を設定し、この仮想位置における仮想データを、仮想位置の設定に係る複数の実測位置における複数の実測データに基づく補間により求め、前記実測位置の計算上のデータが当該実測位置の実測データと同じになるように前記実測データ及び前記仮想データを正規方程式に代入し、最小自乗法による前記磁気源の計算を行うものである。 According to the present invention , a plurality of external magnetic fields are actually measured for a magnetic structure, a virtual magnetic source in the magnetic structure is calculated based on the plurality of measured data, and an arbitrary position other than the actually measured position where the measured data is measured is calculated . In a magnetic model calculation method in a magnetic measurement system for calculating a magnetic field at a position, a virtual position without actual measurement data is set for the actual measurement position, and the virtual data at the virtual position is used as a plurality of actual measurement positions related to the setting of the virtual position. Substituting the actual measurement data and the virtual data into a normal equation so that the calculation data of the actual measurement position is the same as the actual measurement data of the actual measurement position, and by the least square method The calculation of the magnetic source is performed.
仮想データを得る具体例としては、例えば、キール下の磁界を実測データ、左右舷遠方磁界を零とし、その間のデータをスプライン関数等で近似し、仮定する方法を採用する。 As a specific example of obtaining the virtual data, for example, a method of assuming that the magnetic field under the keel is actually measured data, the left and right far field is zero, and the data between them is approximated by a spline function or the like is used.
また、実測データに重みつけするために、仮想データは、実測データの数より少ない個数の値に設定する。 Further, in order to weight the actually measured data, the virtual data is set to a number of values smaller than the number of actually measured data.
また、この発明において、磁気源の計算がマルチダイポール方式である場合、前記磁性構造体の磁性体量から想定され得る最大磁気量を規定値として設定し、前記磁気源の磁気モーメントの磁気量が前記規定値以下となるまで、規定値の再設定、前記実測データ及び前記仮想データの正規方程式への代入並びに最小自乗法による前記磁気源の計算を繰り返すと良い。 In the present invention , when the calculation of the magnetic source is a multi-dipole method, the maximum magnetic amount that can be assumed from the amount of magnetic material of the magnetic structure is set as a specified value, and the magnetic amount of the magnetic moment of the magnetic source is It is preferable to repeat the resetting of the specified value, the substitution of the actual measurement data and the virtual data into the normal equation, and the calculation of the magnetic source by the method of least squares until the value becomes the specified value or less .
請求項１に係る発明によれば、実測位置の他に仮想位置を設定し、この仮想位置の仮想データを仮想位置の設定に係る複数の実測位置における複数の実測データに基づく補間により求め、実測データと仮想データにより磁気モデル計算を行うので、磁気源の発散を抑えることができる。従ってより精度の高い近傍磁界計算値や等磁界曲線を得ることができる。
また、請求項２に係る発明によれば、計算される磁気源の磁気モーメントが均一化し、磁気源の発散を防ぐことができる。
According to the first aspect of the present invention, a virtual position is set in addition to the actually measured position, and virtual data of this virtual position is obtained by interpolation based on a plurality of actually measured data at a plurality of actually measured positions related to the setting of the virtual position. Since magnetic model calculation is performed using data and virtual data , the divergence of the magnetic source can be suppressed. Accordingly, it is possible to obtain a near magnetic field calculation value and an isomagnetic field curve with higher accuracy.
Moreover, according to the invention which concerns on
以下、実施の形態により、この発明をさらに詳細に説明する。図１の（ａ）は、艦艇１の発生磁界を計測し、磁気モデルを計算し、磁気モデルを計算するための、磁気計測システムの設置例を示す概略図である。図１の（ａ）において海底下に直線上に複数個の磁気センサＭＤ１、ＭＤ２、・・・、ＭＤ１６が固定的に設置されている。この磁気センサＭＤ１、ＭＤ２、・・・、ＭＤ１６の設置された領域の海面に、被測定艦艇１が係留されている。
Hereinafter, the present invention will be described in more detail with reference to embodiments. FIG. 1A is a schematic diagram illustrating an installation example of a magnetic measurement system for measuring a magnetic field generated by a
図２は、この実施形態磁気計測システムの構成を示すブロック図である。図２において、磁気計測装置１１は、制御、演算の種々の処理を実行するＣＰＵ１２と、実測データ、仮想データ、設定値等を記憶する各記憶部を存するメモリ１３と、磁気センサＭＤ１、ＭＤ２、・・・、ＭＤＮ（Ｎ＝１６）に接続するための入力ポート１４とを備えている。磁気計測装置１１は、艦艇１の係留地点に近い港湾の施設に設置されるが、係留された艦艇１上に搭載してもよい。
FIG. 2 is a block diagram showing the configuration of the magnetic measurement system of this embodiment. In FIG. 2, a magnetic measurement device 11 includes a
次に、この実施形態磁気測定システムにおいて、各磁気センサＭＤ１、・・・ＭＤ１６による実測データに基づき、磁気モデル計算を行う場合の処理を、図３に示すフロー図を参照して説明する。 Next, in the magnetic measurement system of this embodiment, the processing in the case of performing a magnetic model calculation based on the measured data by each magnetic sensor MD1,... MD16 will be described with reference to the flowchart shown in FIG.
この処理動作が開始されると、ステップＳＴ１において各磁気センサＭＤ１、・・・ＭＤ１６から各点におけるＸ、Ｙ、Ｚの各成分磁界が検知される。つまり、外部磁界が測定され、各Ｘ、Ｙ、Ｚの成分磁界がメモリ１３の実測データ記憶部に記憶される。 When this process operation is started, the magnetic sensors MD1 at step ST1, X from · · · MD16 at each point, Y, each component magnetic field Z is Ru is detected. That is, the external magnetic field is measured, and the component magnetic fields of X, Y, and Z are stored in the actual measurement data storage unit of the memory 13.
この場合、各磁気センサＭＤ１，・・・ＭＤ１６には、それぞれ図１の（ｂ）に示す対応位置の各成分が検知される。ここでＸ成分は艦艇１の船首尾方向の磁界成分であり、船首方向が＋である。Ｚ成分は艦艇１の垂直方向の磁界成分であり、海底方向が＋である。Ｙ成分は左右方向の磁界成分であり、左舷方向が＋である。図１で艦艇１のキール中央下に近い磁気センサＭＤ８は、Ｚ成分が最大値となり、Ｘ成分は、０に近い値となる。続いてステップＳＴ２へ移行する。
In this case, each magnetic sensor MD1,... MD16 detects each component at the corresponding position shown in FIG. Here, the X component is a magnetic field component in the bow-stern direction of the
ステップＳＴ２においては、実測データのない位置の仮想データの設定を行う。図１の磁気センサは、実際に１直線上に、１６個の磁気センサＭＤ１、・・・ＭＤ１６が配置されている。ここでは仮想データを設定する位置として、磁気センサＭＤ８，ＭＤ６、ＭＤ１０の左右方向に一定距離離れた位置Ｐ８a、Ｐ８ｂ、Ｐ６a、Ｐ６ｂ、Ｐ１０a、Ｐ１０ｂとする。 In step ST2, virtual data at a position without actual measurement data is set. In the magnetic sensor shown in FIG. 1, 16 magnetic sensors MD1,... MD16 are actually arranged on one straight line. Here, the positions where virtual data are set are positions P8a, P8b, P6a, P6b, P10a, and P10b that are a fixed distance apart in the left-right direction of the magnetic sensors MD8, MD6, and MD10.
これらの位置における仮想データの設定方法を、磁気センサＭＤ８の左右に設ける位置Ｐ８a、Ｐ８bを例にして説明する。図４の矢印Ａが艦艇１の船首方向とし、Ｐ８が磁気センサＭＤ８の位置、Ｐ８aが左側の仮想データを設定する位置、Ｐ８bが右側の仮想データを設定する位置である。曲線ｍは艦艇１のキール下における磁界の強さのＺ成分を示す波形であり、磁気センサＭＤ８の位置Ｐ８でＺ成分磁界が最大値である。これに対し、左右舷方向の遠方は、指数関数的に磁界が減衰するので、零に収束する。
A method for setting virtual data at these positions will be described using positions P8a and P8b provided on the left and right of the magnetic sensor MD8 as an example. 4 is the bow direction of the
そこで、位置Ｐ８の磁界（Ｘ成分、Ｙ成分、Ｚ成分）を最大とし左右舷方向の遠方の磁界が零になるように例えばスプライン関数を用いて近似計算し、仮想位置Ｐ８a、Ｐ８bのデータを求める。図４においてaは仮想位置Ｐ８aにおけるＺ成分波形、bは仮想位置Ｐ８bにおけるＺ成分波形である。なおここで、使用するスプライン補間式は次式で与えられる。 Therefore, an approximate calculation is performed using, for example, a spline function so that the magnetic field at the position P8 (X component, Y component, Z component) is maximized, and the far magnetic field in the left-right direction becomes zero, and the data at the virtual positions P8a and P8b are obtained. Ask. In FIG. 4, a is a Z component waveform at the virtual position P8a, and b is a Z component waveform at the virtual position P8b. Here, the spline interpolation formula to be used is given by the following formula.
Ｓｊ（ｘ）＝ａｊ（ｘ-ｘｊ）２＋ｂｊ（ｘ-ｘｊ）＋ｃｊ
ｊは１〜３｛センサ位置（１．３は遠方）｝
この仮想データはあくまでも「発散」防止の手段であるため、実測データより少なく個数を設定する。これにより、個数の上から実測データに重みがかかる解が計算される。重み付けについては、ここでは個数による重み付けを使用しているが、他の重み付け方法を用いてもよい。仮想データの設定が終了すると次にステップＳＴ３に移行する。
Sj (x) = aj (x-xj) 2 + bj (x-xj) + cj
j is 1 to 3 {sensor position (1.3 is far away)}
Since this virtual data is only a means for preventing “divergence”, the number is set smaller than the actual measurement data. Thereby, a solution in which the actually measured data is weighted is calculated from the number. For weighting, weighting by number is used here, but other weighting methods may be used. When the virtual data setting is completed, the process proceeds to step ST3.
ステップＳＴ３においては、最大磁気量の計算を行い、規定値の設定を行う。この最大磁気量は、被測定艦艇の磁性体量から想定する。この規定値の設定の具体例を説明する。図５に示すように、艦艇１を６つのブロックに分け、その各部分に、それぞれ磁気モーメントＭ１、Ｍ２、・・・Ｍ６があると仮定する。今、キール下直近の磁気センサＭＤ８で検出したキール下磁界Ｚ成分を例えば２μ・Ｔとすると艦艇１から、磁気センサＭＤ８までの距離を２０mとして磁気モーメントＭ４の垂直方向磁気量は、約８００μ・Ｔ＾３と推定されるが、余裕をみて、その１．５倍の１２００μ・Ｔ＾３とする。これを最大磁気量の規定値として設定し、メモリ１３の設定値記憶部に記憶する。次にステップＳＴ４に移行する。
In step ST3, the maximum magnetic quantity is calculated, and a prescribed value is set. This maximum magnetic quantity is assumed from the magnetic substance quantity of the ship to be measured. A specific example of setting the specified value will be described. As shown in FIG. 5, it is assumed that the
ステップＳＴ４においては、繰り返し計算回数を設定し、メモリ１３の設定値記憶部に記憶する。続いてステップＳＴ５へ移行する。 In step ST4, the number of repeated calculations is set and stored in the set value storage unit of the memory 13. Subsequently, the process proceeds to step ST5.
ステップＳＴ５においては最小自乗法による正規方程式の演算及び磁気源の計算を行う。この計算は、実測データ及び仮想データを併せて行う。 In step ST5, a normal equation is calculated by a least square method and a magnetic source is calculated. This calculation is performed by combining the actual measurement data and the virtual data.
ここで上記例のように、ダイポールの個数を６個とすると、未知数は、６個×３（Ｍx、Ｍy、Ｍz）で１８個であり、仮に、実測点以外の仮想位置も含めた、位置ｐの点数が１００個であるとすれば、未知数１８個を含む１次結合式が１００個できる。これを行列式で表すと次のとおりとなる。 Here, if the number of dipoles is 6 as in the above example, the number of unknowns is 18 by 6 × 3 (Mx, My, Mz), and the position including the virtual position other than the actual measurement point is assumed. If the number of points of p is 100, 100 linear coupling equations including 18 unknowns can be made. This is expressed as a determinant as follows.
これらから、Ｍｍ（ｍ＝1〜１８）の１８元連立方程式が次のように決まる。 From these, the 18-way simultaneous equations of Mm (m = 1 to 18) are determined as follows.
続いてステップＳＴ６へ移行する。ステップＳＴ６においては全ての磁気源が規定値以下か否かを判定する。規定値を超える磁気源が１個でもあると、ステップＳＴ７へ移行する。一方、計算で得た全ての磁気源が規定値以下の場合は処理を終了する。 Subsequently, the process proceeds to step ST6. In step ST6, it is determined whether or not all the magnetic sources are below a specified value. If there is even one magnetic source exceeding the specified value, the process proceeds to step ST7. On the other hand, if all the magnetic sources obtained by the calculation are below the specified value, the process is terminated.
ステップＳＴ７においては、全ての磁気源の計算回数（繰り返し回数）が設定値以上になったか否かを判定する。計算の繰り返し回数が設定値よりも小さい場合は、ステップＳＴ８へ移行する。一方、繰り返し回数が設定値以上の場合は、ここで処理を終了する。 In step ST7, it is determined whether or not the number of calculations (the number of repetitions) of all magnetic sources has reached a set value or more. When the number of repetitions of calculation is smaller than the set value, the process proceeds to step ST8. On the other hand , if the number of repetitions is equal to or greater than the set value, the process ends here.
ステップＳＴ８においては、それまで設定していた最大磁気量の規定値よりも大きな磁気源を新たな規定値として設定する。例えばそれまで最大磁気量を８００μ・Ｔ＾３の１．５倍の１２００μ・Ｔ＾３としていた場合に更に８００μ・Ｔ＾３の１．６倍の１２８０μ・Ｔ＾３とする。そして、ステップＳＴ５へ戻り、再び磁気源の計算を行い、その後、全ての磁気源が規定値以下となるか、繰り返し回数が規定値以上となるまで、ステップＳＴ６〜ステップＳＴ８の処理を繰り返す。 In step ST8, a magnetic source that is larger than the specified value of the maximum magnetic amount that has been set up to that point is set as a new specified value. For example, when the maximum magnetic amount has been set to 1200 μ · T 3 , 1.5 times 800 μ · T 3 , it is further set to 1280 μ · T 3 , 1.6 times 800 μ · T 3 . Then, returning to step ST5, the calculation of the magnetic source is performed again, and thereafter, the processing of step ST6 to step ST8 is repeated until all the magnetic sources become the specified value or less or the number of repetitions becomes the specified value or more.
１、 艦艇
ＭＤ１、・・、ＭＤＮ 磁気センサ
１１、 磁気計測装置
１２、 ＣＰＵ
１３、 メモリ
１４、 入力ポート
1. Ship MD1, .., MDN Magnetic sensor 11,
13,
Claims (2)
前記実測位置に対して実測データのない仮想位置を設定し、この仮想位置における仮想データを、仮想位置の設定に係る複数の実測位置における複数の実測データに基づく補間により求め、前記実測位置の計算上のデータが当該実測位置の実測データと同じになるように前記実測データ及び前記仮想データを正規方程式に代入し、最小自乗法による前記磁気源の計算を行うことを特徴とする磁気モデル計算方法。 A plurality of external magnetic fields are measured for the magnetic structure, a virtual magnetic source in the magnetic structure is calculated based on the plurality of measured data, and a magnetic field at an arbitrary position other than the actually measured position where the measured data is measured is calculated. In the magnetic model calculation method in the magnetic measurement system to calculate,
A virtual position without actual measurement data is set for the actual measurement position, and virtual data at the virtual position is obtained by interpolation based on a plurality of actual measurement data at a plurality of actual measurement positions related to the setting of the virtual position, and calculation of the actual measurement position is performed. A magnetic model calculation method characterized by substituting the actual measurement data and the virtual data into a normal equation so that the above data becomes the same as the actual measurement data of the actual measurement position, and calculating the magnetic source by the method of least squares .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005130776A JP4815859B2 (en) | 2005-04-28 | 2005-04-28 | Magnetic model calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005130776A JP4815859B2 (en) | 2005-04-28 | 2005-04-28 | Magnetic model calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2006308404A JP2006308404A (en) | 2006-11-09 |
| JP4815859B2 true JP4815859B2 (en) | 2011-11-16 |
Family
ID=37475466
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005130776A Expired - Fee Related JP4815859B2 (en) | 2005-04-28 | 2005-04-28 | Magnetic model calculation method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4815859B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105547254A (en) * | 2016-01-30 | 2016-05-04 | 武汉大学 | Normal equation-based quick-built super-large scale block adjustment method and system |
| CN105716580A (en) * | 2016-01-30 | 2016-06-29 | 武汉大学 | Super-large-scale block adjustment method and system based on adjustment parameter fast resolution |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5162761B2 (en) * | 2009-11-28 | 2013-03-13 | 防衛省技術研究本部長 | Method and apparatus for calculating demagnetizing current of ship hull |
| GB201105450D0 (en) * | 2011-03-31 | 2011-05-18 | Babcock Nuclear Ltd | Improvements in and relating to methods and systems for investigating radioactive sources in locations |
| CN109765292B (en) * | 2019-02-18 | 2024-03-26 | 西南石油大学 | Accurate positioning method for pipeline defects |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60250267A (en) * | 1984-05-25 | 1985-12-10 | Shimadzu Corp | Apparatus for magnetic measurement of ship |
| JP2590443B2 (en) * | 1994-09-07 | 1997-03-12 | 防衛庁技術研究本部長 | Ship magnetic reduction device |
| JP4055416B2 (en) * | 2001-12-28 | 2008-03-05 | 株式会社島津製作所 | Magnetic model calculation method |
-
2005
- 2005-04-28 JP JP2005130776A patent/JP4815859B2/en not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105547254A (en) * | 2016-01-30 | 2016-05-04 | 武汉大学 | Normal equation-based quick-built super-large scale block adjustment method and system |
| CN105716580A (en) * | 2016-01-30 | 2016-06-29 | 武汉大学 | Super-large-scale block adjustment method and system based on adjustment parameter fast resolution |
| CN105547254B (en) * | 2016-01-30 | 2018-01-12 | 武汉大学 | Ultra-large block adjustment method and system based on normal equation rapid build |
| CN105716580B (en) * | 2016-01-30 | 2018-01-23 | 武汉大学 | The ultra-large block adjustment method and system quickly resolved based on adjustment parameter |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2006308404A (en) | 2006-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4815859B2 (en) | Magnetic model calculation method | |
| AU2019211746B2 (en) | Devices, methods, and systems for underwater surveying | |
| CN103649783A (en) | Declination compensation for seismic surveys | |
| CN107044854B (en) | Method for positioning remote magnetic target based on scalar magnetometer array | |
| CN102927984B (en) | Method of eliminating influence of carrier magnetizing field on geomagnetic measurement | |
| AU2012284778B2 (en) | Device for measuring location of underwater vehicle and method thereof | |
| JP2013515253A (en) | Method and apparatus for identifying a subset of measurements, method and system for locating an object, and recording medium for these methods | |
| JP3274308B2 (en) | Magnetic exploration device and its magnetic sensor device | |
| JP2007245791A (en) | Demagnetizing coil settling method, demagnetization control method, demagnetization control device and marine vessel and demagnetizing coil settling program | |
| US8362762B2 (en) | Magnetic field measurement device | |
| JP4259327B2 (en) | Magnetic measurement method for ships | |
| JP2014052188A (en) | Method of estimating inertia characteristics | |
| JP5162760B2 (en) | Hull magnetism reduction method and apparatus | |
| JP2590443B2 (en) | Ship magnetic reduction device | |
| US8584586B1 (en) | Roll frequency dependency correction to control magnetic ship signatures | |
| EP0242391B1 (en) | A magnetic self-ranging system for use in the degaussing of ships | |
| JP2007106140A (en) | Demagnetization processing system for vessel | |
| US6512356B1 (en) | Apparatus and method for determining electrical polarization of distribution of a vessel | |
| JP5162761B2 (en) | Method and apparatus for calculating demagnetizing current of ship hull | |
| JP2020527242A (en) | How to calibrate the magnetic detector | |
| JP4635732B2 (en) | Hull magnetic measurement system | |
| JP4737061B2 (en) | Lateral magnetic field adjustment system | |
| JP2012202781A (en) | Moving body position estimation and detection device and program for moving body position estimation and detection | |
| JP2006335255A (en) | Current adjustment method of demagnetization coil for vessel | |
| CN111337870B (en) | Three-axis electric field correction method based on Euler formula |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070706 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20091007 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101012 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101208 |
|
| 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: 20110802 |
|
| 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: 20110815 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140909 Year of fee payment: 3 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 4815859 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: 20140909 Year of fee payment: 3 |
|
| LAPS | Cancellation because of no payment of annual fees |