Micro azimuth-level detector based on micro electro-mechanical systems and a method for determination of attitude

Information

  • Patent Grant
  • 6813584
  • Patent Number
    6,813,584
  • Date Filed
    Tuesday, March 12, 2002
    22 years ago
  • Date Issued
    Tuesday, November 2, 2004
    20 years ago
Abstract
A system and method for determining attitude using triaxial micro electro-mechanical accelerometer-magnetometer sensors. The triaxial accelerometer-magnetometer sensors measure the triaxial components of the local gravity and magnetic field along three sensing axes of the sensors. The misalignment of the components due to nonorthogonal triaxial assembling of the sensors is estimated via performing three optical alignments and is further compensated in the sensing outputs. A micro-cube with three orthogonal coated mirrors is fixed on the circuit board of the system for implementing the optical alignments. The compensated sensing components are processed to figure out the attitude using Orientation Cosine Conversion of the local gravity and magnetic field. The computation for determining the attitude is based on the inverse tangent of the ratio of the processed components, which is more precise than the computation using the absolute values of the components.
Description




FIELD OF THE INVENTION




This present invention is generally subject to the field of instrument & measurement technology based on Micro Electro-mechanical Systems (MEMS), especially relates to a novel design of micro navigation system.




BACKGROUND OF THE INVENTION




Navigation is actually to determine position and attitude for a vehicle by using kinds of sensors. The conventional navigation schemes can be cataloged as: (1) Inertial navigation system, uses gyroscopes and accelerometers to measure rate of rotation and acceleration respectively. However gyroscope is always expensive, has relative big size, slow start-up and long-term drift. (2) Imaging navigation system is composed of camera, units for processing and recognizing, etc. But the imaging processing is complicated and difficult to be implemented in real time. (3) Global position system (GPS), comprises 24 satellites, receiver and signal-processor. It highly depends on the outside information from satellites.




Compared with the conventional methods, the burgeoning MEMS navigation equipments on the basis of micro silicon-based sensors can put miniaturization, practicability and intelligence together. The MEMS-based navigation systems take advantages of low cost, small size, lightweight and adequate performances with respect to the conventional systems.




MEMS is based on integration of micro-electron technology and micro-mechanics technology. A main representative of MEMS-based attitude system is Honeywell Digital Compass HMR3000. It is a tilt compensated compass that uses tilt sensor for enhanced performance up to ±40° tilt range. Its composing and block diagram is shown in FIG.


1


: liquid filled tilt sensors


11


determine inclination angles, the heading (or azimuth) is calculated


13


from the tilt information and the three magnetic field components detected by triaxial magnetometers


12


, the angles then can be communicated with PC


15


through RS-232


14


. HMR3000 measures 0°˜360° heading within ±40° tilt, and the heading accuracy is ±1° for level, but ±2° for incline. The tilt range is only ±40°, and the tilt accuracy is ±0.4° for 0°˜20°, and ±0.6° when exceed ±20°.




The disadvantages of HMR3000 are:




(1) The tilt range is small, only covers ±40°.




(2) The measurement accuracy will decline with bevel rising.




SUMMARY OF THE INVENTION




The present invention has a novel structure and attitude calculation algorithm to avoid the foregoing disadvantages.




The object of the present invention is to design a novel MEMS attitude system with nonrestrictive tilt range and excellent performance. The inventer designs architecture, composing of three-axis orthogonal accelerometers and three-axis orthogonal magnetometers, and schemes out the optimum attitude calculation algorithm for getting higher accuracy and higher reliability.




During sensors assembly, the directions of the tilt sensor (accelerometers) must be carefully aligned with the corresponding magnetic axes, and the three axes in the same orthogonal placement must be mutually orthogonalized strictly. It is critical for determining the measurement accuracy. Aiming at this assembling problem, the present invention constructs an optical alignment approach against the misalignment.




The present integrative system can attain overall attitude (−90°˜90° for pitch, −180°˜180° for roll, 0°˜360° for heading) with high accuracy in any tilt, and possesses characteristics as solid-state structure, high accuracy, small size, lightweight, low consumption, fast start-up and low cost.




The invention is about a novel micro azimuth-level detector based on MEMS, which comprises multi-sensors, A/D converter for converting the analog signals of the sensors to the digital signals, microprocessor for calculation of the digital signals, RS-232 for connecting the microprocessor with PC, and further comprises signal-processing in the microprocessor and Operation-Display software in PC, wherein multi-sensors are triaxial silicon accelerometers with the processing circuit and triaxial silicon magnetometers with the processing circuit.




The foregoing the signal-processing and the Operation-Display software include signal-syncretizing, attitude computation and graphical Operation-Display interface.




For overcoming the misalignment in assembling, a micro-cube with three micro-mirrors in orthogonal directions using for optical alignment is added to the system structure.




The signal-syncretizing is combination of signal-filter and electric/physical (E/P) signal conversion.




For signal-filter, 2nd-order FIR filter is employed:








y


(


n


)=


h




0




x


(


n


)+


h




1




x


(


n−


1)






where x(n),x(n−1) are the output of sensor (input of the filter) and its backward shift, h


0


, h


1


are constant coefficients, y(n) is the output of the filter.




For E/P signal-conversion, that is actually conversion from electrical voltage v(n) to corresponding physical signal u(n):








u


(


n


)=(


v


(


n


)−


v




0


)/


k




v








where v


0


is zero-point, and k


v


is scale factor. v


0


and k


v


are obtained in advance by calibration.




The attitude computation is as following:




Establish two reference coordinate systems: the geographic coordinate system (GCS) shown in FIG.


2


(


a


) and the vehicle coordinate system (VCS) shown in FIG.


2


(


b


). Defining the three orthogonal axes of the GCS as North21-East22-Ground23 (N-E-G), and the three orthogonal axes of the VCS as X-Y-Z (X


31


refers to the forward direction of the vehicle, Y


32


refers to the left-right direction of the vehicle, and Z


33


refers to downward direction of the vehicle).




The rotations of the VCS relative to the GCS represent the attitude status of the vehicle. The terms ROLL (γ), PITCH (θ) and HEADING (ψ) are commonly used in aviation to represent attitude: ROLL refers to the rotation round the X-axis, PITCH refers to the rotation around the Y-axis, and HEADING refers to the rotation round the G-axis. The two representations of the same vector in GCS and VCS can be transformed each other through the Orientation Cosine Matrix C


n




b


.










C
n
b

=

[




cos


(

N
,
X

)





cos


(

E
,
X

)





cos


(

G
,
X

)







cos


(

N
,
Y

)





cos


(

E
,
Y

)





cos


(

G
,
Y

)







cos


(

N
,
Z

)





cos


(

E
,
Z

)





cos


(

G
,
Z

)





]







=


[

T
ij

]


3
×
3















The element T


ij


of C


n




b


is composed of pitch, roll and heading. The three magnetometers are assembled along the X-Y-Z axes to form an orthogonal placement. And the assembling placement is also for the three accelerometers. The independent X, Y, and Z magnetic and gravity readings [x


M


,y


M


,z


M


]


T


and [x


g


,y


g


,z


g


]


T


, which figure the independent X, Y, Z components of the earth's magnetic and gravity field, can be transformed back to the geographic coordinates by applying the Cosine Matrix C


n




b


. Further considering about errors [e


Mx


e


My


e


Mz


]


T


and [e


gx


e


gy


e


gz


]


T


, the following equations are attained (where H is magnetic hall effect, β is magnetic dip angle, f


g


is gravity acceleration):







[




x
M






y
M






z
M




]

=




HC
n
b



[




cos





β





0





sin





β




]


+

[




e
Mx






e
My






e
Mz




]




{








H
×

(



T
11


cos





β

+


T
13


sin





β


)


+

e
Mx


=

x
M









H
×

(



T
21


cos





β

+


T
23


sin





β


)


+

e
My


=

y
M









H
×

(



T
31


cos





β

+


T
33


sin





β


)


+

e
Mz


=

z
M









[




x
g






y
g






z
g




]

=




C
n
b



[



0




0





f
g




]


+

[




e
gx






e
gy






e
gz




]




{





x
g

=



T
13



f
g


+

e
gx









y
g

=



T
23



f
g


+

e
gy









z
g

=



T
33



f
g


+

e
gz





















In the equations above, the errors can be eliminated through signal filter and error compensation.




From the foregoing equations, the three attitude angles are figured out via inverse tangent. The procedures of computation are:




Step one: calculating roll γ according to Table 1 and value of γ


ref








γ


ref


=arctg(


y




g




/z




g


),



















TABLE 1











z


g






y


g






γ




Quadrant













→0




+




  90°








→0









−90°







+




+




γ


ref






(0°, 90°)







+









γ


ref






(−90°, 0°)












+




γ


ref


+180°




(90°, 180°)

















γ


ref


−180°




(−180°, −90°)















Step two: calculating pitch θ according to






θ=−arctg(


x




g


·cosγ/


z




g


),






Step three: and calculating heading ψ according to Table 2, ψ


ref


, Xh and Yh,






Xh
=



x
m


cos





θ

+


(



y
m


sin





γ

+


z
m


cos





γ


)


sin





θ






Yh
=



z
m


sin





γ

-


y
m


cos





γ








ψ
ref

=


arc






tg


(




z
m


sin





γ

-


y
m


cos





γ





x
m


cos





θ

+


(



y
m


sin





γ

+


z
m


cos





γ


)


sin





θ



)



=

arc






tg


(

Yh
Xh

)





,























TABLE 2











Xh




Yh




ψ




Quadrant













→0




+




  90°








→0









270°







+




+




ψ


ref






(0°, 90°)







+









360° +ψ


ref






(270°, 360°)












+




180° +ψ


ref






(90°, 180°)

















180° +ψ


ref






(180°, 270°)















For eliminating the misalignment in assembling, a micro-cube is fixed on the circuit board. The micro-cube is a small right angle cube with three orthogonal coated mirrors, which are used for reflecting beams. It can be bought from Industrial Optics Company. It is a kind of externally coated cubic prism, which has an enhanced aluminum coating on the three orthogonal side-faces. The three beams upright the three orthogonal mirrored surfaces form a 3-D orthogonal beam adjustment X-Y-Z. X-Y-Z is defined as the datum frame of the groupware. The groupware is actually circuit board where three-axis accelerometers and three-axis magnetometers are installed. Define X′-Y′-Z′ as the frame of the accelerometers, and X″-Y″-Z″ as the frame of the magnetometers. The optical alignment is carried out through three steps below:




Step one: Three optical alignments




(a) Align X-Y-Z to North (magnetic north)-East-Ground by beam-based alignment. The independent components of the gravity field along N(X)-E(Y)-G(Z) axis are [0,0,f


g


]


T


, the independent components of the earth's magnetic field along N(X)-E(Y)-G(Z) axis are [H


N


,0,H


G


]


T


. The relation between the three accelerometers' readings [X′


a


,Y′


a


,Z′


a


]


T


and [0,0,f


g


]


T


is:






{





X
a


=


f
g

·

cos


(


X


,
Z

)










Y
a


=


f
g

·

cos


(


Y


,
Z

)











Z
a


=


f
g

·

cos


(


Z


,
Z

)




,




&AutoRightMatch;










The relation between the three magnetometers' readings [X″


a


,Y″


a


,Z″


a


]


T


and [H


N


,0,H


G


]


T


is:






{





X
a


=



H
N

·

cos


(


X


,
X

)



+


H
G

·

cos


(


X


,
Z

)











Y
a


=



H
N

·

cos


(


Y


,
X

)



+


H
G

·

cos


(


Y


,
Z

)













Z
a


=



H
N

·

cos


(


Z


,
X

)



+


H
G

·

cos


(


Z


,
Z

)





;









&AutoRightMatch;










(b) Align X-Y-Z to East-Ground-North(magnetic north) by beam based alignment. Here the relation between the three accelerometers' readings [X′


b


,Y′


b


,Z′


b


]


T


and the gravity XYZ-stage components [0,f


g


,0]


T


is:






{





X
b


=


f
g

·

cos


(


X


,
Y

)










Y
b


=


f
g

·

cos


(


Y


,
Y

)











Z
b


=


f
g

·

cos


(


Z


,
Y

)




,




&AutoRightMatch;










The relation between the three magnetometers' readings [X″


b


,Y″


b


,Z″


b


]


T


and the magnetic XYZ-stage components [0,H


G


,H


N


]


T


is:






{





X
b


=



H
N

·

cos


(


X


,
Z

)



+


H
G

·

cos


(


X


,
Y

)











Y
b


=



H
N

·

cos


(


Y


,
Z

)



+


H
G

·

cos


(


Y


,
Y

)













Z
b


=



H
N

·

cos


(


Z


,
Z

)



+


H
G

·

cos


(


Z


,
Y

)





;









&AutoRightMatch;










(c) Align X-Y-Z to Ground-North(magnetic north)-East by beam based alignment. The relation between the three accelerometers' readings [X′


c


,Y′


c


,Z′


c


]


T


and the gravity XYZ-stage components [f


g


,0,0]


T


is:






{





X
c


=


f
g

·

cos


(


X


,
X

)










Y
c


=


f
g

·

cos


(


Y


,
X

)











Z
c


=


f
g

·

cos


(


Z


,
X

)




,




&AutoRightMatch;










The relation between the three magnetometers' readings [X″


c


,Y″


c


,Z″


c


]


T


and the magnetic XYZ-stage components [H


G


,H


N


,0]


T


is:






{





X
c


=



H
N

·

cos


(


X


,
Y

)



+


H
G

·

cos


(


X


,
X

)











Y
c


=



H
N

·

cos


(


Y


,
Y

)



+


H
G

·

cos


(


Y


,
X

)












Z
c


=



H
N

·

cos


(


Z


,
Y

)



+


H
G

·


cos


(


Z


,
X

)






.












&AutoRightMatch;










Step two: Determining the misaligned angles




Using the relation equations established in step one, figure out the Orientation Cosine Matrix Mg between the datum frame X-Y-Z and the acceleration frame






Mg
=

[




cos


(


X


,
X

)





cos


(


X


,
Y

)





cos


(


X


,
Z

)







cos


(


Y


,
X

)





cos


(


Y


,
Y

)





cos


(


Y


,
Z

)







cos


(


Z


,
X

)





cos


(


Z


,
Y

)





cos


(


Z


,
Z

)





]











and the Orientation Cosine Matrix Mm between X-Y-Z and the magnetic frame X″-Y″-Z″,






Mm
=


[




cos


(


X


,
X

)





cos


(


X


,
Y

)





cos


(


X


,
Z

)







cos


(


Y


,
X

)





cos


(


Y


,
Y

)





cos


(


Y


,
Z

)







cos


(


Z


,
X

)





cos


(


Z


,
Y

)





cos


(


Z


,
Z

)





]

.











Step three: Compensating the misaligned angles




When using the azimuth-level detector to measure attitude, always switch the outputs of the sensors into the orthogonal tuning X-Y-Z using the determined Orientation Cosine Matrix Mg and Mm.








[




x
g






y
g






z
g




]

=


Mg

-
1




[





x


g







y


g







z


g




]



,






[




x
M






y
M






z
M




]

=



Mm

-
1




[





x


M







y


M







z


M




]


.












Where, [x′


g


,y′


g


,z′


g


]


T


, [x′


M


,y′


M


,z′


M


]


T


are respectively readings of three-axis accelerometers and three-axis magnetometers. And [x


g


,y


g


,z


g


]


T


,[x


M


,y


M


,z


M


]


T


are respectively the orthogonal tuning XYZ-stage components of the corresponding sensors.




The characteristics and effects are:




The present invention has largely improved on the composing, and the attitude computation algorithm. Three-axis silicon accelerometers make up of the inclinometer as the accelerometer can provide three independent orthogonal components of the earth's gravity acceleration. The accelerometer-based inclinometer acts on not only the tilt-compensation, but also measurement for any pitch and roll. Through the Orientation Cosine Matrix C


n




b


, the relation equations between GCS and VCS can be established, further export the three attitude angles (pitch, roll, heading). The micro attitude detector of the invention is an integration of multi-MEMS-sensors, such as MEMS accelerometers and MEMS magnetometers. By taking the benefits of MEMS, the system takes on solid-state structure, small size, lightweight, low consumption, fast start-up and low cost.




Through the said inverse tangent algorithm of the invention are obtained high attitude accuracy in the full-scale range and excellent robust for withstanding effect field change. By optical alignment, the misalignment of each sensor axis deflective from the datum frame X-Y-Z are detected and compensated so as to make the attitude accuracy further largely enhanced.




As will be realized, the invention is capable of other and different embodiments, and its several details are capable of modifications in various obvious respects, all without departing from the invention. Accordingly, the drawings and description thereof are to be regarded as illustrative in nature, and not as restrictive.











BRIEF DESCRIPTION OF THE DRAWINGS





FIG. 1

is the composition and block diagram of HMR3000.





FIG. 2

is the view of vehicle coordinate system and geographic coordinate system.





FIG. 3

is the composition and block diagram of the present invention.





FIG. 4

is the software flow chart of the present invention.





FIG. 5

is the circuit diagram of best mode for carrying out the present invention.





FIG. 6

is the schematic view of Operation-Display interface.











BEST MODES FOR CARRYING OUT THE INVENTION




The hardware composition and the block diagram of the present invention are shown in FIG.


3


. The hardware includes: a micro-cube with three orthogonal mirrors


49


, three-axis MEMS accelerometers


42


and the processing circuit, three-axis MEMS magnetometers


43


and the processing circuit, magnetic Set/Reset pulse circuit


48


(reference from Honeywell), voltagestabilizing circuit


41


, A/D converter


44


, microprocessor


45


, RS-232


46


and PC


47


. The micro-cube with three orthogonal coated mirrors is for accomplishing optical alignment to deal with axis-misalignment. The three-axis accelerometers measure the three independent components of gravity acceleration and provide signals indicative of the measurement. The three-axis magnetometers measure the three independent components of the earth's magnetic field and provide signals indicative of the measurement. Magnetic Set/Reset pulse circuit is for eliminating the effects of past magnetic effects. The voltagestabilizing circuit is to convert 7˜12 volt power supply to 5 volt with more stabilization. A/D converter is to convert the analog signals to the digital signals. Then the digital signals can readily interface to microprocessor for signal collection, conversion, processing and attitude computation. The computed attitude information is transferred into PC through RS-232.




The software flow chart is shown in FIG.


4


. The software includes: data collection


51


, A/D conversion


52


, signal-syncretizing


53


(signal-filter, E/P signal-conversion), misalignment compensation


54


, attitude computation


55


, serial output


56


, and Operation-Display interface


50


in PC. The signals flow is as: the sensors' electrical signals are detected, A/D converted, filtered, further converted to physical signals, and through misalignment compensation processing, then computed to get attitude angles (the above steps all proceed in microprocessor). The attitude angles are exported into PC through RS-232. In PC the attitude angles are read from a serial connection port


57


, then conditionally smooth-processed


58


and


59


, that is to apply smoothing to small variations and little or no smoothing for large change in attitude outputs. The smooth-processing


59


is carried out as:






If, |angle(


n


)−angle(


n−


1)|<


ss










Then, attitude(


n


)=


a




0


·angle(


n


)+


a




1


·angle(


n−


1)






Where, angle(n),angle(n−1) are the current attitude outputs and their backward shift, attitude(n) is the smoothed attitude, ss is threshold value, a


0


,a


1


are constant coefficients, 0<a


0


<0.2, a


1


=1−a


0


.




The processed attitude results are finally displayed on Operation-Display interface


60


.





FIG. 5

depicts the circuit diagram of this mode. The power supply should be stabilized through


78


D


05


U


15


before connect to the circuit. The C


17


, C


18


and C


19


are three capacitors around


78


D


05


. The function of C


17


is to stabilize input-voltage, C


18


and C


19


are used for smoothing fluctuation of out-voltage. Put the power supply switch on, then the microprocessor P


89


C


51


RC+U


12


with resistor R


19


, capacitors C


23


, C


24


, C


25


and oscillator Y


1


(Y


1


, C


24


and C


25


provide clock-signal for microprocessor) begins to read the outputs of the accelerometers two ADXL


202


U


1


and U


2


(C


11


and C


14


are filter-capacitors of supply voltage. Capacitors C


12


, C


13


and C


15


are added at output pins to implement low-pass filtering for antialiasing and noise reduction. The period of the power supply is set for the channels of accelerometers by resistors R


14


, R


15


, R


16


and R


17


), and the two-axis magnetometer HMC


1002


U


4


and single-axis magnetometer HMC


1001


U


5


(Resistors R


6


, R


7


, R


8


, R


10


, R


11


and R


12


are used to trim switchpoint), through 12 bit A/D converter MAX


197


U


13


with resistor R


20


, capacitors C


31


, C


32


, C


33


, C


34


and C


35


. C


31


between the CLK pin of U


13


and GND is necessary for selecting internal clock mode of A/D. R


20


, C


32


and C


33


make up of the filter circuit. Thereinto the accelerometer signals are Rail-to-Rail treated through MAX


4044


U


3


. And the magnetometer signals are amplified through amplifiers MAX


4194


U


6


, U


7


, U


8


respectively with gain-setting resistors R


9


, R


13


and R


5


. MAX


6025


U


9


with capacitor C


16


and resistor R


18


is to supply 2.5V reference-voltage for the amplifiers. Wherein C


16


will reduce the amount of overshoot and assist the circuit's transient response, R


18


is to limit the supply current. Data output serially using the RS-232 standard MAX


232


CSE U


14


is transferred to most personal computers. The MAX232CSE has two internal charge-pumps that convert +5V to ±10V for RS-232 driver operation. Wherein the first converter uses capacitor C


29


and C


26


to double the +5V input to +10V, and the second converter uses capacitor C


30


and C


27


to invert +10V to −10V. Capacitor C


28


is used for stabilizing supply-voltage.




In Set/Reset pulse circuit


100


encircled by dashed line (reference from Honeywell), the SET and RESET signals are generated from the microprocessor and are used to control the NPN General Purpose Amplifier


2


N


3904


Y


2


and HEXFET Power MOSFET IRF


7307


U


10


, as well as resistors R


1


, R


2


, R


3


, R


4


, and capacitors C


6


, C


7


, C


8


, to create Set/Reset current pulse. The magnitude of the Set/Reset current pulse depends on the output voltage of MAX


662


U


11


, which regulates 5V power supply to the claimed 12V voltage. MAX


662


uses internal charge pumps and external capacitors C


1


, C


2


to generate +12V. C


3


is the input bypass capacitor and C


5


is output filter capacitor.




The schematic view of Operation-Display interface is shown in FIG.


6


. In the top left side is presented Pitch&Roll graphical display


61


and Heading graphical display


62


, and pitch, roll and heading numerical display


63




64




65


is listed below the graph. Below attitude numerical display, a pair of buttons “Display”


66


and “Stop”


67


are used to activate and stop attitude display; “Update Frequency”


68


is used to set the specific update frequency of attitude display. “Declination Angles”


69


is used to set the declination between the forward direction of the detector and the true geographical north; “Geography Heading”


70


and “Tare Heading”


71


are respectively used to display the true geographic heading and display the relative heading from any initial heading. Furthermore, in the top right side of interface, included “Set/Reset” function region


80


(Set/Reset technique is for eliminating the effects of past magnetic effects, reference from Honeywell), wherein “Start”


81


and “End”


82


buttons are used to activate and finish Set/Reset, “Check”


83


is used to check up if the Set/Reset is success or failure; In the under right side of interface, included Calibration function region


90


for calibrating zero-value v


0


and scale factor k


v


of three-axis magnetometers, “Start”


91


and “End”


92


are used to startup and stop the calibration, the below six windows show calibration parameters (the left three windows


93




94




95


show the zero values and the right three


96




97




98


show the scale factors corresponding to X, Y, and Z magnetic axis), “Default”


99


is to default the calibration parameters to be aforehand-setting values. For the calibration parameters of the accelerometers, the parameters are fixed and stored in microprocessor beforehand as they are not affected by external environment.




The detector followed by this mode can measure a full-scale range of attitude, that is −90°˜+90° for pitch, −180°˜+180° for roll, 0°˜360° for heading. The attitude accuracy is excelled 1° for the full-scale range. It can also be applied in uncertain effect field to detect attitude of vehicle. By optical alignment, the mode can further improve the performances in measurement accuracy that is 0.2° for tilt and 0.5° for heading.




After reading the foregoing specification, one of ordinary skill will be able to affect various changes, substitutions of equivalents and various other aspects of the invention as broadly disclosed herein. It is therefore intended that the protection granted hereon be limited only by the definition contained in the appended claims and equivalents thereof.



Claims
  • 1. A Micro Electro-mechanical System (MEMS) based micro azimuth-level detector comprising:multi-sensors, an analog/digital (A/D) converter, a microprocessor for signal processing, an RS-232 electrical interface, and a PC having an Operation-Display software, wherein the multi-sensors comprise three accelerometers and three magnetometers, the three accelerometers are assembled along three orthogonal axes and having a processing circuit, and the three magnetometers are assembled along three orthogonal axes and have a processing circuit, a micro-cube with three orthogonal coated mirrors for optical alignment, and wherein analogue signals from the multi-sensors are converted into digital signals by the A/D converter, sent to and processed in the microprocessor, and further transmitted via the RS-232 electrical interface to be shown in the PC.
  • 2. The micro azimuth-level detector of claim 1, wherein the digital signals are processed in the microprocessor by signal syncretizing and attitude computation, and sent via the RS232 to be shown in a graphical Operation-Display Software in the PC.
  • 3. A method for determining attitude comprisingobtaining electrical/analog signals from sensors, converting the electrical/analog signals into digital signals by an A/D converter, filtering the digital signals and processing the digital signals by electric/physical conversion, eliminating misalignment in the processed digital signals by misalignment compensation based on optical alignment, computing attitude angles including pitch, roll and heading from the processed signals by Orientation Cosine Conversion of the earth's gravity and magnetic field and inverse tangent method, conditionally smooth-processing the computed attitude angles by applying a large amount of smoothing to small variations and little or no smoothing for large change in attitude outputs, displaying the conditionally smooth-processed attitude on an Operation-Display interface.
  • 4. The method for determining attitude according to claim 3, wherein the processed triaxial gravity signals xg, yg, and zg, and the processed triaxial magnetic signals xm, ym, and zm are obtained through the sensors and processed, and the attitude angles pitch, roll and heading are computed by(1) calculating roll γ according to Table 1 and value of γref γref=arcg(yg/zg), TABLE 1 =zgygγQuadrant→0+  90°→0−−90°++γref(0°, 90°)+−γref(−90°, 0°)−+γref +180°(90°, 180°)−−γref −180°(−180°, −90°)(2) calculating pitch θ according to θ=−arctg(xg·cos γ/zg) (3) and calculating heading ψ according to Table 2, ψref, Xh and Yh, Xh=xm⁢cos⁢ ⁢θ+(ym⁢sin⁢ ⁢γ+zm⁢cos⁢ ⁢γ)⁢sin⁢ ⁢θYh=zm⁢sin⁢ ⁢γ-ym⁢cos⁢ ⁢γψref=arc⁢ ⁢tg⁢(zm⁢sin⁢ ⁢γ-ym⁢cos⁢ ⁢γxm⁢cos⁢ ⁢θ+(ym⁢sin⁢ ⁢γ+zm⁢cos⁢ ⁢γ)⁢sin⁢ ⁢θ)=arc⁢ ⁢tg⁢(YhXh),TABLE 2 =XhYhψQuadrant→0+  90°→0−270°++ψref(0°, 90°)+−360° +ψref(270°, 360°)−+180° +ψref(90°, 180°)−−180° +ψref(180°, 270°).
  • 5. The method for determining attitude according to claim 3, wherein the misalignment is eliminated by(1) defining the orthogonal beam adjustment X-Y-Z of micro-cube as a datum orthogonal frame of groupware, (2) defining X′-Y′-Z′ as a frame of the accelerometers in the orthogonal directions, and X″-Y″-Z″ as a frame of the magnetometers in the orthogonal directions, (3) Aligning the datum frame X-Y-Z to North(magnetic north)-East-Ground by beam-based alignment, a relation between three accelerometers' readings [X′a,Y′a,Z′a]T and independent components [0,0, fg]T of the gravity field along N(X)-E(Y)-G(Z) being {Xa′=fg·cos⁡(X′,Z)Ya′=fg·cos⁡(Y′,Z)Za′=fg·cos⁡(Z′,Z); ⁢and&AutoRightMatch; a relation between three magnetometers' readings [X″a,Y″a,Z″a]T and independent components [HN,0,HG]T of the earth's magnetic field along N(X)-E(Y)-G(Z) being {Xa″=HN·cos⁡(X″,X)+HG·cos⁡(X″,Z)Ya″=HN·cos⁡(Y″,X)+HG·cos⁡(Y″,Z)⁢ Za″=HN·cos⁡(Z″,X)+HG·cos⁡(Z″,Z)⁢ ;⁢ &AutoRightMatch;(4) Aligning X-Y-Z to East-Ground-North(magnetic north); A relation between three accelerometers' readings [X′b,Y′b,Z′b]T and gravity XYZ-stage components [0, fg,0]T being {Xb′=fg·cos⁡(X′,Y)Yb′=fg·cos⁡(Y′,Y)Zb′=fg·cos⁡(Z′,Y); ⁢and&AutoRightMatch; a relation between three magnetometers' readings [X″b,Y″b,Z″b]T and magnetic XYZ-stage components [0, HG,HN]T being {Xb″=HN·cos⁡(X″,Z)+HG·cos⁡(X″,Y)Yb″=HN·cos⁡(Y″,Z)+HG·cos⁡(Y″,Y)⁢ Zb″=HN·cos⁡(Z″,Z)+HG·cos⁡(Z″,Y);&AutoRightMatch;(5) Aligning X-Y-Z to Ground-North(magnetic north)-East; a relation between three accelerometers' readings [X′c,Y′c,Z′c]T and gravity XYZ-stage components [fg,0,0]T being {Xc′=fg·cos⁡(X′,X)Yc′=fg·cos⁡(Y′,X)Zc′=fg·cos⁡(Z′,X); ⁢and&AutoRightMatch; a relation between three magnetometers' readings [X″c,Y″c,Z″c]T and magnetic XYZ-stage components [HG,HN,0]T being {Xc″=HN·cos⁡(X″,Y)+HG·cos⁡(X″,X)Yc″=HN·cos⁡(Y″,Y)+HG·cos⁡(Y″,X)Zc″=HN·cos⁡(Z″,Y)+HG·cos⁡(Z″,X);&AutoRightMatch;(6) computing Orientation Cosine Matrix Mg between the datum frame X-Y-Z and the acceleration frame X′-Y′-Z′, Mg=[cos⁡(X′,X)cos⁡(X′,Y)cos⁡(X′,Z)cos⁡(Y′,X)cos⁡(Y′,Y)cos⁡(Y′,Z)cos⁡(Z′,X)cos⁡(Z′,Y)cos⁡(Z′,Z)],(7) computing Orientation Cosine Matrix Mm between X-Y-Z and the magnetic frame X″-Y″-Z″, Mm=[cos⁡(X″,X)cos⁡(X″,Y)cos⁡(X″,Z)cos⁡(Y″,X)cos⁡(Y″,Y)cos⁡(Y″,Z)cos⁡(Z″,X)cos⁡(Z″,Y)cos⁡(Z″,Z)];(8) switching three accelerometers' readings [x′g,y′g,z′g]T and three magnetometers' readings' [x′M,y′M,z′M]T into the orthogonal datum XYZ-stage components [xg,yg,zg]T and [xM,yM,zM]T using Mg and Mm when detecting attitude [xgygzg]=Mg-1⁡[xg′yg′zg′],⁢[xMyMzM]=Mm-1⁡[xM′yM′zM′].
  • 6. The method for determining attitude according to claim 3, wherein conditionally smooth-processing to the computed attitude angles is performed by applying smoothing to small variations and little or no smoothing for large change in attitude outputs as:If, |angle(n)−angle(n−1)|<ss Then, attitude(n)=a0·angle(n)+a1·angle(n−1) Where, angle(n), angle(n−1) are the current attitude outputs and their backward shift, attitude(n) is the smoothed attitude, ss is threshold value, a0, a1 are constant coefficients.
Priority Claims (1)
Number Date Country Kind
01110135 A Mar 2001 CN
US Referenced Citations (10)
Number Name Date Kind
5090231 Gallagher Feb 1992 A
5333110 Al-Attar Jul 1994 A
5854843 Jacknin et al. Dec 1998 A
6208936 Minor et al. Mar 2001 B1
6243660 Hsu et al. Jun 2001 B1
6421622 Horton et al. Jul 2002 B1
6476610 Wiegert et al. Nov 2002 B1
6543146 Smith et al. Apr 2003 B2
6647352 Horton Nov 2003 B1
6725173 An et al. Apr 2004 B2
Non-Patent Literature Citations (7)
Entry
M.J. Caruso, “Applications of Magnetic Sensors for Low Cost Compass Systems”, http//www.ssec.honeywell.com/magnetic/datasheets/lowcost.pdf, Mar. 2000.*
60/230,567, An et al., “Digital signal processing method and system thereof for precision orientation measurements”, Sep. 2, 2000.*
“Honeywell Sensor Products: Digital Compass Module HMR 3000,” http://www.ssec.honeywell.com/magnetic/datasheets/hmr3000.pdf, Dec. 1999.
“Set/Reset Pulse Circuits For Magnetic Sensors AN-201,” http://www.ssec.honeywell.com/magnetic/datasheets/an201.pdf, May 1996.
“Magnetic Sensor Hybrid Application Circuit AN-202,” http://www.ssec.honeywell.com/magnetic/datasheets/an202.pdf, Sep. 1995.
“Compass Heading Using Magnetometers AN-203,” http://www.ssec.honeywell.com/magnetic/datasheets/an203.pdf, Jul. 1995.
“Getting Stable Output From HMR3000 Compass Module AN-208,” http://www.ssec.honeywell.com/magnetic/datasheets/an208.pdf, Nov. 1998.