KR101698682B1 - 지자기 센서의 출력값을 보정하는 방법 및 장치 - Google Patents

지자기 센서의 출력값을 보정하는 방법 및 장치 Download PDF

Info

Publication number
KR101698682B1
KR101698682B1 KR1020150120477A KR20150120477A KR101698682B1 KR 101698682 B1 KR101698682 B1 KR 101698682B1 KR 1020150120477 A KR1020150120477 A KR 1020150120477A KR 20150120477 A KR20150120477 A KR 20150120477A KR 101698682 B1 KR101698682 B1 KR 101698682B1
Authority
KR
South Korea
Prior art keywords
estimated
gain vector
output value
kth
error covariance
Prior art date
Application number
KR1020150120477A
Other languages
English (en)
Inventor
김동윤
김병기
박승환
김영준
나조안
최연수
신형식
기경석
Original Assignee
매그나칩 반도체 유한회사
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 매그나칩 반도체 유한회사 filed Critical 매그나칩 반도체 유한회사
Priority to KR1020150120477A priority Critical patent/KR101698682B1/ko
Priority to US15/219,993 priority patent/US10060743B2/en
Application granted granted Critical
Publication of KR101698682B1 publication Critical patent/KR101698682B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C17/00—Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
    • G01C17/38—Testing, calibrating, or compensating of compasses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C17/00—Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
    • G01C17/02—Magnetic compasses
    • G01C17/28—Electromagnetic compasses
    • G01C17/30—Earth-inductor compasses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Measuring Magnetic Variables (AREA)
  • Mathematical Physics (AREA)
  • Electromagnetism (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Manufacturing & Machinery (AREA)

Abstract

본 발명은 지자기 센서의 출력값을 보정하는 방법 및 그 장치에 관한 것으로, 매번 지자기 센서의 출력값 측정시 추정 오차(오프셋 및 반지름 값)를 추적함으로써 지자기 센서의 출력값에 대한 보정의 정확도를 향상시킬 수 있는 방법 및 그 장치에 관한 것이다.

Description

지자기 센서의 출력값을 보정하는 방법 및 장치{Method and Apparatus of correcting output value of terrestrial magnetism sensor}
본 발명은 지자기 센서의 출력값을 보정하는 방법 및 그 장치에 관한 것으로, 매번 지자기 센서의 출력값 측정시 추정 오차(오프셋 및 반지름 값)를 추적함으로써 지자기 센서의 출력값에 대한 보정의 정확도를 향상시킬 수 있는 방법 및 그 장치에 관한 것이다.
지자기센서를 이용하여 측정되는 데이터에는 지구 자기장 외에 바이어스가 포함되어 있다. 이 바이어스는 센서 주변에서 고정적으로 자기를 발산 혹은 상쇄하는 물리적인 성질의 강철효과와, 센서의 측정 오차, 각 축의 직교 오차, 연철 효과 등 다양한 오류 가능성을 내포하고 있다. 지자기 센서의 데이터 보정은 지자기센서를 통해 측정된 데이터로부터 예상되는 바이어스를 제거함으로써 실제 지구 자기장의 성분만을 분리하는 것을 의미하고, 보정 성능은 정확한 바이어스를 추정해 내는 능력이다.
2축 혹은 3축 지자기센서로부터 측정된 데이터는 원, 혹은 구의 표면에 분포하는 성질이 있으며 바이어스에 따라 원 혹은 구의 중점과 반지름의 추정치가 달라진다. 바이어스를 데이터 관점에서 오프셋이라고 가정할 때, 원 혹은 구의 중점은 원점으로부터 오프셋만큼 떨어진 곳에 위치하게 된다. 측정된 지자기 데이터로부터 오프셋을 추정하는 방법은 일반적으로 최소제곱법을 이용한다.
최소제곱법을 이용하여 오프셋과 반지름을 추정할 때, 2축의 경우에는 3개 이상, 3축의 경우에는 4개 이상의 측정데이터로 추정이 가능하다. 이상적으로 측정오차 없이 고정적인 강철효과만 존재한다면 측정 데이터 4개로도 충분히 정확도 높은 오프셋을 추정할 수 있다. 하지만 측정데이터는 다양한 측정노이즈를 동반하므로 추정오류 발생 가능성을 배제할 수 없다. 최소제곱법을 이용하여 지자기 센서 데이터를 보정하고자 할 때, 측정 데이터의 양이 많으면 측정노이즈의 비중이 적어 종래기술에서는 측정노이즈를 고려하지 않는다. 하지만 실제 지자기센서를 이용하는 시스템에서는 측정 데이터의 양을 늘리면 연산 부하에 따른 처리 지연을 가져올 수 있고, 저장 공간의 부담이 있을 수 있어 측정 데이터를 일정 양 만큼 나누어 처리해야 한다. 일정 양 만큼의 측정 데이터 집합으로부터 오프셋을 추정할 때, 집합 간의 추정한 값의 산포가 발생하며, 결국 추정오차가 클 경우 정확도가 낮아지는 문제점이 있었다.
US 7,177,779 EP 00451839 US 5,151,872
본 발명은 상기와 같은 문제점을 해결하기 위해 제안된 발명으로써, 종래 기술과 달리 현재 시점의 상태 정보(오프셋, 반지름 등)를 재귀적으로(reculsive) 추정함으로써 높은 보정 정확도를 갖는 지자기 센서의 출력값의 보정 방법 및 장치를 제공하고자 한다.
본 발명의 일 예에 따라 순차적으로 획득되는 지자기 센서의 출력값에 기반하여 추정 상태, 추정게인벡터 및 추정오차공분산을 재귀적으로 산출하고, 상기 추정 상태에 대한 값에 기반하여 지자기 센서의 출력값을 보정하는 방법은, k번째 지자기 센서의 출력값을 획득하는 단계; k-1번째 추정오차공분산 및 k번째 출력값에 기반하여 k번째 추정게인벡터를 산출하는 단계; 산출된 k번째 추정게인벡터, k번째 출력값 및 k-1번째 추정 상태에 기반하여 k번째 추정 상태를 산출하는 단계; 상기 k번째 추정게인벡터 및 k-1번째 추정오차공분산을 이용하여 k번째 추정오차공분산을 산출하는 단계; 및 상기 k번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정하는 단계;를 포함한다. (단, k는 1보다 큰 자연수)
상기 지자기 센서의 출력값은, 2축 지자기 센서의 출력값 또는 3축 지자기 센서의 출력값이 적용될 수 있다.
상기 k번째 추정게인벡터를 산출하는 단계는, 아래의 수학식에 따라 k번째 추정게인벡터를 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00001
(여기서, Kk는 k번째 추정게인벡터, Pk -1은 k-1번째 추정오차공분산, Hk는 k번째 측정 행렬(measurement matrix), Rk는 k번째 측정 노이즈의 공분산)
상기 k번째 추정 상태를 산출하는 단계는, 아래의 수학식에 따라 k번째 추정 상태를 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00002
(여기서, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix), yk는 k번째 지자기 센서 측정값)
상기 k번째 추정오차공분산을 산출하는 단계는, 아래의 수학식에 따라 k번째 추정오차공분산을 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00003
(여기서, Pk -1은 k-1번째 추정오차공분산, I는 단위 행렬, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix))
본 발명의 다른 측면에 따라 순차적으로 획득되는 지자기 센서의 출력값에 기반하여 추정 상태, 추정게인벡터 및 추정오차공분산을 재귀적으로 산출하고, 상기 추정 상태에 대한 값에 기반하여 지자기 센서의 출력값을 보정하는 장치는, k-1번째 추정 상태, k-1번째 추정게인벡터, k-1번째 추정오차공분산 정보를 저장하는 저장부; k번째 지자기 센서의 출력값을 획득하는 센서부; 상기 저장부에 저장된 k-1번째 추정오차공분산 및 상기 센서부를 통해 획득된 k번째 출력값에 기반하여 k번째 추정게인벡터를 산출하여 저장부에 저장하는 추정게인벡터 산출부; 상기 k번째 추정게인벡터, k번째 출력값 및 k-1번째 추정 상태에 기반하여 k번째 추정 상태를 산출하여 저장부에 저장하는 추정 상태 산출부; 상기 k번째 추정게인벡터 및 k-1번째 추정오차공분산을 이용하여 k번째 추정오차공분산을 산출하는 추정오차공분산 산출부; 및 상기 k번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정하는 출력값 보정부;를 포함한다. (단, k는 1보다 큰 자연수)
상기 지자기 센서의 출력값은, 2축 지자기 센서의 출력값 또는 3축 지자기 센서의 출력값이 적용될 수 있다.
상기 추정게인벡터 산출부는, 아래의 수학식에 따라 k번째 추정게인벡터를 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00004
(여기서, Kk는 k번째 추정게인벡터, Pk -1은 k-1번째 추정오차공분산, Hk는 k번째 측정 행렬(measurement matrix), Rk는 k번째 측정 노이즈의 공분산)
상기 추정 상태 산출부는, 아래의 수학식에 따라 k번째 추정 상태를 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00005
(여기서, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix), yk는 k번째 지자기 센서 측정값)
상기 추정오차공분산 산출부는, 아래의 수학식에 따라 k번째 추정오차공분산을 산출할 수 있다.
[수학식]
Figure 112015083023916-pat00006
(여기서, Pk -1은 k-1번째 추정오차공분산, I는 단위 행렬, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix))
상기 저장부는, 백업 추정 상태, 백업 추정게인벡터 및 백업 추정오차공분산을 추가적으로 더 저장하고,
상기 센서부가 1번째 지자기 출력 값을 획득하면,
상기 추정게인벡터 산출부는 상기 백업 추정오차공분산 및 1번째 출력값에 기반하여 1번째 추정게인벡터를 산출하여 저장부에 저장하고,
상기 추정 상태 산출부는 1번째 추정게인벡터, 1번째 출력값 및 백업 추정 상태에 기반하여 1번째 추정 상태를 산출하여 저장부에 저장하고,
상기 추정오차공분산 산출부는 1번째 추정게인벡터 및 백업 추정오차공분산을 이용하여 1번째 추정오차공분산을 산출하고,
상기 출력값 보정부는 상기 1번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정할 수 있다.
이때, 상기 백업 추정 상태, 백업 추정게인벡터 및 백업 추정오차공분산은, 각각 상기 지자기 센서의 동작이 멈추는 시점의 추정 상태, 추정게인벡터, 추정오차공분산 값이 적용될 수 있다.
본 발명에 따른 지자기 센서의 출력값을 보정하는 방법 및 장치는 현재 시점의 사태 정보를 재귀적으로 추정하기 위해 칼만필터의 동적 시스템 모델을 적용함으로써 높은 보정 정확도로 지자기 센서의 출력값의 보정할 수 있다는 효과가 있다.
또한, 본 발명과 종래기술간 시뮬레이션 결과에서 확인할 수 있듯이, 본 발명은 종래 기술보다 높은 정확도로 레퍼렌스의 오프셋과 반지름 값을 추정할 수 있다는 효과가 있다.
도 1은 본 발명의 일 예에 따른 지자기 센서의 출력값을 보정하는 방법을 나타낸 도면,
도 2는 본 발명의 다른 예에 따른 지자기 센서의 출력값을 보정하는 장치를 나타낸 도면,
도 3은 종래 기술과 본 발명에 따라 출력되는 지자기 센서의 출력값을 비교한 시뮬레이션 결과이다.
본 발명은 다양한 변환을 가할 수 있고 여러 가지 실시 예를 가질 수 있는 바, 특정 실시 예들을 도면에 예시하고 상세한 설명에 상세하게 설명하고자 한다. 그러나, 이는 본 발명의 특정한 실시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변환, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다. 본 발명을 설명함에 있어서 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 흐릴 수 있다고 판단되는 경우 그 상세한 설명을 생략한다.
제1, 제2 등의 용어는 다양한 구성요소들을 설명하는데 사용될 수 있지만, 상기 구성요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다.
본 출원에서 사용한 용어는 단지 특정한 실시 예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함한다. 본 출원에서, "포함하다" 또는 "가지다" 등의 용어는 명세서상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다.
공간적으로 상대적인 용어인 아래(below, beneath, lower), 위(above, upper) 등은 도면에 도시되어 있는 바와 같이 하나의 소자 또는 구성 요소들과 다른 소자 또는 구성 요소들과의 상관 관계를 용이하게 기술하기 위해 사용될 수 있다. 공간적으로 상대적인 용어는 도면에 도시되어 있는 방향에 더하여 사용시 또는 동작시 소자의 서로 다른 방향을 포함하는 용어로 이해되어야 한다. 예를 들면, 도면에 도시되어 있는 소자를 뒤집을 경우, 다른 소자의 아래(below, beneath)로 기술된 소자는 다른 소자의 위(above, upper)에 놓여질 수 있다. 따라서, 예시적인 용어인 아래는 아래와 위의 방향을 모두 포함할 수 있다. 소자는 다른 방향으로도 배향될 수 있고, 이에 따라 공간적으로 상대적인 용어들은 배향에 따라 해석될 수 있다.
이하, 본 발명의 실시예를 첨부한 도면들을 참조하여 상세히 설명하기로 한다.
도 1은 본 발명의 일 예에 따른 지자기 센서의 출력값을 보정하는 방법을 나타낸 도면이다.
본 발명에 따른 지자기 센서의 출력값을 보정하는 방법은 현재까지 측정된 데이터 중 일부 개수의 데이터를 선별하여 실제 지자기 값을 추정하는 종래 기술과 컨셉을 달리한다. 본 발명에서는 이전 시점에서의 실제 지자기 값에서 현재 시점의 추정오류를 예측하여 현재 시점에서의 실제 지자기 값을 재귀적(recursive)으로 추정하는 것을 특징으로 한다. 이를 위해, 본 발명에서는 칼만필터의 동적 시스템 모델을 적용한다.
이를 위해 본 발명에서는 순차적으로 획득되는 지자기 센서의 출력값에 기반하여 추정 상태, 추정게인벡터 및 추정오차공분산을 재귀적으로 산출하고, 상기 값에 기반하여 지자기 센서의 출력값을 보정한다. 이에 대한 상세한 설명 및 청구항에서는 현재 시점과 바로 직전인 이전 시점의 추정 상태, 추정게인벡터 및 추정오차공분산을 각각 k번째, k-1번째로 표시함으로써 현재 시점과 이전 시점을 구분한다.
먼저, 본 발명에 따른 지자기 센서의 출력값을 보정하기 위해 현재 시점의 지자기 센서의 출력값을 획득한다(S10).
이때, 상기 지자기 센서의 출력값은 지자기 센서의 축의 개수에 따라 결정될 수 있다. 일 예로, 2축 지자기 센서의 경우에는 (hmx, hmy), 3축 지자기 센서의 경우에는 (hmx, hmy, hmz) 과 같은 지자기 센서의 출력값을 획득할 수 있다. 본 발명에 있어 상기 S10 단계를 통해 획득되는 지자기 센서의 출력값의 형태는 중요하지 않으며, 이는 다양한 형태의 지자기 센서의 출력값에도 적용 가능하다.
이하, 본 발명에 대해서는 3축 지자기 센서의 경우로 구체화하여 설명하나, 본 발명은 이에 한정되지 않는다.
이어, 현재 시점의 추정게인벡터를 산출한다(S20). 보다 구체적으로, 상기 S20 단계에서는 이전 시점의 추정오차공분산, 현재 시점의 지자기 센서의 출력값에 기반하여 현재 시점의 추정게인벡터를 산출하며, 하기의 수학식 1에 따라 산출될 수 있다.
Figure 112015083023916-pat00007
(여기서, Kk는 k번째 추정게인벡터, Pk -1은 k-1번째 추정오차공분산, Hk는 k번째 측정 행렬(measurement matrix), Rk는 k번째 측정 노이즈의 공분산)
상기 수학식 1에 있어, 현재 시점의 측정 행렬(measurement matrix)은 하기의 수학식 2와 같이 획득할 수 있다.
Figure 112015083023916-pat00008
상기 S20 단계를 통해 산출된 현재 시점의 추정게인벡터를 활용하여 현재 시점의 추정 상태를 산출한다(S30). 상기 S30 단계에 있어 추정 상태란 특정 시점의 추정 오프셋(offset) 값으로써 상기 S10 단계를 통해 획득된 지자기 센서의 출력값을 보정하기 위한 값으로 활용된다.
보다 구체적으로, 상기 S30 단계에서는 하기의 수학식 3에 따라 현재 시점의 추정 상태를 산출한다.
Figure 112015083023916-pat00009
(여기서, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix), yk는 k번째 지자기 센서 측정값)
상기 수학식 3을 참고하면, 현재 시점의 추정 상태(오프셋 값)를 산출하기 위하여 이전 시점의 추정상태, 현재 시점의 추정게인벡터, 현재시점의 지자기 센서의 측정값을 이용한다.
이와 같이, 본 발명에서는 현재 시점의 추정 상태를 재귀적으로 추정하는 칼만필터의 동적 시스템 모델을 적용하였다. 이는 실시간으로 오프셋 값을 추정할 때, 현재 시점에서의 측정 노이즈 값의 비중이 커져 이를 간과할 수 없기 때문이다.
상기 수학식 3에서 yk-Hkxk -1 은 Correction term이 된다. 즉, 현재 시점의 추정 상태 값은 이전 시점의 추정 상태 값, 현재 시점의 추정게인벡터를 통해 측정 노이즈에 대한 보정을 예측하여 적용했기 때문에 측정노이즈는 이미 고려되어 졌다고 볼 수 있다.
이와 같은 과정을 통해 획득된 현재 시점의 추정 상태로부터 현재 시점의 추정 오프셋값 및 추정 반지름 값을 획득할 수 있다. 일 실시예로, 하기와 같은 수학식4에 따라 각 추정 오프셋값 및 추정 반지름 값을 산출할 수 있다.
Figure 112015083023916-pat00010
(여기서, Hrx, Hry, Hrz는 각 축별 추정 오프셋, xk는 k번째 추정 상태, r은 현재 시점의 추정 반지름)
또한, 현재 시점의 추정게인벡터, 이전 시점의 추정오차공분산을 이용하여 현재 시점의 추정오차공분산을 산출한다(S40).
상기 S40 단계에 있어, 현재 시점의 추정오차공분산은 하기의 수학식 5에 의해 산출될 수 있다.
Figure 112015083023916-pat00011
(여기서, Pk -1은 k-1번째 추정오차공분산, I는 단위 행렬, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix))
상기 S40 단계를 통해 획득된 현재 시점의 추정오차공분산은 이후 시점의 추정게인벡터 및 추정오차공분산을 획득하는데 활용 가능하다.
최종적으로, 현재 시점의 추정 상태에 기반하여 지자기 센서의 출력값을 보정한다(S50).
일 예로, 3축의 지자기 센서의 출력값인 경우에는 하기의 수학식 6과 같이 지자기 데이터 값을 보정할 수 있다.
Figure 112015083023916-pat00012
(여기서, hmx, hmy, hmz는 각 축별 지자기 센서의 출력값, Hrx, Hry, Hrz는 각 축별 추정 오프셋)
이하, 본 발명에 따른 지자기 센서의 출력값을 보정하는 방법에 적용되는 각각의 수학식이 도출된 과정에 대하여 추가적으로 설명한다.
상기 수학식 3을 참고할 때, 현재 시점에서의 추정오차(εk)는 하기의 수학식 7과 같이 표현할 수 있다.
Figure 112015083023916-pat00013
(여기서, εk는 k번째 추정오차, x는 알 수 없는 벡터 상수, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, vk는 k번째 추정 노이즈 크기, I는 단위행렬)
이때, E(εk)는 하기의 수학식 8과 같이 정의할 수 있다.
Figure 112015083023916-pat00014
또한, 현재 시점에서 최소제곱법을 이용하여 추정오차를 최소제곱법을 이용하여 추정오차를 최소화 하기 위한 Jk를 하기의 수학식 9와 같이 정의할 수 있다.
Figure 112015083023916-pat00015
(여기서, tr(Pk)는 k번째 추적 연산자, Pk는 k번째 추정오차공분산(estimation-error covariation))
이전 시점에서의 추정오차 εk - 1는 현재 시점에서의 측정노이즈 vk와는 독립적이다. 이에, Pk는 하기의 수학식 10과 같이 다시 정리할 수 있다.
Figure 112015083023916-pat00016
(여기서, Rk는 E(vkvk T)로 vk의 공분산)
추정게인벡터를 구하기 위해 Jk를 Kk로 편미분하게 된다. 이에 따라 추정게인벡터는 하기의 수학식 11과 같이 구할 수 있다.
Figure 112015083023916-pat00017
여기서,
Figure 112015083023916-pat00018
라고 할 때, Kk와 Pk를 하기의 수학식 12 및 13과 같이 정리할 수 있다.
Figure 112015083023916-pat00019
Figure 112015083023916-pat00020
상기와 같이 본 발명에서는 측정오차가 발생하더라도 실시간적으로 추정게인벡터 Kk와 추정오차공분산 Pk를 산출함으로써 현재시점에서의 x에 근접한 xk를 추정할 수 있다.
도 2는 본 발명의 다른 예에 따른 지자기 센서의 출력값을 보정하는 장치를 나타낸 도면이다.
도 2에 도시된 바와 같이, 본 발명에 따른 지자기 센서의 출력값 보정 장치(1)는 저장부(10), 센서부(20), 추정게인벡터 산출부(30), 추정 상태 산출부(40), 추정오차공분산 산출부(50), 출력값 보정부(60)를 포함한다.
이하, 도 1 및 이에 대한 상세한 설명을 통해 상세히 설명한 내용에 대해서는 간략히 설명한다.
저장부(10)는 이전 시점의 추정 상태, 추정게인벡터, 추정오차공분산 정보를 저장한다. 이때, 상기 저장부(10)는 가장 최근의 추정 상태, 추정게인벡터, 추정오차공분산 정보를 저장할 수도 있으며, 현 시점까지의 추정 상태, 추정게인벡터, 추정오차공분산 정보를 시점별로 구분하여 모두 저장할 수도 있다. 이는 사용자(관리자)의 설정에 따라 달리 설정될 수 있다.
또한, 상기 저장부(10)는 현재까지 획득된 지자기 센서의 출력값을 저장할 수 있다. 이를 통해, 추정게인벡터 산출부(30), 추정 상태 산출부(40)는 현재 시점의 지자기 센서의 출력값을 센서부(20)가 아닌 저장부(10)로부터 제공받을 수 있다.
바람직한 실시예에서, 상기 저장부(10)는 백업 추정 상태, 백업 추정게인벡터 및 백업 추정오차공분산을 추가적으로 더 저장할 수 있다.
이때, 백업 추정 상태, 백업 추정게인벡터, 백업 추정오차공분산 정보는 지자기 센서의 사용을 멈추는 시점의 대응되는 정보가 적용될 수 있다.
즉, 저장부(10)는 지자기 센서의 사용(측정)이 멈추는 시점을 기준으로 가장 최근의 추정 상태, 추정게인벡터, 추정오차공분산 정보를 각각 백업 추정 상태, 백업 추정게인벡터, 백업 추정오차공분산으로 저장할 수 있다.
차후, 지자기 센서의 사용(측정 동작)이 수행되게 되면, 추정게인벡터 산출부(30), 추정 상태 산출부(40), 추정오차공분산 산출부(50)는 상기 저장부(10)에 미리 저장되어 있는 백업 추정 상태, 백업 추정게인벡터, 백업 추정오차공분산 정보를 이용하여 1번째 추정게인벡터, 추정 상태, 추정오차공분산을 산출할 수 있다.
이와 같은 방법을 통해, 본 발명에 따른 지자기 센서의 출력값 보정 장치(1)는 도 3b에 도시된 초기에 안정화되지 않는 문제점을 극복할 수 있다.
센서부(20)는 현재 시점의 지자기 센서의 출력값을 획득한다. 보다 구체적으로는, 도 1의 S10 단계와 같이 상기 센서부(20)는 지자기 센서의 출력값을 획득한다. 이하 상세한 설명은 생략한다.
추정게인벡터 산출부(30)는 이전 시점의 추정오차공분산, 현재 시점의 지자기 센서의 출력값을 이용하여 현재 시점의 추정게인벡터를 산출한다. 이어, 산출된 현재 시점의 추정게인벡터를 상기 저장부(10)에 저장한다.
본 발명에 있어, 각 시점의 추정오차공분산, 추정 상태, 추정게인벡터 정보는 모두 저장부(10)에 저장될 수 있으며, 추정게인벡터 산출부(30) 등은 상기 저장부(10)에 저장된 정보를 이용하여 현재 시점의 추정게인벡터 등을 산출할 수 있다.
추정 상태 산출부(40), 추정오차공분산 산출부(50)는, 상기 추정게인벡터 산출부(30)와 유사하게, 상기 저장부(10)에 저장된 정보들을 이용하여 각각 현재 시점의 추정 상태, 추정오차공분산을 산출할 수 있다.
이하, 추정게인벡터 산출부(30), 추정 상태 산출부(40), 추정오차공분산 산출부(50)의 동작에 대해서는 도 1의 S20, S30, S40 단계에 대한 설명을 통해 상세히 설명한 바 생략한다.
출력값 보정부(60)는 상기 추정 상태 산출부(40)를 통해 산출된 현재 시점의 추정 상태 정보에 기반하여 지자기 센서의 출력값을 보정한다. 보다 구체적인 동작 사항은 도 1의 S50 단계를 통해 상세히 설명한 바 생략한다.
도 3은 종래 기술과 본 발명에 따라 출력되는 지자기 센서의 출력값을 비교한 시뮬레이션 결과이다. X축은 측정 데이터 수를, Y축은 추정된 오프셋 값과 반지름 값 및 각각의 레퍼런스 값을 나타낸다.
도 3a에 도시된 종래 기술의 경우, 측정 데이터를 20개 단위로 구분하여 총 20개의 그룹으로 구분하였으며, 각각의 그룹으로부터 20개의 오프셋과 반지름을 추정하였다. 이때, 전체 측정 데이터 400개를 1개의 그룹으로 구분하였을 때 추정된 오프셋과 반지름을 레퍼런스(Ref)로 하였다.
이때, 도 3a에 도시된 바와 같이, 종래 기술에 따라 산출된 각 그룹의 오프셋(offset.X, offset.Y, offset.Z)과 반지름(radius) 값은 레퍼런스의 오프셋(Ref.offset.X, Ref.offset.Y, Ref.offset.Z)과 레퍼런스의 반지름(Ref.radius) 값을 기준으로 큰 산포를 보이고 있다. 또한, 각 그룹의 측정 데이터에는 측정 노이즈가 그룹별로 다르게 발생하였으며, 그룹간의 편차도 매우 크게 나타나고 있다.
반면, 도 3b에 도시된 본 발명의 경우, 측정 데이터 1개 단위로 오프셋(offset.X, offset.Y, offset.Z)과 반지름(radius)을 추정하였으며, 첫번째 측정 데이터에서는 이전 데이터가 없었기 때문에 추정이 생략되어 총 399번의 추정 결과를 얻었다.
이때, 도 3b에 도시된 바와 같이, 추정 초기에는 이전 시점의 추정 상태가 안정화되지 않았으므로 큰 편차를 보이지만 시간이 갈수록 레퍼런스의 오프셋(Ref.offset.X, Ref.offset.Y, Ref.offset.Z)과 레퍼런스의 반지름(Ref.radius) 값에 근접한 값을 추정하고 있음을 확인할 수 있다.
한편, 다른 실시예에서, 이전 시점의 추정 상태, 추정오차공분산, 추정게인벡터 값은 지자기 센서의 사용을 멈추는 시점을 기준으로 백업을 해두었다가 지자기 센서를 다시 사용할 때 백업된 데이터를 로딩함으로써 초기에 안정화되지 않는 문제점을 극복할 수 있다.
이제까지 본 발명에 대하여 그 바람직한 실시예들을 중심으로 살펴보았다. 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자는 본 발명이 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현될 수 있음을 이해할 수 있을 것이다. 그러므로 개시된 실시예들은 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 한다. 본 발명의 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다.
1: 지자기 센서의 출력값 보정 장치
10: 저장부 20: 센서부
30: 추정게인벡터 산출부 40: 추정 상태 산출부
50: 추정오차공분산 산출부 60: 출력값 보정부

Claims (12)

  1. 순차적으로 획득되는 지자기 센서의 출력값에 기반하여 추정 상태, 추정게인벡터 및 추정오차공분산을 재귀적으로 산출하고, 상기 추정 상태에 대한 값에 기반하여 지자기 센서의 출력값을 보정하는 방법에 있어서,
    k번째 지자기 센서의 출력값을 획득하는 단계;
    k-1번째 추정오차공분산 및 k번째 출력값에 기반하여 k번째 추정게인벡터를 산출하는 단계;
    산출된 k번째 추정게인벡터, k번째 출력값 및 k-1번째 추정 상태에 기반하여 k번째 추정 상태를 산출하는 단계;
    상기 k번째 추정게인벡터 및 k-1번째 추정오차공분산을 이용하여 k번째 추정오차공분산을 산출하는 단계; 및
    상기 k번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정하는 단계;를 포함하는 지자기 센서의 출력값을 보정하는 방법.
    (단, k는 1보다 큰 자연수)
  2. 제 1항에 있어서,
    상기 지자기 센서의 출력값은,
    2축 지자기 센서의 출력값 또는 3축 지자기 센서의 출력값인 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 방법.
  3. 제 1항에 있어서,
    상기 k번째 추정게인벡터를 산출하는 단계는,
    아래의 수학식에 따라 k번째 추정게인벡터를 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 방법.
    [수학식]
    Figure 112015083023916-pat00021

    (여기서, Kk는 k번째 추정게인벡터, Pk -1은 k-1번째 추정오차공분산, Hk는 k번째 측정 행렬(measurement matrix), Rk는 k번째 측정 노이즈의 공분산)
  4. 제 1항에 있어서,
    상기 k번째 추정 상태를 산출하는 단계는,
    아래의 수학식에 따라 k번째 추정 상태를 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 방법.
    [수학식]
    Figure 112015083023916-pat00022

    (여기서, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix), yk는 k번째 지자기 센서 측정값)
  5. 제 1항에 있어서,
    상기 k번째 추정오차공분산을 산출하는 단계는,
    아래의 수학식에 따라 k번째 추정오차공분산을 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 방법.
    [수학식]
    Figure 112015083023916-pat00023

    (여기서, Pk -1은 k-1번째 추정오차공분산, I는 단위 행렬, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix))
  6. 순차적으로 획득되는 지자기 센서의 출력값에 기반하여 추정 상태, 추정게인벡터 및 추정오차공분산을 재귀적으로 산출하고, 상기 추정 상태에 대한 값에 기반하여 지자기 센서의 출력값을 보정하는 장치에 있어서,
    k-1번째 추정 상태, k-1번째 추정게인벡터, k-1번째 추정오차공분산 정보를 저장하는 저장부;
    k번째 지자기 센서의 출력값을 획득하는 센서부;
    상기 저장부에 저장된 k-1번째 추정오차공분산 및 상기 센서부를 통해 획득된 k번째 출력값에 기반하여 k번째 추정게인벡터를 산출하여 저장부에 저장하는 추정게인벡터 산출부;
    상기 k번째 추정게인벡터, k번째 출력값 및 k-1번째 추정 상태에 기반하여 k번째 추정 상태를 산출하여 저장부에 저장하는 추정 상태 산출부;
    상기 k번째 추정게인벡터 및 k-1번째 추정오차공분산을 이용하여 k번째 추정오차공분산을 산출하는 추정오차공분산 산출부; 및
    상기 k번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정하는 출력값 보정부;를 포함하는 지자기 센서의 출력값을 보정하는 장치.
    (단, k는 1보다 큰 자연수)
  7. 제 6항에 있어서,
    상기 지자기 센서의 출력값은,
    2축 지자기 센서의 출력값 또는 3축 지자기 센서의 출력값인 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
  8. 제 6항에 있어서,
    상기 추정게인벡터 산출부는,
    아래의 수학식에 따라 k번째 추정게인벡터를 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
    [수학식]
    Figure 112015083023916-pat00024

    (여기서, Kk는 k번째 추정게인벡터, Pk -1은 k-1번째 추정오차공분산, Hk는 k번째 측정 행렬(measurement matrix), Rk는 k번째 측정 노이즈의 공분산)
  9. 제 6항에 있어서,
    상기 추정 상태 산출부는,
    아래의 수학식에 따라 k번째 추정 상태를 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
    [수학식]
    Figure 112015083023916-pat00025

    (여기서, xk는 k번째 추정 상태, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix), yk는 k번째 지자기 센서 측정값)
  10. 제 6항에 있어서,
    상기 추정오차공분산 산출부는,
    아래의 수학식에 따라 k번째 추정오차공분산을 산출하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
    [수학식]
    Figure 112015083023916-pat00026

    (여기서, Pk -1은 k-1번째 추정오차공분산, I는 단위 행렬, Kk는 k번째 추정게인벡터, Hk는 k번째 측정 행렬(measurement matrix))
  11. 제 6항에 있어서,
    상기 저장부는,
    백업 추정 상태, 백업 추정게인벡터 및 백업 추정오차공분산을 추가적으로 더 저장하고,
    상기 센서부가 1번째 지자기 출력 값을 획득하면,
    상기 추정게인벡터 산출부는 상기 백업 추정오차공분산 및 1번째 출력값에 기반하여 1번째 추정게인벡터를 산출하여 저장부에 저장하고,
    상기 추정 상태 산출부는 1번째 추정게인벡터, 1번째 출력값 및 백업 추정 상태에 기반하여 1번째 추정 상태를 산출하여 저장부에 저장하고,
    상기 추정오차공분산 산출부는 1번째 추정게인벡터 및 백업 추정오차공분산을 이용하여 1번째 추정오차공분산을 산출하고,
    상기 출력값 보정부는 상기 1번째 추정 상태에 기반하여 지자기 센서의 출력값을 보정하는 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
  12. 제 11항에 있어서,
    상기 백업 추정 상태, 백업 추정게인벡터 및 백업 추정오차공분산은,
    각각 상기 지자기 센서의 동작이 멈추는 시점의 추정 상태, 추정게인벡터, 추정오차공분산 값인 것을 특징으로 하는 지자기 센서의 출력값을 보정하는 장치.
KR1020150120477A 2015-08-26 2015-08-26 지자기 센서의 출력값을 보정하는 방법 및 장치 KR101698682B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020150120477A KR101698682B1 (ko) 2015-08-26 2015-08-26 지자기 센서의 출력값을 보정하는 방법 및 장치
US15/219,993 US10060743B2 (en) 2015-08-26 2016-07-26 Method and apparatus of correcting output value of geomagnetic sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150120477A KR101698682B1 (ko) 2015-08-26 2015-08-26 지자기 센서의 출력값을 보정하는 방법 및 장치

Publications (1)

Publication Number Publication Date
KR101698682B1 true KR101698682B1 (ko) 2017-01-23

Family

ID=57989793

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150120477A KR101698682B1 (ko) 2015-08-26 2015-08-26 지자기 센서의 출력값을 보정하는 방법 및 장치

Country Status (2)

Country Link
US (1) US10060743B2 (ko)
KR (1) KR101698682B1 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101938564B1 (ko) * 2017-02-15 2019-01-15 주식회사 엠앤디 센서 교정 유무 판단 방법
KR101985458B1 (ko) * 2018-11-15 2019-06-03 주식회사 미래와도전 지진 센서의 오차를 예측하는 지진 분석 장치 및 오차 예측 방법
KR102357109B1 (ko) * 2021-02-17 2022-02-08 주식회사 비엘 시공중 터널 막장면 상태평가 시스템

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017111926A1 (de) 2017-05-31 2018-12-06 Infineon Technologies Ag Prozesssteuerschaltung und Verfahren zum Steuern einer Bearbeitungsanordnung
CN108490472B (zh) * 2018-01-29 2021-12-03 哈尔滨工程大学 一种基于模糊自适应滤波的无人艇组合导航方法
CN108897053B (zh) * 2018-04-02 2021-03-30 广州吉欧电子科技有限公司 一种基于设备磁环境变化的磁干扰检测方法及装置
CN115200607A (zh) * 2021-04-09 2022-10-18 广东小天才科技有限公司 一种误差系数计算方法及电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0451839A2 (en) 1990-04-12 1991-10-16 Matsushita Electric Industrial Co., Ltd. Direction finding apparatus
US5151872A (en) 1989-10-17 1992-09-29 Pioneer Electronic Corporation Method and apparatus for correcting the output of an onboard vehicle terrestrial magnetism sensor
US7177779B2 (en) 2002-07-01 2007-02-13 Asahi Kasei Emd Corporation Azimuth measuring device and azimuth measuring method
KR20070049419A (ko) * 2005-11-08 2007-05-11 한국전자통신연구원 2축 지자계 센서의 경사각 오차를 보상하기 위한 가상축지자계 데이터 추정 장치 및 그 방법과, 그를 이용한방위각 산출 장치
KR20140093111A (ko) * 2013-01-17 2014-07-25 삼성전자주식회사 지자기 센서 오차 보정 장치 및 지자기 센서 오차 보정 방법
JP5569681B2 (ja) * 2010-04-23 2014-08-13 国立大学法人 東京大学 慣性センサ，磁気センサおよび速度計を用いた移動体の姿勢推定装置および姿勢推定方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020008661A1 (en) * 2000-07-20 2002-01-24 Mccall Hiram Micro integrated global positioning system/inertial measurement unit system
US6697736B2 (en) * 2002-02-06 2004-02-24 American Gnc Corporation Positioning and navigation method and system thereof
US8615253B2 (en) * 2011-06-03 2013-12-24 Apple Inc. State estimation using motion context and multiple input observation types
US20130060516A1 (en) * 2011-09-06 2013-03-07 Yin-Chen CHANG Trace-generating devices and methods thereof
JP2013064695A (ja) * 2011-09-20 2013-04-11 Yamaha Corp 状態推定装置、オフセット更新方法およびオフセット更新プログラム
JP2013096724A (ja) * 2011-10-28 2013-05-20 Yamaha Corp 状態推定装置
JP2014089113A (ja) * 2012-10-30 2014-05-15 Yamaha Corp 姿勢推定装置及びプログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5151872A (en) 1989-10-17 1992-09-29 Pioneer Electronic Corporation Method and apparatus for correcting the output of an onboard vehicle terrestrial magnetism sensor
EP0451839A2 (en) 1990-04-12 1991-10-16 Matsushita Electric Industrial Co., Ltd. Direction finding apparatus
US7177779B2 (en) 2002-07-01 2007-02-13 Asahi Kasei Emd Corporation Azimuth measuring device and azimuth measuring method
KR20070049419A (ko) * 2005-11-08 2007-05-11 한국전자통신연구원 2축 지자계 센서의 경사각 오차를 보상하기 위한 가상축지자계 데이터 추정 장치 및 그 방법과, 그를 이용한방위각 산출 장치
JP5569681B2 (ja) * 2010-04-23 2014-08-13 国立大学法人 東京大学 慣性センサ，磁気センサおよび速度計を用いた移動体の姿勢推定装置および姿勢推定方法
KR20140093111A (ko) * 2013-01-17 2014-07-25 삼성전자주식회사 지자기 센서 오차 보정 장치 및 지자기 센서 오차 보정 방법

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101938564B1 (ko) * 2017-02-15 2019-01-15 주식회사 엠앤디 센서 교정 유무 판단 방법
KR101985458B1 (ko) * 2018-11-15 2019-06-03 주식회사 미래와도전 지진 센서의 오차를 예측하는 지진 분석 장치 및 오차 예측 방법
KR102357109B1 (ko) * 2021-02-17 2022-02-08 주식회사 비엘 시공중 터널 막장면 상태평가 시스템

Also Published As

Publication number Publication date
US10060743B2 (en) 2018-08-28
US20170059318A1 (en) 2017-03-02

Similar Documents

Publication Publication Date Title
KR101698682B1 (ko) 지자기 센서의 출력값을 보정하는 방법 및 장치
Tedaldi et al. A robust and easy to implement method for IMU calibration without external equipments
JP5237723B2 (ja) 動的に較正されるセンサデータと、ナビゲーションシステム内の繰り返し拡張カルマンフィルタとを使用する、ジャイロコンパスの整合用のシステム及び方法
US9031805B2 (en) Geomagnetic field measurement device, offset determination method, and computer readable recording medium therefor
KR101485142B1 (ko) 자체 보정 멀티-자기력계 플랫폼을 위한 시스템 및 방법
CN103941309B (zh) 地磁传感器校准设备及其方法
CN108426571B (zh) 一种电子罗盘本地实时校准方法及装置
US11248932B2 (en) Method for calibrating magnetometers of an object
CN103575293B (zh) 一种磁力计方向角校正方法及磁力计
US20140214364A1 (en) Method of estimating offset of magnetic sensor
KR20130143576A (ko) 중력 기준계에서의 장치에 부착된 움직임 센서 및 자력계의 측정치를 사용하여 장치의 요각을 추정하는 장치 및 방법
KR101106048B1 (ko) 센서오차의 작동 중 자동교정 방법과 이를 이용한 관성항법장치
CN107870001A (zh) 一种基于椭球拟合的磁力计校正方法
US20180313933A1 (en) Method of tracking a position of a magnet by differential measurement
CN110319851B (zh) 传感器的校正方法、装置、设备及存储介质
JP2017207456A (ja) 姿勢推定装置、姿勢推定方法、制御プログラム、および記録媒体
KR20140025319A (ko) 근 자기장의 동적 추적 및 보상을 위한 장치 및 방법
US10325005B2 (en) Apparatuses and methods for calibrating magnetometer attitude-independent parameters
US9863867B2 (en) Automatically updating hard iron and soft iron coefficients of a magnetic sensor
US11709056B2 (en) Method and device for magnetic field measurement by magnetometers
US20150260543A1 (en) Background calibration
CN107356786B (zh) 加速度计的校准方法和装置、计算机可读存储介质
CN106931965B (zh) 一种确定终端姿态的方法及装置
JP2005195376A (ja) 方位角計測方法および方位角計測装置
KR101928231B1 (ko) 오프셋 산출 장치 및 이것을 사용한 방위각 센서

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20191203

Year of fee payment: 4