CN107290057A - A kind of calibrating method for playing light modulation Fourier trasform spectroscopy - Google Patents
A kind of calibrating method for playing light modulation Fourier trasform spectroscopy Download PDFInfo
- Publication number
- CN107290057A CN107290057A CN201710572176.8A CN201710572176A CN107290057A CN 107290057 A CN107290057 A CN 107290057A CN 201710572176 A CN201710572176 A CN 201710572176A CN 107290057 A CN107290057 A CN 107290057A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- light modulation
- refractive index
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004611 spectroscopical analysis Methods 0.000 title claims description 12
- 239000013078 crystal Substances 0.000 claims abstract description 55
- 230000003287 optical effect Effects 0.000 claims abstract description 53
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 238000012937 correction Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 6
- 238000009966 trimming Methods 0.000 claims abstract description 5
- SBIBMFFZSBJNJF-UHFFFAOYSA-N selenium;zinc Chemical compound [Se]=[Zn] SBIBMFFZSBJNJF-UHFFFAOYSA-N 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 9
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- PFNQVRZLDWYSCW-UHFFFAOYSA-N (fluoren-9-ylideneamino) n-naphthalen-1-ylcarbamate Chemical compound C12=CC=CC=C2C2=CC=CC=C2C1=NOC(=O)NC1=CC=CC2=CC=CC=C12 PFNQVRZLDWYSCW-UHFFFAOYSA-N 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 239000002178 crystalline material Substances 0.000 description 2
- 239000010453 quartz Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 230000010748 Photoabsorption Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009615 fourier-transform spectroscopy Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 208000014733 refractive error Diseases 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- GTLQJUQHDTWYJC-UHFFFAOYSA-N zinc;selenium(2-) Chemical compound [Zn+2].[Se-2] GTLQJUQHDTWYJC-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/45—Interferometric spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
The present invention relates to playing in light modulation Fourier transform spectrometer, light modulation interference signal is played after discrete Fourier transform, the calibrating method of frequency spectrum, including：In light modulation Fourier transform spectrometer, is played, using one-wavelength laser as reference, the maximum optical path difference of laser interference figure is measured, using this maximum optical path difference as parameter, the frequency spectrum calibration for rebuilding spectrum is realized；But surveyed maximum optical path difference is related to the refractive index of the wavelength of incident light, crystal.Therefore, it is necessary to be corrected during the maximum optical path difference realization reconstruction spectral calibration measured with reference laser；Timing is usedFormula, the refractive index to various discrete point is corrected.Due to the n in formula0(N)、fc(N) it is unknown quantity, in trimming process, f is replaced using f (N)c(N) f, is calculatedc(N) refractive index corresponding to, then realizes the correction of frequency coordinate.Frequency after correction is again by testing the secondary correction for carrying out introducing error.The discrete frequency values of rebuild spectral intensity values and correction are finally based on, the calibration for rebuilding spectrum is realized.
Description
Technical field
The present invention relates to light modulation Fourier Transform Spectroscopy field is played, more specifically, it is related to a kind of bullet light modulation
The calibrating method of Fourier trasform spectroscopy, the calibrating method of the spectrum is related to the maximum optical path difference and folding that play light modulation interference pattern
Penetrate rate, the relation of wavelength；Various discrete frequency values after Fourier transformation are carried out with maximum optical path difference parameter acquiring with reference light
The problem of calibration needs to carry out compensation of refractive index；And the frequency compensation method based on refractive index.
Background technology
Have that measuring speed is fast, spectrum model to play the bullet light modulation Fourier transform spectrometer, of optical modulator as interferometer
The advantages of wide, anti-seismic performance is good is enclosed, there is potential using value in transient state spectrographic detection field.Optical modulator is played to adjust as bullet light
The critical piece of interferometer processed, is made up of bullet luminescent crystal and piezoelectric quartz crystal.Wherein, luminescent crystal is played as thang-kng crystal,
It realizes the regulation of light beam based on elasto-optical effect, produces interference pattern.
In order to based on light modulation interference pattern is played, rebuild the spectral information of incoming signal, it is necessary to detect the maximum light of interference pattern
Path difference parameter.In conventional time-modulation type Fourier transform spectrometer, generally using reference laser measurement interference pattern most
Big optical path difference.In Fourier transformation, the frequency resolution of Fourier transformation is described with frequency sampling interval delta f, represents spectrum point
The minimum interval of discernmible two spectrum components in analysis, and this interval and continuous signal time domain duration T into anti-
Than.Continuous signal is reflected in interference pattern in the duration of time domain, then is the maximum optical path difference of interference pattern.Play optical modulator
The refractive index for playing luminescent crystal is not steady state value, and it is the variable relevant with the wavelength of incident light, temperature.This is allowed in identical
Under strained condition, for different incident optical signals, it is different to play the phase-delay quantity of optical modulator generation.
Typically with 808nm laser as reference light source, measurement plays the maximum optical path difference of light modulation interference pattern, and is based on
This parameter, the discrete interval Frequency point of setting frequency domain will bring frequency spectrum calibrated error.In order to accurately rebuild bullet light modulation interference
The spectrogram of figure, realizes the accurate calibration of frequency coordinate, it is necessary to which the calibrating method for playing photoabsorption modulation light spectrogram is improved.
The content of the invention
In order in the presence of overcoming the shortcomings of prior art, the present invention, which is provided, a kind of plays light modulation Fourier trasform spectroscopy
Calibrating method, this method can realize that rebuilding spectrum and frequency spectrum based on the Fourier transformation interference pattern for playing optical modulator are accurately fixed
Mark.
In order to solve the above-mentioned technical problem, the technical solution adopted in the present invention is：
A kind of calibrating method for playing light modulation Fourier trasform spectroscopy, is followed the steps below：
Step one：One-wavelength laser is used for light source, bullet optical modulator is thang-kng crystal, and bullet light is calculated using following formula
The maximum optical path difference of interferometric modulator figure,
In above formula, d is plays the thickness of optical modulator optical direction, and λ is the wavelength of incident light, nr0(λ) is crystal in the ripple
Refractive index under long, δ0For maximum modulation stress, π12,π11For the stress elasto-optical coefficient tensor of different optical directions；
Step 2：Using maximum optical path difference obtained above as parameter, the mark of discrete Fourier transform spectrum abscissa is obtained
Determine parameter；
Step 3：Each Frequency point of discrete Fourier transform spectrum abscissa, is the static state folding with the discrete point in frequency
Penetrate rate and acquisition is corrected to the calibrating parameters of discrete Fourier transform spectrum abscissa.
Step 4：Frequency point coordinates based on correction, realizes the frequency spectrum calibration for rebuilding spectrum.
In the step 2, in discrete Fourier transform, each frequency values calculated are：
In above formula, N is discrete Frequency point, and K is the points of discrete Fourier transform, and f (N) is discrete frequency values.
In the step 3, the error based on refractive index, the correction of various discrete frequency values is public after discrete Fourier transform
Formula is：
In above formula, nr0=2.5242, the value is refractive index of the ZnSe crystal in 808nm wavelength in the environment of 23 DEG C,
fc(N) it is the discrete frequency values after correction；N (N) is the refractive index of n-th discrete point in frequency after correction.
In formula, n0(N)、fc(N) it is unknown quantity, in trimming process, f is replaced using f (N)c(N), by following
Formula calculates fc(N) refractive index corresponding to：
After power correction, then using the monochromatic light of different wave length as signal source, two are carried out to the error introduced in algorithm
Secondary correction, modified index of refraction parameter.
The bullet optical modulator is made up of piezo-electric crystal and ZnSe crystal, and ZnSe crystal, as thang-kng element, is horizontal bullet
Luminous effect.
The one-wavelength laser uses the narrow band laser of 808nm wavelength.
Compared with prior art, the advantageous effect of present invention is that：
Static refractive index of the present invention based on bullet luminescent crystal is function of wavelength, and the maximum light path measured with reference laser
Difference, realizes the abscissa calibration of the Frequency point of discrete Fourier transform.Using this method, it can be achieved to play light modulation interference pattern non-
After uniform Fast Fourier Transform (FFT), frequency scaling is carried out.
Brief description of the drawings
Fig. 1 is the compound bullet light modulation interference pattern of multiple laser of narrowband of the present invention；
The index of refraction diagram of Fig. 2 Znse crystal at different wavelengths；
Fig. 3 is the reconstruction spectrogram demarcated with reference laser of the invention；
Fig. 4 is the calibration spectrogram after the power correction of the present invention；
Fig. 5 is the frequency values and the comparison figure of true value after the correction of the present invention；
Fig. 6 is the calibration spectrogram after the error compensation of the present invention.
Embodiment：
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
The present invention bullet optical modulator be made up of piezo-electric crystal and ZnSe crystal, ZnSe crystal as thang-kng element,
It is horizontal elasto-optical effect.Under stress drive, play luminescent crystal and produce birefringence effect, its refringence is and answering for being subject to
Power, static refractive index n0, stress elasto-optical coefficient tensor it is relevant.And after crystalline material determination, stress elasto-optical coefficient tensor is solid
Determine parameter.
The maximum optical path difference of the bullet light modulation interference pattern of the present invention is measured by 808nm laser interference figure, and is applied
The maximum optical path difference of measurement as discrete Fourier transform spectrum abscissa calibrating parameters.
Each Frequency point of the abscissa of discrete Fourier transform of the present invention, should be using the static refringence of the frequency as ginseng
Number, measures maximum optical path difference, rather than with the static refractive index of 808nm reference wavelength, calculates maximum optical path difference.Therefore, need
The static refractive index for influenceing maximum optical path difference parameter is corrected.
The present invention each discrete point in frequency updating formula be：
In formula, nr=2.5242, it is refractive index of the ZnSe crystal in 808nm wavelength under 23 DEG C of environment；n0(N) it is each
The static refractive index of discrete point in frequency.
In above formula, n0(N)、fc(N) it is unknown quantity, in trimming process, f is replaced using f (N)c(N) f, is calculatedc
(N) refractive index corresponding to.The relation of refractive index and frequency is：
After correction, then by testing progress error correction.
To make the purpose of the present invention, technical scheme, feature more protrude, the specific implementation below in conjunction with accompanying drawing to invention
It is illustrated.
In light modulation Fourier transform spectrometer, is played, optical modulator is played as its static interferometer, when radiating light source warp
The polarizer, produces two beam o light and e light, and this two-beam produces phase delay, realizes two-beam interference, produce when by modulator
It is raw to play light modulation interference pattern.
It is the Primary Component for producing interference pattern to play optical modulator, and the analysis to its characteristic helps to understand that playing light modulation interferes
Figure.
It is the resonating device being made up of bullet luminescent crystal ZnSe and quartz crystal to play optical modulator.Its operation principle is：Piezoelectricity
Crystal produces strain in the presence of extraneous electric signal based on inverse piezoelectric effect, and the strain causes the inverse dielectric for playing luminescent crystal
Measure ηijMicrovariations.
ZnSe crystal belongs to cubic systemCrystal class, for designed bullet optical modulator, its machinery being subject to
Stress is single direction, and physical limit, therefore the bullet light is not present in crystal at free vibrational flexing pattern, crystal edge
The stress elasto-optical coefficient tensor of crystal only has π11π12π13π44, wherein,
The stress elasto-optical coefficient matrix of zinc selenide crystal can be written as：
It is assumed that to zinc selenide crystal along x1Direction of principal axis produces strain σ1, the relation of dielectric isolation rate tensor and stress tensor is：
Δβm=βm-β0 m=πmnσn, (m, n=1,2 ..., 6) (2)
By formula (1) and (2), it can obtain：
Because zinc selenide crystal belongs toCrystal class, as cubic crystal, thereforeBy formula (3) band
Entering index ellipsoid equation has：
The formula shows, crystal refractive index ellipsoid is changed into three-axial ellipsoid from original spheroid, that is, crystal from optics it is each to
Same sex body becomes biaxal crystal, but three major axes orientations do not change.
Can solve three principal refractive indexs using formula (4) is respectively
If along x1Direction of principal axis thang-kng (longitudinal effect), then phase difference is
Wherein, d is the thickness of crystal optical direction, and λ is the wavelength of incident light；
If along x2And x3Direction of principal axis thang-kng (transversal effect), phase difference is respectively
There is π for homogeneous material zinc selenide12=π13, thus longitudinal elasto-optical effect is not present in it.
In this programme, from x3Direction of principal axis thang-kng, final refringence and position difference are：
For polygon-octagonal bullet ray crystal structure, its resonance mode is two-dimentional in plane vibration, and its strain is one-dimensional vibration
Twice, thus refringence and position difference are twice of one-dimensional case, i.e.,：
(10), (11) provide refringence and position difference and the relation of stress.
It was found from formula (10), the refringence for playing optical modulator is stress, the refractive index n with being subject to0, stress elasto-optical coefficient
Tensor is relevant.
Meanwhile, the refractive index of crystal is the function n of wavelength0(λ), for the incident light of different wave length, its refractive index is not
Identical；And after crystalline material is fixed, stress elasto-optical coefficient tensor is fixed parameter.
For one-dimensional bullet optical modulator, its stress distribution is：
δ (x, t)=δ0cos(πx/L)sin(ωt-ψ) (12)
Wherein, δ0It is maximum modulation stress, t is the time, and ψ is initial phase, and L is the length for playing luminescent crystal, space variable x's
Scope is from 0 to L.
The central point of luminescent crystal is being played, the strain minimum of generation, stress are maximum.If initial phase is zero, then have：
δ (t)=δ0sin(ωt)
Above formula is substituted into (10), (11), while considering that refractive index is the function of wavelength, then had：
As can be seen that playing the timing of thickness one of luminescent crystal from phase delay formula (11), the transient state of optical modulator is played
Birefringence difference Δ n is bigger, and phase-delay quantity is bigger.The maximum optical path difference for producing interference pattern is also bigger, so that rebuilding
The spectral resolution of spectrum is improved.
Monochromatic light is after playing the modulation of light modulation interferometer, and the interference pattern of generation can be expressed as with equation below：
When incident light is polychromatic light, through playing light modulation interferometer, produces interference pattern and be represented by：
From formula (16), it is assumed that refractive index is definite value, then interference pattern only has with incident optical signal, driving stress signal
Close.
Using 808nm (wave number 12376) laser as reference, by wave number be 2020,3000,5000,8000,15000 six
Plant narrow band light to be combined, interference pattern such as Fig. 1 of generation.
In order to rebuild the spectrogram of radiation source by interference pattern, it need to be realized using inverse discrete Fourier transform algorithm.
In Fourier transformation, the frequency resolution of Fourier transformation is described with frequency sampling interval delta f, represents analysis of spectrum
In discernmible two spectrum components minimum interval, and the duration T in time domain of this interval and continuous signal is into anti-
Than.Continuous signal is reflected in interference pattern in the duration of time domain, then is the maximum optical path difference of interference pattern.
In order to realize the frequency demarcating of abscissa, then need to measure the maximum optical path difference of interference pattern.
In spectrum calibration method of the present invention, it is 808nm narrow band laser as reference light source to use wavelength,
Measure the transient state maximum optical path difference of interference pattern.
Play light modulation interference pattern transient state optical path difference be：
D is the thickness of crystal optical direction in formula, and λ is the wavelength of incident light；n0(λ) is refraction of the crystal in the wavelength
Rate；δ0It is maximum modulation stress；π12,π11The stress elasto-optical coefficient tensor of different optical directions.
Play light modulation interference pattern maximum optical path difference be：
N in formular0(λ) is refractive index of the reference laser in the crystal.
In Fourier transformation frequency spectrum, its frequency resolution is：
In discrete Fourier transform, the index of each Frequency point is：
In formula, N is discrete Frequency point, and K is the points of discrete Fourier transform, and f (N) is discrete frequency values.
From formula (20), under driving voltage, crystalline size certain condition, the coordinate of each Frequency point is and different ripples
The cube of refractive index is related under long (wave number).And the refractive index of crystal at different wavelengths is a variable, such as Fig. 2
It is shown.
In light modulation Fourier transform spectrometer, is played, wavelength is used for 808nm laser measurement maximum optical path difference, and in fact
Existing frequency spectrum calibration.Fig. 3 is with the spectrogram of 808nm laser scaling.Because introducing refractive error in this calibration, so weight
The spectrogram wave number built is wrong.
In order to overcome the abscissa error that refractive index is introduced, it need to be corrected using the refractive index parameter of each Frequency point.
ZnSe crystal is as thang-kng crystal, and its spectral window is 0.54~18.2 μm (wave number is 549.45~18518), folding
The relation for penetrating rate and wavelength is：
In formula, α1=4.45813734, α2=0.467216334, α3=2.89566290, β1=0.200859853, β2=
0.391371166, β3=47.1362108
The relation for being converted to refractive index and wave number is：
In formula, k=10000.According to formula (22), the refractive index under different wave numbers can be calculated.
From formula (12), each frequency coordinate of discrete Fourier transform is inversely proportional with the cube of refractive index,
In order to eliminate the error introduced by reference laser refractive index, it can be corrected using formula (24).
In formula, nr0=2.5242, it is refractive index of the ZnSe crystal in 808nm wavelength under 23 DEG C of environment；fc(N) it is
Discrete frequency values after correction；N (N) is the refractive index of n-th discrete point in frequency after correction.
In formula (24), n (N), fc(N) it is unknown quantity, in view of refractive index is gradual process, in trimming process
In, using f (N), its corresponding refractive index is calculated using formula (23).Instead of the refractive index n (N) of the check point in formula.
Based on the bearing calibration of formula (24), the spectrogram after being corrected, as shown in Figure 4.Because there are two unknown numbers,
The error of some calibrations can be introduced, as shown in Figure 5.Error after near infrared band, correction is smaller, meets bullet light tune
The index request of Fourier transform spectrometer, processed；And it is big in visible light wave range application condition.Need to carry out secondary correction.
Secondary correction, using the method for experiment correction, refractive index is modified so that whole service band is met will
Ask, as shown in Figure 6.
The foregoing is only the present invention is preferable to carry out technology and data processing method, is not limited in the present invention.To this
For the technical staff in field, on the premise of the implementation principle of the present invention is not departed from, improvement, the polishing made fall within this
The protection domain of invention.
Claims (7)
1. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy, it is characterised in that：
Step one：One-wavelength laser is used for light source, bullet optical modulator is thang-kng crystal, and bullet light modulation is calculated using following formula
The maximum optical path difference of interference pattern,
<mrow>
<mi>L</mi>
<mo>=</mo>
<msubsup>
<mi>dn</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
<mn>3</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;</mi>
<mn>12</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&pi;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mn>0</mn>
</msub>
</mrow>
In above formula, d is plays the thickness of optical modulator optical direction, and λ is the wavelength of incident light, nr0(λ) be crystal at that wavelength
Refractive index, δ0For maximum modulation stress, π12,π11For the stress elasto-optical coefficient tensor of different optical directions；
Step 2：Using maximum optical path difference obtained above as parameter, the demarcation ginseng of discrete Fourier transform spectrum abscissa is obtained
Number；
Step 3：Each Frequency point of discrete Fourier transform spectrum abscissa, is the static refractive index with the discrete point in frequency
Acquisition is corrected to the calibrating parameters of discrete Fourier transform spectrum abscissa.
Step 4：Frequency point coordinates based on correction, realizes the frequency spectrum calibration for rebuilding spectrum.
2. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 1, it is characterised in that：It is described
In step 2, in discrete Fourier transform, each frequency values calculated are：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>N</mi>
<mrow>
<mn>2</mn>
<msubsup>
<mi>dn</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
<mn>3</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;</mi>
<mn>12</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&pi;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&delta;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mn>0</mn>
<mo>:</mo>
<mn>1</mn>
<mo>:</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
In above formula, N is discrete Frequency point, and K is the points of discrete Fourier transform, and f (N) is discrete frequency values.
3. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 1, it is characterised in that：It is described
In step 3, the error based on refractive index, the updating formula of various discrete frequency values is after discrete Fourier transform：
<mrow>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>n</mi>
<mn>0</mn>
</msub>
<mn>3</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mi>n</mi>
<mrow>
<mi>r</mi>
<mn>0</mn>
</mrow>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
</mrow>
In above formula, nr0=2.5242, the value is refractive index of the ZnSe crystal in 808nm wavelength, f in the environment of 23 DEG Cc
(N) it is the discrete frequency values after correction；N (N) is the refractive index of n-th discrete point in frequency after correction.
4. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 3, it is characterised in that：In public affairs
In formula, n0(N)、fc(N) it is unknown quantity, in trimming process, f is replaced using f (N)c(N) f, is calculated by below equationc
(N) refractive index corresponding to：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<mi>f</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<mi>f</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&beta;</mi>
<mn>3</mn>
</msub>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<mi>f</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msqrt>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mn>0</mn>
<mo>:</mo>
<mn>1</mn>
<mo>:</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
5. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 4, it is characterised in that：In folding
Penetrate after rate correction, then using the monochromatic light of different wave length as signal source, secondary correction, amendment folding are carried out to the error introduced in algorithm
Penetrate rate parameter.
6. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 1, it is characterised in that：It is described
Play optical modulator to be made up of piezo-electric crystal and ZnSe crystal, ZnSe crystal, as thang-kng element, is horizontal elasto-optical effect.
7. a kind of calibrating method for playing light modulation Fourier trasform spectroscopy according to claim 1, it is characterised in that：It is described
One-wavelength laser uses the narrow band laser of 808nm wavelength.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710572176.8A CN107290057B (en) | 2017-07-13 | 2017-07-13 | A kind of calibrating method playing light modulation Fourier trasform spectroscopy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710572176.8A CN107290057B (en) | 2017-07-13 | 2017-07-13 | A kind of calibrating method playing light modulation Fourier trasform spectroscopy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107290057A true CN107290057A (en) | 2017-10-24 |
| CN107290057B CN107290057B (en) | 2018-08-31 |
Family
ID=60101253
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710572176.8A Expired - Fee Related CN107290057B (en) | 2017-07-13 | 2017-07-13 | A kind of calibrating method playing light modulation Fourier trasform spectroscopy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107290057B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113544491A (en) * | 2019-03-15 | 2021-10-22 | 布鲁克光学有限公司 | Method for determining a correction quantity function and method for generating a frequency-corrected hyperspectral image |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007061436A1 (en) * | 2005-11-28 | 2007-05-31 | University Of South Carolina | Self calibration methods for optical analysis system |
| CN103162834A (en) * | 2013-03-05 | 2013-06-19 | 中北大学 | Multi-reflection single elastic light modulation Fourier transformation spectrometer light path structure |
| CN103234635A (en) * | 2013-03-22 | 2013-08-07 | 中北大学 | Photoelastic-modulation Fourier transform interference imaging spectrometer |
| CN203719771U (en) * | 2014-01-03 | 2014-07-16 | 南京邮电大学 | Spectral measurement apparatus based on elasto-optical effect |
-
2017
- 2017-07-13 CN CN201710572176.8A patent/CN107290057B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007061436A1 (en) * | 2005-11-28 | 2007-05-31 | University Of South Carolina | Self calibration methods for optical analysis system |
| CN103162834A (en) * | 2013-03-05 | 2013-06-19 | 中北大学 | Multi-reflection single elastic light modulation Fourier transformation spectrometer light path structure |
| CN103234635A (en) * | 2013-03-22 | 2013-08-07 | 中北大学 | Photoelastic-modulation Fourier transform interference imaging spectrometer |
| CN203719771U (en) * | 2014-01-03 | 2014-07-16 | 南京邮电大学 | Spectral measurement apparatus based on elasto-optical effect |
Non-Patent Citations (1)
| Title |
|---|
| 张敏娟: "《弹光调制傅里叶变换光谱复原高速数据处理技术研究》", 《中国博士学位论文全文数据库》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113544491A (en) * | 2019-03-15 | 2021-10-22 | 布鲁克光学有限公司 | Method for determining a correction quantity function and method for generating a frequency-corrected hyperspectral image |
| US11994430B2 (en) | 2019-03-15 | 2024-05-28 | Bruker Optics Gmbh & Co. Kg | Method for determining a correction value function and method for generating a frequency-corrected hyperspectral image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107290057B (en) | 2018-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110261317B (en) | Measurement system and method for Mueller matrix spectrum | |
| CN104535192B (en) | Stokes polarimeter error calibration based on rotating wave plate method and compensation method | |
| CN102426058B (en) | Static interference imaging polarizer and method for acquiring polarization information of target | |
| CN106706521A (en) | Optical film ultra wide band optical constant testing method | |
| CN104677834A (en) | Method for carrying out optical measurement by using full-Mueller matrix ellipsometer | |
| CN111443045B (en) | Spectral imaging analysis system and spectral imaging analysis method | |
| CN104677838A (en) | Self-calibration type full-muller-matrix ellipsometer measurement system | |
| CN104568391A (en) | Performance testing method and device for dual optical path switching mutual reference high precession AOTF | |
| KR20080091002A (en) | Apparatus for measuring phase difference using spectrometer | |
| CN114136894B (en) | Error calibration method and device of polarization detection system based on vortex wave plate | |
| CN104677833A (en) | Method for carrying out optical measurement by using full-Mueller matrix ellipsometer | |
| CN107290057B (en) | A kind of calibrating method playing light modulation Fourier trasform spectroscopy | |
| CN204855372U (en) | Heavy -calibre uniaxial crystal refracting index homogeneity measuring device | |
| CN204679246U (en) | Double light path switches mutually with reference to high-precision A OTF performance testing device | |
| CN104677837A (en) | Full-Mueller matrix ellipsometer calibration method | |
| CN103411890A (en) | Rotating compensator type ellipsometer system error assessment and elimination method thereof | |
| CN110567883B (en) | Method for measuring Mueller matrix spectrum in real time | |
| Zhai et al. | Study on high sensitivity measurement of seawater temperature based on bow tie fiber | |
| JP5041508B2 (en) | Optical characteristic measuring apparatus and method | |
| CN104677835A (en) | Full-Mueller matrix ellipsometer calibration method | |
| CN105158209A (en) | Large aperture uniaxial crystal refractive index uniformity measurement apparatus and method thereof | |
| Naciri et al. | Fixed polarizer, rotating-polarizer and fixed analyzer spectroscopic ellipsometer: accurate calibration method, effect of errors and testing | |
| CN108801604A (en) | It is a kind of play optical modulator phase prolong amplitude calibration with closed-loop control device and method | |
| CN107314839A (en) | Stress detection device and method based on Muller matrix | |
| CN111812032A (en) | Stress measuring instrument for wide-spectrum optical material and element |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Hao Jian Inventor after: Zhang Minjuan Inventor after: Li Jinhua Inventor after: Jing Ning Inventor after: Wang Zhibin Inventor after: Zhang Rui Inventor after: Li Kewu Inventor after: Su Limei Inventor before: Zhang Minjuan Inventor before: Hao Jian Inventor before: Li Jinhua Inventor before: Jing Ning Inventor before: Wang Zhibin Inventor before: Zhang Rui Inventor before: Li Kewu Inventor before: Su Limei |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180831 Termination date: 20190713 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |