CN105044798A - 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法 - Google Patents

Abstract

本发明公开了一种旋转加速度计重力梯度仪加速度计标度因子反馈调整方法，重力梯度仪旋转圆盘上四只加速度计信号经过加法、减法放大器实现信号的加减操作，加法放大器的两组输出信号和减法放大器的输出信号分别经过三个带通滤波器，对滤波器输出的三组信号分别进行A/D转换，输出的数字信号分别被一倍圆盘旋转频率的正余弦参考信号、圆盘摇摆频率的正弦参考信号解调，解调得到的信号分别经过平滑滤波、PID控制器和D/A转换，输出的模拟信号反馈到相应加速度计的标度因数调整输入端，实现加速度计标度因数实时在线调整。本发明可以进一步提高旋转加速度计梯度仪的加速度计标度因数在线调整的一致性精度，实现高精度的梯度信号解调。

Description

旋转加速度计重力梯度仪加速度计标度因子反馈调整方法

技术领域

本发明涉及一种旋转加速度计重力梯度仪加速度计标度因子调整方法，尤其涉及一种旋转加速度计重力梯度仪加速度计标度因子一致性实时在线反馈调整方法。

背景技术

地球重力场反映了地球表层和地球内部的密度分布和物质运动状态，重力梯度反映的是重力场在空间的变化率，具有比重力更高的空间分辨率，能够更加精确的反映场源体的细节。高精度重力梯度测量对于空间科学、地球科学、地质科学的发展以及在惯性导航等方面具有非常重要的意义，同时重力梯度测量被认为是资源勘探的最有效的手段之一，对于基础地质调查、基础地质研究、大地测量、油气矿藏资源勘查等领域具有重要的应用前景，航空重力梯度测量对山区、无人区、海岛和沿海大陆架部分的基础数据获取具有重要应用价值。

重力梯度信号非常微弱，且受到地形、地质噪声、电气噪声和机械噪声等各种噪声的影响，特别是移动式重力梯度仪面临严峻的线加速度计和角运动加速度等外部环境的影响，使得原本微弱的重力梯度信号更加难以获取，如何消除或降低这些干扰和噪声给科技工作者带来了很大挑战。

旋转加速度计重力梯度仪测量原理如图1所示，四个加速度计相互正交的安装在旋转圆盘上，加速度计敏感轴正切于圆盘边沿，相对的两只加速度计敏感轴方向相反。四只加速度计输出信号经过加法-减法放大电路进行信号组合和信号放大，最后经过滤波/梯度解调进行处理，最后得到此时的重力梯度分量；重力梯度仪在以角频率为ω旋转状态下，四只加速度计输出信号组合基本输出信号形式为：

$\begin{array}{c}(A_1+A_2)-(A_3+A_4)\\ =f\left(c\mathrm{\ω}\right)+f\left(d\mathrm{\ω}\right)+(K_1+K_2-K_3-K_4)R\stackrel{\·}{\mathrm{\ω}}+\\ \[a(K_1-K_2)+b(K_3-K_4)\]\sin \mathrm{\ω}t+\[b(K_1-K_2)+a(K_3-K_4)\]\cos \mathrm{\ω}t+\\ (K_1+K_2+K_3+K_4)R\[\frac{1}{2}({\mathrm{\Γ}}_{xx}-{\mathrm{\Γ}}_{yy})\sin 2\mathrm{\ω}t+(-{\mathrm{\Γ}}_{xy})\cos 2\mathrm{\ω}t\]\end{array}$

其中，K₁、K₂、K₃、K₄分别为加速度计A₁、A₂、A₃、A₄的标度因数，Γ_xx、Γ_yy、Γ_xy为圆盘中心处的三个重力梯度张量分量，R为圆盘中心到加速度计质量中心的距离，为圆盘旋转角加速度，[a(K₁-K₂)+b(K₃-K₄)]和[b(K₁-K₂)+a(K₃-K₄)]分别为圆盘转速正余弦信号幅值，a,b为重力梯度仪与外部环境相关系数，f(cω)为与小于圆盘旋转频率ω的低频信号成分，c<1，f(dω)为与大于圆盘旋转频率ω的高频信号成分，d>1，由于转台不平稳产生的旋转角加速度会带来的二倍频信号干扰，此干扰信号频率与梯度信号频率一致，会严重影响梯度信号的精确性，只有通过调整加速度计标度因数，使得K₁＝K₂＝K₃＝K₄，消除圆盘旋转角加速度所引入的二倍频干扰信号，最终解调出高精度的重力梯度分量(Γ_xx–Γ_yy)和(–Γ_xy)。

发明内容

技术问题：本发明的目的在于提高旋转加速度计重力梯度仪圆盘上四只加速度计标度因子的一致性，提出一种设计合理且能够实现加速度计标度因子实时在线反馈方法，可以快速实现四只加速度计标度因子的一致性，有利于提高重力梯度解调精度的旋转加速度计重力梯度仪加速度计标度因子反馈调整方法。

技术方案：本发明的旋转加速度计重力梯度仪加速度计标度因子反馈调整方法，包括以下步骤：

1)将四只高精度加速度计均匀正交放置在旋转圆盘上，相对两只加速度计的输入轴方向相反，相邻两只加速度计的输入轴相互垂直，且四只加速度计的输入轴方向正切于圆盘边沿，其中加速度计A₁和A₂为一对，A₃和A₄为一对；将圆盘坐标系设为东北天地理坐标系，记为OXYZ，原点为O，圆盘圆心位于坐标系原点，圆盘旋转轴垂直于圆盘平面，旋转加速度计重力梯度仪的圆盘旋转角频率为ω；

2)将加速度计A₁、A₂的输出信号和加速度计A₃、A₄的输出信号分别经过加法放大器实现信号相加放大操作，将加法放大器的输出信号V₁₂和V₃₄经过减法放大器实现信号相减放大操作，得到输出信号E，所述加法放大器的输出信号V₁₂与加速度计A₁、A₂对应，输出信号V₃₄与加速度计A₃、A₄对应；

3)将所述步骤2)中得到的信号V₁₂和V₃₄分别经过第一带通滤波器BPF1和第二带通滤波器BPF2实现信号滤波处理，然后将两个带通滤波器的输出信号分别转换为数字信号，以参考信号sinωt为解调基准信号，将第一带通滤波器BPF1的数字信号解调得到幅值1，以参考信号cosωt为解调基准信号，将第二带通滤波器BPF2的数字信号解调得到幅值2，其中t为时间；

同时将所述步骤2)中得到的信号E经过第三带通滤波器BPF3滤波处理后，将BPF3的输出信号转换为数字信号，以参考信号sinω_st为解调基准信号，对所述的数字信号进行幅值解调得到幅值3，其中ω_s为旋转加速度计重力梯度仪的圆盘摇摆角频率；

4)将所述步骤3)中得到的三组解调后的幅值信号，分别用平滑滤波器进行平滑处理；

5)将所述步骤4)中得到的三组平滑处理后的信号分别经过比例、微分积分控制器PID1、PID2和PID3进行调节控制，得到三组调整加速度计标度因子的数字信号；

6)将所述的三组调整加速度计标度因子的数字信号分别经过数模转换，得到与所述三组幅值信号对应的三组模拟电流信号；

7)将所述步骤6)中得到的与幅值1对应的模拟电流信号反馈到第一加速度计A₁的标度因数调整输入端口；将与幅值2对应的模拟电流信号反馈到第三加速度计A₃的标度因数调整输入端口；将与幅值3对应的模拟电流信号反馈到位于第四加速度计A₄的标度因数调整输入端口，所述的四个加速度计相互正交安装在旋转圆盘上，第一加速度计A₁和第二加速度计A₂关于圆盘中心对称设置，第三加速度计A₃和第四加速度计A₄关于圆盘中心对称设置，加速度计敏感轴正切于圆盘边沿，相对的两只加速度计敏感轴方向相反。

