CN109141827A - A kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device - Google Patents
A kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device Download PDFInfo
- Publication number
- CN109141827A CN109141827A CN201710503922.8A CN201710503922A CN109141827A CN 109141827 A CN109141827 A CN 109141827A CN 201710503922 A CN201710503922 A CN 201710503922A CN 109141827 A CN109141827 A CN 109141827A
- Authority
- CN
- China
- Prior art keywords
- maximum value
- peak
- sequence
- value position
- optic filter
- 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.)
- Granted
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
Landscapes
- Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The invention discloses a kind of spectrum Peak Search Methods applied to tunable optic filter optical property detection device.It realizes that steps are as follows: data is grouped；Maximum value position is searched in every part of data；Maximum value position inceptive filtering；It constructs position difference sequence and carries out linear fit；Calculate filtration parameter array；Maximum value position filtering；Maximum value position stable convergence determines；Obtain peak position.The invention has the benefit that solving the problems, such as that existing fixed wave crest threshold value, fixed corrugation pitch Peak Search Method can not accurately be not fixed peak, peak period is not fixed, peak repressed spectrum in low light region carries out peak-seeking；Calculating process is simplified, is compared using simple numerical value, peak-seeking process can be completed in linear fit, has the advantages that peak-seeking is accurate, noise resisting ability is strong, speed is fast, adaptation is strong, convenient for the realization of circuit.
Description
Technical field
The present invention relates to fiber optic sensor technology fields, more particularly to a kind of tunable optic filter optical property that is applied to detect
The spectrum Peak Search Method of device.
Background technique
Using the optical property detection device of tunable optic filter in the process for being heated or being cooled down to tunable optic filter
In its can move through peak position, t is recorded using PD0To tnM- optical power is bent when the optical power value at moment can be obtained by
Line.Wavelength-optical power curve in order to obtain, m- wavelength curve, i.e., calibrate wavelength when needing to obtain.
In the wavelength scaling device using SLED light source+WDM+ etalon, SLED light source is penetrated through tunable optic filter
Peak acts on WDM and etalon, and when tunable optic filter occurs mobile through peak position, spectrum is modulated into tool by etalon
There is the curve of wave crest form, during tunable optic filter moves on to stopband from the passband of WDM through peak position, optical power is bent
There is low light region in line.Pass through the low light region position of appearance and crest location and WDM pass-band performance and etalon cyclophysis
Wavelength can be calibrated, m- wavelength curve when being obtained by fitting or interpolation, and then obtain wavelength-optical power curve.
It can be seen that accurate peak-seeking is the key that guarantee high-accuracy wavelength calibration.
In practical applications, since the Insertion Loss of tunable optic filter can change with temperature, in the process for heating or cooling down
In it is difficult to ensure that tunable optic filter Insertion Loss is constant or at the uniform velocity change, simultaneously because it is difficult to ensure that tunable filtering during temperature control
Device temperature at the uniform velocity changes, and leads to that optical power curve peak is not fixed, peak period is not fixed, this brings to accurate peak-seeking
It is difficult.Simultaneously because the presence of WDM, during tunable optic filter penetrates peak position from WDM passband to stopband, wave crest peak
Value can be suppressed, this increases difficulty to peak-seeking.
Summary of the invention
The purpose of the present invention is to provide a kind of spectrum peak-seeking sides applied to tunable optic filter optical property detection device
Method can be not fixed from spectrum peak, peak period is not fixed, is carried out in the repressed curve of low light region peak
The method of accurate peak-seeking, to obtain wavelength-optical power curve.
In order to achieve the above object, the technical scheme is that a kind of be applied to tunable optic filter optical property detection dress
The spectrum Peak Search Method set, it is characterised in that steps are as follows for the realization of this method:
Step 1: data are grouped, according to the cyclophysis of standard and tunable optic filter through peak movement speed and
The PD sampling period determines etc. that point interval I0, I0 are less than the half of minimum corrugation pitch, according to this equal part interval I0 by data point
At N parts, if data length cannot become portion by I0 equal part, remainder；
Step 2: its maximum value position L1, L2 ..., LN are searched in the every part of data generated in step 1；
Step 3: the interval I0 determined with step 1 the L1 to step 2, L2 ..., LN carry out inceptive filtering, filtering rule are as follows: each
Maximum value position L ' should be maximum value within the scope of the I0 of its left and right, if not then rejecting, obtain maximum value position by step 3
L1, L2 ..., LM；
Step 4: obtain D1 with the L2-L1 that step 3 obtains maximum value position, L3-L2 obtains D2, and so on obtain adjacent most
Big value alternate position spike sequence D 1, D2 ..., DM-1；With L1, L2 ..., LM-1 and D1, D2 ..., DM-1 carry out linear fit and obtain letter
Number D=f (L);
Step 5: by the L1 obtained function of substitution step 4 and divided by a coefficient n, obtaining f (L1)/n, rounding is denoted as I1, coefficient n
Determine the filter area of maximum value position, similarly, to maximum value position L2 ..., LM is handled obtain parameter in the same way
I2 ..., IM；
Step 6: the parameter I1 obtained according to step 5 is filtered L1, filtering rule are as follows: L1 is answered within the scope of the I1 of its left and right
For maximum value, if not then rejecting；Similarly, maximum value position L2 ..., LM are filtered according to parameter I2 ..., IM；It is logical
It crosses step 6 and obtains new maximum value sequence L1, L2 ..., LP；
Step 7: judge maximum value position sequence whether stable convergence, decision rule be by maximum value position sequence L1, L2 ...,
LP ' is denoted as L1, L2 ..., LM, substitutes into step 4, and 5,6 obtain new maximum value sequence, if new maximum value sequence and substitution step
4,5,6 sequence unanimously then thinks maximum value sequence stable convergence, enters step 8, if inconsistent, then it is assumed that maximum value sequence is not
Obtained new maximum value sequence is substituted into step 4 by stable convergence again, and 5,6 are calculated, until maximum value position sequence is steady
Fixed convergence, enters step 8；
Step 8: the stable convergence maximum value sequence obtained according to step 7 is spectral peak position, is denoted as P1, P2 ..., PQ.
Further, the interval I0 in the step 1 has just been fixed when system designs.
Further, the coefficient n in the step 5 has just been fixed when system designs.
The invention has the benefit that
1) present invention solves existing fixed wave crest threshold value, fixation corrugation pitch Peak Search Method and can not consolidate to peak accurately
The problem of fixed, peak period is not fixed, peak repressed spectrum in low light region carries out peak-seeking；
2) Peak Search Method of the invention simplifies calculating process, is compared using simple numerical value, peak-seeking can be completed in linear fit
Process has the advantages that peak-seeking is accurate, noise resisting ability is strong, speed is fast, adaptation is strong, convenient for the realization of circuit.
Detailed description of the invention
Fig. 1 is algorithm flow chart of the invention；
Fig. 2 is the spectrogram of present example；
Fig. 3 is maximum value position schematic diagram after present example spectrum array equal part；
Fig. 4 is the schematic diagram after the inceptive filtering of present example spectral maximum position；
Fig. 5 is present example spectral peak position schematic diagram.
Specific embodiment
With reference to the accompanying drawings and detailed description, the present invention will be further described.
Flow chart of the invention as shown in Figure 1,101 be data grouping module, and 102 be data maximums searching module, 103
It is maximum value position inceptive filtering module, 104 be that building position difference sequence carries out linear fit module, and 105 be filtration parameter number
Group computing module, 106 be maximum value position filtering module, and 107 be maximum value position stable convergence judgment module, and 108 be peak value
Assignment module.
As shown in Figure 2 with the spectrogram of 7500 point datas, it can be seen from the figure that spectrum has peak not solid
Fixed, peak period is not fixed, the repressed characteristic of low light region peak, traditional to fix wave crest threshold value, fixed wave crest
The peak-seeking mode of spacing is not suitable for the present embodiment.The present invention is as follows to the processing step of embodiment.
Step 1: data being grouped, peak movement speed is penetrated according to the cyclophysis of standard and tunable optic filter
And PD sampling period determination etc. point interval I0=25,7500 spectroscopic datas are divided into N=300 part；
Step 2: its maximum value position L1=11, L2=50 ..., L300=7477 are searched in the every part of data generated in step 1；
As shown in figure 3, with " * " mark at the maximum value position found；
Step 3: with the L1 for the I0=25 pair step 2 that step 1 determines, L2 ..., L300 carry out inceptive filtering, filtering rule are as follows: every
A maximum value position L ' should be maximum value in the range of I0=25 of its left and right, if not then rejecting.Step 2 is obtained by step 3
To 300 maximum value positions be filtered into 37, be denoted as L1=55, L2=127 ..., L37=7181, with " * " mark in Fig. 4；
Step 4: obtain D1=72 with the L2-L1 for obtaining maximum value position, L3-L2 obtains D2=89, and so on obtain adjacent
Maximum value position difference sequence D1=72, D2=89 ..., D36=373；With L1, L2 ..., L36 and D1, D2 ..., D36 carry out Linear Quasi
Conjunction obtains function D=0.03699 × L+95.6852;
Step 5: by the L1 obtained function of substitution step 4 and divided by a coefficient n, here according to light path system feature take n=
2.5, it obtains f (L1)/2.5, rounding is denoted as I1=39；Similarly, maximum value position L2 ..., L37 are handled obtain in the same way
Parameter I2=40 ..., I37=145；
Step 6: the parameter I1 obtained according to step 5 is filtered L1, filtering rule are as follows: L1 is answered within the scope of the I1 of its left and right
For maximum value, if not then rejecting；Similarly, maximum value position L2 ..., L37 are filtered according to parameter I2 ..., I37；
27 are filtered by 37 maximum value positions that step 6 obtains step 5, is denoted as L1=55, L2=127 ..., L27=7181,
With " * " mark in Fig. 5；
Step 7: judge maximum value position sequence whether stable convergence: by sequence L1=55, L2=127 ..., L27=7181 substitutes into step
Rapid 4,5,6 obtain new maximum value sequence L1 '=55, L2 '=127 ..., and L27 '=7181 determine new maximum value sequence L1 ', L2
' ..., L27 ' and L1, L2 ..., L27 is consistent, and maximum value sequence stable convergence enters step 8；
Step 8: the stable convergence maximum value sequence obtained according to step 7 is spectral peak position, is denoted as P1=55, and P2=
There are 27 wave crests in 127 ..., P27=7181 in spectrum.
Although specifically showing and describing the present invention in conjunction with preferred embodiment, those skilled in the art should be bright
It is white, it is not departing from the spirit and scope of the present invention defined by the appended claims, in the form and details to this hair
The bright various change made, is protection scope of the present invention.
Claims (3)
1. a kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device, which is characterized in that this method
It comprises the following steps that
Step 1: data are grouped, according to the cyclophysis of standard and tunable optic filter through peak movement speed and
The PD sampling period determines etc. that point interval I0, I0 are less than the half of minimum corrugation pitch, according to this equal part interval I0 by data point
At N parts, if data length cannot become portion by I0 equal part, remainder；
Step 2: its maximum value position L1, L2 ..., LN are searched in the every part of data generated in step 1；
Step 3: the interval I0 determined with step 1 the L1 to step 2, L2 ..., LN carry out inceptive filtering, filtering rule are as follows: each
Maximum value position L ' should be maximum value within the scope of the I0 of its left and right, if not then rejecting, obtain maximum value position by step 3
L1, L2 ..., LM；
Step 4: obtain D1 with the L2-L1 that step 3 obtains maximum value position, L3-L2 obtains D2, and so on obtain adjacent most
Big value alternate position spike sequence D 1, D2 ..., DM-1；With L1, L2 ..., LM-1 and D1, D2 ..., DM-1 carry out linear fit and obtain letter
Number D=f (L);
Step 5: by the L1 obtained function of substitution step 4 and divided by a coefficient n, obtaining f (L1)/n, rounding is denoted as I1, coefficient n
Determine the filter area of maximum value position, similarly, to maximum value position L2 ..., LM is handled obtain parameter in the same way
I2 ..., IM；
Step 6: the parameter I1 obtained according to step 5 is filtered L1, filtering rule are as follows: L1 is answered within the scope of the I1 of its left and right
For maximum value, if not then rejecting；Similarly, maximum value position L2 ..., LM are filtered according to parameter I2 ..., IM；It is logical
It crosses step 6 and obtains new maximum value sequence L1, L2 ..., LP；
Step 7: judge maximum value position sequence whether stable convergence, decision rule be by maximum value position sequence L1, L2 ...,
LP ' is denoted as L1, L2 ..., LM, substitutes into step 4, and 5,6 obtain new maximum value sequence, if new maximum value sequence and substitution step
4,5,6 sequence unanimously then thinks maximum value sequence stable convergence, enters step 8, if inconsistent, then it is assumed that maximum value sequence is not
Obtained new maximum value sequence is substituted into step 4 by stable convergence again, and 5,6 are calculated, until maximum value position sequence is steady
Fixed convergence, enters step 8；
Step 8: the stable convergence maximum value sequence obtained according to step 7 is spectral peak position, is denoted as P1, P2 ..., PQ.
2. a kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device as described in claim 1, special
Sign is that the interval I0 in the step 1 has just been fixed when system designs.
3. a kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device as described in claim 1, special
Sign is that the coefficient n in the step 5 has just been fixed when system designs.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710503922.8A CN109141827B (en) | 2017-06-28 | 2017-06-28 | Spectrum peak searching method applied to tunable filter optical performance detection device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710503922.8A CN109141827B (en) | 2017-06-28 | 2017-06-28 | Spectrum peak searching method applied to tunable filter optical performance detection device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109141827A true CN109141827A (en) | 2019-01-04 |
| CN109141827B CN109141827B (en) | 2020-07-07 |
Family
ID=64805431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710503922.8A Expired - Fee Related CN109141827B (en) | 2017-06-28 | 2017-06-28 | Spectrum peak searching method applied to tunable filter optical performance detection device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109141827B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112113673A (en) * | 2020-08-11 | 2020-12-22 | 福建华科光电有限公司 | Spectrum calibration method of optical tunable filter, readable storage medium and device |
| CN114047160A (en) * | 2021-10-28 | 2022-02-15 | 中南大学 | Second harmonic threading peak-searching method |
| CN114414048A (en) * | 2021-12-23 | 2022-04-29 | 电子科技大学 | Device and method for improving modulation spectrum measurement precision |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4375037A (en) * | 1980-01-07 | 1983-02-22 | Hitachi, Ltd. | Receiving circuit |
| EP0621686A3 (en) * | 1993-04-19 | 1995-10-18 | Nippon Telegraph & Telephone | Amplifier. |
| CN102435214A (en) * | 2011-09-20 | 2012-05-02 | 华北电力大学（保定） | High-precision and fast optical fiber grating spectrum peak-searching method |
| CN102706372A (en) * | 2012-03-28 | 2012-10-03 | 中国航空工业集团公司北京长城计量测试技术研究所 | Optical fiber wavelength demodulating light spectrum peak value positioning method |
| CN103578086A (en) * | 2013-11-14 | 2014-02-12 | 北京航空航天大学 | Method for restoring interference pattern data spectra based on wavelet analysis |
| CN103954368A (en) * | 2014-05-21 | 2014-07-30 | 北京遥测技术研究所 | Narrow-band light demodulating system and method based on photoelectric detection array |
| CN105823497A (en) * | 2016-05-24 | 2016-08-03 | 北京信息科技大学 | Fiber grating reflection spectrum demodulation algorithm based on signal autocorrelation matching |
-
2017
- 2017-06-28 CN CN201710503922.8A patent/CN109141827B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4375037A (en) * | 1980-01-07 | 1983-02-22 | Hitachi, Ltd. | Receiving circuit |
| EP0621686A3 (en) * | 1993-04-19 | 1995-10-18 | Nippon Telegraph & Telephone | Amplifier. |
| CN102435214A (en) * | 2011-09-20 | 2012-05-02 | 华北电力大学（保定） | High-precision and fast optical fiber grating spectrum peak-searching method |
| CN102706372A (en) * | 2012-03-28 | 2012-10-03 | 中国航空工业集团公司北京长城计量测试技术研究所 | Optical fiber wavelength demodulating light spectrum peak value positioning method |
| CN103578086A (en) * | 2013-11-14 | 2014-02-12 | 北京航空航天大学 | Method for restoring interference pattern data spectra based on wavelet analysis |
| CN103954368A (en) * | 2014-05-21 | 2014-07-30 | 北京遥测技术研究所 | Narrow-band light demodulating system and method based on photoelectric detection array |
| CN105823497A (en) * | 2016-05-24 | 2016-08-03 | 北京信息科技大学 | Fiber grating reflection spectrum demodulation algorithm based on signal autocorrelation matching |
Non-Patent Citations (1)
| Title |
|---|
| 李乔艺: ""可调谐DFB激光器动态扫描的FBG解调系统研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112113673A (en) * | 2020-08-11 | 2020-12-22 | 福建华科光电有限公司 | Spectrum calibration method of optical tunable filter, readable storage medium and device |
| CN114047160A (en) * | 2021-10-28 | 2022-02-15 | 中南大学 | Second harmonic threading peak-searching method |
| CN114414048A (en) * | 2021-12-23 | 2022-04-29 | 电子科技大学 | Device and method for improving modulation spectrum measurement precision |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109141827B (en) | 2020-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109141827A (en) | A kind of spectrum Peak Search Method applied to tunable optic filter optical property detection device | |
| CN102706372B (en) | Optical fiber wavelength demodulating light spectrum peak value positioning method | |
| CN103090961B (en) | Disturbance source positioning method of distributed type optical fiber sensing system | |
| CN103488751B (en) | Spectral matching method based on curve of spectrum wave-form similarity | |
| CN104573248A (en) | EMD based fiber-optic gyroscope temperature drift multi-scale extreme learning machine training method | |
| CN105823494A (en) | System and method of intelligently recognizing optical cable | |
| CN106027157A (en) | Current compensation method for optical module and optical module | |
| CN104111082A (en) | High-precision FBG sensing signal peak searching method | |
| CN102928003A (en) | Fiber grating demodulating system with real-time reference | |
| CN110160980A (en) | Analysis method, analytical equipment and the Systems for optical inspection of sample absorbance change rate | |
| CN105826811A (en) | Characterization method and device of tunable laser | |
| CN104579458A (en) | Light power compensation method and circuit | |
| CN106253973A (en) | A kind of distance less fundamental mode optical fibre characteristic measurement method and device | |
| CN106225816B (en) | A kind of grating sensing apparatus and method based on Brillouin's wave filter | |
| CN110987038B (en) | Distributed optical fiber disturbance detection device | |
| CN1256570C (en) | Method for realizing multi-channel optical fibre raster sensor high-sensitivity measurement | |
| CN108051018A (en) | Distributed fiber grating transducing signal peak-seeking control system based on FPGA | |
| CN103823175B (en) | A kind of photodetection circuit frequency response function test method based on OTDR | |
| CN103033205A (en) | Fiber bragg grating interrogator based on digitization tunable light sources and method thereof | |
| CN107543537A (en) | A kind of method for improving optic fiber gyroscope graduation factor stability | |
| CN108254000A (en) | A kind of FBG Sensor Network spectrum demodulation methods based on simulated annealing improved | |
| CN104462735A (en) | Method for simulating coaxially-packaged TOSA temperature distribution | |
| CN110048768A (en) | Identification system and identification method for same-wavelength optical fiber coding group | |
| CN107037583A (en) | Demodulate the methods, devices and systems of phase-shifted grating centre wavelength and phase-shift phase | |
| CN113934256A (en) | Variable-coefficient differential equation photon calculation solving system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200707 Termination date: 20210628 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |