CN101713685B - 一种无源零差正交解调技术的相位补偿方法 - Google Patents
一种无源零差正交解调技术的相位补偿方法 Download PDFInfo
- Publication number
- CN101713685B CN101713685B CN2009102240600A CN200910224060A CN101713685B CN 101713685 B CN101713685 B CN 101713685B CN 2009102240600 A CN2009102240600 A CN 2009102240600A CN 200910224060 A CN200910224060 A CN 200910224060A CN 101713685 B CN101713685 B CN 101713685B
- Authority
- CN
- China
- Prior art keywords
- interference signal
- signal
- way
- quadrature
- sin
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
- Measuring Phase Differences (AREA)
Abstract
本发明涉及光纤传感技术领域,具体涉及一种无源零差正交解调的相位补偿方法。该方法先将这两路干涉信号通过去除直流信号、归一化处理后,再相乘、低通滤波得到直流信号,对获得的直流信号进行反正弦函数运算,从而获得两路干涉信号偏离正交的角度;利用此偏离角度得到两路正交的干涉信号,再利用正交信号的相位解调技术,就可以解调出调制信号。本发明直接对两路原始的干涉信号进行处理和计算,能够有效地避免两路干涉信号由于偏离正交所造成的调制信号的失真。
Description
技术领域
本发明涉及光纤传感技术领域,具体涉及一种无源零差正交解调技术的相位补偿方法。
背景技术
干涉型光纤传感器按照是否需要主动控制可以分为有源型和无源型。无源型光纤干涉仪由于不需要主动控制,可以达到的测量频率很高,测量频率只受电子线路工作频率的限制,因此广泛应用到动态信号(如振动、声)的测量中。
无源零差正交解调技术广泛应用于干涉型光纤传感器的相位调制恢复过程,在测量动态信号时具有测量频率高,相位分辨率高,测量范围大等特点。该技术的关键所在是获得两路正交的干涉信号。这两路干涉信号可以是两个不同光程的干涉仪产生的干涉信号,如:[K.A.Murphy,M.F.Gunther,A.M.Vengsarkar,and R.O.Claus,“Quadrature phase-shifted,extrinsic Fabry-Perot optical fiber sensors”,Optics Letters,1991,16(4),273-275];还可以是两束不同波长的光通过同一光纤干涉仪的干涉信号,如:[O.B.Wright,“Stabilizeddual-wavelength fiber-optic interferometer for vibration measurement”,Optics Letter,1991,16,56-58];也可以利用基于3×3耦合器的干涉仪产生的三路干涉信号(如g1,g2和g3)得到两路正交的干涉信号,如用g2-g3和g2+g3,如:[S.K.Sheem,“Optical fiber interferometers with3×3 directional coupler:Analysis”,Journal of Appllied Physics.,1981,52(6),3865-3872]或者g2-g3和2g1-g2-g3[R.G.Priest,“Analysis of fiberinterferometer utilizing 3×3 fiber coupler”,IEEE Transactions onMicrowave Theory and Techniques,1982,30(10),1589-1591]。这些方法都可以得到两路正交的干涉信号。在获得两路正交信号后,再利用反正切算法或交叉相乘法,就可以解调出干涉信号中的相位调制。
在无源零差正交解调技术中,要正确地解调出调制信号,就必须要求两路干涉信号完全正交。但是,在实际应用中,由于加工工艺所限,以及工作环境的影响,两路干涉信号常常并不完全正交,而是偏离正交状态一定的相位角,这样会引起测量误差,并造成解调出的信号失真。已经有其它研究人员对两路干涉信号的相位差偏离90°对正交解调产生的影响做了研究,例如Y.L.Lo等人通过测量两路干涉信号的李萨如图形与坐标轴的交点位置,来测量两路干涉信号的相位差偏离90°的相位误差,并进一步对其进行补偿[Y.L.Lo,J.S.Sirkis,andC.C.Chang,“Passive signal processing of in-line fiber etalon sensor forhigh strain-rate loading”,IEEE Journal of Lightwave Technology,15(8),1578-1586]。但是这种方法需要人工操作,且易受李萨如图形线宽影响,准确性不高,难以在实际中应用。本发明直接对原始的两路干涉信号进行处理和计算,获得的两路干涉信号偏离90°的相位差,再产生一路补偿掉该相位差的新的干涉信号,与其中一路原始信号形成两路完全正交的两路干涉信号,最终能正确地解调出调制信号。该方法简单易行,便于应用,测量精度高。
发明内容
本发明的目的是克服现有技术中的不足,提出一种无源零差正交解调技术的相位补偿方法,应用于干涉型光纤传感器系统中;该方法对两路原始的干涉信号的正交性要求降低,抗干扰能力强,测量精度高。
本发明的目的是通过下述技术方案实现的。
本发明的一种无源零差正交解调技术的相位补偿方法,具体实现步骤如下:
(1)定义两路不完全正交的干涉信号I1和I2为I1=B1+C1cos(α+Δφ)和I2=B2+C2sinα,其中B1和B2分别是两路干涉信号的背景光,C1和C2分别是两路干涉信号的对比度,Δφ为两路干涉信号偏离正交的相位差,α为需要测量的调制信号。
(2)将上述两路信号分别经过信号的直流滤除来去除直流分量B1和B2,然后对信号的振幅进行归一化处理之后,两路干涉信号表示为I1′=cos(α+Δφ)和I2′=sinα。
(3)将步骤(2)得到的不完全正交的这两路干涉信号相乘可得
上式中等号右边第一项正弦信号的频率是干涉信号的频率的两倍,是交流信号,第二项中Δφ为常量,所以第二项为直流分量;滤出I1′×I2′中的直流分量得到的值,再对该直流分量乘以-2计算反正弦函数得到Δφ的值;而与I2′相正交的新生成的干涉信号可以表示为
I1 new=cosα=[cos(α+Δφ)+sinαsinΔφ]/cosΔφ,
化简后为
I1 new=[I1′+I2′sinΔφ]/cosΔφ。
最后得到新生成的两路正交的干涉信号是I1 new和I2’。再利用已有的正交信号解调方法(如反正切法,或交叉相乘法等)就可以解调出更加准确地调制信号α。
有益效果
本发明所公开的方法经过实验的验证,简单实用,既不需改变硬件设置又无需人工干预,提高了测量精度,效果显著。
附图说明
图1为未校正前信号数据
其中,a-实际测量获得的原始两路不正交干涉信号、b-为这两路干涉信号的李萨如图、c-为应用不正交的原始的两路干涉信号直接解调出的波形图;
图2为应用本方法所得出的两路干涉信号偏移角度连续测量图;
图3为应用本方法校正干涉信号后所得到的信号数据;
其中,a-应用本方法校正后的新生成的两路正交的干涉信号,b-为新生成的两路正交的干涉信号的李萨如图,c-为新生成的两路正交的干涉信号解调出的波形图;
图4应用本方法前后解调出的信号峰峰值随相位差偏移的关系
其中,a-原始的两路干涉信号解调出的信号的峰峰值随相位差偏移的关系,b-应用本方法后解调出的信号峰峰值随相位差偏移的关系。
具体实施方式
下面将结合附图和实施例对本发明做详细说明。
将输入的两路不完全正交的干涉信号I1和I2表示为
I1=B1+C1cos(α+Δφ) (1)
I2=B2+C2sinα (2)
其中B1和B2分别是两路干涉信号的背景光,C1和C2分别是两路干涉信号的对比度,Δφ为两路干涉信号偏离正交的相位差,α为需要测量的调制信号。经过去直流和归一化处理后,原始的干涉信号表示为
I1′=cos(α+Δφ) (3)
I2′=sinα (4)
将这两路干涉信号相乘,利用积化和差公式可得
式(5)中等号右边第一项正弦信号的频率是干涉信号频率的两倍,是交流信号,第二项中Δφ为常量,所以第二项为直流分量,滤出直流分量就可以得到的值,再将该直流分量乘以-2,计算反正弦函数就可以得到Δφ的值;而与I2相正交的新生成的的干涉信号I1 new即:
I1 new=cosα (6)
由三角函数关系cos(α+Δφ)=cosαcosΔφ-sinαsinΔφ有
I1 new=cosα=[cos(α+Δφ)+sinαsinΔφ]/cosΔφ (7)
化简得
I1 new=[I1′+I2′sinΔφ]/cosΔφ (8)
得到新生成的两路正交的干涉信号I1 new和I2’。再利用已有的正交信号解调方法就可以解调出调制信号α。
下面是该方法的实用例子:
图1(a)中为实际测量获得的原始两路不正交干涉信号;图1((b)为这两路干涉信号的李萨如图,从李萨如图中可以很明显地看出原始的两路干涉信号并不正交,并且相位差偏离90°较大;图1(c)为应用不正交的原始的两路干涉信号直接解调出的波形,解调出的信号已经失真。可见原始的两路干涉信号是否正交对解调所得到信号的质量影响很大。
接下来用本发明的方法处理原始的干涉信号。首先应用本发明的方法来测量两路干涉信号偏离90°的相位差,图2是连续测量的结果,从图中可以看出偏离正交的相位差约为52°。利用得到的这个相位差将原始的两路不正交的干涉信号转换为新的两路正交的干涉信号,如图3(a)所示。图3(b)为新产生的两路正交的干涉信号的李萨如图,可以看出,这是一个正圆,表明这两路新生成的干涉信号的相位差为90°。用新生成的两路正交的干涉信号进行正交解调,得到的解调后的信号为图3(c),与原来的解调结果图1(c)相比较,可以看出,应用本发明来解调出的信号没有失真。
另外,利用计算机模拟了当调制信号的峰峰值为8弧度时,两路干涉信号偏离正交的角度对解调出来的信号的影响。图4(a)是应用原始的两路干涉信号解调后的峰峰值随相位差偏移的关系,随着两路干涉信号偏离正交角度的增大,解调出来信号的峰峰值与输入信号峰峰值的差距也在增大。用本发明所述方法处理后,虽然两路干涉信号偏离正交的角度在改变,但解调信号的峰峰值仍然稳定在8弧度周围,最大误差为0.006弧度,误差极小,如图4(b)所示。
本发明所公开的方法经过实验的验证,简单实用,既不需改变硬件设置又无需人工干预,提高了测量精度,效果显著。
Claims (1)
1.一种无源零差正交解调技术的相位补偿方法,其特征在于具体实现步骤如下:
步骤(1)定义两路不完全正交的干涉信号I1和I2为I1=B1+C1cos(α+Δφ)和I2=B2+C2sinα,其中B1和B2分别是两路干涉信号的背景光,C1和C2分别是两路干涉信号的对比度,Δφ为两路干涉信号偏离正交的相位差,α为需要测量的调制信号;
步骤(2)将上述两路信号分别经过信号的直流滤除来去除直流分量B1和B2,然后对信号的振幅进行归一化处理之后,两路干涉信号表示为I1′=cos(α+Δφ)和I2′=sinα;
步骤(3)将步骤(2)得到的不完全正交的这两路干涉信号相乘可得
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102240600A CN101713685B (zh) | 2009-12-03 | 2009-12-03 | 一种无源零差正交解调技术的相位补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102240600A CN101713685B (zh) | 2009-12-03 | 2009-12-03 | 一种无源零差正交解调技术的相位补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101713685A CN101713685A (zh) | 2010-05-26 |
CN101713685B true CN101713685B (zh) | 2011-07-27 |
Family
ID=42417510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102240600A Expired - Fee Related CN101713685B (zh) | 2009-12-03 | 2009-12-03 | 一种无源零差正交解调技术的相位补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101713685B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630164B (zh) * | 2013-12-23 | 2016-11-16 | 中国计量科学研究院 | 双路干涉信号正交标准化方法及装置 |
CN107421435B (zh) * | 2017-03-24 | 2019-11-12 | 浙江大学 | 基于软件移相的零差正交激光干涉信号生成方法 |
CN109660305A (zh) * | 2018-11-16 | 2019-04-19 | 上海无线电设备研究所 | 一种对多通道接收机幅度、相位校正及补偿方法 |
CN111006611B (zh) * | 2019-12-20 | 2021-10-22 | 哈尔滨工业大学(深圳) | 一种基于两步非线性移相的快速高精度相位提取方法 |
CN112344974B (zh) * | 2020-10-21 | 2022-05-20 | 北京自动化控制设备研究所 | 用于干涉型光纤传感器的相位解调系统及解调方法 |
CN113405578A (zh) * | 2021-05-25 | 2021-09-17 | 天津大学 | 基于偏振干涉和dcm算法的高稳定性动态相位解调补偿方法 |
CN114858062B (zh) * | 2022-06-13 | 2022-09-16 | 探维科技(北京)有限公司 | 探测信号处理方法、装置、介质、设备及光栅尺 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6556509B1 (en) * | 1999-11-08 | 2003-04-29 | Optiphase, Inc. | Demodulator and method for interferometric outputs of increased accuracy |
CN1752836A (zh) * | 2005-09-27 | 2006-03-29 | 北京理工大学 | 一种光纤干涉信号解调方法 |
-
2009
- 2009-12-03 CN CN2009102240600A patent/CN101713685B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6556509B1 (en) * | 1999-11-08 | 2003-04-29 | Optiphase, Inc. | Demodulator and method for interferometric outputs of increased accuracy |
CN1752836A (zh) * | 2005-09-27 | 2006-03-29 | 北京理工大学 | 一种光纤干涉信号解调方法 |
Non-Patent Citations (2)
Title |
---|
葛云生.干涉型光纤传感器的正交解调中调制初相位影响及校正算法研究.《仪器仪表学报》.2009,第30卷(第10期),229-232. * |
葛云生.干涉型光纤传感器的正交解调算法分析.《仪器仪表学报》.2008,第29卷(第8期),121-123. * |
Also Published As
Publication number | Publication date |
---|---|
CN101713685A (zh) | 2010-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101713685B (zh) | 一种无源零差正交解调技术的相位补偿方法 | |
US10641600B2 (en) | Real-time normalization apparatus and method of phase generated carrier demodulation in sinusoidal phase modulation interferometer | |
CN109459070B (zh) | 一种pgc相位解调法中相位延迟提取与补偿方法 | |
US20030030816A1 (en) | Nonlinearity error correcting method and phase angle measuring method for displacement measurement in two-freqency laser interferometer and displacement measurement system using the same | |
EP1985967B1 (en) | Fiber interferometric sensor and phase compensation method of PGC demodulator | |
US10162245B2 (en) | Distributed acoustic sensing system based on delayed optical hybrid phase demodulator | |
Yu et al. | High stability and low harmonic distortion PGC demodulation technique for interferometric optical fiber sensors | |
US20070041020A1 (en) | Method and apparatus for direct detection of signals from a differential delay heterodyne interferometric system | |
US9658052B2 (en) | Method for reducing interference from scattered light/reflected light of interference path by generating carrier through phase | |
US11469816B2 (en) | Distributed sensing apparatus | |
CN101718563A (zh) | 基于3×3光纤耦合器的相移白光干涉测量方法 | |
CN110285843B (zh) | 一种大振幅信号分布式弱光栅阵列传感系统及解调方法 | |
CN102997945A (zh) | 光纤分布式扰动传感器的多点扰动定位方法 | |
CN106323478B (zh) | 抗偏振衰落的光纤干涉型传感器相位生成载波调制解调系统 | |
Yan et al. | Phase-modulated dual-homodyne interferometer without periodic nonlinearity | |
US20010028679A1 (en) | Phase demodulator, phase difference detector, and interferometric system using the phase difference detector | |
US6600586B1 (en) | Normalization method for acquiring interferometer phase shift from frequency division multiplexed fiber optic sensor arrays | |
CN103983209A (zh) | 用于光纤干涉条纹投射三维形貌测量的条纹相位稳定方法 | |
CN110118594B (zh) | 一种基于偏振分极接收的光学相位解调方法和系统 | |
US11885607B2 (en) | Device for interferometric distance measurement | |
CN108007340A (zh) | 相位生成载波反正切解调中非线性误差的实时计算方法 | |
US20230296367A1 (en) | Distortion-corrected phase generated carrier demodulation method using multitone mixing | |
Jiang et al. | Note: Phase compensation in the fiber optical quadrature passive demodulation scheme | |
Chawah et al. | Real time and adaptive Kalman filter for joint nanometric displacement estimation, parameters tracking and drift correction of EFFPI sensor systems | |
JP3108866B2 (ja) | レーザ振動変位測定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110727 Termination date: 20131203 |