TWI407102B - 三維空間運動感測方法 - Google Patents
三維空間運動感測方法 Download PDFInfo
- Publication number
- TWI407102B TWI407102B TW99102073A TW99102073A TWI407102B TW I407102 B TWI407102 B TW I407102B TW 99102073 A TW99102073 A TW 99102073A TW 99102073 A TW99102073 A TW 99102073A TW I407102 B TWI407102 B TW I407102B
- Authority
- TW
- Taiwan
- Prior art keywords
- angular velocity
- coordinate directions
- readings
- motion
- acceleration
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/03—Arrangements for converting the position or the displacement of a member into a coded form
- G06F3/033—Pointing devices displaced or positioned by the user, e.g. mice, trackballs, pens or joysticks; Accessories therefor
- G06F3/0346—Pointing devices displaced or positioned by the user, e.g. mice, trackballs, pens or joysticks; Accessories therefor with detection of the device orientation or free movement in a 3D space, e.g. 3D mice, 6-DOF [six degrees of freedom] pointers using gyroscopes, accelerometers or tilt-sensors
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Input By Displaying (AREA)
- User Interface Of Digital Computer (AREA)
Description
本發明係有關於一種感測裝置在三維空間運動的方法,尤指一種可除去裝置於被移動時伴隨之轉動的影響,並將裝置座標轉換至相對於重力的座標,以感測出裝置在三維空間運動的方法。
由於當如3D滑鼠之裝置被移動時,即使使用者並無意要轉動3D滑鼠,移動3D滑鼠的手掌仍會不自覺地些微轉動3D滑鼠,然而3D滑鼠的三軸加速度計非常靈敏,其加速度讀值包含受重力加速度影響的分量以及手部運動所產生的加速度分量,因此若只使用傳統的座標演算法計算三軸加速度計相對於重力座標的加速度讀值,將因為三軸加速度計讀值所包含的手部運動所產生的加速度分量而無法得到3D滑鼠在空間中正確的運動軌跡。因為傳統的演算法是假設三軸加速度計之讀值只受到重力的影響而沒有考慮到讀值中另包含手部運動所產生的加速度分量,因此當3D滑鼠並不如預期為靜態或是等速運動時,藉由傳統演算法所計算出的三軸加速度計相對於重力座標的讀值,將會因非重力的加速度的影響而產生誤差,導致無法精準計算出3D滑鼠相對於重力的運動軌跡。
請參照第1圖,第1圖係先前技術的一裝置102在螢幕上所呈現的指標軌跡。請參照第1圖,裝置102向左移動兼轉動,因為裝置102未將所產生之運動軌跡由裝置座標轉換至重力座標,導致裝置102射向螢幕100的指標軌跡會隨著裝置102的轉動而改變其移動方向,也就是說,指標軌跡並非只是向左移動而是向左上方移動。
台灣專利第200639406號“慣性感測輸入裝置”揭露一種慣性感測輸入裝置,此篇專利所描述之運動裝置與方法,無論於靜止或運動中皆無法進行相對於重力之座標轉換,使其輸出值為相對於裝置座標之運動量,非相對於重力之移動量,導致指標之移動方向會隨著裝置相對於重力之方位而改變。
台灣專利第200538751號“自由空間定位裝置和方法”揭露一種自由空間定位裝置和方法,此篇專利所描述之運動裝置與方法,只能在靜態或慢速移動的條件下,計算從本體參考坐標系到用戶參考坐標系之座標轉換,由於缺乏動態下的補償修正機制,而不可避免的無法進行運動過程中之座標系轉換工作,導致若在運動過程中出現裝置座標與自由空間座標之相對改變,將造成不正確之運動裝置輸出。
美國專利第5,898,421號“陀螺儀指針和方法(Gyroscopic pointer and method)”揭露一種陀螺儀指針和方法,此篇專利未能計算加速度計與角度感測器的相對關係,因此需要兩顆或以上之雙軸角度感測器才能運算座標轉換時的所需之資訊。
本發明揭露一種三維空間運動感測方法,包含讀取至少一加速度計在三個座標方向的複數個讀值;讀取至少一角速度感測器在二個座標方向的複數個讀值;利用該至少一加速度計在該三個座標方向的複數個讀值及該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置的運動旋轉半徑;利用該裝置的運動旋轉半徑與該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置在該二個座標方向的切線加速度;及利用該至少一加速度計在該二個座標方向的複數個讀值,該裝置在該二個座標方向的切線加速度,該至少一角速度感測器在該二個座標方向的複數個讀值,及該裝置的運動旋轉半徑,產生該裝置的運動輸出。
第2圖係說明一裝置202的示意圖。在裝置202內配置三軸加速度計204與二軸角速度感測器206。三軸加速度計204與二軸角速度感測器206在裝置202於三維空間中運動時,可記錄裝置202運動時的加速度與角速度。
請參照第3圖和第4圖,第3圖係本發明的一實施例所揭露的一裝置在三維空間運動感測的方法之流程圖,第4圖係說明裝置202的座標軸,而三軸加速度計204和二軸角速度感測器206座標軸和裝置202的座標軸相同。第3圖之方法係利用第2圖的裝置202說明,詳細步驟如下:步驟300:開始;步驟302:讀取三軸加速度計204在三個座標方向的複數個讀值A x
、A y
、A z
,以及讀取二軸角速度感測器206在二個座標方向的複數個讀值w x
、w z
;步驟304:利用三軸加速度計204在三個座標方向的複數個讀值A x
、A y
、A z
及二軸角速度感測器206在二個座標方向的複數個讀值w x
、w z
,產生裝置202的運動旋轉半徑 r
;步驟306:利用裝置202的運動旋轉半徑r
與二軸角速度感測器206在二個座標方向的複數個讀值w x
、w z
,產生裝置202的法線加速度A ry
及裝置202在二個座標方向的切線加速度A rx
、A rz
;步驟308:利用三軸加速度計204在二個座標方向的複數個讀值A x
、A z
及裝置202在二個座標方向的切線加速度A rx
、A rz
,產生重力加速度在二個座標方向的加速度分量G x
、G z
;步驟310:利用二個座標方向的重力加速度分量G x
、G z
,產生x
軸方向和重力加速度合成分量G xz
的夾角θ,其中步驟312:利用二軸角速度感測器206在二個座標方向的複數個讀值w x
、w z
,裝置202的運動旋轉半徑r
及夾角θ進行座標轉換,以產生相對於重力座標的橫向運動輸出與縱向運動輸出;步驟314:結束。
在步驟302,其中A x
為三軸加速度計204中的x
軸加速度計偵測到之加速度分量、A y
為三軸加速度計204中的y
軸加速度計偵測到之加速度分量、A z
為三軸加速度計204中的z
軸加速度計偵測到之加速度分量、w z
為二軸角速度感測器206中之z
軸角速度感測器偵測到之繞z
軸旋轉之角速度以及w x
為二軸角速度感測器206中之x
軸角速度感測器偵測到繞x
軸旋轉之角速度。
在步驟304,裝置202在以y
軸為法線之旋轉運動中,三軸加速度計204在x
、z
軸方向的讀值A x
、A z
會包含裝置202在切線方向的加速度及重力加速度的成分,三軸加速度計204在y軸方向的讀值A y
會包含裝置202在法線方向的加速度及重力加速度的成分;由x
、z
軸方向的兩個角加速度分別乘上裝置202的旋轉半徑r
便可得到裝置202在切線方向的加速度A rx
、A rz
。而藉由三軸加速度計204在三軸方向的讀值A x
、A y
、A z
、裝置202在切線方向的加速度A rx
、A rz
,裝置202在法線方向的加速度A ry
及重力加速度G
的關係便可取得裝置202的旋轉半徑r
如下:
式中Δt
為三軸加速度計204和二軸角速度感測器206的取樣週期;Δw z
係為在繞z
軸旋轉之角速度於Δt
的差值;Δw x
係為在繞x
軸旋轉之角速度於Δt
的差值;三軸加速度計204讀值A x
、A y
、A z
分別扣除切線加速度A rx
、A rz
與法線加速度A ry
後之平方總和等於重力加速度G
的平方。
(A x
-A rx
)2
+(A z
-A rz
)2
+(A y
-A ry
)2
=G2
(4)
將式(1)、(2)、(3)代入式(4)可得到式(5),
式(5)中除了裝置202的旋轉半徑r
之外皆為已知參數,故可得裝置202的旋轉半徑r
。
步驟306,利用步驟304得到的旋轉半徑r
與二軸角速度感測器206在二個座標方向的複數個讀值w x
、w z
,代入式(6)、(7)可得裝置202在二個座標方向的切線加速度A rx
、A rz
:
步驟308,裝置202於運動中,由於三軸加速度計204在x
軸及z
軸方向的讀值A x
、A z
包含裝置202在切線方向的加速度及重力加速度的成分,因此將三軸加速度計204之輸出數值A x 、A z
扣除切線加速度A rx
、A rz
後,便可得到重力加速度在x
軸及z
軸方向的分量G x
、G z
。G x
、G z
可由式(8)、(9)得到。
G x
=A x
-A rx
(8)
G z
=A z
-A rz
(9)
其中G x
為x
軸方向的重力加速度分量,G z
為z
軸方向的重力加速度分量。
請參照第5圖,第5圖顯示重力加速度分量G x
、G z
與重力加速度合成分量G xz
的關係。步驟310,利用x
軸及z
軸方向的重力加速度分量G x
、G z
,產生於裝置202運動過程中和重力加速度合成分量G xz
的夾角θ,其中
在步驟312中,利用夾角θ進行座標轉換成相對於重力之座標的橫向運動輸出與縱向運動輸出,而橫向運動輸出與縱向運動輸出可以位移的方式輸出、以速度的方式輸出、及/或位移與速度混合的方式輸出。式(10)、(11)所示係橫向運動輸出與縱向運動輸出為位移的狀況,式(12)、(13)所示係橫向運動輸出與縱向運動輸出為速度的狀況。
橫向運動(位移)輸出=ωz
×r
×Δt
×sinθ+ω x
×r
×Δt
×cosθ (10)
縱向運動(位移)輸出=ωz
×r
×Δt
×cosθ+ω x
×r
×Δt
×sinθ (11)
橫向運動(速度)輸出=ωz
×r
×sinθ+ω x
×r
×cosθ (12)
縱向運動(速度)輸出=ωz
×r
×cosθ+ω x
×r
×sinθ (13)
其中式(10)、(11)及/或式(12)、(13)的橫向運動輸出與縱向運動輸出係為三軸加速度計204的一個取樣週期內的運動輸出。
請參照第6圖,第6圖係說明利用第2圖的裝置202和第3圖揭露的方法在螢幕100上所呈現的指標軌跡之示意圖。請參照第6圖,三軸加速度計204與二軸角速度感測器206可在裝置202於三維空間中運動時,偵測裝置202的加速度與角速度,然後利用第3圖揭露的方法校正手掌不自覺地轉動裝置202而產生的誤差,完成校正後的裝置202可使發射至螢幕100上的指標隨著裝置202相對於重力座標系的運動方向而移動,以克服在螢幕100上的指標會隨著使用者的手腕不自覺轉動而錯誤地判斷裝置202之移動方向的現象,以增加使用上的準確度。除此之外,三軸加速度計204可為三個可輸出單軸方向加速度讀值的加速度計、一個可輸出二軸方向加速度讀值的加速度計搭配一個可輸出單軸方向加速度讀值的加速度計、及/或一個可輸出三軸方向加速度讀值的加速度計;二軸角速度感測器206可為二個可輸出單軸方向角速度讀值的角速度感測器、及/或一可輸出二軸方向角速度讀值的角速度感測器。
相較於先前技術,本發明在於能隨時計算裝置座標相對於重力座標的關係,在自由空間慣性運動裝置的應用上,可使其於運動中輸出相對於重力座標的運動方向,而非相對於裝置自身座標的運動方向,以避免運動中與重力方向之相對座標改變而造成不正常之運動輸出。再者,本發明能藉由加速度計與角度感測器間之相關性,修正手掌不自覺地轉動裝置而產生的誤差,因此只需要一顆三軸加速度計與一顆二軸角度感測器即可運算座標轉換所需之資訊,進而輸出裝置相對於重力之運動方向。
以上所述僅為本發明之較佳實施例,凡依本發明申請專利範圍所做之均等變化與修飾,皆應屬本發明之涵蓋範圍。
100...螢幕
102、202...裝置
204...三軸加速度計
206...二軸角速度感測器
300-314...步驟
G xz
...重力加速度於x
、z
平面之分量
G x
...於x
軸方向之重力加速度分量
G z
...於z
軸方向之重力加速度分量
θ...夾角
第1圖係說明先前技術的裝置在螢幕上所呈現的指標軌跡。
第2圖係本發明裝置的示意圖。
第3圖係第2圖裝置在三維空間運動感測的方法之流程圖。
第4圖係說明第2圖裝置之座標軸。
第5圖係說明第2圖裝置的二個座標方向和重力方向的夾角。
第6圖係說明利用第2圖的裝置和第3圖揭露的方法在螢幕上所呈現的指標軌跡。
300-314...步驟
Claims (8)
- 一種三維空間運動感測方法,包含:讀取至少一加速度計在三個座標方向的複數個讀值;讀取至少一角速度感測器在二個座標方向的複數個讀值;利用該至少一加速度計在該三個座標方向的複數個讀值及該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置的運動旋轉半徑;利用該裝置的運動旋轉半徑與該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置在該二個座標方向的切線加速度;及利用該至少一加速度計在該二個座標方向的複數個讀值,該裝置在該二個座標方向的切線加速度,該至少一角速度感測器在該二個座標方向的複數個讀值,及該裝置的運動旋轉半徑,產生該裝置的運動輸出。
- 如請求項1所述之方法,其中利用該至少一加速度計在該二個座標方向的複數個讀值,該裝置在該二個座標方向的切線加速度,該至少一角速度感測器在該二個座標方向的複數個讀值,及該裝置的運動旋轉半徑,產生該裝置的運動輸出包含:利用該至少一加速度計在該二個座標方向的複數個讀值及該裝置在該二個座標方向的切線加速度,產生該裝置在該二個座標方向的重力加速度分量;利用該二個座標方向的重力加速度分量,得到該二個座標方向的重力加速度分量之重力加速度合成分量與該二個座標方向中一座標方向的夾角;及利用該至少一角速度感測器在該二個座標方向的複數個讀值、該裝置的運動旋轉半徑及該夾角,產生該裝置的橫向運動輸出與縱向運動輸出。
- 如請求項2所述之方法,其中利用該至少一加速度計在該三個座標方向的複數個讀值及該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置的運動旋轉半徑係利用下列方程式:
- 如請求項3所述之方法,其中利用該裝置的運動旋轉半徑與該至少一角速度感測器在該二個座標方向的複數個讀值,產生該裝置在該二個座標方向的切線加速度係利用下列方程式:
- 如請求項4所述之方法,其中利用該至少一加速度計在該二個座標方向的複數個讀值及該裝置在該二個座標方向的切線加速度,產生該裝置在該二個座標方向的重力加速度分量係利用下列方程式:G x =A x -A rx ;及G z =A z -A rz ;其中:G x 為x 軸方向的重力加速度分量;及G z 為z 軸方向的重力加速度分量。
- 如請求項5所述之方法,其中利用該二個座標方向的重力加速度分量,得到該二個座標方向的重力加速度分量之重力加速度合成分量與該二個座標方向中一座標方向的夾角係利用下列方程式:
- 如請求項6所述之方法,其中利用該至少一角速度感測器在該二個座標方向的複數個讀值、該裝置的運動旋轉半徑及該夾角,產生該裝置的該橫向運動輸出與該縱向運動輸出包含下列方程式:該橫向運動輸出=ω z ×r ×Δt ×sinθ+ω x ×r ×Δt ×cosθ;及該縱向運動輸出=ω z ×r ×Δt ×cosθ+ω x ×r ×Δt ×sinθ。
- 如請求項6所述之方法,其中利用該至少一角速度感測器在該二個座標方向的複數個讀值、該裝置的運動旋轉半徑及該夾角,產生該裝置的該橫向運動輸出與該縱向運動輸出包含下列方程式:該橫向運動輸出=ω z ×r ×sinθ+ω x ×r ×cosθ;及該縱向運動輸出=ω z ×r ×cosθ+ω x ×r ×Δt ×sinθ。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW99102073A TWI407102B (zh) | 2010-01-26 | 2010-01-26 | 三維空間運動感測方法 |
US12/748,461 US8797262B2 (en) | 2010-01-26 | 2010-03-29 | Method of sensing motion in three-dimensional space |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW99102073A TWI407102B (zh) | 2010-01-26 | 2010-01-26 | 三維空間運動感測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201126166A TW201126166A (en) | 2011-08-01 |
TWI407102B true TWI407102B (zh) | 2013-09-01 |
Family
ID=44308580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW99102073A TWI407102B (zh) | 2010-01-26 | 2010-01-26 | 三維空間運動感測方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8797262B2 (zh) |
TW (1) | TWI407102B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI637344B (zh) * | 2016-02-18 | 2018-10-01 | 緯創資通股份有限公司 | 空間繪畫計分方法、用於進行空間繪畫計分的裝置及系統 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8555205B2 (en) * | 2010-10-08 | 2013-10-08 | Cywee Group Limited | System and method utilized for human and machine interface |
KR20130127991A (ko) * | 2010-12-01 | 2013-11-25 | 꼼미사리아 아 레네르지 아토미끄 에뜨 옥스 에너지스 앨터네이티브즈 | 이동체의 경로를 추정하는 방법 및 시스템 |
TWI447356B (zh) * | 2011-05-27 | 2014-08-01 | Sitronix Technology Corp | Direction sensing device |
US10652696B2 (en) * | 2014-07-30 | 2020-05-12 | Trusted Positioning, Inc. | Method and apparatus for categorizing device use case for on foot motion using motion sensor data |
WO2017132563A1 (en) * | 2016-01-29 | 2017-08-03 | Baylor Research Institute | Joint disorder diagnosis with 3d motion capture |
CN116048296B (zh) * | 2022-08-11 | 2024-04-02 | 荣耀终端有限公司 | 电子设备及其轨迹生成方法和介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200538751A (en) * | 2004-04-30 | 2005-12-01 | Hillcrest Lab Inc | Free space pointing devices and methods |
TW200639406A (en) * | 2005-05-13 | 2006-11-16 | Ind Tech Res Inst | Initial sensing input apparatus |
US20070214886A1 (en) * | 2006-03-15 | 2007-09-20 | Qualcomm Incorporated | Sensor-based orientation system |
TW200825867A (en) * | 2006-12-12 | 2008-06-16 | Ind Tech Res Inst | Inertial input apparatus with six-axial detection ability and the opearting method thereof |
US20090303204A1 (en) * | 2007-01-05 | 2009-12-10 | Invensense Inc. | Controlling and accessing content using motion processing on mobile devices |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5440326A (en) | 1990-03-21 | 1995-08-08 | Gyration, Inc. | Gyroscopic pointer |
DE112007000074T5 (de) * | 2006-07-10 | 2009-04-02 | Memsic Inc., Andover | System zum Erfassen einer Gierrate unter Verwendung eines Magnetfeldsensors sowie tragbare elektronische Vorrichtungen unter Verwendung hiervon |
WO2008128087A1 (en) * | 2007-04-13 | 2008-10-23 | Keynetik, Inc. | A force sensing apparatus and method to determine the radius of rotation of a moving object |
-
2010
- 2010-01-26 TW TW99102073A patent/TWI407102B/zh active
- 2010-03-29 US US12/748,461 patent/US8797262B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200538751A (en) * | 2004-04-30 | 2005-12-01 | Hillcrest Lab Inc | Free space pointing devices and methods |
TW200639406A (en) * | 2005-05-13 | 2006-11-16 | Ind Tech Res Inst | Initial sensing input apparatus |
US20070214886A1 (en) * | 2006-03-15 | 2007-09-20 | Qualcomm Incorporated | Sensor-based orientation system |
TW200825867A (en) * | 2006-12-12 | 2008-06-16 | Ind Tech Res Inst | Inertial input apparatus with six-axial detection ability and the opearting method thereof |
US20090303204A1 (en) * | 2007-01-05 | 2009-12-10 | Invensense Inc. | Controlling and accessing content using motion processing on mobile devices |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI637344B (zh) * | 2016-02-18 | 2018-10-01 | 緯創資通股份有限公司 | 空間繪畫計分方法、用於進行空間繪畫計分的裝置及系統 |
Also Published As
Publication number | Publication date |
---|---|
US20110181505A1 (en) | 2011-07-28 |
TW201126166A (en) | 2011-08-01 |
US8797262B2 (en) | 2014-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI407102B (zh) | 三維空間運動感測方法 | |
JP4980541B2 (ja) | 座標測定マシン(cmm)の振動に起因した座標測定誤差を補正する方法および装置 | |
CN102289306B (zh) | 姿态感知设备及其定位、鼠标指针的控制方法和装置 | |
Favre et al. | Quaternion-based fusion of gyroscopes and accelerometers to improve 3D angle measurement | |
TWI512548B (zh) | 移動軌跡產生方法 | |
KR100533106B1 (ko) | 지자계 센서의 자세 오차 보상장치 및 방법 | |
CN107655493B (zh) | 一种光纤陀螺sins六位置***级标定方法 | |
US8396684B2 (en) | Method and system for motion tracking | |
CN111678538A (zh) | 一种基于速度匹配的动态水平仪误差补偿方法 | |
US20130151195A1 (en) | System and method for compensating orientation of a portable device | |
Cismas et al. | Crash detection using imu sensors | |
JP6739662B2 (ja) | 空間的な向きの決定 | |
EP2930467A1 (en) | A system and method for sensing the inclination of a moving platform with respect to gravity | |
CN112567097B (zh) | 用于确定机器的作业设备的角度的方法 | |
CN104570142B (zh) | 一种重力梯度仪重力梯度测量信号的解调方法 | |
CN116972875B (zh) | 一种基于陀螺仪的物体运动轨迹监测方法 | |
CN107576334B (zh) | 惯性测量单元的标定方法及装置 | |
CN102141845B (zh) | 三维空间运动感测方法 | |
CN107402004B (zh) | 传感器的姿态信息获取方法及装置 | |
CN112254717B (zh) | 一种基于冷原子干涉陀螺仪的惯性导航装置及方法 | |
CN117782154A (zh) | 一种地面旋转状态下的滚转角测量精度评估方法 | |
CN113959464B (zh) | 一种陀螺仪辅助的加速度计现场校准方法和*** | |
CN108716925A (zh) | 一种九轴传感器的标定方法及装置 | |
Žumer et al. | An advanced nonlinear model of a low-g MEMS accelerometer for a computer pen | |
CN113227714B (zh) | 用于表征惯性测量单元的方法 |