进一步的，所述步骤3)中的第一带通滤波器BPF1、第二带通滤波器BPF2和第三带通滤波器BPF3均为四阶巴特沃斯模拟滤波器，其中带通滤波器BPF1、BPF2的中心频率为ω/(2π)，过渡带带宽为2ω/(5π)，通带带宽为ω/(5π)；带通滤波器BPF3的中心频率为ω_s/(2π)，过渡带带宽为2ω_s/(5π)，通带带宽为ω_s/(5π)。

进一步的，所述步骤3)中的参考信号sinωt、cosωt和sinω_st均由数字信号处理器芯片产生。

本发明采用三组带通模拟滤波器分别对两组加法器输出信号和一组减法器输出信号进行滤波处理，然后对分别对滤波后的信号进行A/D转换、平滑滤波、PID调整和DAC模拟反馈信号输出，实现加速度计标度因数实时在线反馈调整，最终达到四只加速度计标度因数一致。

有益效果：本发明与现有技术相比，具有以下优点：

本发明是采用三组模拟带通滤波器分别对两组加法器输出信号和一组减法器输出信号进行滤波处理，采用带通滤波器可以有效过滤掉低频漂移和高频干扰，比传统的采用低通滤波器效果更优。带通滤波器输出后的梯度信号经过A/D转换、幅值解调、数据平滑、PID调整和D/A转换电路处理后，DAC输出的模拟电流信号反馈到相应的加速度计标度因数调整输入端，实现旋转加速度计重力梯度仪上加速度计标度因数在线实时调整，最后使得四只加速度计标度因数趋于一致，为高精度梯度信号解调提供保障。

附图说明

图1为旋转加速度计重力梯度仪原理图。

图2为加速度计标度因数在线反馈调整示意图。

图3为带通滤波器电路原理图。

具体实施方式

下面结合实施例和说明书附图对本发明作进一步地说明。

本发明的旋转加速度计重力梯度仪加速度计标度因子调整方法具体流程步骤如下：

1)将四只高精度加速度计均匀正交放置在旋转圆盘上，相对两只加速度计的输入轴方向相反，相邻两只加速度计的输入轴相互垂直，且四只加速度计的输入轴方向正切于圆盘边沿，其中加速度计A₁和A₂为一对，A₃和A₄为一对；将圆盘坐标系设为东北天地理坐标系，记为OXYZ，原点为O，圆盘圆心位于坐标系原点，圆盘旋转轴垂直于圆盘平面，旋转加速度计重力梯度仪的圆盘旋转角频率为ω；

2)将加速度计A₁、A₂的输出信号和加速度计A₃、A₄的输出信号分别经过加法放大器实现信号相加放大操作，将加法放大器的输出信号V₁₂和V₃₄经过减法放大器实现信号相减放大操作，得到输出信号E，所述加法放大器的输出信号V₁₂与加速度计A₁、A₂对应，输出信号V₃₄为与加速度计A₃、A₄对应；

3)将所述步骤2)中得到的信号V₁₂和V₃₄分别经过第一带通滤波器BPF1和第二带通滤波器BPF2实现信号滤波处理，然后将两个带通滤波器的输出信号分别转换为数字信号，以参考信号sinωt为解调基准信号，将第一带通滤波器BPF1的数字信号解调得到幅值1，以参考信号cosωt为解调基准信号，将第二带通滤波器BPF2的数字信号解调得到幅值2；

同时将所述步骤2)中得到的信号E经过第三带通滤波器BPF3滤波处理，然后将BPF3的输出信号转换为数字信号，以参考信号sinω_st为解调基准信号，对所述的数字信号进行幅值解调得到幅值3，其中ω_s为旋转加速度计重力梯度仪的圆盘摇摆角频率；

4)将所述步骤3)中得到的三组解调后的幅值信号，分别用平滑滤波器进行平滑处理；

5)将所述步骤4)中得到的三组平滑处理后的信号分别经过比例、微分积分控制器PID1、PID2和PID3进行调节控制，得到三组调整加速度计标度因子的数字信号；

6)将所述的三组调整加速度计标度因子的数字信号分别经过数模转换，得到与所述三组幅值信号对应的三组模拟电流信号；

7)将所述步骤6)中得到的与幅值1对应的模拟电流信号反馈到第一加速度计A₁的标度因数调整输入端口；将与幅值2对应的模拟电流信号反馈到第三加速度计A₃的标度因数调整输入端口；将与幅值3对应的模拟电流信号反馈到位于第四加速度计A₄的标度因数调整输入端口，所述的四个加速度计相互正交安装在旋转圆盘上，第一加速度计A₁和第二加速度计A₂关于圆盘中心对称设置，第三加速度计A₃和第四加速度计A₄关于圆盘中心对称设置，加速度计敏感轴正切于圆盘边沿，相对的两只加速度计敏感轴方向相反。

实施例：

如图2所示，具体实施方式如下：

所述加速度计标度因子反馈调整方法由三个反馈回路组成。第一反馈回路为：加速度计A₁和加速度计A₂输出信号经过加法放大器实现两信号相加操作，加法器1的输出信号V₁₂依次经过BPF1、ADC、1倍频幅值解调、平滑滤波、PID1调整和DAC，DAC输出的模拟电流信号反馈给加速度计A1的标定因子调整端口；实现方法为：

相对加速度计A₁、A₂经过加法器1相加后的输出信号形式为：

$\begin{array}{c}{V}_{12}=K(A_1+A_2)\\ =f\left(c_1\mathrm{\&omega;}\right)+f\left(d_1\mathrm{\&omega;}\right)+K(K_1-K_2)(a\sin \mathrm{\&omega;}t+b\cos \mathrm{\&omega;}t)+\\ K(K_1+K_2)R\&lsqb;\frac{1}{2}({\mathrm{\&Gamma;}}_{xx}-{\mathrm{\&Gamma;}}_{yy})\sin 2\mathrm{\&omega;}t+(-{\mathrm{\&Gamma;}}_{xy})\cos 2\mathrm{\&omega;}t+\stackrel{\&CenterDot;}{\mathrm{\&omega;}}\&rsqb;\end{array}---\left(1\right)$

其中，f(c₁ω)为与小于圆盘旋转频率ω的低频信号成分，c₁<1，f(d₁ω)为与大于圆盘旋转频率ω的高频信号成分，d₁>1，a，b为一倍频信号幅值，K为信号放大倍数。输出信号V₁₂经过带通滤波器，把低于圆盘旋转频率的低频信号和高于圆盘旋转频率的高频信号成分过滤掉，输出与圆盘旋转频率相同的信号，即：

$V_{12}^{BPF}=K(K_1-K_2)(asin\mathrm{\&omega;}t+bcos\mathrm{\&omega;}t)---\left(2\right)$

然后带通滤波器输出信号经过A/D转换，转换成数字信号D₁₂，对采样到的数字信号D₁₂进行幅值解调，解调公式为：

$D_{12}^1=\frac{2\underset{i=1}{\overset{N}{\&Sigma;}}{D}_{12}\left(i\right)\sin \mathrm{\&omega;}\mathrm{\&Delta;}t(i-1)}{NK}=(K_1-K_2)a---\left(3\right)$

解调得到的幅值信号经过平滑、PID控制器调整，变为较为稳定、精确的数值，然后经过D/A转换器，转换为模拟电流信号A₁₂，模拟信号A₁₂反馈给加速度计A₁的标度因子调整输入端，使加速度计A₁的标度因子跟踪加速度计A₂的标度因子，最终实现加速度计A₁与加速度计A₂的标度因子一致，即：K₁＝K₂；

第二反馈回路为：加速度计A₃和加速度计A₄输出信号经过加法放大器2实现两信号相加操作，加法器2的输出信号V₃₄依次经过BPF2、ADC、1倍频幅值解调、平滑滤波、PID2调整和DAC，DAC输出的模拟电流信号反馈给加速度计A₃的标定因子调整端口；实现方法为：

相对加速度计A₃、A₄经过加法器2相加后的输出信号形式为：

$\begin{array}{c}{V}_{34}=K(A_3+A_4)\\ =f\left(c_2\mathrm{\&omega;}\right)+f\left(d_2\mathrm{\&omega;}\right)+K(K_3-K_4)(a\cos \mathrm{\&omega;}t+b\sin \mathrm{\&omega;}t)-\\ K(K_3+K_4)R\&lsqb;\frac{1}{2}({\mathrm{\&Gamma;}}_{xx}-{\mathrm{\&Gamma;}}_{yy})\sin 2\mathrm{\&omega;}t+(-{\mathrm{\&Gamma;}}_{xy})\cos 2\mathrm{\&omega;}t-\stackrel{\&CenterDot;}{\mathrm{\&omega;}}\&rsqb;\end{array}---\left(4\right)$

其中，其中，f(c₂ω)为与小于圆盘旋转频率ω的低频信号成分，c₂<1，f(d₂ω)为与大于圆盘旋转频率ω的高频信号成分，d₂>1，对加法器A₃、A₄输出信号V₃₄进行带通滤波处理，把低于圆盘旋转频率的低频信号和高于圆盘旋转频率的高频信号成分过滤掉，输出与圆盘旋转频率相同的信号，即：

$V_{34}^{BPF}=K(K_3-K_4)(acos\mathrm{\&omega;}t+bsin\mathrm{\&omega;}t)---\left(5\right)$

然后对滤波器输出信号进行A/D转换，转换成数字信号D₃₄，对采样到的数字信号D₃₄进行幅值解调，解调公式为：

$D_{34}^2=\frac{2\underset{i=1}{\overset{N}{\&Sigma;}}{D}_{34}\left(i\right)cos\mathrm{\&omega;}\mathrm{\&Delta;}t(i-1)}{NK}=(K_3-K_4)a---\left(6\right)$

解调得到的幅值信号经过平滑、PID控制器调整，变为较为稳定、精确的数值，然后再经过D/A转换器，转换为模拟电流信号A₃₄，模拟信号A₃₄反馈给加速度计A₃的标度因子调整输入端，使加速度计A₃的标度因子跟踪加速度计A₄的标度因子，最终实现加速度计A₃与加速度计A₄的标度因子一致，即：K₃＝K₄；

第三反馈回路为：加法信号V₁₂和V₃₄经过减法放大器实现两信号相减操作，减法器的输出信号E依次经过BPF3、ADC、圆盘摇摆信号幅值解调、平滑滤波、PID3调整和DAC，DAC输出的模拟电流信号反馈给加速度计A₄的标定因子调整端口；

四只加速度计经过加法器、减法器后的输出信号形式为：

$\begin{array}{c}E=K\&lsqb;(A_1+A_2)-(A_3+A_4)\&rsqb;\\ =f\left(c\mathrm{\&omega;}\right)+f\left(\mathrm{\&omega;}\right)+f\left(d\mathrm{\&omega;}\right)+K(K_1+K_2-K_3-K_4)R\stackrel{\&CenterDot;}{\mathrm{\&omega;}}+\\ K(K_1+K_2+K_3+K_4)R\&lsqb;\frac{1}{2}({\mathrm{\&Gamma;}}_{xx}-{\mathrm{\&Gamma;}}_{yy})\sin 2\mathrm{\&omega;}t+(-{\mathrm{\&Gamma;}}_{xy})\cos 2\mathrm{\&omega;}t\&rsqb;\end{array}---\left(7\right)$

其中，f(ω)为与圆盘旋转频率ω相同的信号成分，从上式表达式可以看出，减法器输出信号中含有圆盘旋转角加速度信号但角加速度信号是一种不平稳信号，里面含有各种频率成分信息，若含有两倍圆盘旋转频率信息且幅值较大，会严重干扰梯度信息的精确解调，必须消除此影响。而角加速度信号幅度是由四只加速度计标度因数组合而成，只要实现K₄＝K₂，再结合第一、第二反馈回路标度因数调整方法，就可以实现四只加速度计标度因数一致，实现方法是给圆盘施加小幅度高频摇摆信号，设定圆盘摇摆频率为ω_s，摇摆振幅幅值为δ，即：δsinω_st，此时圆盘旋转角加速度变为：带入减法器输出信号表达式(7)，变为：

$\begin{array}{c}E=K\&lsqb;(A_1+A_2)-(A_3+A_4)\&rsqb;\\ =f\left(c\mathrm{\&omega;}\right)+f\left(\mathrm{\&omega;}\right)+f\left(d\mathrm{\&omega;}\right)+K(K_1+K_2-K_3-K_4)R(2\stackrel{\&CenterDot;}{\mathrm{\&omega;}}-{\mathrm{\&delta;\&omega;}}_s^2{\mathrm{sin\&omega;}}_{s}t)+\\ K(K_1+K_2+K_3+K_4)R\&lsqb;\frac{1}{2}({\mathrm{\&Gamma;}}_{xx}-{\mathrm{\&Gamma;}}_{yy})\sin 2\mathrm{\&omega;}t+(-{\mathrm{\&Gamma;}}_{xy})\cos 2\mathrm{\&omega;}t\&rsqb;\end{array}---\left(8\right)$

减法器输出信号经过带通滤波器，把低于摇摆频率的低频信号和高于摇摆频率的高频信号成分滤掉，剩下与摇摆频率相同的频率成分，即：

$E^{BPF}=K(K_1+K_2-K_3-K_4){\mathrm{R\&delta;\&omega;}}_s^2{\mathrm{sin\&omega;}}_{s}t---\left(9\right)$

然后对滤波后的信号进行A/D转换，得到数字信号D，然后对其进行幅值解调，解调公式为：

$D_{1234}^3=\frac{2\underset{i=1}{\overset{N}{\&Sigma;}}D\left(i\right){\mathrm{sin\&omega;}}_s\mathrm{\&Delta;}t(i-1)}{{\mathrm{NKR\&delta;\&omega;}}_s^2}=(K_1+K_2-K_3-K_4)---\left(10\right)$

解调得到的幅值信号经过平滑、PID控制器调整，变为较为稳定、精确的数值，然后再经过D/A转换器，转换为模拟电流信号A₁₂₃₄，模拟信号A₁₂₃₄反馈给加速度计A₄的标度因子调整输入端，使加速度计A₄的标度因子跟踪加速度计A₂的标度因子，最终实现加速度计A₄与加速度计A₂的标度因子一致，即：K₄＝K₂；从第一、二反馈调整回路中实现了K₁＝K₂，K₃＝K₄，而第三反馈回路实现K₄＝K₂，最终实现四只加速度计标度因子一致。

所述步骤3)中的参考信号sinωt、cosωt均由DSP芯片产生，信号幅值解调参考信号sinωt、cosωt和sinω_st由DSP芯片发出，但该数字解调功能并不局限于DSP芯片，利用其他数字处理器也可以实现数字信号幅值解调功能。

所述步骤3)中的第一带通滤波器BPF1、第二带通滤波器BPF2和第三带通滤波器BPF3均为四阶巴特沃斯模拟滤波器，其中带通滤波器BPF1、BPF2的中心频率为ω/(2π)，过渡带带宽为2ω/(5π)，通带带宽为ω/(5π)；带通滤波器BPF3的中心频率为ω_s/(2π)，过渡带带宽为2ω_s/(5π)，通带带宽为ω_s/(5π)。带通滤波电路由两个两阶巴特沃斯带通滤波器级联组成一个四阶巴特沃斯带通滤波器，如图3所示；第一级二阶带通滤波电路由运算放大器A1，电阻R1、R2、R3，电容C1、C2构成；第二级二阶带通滤波电路由运算放大器A2，电阻R4、R5、R6，电容C3、C4构成。

利用本发明的加速度计标度因数在线反馈调整方法，可以精确实现实时在线高稳定性、高精度的加速度计标度因数调整，实现四只加速度计标度因数高度一致性，实现高精度的梯度解调。

应理解上述实施例仅用于说明本发明技术方案的具体实施方式，而不用于限制本发明的范围。在阅读了本发明之后，本领域技术人员对本发明的各种等同形式的修改和替换均落于本申请权利要求所限定的保护范围。

Claims (3)

1.一种旋转加速度计重力梯度仪加速度计标度因子反馈调整方法，其特征在于，该方法包括以下步骤：

1）将四只高精度加速度计均匀正交放置在旋转圆盘上，相对两只加速度计的输入轴方向相反，相邻两只加速度计的输入轴相互垂直，且四只加速度计的输入轴方向正切于圆盘边沿，其中加速度计A₁和A₂为一对，A₃和A₄为一对；将圆盘坐标系设为东北天地理坐标系，记为OXYZ，原点为O，圆盘圆心位于坐标系原点，圆盘旋转轴垂直于圆盘平面，旋转加速度计重力梯度仪的圆盘旋转角频率为ω；

2）将加速度计A₁、A₂的输出信号和加速度计A₃、A₄的输出信号分别经过加法放大器实现信号相加放大操作，将加法放大器的输出信号V ₁₂和V ₃₄经过减法放大器实现信号相减放大操作，得到输出信号E，所述加法放大器的输出信号V ₁₂与加速度计A₁、A₂对应，输出信号V ₃₄为与加速度计A₃、A₄对应；

3）将所述步骤2）中得到的信号V ₁₂和V ₃₄分别经过第一带通滤波器BPF1和第二带通滤波器BPF2实现信号滤波处理，然后将两个带通滤波器的输出信号分别转换为数字信号，以参考信号sinωt为解调基准信号，将第一带通滤波器BPF1的数字信号解调得到幅值1，以参考信号cosωt为解调基准信号，将第二带通滤波器BPF2的数字信号解调得到幅值2，其中t为时间；

同时将所述步骤2）中得到的信号E经过第三带通滤波器BPF3滤波处理，然后将BPF3的输出信号转换为数字信号，以参考信号sinω _s t为解调基准信号，对所述的数字信号进行幅值解调得到幅值3，其中ω _s 为旋转加速度计重力梯度仪的圆盘摇摆角频率；

4）将所述步骤3）中得到的三组解调后的幅值信号，分别用平滑滤波器进行平滑处理；

5）将所述步骤4）中得到的三组平滑处理后的信号分别经过比例、微分积分控制器PID1、PID2和PID3进行调节控制，得到三组调整加速度计标度因子的数字信号；

6）将所述的三组调整加速度计标度因子的数字信号分别经过数模转换，得到与所述三组幅值信号对应的三组模拟电流信号；

7）将所述步骤6）中得到的与幅值1对应的模拟电流信号反馈到第一加速度计A₁的标度因数调整输入端口；将与幅值2对应的模拟电流信号反馈到第三加速度计A₃的标度因数调整输入端口；将与幅值3对应的模拟电流信号反馈到位于第四加速度计A₄的标度因数调整输入端口，所述的四个加速度计相互正交安装在旋转圆盘上，第一加速度计A₁和第二加速度计A₂关于圆盘中心对称设置，第三加速度计A₃和第四加速度计A₄关于圆盘中心对称设置，加速度计敏感轴正切于圆盘边沿，相对的两只加速度计敏感轴方向相反。

2.根据权利要求1所述的旋转加速度计重力梯度仪加速度计标度因子反馈调整方法，其特征在于，所述步骤3）中的第一带通滤波器BPF1、第二带通滤波器BPF2和第三带通滤波器BPF3均为四阶巴特沃斯模拟滤波器，其中带通滤波器BPF1、BPF2的中心频率为ω/(2π)，过渡带带宽为2ω/(5π)，通带带宽为ω/(5π)；带通滤波器BPF3的中心频率为ω _s /(2π)，过渡带带宽为2ω _s /(5π)，通带带宽为ω _s /(5π)。

3.根据权利要求1所述的旋转加速度计重力梯度仪加速度计标度因子反馈调整方法，其特征在于，所述步骤3）中的参考信号sinωt、cosωt和sinω _s t均由数字信号处理器芯片产生。