CN103616020B - 基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 - Google Patents
基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 Download PDFInfo
- Publication number
- CN103616020B CN103616020B CN201310610990.6A CN201310610990A CN103616020B CN 103616020 B CN103616020 B CN 103616020B CN 201310610990 A CN201310610990 A CN 201310610990A CN 103616020 B CN103616020 B CN 103616020B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- fiber optic
- optic loop
- sine wave
- eigenfrequency
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C19/00—Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
- G01C19/58—Turn-sensitive devices without moving masses
- G01C19/64—Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams
- G01C19/72—Gyrometers using the Sagnac effect, i.e. rotation-induced shifts between counter-rotating electromagnetic beams with counter-rotating light beams in a passive ring, e.g. fibre laser gyrometers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optics & Photonics (AREA)
- Electromagnetism (AREA)
- Power Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Manufacturing & Machinery (AREA)
- Gyroscopes (AREA)
Abstract
本发明公开了一种基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法,该装置及方法采用正弦波施加调制,通过检测光纤环干涉输出信号中的一次谐波幅度来测得本征频率,其技术效果在于使用正弦波而不是方波施加调制,克服了方波调制带来的多种缺陷,并且检测的是输出信号一次谐波的幅度,可采用相关检测技术获得高的信噪比。
Description
技术领域
本发明涉及一种光纤环本征频率的测量装置及方法,特别涉及一种基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法。
背景技术
光纤环是光纤陀螺中敏感角速度的关键器件,进入光纤环相反方向传播的两光波沿光纤环传播一周后干涉产生Sagnac效应,通过Sagnac效应实现对角速度的测量。光波在光纤环中传播时间τ的倒数的二分之一为光纤环的本征频率本征频率是光纤环的重要参数之一,特别是对于闭环光纤陀螺,其工作时调制频率就应等于光纤环的本征频率,当调制频率偏离本征频率时,会产生调制误差,并影响光纤陀螺零偏。因此精确测量光纤环的本征频率对提高光纤陀螺的性能意义重大。
测量光纤环本征频率多采用的是基于方波调制的方法。由于实际应用中的调制方波都有一定的上升沿与下降沿,这将使输出信号必然存在无效调制区,即调制频率即使等于本征频率,输出信号依然存在一定尖峰脉冲,从而影响本征频率的测量精度。
也有人对基于方波调制的光纤环本征频率测量方法进行了改进,诸如采用1/2本征频率的方波、不对称方波等。这些方法在一定程度上能够提高本征频率的精度,但由于需要测量输出方波的占空比或脉冲宽度,对采集设备的要求较高,致使测量成本较高。
发明内容
为了解决目前光纤环本征频率测量精度不足及测量成本高的问题,本发明提出一种能以较低成本及较高精度测得光纤环本征频率的基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法。
为了实现上述技术目的,本发明的技术方案是,一种基于正弦波调制及一次谐波检测的光纤环本征频率测量装置,包括光源、耦合器、调制器、光纤环、D/A及滤波放大电路、波形调制输出及一次谐波检测电路、放大采样电路和探测器,所述的光源连接耦合器,所述的耦合器分别连接至探测器和调制器,光纤环连接调制器,探测器串联放大采样电路后连接至波形调制输出及一次谐波检测电路,波形调制输出及一次谐波检测电路串联D/A及滤波放大电路后连接至调制器,其特征在于,所述的波形调制输出及一次谐波检测电路为正弦波调制输出及一次谐波检测电路。
所述的装置,所述的调制器为Y波导,D/A及滤波放大电路输出的正弦信号直接施加于Y波导上。
所述的装置,所述的调制器包括压电换能器相位调制器、调制耦合器和偏振器,所述的压电换能器相位调制器设置于光纤环的一个输出端上,所述的调制耦合器分别连接偏振器、压电换能器相位调制器和光纤环的另一个输出端,D/A及滤波放大电路输出的正弦信号施加于压电换能器相位调制器上,所述的偏振器设置于耦合器与调制耦合器之间。
一种基于正弦波调制及一次谐波检测的光纤环本征频率测量方法,采用所述的装置,步骤包括:
第一步:将光纤环转动或摇摆,使光纤环的Sagnac相位差值与π/2之间差值小于0.2,由光纤环的长度L0计算出光在光纤环中传播时间的估计值τ0,
其中n为光纤的折射率,c为真空中的光速,进而得到光纤环本征频率的估计值fe0,
则调制频率初值fm0=2kfe0,其中k取使调制频率不超过调制器可用频率及波形调制输出及一次谐波检测电路、放大采样电路和探测器的可检测频率上限前提下的最大正整数,正弦波调制输出及一次谐波检测电路输出调制频率初值fm0的数字正弦波;
第二步:数字正弦波经D/A及滤波放大电路后,在调制器上施加正弦波调制;
第三步:探测器接收经正弦波调制后的干涉输出,经放大采样电路后,信号输出至正弦波调制输出及一次谐波检测电路,由正弦波调制输出及一次谐波检测电路检测出信号中的一次谐波幅度;
第四步:在大于(2k-1)fe0而小于(2k+1)fe0的范围内改变调制频率,重复上述的第二步、第三步,测得10至20个不同调制频率下的干涉输出中的一次谐波幅度;
第五步:在本征频率的2k倍两侧将调制频率与一次谐波幅度做线性拟合,两直线与横轴交点的中点认为是拟合曲线中一次谐波幅度最小的调制频率fm,频率fm为本征频率fe的2k倍,即fe=fm/2k,从而得到光纤环的本征频率fe。
本发明的技术效果在于:
(1)使用正弦波而不是方波施加调制,克服了方波调制带来的多种缺陷,即不存在方波调制中上升沿与下降沿所导致的无效调制区问题,且无需方波调制中的精确测量占空比或脉冲宽度的设备,即在保证测量精度的前提下并没有增加测量成本。并且检测的是输出信号一次谐波的幅度,可采用相关检测技术获得高的信噪比。
(2)由于使用正弦波施加调制,调制器除可使用Y波导外,还可使用PZT施加调制,从而降低成本。
(3)工作时调制频率在本征频率的偶数倍附近(即fm=2kfe),测量误差与k成反比,即测量时可根据调制器可用频率及检测设备的可检测频率上限选择最大的k,以获得最小的测量误差,从而提高测量精度。
(4)测量误差与正弦波调制信号幅度成反比,即测量时可使用调制器能承受的最大调制幅度,以获得更小的测量误差,从而提高测量精度。而对于方波调制,调制信号产生的光相位变化超过π后再增大调制信号幅度无实际意义,而且幅度越大方波上升沿与下降沿产生的无效调制的影响越大。
下面结合附图对本发明作进一步说明。
附图说明
图1为本发明的测量方法所使用系统的结构示意图;
图2为使用Y波导的调制结构示意图;
图3为使用压电换能器相位调制器的调制结构示意图;
图4为调制频率与一次谐波幅度关系曲线。
其中1为光源,2为探测器,3为耦合器,4为调制器,41为Y波导,42为偏振器,43为调制耦合器,44为压电换能器相位调制器,5为光纤环,6为放大采样电路,7为正弦波调制输出及一次谐波检测电路,8为D/A及滤波放大电路。
具体实施方式
参见图1、图2、图3,本发明装置包括光源、耦合器、调制器、光纤环、D/A及滤波放大电路、波形调制输出及一次谐波检测电路、放大采样电路和探测器,光源连接耦合器,耦合器分别连接至探测器和调制器,光纤环连接调制器,探测器串联放大采样电路后连接至波形调制输出及一次谐波检测电路,波形调制输出及一次谐波检测电路串联D/A及滤波放大电路后连接至调制器,波形调制输出及一次谐波检测电路为正弦波调制输出及一次谐波检测电路,调制器为Y波导或为PZT、调制耦合器和偏振器,PZT设置于光纤环的一个输出端上,调制耦合器分别连接偏振器、压电换能器相位调制器和光纤环的另一个输出端,偏振器设置于耦合器与调制耦合器之间。调制器如采用Y波导,则可采用较高的调制频率(即可选用较大的k),从而获得较高的测量精度,而调制器采用PZT则可降低系统成本,应用中可根据实际情况选择Y波导或PZT。
本发明测量方法的测量步骤为:
(1)测量系统放置于转台或摇摆台上,转台或摇摆台以一定角速度转动或摇摆,使光纤环的Sagnac相位差值与π/2之间差值小于0.2。由光纤环本征频率的估计值fe0及由调制器可用频率及检测设备的可检测频率上限选取出的k确定调制频率初值fm0=2kfe0,正弦波调制输出及一次谐波检测电路7输出该频率的数字正弦波。
(2)数字正弦波经D/A及滤波放大电路8后,在调制器4上施加正弦波调制,正弦波的幅度可使用调制器能承受的最大调制幅度,以此降低测量误差。
(3)探测器2接收干涉输出,经放大采样电路6后,信号输出至正弦波调制输出及一次谐波检测电路7,由正弦波调制输出及一次谐波检测电路7检测出信号中的一次谐波幅度,一次谐波幅度的检测采用相关检测技术以提高信噪比。
(4)在一定范围内改变调制频率,重复上述的(2)~(3)步骤,测得一系列调制频率下的干涉输出中的一次谐波幅度。
(5)在本征频率的2k倍两侧将调制频率与一次谐波幅度做线性拟合,即拟合曲线V1=af+b,其中V2为一次谐波幅度,f为调制频率,a、b为拟合参数,两直线与横轴交点的中点认为是拟合曲线中一次谐波幅度最小的调制频率,该频率fm为本征频率的2k倍,即fe=fm/2k,从而测得光纤环的本征频率。
本发明的原理为:
设调制器上所加的正弦波信号产生的光相位变化为:
式中为光相位变化,ωm=2πfm为调制角频率,fm为调制频率,为调制幅度,t为时间,由于光纤环中相反方向传播的两束光通过调制器的时间不同,其时间差恰好为光波在光纤环中的传播时间τ,则可得调制产生的正反传播光的相位差为:
式中fe为本征频率,令称为正弦波调制时的有效调制深度,干涉输出信号为:
上式中为Sagnac相位差。上式可展开为:
其中Jn为第一类n阶贝塞尔函数。则可得输出信号的一次谐波分量为:
由上式可得一次谐波的幅度为由于系统放置于转动的转台上,与π/2之间差值小于0.2,不妨设当调制深度接近0时,J1(φ0)≈φ0/2,即此时一次谐波的幅度为V1=V0|φ|。当调制频率在本征频率的偶数倍附近时,令fm=2kfe(1+εfe),2k为倍频数,εfe为本征频率测量的相对误差,则可得:
即一次谐波幅度是fm/2k的V形折线函数,如图4所示。
设系统可检测到的一次谐波最小幅度为ΔV,即测量时2V0|J1(φ0)|≥△V,令εV=ΔV/2V0,称其为谐波相对探测下限,由此得测量误差:
可见测量误差与调制信号幅度成反比,测量时可使用调制器能承受的最大调制幅度,以获得更小的测量误差。测量误差与k成反比,即测量时可根据检测设备的可检测频率上限选择最大的k,以获得最小的测量误差。由式(6)知一次谐波幅度随调制频率的变化是V形折线函数,V形折线的顶点对应本征频率的偶数倍,因此可通过测量不同调制频率下的一次谐波幅度,然后通过在V形折线顶点两侧分别拟合直线来获得V形折线顶点的方式获得本征频率,这样可以使误差进一步降低。
因此本发明可以较低成本及较高精度测得光纤环本征频率。
Claims (1)
1.一种基于正弦波调制及一次谐波检测的光纤环本征频率测量方法,采用一种基于正弦波调制及一次谐波检测的光纤环本征频率测量装置,包括光源、耦合器、调制器、光纤环、D/A及滤波放大电路、波形调制输出及一次谐波检测电路、放大采样电路和探测器,所述的光源连接耦合器,所述的耦合器分别连接至探测器和调制器,光纤环连接调制器,探测器串联放大采样电路后连接至波形调制输出及一次谐波检测电路,波形调制输出及一次谐波检测电路串联D/A及滤波放大电路后连接至调制器,所述的波形调制输出及一次谐波检测电路为正弦波调制输出及一次谐波检测电路,其特征在于,步骤包括:
第一步:将光纤环转动或摇摆,使光纤环的Sagnac相位差值与π/2之间差值小于0.2,由光纤环的长度L0计算出光在光纤环中传播时间的估计值τ0,
其中n为光纤的折射率,c为真空中的光速,进而得到光纤环本征频率的估计值fe0,
则调制频率初值fm0=2kfe0,其中k取使调制频率不超过调制器可用频率及波形调制输出及一次谐波检测电路、放大采样电路和探测器的可检测频率上限前提下的最大正整数,正弦波调制输出及一次谐波检测电路输出调制频率初值fm0的数字正弦波;
第二步:数字正弦波经D/A及滤波放大电路后,在调制器上施加正弦波调制;
第三步:探测器接收经正弦波调制后的干涉输出,经放大采样电路后,信号输出至正弦波调制输出及一次谐波检测电路,由正弦波调制输出及一次谐波检测电路检测出信号中的一次谐波幅度;
第四步:在大于(2k-1)fe0而小于(2k+1)fe0的范围内改变调制频率,重复上述的第二步、第三步,测得10至20个不同调制频率下的干涉输出中的一次谐波幅度;
第五步:在本征频率的2k倍两侧将调制频率与一次谐波幅度做线性拟合,两直线与横轴交点的中点认为是拟合曲线中一次谐波幅度最小的调制频率fm,频率fm为本征频率fe的2k倍,即fe=fm/2k,从而得到光纤环的本征频率fe。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310610990.6A CN103616020B (zh) | 2013-11-26 | 2013-11-26 | 基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310610990.6A CN103616020B (zh) | 2013-11-26 | 2013-11-26 | 基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103616020A CN103616020A (zh) | 2014-03-05 |
CN103616020B true CN103616020B (zh) | 2016-01-06 |
Family
ID=50166727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310610990.6A Expired - Fee Related CN103616020B (zh) | 2013-11-26 | 2013-11-26 | 基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103616020B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973219B (zh) * | 2016-04-29 | 2018-08-17 | 重庆华渝电气集团有限公司 | 光纤陀螺调制解调方法 |
CN107389097B (zh) * | 2017-07-25 | 2020-03-10 | 北京航空航天大学 | 光纤陀螺Sagnac光纤环本征频率跟踪测量方法 |
CN107356266B (zh) * | 2017-07-25 | 2020-03-10 | 北京航空航天大学 | 基于偶倍本征频率锯齿波调制的光纤陀螺本征频率测量方法 |
CN112833874A (zh) * | 2021-01-08 | 2021-05-25 | 中国船舶重工集团公司第七0七研究所 | 一种适用于光纤陀螺的新型四态调制解调方法 |
CN114264318B (zh) * | 2021-12-06 | 2024-06-11 | 河北汉光重工有限责任公司 | 一种闭环光纤陀螺固有频率测试方法和装置 |
CN116045957B (zh) * | 2023-03-31 | 2023-06-16 | 中国船舶集团有限公司第七〇七研究所 | 一种基于光纤陀螺扩频采样的误差消除方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101482446A (zh) * | 2009-02-24 | 2009-07-15 | 北京航天时代光电科技有限公司 | 光纤陀螺用y波导相位调制线性度测试方法 |
CN101696882A (zh) * | 2009-10-26 | 2010-04-21 | 浙江大学 | 光纤陀螺仪开环信号调制解调电路 |
CN102183249A (zh) * | 2011-03-15 | 2011-09-14 | 北京大学 | 一种光纤陀螺仪的萨格奈克相移跟踪方法 |
CN102435186A (zh) * | 2011-11-28 | 2012-05-02 | 北京大学 | 一种光纤陀螺的数字信号处理方法、装置及光纤陀螺仪 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7817284B2 (en) * | 2007-08-08 | 2010-10-19 | The Charles Stark Draper Laboratory, Inc. | Interferometric fiber optic gyroscope with off-frequency modulation signals |
-
2013
- 2013-11-26 CN CN201310610990.6A patent/CN103616020B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101482446A (zh) * | 2009-02-24 | 2009-07-15 | 北京航天时代光电科技有限公司 | 光纤陀螺用y波导相位调制线性度测试方法 |
CN101696882A (zh) * | 2009-10-26 | 2010-04-21 | 浙江大学 | 光纤陀螺仪开环信号调制解调电路 |
CN102183249A (zh) * | 2011-03-15 | 2011-09-14 | 北京大学 | 一种光纤陀螺仪的萨格奈克相移跟踪方法 |
CN102435186A (zh) * | 2011-11-28 | 2012-05-02 | 北京大学 | 一种光纤陀螺的数字信号处理方法、装置及光纤陀螺仪 |
Non-Patent Citations (3)
Title |
---|
一种低成本的光纤陀螺信号处理方法;邢怀飞等;《光学技术》;20050930;第31卷(第5期);790-792 * |
全数字开环光纤陀螺工作点实时控制;宋章启等;《光电工程》;20060531;第33卷(第5期);17-20 * |
开环光纤陀螺数字解调的实现与探讨;牟旭东;《压电与声光》;20001231;第22卷(第6期);360-362 * |
Also Published As
Publication number | Publication date |
---|---|
CN103616020A (zh) | 2014-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103604444B (zh) | 基于正弦波调制及二次谐波检测的光纤环本征频率测量装置及方法 | |
CN103616020B (zh) | 基于正弦波调制及一次谐波检测的光纤环本征频率测量装置及方法 | |
CN102840869B (zh) | 一种光纤陀螺本征频率的测量方法 | |
CN101339093B (zh) | 光纤陀螺用光纤环质量的测量方法及其装置 | |
CN102692314B (zh) | 基于光纤谐振腔测试激光器频率噪声功率谱密度的装置及方法 | |
CN102538822B (zh) | 光纤陀螺仪的动态特性快速测试标定方法 | |
CN107957276A (zh) | 基于频率漂移补偿的相位敏感光时域反射计及其测量方法 | |
CN102788595A (zh) | 基于Faraday效应的光纤陀螺频率特性评估方法与装置 | |
CN100458367C (zh) | 一种对光纤陀螺2π电压进行自动跟踪的四态调制及解调方法 | |
CN107356266B (zh) | 基于偶倍本征频率锯齿波调制的光纤陀螺本征频率测量方法 | |
CN100541127C (zh) | 采用不对称方波调制方法测量干涉式光纤陀螺本征频率的方法 | |
CN102506896B (zh) | 谐振腔技术测试谐振式光学陀螺中背散射噪声的方法 | |
CN101713701B (zh) | 基于外差技术测试相位调制器半波电压的装置及方法 | |
CN115077567B (zh) | 一种基于波导复位误差的标度因数补偿系统及方法 | |
CN104330104A (zh) | 一种干涉型传感器臂长差的测量装置 | |
CN100533064C (zh) | 光纤陀螺渡越时间测试方法 | |
CN102607547B (zh) | 光纤陀螺仪的数字域平衡探测方法及装置 | |
CN102183249B (zh) | 一种光纤陀螺仪的萨格奈克相移跟踪方法 | |
US5018859A (en) | Fiber optic gyroscope balanced plural serrodyne modulators phase difference control | |
Li et al. | Sagnac vibration sensing system with nested pulse method | |
CN113310483B (zh) | 一种数字闭环光纤陀螺本征频率实时跟踪装置及方法 | |
CN102901515A (zh) | 一种光纤陀螺渡越时间的在线快速测量方法 | |
CN1945209B (zh) | 一种抑制光纤陀螺交叉干扰的组合调制及解调方法 | |
CN102313558A (zh) | 基于Sagnac干涉的集成光学相位调制器直流漂移检测方法 | |
CN100386601C (zh) | Sagnac干涉型光纤传感器的一种开环信号处理方法和电路 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160106 Termination date: 20171126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |