[發明專利]一種全自動重力測試儀高精度動態數據處理方法

有效

• 申請號: 201811553674.9
• 申請日: 2018-12-19
• 公開（公告）號: CN109471192B
• 公開（公告）日: 2021-09-14
• 發明（設計）人: 金瑩;陳家俊
• 申請（專利權）人: 湖南航天機電設備與特種材料研究所
• 主分類號: G01V7/06
• 分類號: G01V7/06
• 代理公司: 長沙正奇專利事務所有限責任公司 43113
• 代理人: 郭立中;李美麗

1.一種全自動重力測試儀高精度動態數據處理方法，其特征在于，包括以下步驟：

步驟A.獲得重力測量儀的輸出信號y(k)＝s(k)+v(k)，其中，s(k)為重力異常特征信號，k為采樣時刻且k＝1,2,…,N，v(k)為k時刻測量噪聲；N為采樣點的數目；

步驟B.設根據時間序列分析方法建立的重力測試儀隨機誤差p階AR模型為

其中，{e(k)}為干擾噪聲序列，v(k)⊥e(k)；

將系統方程表示為：

y(k)＝H(k)·x(k)+v(k)，

其中，x(k+1)為k+1時刻系統待估計的狀態變量；y(k)為k時刻重力測試儀的觀測值；為k時刻到k+1時刻系統狀態一步轉移矩陣；w(k)為k時刻系統噪聲向量；Γ(k+1,k)為系統噪聲矩陣；H(k)為k時刻系統量測轉移矩陣；

步驟C.建立狀態一步預測方程：

一步預測均方誤差：

誤差序列：v(k)＝y(k)-H(k)·x(k/k-1)，

估計量測噪聲：

R(k)＝R(k-1)+(v(k)·v^T(k)-H(k)·P(k/k-1)·H^T(k))，

濾波增益：

K(k)＝P(k/k-1)·H^T(k)·(H(k)·P(k/k-1)·H^T(k)+R(k))^-1，

狀態估計：

估計均方誤差：P(k)＝(I-K(k)·H(k))·P(k/k-1)。