CN106895958A - 一种用于高精度光纤测量的快速方法 - Google Patents

一种用于高精度光纤测量的快速方法 Download PDF

Info

Publication number
CN106895958A
CN106895958A CN201611270815.7A CN201611270815A CN106895958A CN 106895958 A CN106895958 A CN 106895958A CN 201611270815 A CN201611270815 A CN 201611270815A CN 106895958 A CN106895958 A CN 106895958A
Authority
CN
China
Prior art keywords
signal
dft
measurement
optical fiber
precision optical
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.)
Pending
Application number
CN201611270815.7A
Other languages
English (en)
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.)
Jiangsu Chun Long Photoelectric Polytron Technologies Inc
Original Assignee
Jiangsu Chun Long Photoelectric Polytron Technologies Inc
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 Jiangsu Chun Long Photoelectric Polytron Technologies Inc filed Critical Jiangsu Chun Long Photoelectric Polytron Technologies Inc
Priority to CN201611270815.7A priority Critical patent/CN106895958A/zh
Publication of CN106895958A publication Critical patent/CN106895958A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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/30—Testing of optical devices, constituted by fibre optics or optical waveguides
    • G01M11/31—Testing of optical devices, constituted by fibre optics or optical waveguides with a light emitter and a light receiver being disposed at the same side of a fibre or waveguide end-face, e.g. reflectometers
    • G01M11/3109—Reflectometers detecting the back-scattered light in the time-domain, e.g. OTDR
    • G01M11/3145—Details of the optoelectronics or data analysis

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种用于高精度光纤测量的快速方法，步骤包括：对待测光纤进行采样时增加采样点数；由辅助参考数据获得相位噪声信息；对测量数据去除相位噪声；对测量数据的频谱图进行分析，获得反射点信号对应的尖峰点。本发明在不影响精度的情况下，测量耗时减少了80％，将原来需要一分钟的测试时间缩短为十几秒。

Description

一种用于高精度光纤测量的快速方法
技术领域
本发明涉及光缆检测技术领域，特别是一种用于高精度光纤测量的快速方法。
背景技术
光频域反射技术的原始测量算法，为了能达到毫米级的测量精度，进行了多次测量的采样及运行算法，非常耗时，得出结果需要数分钟，限制了光频域反射仪的应用场景。在不影响测量精度的前提下，有效缩短测量时间，提高测试效率，让光频域反射仪能够满足不同的应用要求，本发明提出一种用于高精度光纤测量的快速算法，不影响精度的情况下单次测量即可获得原本对此测量的高精度测试结果。
发明内容
针对现有技术中存在的问题，本发明提供了一种在不影响精度的情况下，减少测量耗时的用于高精度光纤测量的快速方法。
本发明的目的通过以下技术方案实现。
一种用于高精度光纤测量的快速方法，步骤包括：
1)对待测光纤进行采样时增加采样点数；
2)由辅助参考数据获得相位噪声信息；
3)对测量数据去除相位噪声；
4)对测量数据的频谱图进行分析，获得反射点信号对应的尖峰点。
进一步的，步骤3)通过快速傅里叶变换：利用离散信号x(nT)来计算信号x(t)的频谱，有限长离散信号x(n)，n＝0，1，…，N-1的DFT定义为：
非周期性连续时间信号x(t)的傅里叶变换可以表示为：式中计算出来的是信号x(t)的连续频谱，将x(n)分解为偶数与奇数的两个序列之和，即x(n)＝x1(n)+x2(n)x1(n)和x2(n)的长度都是N/2，x1(n)是偶数序列，x2(n)是奇数序列，则
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT。由于X1(k)和X2(k)均以N/2为周期，且WN k+N/2＝-WN k，所以X(k)又可表示为：
依此类推，经过m-1次分解，最后将N点DFT分解为N/2个两点DFT。
相比于现有技术，本发明的优点在于：本发明在不影响精度的情况下，测量耗时减少了80％，将原来需要一分钟的测试时间缩短为十几秒。
经分析主要耗时的步骤部分为去除相位噪声，这部分改为使用快速傅里叶变换，时间大幅降低。将部分算法放进FPGA处理，也可加快数据处理速度。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合说明书附图和具体的实施例，对本发明作详细描述。
如图1所示，一种用于高精度光纤测量的快速方法，步骤包括：
1)对待测光纤进行采样时增加采样点数；
2)由辅助参考数据获得相位噪声信息；
3)对测量数据去除相位噪声；
4)对测量数据的频谱图进行分析，获得反射点信号对应的尖峰点。
步骤3)通过快速傅里叶变换：在实际的控制系统中能够得到的是连续信号x(t)的离散采样值x(nT)。因此需要利用离散信号x(nT)来计算信号x(t)的频谱。利用离散信号x(nT)来计算信号x(t)的频谱，有限长离散信号x(n)，n＝0，1，…，N-1的DFT定义为：
非周期性连续时间信号x(t)的傅里叶变换可以表示为：式中计算出来的是信号x(t)的连续频谱，可以看出，DFT需要计算大约N2次乘法和N2次加法。当N较大时，这个计算量是很大的。利用WN的对称性和周期性，将N点DFT分解为两个N/2点的DFT，这样两个N/2点DFT总的计算量只是原来的一半，即(N/2)2+(N/2)2＝N2/2，这样可以继续分解下去，将N/2再分解为N/4点DFT等。对于N＝2m点的DFT都可以分解为2点的DFT，这样其计算量可以减少为(N/2)log2N次乘法和Nlog2N次加法。
将x(n)分解为偶数与奇数的两个序列之和，即x(n)＝x1(n)+x2(n)x1(n)和x2(n)的长度都是N/2，x1(n)是偶数序列，x2(n)是奇数序列，则
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT。由于X1(k)和X2(k)均以N/2为周期，且WN k+N/2＝-WN k，所以X(k)又可表示为：
根据其形状称之为蝶形运算。依此类推，经过m-1次分解，最后将N点DFT分解为N/2个两点DFT。

Claims (2)

1.一种用于高精度光纤测量的快速方法，其特征在于步骤包括：
1)对待测光纤进行采样时增加采样点数；
2)由辅助参考数据获得相位噪声信息；
3)对测量数据去除相位噪声；
4)对测量数据的频谱图进行分析，获得反射点信号对应的尖峰点。
2.根据权利要求1所述的一种用于高精度光纤测量的快速方法，其特征在于步骤3)通过快速傅里叶变换：利用离散信号x(nT)来计算信号x(t)的频谱，有限长离散信号x(n)，n＝0，1，…，N-1的DFT定义为：非周期性连续时间信号x(t)的傅里叶变换可以表示为：式中计算出来的是信号x(t)的连续频谱，将x(n)分解为偶数与奇数的两个序列之和，即x(n)＝x1(n)+x2(n)x1(n)和x2(n)的长度都是N/2，x1(n)是偶数序列，x2(n)是奇数序列，则
由于则
其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT。由于X1(k)和X2(k)均以N/2为周期，且WN k+N/2＝-WNk，所以X(k)又可表示为：
依此类推，经过m-1次分解，最后将N点DFT分解为N/2个两点DFT。
CN201611270815.7A 2016-12-30 2016-12-30 一种用于高精度光纤测量的快速方法 Pending CN106895958A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611270815.7A CN106895958A (zh) 2016-12-30 2016-12-30 一种用于高精度光纤测量的快速方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611270815.7A CN106895958A (zh) 2016-12-30 2016-12-30 一种用于高精度光纤测量的快速方法

Publications (1)

Publication Number Publication Date
CN106895958A true CN106895958A (zh) 2017-06-27

Family

ID=59197791

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611270815.7A Pending CN106895958A (zh) 2016-12-30 2016-12-30 一种用于高精度光纤测量的快速方法

Country Status (1)

Country Link
CN (1) CN106895958A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101937424A (zh) * 2010-09-10 2011-01-05 南京航空航天大学 基于fpga实现高速fft处理的方法
CN102830282A (zh) * 2012-09-04 2012-12-19 国电南京自动化股份有限公司 一种2560点分组基2快速傅里叶变换方法
CN103856431A (zh) * 2014-03-04 2014-06-11 中国人民解放军理工大学 基于快速傅里叶变换的mpsk调制下开环载波同步算法
CN103852189A (zh) * 2012-11-29 2014-06-11 西安上尚机电有限公司 一种利用光纤测量配电柜温度的系统
CN105806594A (zh) * 2016-03-25 2016-07-27 江苏骏龙电力科技股份有限公司 一种毫米级分辨率的光纤检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101937424A (zh) * 2010-09-10 2011-01-05 南京航空航天大学 基于fpga实现高速fft处理的方法
CN102830282A (zh) * 2012-09-04 2012-12-19 国电南京自动化股份有限公司 一种2560点分组基2快速傅里叶变换方法
CN103852189A (zh) * 2012-11-29 2014-06-11 西安上尚机电有限公司 一种利用光纤测量配电柜温度的系统
CN103856431A (zh) * 2014-03-04 2014-06-11 中国人民解放军理工大学 基于快速傅里叶变换的mpsk调制下开环载波同步算法
CN105806594A (zh) * 2016-03-25 2016-07-27 江苏骏龙电力科技股份有限公司 一种毫米级分辨率的光纤检测方法

Similar Documents

Publication Publication Date Title
CN101701984B (zh) 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法
RU2668336C1 (ru) Способ определения места короткого замыкания на линиях электропередач
CN103956756A (zh) 一种电力系统低频振荡模态辨识方法
Lahiri et al. Estimating the parameters of multiple chirp signals
Yao et al. Improved discrete Fourier transform algorithm for harmonic analysis of rotor system
CN103543331B (zh) 一种计算电信号谐波和间谐波的方法
CN102269803A (zh) 基于时间延迟的离散频谱低频成分的校正方法
Li et al. Frequency estimation based on symmetric discrete Fourier transform
Ponomarev et al. Method for effective measurement of a sliding parametric Fourier spectrum
CN105372492B (zh) 基于三条dft复数谱线的信号频率测量方法
CN106895958A (zh) 一种用于高精度光纤测量的快速方法
Xu et al. xASTNN: Improved Code Representations for Industrial Practice
CN104849551B (zh) 一种谐相角分析方法
CN113688351A (zh) 用于检测物品重量的方法、装置、电子设备和可读介质
RU2551400C1 (ru) Способ гармонического анализа периодического многочастотного сигнала на фоне шума
Wang et al. Parameters estimation algorithm for the exponential signal by the interpolated all-phase DFT approach
CN113592010B (zh) 基于局部omp的快速卷积稀疏字典学习方法
Wu et al. A faster method for accurate spectral testing without requiring coherent sampling
Nur et al. High Performance Matched Filtering Algorithm Based on CUDA
CN113884761B (zh) 一种提高测频精度的方法、装置、设备及介质
CN104462017B (zh) I/o个数不相等的快速非均匀离散傅里叶变换方法及系统
US9471543B1 (en) Method and associated algorithm and firmware for extracting magnetometer signals from noise utilizing digital signal processing (DSP)
CN116165432A (zh) 一种数字测频方法、系统、存储介质和电子设备
Gillich et al. ABOUT THE ACCURACY OF ESTIMATED FREQUENCIES WITH THE PyFEST SOFTWARE
Li et al. Stress wave propagation analysis on vortex-induced vibration of marine risers

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170627