CN107421435B - 基于软件移相的零差正交激光干涉信号生成方法 - Google Patents

基于软件移相的零差正交激光干涉信号生成方法 Download PDF

Info

Publication number
CN107421435B
CN107421435B CN201710183123.7A CN201710183123A CN107421435B CN 107421435 B CN107421435 B CN 107421435B CN 201710183123 A CN201710183123 A CN 201710183123A CN 107421435 B CN107421435 B CN 107421435B
Authority
CN
China
Prior art keywords
signal
phase
photosignal
revise
period
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710183123.7A
Other languages
English (en)
Other versions
CN107421435A (zh
Inventor
何闻
谢智郎
郝凌凌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710183123.7A priority Critical patent/CN107421435B/zh
Publication of CN107421435A publication Critical patent/CN107421435A/zh
Application granted granted Critical
Publication of CN107421435B publication Critical patent/CN107421435B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/0201Interferometers characterised by controlling or generating intrinsic radiation properties using temporal phase variation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Communication System (AREA)

Abstract

基于软件移相的零差正交激光干涉信号生成方法,包括从激光干涉仪按整周期采样获取单路光电信号;将光电信号的幅值按比例缩放到峰值为1V,同时去除直流量得到幅值均值为0的归一化信号;对归一化信号进行自适应修正以获得初始值位于相位反转点,且整个信号的满足整周期采样特性的修正信号;对修正信号以半周期为单位进行+90°和‑90°交替移相,获得移相信号;移相信号和来自激光干涉仪的单路光电信号形成两路正交光电信号。本发明的优点在于无需使用硬件移相,且两路光电信号的相位始终是相差+90°或者‑90°。

Description

基于软件移相的零差正交激光干涉信号生成方法
技术领域
本发明涉及一种基于软件移相的零差正交激光干涉信号生成方法,可实现利用单路激光干涉信号生成两路零差正交激光干涉信号,用于振动传感器幅频特性和相频特性绝对法测量。
技术背景
现在的国家标准中,对于振动传感器幅频特性和相频特性绝对法测量推荐使用的是零差正交激光干涉仪,如中国计量科学研究院的基于正交输出的零差激光干涉仪。但是由于上述激光干涉仪需要两路正交光电信号输出,这就需要用到四分之一波片等光学移相元器件,但是由于零差正交激光干涉仪中四分之一波片、偏振分光棱镜的安装调整位置等因素的影响可能会导致其输出的两路正交光电信号的相位差与90°会产生偏差,当与90°名义角度的偏差超过±20°的时候两路正交光电信号所组成的李萨如图则是一个离心率较大的椭圆,这样的两路光电信号不能正确地进行后续的解调处理。并且,从硬件上实现移相,再加上附加的一套光电转换设备及高速数据采集系统,造价昂贵。
发明内容
本发明的目的在于提供一种无需使用光学移相元器件,并能保证两路正交光电信号不存在正交相位误差的基于软件移相的零差正交激光干涉信号生成方法。
基于软件移相的零差正交激光干涉信号生成方法,包括以下步骤:
步骤1:从激光干涉仪按被测振动信号的周期做整周期采样获取单路光电信号;
步骤2:将光电信号的幅值按比例缩放到峰值为1V,同时去除直流量得到幅值均值为0的归一化信号;
步骤3:对归一化信号进行自适应修正以获得初始值位于相位反转点,且满足整周期采样特性的修正信号;
步骤4:对修正信号以半周期为单位进行+90°和-90°交替移相,获得移相信号;
步骤5:移相信号和来自激光干涉仪的单路光电信号形成两路正交光电信号。根据得到的两路正交光电信号,可以解调得到振动信号波形,经过参数辨识后即可得到振动信号的幅值和相位,从而可以测量得到振动传感器的幅相特性。
进一步,步骤3中修正信号的获取包括以下步骤:
步骤3.1:判断归一化信号的初始值是否位于相位反转点,若是,则将归一化信号作为修正信号;若否,则进入步骤3.2;
步骤3.2:对归一化信号进行差分运算,对差分运算后的信号求上包络曲线,搜索获得包络线的第一个极小值,第一个极小值对应所要求解的相位反转点,将第一个极小值点作为初始相位点;
步骤3.3:将初始相位点之前的信号作为第一修正部,将第一修正部删除获得初步修正信号;
步骤3.4:按周期遍历初步修正信号,直到出现不满一个周期的信号部分,将该不满一个周期的信号部分作为第二修正部,将第二修正部删除。
或者:步骤3.4替换为:将第一修正部的起始点与初始信号的最后一点重合,从第一修正部的终点起向前删掉最后一个周期的信号。
进一步,为了降低噪声的影响,得到准确的极小值,步骤3.2中,对上包络曲线进行低通滤波后,再搜索极小值;
进一步,步骤4中,用傅里叶变换对修正信号移相+90°,用希尔伯特变换对修正信号移相-90°,傅里叶变换和希尔伯特变换交替进行。如:第一个单位进行一次傅里叶变换将其变换+90°,第二个单位进行一次希尔伯特变换将其变换-90°,第三个单位则又进行傅里叶变换,如此交替循环进行,直到最后结束移相处理,这样即可达到零差正交激光干涉仪两路零差正交光电信号的输出效果。
本发明的优点在于:1、用傅里叶变换和希尔伯特变换实现了对光电信号的正交移相,从而使移相后的信号与来自激光干涉仪的信号形成两路正交光电信号,节省了光学移相与其连带的硬件成本和占用空间。2、通过移相算法处理得到的两路光电信号的相位始终是相差+90°或者-90°,不会存在正交相位误差的问题。
附图说明
图1为本发明的流程图。
图2为光电信号自适应修正方法流程图。
图3为交替使用傅里叶变换和希尔伯特变换移相后的光电信号。
图4为软件仿真的两路标准的零差正交光电信号。
图5为基于本发明解调得到的加速度波形。
图6为基于零差激光干涉仪解调得到的加速度波形。
具体实施方式
实施例1:
如图1,基于软件移相的零差正交激光干涉信号生成方法,包括以下步骤:
步骤1:获取单路信号:由数据采集卡读取从迈克尔逊激光干涉仪输出的单路光电信号,该单路光电信号按被测振动信号的周期做整周期采样。
步骤2:归一化处理:由于光电转换及放大电路增益与零点漂移的影响导致输出的光电信号可能出现峰值电压不为1 V或者均值不为0的现象,为方便后续的交替移相算法处理需要对光电信号进行归一化处理,即将光电信号的幅值按比例缩放至峰值为1 V,同时去除直流量得到均值为0的信号。
步骤3:对归一化信号进行自适应修正以获得初始值位于相位反转点,且满足整周期采样特性的修正信号:寻找振动激励信号周期的反向点,即振动速度的反向点,亦称为半周期分界点,对应于光电信号即为光电信号的相位反转点。如果采集得到的光电信号的初始值正好位于相位反转点,则可以直接进入到下一步的移相处理;如果采集得到的光电信号的初始值不是相位反转点,则需要对光电信号进行修正,如图2所示,以保证进入后续的移相处理过程的光电信号的初始值为相位反转点;
在移相算法中所采集到的光电信号所携带的振动信号的初相位不一定是0°或者180°,也就是说不一定正好是从振动激励装置反向运动开始的,那么±90°交替移相就不能够从半周期单位开始了,这也就意味着疏-密-疏的半周期单位被破坏了,那么此时依然从开始来计算半个周期采样点就不能正确地得到与初始光电信号正交的光电信号,所以需要对归一化信号进行自适应修正以获得初始值位于相位反转点。
步骤3.1:判断归一化信号的初始值是否位于相位反转点,若是,则将归一化信号作为修正信号;若否,则进入步骤3.2。
步骤3.2:对归一化信号进行差分运算,对差分运算后的信号求上包络曲线,搜索获得包络线的第一个极小值,第一个极小值对应所要求解的相位反转点,将第一个极小值点作为初始相位点;
为了降低噪声的影响,得到准确的极小值点,对上包络曲线进行适当的低通滤波后,再搜索极小值。当原始信号的噪声影响可以忽略时,可以不使用低通滤波。
步骤3.3:将初始相位点之前的信号作为第一修正部,将第一修正部删除获得初步修正信号。
步骤3.4:按周期遍历初步修正信号,直到出现不满一个周期的信号部分,将该不满一个周期的信号部分作为第二修正部,将第二修正部删除。满足整周期采样特性的修正信号为步骤4的移相做准备。
步骤4:移相处理:对修正信号以半周期为单位进行+90°和-90°交替移相;用傅里叶变换对修正信号移相+90°,用希尔伯特变换对修正信号移相-90°,如此交替循环进行,直到最后移相处理,即可完成用软件来生成另一路正交光电信号,从而获得两路正交光电信号,如图3所示。图4为软件仿真的两路标准的零差正交光电信号。
这是因为零差正交激光干涉仪所输出的两路光电信号的相位差并不都是90°,而是当振动激励装置反向的时候相位差也反相,即两路信号的相位差为+90°,当振动激励装置反向时两路信号的相位差就变为-90°,当振动激励装置振动再次反向时两路信号的相位差就又变为+90°,因此在移相处理的过程中考虑对信号进行交替移相;而振动信号在一个周期内会发生两次反向,因此振动信号在半个周期内会发生一次反向运动,故而以半个周期为单位。
步骤5:移相信号和来自激光干涉仪的单路光电信号形成两路正交光电信号。根据得到的两路正交光电信号,可以解调得到振动信号波形,如图5所示(图6为基于零差激光干涉仪解调得到的加速度波形)。经过参数辨识后即可得到振动信号的幅值和相位,从而可以测量得到振动传感器的幅相特性。
实施例2
本实施例与实施例1的区别在于:步骤3中获得修正信号的具体方法不同,除了获取修正信号的步骤以外,其余步骤均与实施例1相同。
步骤3中修正信号的获取包括以下步骤:
步骤3.1:判断归一化信号的初始值是否位于相位反转点,若是,则将归一化信号作为修正信号;若否,则进入步骤3.2;
步骤3.2:对归一化信号进行差分运算,对差分运算后的信号求上包络曲线;为了降低噪声的影响,得到准确的极小值点,对上包络曲线进行适当的低通滤波后,再搜索极小值。当原始信号的噪声影响可以忽略时,可以不使用低通滤波。
步骤3.3:将初始相位点之前的信号作为第一修正部,将第一修正部截取,同时获得初步修正信号;
步骤3.4:将第一修正部的起始点与初步修正信号的最后一点重合,从第一修正部的终点起向前删掉最后一个周期的信号,获得所需的修正信号。
这种方法的优势在于:无需按周期遍历光电信号,节约运算资源。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (3)

1.基于软件移相的零差正交激光干涉信号生成方法,包括以下步骤:
步骤1:从激光干涉仪按被测振动信号的周期做整周期采样获取单路光电信号;
步骤2:将光电信号的幅值按比例缩放到峰值为1V,同时去除直流量得到幅值均值为0的归一化信号;
步骤3:对归一化信号进行自适应修正以获得初始值位于相位反转点,且满足整周期采样特性的修正信号;修正信号的获取包括以下步骤:
步骤3.1:判断归一化信号的初始值是否位于相位反转点,若是,则将归一化信号作为修正信号;若否,则进入步骤3.2;
步骤3.2:对归一化信号进行差分运算,对差分运算后的信号求上包络曲线,搜索获得包络线的第一个极小值,第一个极小值对应所要求解的相位反转点,将第一个极小值点作为初始相位点;
步骤3.3:将初始相位点之前的信号作为第一修正部,将第一修正部删除获得初步修正信号;
步骤3.4:按周期遍历初步修正信号,直到出现不满一个周期的信号部分,将该不满一个周期的信号部分作为第二修正部,将第二修正部删除。
步骤4:对修正信号以半周期为单位进行+90°和-90°交替移相,获得移相信号;用傅里叶变换对修正信号移相+90°,用希尔伯特变换对修正信号移相-90°,傅里叶变换和希尔伯特变换交替进行。
步骤5:移相信号和来自激光干涉仪的单路光电信号形成两路正交光电信号。
2.如权利要求1所述的基于软件移相的零差正交激光干涉信号生成方法,其特征在于:步骤3.1:判断归一化信号的初始值是否位于相位反转点,若是,则将归一化信号作为修正信号;若否,则进入步骤3.2;
步骤3.2:对归一化信号进行差分运算,对差分运算后的信号求上包络曲线,搜索获得包络线的第一个极小值,第一个极小值对应所要求解的相位反转点,将第一个极小值点作为初始相位点;
步骤3.3:将初始相位点之前的信号作为第一修正部,将第一修正部删除获得初步修正信号;
步骤3.4:按周期遍历初步修正信号,直到出现不满一个周期的信号部分,将该不满一个周期的信号部分作为第二修正部,将第二修正部删除。
3.如权利要求1或2所述的基于软件移相的零差正交激光干涉信号生成方法,其特征在于:步骤3.2中,对上包络曲线进行低通滤波后,再搜索极小值。
CN201710183123.7A 2017-03-24 2017-03-24 基于软件移相的零差正交激光干涉信号生成方法 Active CN107421435B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710183123.7A CN107421435B (zh) 2017-03-24 2017-03-24 基于软件移相的零差正交激光干涉信号生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710183123.7A CN107421435B (zh) 2017-03-24 2017-03-24 基于软件移相的零差正交激光干涉信号生成方法

Publications (2)

Publication Number Publication Date
CN107421435A CN107421435A (zh) 2017-12-01
CN107421435B true CN107421435B (zh) 2019-11-12

Family

ID=60424126

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710183123.7A Active CN107421435B (zh) 2017-03-24 2017-03-24 基于软件移相的零差正交激光干涉信号生成方法

Country Status (1)

Country Link
CN (1) CN107421435B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063668B (zh) * 2018-08-14 2021-09-14 哈工大机器人(合肥)国际创新研究院 一种基于峰值保留降采样的冲击信号包络解调方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7149029B1 (en) * 2005-01-11 2006-12-12 Hrl Laboratories, Llc Interferometric PPM demodulators based on semiconductor optical amplifiers
KR100830694B1 (ko) * 2007-05-22 2008-05-19 광주과학기술원 광학식 호모다인 간섭계의 전자식 헤테로다인 간섭계로의변환방법 및 이를 이용한 간섭계
CN101713685A (zh) * 2009-12-03 2010-05-26 北京理工大学 一种无源零差正交解调技术的相位补偿方法
DE112010000705B4 (de) * 2009-01-28 2014-09-25 Kabushiki Kaisha Kobe Seiko Sho Formbestimmungsvorrichtung
CN105043243A (zh) * 2015-05-29 2015-11-11 南京师范大学 正交型零差激光干涉仪及其测量方法
CN105121998A (zh) * 2012-12-13 2015-12-02 哈德斯菲尔德大学 干涉设备及使用这种设备的样品特征确定设备
CN105606194A (zh) * 2015-12-29 2016-05-25 哈尔滨工业大学 一种基于激光正交偏振干涉技术的水下声信号实时提取方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7149029B1 (en) * 2005-01-11 2006-12-12 Hrl Laboratories, Llc Interferometric PPM demodulators based on semiconductor optical amplifiers
KR100830694B1 (ko) * 2007-05-22 2008-05-19 광주과학기술원 광학식 호모다인 간섭계의 전자식 헤테로다인 간섭계로의변환방법 및 이를 이용한 간섭계
DE112010000705B4 (de) * 2009-01-28 2014-09-25 Kabushiki Kaisha Kobe Seiko Sho Formbestimmungsvorrichtung
CN101713685A (zh) * 2009-12-03 2010-05-26 北京理工大学 一种无源零差正交解调技术的相位补偿方法
CN105121998A (zh) * 2012-12-13 2015-12-02 哈德斯菲尔德大学 干涉设备及使用这种设备的样品特征确定设备
CN105043243A (zh) * 2015-05-29 2015-11-11 南京师范大学 正交型零差激光干涉仪及其测量方法
CN105606194A (zh) * 2015-12-29 2016-05-25 哈尔滨工业大学 一种基于激光正交偏振干涉技术的水下声信号实时提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于相干光传输系统的连续相位调制方法的研究;邹宁;《中国博士学位论文全文数据库信息科技辑》;20140215;第I136-26 *
激光测振在振动计量中的若干技术问题的研究;李民;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20030315;第C030-55 *

Also Published As

Publication number Publication date
CN107421435A (zh) 2017-12-01

Similar Documents

Publication Publication Date Title
Kang et al. Phase difference correction method for phase and frequency in spectral analysis
JP2019525181A (ja) 分散型光ファイバセンシングシステム及びその振動検知位置決め方法
CN106992818B (zh) 一种基于相位程控的相位生成载波解调装置及方法
CN205427032U (zh) 一种基于锁相环的弱信号检测装置
TWI429295B (zh) 用來偵測一數位資料串流的電路及相關方法
CN106895904B (zh) 一种高精度的外差式激光振动校准灵敏度相位测量方法
CN107505497B (zh) 一种无源磁电转速传感器信号峰峰值时域测量方法
CN106813564A (zh) 一种lvdt位移传感器数字化处理方法及装置
CN105157733A (zh) 一种改进的生成载波相位pgc解调方法
CN104122439B (zh) 提高相位校正精度的电能表
CN107421435B (zh) 基于软件移相的零差正交激光干涉信号生成方法
CN106872017A (zh) 一种磁电式低频振动传感器的幅相补偿方法
CN111766451B (zh) 一种高精度电容参数测试的系统及方法
US6794857B2 (en) Apparatus and method for measuring a phase delay characteristic
CN107561306A (zh) 基于转子电流的绕线式电机转速测量方法
CN107579734B (zh) 信号处理方法、装置、存储介质和处理器
CN115128159A (zh) 一种基于信号脉宽捕获的轨道缺陷快速检测方法
JP3498908B2 (ja) 交流信号の振幅サンプリング方法及び振幅検出回路
JP4225651B2 (ja) 回路素子測定器の位相誤差補正方法
JP3252641B2 (ja) 位相差測定装置
CN106644031B (zh) 一种消除光强扰动的高稳定性光纤传感装置及解调方法
CN105004259A (zh) 一种电容式mems传感器检测电路
CN212059293U (zh) 一种桥梁质量变频信号检测装置和检测系统
CN111130660B (zh) 一种相位噪声测量探头电路
JP2005517929A (ja) ゼロミキシングスペクトル分析の方法及び装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant