CN105486331A - 一种具有高精度的光学信号相位解调系统及解调方法 - Google Patents

Abstract

本发明公开了一种具有高精度的光学信号相位解调系统及解调方法。包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块。本发明能够进行独立的不相关解调，最终解调结果中，光路中实相位变化成相关性，解调过程及外界引入噪声没有相关性，可以通过累加平均降噪的方法获得更高的信噪比；本发明有效提高光学信号相位解调精度，可广泛用于高精度光纤测量和光纤传感等领域。

Description

一种具有高精度的光学信号相位解调系统及解调方法

技术领域

本发明属于光学(纤)干涉仪测量领域，尤其涉及一种具有高精度的光学信号相位解调系统及解调方法。

背景技术

光纤传感器在振动信号探测方面具有灵敏度高，动态范围大，布设方便等优点，其中，相位型光纤传感器广泛应用于水声信号监测，地震信号记录等领域内。相位型光纤传感器的工作原理即通过测量应力，应变，位移以及加速度等物理量所引起的光路中的相位变化来标定对应的物理量，如何准确测得相位信息成为了传感器的关键所在。1982年美国海军实验室KIRKENDALLCK，DANDRIDGEA.等人提出了一种相位生成载波(PGC)的方法计算干涉光路的相位变化，该方法实现简单，计算精度高，广泛应用于干涉型光纤传感器中，但是又有着严重的问题。传统的PGC算法受环境影响很大，包括光源的抖动，调制器状态的漂移等问题。在此基础上，人们致力于如何抑制外界的环境噪声，保证获取更为精准的测量信号。

清华大学在如何抑制噪声这方面做了很多研究，施清平发表的“一种消除伴生调幅的光源调频型相位生成载波解调方法”提出了利用3×3耦合器构造辅助光路的办法，引入额外的参考信号消除掉伴生调幅影响；张敏等人提出了一种去相关的方法(CN201110191719.4)，通过在PGC系统中引入一个3×2耦合器，利用另外两路参考信号之间存在固定相位差的特点实现消除同源噪声；美国DavidB.Hall,LaCrescenta,Calif等人提出了一种多路归一化的方法(US006134015A)，通过同时测量6路相同的干涉仪数据，然后将数据做归一化处理得到更好的噪声结果。以上几种方法核心思想是利用多路信号之间的相关性，将不相关的噪声进行剔除，从而获得更好的输出结果，但是在构造相关信号方面都增加了光信号路数，无疑会带来额外的环境影响因素；

另外美国的AntonioL.Deus，GeraldL.Assard等人提出一种改进型交叉微分相乘方法(US6580314B1)，该方法不必额外增加光路结构，利用一路传感信号中的基频分量，倍频分量开方结果引入求解过程，消除部分噪声，但是这种方法容易受伴生调幅影响，当基频分量，倍频分量波动很大时，系统输出质量会下降；哈尔滨工程大学杨军等人也同时提出一种利用椭圆拟合矫正的方法(CN201510293443.9)以及一种增大动态范围、增强稳定性的光学干涉相位解调方法(201510293443.3)，来修正光源伴生调幅影响。以上这两种方法都在光路结构上同时引入了3×3耦合器，利用另外一路或多路输出作为辅助计算信号，进行相关运算，抑制噪声。

发明内容

本发明的目的是提供一种能够进行独立的不相关解调的，一种具有高精度的光学信号相位解调系统。本发明的目的还包括提供一种能够增强稳定性和精度的，一种具有高精度的光学信号相位解调方法。

一种具有高精度的光学信号相位解调系统，包括前端采集模块、多路锁相模块和数据解调模块，

前端采集模块包括数据采集子模块和调相波子模块，调相波子模块输出N路信号，包括基频信号cos(ω₀t)、倍频信号cos(2ω₀t)……N阶频信号cos(Nω₀t)，数据采集子模块采集并输出PGC干涉信号；

多路锁相模块包括N个乘法器和N个滤波器，PGC干涉信号分别与调相波子模块输出的N路信号中的一路信号经过一个乘法器后传入一个滤波器中，得到N路信号的谐波成分；

数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块，N路信号的谐波成分分别经过一个希尔伯特变换子模块后传入一个相位累加子模块，N个相位累加子模块将输出信号传送给降采样平滑滤波子模块，得到最终相位进行输出。

一种具有高精度的光学信号相位解调方法，包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块；

a、调相波子模块生成输出N路信号，包括基频信号cos(ω₀t)、倍频信号cos(2ω₀t)……N阶频信号cos(Nω₀t)，数据采集子模块采集并输出PGC干涉信号：

其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω₀为调制信号频率，为待测量相位信号；

b、将调相波子模块生成输出的N路信号分别与PGC干涉信号进行相乘，再进行滤波操作，得到锁相后N路信号的谐波成分，包括基波成分、二次谐波成分……N次谐波成分；

c、调整调制深度C，得到调整后的N路信号的谐波成分；

d、对调整后的N路信号的谐波成分进行希尔伯特变换操作；

e、求取N路信号被测相位值；

f、对求得的相位值使用降采样平滑的方法得到最终相位输出。

本发明一种具有高精度的光学信号相位解调方法，还可以包括：

1、步骤b中调整后的基波成分为：

二次谐波成分为：

N次谐波成分分别为：

其中，J_k(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(ω₀t)至光源调制器223，同时生成本地载波信号cos(ω₀t)即基频信号103，倍频载波cos(2ω₀t)即倍频信号104，以此规律生成信号cos(Nω₀t)即N阶频信号105，共计N路；

2)数据采集子模块102采集放大器222输出，其形式为：

即PGC干涉信号106，其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω₀为调制信号频率，为待测量相位信号；将公式(1)所示进行贝赛尔展开，可以得到PGC干涉信号106的频谱成分为：

其中J_k(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)，对应输出结果为有如下关系：

$\hat{x}\left(t\right)=\frac{1}{\mathrm{\π}}{\∫}_{-\mathrm{\∞}}^{\mathrm{\∞}}\frac{x\left(\mathrm{\τ}\right)}{t-\mathrm{\τ}}d\mathrm{\τ}=\frac{1}{\mathrm{\π}}{\∫}_{-\mathrm{\∞}}^{\mathrm{\∞}}\frac{x(t-\mathrm{\τ})}{\mathrm{\τ}}d\mathrm{\τ}=x\left(t\right)*\frac{1}{\mathrm{\π}t}---\left(9\right)$

所以在公式(6)(7)(8)输入后即可得到对应的正交分量，使用反正切求解的方法可以分别求得N路信号被测相位值：

其中分别为第一路被测相位信号3111，第二路被测相位信号3112……第N路相位被测信号3113，H为希尔伯特变换过程。最后在输出末端使用降采样平滑的方法得到最终输出结果：

2.所述的前端采集模块10，多路锁相模块11，数据解调模块12，其内部参数设定：

1)调相波子模块101，根据光源调制器223增益倍数输出幅度为1rad-10rad，频率ω₀为20kHz-2MHz的正弦信号进行载波调制；调相波子模块101输出N路本地调相波信号，信号频率为ω₀，2ω₀……Nω₀不等，根据系统采样率取N＝1-10⁴皆可；

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)

其希尔伯特变换的幅度，相位特性为

$H\left(j\mathrm{\&Omega;}\right)=-j\mathrm{sgn}\left(\mathrm{\&Omega;}\right)=\left\{\begin{array}{cc}-j& \mathrm{\&Omega;}>0\\ j& \mathrm{\&Omega;}<0\end{array}\right.---\left(15\right)$

即通过希尔伯特变换为幅频特性为1的全通滤波器。任意信号通过希尔伯特变换器后，其负频率成分作+90°相移，而正频率成分作－90°相移。

基波变换分量311与基波成分117同时经过第一除法器314得到被测相位信号的正切值；二次谐波变换分量312与二次谐波成分118同时通过第二除法器315也可得到被测信号的正切值，依次类推，有N路信号。这几路信号分别由第一反正切子模块317，第二反正切子模块318，第N反正切子模块319求解得到对应的第一路被测相位信号3111与第二路被测相位信号3112，第N路被测相位信号3113；

●多路降噪取平均原理

当系统共有N路信号参与解调时，即有N路相关的信号，但是每路信号分别带有不相关的噪声值，附图4说明N路取平均降噪结果，以下是证明过程：

设任意信号x(t)，其对应的傅立叶变换为：

$F\left(X\right(t\left)\right)={\&Integral;}_{-\mathrm{\&infin;}}^{\mathrm{\&infin;}}x\left(t\right)e^{-j\mathrm{\&omega;}t}dt---\left(16\right)$

当取一定长度信号做傅里叶变换，得到对应频谱的幅度值为：

$F\left(X\right(t\left)\right)={\&Integral;}_0^{T}x\left(t\right)e^{-j\mathrm{\&omega;}t}dt---\left(17\right)$

若N路信号相关，则将N路信号相加，可以得到N路信号累加后的频谱为：

$F\left(X\right(t\left)\right)=N{\&Integral;}_0^{T}x\left(t\right)e^{-j\mathrm{\&omega;}t}dt---\left(18\right)$

此时，求该信号的功率谱为傅里叶变换模值的平方：

$P_N={\left(N{\&Integral;}_0^{T}x\left(t\right)e^{-j\mathrm{\&omega;}t}dt\right)}^2=N^2{P}_0---\left(19\right)$

即N路信号叠加，信号功率为原单一信号功率的N²倍；

但是对于噪声信号而言，N路噪声信号能量不同，根据Parseval定理，信号能量守恒，时域能量等于频域能量，不会因为变换而发生改变，即函数平方的和(或积分)等于其傅里叶转换式平方之和(或者积分)。此时，我们已知噪声为随机信号，服从概率分布，所以噪声功率为：σ²；此时可以计算出其对应的频域功率大小为：

$P_{noise}={\&Integral;}_0^T{\mathrm{\&sigma;}}^{2}dt---\left(20\right)$

当N路随进信号叠加，其噪声功率和为

$P_{Nnoise}=N{\&Integral;}_0^T{\mathrm{\&sigma;}}^{2}dt---\left(21\right)$

所以当信号路数叠加至N路情况下，其信噪比为：

$\frac{P_N}{P_{Nnoise}}=\frac{N^2{P}_0}{{\mathrm{NP}}_{noise}}=N\frac{P_0}{P_{noise}}---\left(22\right)$

由此可以得出结论，在信号路数经过希尔伯特变换被拓展为N路情况下，其信噪比为原来的N倍。

一种提高光学信号相位解调精度的方法，光纤干涉仪输出信号经光电探测器转换进入光学信号相位解调处理系统内，采集模块通过合理调节光源调制器增益倍数，抑制伴生调幅现象的影响；相位计算过程中通过多通道，多频次锁相的办法构造出多路测量信号，每一路信号分别进行希尔伯特变换得到自身的正交分量，各路信号通过与自身正交分量进行反正切查表计算分别得到相位变化量，这些相位变化量呈相关性，且带有一定的共模噪声，通过平均平滑的方法抑制共模噪声，使输出获得更好的信噪比；本方法有效提高光学信号相位解调精度，可广泛用于高精度光纤测量和光纤传感等领域。

一种提高光学信号相位解调精度的方法，包含前端采集模块10，其特征是多路锁相模块11，数据解调模块12。

1)调相波子模块101输出调相波信号cos(ω₀t)至光源调制器223，同时生成本地载波信号cos(ω₀t)即基频信号103，倍频载波cos(2ω₀t)即倍频信号104，以此规律生成信号cos(Nω₀t)即N阶频信号105，共计N路；

2)数据采集子模块102采集放大器222输出，其形式为：

即PGC干涉信号106，其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω₀为调制信号频率，为待测量相位信号；将公式(23)所示进行贝赛尔展开，可以得到PGC干涉信号106的频谱成分为：

其中J_k(C)为贝塞尔函数系数，k为信号高阶分量。

3)将基频信号cos(ω₀t)103，倍频信号cos(2ω₀t)104……N阶频信号cos(Nω₀t)105分别与PGC干涉信号106进行相乘运算，之后进行滤波运算可以得到锁相后的基波成分117，二次谐波成分118……N次谐波成分119，其形式为：

其中A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。

4)调整调制深度C，使J₁(C)≈J₃(C)、J₀(C)≈J₂(C)，J_p(C)≈J_q(C)就可以得到修正后的基波成分117，二次谐波成分118……N次谐波成分119，再经过高通滤波处理即消除伴生调幅系数m对解调结果的影响，得到

