CN105486331B - 一种具有高精度的光学信号相位解调系统及解调方法 - Google Patents
一种具有高精度的光学信号相位解调系统及解调方法 Download PDFInfo
- Publication number
- CN105486331B CN105486331B CN201510869444.3A CN201510869444A CN105486331B CN 105486331 B CN105486331 B CN 105486331B CN 201510869444 A CN201510869444 A CN 201510869444A CN 105486331 B CN105486331 B CN 105486331B
- Authority
- CN
- China
- Prior art keywords
- phase
- signal
- submodule
- roads
- module
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000003287 optical effect Effects 0.000 title claims abstract description 37
- 230000011664 signaling Effects 0.000 title claims abstract description 22
- 238000001914 filtration Methods 0.000 claims abstract description 12
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 239000000835 fiber Substances 0.000 abstract description 13
- 230000008569 process Effects 0.000 abstract description 8
- 230000009467 reduction Effects 0.000 abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 5
- 239000013307 optical fiber Substances 0.000 abstract description 4
- 230000001186 cumulative effect Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 5
- 239000000919 ceramic Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000007667 floating Methods 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000004579 marble Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000002463 transducing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35306—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了一种具有高精度的光学信号相位解调系统及解调方法。包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块。本发明能够进行独立的不相关解调，最终解调结果中，光路中实相位变化成相关性，解调过程及外界引入噪声没有相关性，可以通过累加平均降噪的方法获得更高的信噪比；本发明有效提高光学信号相位解调精度，可广泛用于高精度光纤测量和光纤传感等领域。
Description
技术领域
本发明属于光学(纤)干涉仪测量领域，尤其涉及一种具有高精度的光学信号相位解调系统及解调方法。
背景技术
光纤传感器在振动信号探测方面具有灵敏度高，动态范围大，布设方便等优点，其中，相位型光纤传感器广泛应用于水声信号监测，地震信号记录等领域内。相位型光纤传感器的工作原理即通过测量应力，应变，位移以及加速度等物理量所引起的光路中的相位变化来标定对应的物理量，如何准确测得相位信息成为了传感器的关键所在。1982年美国海军实验室KIRKENDALL C K，DANDRIDGE A.等人提出了一种相位生成载波(PGC)的方法计算干涉光路的相位变化，该方法实现简单，计算精度高，广泛应用于干涉型光纤传感器中，但是又有着严重的问题。传统的PGC算法受环境影响很大，包括光源的抖动，调制器状态的漂移等问题。在此基础上，人们致力于如何抑制外界的环境噪声，保证获取更为精准的测量信号。
清华大学在如何抑制噪声这方面做了很多研究，施清平发表的“一种消除伴生调幅的光源调频型相位生成载波解调方法”提出了利用3×3耦合器构造辅助光路的办法，引入额外的参考信号消除掉伴生调幅影响；张敏等人提出了一种去相关的方法(CN201110191719.4)，通过在PGC系统中引入一个3×2耦合器，利用另外两路参考信号之间存在固定相位差的特点实现消除同源噪声；美国David B.Hall,La Crescenta,Calif等人提出了一种多路归一化的方法(US006134015A)，通过同时测量6路相同的干涉仪数据，然后将数据做归一化处理得到更好的噪声结果。以上几种方法核心思想是利用多路信号之间的相关性，将不相关的噪声进行剔除，从而获得更好的输出结果，但是在构造相关信号方面都增加了光信号路数，无疑会带来额外的环境影响因素；
另外美国的Antonio L.Deus，Gerald L.Assard等人提出一种改进型交叉微分相乘方法(US6580314B1)，该方法不必额外增加光路结构，利用一路传感信号中的基频分量，倍频分量开方结果引入求解过程，消除部分噪声，但是这种方法容易受伴生调幅影响，当基频分量，倍频分量波动很大时，系统输出质量会下降；哈尔滨工程大学杨军等人也同时提出一种利用椭圆拟合矫正的方法(CN 201510293443.9)以及一种增大动态范围、增强稳定性的光学干涉相位解调方法(201510293443.3)，来修正光源伴生调幅影响。以上这两种方法都在光路结构上同时引入了3×3耦合器，利用另外一路或多路输出作为辅助计算信号，进行相关运算，抑制噪声。
发明内容
本发明的目的是提供一种能够进行独立的不相关解调的，一种具有高精度的光学信号相位解调系统。本发明的目的还包括提供一种能够增强稳定性和精度的，一种具有高精度的光学信号相位解调方法。
一种具有高精度的光学信号相位解调系统，包括前端采集模块、多路锁相模块和数据解调模块，
前端采集模块包括数据采集子模块和调相波子模块，调相波子模块输出N路信号，包括基频信号cos(ω0t)、倍频信号cos(2ω0t)……N阶频信号cos(Nω0t)，数据采集子模块采集并输出PGC干涉信号；
多路锁相模块包括N个乘法器和N个滤波器，PGC干涉信号分别与调相波子模块输出的N路信号中的一路信号经过一个乘法器后传入一个滤波器中，得到N路信号的谐波成分；
数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块，N路信号的谐波成分分别经过一个希尔伯特变换子模块后传入一个相位累加子模块，N个相位累加子模块将输出信号传送给降采样平滑滤波子模块，得到最终相位进行输出。
一种具有高精度的光学信号相位解调方法，包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块；
a、调相波子模块生成输出N路信号，包括基频信号cos(ω0t)、倍频信号cos(2ω0t)……N阶频信号cos(Nω0t)，数据采集子模块采集并输出PGC干涉信号：
其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω0为调制信号频率，为待测量相位信号；
b、将调相波子模块生成输出的N路信号分别与PGC干涉信号进行相乘，再进行滤波操作，得到锁相后N路信号的谐波成分，包括基波成分、二次谐波成分……N次谐波成分；
c、调整调制深度C，得到调整后的N路信号的谐波成分；
d、对调整后的N路信号的谐波成分进行希尔伯特变换操作；
e、求取N路信号被测相位值；
f、对求得的相位值使用降采样平滑的方法得到最终相位输出。
本发明一种具有高精度的光学信号相位解调方法，还可以包括：
1、步骤b中调整后的基波成分为：
二次谐波成分为：
N次谐波成分分别为：
其中，Jk(C)为贝塞尔函数系数，k为信号高阶分量，A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。
2、步骤c中调整后的基波成分为：
二次谐波成分为：
N次谐波成分为：
3、步骤e中N路信号被测相位值的第一路被测相位信号为：
第二路被测相位信号为：
第N路相位被测信号为：
4、步骤f中，最终相位为：
有益效果：
本发明公开一种增强稳定性，抑制系统噪声的光学测量方法，该方法针对以上方法的不足，不需要额外增加参考光路，使用数学变换的方法引入相关信号，进行相关去噪运算，在控制传感器制作成本以及体积外形基础上，提升了解调信噪比；该方法包括前端采集模块，其作用是采集光路信号至解调装置内并输入至解调算法；多路锁相模块，其作用是根据干涉信号频率分量的不同提取各个频率成分，构造多路相位解调信号；数据解调模块，其作用是利用希尔伯特变换构造正交分量进行相位解调。整个过程为，光纤干涉仪输出信号经光电探测器转换进入光学测量系统内，将该路信号作为原始基准，进行多路不同频次分量锁相，得到各种频次信号分量，每一路分量可根据希尔伯特变换得到自身正交信号，通过反正切查表计算的方法，进行独立的不相关解调；最终解调结果中，光路中实相位变化成相关性，解调过程及外界引入噪声没有相关性，可以通过累加平均降噪的方法获得更高的信噪比；本方法有效提高光学信号相位解调精度，可广泛用于高精度光纤测量和光纤传感等领域。
与现有技术相比，本发明的优点在于：
1)在不改变原有的光路结构、不提升硬件性能的情况下，有效的消除伴生调幅，抑制系统噪声，提高信噪比N倍。
2)提高系统稳定性与相位分辨率，在原有的光路基础上可以有效消除光源调制引起的系统状态抖动。
3)灵活性高，根据传感器需求修改算法层，即可完成对传感器信号处理性能的优化。
附图说明
图1为一种改进相位解调算法的流程图；
图2为改进相位解调方法实验装置图；
图3为希尔伯特变换及信号求解过程详细流程图；
图4为N路叠加降噪效果图，图4(a)为N＝10时的N路叠加降噪效果图，图4(b)为N＝100时的N路叠加降噪效果图。
具体实施方式
下面将结合附图对本发明做进一步详细说明。
1.一种提高光学信号相位解调精度的方法，包含前端采集模块10，多路锁相模块11，数据解调模块12。
1)调相波子模块101输出调相波信号cos(ω0t)至光源调制器223，同时生成本地载波信号cos(ω0t)即基频信号103，倍频载波cos(2ω0t)即倍频信号104，以此规律生成信号cos(Nω0t)即N阶频信号105，共计N路；
2)数据采集子模块102采集放大器222输出，其形式为：
即PGC干涉信号106，其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω0为调制信号频率，为待测量相位信号；将公式(1)所示进行贝赛尔展开，可以得到PGC干涉信号106的频谱成分为：
其中Jk(C)为贝塞尔函数系数，k为信号高阶分量。
3)将基频信号103，倍频信号104……N阶频信号105分别与PGC干涉信号106进行相乘运算，之后进行滤波运算可以得到锁相后的基波成分117，二次谐波成分118……N次谐波成分119，其形式为：
其中A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。
4)调整调制深度C，使J1(C)≈J3(C)、J0(C)≈J2(C)，Jp(C)≈Jq(C)就可以得到修正基波成分117，二次谐波成分118……N次谐波成分119，再经过高通滤波处理即消除伴生调幅系数m对解调结果的影响，得到
5)对调整后的基波成分117即公式(6)，二次谐波成分118即公式(7)……N次谐波成分119即公式(8)进行希尔伯特变换操作，设任意需要进行变换的信号为x(t)，对应输出结果为有如下关系：
所以在公式(6)(7)(8)输入后即可得到对应的正交分量，使用反正切求解的方法可以分别求得N路信号被测相位值：
其中分别为第一路被测相位信号3111，第二路被测相位信号3112……第N路相位被测信号3113，H为希尔伯特变换过程。最后在输出末端使用降采样平滑的方法得到最终输出结果：
2.所述的前端采集模块10，多路锁相模块11，数据解调模块12，其内部参数设定：
1)调相波子模块101，根据光源调制器223增益倍数输出幅度为1rad-10rad，频率ω0为20kHz-2MHz的正弦信号进行载波调制；调相波子模块101输出N路本地调相波信号，信号频率为ω0，2ω0……Nω0不等，根据系统采样率取N＝1-104皆可；
2)数据采集子模块102根据调相波信号频率范围，设定采样率为1MHz-100MHz；
3)多路锁相模块11中第一滤波器114，第二滤波器115以及第N滤波器116为低通滤波器，截止频率根据调相波信号频率范围选定为10kHz-2MHz不等，衰减速度至少为-80dB/十倍频。
4)数据解调模块12中第一希尔伯特变换子模块121，第二希尔伯特变换子模块122以及第N希尔伯特变换子模块123根据采样率匹配运行速率为10MHz-500MHz。
●发明原理：
本算法用于解决干涉仪的相位解调，干涉仪基本结构c示，主要包括以下几个模块：光源模块20，干涉仪21，探测及控制模块22，数字解调模块23。
系统工作开始，光源201输出调制后的稳定干涉信号，经过隔离器202并通过衰减器203注入到干涉仪21中。其中隔离器202用来防止光路中的回光反射对光源造成伤害，衰减器203用于调节光路的输出功率，控制光路的输出光强在被采样的时候能有更好的信噪比。光源模块20输出光射入环形器211中，从环形器A口输入通过B口注入到2×2耦合器212的一臂中，经过2×2耦合器212分成两路；第一路光经过光纤环213到达第一法拉第旋镜214处，被法拉第旋镜反射回2×2耦合器212中；第二路光经过压电陶瓷环216到达第二法拉第旋镜215处，同样被反射回2×2耦合器212中，并与第一路光在此发生干涉。这两个法拉第旋镜在此作用是消除光路中的偏振衰落现象。干涉信号分别由2×2耦合器两路输出，一路光经过环形器211的B口并由C口输出至差分探测器221处；另一路信号直接输出至差分探测器211处；由于这两路信号具有半个周期的相位差，在差分后可以去掉一部分共光路的噪声。
差分探测器221的输出连接至放大器222，将前端光信号调理至采集系统动态范围内的电压信号。放大器222的增益倍数由增益控制器224控制，增益控制器根据采集模块231的动态范围与实时信号值大小调节光信号增益倍数。放大器222输出连接至采集模块231，采集到的电压值送入计算机232进行处理，同时计算机231将调相波信号送入采集模块231中。为了能得到稳定可靠的高信噪比信号，需要设定调制深度在2-10rad的一个合理范围内，即通过光源调制器223，使光路的调相波在光路中的相位变化幅度为2-10rad。由于光路输出带有调相波信号，在放大器222的输出端得到干涉信号如公式(1)所示；经过采集模块231将数据送入计算机232中进行算法解调。
●算法流程结构
该算法由前端采集模块10，多路锁相模块11，数据解调模块12组成，其中前端采集模块10采集干涉信号，经过多路锁相模块11与数据解调模块12得到系统输出。
前端采集模块10中调相波子模块101用于产生调相波与本地载波，数据采集子模块102用于采集PGC干涉信号106。
多路锁相模块11完成对PGC干涉信号106不同信号频率分量提取功能，经过第一乘法器111，第二乘法器112……第N乘法器113等完成信号混频过程；经过第一滤波器114，第二滤波器115……第N滤波器116等实现PGC干涉信号106不同信号频率分量提取。
数据解调模块12完成希尔伯特变换功能；通过第一希尔伯特变换子模块121，第二希尔伯特变换子模块122……第N希尔伯特变换子模块123对基波成分117，二次谐波成分118……N次谐波成分119进行希尔伯特变换操作。
●希尔伯特变换解调原理
基波成分117，二次谐波成分118，N次谐波成分119经过修正后，通过第一高通滤波器301、第二高通滤波器303以及第N高通滤波器306得到N组交流信号；希尔伯特变换过程为自身与希尔伯特变换常量305作卷积运算，根据希尔伯特变换原理如公式(9)所示。设原始信号为h(t)则
h(t)＝1/πt (14)
其希尔伯特变换的幅度，相位特性为
即通过希尔伯特变换为幅频特性为1的全通滤波器。任意信号通过希尔伯特变换器后，其负频率成分作+90°相移，而正频率成分作－90°相移。
基波变换分量311与基波成分117同时经过第一除法器314得到被测相位信号的正切值；二次谐波变换分量312与二次谐波成分118同时通过第二除法器315也可得到被测信号的正切值，依次类推，有N路信号。这几路信号分别由第一反正切子模块317，第二反正切子模块318，第N反正切子模块319求解得到对应的第一路被测相位信号3111与第二路被测相位信号3112，第N路被测相位信号3113；
●多路降噪取平均原理
当系统共有N路信号参与解调时，即有N路相关的信号，但是每路信号分别带有不相关的噪声值，附图4说明N路取平均降噪结果，以下是证明过程：
设任意信号x(t)，其对应的傅立叶变换为：
当取一定长度信号做傅里叶变换，得到对应频谱的幅度值为：
若N路信号相关，则将N路信号相加，可以得到N路信号累加后的频谱为：
此时，求该信号的功率谱为傅里叶变换模值的平方：
即N路信号叠加，信号功率为原单一信号功率的N2倍；
但是对于噪声信号而言，N路噪声信号能量不同，根据Parseval定理，信号能量守恒，时域能量等于频域能量，不会因为变换而发生改变，即函数平方的和(或积分)等于其傅里叶转换式平方之和(或者积分)。此时，我们已知噪声为随机信号，服从概率分布，所以噪声功率为：σ2；此时可以计算出其对应的频域功率大小为：
当N路随进信号叠加，其噪声功率和为
所以当信号路数叠加至N路情况下，其信噪比为：
由此可以得出结论，在信号路数经过希尔伯特变换被拓展为N路情况下，其信噪比为原来的N倍。
一种提高光学信号相位解调精度的方法，光纤干涉仪输出信号经光电探测器转换进入光学信号相位解调处理系统内，采集模块通过合理调节光源调制器增益倍数，抑制伴生调幅现象的影响；相位计算过程中通过多通道，多频次锁相的办法构造出多路测量信号，每一路信号分别进行希尔伯特变换得到自身的正交分量，各路信号通过与自身正交分量进行反正切查表计算分别得到相位变化量，这些相位变化量呈相关性，且带有一定的共模噪声，通过平均平滑的方法抑制共模噪声，使输出获得更好的信噪比；本方法有效提高光学信号相位解调精度，可广泛用于高精度光纤测量和光纤传感等领域。
一种提高光学信号相位解调精度的方法，包含前端采集模块10，其特征是多路锁相模块11，数据解调模块12。
1)调相波子模块101输出调相波信号cos(ω0t)至光源调制器223，同时生成本地载波信号cos(ω0t)即基频信号103，倍频载波cos(2ω0t)即倍频信号104，以此规律生成信号cos(Nω0t)即N阶频信号105，共计N路；
2)数据采集子模块102采集放大器222输出，其形式为：
即PGC干涉信号106，其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω0为调制信号频率，为待测量相位信号；将公式(23)所示进行贝赛尔展开，可以得到PGC干涉信号106的频谱成分为：
其中Jk(C)为贝塞尔函数系数，k为信号高阶分量。
3)将基频信号cos(ω0t)103，倍频信号cos(2ω0t)104……N阶频信号cos(Nω0t)105分别与PGC干涉信号106进行相乘运算，之后进行滤波运算可以得到锁相后的基波成分117，二次谐波成分118……N次谐波成分119，其形式为：
其中A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。
4)调整调制深度C，使J1(C)≈J3(C)、J0(C)≈J2(C)，Jp(C)≈Jq(C)就可以得到修正后的基波成分117，二次谐波成分118……N次谐波成分119，再经过高通滤波处理即消除伴生调幅系数m对解调结果的影响，得到
5)对调整后的基波成分117即公式(28)，二次谐波成分118即公式(29)……N次谐波成分119即公式(30)进行希尔伯特变换操作，设任意需要进行变换的信号为x(t)，对应输出结果为有如下关系：
所以在公式(28)(29)(30)输入后即可得到对应的正交分量，使用反正切求解的方法可以分别求得N路信号被测相位值：
其中分别为第一路被测相位信号3111，第二路被测相位信号3112……第N路相位被测信号3113，H为希尔伯特变换过程。最后在输出末端使用降采样平滑的方法得到最终输出结果：
前端采集模块10，多路锁相模块11，数据解调模块12，其内部参数是：
1)调相波子模块101，根据光源调制器223增益倍数输出幅度为1rad-10rad，频率ω0为20kHz-2MHz的正弦信号进行载波调制；调相波子模块101输出N路本地调相波信号，信号频率为ω0，2ω0……Nω0不等，根据系统采样率取N＝1-104皆可；
2)数据采集子模块102根据调相波信号频率范围，设定采样率为1MHz-100MHz；
3)多路锁相模块11中第一滤波器114，第二滤波器115以及第N滤波器116为低通滤波器，截止频率根据调相波信号频率范围选定为10kHz-2MHz不等，衰减速度至少为-80dB/十倍频。
4)数据解调模块12中第一希尔伯特变换子模块121，第二希尔伯特变换子模块122以及第N希尔伯特变换子模块123根据采样率匹配运行速率为10MHz-500MHz。
实施例——增加N路抑制系统噪声
干涉仪装置如图2所示，干涉仪测量装置的器件选择与参数如下：
1.光源201的中心波长1550nm、半谱宽度大于45nm，出纤功率大于1～10mW；
2.光纤隔离器202工作波长1550nm±5nm，插入损耗≤1.0dB(23℃工作温度时)，回波损耗≥55dB；
3.环形器211工作波长1550nm&1310nm，插入损耗1.0dB，隔离度28dB，方向性50dB，工作温度0～70℃，回波损耗45dB；
4.第一法拉第旋镜214、第二法拉第旋镜215工作波长1550nm±5nm，插入损耗0.6dB，法拉第旋转角度90°，旋转角误差@23℃为±1°，最大光源承受能力1W，工作温度-40至85℃；
5.用于加载标定信号的压电陶瓷尺寸为24mm，电容22nF，耐压幅度0～120V；
6.2×2耦合器212工作波长1550nm，使用2×2耦合器的1端口与2端口输入，1端口输入对应2路输出分光比为49.8％，50.2％，2端口输入对应2路输出分光比为49.6％，50.4％；
7.差分探测器221为InGaAs型光电探测器，连接模式属于尾纤式FC/PC，工作波长为1100nm～1650nm，光强响应度R＝0.85A/W，电容为0.35pF；
8.放大器222用于放大转换后的光电压信号，工作带宽为200kHz，工作过程包括使用MSP430单片机采集信号幅度，调节信号增益，保证信号幅度并不会饱和；
9.采集模块231为NI-6366采集卡，采样率为2Mbps，同步采集，输入电压幅度±10V，采样时钟为采集卡内部时钟，三路同步误差小于10ns，输入电阻20kΩ；
10.压电陶瓷驱动器233为功率放大器，使用AD公司的AD8040轨对轨功率放大器，工作电压2.7V～12V，工作带宽125MHz，最大输出电流200mA，负载电容15pF；
相位解调系统的调制深度稳定测试具体流程如附图1和图3所示：
1.系统启动，信号采集模块10产生载波调制光源，采样率为2Mbps，载波频率为20kHz，调制波形中混入调制深度C值得抖动与变化，将调制深度变化波动量设置为±10％，压电陶瓷产生标定信号，频率为10Hz，随着调制电压增加，产生光路相位变化10-5rad～105rad；
2.经过数据采集子模块102得到PGC干涉信号106，调相波子模块101产生基频信号103，倍频信号104，设置本地调相波信号幅度为1V，频率为20kHz，采样率为2MHz；
3.锁相滤波模块11中，设置第一、第二、第N滤波器为FIR等纹波滤波器，参数为通带10kHz，截阻带12kHz，衰减-120dB，通带纹波为0.01dB，阶数为764阶，数据经过滤波器后得到两路正交信号。
4.设置第一，第二、第N乘法器为有符号同步乘法器，输入数据为16bit单精度有符号浮点数据，输出数据为32bit有符号双精度浮点类型数据。
5.设置第一，第二、第N高通滤波器为FIR等纹波滤波器，参数为通带1Hz，截阻带1Hz，衰减-120dB，通带纹波为0.01dB，阶数为953阶，数据经过滤波器后得到两路交流信号。
6.设置第一，第二、第N卷积器为串行工作方式，卷积常量根据实际信号幅度大小调节，卷积器工作时钟频率与乘法器相同。
7.设置第一，第二、第N除法器为有符号除法器，输入数据为16bit单精度有符号浮点数据，输出数据为32bit有符号双精度浮点类型数据，信号经过除法器得到消掉光强波动与调制量波动的信号比值。
8.设置第一，第二、第N相位累加器为串行工作模式，相位累计器内部采用二分查表法与8分圆区间工作模式，保证信号每次变换不会超过系统解调的动态范围。
9.设置平均降噪子模块为并行数据输出模式，对N路输入信号取平均值得到去掉共模噪声的输出结果。
10.经过改进后的光学干涉相位解调方法，当变换路数N＝10时，信噪比为26.9dB，当拓展计算路数N＝100时，信噪比为38.3dB，增加10dB左右。
Claims (6)
1.一种具有高精度的光学信号相位解调系统，其特征在于：包括前端采集模块、多路锁相模块和数据解调模块，
前端采集模块包括数据采集子模块和调相波子模块，调相波子模块输出N路信号，包括基频信号cos(ω0t)、倍频信号cos(2ω0t)……N阶频信号cos(Nω0t)，数据采集子模块采集并输出PGC干涉信号；
多路锁相模块包括N个乘法器和N个滤波器，PGC干涉信号分别与调相波子模块输出的N路信号中的其中一路信号经过一个乘法器后传入一个滤波器中，得到N路信号的谐波成分；
数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块，N路信号的谐波成分分别经过一个希尔伯特变换子模块后传入一个相位累加子模块，N个相位累加子模块将输出信号传送给降采样平滑滤波子模块，得到最终相位进行输出。
2.一种具有高精度的光学信号相位解调方法，其特征在于：包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块；
a、调相波子模块生成输出N路信号，包括基频信号cos(ω0t)、倍频信号cos(2ω0t)……N阶频信号cos(Nω0t)，数据采集子模块采集并输出PGC干涉信号：
其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω0为调制信号频率，为待测量相位信号；
b、将调相波子模块生成输出的N路信号分别与PGC干涉信号进行相乘，再进行滤波操作，得到锁相后N路信号的谐波成分，包括基波成分、二次谐波成分……N次谐波成分；
c、调整调制深度C，得到调整后的N路信号的谐波成分；
d、对调整后的N路信号的谐波成分进行希尔伯特变换操作；
e、求取N路信号被测相位值；
f、对求得的相位值使用降采样平滑的方法得到最终相位输出。
3.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤b中调整后的基波成分为：
二次谐波成分为：
N次谐波成分分别为：
其中，Jk(C)为贝塞尔函数系数，k为信号高阶分量，A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。
4.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤c中调整后的基波成分为：
二次谐波成分为：
N次谐波成分为：
5.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤e中N路信号被测相位值的第一路被测相位信号为：
第二路被测相位信号为：
第N路相位被测信号为：
6.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤f中，最终相位为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510869444.3A CN105486331B (zh) | 2015-12-01 | 2015-12-01 | 一种具有高精度的光学信号相位解调系统及解调方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510869444.3A CN105486331B (zh) | 2015-12-01 | 2015-12-01 | 一种具有高精度的光学信号相位解调系统及解调方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105486331A CN105486331A (zh) | 2016-04-13 |
| CN105486331B true CN105486331B (zh) | 2017-12-19 |
Family
ID=55673465
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510869444.3A Active CN105486331B (zh) | 2015-12-01 | 2015-12-01 | 一种具有高精度的光学信号相位解调系统及解调方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105486331B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107104740A (zh) * | 2017-01-16 | 2017-08-29 | 中国科学院声学研究所 | 一种信号处理方法及系统 |
| CN108007340B (zh) * | 2017-10-24 | 2019-12-06 | 浙江理工大学 | 相位生成载波反正切解调中非线性误差的实时计算方法 |
| CN110530496B (zh) * | 2018-05-25 | 2021-05-07 | 南京理工大学 | 一种基于数字滤波和二元脉冲调制的光纤传感解调方法 |
| CN109115722B (zh) * | 2018-06-22 | 2023-10-31 | 山东建筑大学 | 一种高灵敏相位解调光纤气体传感系统 |
| JP7020385B2 (ja) * | 2018-11-29 | 2022-02-16 | 日本電信電話株式会社 | 位相測定方法、信号処理装置、およびプログラム |
| CN110319826B (zh) * | 2019-06-24 | 2021-06-11 | 北京控制工程研究所 | 基于自适应滤波的光纤陀螺阶梯波串扰抑制方法 |
| CN111561953B (zh) * | 2020-05-14 | 2021-11-05 | 上海科技大学 | 基于波分复用与平衡探测的片上集成光学矩阵矢量乘法器 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6580314B1 (en) * | 2000-10-11 | 2003-06-17 | The United States Of America As Represented By The Secretary Of The Navy | Demodulation system and method for recovering a signal of interest from a modulated carrier sampled at two times the phase generated carrier frequency |
| CN102169273A (zh) * | 2011-04-15 | 2011-08-31 | 天津工业大学 | 非对称性3×3耦合器输出信号的无载波解调方法 |
| WO2012147324A1 (ja) * | 2011-04-28 | 2012-11-01 | 日本電気株式会社 | 光周波数変換方法および装置 |
| CN202267808U (zh) * | 2011-06-20 | 2012-06-06 | 深圳职业技术学院 | 干涉型光纤传感器数字解调装置 |
| CN102353393B (zh) * | 2011-07-01 | 2014-08-20 | 浙江大学 | 基于π/2相位调制的干涉型光传感器的正交解调装置 |
| CN102359797B (zh) * | 2011-07-08 | 2014-05-07 | 清华大学 | 用于光纤水听器中pgc系统降噪的系统及方法 |
| CN102624456B (zh) * | 2012-02-21 | 2014-07-23 | 中国科学院半导体研究所 | 一种光纤干涉型传感信号解调的系统及方法 |
| CN102620757A (zh) * | 2012-03-27 | 2012-08-01 | 中国科学院半导体研究所 | 一种对光纤干涉型传感信号进行解调的系统及方法 |
| CN103592805B (zh) * | 2013-11-28 | 2015-11-11 | 武汉光迅科技股份有限公司 | 一种相位调制信号的全光波长转换装置 |
| CN104197966B (zh) * | 2014-08-20 | 2017-03-01 | 中国科学院半导体研究所 | 一种抗调制深度漂移的光纤光栅的解调方法 |
| CN105067017B (zh) * | 2015-06-02 | 2017-11-28 | 哈尔滨工程大学 | 一种改进的生成载波相位pgc解调方法 |
-
2015
- 2015-12-01 CN CN201510869444.3A patent/CN105486331B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105486331A (zh) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105486331B (zh) | 一种具有高精度的光学信号相位解调系统及解调方法 | |
| CN105486225B (zh) | 一种抑制光强波动噪声的相位解调装置及解调方法 | |
| CN105067017B (zh) | 一种改进的生成载波相位pgc解调方法 | |
| CN105157733B (zh) | 一种改进的生成载波相位pgc解调方法 | |
| CN112097813B (zh) | 一种基于光频调制的3×3耦合器光电检测方法与装置 | |
| CN102353393B (zh) | 基于π/2相位调制的干涉型光传感器的正交解调装置 | |
| CN107340050B (zh) | 一种光纤分布式振动传感系统及鉴相非线性误差修正方法 | |
| CN109459070A (zh) | 一种pgc相位解调法中相位延迟提取与补偿方法 | |
| CN106992818A (zh) | 一种基于相位程控的相位生成载波解调装置及方法 | |
| CN108007550A (zh) | 一种改进的pgc调制解调检测方法 | |
| CN105203136B (zh) | 一种基于差分放大技术的分布式传感系统 | |
| CN109141490A (zh) | 一种扰动波形和位置同时测量的光纤传感装置和解调方法 | |
| CN103117973A (zh) | 基于正交检波的多路并行解调系统 | |
| CN112698253B (zh) | 一种数字化三轴tmr磁传感系统 | |
| CN110411334B (zh) | 一种改进的相位载波pgc解调方法及系统 | |
| CN101799610A (zh) | 外差式相位干涉型光纤传感器的正交解调装置 | |
| CN108760021A (zh) | 基于双折射晶体的光纤法珀声振动传感装置及解调方法 | |
| CN104729493A (zh) | 一种光纤陀螺仪的新型探测方法 | |
| CN109361636A (zh) | 一种基于相位生成载波的解调装置及方法 | |
| CN110332980A (zh) | 一种双路正交相位光纤声振动传感器解调装置及解调方法 | |
| CN209002004U (zh) | 一种基于相位生成载波的解调装置 | |
| CN111366179A (zh) | 一种自微分相除的相位生成载波解调方法 | |
| CN103322926B (zh) | 信号传输过程中的周期性非线性误差或干扰消除法 | |
| CN100386601C (zh) | Sagnac干涉型光纤传感器的一种开环信号处理方法和电路 | |
| CN206850776U (zh) | 一种基于相位程控的相位生成载波解调装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |