一种用于外差干涉测振系统的信号解调装置和方法。该信号解调方法包括：产生载波调制模拟信号s(t)和AOM驱动信号fAOM；将模拟信号s(t)转换为数字信号s(n)；产生单位正交信号I0和Q0，并与s(n)混频，得到正交信号IS和QS；对混频后的正交信号IS和Qs信号进行抽取和滤波；对IS和QS整形和进一步滤波，只保留含有基带信号的正交信号IB和QB；根据坐标旋转原理进行相位计算得到相位对超出[‑π，π]量程的相位进行补偿，得到真实相位将数字相位信号转换为模拟位移信号d(t)输出。本发明在FPGA中利用数字方法实现了核心模块的设计，提高了系统的可移植性和修改性，算法抗干扰性强，准确度高。
1.一种用于外差干涉测振系统的信号解调装置，其特征在于，包括：外差干涉光路，用于产生经过振动信号调制的载波调制信号 ，并产生AOM驱动参考信号 ；
信号采集模块，用于将所述载波调制模拟信号 转换为载波调制数字信号 ，并将所述AOM驱动参考信号 转换为数字信号；
本振信号模块，用于基于数字 信号产生两路单位正交信号 和 ，其中所述两路单位正交信号 和 的信号频率与载波调制数字信号 的信号频率相同，数据速率与信号采集模块的采样率 相同；
混频模块，用于将载波调制数字信号 与所述两路单位正交信号 和 进行混频，得到两路正交信号 和 ，所述正交信号 和 同时保留了倍频项和基带成分；
数据抽取模块，用于对所述正交信号 和 降低数据速率和滤波，滤除其部分倍频成分，保留基带成分；
滤波模块，用于对所述正交信号 和 进一步滤波，经过滤波模块处理后，得到信号 和，所述正交信号只保留了基带信号；
Cordic相位计算模块，用于对滤波模块处理后的所述正交信号 和 的相位进行计算；
相位补偿模块，针对目标相位信号 超出上述Cordic相位计算模块量程 的情况下，对Cordic相位计算结果 进行补偿，得到具有真实相位的相位信号 ；
数据输出模块，用于将相位补偿模块输出的数字相位信号转换为模拟位移信号后输出。
2.如权利要 求1所述的 装置 ，其特 征在于 ，所述载波 调制模拟 信号，其中A为载波幅度，为载波频率， 为振动信号调制相位，所述载波调制模拟信号 中包含有振动调制信息 。
3.如权利要求1所述的装置，其特征在于，信号采集模块将载波调制模拟信号 进行离散化处理，离散化后的载波调制数字信号 的表达式为 ，为数字化后的载波幅度，为数字化后的载波频率， 为振动信号调制相位。
4.如权利要求3所述的装置，其特征在于，所述数据采样模块的采样率 满足，为模拟载波频率，则离散化后的载波频率为 。
5.如权利要求3所述的装置，其特征在于，离散化的信号 在混频模块与本振模块产生的单位正交信号 及 混频，混频公式如下：路：
混频后，得到正交化的 、 两路信号，包括倍频项与基带项。
6.如权利要求5所述的装置，其特征在于，数据抽取模块滤波之后，所述正交信号 和消除了倍频项，得到正交信号 和 ，和 只保留基带项，表达式为：路： ，
7.如权利要求3所述的装置，其特征在于，Cordic相位计算模块依据坐标旋转原理得到目标相位信号 ， 与反正切算法 得到的值相同，Cordic算法相位 计算范围为 ，采用Cordic相位计算模块避免反正切算法中分母 为零时无法计算的情况。
8.如权利要求7所述的装置，其特征在于，相位补偿模块基于Cordic相位计算模块对产生 之外的 信号进行补偿，根据产生折叠时相位值 与真实相位值 相差 的关系，对相位 进行补偿，得到补偿后的相位信号 ，其中k为整数，，然后数据输出模块将相位补偿模块输出的数字相位信号转换为模拟信号，并且滤除噪声，得到准确的目标信号，通过振动位移与目标相位的关系， ，得到物体振动的位移d，其中 为激光波长。
9.一种用于外差干涉测振系统的信号解调方法，其特征在于，包括以下步骤：步骤1，将外差干涉光路中的光信号转换成电信号，产生载波调制模拟信号 ，所述载波调制模拟信号 中含有目标相位信息；产生AOM驱动信号 ；
步骤2，将所述载波调制信号 转换为数字信号 ；将所述AOM驱动信号 离散化成数字信号；
步骤3，本振信号模块基于所述AOM驱动信号 产生单位正交信号 和 ，然后单位正交信号 和 与载波调制信号 进行混频，得到两路正交信号 和 ，所述正交信号 和中同时保留倍频与基带成分；
步骤4，数据抽取模块对混频后的两路正交信号IS和QS信号进行抽取和滤波，滤除混频产生的部分倍频成分，保留其基带信号；
步骤5，滤波模块对上述得到的正交信号 和 整形和进一步滤波，得到正交信号 和，和 只保留基带信号；
步骤6，Coridic相位计算模块对上述得到的正交信号 和 根据坐标旋转原理进行相位计算，得到原始相位 ；
步骤7，相位补偿模块对实际超出Cordic相位计算模块量程的信号 进行补偿,得到真实相位 ；
步骤8，将上述得到的相位信号 转换为模拟位移信号 输出。
10.如权利要求9所述的方法，其特征在于，步骤（1）中所述载波调制模拟信号，为数字化后的载波幅度，为数字化后的载波频率， 为振动信号调制相位。
11.如权利要求9所述的方法，其特征在于，步骤（2）中信号采集模块对载波调制模拟信号 进行离散化处理，离散化后得到的载波调制数字信号 ，为数字化后的载波幅度，为数字化后的载波频率， 为振动信号调制相位。
12.如权利要求11所述的方法，其特征在于，离散化的信号 在混频模块与本振模块产生的单位正交信号 及 混频，混频公式如下：路：
混频后，得到正交化的 、 两路信号，包括倍频项与基带项。
13.如权利要求12所述的方法，其特征在于，数据抽取模块之后，所述正交信号 和 消除了部分倍频项；再经过滤波模块，得到正交信号 和 ，和 只保留基带项，表达式为：路： ，
14.如权利要求11所述的方法，其特征在于，采用Cordic相位计算模块计算 和 两路信号的相位 ， 与反正切算法 得到的值相同，Cordic算法相位 计算范围为 ，采用Cordic相位计算模块避免反正切算法中分母 为零时无法计算的情况。
15.如权利要求12所述的方法，其特征在于，相位补偿模块基于Cordic相位计算模块对产生 之外的 信号进行补偿，根据产生折叠时相位值 与真实相位值 相差 的关系，对相位 进行补偿，得到补偿后的相位信号 ，其中k为整数，，然后数据输出模块将相位补偿模块输出的数字相位信号 转换为模拟位移信号 ，并且滤除噪声，得到准确的目标信号；通过振动位移与目标的关系， ，得到物体振动的位移d，其中 为激光波长。
用于外差干涉测振系统的信号解调装置和方法
技术领域
[0001] 本发明涉及振动检测领域，尤其涉及一种用于外差干涉光路的振动测量系统(测振系统)的信号解调装置和方法。
背景技术
[0002] 激光多普勒测振技术具有非接触、远距离、灵敏度高、精度高的特点，已经广泛应用于民用、工业生产、军事及航空等领域。
[0003] 如图1所示，激光多普勒测振系统包括干涉光路和电路信号处理两部分，由于外差干涉光路抗干扰性强与精度高的特点，目前测振系统中干涉光路部分通常采用外差光路法，外差干涉光路典型特点就是光路中有频移器(EOM或者AOM)使信号光或参考光产生一定的频移，频移量就是干涉后载波的频率。
[0004] 激光多普勒测振系统的核心就是光学干涉系统和信号解调单元，系统的高精度主要依赖于这两部分的性能，外差光路得到干涉光信号后，系统进行光电信号转换、信号调理和信号解调，信号解调的精度直接决定了系统精度。
[0005] 常见的基于外差干涉光路的振动测量系统中，所采用的频移器件(AOM或者EOM)，最小的频移量为40MHz，如此高的载波如果按照香侬采样定理的要求进行采样，则采样时钟速率要达到几百MHz甚至GHz，由此，这样采用高频载波的好处就是系统抗干扰性强，但缺点是载波越高要求后端采样率也越高，对信号采集模块和信号处理系统的设计带来一定难度，使信号采集电路的设计变得困难，后续信号的处理复杂度和难度均大幅增加。
发明内容
[0006] 有鉴于此，本发明的主要目的在于提供一种用于外差干涉测振系统的信号解调方法，以期至少部分地解决上述技术问题。
[0007] 为了实现上述目的，作为本发明的一个方面，提供了一种用于外差干涉测振系统的信号解调装置，其特征在于，包括：
[0008] 外差干涉光路，用于产生经过振动信号调制的载波调制信号s(t)，并产生AOM驱动参考信号fAOM；
[0009] 信号采集模块，用于将所述载波调制模拟信号s(t)转换为载波调制数字信号s(n)，并将所述AOM驱动参考信号fAOM转换为数字信号；
[0010] 本振信号模块，用于基于数字fAOM信号产生两路单位正交信号I0和Q0，其中所述两路单位正交信号I0和Q0的信号频率与载波调制数字信号s(n)的信号频率相同，数据速率与信号采集模块的采样率相同；
[0011] 混频模块，用于将载波调制数字信号s(n)与所述两路单位正交信号I0和Q0进行混频，得到两路正交信号IS和QS，所述正交信号IS和QS同时保留了倍频和基带成分；
[0012] 数据抽取模块，用于对所述正交信号IS和QS降低数据速率和滤波，滤除其部分倍频成分，保留基带成分；
[0013] 滤波模块，用于对所述正交信号IS和QS进一步滤波，经过滤波模块处理后，消除倍频成分，得到正交信号IB和QB，所述正交信号只保留了基带信号；
[0014] Cordic相位计算模块，用于对滤波模块处理后的所述正交信号IB和QB的相位进行计算；
[0015] 相位补偿模块，针对目标相位信号 超出上述Cordic相位计算模块的量程[-π，π]情况下，对计算结果 进行补偿，得到具有真实相位的相位信号
[0016] 数据输出模块，用于将相位补偿模块输出的数字相位信号转换为模拟位移信号后输出。
[0017] 作为本发明的另一个方面，还提供了一种用于外差干涉测振系统的信号解调方法，其特征在于，包括以下步骤：
[0018] 步骤1，将外差干涉光路中的光信号转换成电信号，产生载波调制模拟信号s(t)，所述载波调制模拟信号s(t)中含有目标相位信息；产生AOM驱动信号fAOM；
[0019] 步骤2，将所述载波调制信号s(t)转换为数字信号s(n)；将所述AOM驱动信号fAOM离散化成数字信号；
[0020] 步骤3，本振信号模块基于所述AOM驱动信号fAOM产生单位正交信号I0和Q0，然后单位正交信号I0和Q0与载波调制信号s(n)进行混频，得到两路正交信号IS和QS，所述正交信号IS和QS中同时保留倍频与基带成分；
[0021] 步骤4，数据抽取模块对混频后的两路正交信号IS和QS信号进行抽取和滤波，滤除混频产生的部分高频成分，其保留基带信号；
[0022] 步骤5，滤波模块对上述得到的正交信号Is和QS整形和进一步滤波，得到正交信号IB和QB，IB和QB只保留基带信号；
[0023] 步骤6，Coridic相位计算模块对上述得到的正交信号IB和QB根据坐标旋转原理进行相位计算，得到原始相位
[0024] 步骤7，相位补偿模块对超出Cordic相位计算模块量程的相位计算结果 进行补偿，得到真实相位
[0025] 步骤8，将上述得到的相位信号 转换为模拟位移信号d(t)输出。
[0026] 基于上述技术方案可知，本发明的信号解调方法相对于现有技术具有如下优点：
[0027] (1)针对外差干涉测振系统中40MHz载波信号的振动调制信号，利用欠采样采样技术，用低于Nyquist频率以下的采样率进行采样，数据率的降低能够有效的降低信号后端处理复杂度；
[0028] (2)在FPGA中实现数字混频、数字抽取、数字滤波、相位计算功能，实现了目标信号的解调，得到了物体振动信息，通过在FPGA中利用数字方法实现了核心模块的设计，提高了系统的可移植性和修改性；
[0029] (3)本发明利用模/数转换和数/模转换，将信号变化到数字域处理，最后将目标信号转换成模拟信号输出，这样做的优点是算法的抗干扰性强，准确度高，保证了外差干涉测振系统的输出精度。
附图说明
[0030] 图1为现有技术的外差干涉测振系统的组成方框示意图；
[0031] 图2为本发明的外差干涉测振系统中振动信号解调的流程框架结构图；
[0032] 图3为本发明中信号变换的流程图；
[0033] 图4为本发明的外差干涉测振系统的信号解调方法的流程图。
具体实施方式
[0034] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。
[0035] 本发明公开了一种带通采样方法来实现振动信号的解调，在FPGA中实现，这种解调方法将采样率降低到Nyquist频率以下，该方案具有采样率低的特点。
[0036] 具体地，本发明公开了一种用于外差干涉测振系统的信号解调装置，包括：
[0037] 外差干涉光路，用于产生经过振动信号调制的载波调制信号s(t)，并产生AOM驱动参考信号fAOM；
[0038] 信号采集模块，用于将所述载波调制模拟信号s(t)转换为载波调制数字信号s(n)，并将所述AOM驱动参考信号fAOM转换为数字信号；
[0039] 本振信号模块，用于基于数字fAOM信号产生两路单位正交信号I0和Q0，其中所述两路单位正交信号I0和Q0的信号频率与载波调制数字信号s(n)的信号频率相同，数据速率与信号采集模块的采样率相同；
[0040] 混频模块，用于将载波调制数字信号s(n)与所述两路单位正交信号I0和Q0进行混频，得到两路正交信号IS和QS，所述正交信号IS和QS同时保留了倍频和基带成分；
[0041] 数据抽取模块，用于对所述正交信号IS和QS降低数据速率和滤波，滤除其部分倍频成分，保留基带成分；
[0042] 滤波模块，用于对所述正交信号IS和QS进一步滤波，经过滤波模块处理后，得到信号IB和QB，正交信号IB和QB正交信号只保留了基带信号；
[0043] Cordic相位计算模块，用于对滤波模块处理后的所述正交信号IB和QB的相位进行计算；
[0044] 相位补偿模块，针对目标相位信号 超出上述Cordic相位计算模块的量程[-π，π]情况下，对 进行补偿，得到具有真实相位的相位信号
[0045] 数据输出模块，用于将相位补偿模块输出的数字相位信号 转换为模拟位移信号d(t)后输出。
[0046] 作为优选，所述载波调制模拟信号 其中A为载波幅度，f0为载波频率， 为振动信号调制相位，所述载波调制模拟信号s(t)中包含有振动调制信启
[0047] 作为优选，信号采集模块将载波调制模拟信号s(t)进行离散化处理，离散化后的载波调制数字信号s(n)的表达式为 A(n)为数字化后的载波幅度，fC为数字化后的载波频率， 为振动信号调制相位；
[0048] 作为优选，所述数据采样模块的采样率FS满足f0＜FS＜2f0(f0为模拟载波频率)，则离散化后的载波频率为fC＝FS-f0。
[0049] 作为优选，离散化的信号s(n)在混频模块与本振模块产生的单位正交信号I0及Q0混频，混频公式如下：
[0050] IS路：
[0051]
[0052] QS路：
[0053]
[0054] 混频后，得到正交化的IS、QS两路信号，包括倍频项与基带项；
[0055] 作为优选，数据抽取模块滤波之后，所述正交信号IS和QS消除了部分倍频项，经过进一步滤波后，得到正交信号IB和QB，IB和QB只保留基带项，表达式为：
[0056] IB路：
[0057] QB路：
[0058] 作为优选，Cordic相位计算模块依据坐标旋转原理得到目标相位信号与反正切算法arctan(QB/IB)得到的值相同，Cordic算法相位 计算范围为[-π，π]，采用Cordic相位计算模块可以避免反正切算法中分母IB为零时无法计算的情况；
[0059] 作为优选，相位补偿模块基于Cordic相位计算模块对可能产生[-π，π]之外的信号进行补偿，根据产生折叠时相位值 与真实相位值 相差2π的关系，对相位 进行补偿，得到补偿后的相位信号 其中k为整数，k＝0，±1，±2，…)，然后数据输出模块将相位补偿模块输出的数字相位信号 转换为模拟位移信号d(t)，并且滤除噪声，得到准确的目标信号；通过振动位移与目标的关系， 得到物体振动的位移d，其中λ为激光波长。
[0060] 本发明还公开了一种用于外差干涉测振系统的信号解调方法，包括以下步骤：
[0061] 步骤1，将外差干涉光路中的光信号转换成电信号，产生载波调制模拟信号s(t)，所述载波调制模拟信号s(t)中含有目标相位信息；产生AOM驱动信号fAOM；
[0062] 步骤2，将所述载波调制信号s(t)转换为数字信号s(n)；将所述AOM驱动信号fAOM离散化成数字信号；
[0063] 步骤3，本振信号模块基于所述AOM驱动信号fAOM产生单位正交信号I0和Q0，然后单位正交信号I0和Q0与载波调制信号s(n)进行混频，得到两路正交信号IS和QS，所述正交信号IS和QS中同时保留倍频与基带成分；
[0064] 步骤4，数据抽取模块对混频后的两路正交信号IS和Qs信号进行抽取和滤波，滤除混频产生的部分高频成分，其保留基带信号；
[0065] 步骤5，滤波模块对上述得到的正交信号IS和QS整形和进一步滤波，得到正交信号IB和QB，IB和QB只保留基带信号；
[0066] 步骤6，Coridic相位计算模块对上述得到的正交信号IB和QB根据坐标旋转原理进行相位计算，得到原始相位
[0067] 步骤7，相位补偿模块对超出Cordic相位计算模块量程的信号进行补偿，得到真实相位
[0068] 步骤8，将上述得到的相位信号 转换为模拟位移信号d(t)输出。
[0069] 作为优选，步骤(1)中所述载波调制模拟信号 A为数字化后的载波幅度，f0为数字化后的载波频率， 为振动信号调制相位；
[0070] 作为优选，步骤(2)中信号采集模块将载波调制模拟信号，s(t)进行离散化处理，离散化后得到的载波调制数字信号 A(n)为数字化后的载波幅度，fC为数字化后的载波频率fC＝FS-f0， 为振动信号调制相位；
[0071] 作为优选，离散化的信号s(n)在混频模块与本振模块产生的单位正交信号I0及Q0混频，混频公式如下：
[0072] IS路：
[0073]
[0074] QS路：
[0075]
[0076] 混频后，得到正交化的IS、QS两路信号，包括倍频项与基带项；
[0077] 作为优选，数据抽取模块滤波之后，所述正交信号IS和Qs消除了倍频项，得到正交信号IB和QB，IB和QB只保留基带项，表达式为：
[0078] IB路：
[0079] QB路：
[0080] 作为优选，采用Cordic相位计算模块计算IB和QB两路信号的相位 与反正切算法arctan(QB/IB)得到的值相同，Cordic算法相位 计算范围为[-π，π]，采用Cordic相位计算模块可以避免反正切算法中分母IB为零时无法计算的情况；
[0081] 作为优选，相位补偿模块基于Cordic相位计算模块对可能产生[-π，π]之外的 信号进行补偿，根据产生折叠时相位值 与真实相位值 相差2π的关系，对相位 进行补偿，得到补偿后的相位信号 其中k为整数，k＝0，±1，±2，…，然后数据输出模块将相位补偿模块输出的数字相位信号 转换为模拟位移信号d(t)，并且滤除噪声，得到准确的目标信号；通过振动位移与目标的关系， 得到物体振动的位移d，其中λ为激光波长。
[0082] 在一个具体实施方式中，本发明的用于外差干涉测振系统的信号解调装置，包括：
[0083] 外差干涉光路，产生AOM输出的参考频率信号fAOM，产生载波调制信号s(t)，载波调制信号含有目标信号；
[0084] 信号采集模块，用于将模拟的载波调制信号s(t)转换为数字信号s(n)，便于后续数字信号处理模块处理；
[0085] 本振信号模块，用于产生特定频率的I0、Q0正交化正弦波；混频模块，用于将本振信号模块产生的I0和Q0信号和载波信号s(n)的混频，实现载波信号s(t)的正交化，得到IS和QS两路信号；
[0086] 数据抽取模块，用于对混频后的IS和QS两路信号进行抽取，降低IS和QS的数据速率；
[0087] 滤波模块，用于消除混频产生IS和QS的倍频项，保留含有振动信息的基带信号IB和QB；
[0088] Cordic相位计算模块，计算基带信号IB和QB两路信号的原始相位 相位补偿模块用于补偿超出Cordic相位计算模块量程的相位值，得到真实的相位
[0089] 信号输出模块，将数字化的目标相位信号 转换成模拟位移信号d(t)，并且滤波输出。
[0090] 在一个具体实施方法中，如图2所示，本发明的用于外差干涉信号的解调系统包括：
[0091] 外差干涉光路，提供载波调制信号和本振信号模块的参考输入，作为模拟信号向数字信号转换的信号采集模块，该信号采集模块核心单元是一个N bits的A/D芯片，采样率为FS(f0＜FS＜2f0)将f0＝40MHz载波离散成Nbits的离散信号，离散后的信号频率fC＝FS-f0；
[0092] 干涉光路中AOM频移器输出的参考频率信号fAOM，通过数据采集模块将其离散化，然后作为本振信号模块的参考输入信号，本振信号模块产生两路单位正交信号I0和Q0，其频率为fC，数据速率为FS(FS＞2fC)，用于和离散化的外差干涉信号进行混频；
[0093] 离散化的外差干涉信号s(n)和本振信号模块产生的两路正交信号I0和Q0在混频模块进行混频，得到含有2fC的高频项和相位调制的基带项，此时数据速率为FS；
[0094] 经过混频后的信号进入数据抽取模块，该模块核心单元为CIC滤波器，经过1/R的抽取，得到的数据速率为FR＝FS/R，FR越小，后续的滤波器的设计越容易，还要考虑到Nyquist定律，保证波形的平滑；在数据抽取的同时，2fC的高频成分得到滤除，只保留了相位调制的基带项；
[0095] 滤波模块，采用FIR滤波器设计，经过CIC模块抽取和滤波后的IS、QS两路信号进入FIR模块进行进一步滤波，得到平滑的低频波形，FIR通带截止频率为4kHz，这有效的滤除了目标频带信号引起的振动分量，保留基带信号IB和QB；
[0096] 经过抽取和滤波的两路正交信号进入相位计算模块，相位计算模块是是基于坐标旋转原理计算正交两路信号相位 的模块，其基于Cordic旋转坐标系的原理实现，利用反正切函数arctan(QB/IB)计算目标相位有两个缺点：在分母IB处于零点位置时不能有效的计算，而基于坐标旋转的方法很好的解决了这个问题；反正切运算arctan(QB/IB)的相位计算范围 而Cordic相位计算模块的相位计算范围为[-π，π]；
[0097] 上述方案中，所述的信号采集模块包括A/D采集芯片及其外围电路部分，混频模块，本振信号模块，CIC模块，FIR模块，相位计算模块和相位补偿模块是在FPGA中实现的，这些模块采用Verilog语言实现，数据输出模块由D/A转换芯片及外围滤波电路组成。
[0098] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合附图，做出更详细的说明：
[0099] 如图1所示，本发明的外差干涉振动测量系统有光学部分和电学部分，电学部分包括光电探测、信号调理和信号解调部分，光学部分实现光信号的分频、偏转、调制、发射、接收和干涉功能，然后经过光电转换输出载波调制信号，电学部分实现了模拟信号调理和信号解调部分。如图2所示，在一个优选实施例中，本发明提供的振动信号的解调方法是针对载波频率为40MHz的振动调制信号，目标振动频率为4kHz以下的物体振动信号的解调；
[0100] 如图2所示，载波调制信号经过信号采集系统以采样频率FS进行采样，得到数字化的信号s(n)，此时载波频率由f0＝40MHz转变为fC，fC＝FS-f0；
[0101] 将AOM频移器输出的参考信号离散化，将其作为本振信号模块的输入参考信号，本振信号模块产生两路正交信号I0和Q0，这两路正交信号频率为fC，数据速率为FS，这样才可以正确地与离散化的载波调制信号混频；
[0102] 离散化的载波调制信号s(n)和本振信号模块产生的I0和Q0两路信号混频，得到IS和QS两路正交信号，其中IS和QS含有倍频2fC信号和振动调制相位基带项[0103] 经过混频后的信号，由于数据速率过高，直接进入低通滤波器会给滤波器的设计带来一定难度，因此需要降低数据速率，将CIC滤波器抽取率设R，经过抽取后数据速率为FR＝FS/R，目标信号最大频率为4kHz；
[0104] 抽取后的信号，进入到FIR滤波器，FIR有效滤除含有2fC的倍频项，只保留了低频振动调制相位项，得到含有基带信息的正交项IB和QB；
[0105] 滤波后的两路基带信号IB和QB进入到Cordic相位计算模块，Cordic相位计算模块通过坐标旋转方法将目标相位 计算出来，由于低频干扰和Cordic相位计算模块量程[-π，π]的限制，当相位 超出这个范围时会产生折叠效应，计算出的相位依然落在[-π，π]范围内；
[0106] 针对超出上述相位计算范围的相位值，需要相位补偿模块进行相位补偿，通过折叠值 与真实值 相差2π的关系，通过 (k为整数，k＝0，±1，±2，…)进行补偿，得到真实的振动信号相位值 经过公式 得到振动
的位移值，这个位移值d(n)包含了目标振动信号和低频干扰。
[0107] 数据输出模块将解调后的数字信号d(n)转换成模拟信号d(t)输出，并且进行滤波，滤除位移夹杂的噪声或异常点，得到平滑的输出波形；至此，完成整个算法的实现。
[0108] 本发明的上述算法可以通过FPGA编程实现，也可以根据其原理移植到其它通用平台上，只要其中的主要步骤与本申请技术方案一致，均应该包含在本发明的保护范围之内。
[0109] 也即是说，以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