5)对调整后的基波成分117即公式(28)，二次谐波成分118即公式(29)……N次谐波成分119即公式(30)进行希尔伯特变换操作，设任意需要进行变换的信号为x(t)，对应输出结果为有如下关系：

$\hat{x}\left(t\right)=\frac{1}{\mathrm{\&pi;}}{\&Integral;}_{-\mathrm{\&infin;}}^{\mathrm{\&infin;}}\frac{x\left(\mathrm{\&tau;}\right)}{t-\mathrm{\&tau;}}d\mathrm{\&tau;}=\frac{1}{\mathrm{\&pi;}}{\&Integral;}_{-\mathrm{\&infin;}}^{\mathrm{\&infin;}}\frac{x(t-\mathrm{\&tau;})}{\mathrm{\&tau;}}d\mathrm{\&tau;}=x\left(t\right)*\frac{1}{\mathrm{\&pi;}t}---\left(31\right)$

所以在公式(28)(29)(30)输入后即可得到对应的正交分量，使用反正切求解的方法可以分别求得N路信号被测相位值：

其中分别为第一路被测相位信号3111，第二路被测相位信号3112……第N路相位被测信号3113，H为希尔伯特变换过程。最后在输出末端使用降采样平滑的方法得到最终输出结果：

前端采集模块10，多路锁相模块11，数据解调模块12，其内部参数是：

1)调相波子模块101，根据光源调制器223增益倍数输出幅度为1rad-10rad，频率ω₀为20kHz-2MHz的正弦信号进行载波调制；调相波子模块101输出N路本地调相波信号，信号频率为ω₀，2ω₀……Nω₀不等，根据系统采样率取N＝1-10⁴皆可；

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～10⁵rad；

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(ω₀t)、倍频信号cos(2ω₀t)……N阶频信号cos(Nω₀t)，数据采集子模块采集并输出PGC干涉信号；

多路锁相模块包括N个乘法器和N个滤波器，PGC干涉信号分别与调相波子模块输出的N路信号中的一路信号经过一个乘法器后传入一个滤波器中，得到N路信号的谐波成分；

数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块，N路信号的谐波成分分别经过一个希尔伯特变换子模块后传入一个相位累加子模块，N个相位累加子模块将输出信号传送给降采样平滑滤波子模块，得到最终相位进行输出。

2.一种具有高精度的光学信号相位解调方法，其特征在于：包括前端采集模块、多路锁相模块和数据解调模块，前端采集模块包括数据采集子模块和调相波子模块，多路锁相模块包括N个乘法器和N个滤波器，数据解调模块包括N个希尔伯特变换子模块、N个相位累加子模块和降采样平滑滤波子模块；

a、调相波子模块生成输出N路信号，包括基频信号cos(ω₀t)、倍频信号cos(2ω₀t)……N阶频信号cos(Nω₀t)，数据采集子模块采集并输出PGC干涉信号：

其中P为信号幅度，A为光源直流项，B为光源交流项，C为调制深度，ω₀为调制信号频率，为待测量相位信号；

b、将调相波子模块生成输出的N路信号分别与PGC干涉信号进行相乘，再进行滤波操作，得到锁相后N路信号的谐波成分，包括基波成分、二次谐波成分……N次谐波成分；

c、调整调制深度C，得到调整后的N路信号的谐波成分；

d、对调整后的N路信号的谐波成分进行希尔伯特变换操作；

e、求取N路信号被测相位值；

f、对求得的相位值使用降采样平滑的方法得到最终相位输出。

3.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤b中调整后的基波成分为：

二次谐波成分为：

N次谐波成分分别为：

其中，J_k(C)为贝塞尔函数系数，k为信号高阶分量，A为光强直流分量，B为交流幅度，p，q代表不同阶数；m为伴生调幅系数。

4.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤c中调整后的基波成分为：

二次谐波成分为：

N次谐波成分为：

5.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤e中N路信号被测相位值的第一路被测相位信号为：

第二路被测相位信号为：

第N路相位被测信号为：

6.根据权利要求2所述的一种具有高精度的光学信号相位解调方法，其特征在于：所述的步骤f中，最终相位为：