CN105823559B - A kind of adaptive double light comb spectrally compensating method for extracting signal - Google Patents
A kind of adaptive double light comb spectrally compensating method for extracting signal Download PDFInfo
- Publication number
- CN105823559B CN105823559B CN201610308996.1A CN201610308996A CN105823559B CN 105823559 B CN105823559 B CN 105823559B CN 201610308996 A CN201610308996 A CN 201610308996A CN 105823559 B CN105823559 B CN 105823559B
- Authority
- CN
- China
- Prior art keywords
- frequency
- optical
- comb
- optical signals
- signals
- 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.)
- Active
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title abstract description 8
- 230000003595 spectral effect Effects 0.000 claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 20
- 238000000605 extraction Methods 0.000 claims abstract description 13
- 230000003287 optical effect Effects 0.000 claims description 113
- 239000000835 fiber Substances 0.000 claims description 41
- 238000001914 filtration Methods 0.000 claims description 22
- 239000013078 crystal Substances 0.000 claims description 21
- 239000013307 optical fiber Substances 0.000 claims description 8
- 230000035559 beat frequency Effects 0.000 abstract description 13
- 238000004611 spectroscopical analysis Methods 0.000 abstract description 11
- 238000005259 measurement Methods 0.000 abstract description 10
- 239000000284 extract Substances 0.000 abstract description 2
- 230000009977 dual effect Effects 0.000 description 12
- 238000001514 detection method Methods 0.000 description 11
- 210000001520 comb Anatomy 0.000 description 3
- 230000007774 longterm Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 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)
Abstract
本发明公开了一种自适应双光梳光谱补偿信号提取系统,其特点是迈克尔逊干涉仪和电学拍频获得自适应双光梳光谱补偿信号,采用时域延时和频域选择的方式,有效利用迈克尔逊干涉仪相干探测一定延时的前后脉冲,提取出表征单个脉冲激光器自身重复频率抖动和载波包络相位抖动的信号,最后通过电学倍频及混频过程,获得表征两个脉冲激光器相对重复频率抖动和相对载波包络相位抖动的补偿信号,可直接用于自适应双光梳光谱测量。本发明与现有技术相比完全避免了自适应双光梳光谱中需要引入两台窄线宽连续激光器提取补偿信号的方式,补偿信号受连续光波长飘移干扰小,信号稳定可靠,进一步提高了自适应双光梳光谱系统的稳定性和可靠性。
The invention discloses an adaptive dual-comb spectral compensation signal extraction system, which is characterized in that the adaptive dual-comb spectral compensation signal is obtained by Michelson interferometer and electrical beat frequency, and adopts the method of time domain delay and frequency domain selection, Effectively use the Michelson interferometer to coherently detect the front and rear pulses with a certain delay, and extract the signals representing the repetition frequency jitter and carrier envelope phase jitter of a single pulse laser, and finally obtain the two pulse lasers through the process of electrical frequency doubling and frequency mixing Compensation signals for relative repetition frequency jitter and relative carrier envelope phase jitter can be directly used for adaptive dual-comb spectroscopy measurement. Compared with the prior art, the present invention completely avoids the need to introduce two narrow-linewidth continuous lasers to extract the compensation signal in the adaptive dual-comb spectrum. The compensation signal is less disturbed by the continuous light wavelength drift, and the signal is stable and reliable, which further improves the Stability and reliability of an adaptive dual-comb spectroscopy system.
Description
技术领域technical field
本发明涉及光电探测技术领域,具体地说是一种基于基于迈克尔逊干涉仪的自适应双光梳光谱补偿信号提取方法。The invention relates to the technical field of photoelectric detection, in particular to an adaptive dual-comb spectral compensation signal extraction method based on Michelson interferometer.
背景技术Background technique
双光梳光谱技术作为近年来科研领域的前沿课题,能极大的提高光谱探测的精度。双光梳光谱探测技术相比传统技术有两个改进,一是采用光梳替代传统光源,作为标准的频率谱线位置更稳定,线宽更窄,精度更高;另一个是采用两个重复频率略有差异的光梳同时测量,一个作为参考光梳,一个作为激发光梳,类似于游标卡尺的主尺和副尺,通过频率错位和拍频探测,进一步提高频率测量的精度。近期发展起来的自适应双光梳光谱技术,采用两台窄线宽连续激光器和两台超短脉冲激光器拍频,再经过电路混频、倍频、放大、滤波处理,获得表征两台超短脉冲激光器相对重复频率抖动和相对载波包络相位抖动的补偿信号,表征相对重复频率抖动的信号用于光谱探测过程中的异步采样,表征相对载波包络相位抖动的信号用于光谱干涉信号混频,最终在信号处理阶段消除抖动的影响,实现高精度光谱探测。As a frontier topic in the field of scientific research in recent years, dual-comb spectroscopy technology can greatly improve the accuracy of spectral detection. Compared with the traditional technology, the dual optical comb spectral detection technology has two improvements. One is to use an optical comb to replace the traditional light source. As a standard, the position of the frequency spectral line is more stable, the line width is narrower, and the accuracy is higher; Optical combs with slightly different frequencies are measured at the same time. One is used as a reference comb and the other is used as an excitation comb, which is similar to the main scale and auxiliary scale of a vernier caliper. Through frequency misalignment and beat frequency detection, the accuracy of frequency measurement is further improved. The recently developed adaptive dual-comb spectroscopy technology uses two narrow-linewidth CW lasers and two ultrashort pulse lasers to beat frequency, and then undergoes circuit frequency mixing, frequency doubling, amplification, and filtering to obtain and characterize the two ultrashort pulse lasers. Compensation signals for relative repetition frequency jitter and relative carrier envelope phase jitter of pulsed lasers, signals representing relative repetition frequency jitter are used for asynchronous sampling in the process of spectral detection, signals representing relative carrier envelope phase jitter are used for spectral interference signal mixing , and finally eliminate the influence of jitter in the signal processing stage to achieve high-precision spectral detection.
目前的自适应双光梳光谱系统中,补偿信号的源头来自两台窄线宽连续激光器和两台超短脉冲激光器的相互光学拍频,由于环境温度、振动变化,每台窄线宽激光器或者超短脉冲激光器的光频会随时间随机飘移,因此,光学拍频信号的中心频率也会随时间随机飘移,并且飘移范围非常大,很容易超出电路滤波器的带宽,导致拍频信号丢失,电路处理模块失效,自适应补偿系统工作异常。为控制窄线宽连续激光器和超短脉冲激光器的光学拍频信号飘移范围,可以通过建立负反馈环路,通过调节连续激光器的温度、泵浦功率或者脉冲激光器的重复频率等参数,将两者的拍频飘移控制在后续电路滤波器的带宽之内,但是这种负反馈的补偿方式必然会对自适应光梳光谱造成干扰,降低系统的测量精度。In the current adaptive dual-comb spectroscopy system, the source of the compensation signal comes from the mutual optical beat frequency of two narrow linewidth continuous lasers and two ultrashort pulse lasers. Due to environmental temperature and vibration changes, each narrow linewidth laser or The optical frequency of an ultrashort pulse laser will drift randomly with time. Therefore, the center frequency of the optical beat signal will also drift randomly with time, and the drift range is very large, which can easily exceed the bandwidth of the circuit filter, resulting in the loss of the beat signal. The circuit processing module fails, and the adaptive compensation system works abnormally. In order to control the drift range of optical beat frequency signals of narrow-linewidth CW lasers and ultrashort pulse lasers, a negative feedback loop can be established, and parameters such as the temperature of the CW laser, the pump power, or the repetition frequency of the pulsed laser can be adjusted to combine the two The beat frequency drift of the system is controlled within the bandwidth of the follow-up circuit filter, but this negative feedback compensation method will inevitably cause interference to the adaptive optical comb spectrum and reduce the measurement accuracy of the system.
现有技术的自适应双光梳光谱系统,其补偿信号提取方式必须首先使用窄线宽连续激光器与超短脉冲激光器光学拍频,存在拍频信号丢失的重大隐患,影响系统长期运转的稳定性、可靠性。In the adaptive dual-comb spectroscopy system of the prior art, the compensation signal extraction method must first use a narrow linewidth continuous laser and an ultrashort pulse laser to optically beat the frequency, and there is a major hidden danger of the loss of the beat frequency signal, which affects the stability of the long-term operation of the system ,reliability.
发明内容Contents of the invention
本发明的目的是针对现有技术的不足而提供的一种自适应双光梳光谱补偿信号提取方法,采用两个独立的迈克尔逊干涉仪分别测量两台超短脉冲激光器的频率抖动,利用时域延时和频域选择的方式,提取出表征两台脉冲激光器相对重复频率抖动和相对载波包络相位抖动的信号,实现双光梳光谱测量,无需对脉冲激光器进行重复频率或者载波包络相位的主动控制,避免了传统自适应双光梳光谱补偿信号提取方式中窄线宽连续光与超短脉冲激光的光频飘移,保证拍频信号长期稳定存在,较好的解决了光学拍频信号飘移导致的补偿信号丢失问题,系统简单,操作方便,大大提高了自适应双光梳光谱测量系统的稳定性。The purpose of the present invention is to provide a kind of self-adaptive dual-comb spectral compensation signal extraction method for the deficiencies in the prior art, adopt two independent Michelson interferometers to measure the frequency jitter of two ultrashort pulse lasers respectively, utilize time The method of domain delay and frequency domain selection extracts the signals representing the relative repetition frequency jitter and relative carrier envelope phase jitter of two pulsed lasers, and realizes the dual optical comb spectrum measurement without the repetition frequency or carrier envelope phase of the pulsed laser. Active control avoids the optical frequency drift of narrow-linewidth continuous light and ultrashort pulse laser in the traditional adaptive dual-comb spectral compensation signal extraction method, ensures the long-term stable existence of beat frequency signals, and better solves the problem of optical beat frequency signals. The problem of compensation signal loss caused by drift is simple, the operation is convenient, and the stability of the adaptive dual-comb spectral measurement system is greatly improved.
本发明的目的是这样实现的:一种自适应双光梳光谱补偿信号提取方法,包括超短脉冲激光器、光电探测器、电学倍频及滤波单元和双光梳光谱数据处理单元,其特点是采用两个独立的迈克尔逊干涉仪分别测量两台超短脉冲激光器的频率抖动,然后将两迈克尔逊干涉仪的输出光分别进入两光学分束器,两光学分束器将输出光分别耦合到四个光学滤波器,滤波后的四个光信号分别由四个光电探测器进行探测,每个光电探测器将探测的光信号平均分为两路输出,一路输出光经电学倍频及滤波单元后与另一路输出光依次分别进入一级混频和二级混频,然后将二级混频处理后的光信号作为自适应双光梳光谱的补偿信号接入双光梳光谱数据处理单元;所述每个光电探测器将探测到的光信号平均分为两路输出为四个光电探测器分别探测得到①、②、③和④四个光信号,每个光信号平均分为两路输出,两光电探测器一路输出的①和②两个光信号分别经两个电学倍频及滤波单元后进入一级混频,两光电探测器另一路输出的①和②两个光信号分别直接进入一级混频;另两光电探测器一路输出的③和④两个光信号分别经另两个电学倍频及滤波单元后进入另一个一级混频,另两光电探测器另一路输出的③和④两个光信号分别直接进入另一个一级混频;所述四个电学倍频及滤波单元中的两个电学倍频及滤波单元分别对①和③光信号作(p-1)倍的倍频处理,另两个电学倍频及滤波单元分别对②和④光信号作p倍的倍频处理,分别得到①×(p-1)、②×p、③×(p-1)和④×p四个倍频后的信号,将倍频后的①×(p-1)、②×p和未倍频的①和②光信号进入一级混频,一级混频将①×(p-1)与①、②×p与②的光信号混频后,分别得到①-②和②×p-①×(p-1)两个混频的光信号;将倍频后的③×(p-1)、④×p和未倍频的③、④光信号进入另一个一级混频,另一个一级混频将③×(p-1)与③、④×p与④的光信号混频后,分别得到③-④和④×p-③×(p-1)两个混频的光信号,然后将①-②、②×p-①×(p-1)两个混频的光信号和③-④和④×p-③×(p-1)两个混频的光信号分别进入二级混频,经二级混频处理后分别得到(①-②)-(③-④)和[②×p-①×(p-1)]-[④×p-③×(p-1)]两个作为自适应双光梳光谱的补偿信号接入双光梳光谱数据处理单元;所述两个自适应双光梳光谱的补偿信号分别为(n-m)△fr和△f0+q△fr,其中:△fr为两台超短脉冲激光器的相对重复频率抖动;△f0为两台超短脉冲激光器的相对载波包络相位抖动;p、n和m均为正整数。The object of the present invention is achieved like this: a kind of self-adaptive dual optical comb spectral compensation signal extraction method comprises ultrashort pulse laser, photodetector, electrical frequency doubling and filtering unit and dual optical comb spectral data processing unit, and its characteristic is Two independent Michelson interferometers are used to measure the frequency jitter of two ultrashort pulse lasers respectively, and then the output light of the two Michelson interferometers enters two optical beam splitters respectively, and the two optical beam splitters couple the output light to the Four optical filters, the filtered four optical signals are respectively detected by four photodetectors, each photodetector divides the detected optical signal into two outputs on average, and one output light passes through the electrical frequency doubling and filtering unit Afterwards, the output light of the other channel enters the first-stage frequency mixing and the second-stage frequency mixing in turn, and then the optical signal processed by the second-stage frequency mixing is connected to the dual-comb spectrum data processing unit as a compensation signal for adaptive dual-comb spectrum; Each of the photodetectors divides the detected optical signal into two outputs on average, and four photodetectors respectively detect and obtain four optical signals ①, ②, ③ and ④, and each optical signal is divided into two outputs on average. , the two optical signals ① and ② output by one of the two photodetectors respectively enter the first-stage frequency mixing after passing through two electrical frequency multiplication and filtering units, and the two optical signals ① and ② output by the other of the two photodetectors directly enter the First-stage frequency mixing; the two optical signals ③ and ④ output by the other two photodetectors pass through the other two electrical frequency multiplication and filtering units respectively and then enter another first-stage frequency mixing, and the other two optical signals output by the other two photodetectors ③ and ④ two optical signals directly enter another first-level frequency mixing respectively; two electrical frequency multiplication and filtering units in the four electrical frequency multiplication and filtering units perform (p-1) times on ① and ③ optical signals respectively The other two electrical frequency multiplication and filtering units perform p-fold frequency multiplication processing on the ② and ④ optical signals respectively, and respectively obtain ①×(p-1), ②×p, ③×(p-1) and ④×p four frequency-multiplied signals, put the frequency-multiplied ①×(p-1), ②×p and the unmultiplied ① and ② optical signals into the first-level mixing, and the first-level mixing will ① After ×(p-1) is mixed with the optical signals of ①, ②×p and ②, two mixed optical signals of ①-② and ②×p-①×(p-1) are obtained respectively; The ③×(p-1), ④×p and the unmultiplied ③, ④ optical signals enter another level of mixing, and another level of mixing combines ③×(p-1) with ③, ④×p After mixing with the optical signal of ④, two mixed optical signals of ③-④ and ④×p-③×(p-1) are obtained respectively, and then ①-②, ②×p-①×(p-1 ) two mixed optical signals and ③-④ and ④×p-③×(p-1) two mixed optical signals respectively enter the secondary frequency mixing, and after the secondary frequency mixing processing, respectively obtain (①- ②)-(③-④) and [②×p-①×(p-1)]-[④×p-③×(p-1)] are connected as compensation signals for adaptive dual-comb spectrum Dual-comb spectrum data processing unit; the compensation signals of the two adaptive dual-comb spectra are (nm) △ f r and △ f 0 +q △ f r , wherein: △ f r is two ultrashort pulse Δf 0 is the relative carrier envelope phase jitter of two ultrashort pulse lasers; p, n and m are all positive integers.
所述迈克尔逊干涉仪包括分束片或光纤耦合器、两个法拉第反射镜、延时晶体或延时光纤和带驱动的声光调制器,超短脉冲激光器的输出光经过分束片或光纤耦合器按功率比为1:1分成两束光,一束光直接经过第一法拉第反射镜反射,回到分束片或光纤耦合器,另一束光先经过延时晶体或延时光纤,再经过带驱动的声光调制器,最后经过第二法拉第反射镜反射,依次经带驱动的声光调制器和延时晶体或延时光纤回到分束片或光纤耦合器,两束反射光在分束片或光纤耦合器上合成一束从迈克尔逊干涉仪输出。The Michelson interferometer includes a beam splitter or a fiber coupler, two Faraday mirrors, a time-delay crystal or a time-delay fiber and a band-driven acousto-optic modulator, and the output light of the ultrashort pulse laser passes through the beam splitter or fiber The coupler is divided into two beams of light according to the power ratio of 1:1. One beam of light is directly reflected by the first Faraday reflector and returns to the beam splitter or fiber coupler. The other beam of light first passes through the delay crystal or delay fiber. After passing through the belt-driven acousto-optic modulator, and finally reflected by the second Faraday reflector, the belt-driven acousto-optic modulator and time-delay crystal or time-delay fiber return to the beam splitter or fiber coupler in turn, and the two reflected lights Combine a beam on the beam splitter or fiber coupler to output from the Michelson interferometer.
所述光学分束器为1:1分光比的半透半反镜片或光纤耦合器。The optical beam splitter is a semi-transparent and semi-reflective mirror with a splitting ratio of 1:1 or a fiber optic coupler.
所述光学滤波器为允许特定波长的激光通过,隔绝其他波长激光通过的窄带滤波镜片、光纤光栅或光纤滤波器。The optical filter is a narrow-band filter mirror, optical fiber grating or optical fiber filter that allows laser light of a specific wavelength to pass through and prevents laser light of other wavelengths from passing through.
本发明与现有技术相比具有以下优点:Compared with the prior art, the present invention has the following advantages:
⑴、不需要窄线宽连续激光器以及连续激光器和脉冲激光器的光学拍频,因此避免了光学拍频信号的随机抖动,提高自适应双光梳光谱补偿系统的稳定性和可靠性。(1) The optical beat frequency of narrow-linewidth continuous lasers and continuous lasers and pulsed lasers is not required, so random jitter of optical beat frequency signals is avoided, and the stability and reliability of the adaptive dual optical comb spectral compensation system are improved.
⑵、采用两个独立的迈克尔逊干涉仪分别测量两台超短脉冲激光器的频率抖动,彼此之间无干扰,探测噪声低。⑵. Two independent Michelson interferometers are used to measure the frequency jitter of two ultrashort pulse lasers respectively, without interference between each other and low detection noise.
⑶、采用频域选择的方式,在两个不同的频率窗口相干探测脉冲的拍频信号,通过灵活选择频率间隔和窗口大小,能降低频域噪声,提高探测灵敏度。⑶. Using the frequency domain selection method, coherently detect the beat frequency signal of the pulse in two different frequency windows. By flexibly selecting the frequency interval and window size, the frequency domain noise can be reduced and the detection sensitivity can be improved.
⑷、迈克尔逊干涉仪中加入带驱动的声光调制器,将拍频信号的中心频率偏移到声光调制器的中心频率,降低零频附近的噪声。⑷. An acousto-optic modulator with a drive is added to the Michelson interferometer to shift the center frequency of the beat frequency signal to the center frequency of the acousto-optic modulator to reduce the noise near zero frequency.
⑸、无需额外建立负反馈电路补偿信号的飘移,不会引入电子线路噪声,能实现高精度的自适应双光梳光谱。(5) There is no need to establish an additional negative feedback circuit to compensate the drift of the signal, and no electronic circuit noise will be introduced, and high-precision adaptive dual-comb spectroscopy can be realized.
⑹、采用自由运转的两台脉冲激光器和自适应控制技术就能实现双光梳光谱测量,无需对脉冲激光器进行重复频率或者载波包络相位的主动控制,系统简单,操作方便。⑹. Using two free-running pulse lasers and adaptive control technology can realize dual-comb spectrum measurement without active control of repetition frequency or carrier envelope phase of pulse lasers. The system is simple and easy to operate.
附图说明Description of drawings
图1为本发明结构示意图;Fig. 1 is a structural representation of the present invention;
图2为空间结构的实施例结构示意图;Fig. 2 is the structural representation of the embodiment of spatial structure;
图3为光纤结构的实施例结构示意图;Fig. 3 is the structural representation of the embodiment of optical fiber structure;
图4为本发明具体运用示意图。Fig. 4 is a schematic diagram of the specific application of the present invention.
具体实施方式detailed description
参阅附图1,本发明由两个超短脉冲激光器1、两个迈克尔逊干涉仪2、两个光学分束器3、四个光学滤波器4、四个光电探测器5、四个电学倍频及滤波单元6、两个一级混频7、一个二级混频8和双光梳光谱数据处理单元9组成两补偿信号的光路,所述一超短脉冲激光器1的输出光入射到迈克尔逊干涉仪2,迈克尔逊干涉仪2中光程较长的干涉臂包含一个带驱动的声光调制器,调制频率为fa,迈克尔逊干涉仪2的输出光经过一光学分束器3,分别耦合到两光学滤波器4,采用两个光电探测器5分别测量滤波后的光学信号,探测到的信号分别记为①和②,①信号可以表示为2fa+△(nfr1+f01+2fa),n为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器1的重复频率飘移△fr1、载波包络相位飘移△f01和声光调制器驱动频率的飘移△fa,②信号可以表示为2fa+△(mfr1+f01+2fa),m为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器1的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移;另一超短脉冲激光器1的输出光入射到另一迈克尔逊干涉仪2,另一迈克尔逊干涉仪2中光程较长的干涉臂包含一个带驱动的声光调制器,调制频率为fb,另一迈克尔逊干涉仪2的输出光经另一光学分束器2,分别耦合到另两个光学滤波器4,采用另两个光电探测器4分别测量滤波后的光学信号,探测到的信号分别记为③和④,③信号可以表示为2fb+△(nfr2+f02+2fb),n为正整数,其中心频率在2fb,信号的飘移包括另一超短脉冲激光器1的重复频率飘移△fr2、载波包络相位飘移△f02和声光调制器驱动频率的飘移△fb,④信号可以表示为2fb+△(mfr2+f02+2fb),m为正整数,其中心频率在2fb,信号的飘移包括另一超短脉冲激光器1的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移;每个光电探测器5探测到的信号均平均分为两路,一路不倍频,一路采用电学倍频及滤波单元6分别对上述①、②、③和④信号作倍频处理,对①和③信号作(p-1)倍的倍频处理,p为正整数,对②和④信号作p倍的倍频处理,四个电学倍频及滤波单元6得到四个倍频后的信号,分别为①×(p-1)、②×p、③×(p-1)和④×p;采用两个一级混频7,对倍频后的信号及不倍频的信号作第一级混频处理,得到四个第一级混频后的信号,分别为:①-②=(n-m)△fr1、②×p-①×(p-1)=2fa+△f01+(pgm-pgn+n)△fr1+2△fa、③-④=(n-m)△fr2、④×p-③×(p-1)=2fb+△f02+(pgm-pgn+n)△fr2+2△fb;采用二级混频8对第一级混频后的信号作第二级混频处理,得到两个第二级混频后的信号,分别为:(①-②)-③-④)=(n-m)△fr、[②×p-①×(p-1)]-[④×p-③×(p-1)]=2(fa-fb)+△f0+(pgm-pgn+n)△fr+2△(fa-fb);△f0+q△fr,上述两个信号实质上分别等于(n-m)△fr,△f0+q△fr,其中n、m和q为正整数,△fr为两台超短脉冲激光器1的相对重复频率抖动,△f0为两台超短脉冲激光器1的相对载波包络相位抖动,这两个信号就可以直接作为自适应双光梳光谱的补偿信号。其中,迈克尔逊干涉仪2包括分束片或光纤耦合器、第一法拉第反射镜、第二法拉第反射镜、延时晶体或者延时光纤和带驱动的声光调制器。所述脉冲激光器1的输出光经过分束片或光纤耦合器,按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜反射回到分束片或光纤耦合器,另一束光先经过延时晶体或者延时光纤,经带驱动的声光调制器和第二法拉第反射镜反射,再依次经过带驱动的声光调制器和延时晶体或者延时光纤,回到分束片或光纤耦合器,两束反射光在分束片或光纤耦合器上合成一束从迈克尔逊干涉仪2输出。Referring to accompanying drawing 1, the present invention is made of two ultrashort pulse lasers 1, two Michelson interferometers 2, two optical beam splitters 3, four optical filters 4, four photodetectors 5, four electrical times Frequency and filter unit 6, two first-stage frequency mixers 7, one second-stage frequency mixer 8 and dual optical comb spectral data processing unit 9 form the optical path of two compensation signals, and the output light of the ultrashort pulse laser 1 is incident on Michael The Michelson interferometer 2, the longer interferometric arm of the Michelson interferometer 2 includes a band-driven acousto-optic modulator, the modulation frequency is fa, and the output light of the Michelson interferometer 2 passes through an optical beam splitter 3, respectively Coupled to two optical filters 4, two photodetectors 5 are used to measure the filtered optical signals respectively, and the detected signals are respectively recorded as ① and ②, and the signal ① can be expressed as 2f a + △ (nf r1 +f 01 + 2f a ), n is a positive integer, its center frequency is at 2fa, and the signal drift includes the repetition frequency drift △f r1 of the ultrashort pulse laser 1, the carrier envelope phase drift △f 01 and the driving frequency drift of the acousto-optic modulator △ f a , ②The signal can be expressed as 2f a +△(mf r1 +f 01 +2f a ), m is a positive integer, its center frequency is at 2fa, and the signal drift includes repetition frequency drift of ultrashort pulse laser 1, carrier packet network phase drift and the drift of the drive frequency of the AOM; the output light of another ultrashort pulse laser 1 is incident on another Michelson interferometer 2, and the interference arm with a longer optical path in the other Michelson interferometer 2 contains a Acousto-optic modulator with drive, the modulation frequency is fb, the output light of another Michelson interferometer 2 is coupled to the other two optical filters 4 through another optical beam splitter 2, and the other two photodetectors are used 4Measure the filtered optical signals respectively, and the detected signals are recorded as ③ and ④ respectively. The ③ signal can be expressed as 2f b + △(nf r2 + f 02 + 2f b ), n is a positive integer, and its center frequency is at 2fb , the drift of the signal includes the repetition frequency drift △f r2 of another ultrashort pulse laser 1, the carrier envelope phase drift △f 02 and the driving frequency drift △f b of the acousto-optic modulator, ④The signal can be expressed as 2f b + △ (mf r2 +f 02 +2f b ), m is a positive integer, its center frequency is at 2fb, and the drift of the signal includes the repetition frequency drift of another ultrashort pulse laser 1, the carrier envelope phase drift and the driving frequency of the acousto-optic modulator drift; the signal detected by each photodetector 5 is evenly divided into two paths, one path is not frequency multiplied, and the other path adopts electrical frequency multiplication and filtering unit 6 to perform frequency multiplication processing on the above-mentioned ①, ②, ③ and ④ signals respectively, Perform (p-1) times frequency multiplication processing on ① and ③ signals, p is a positive integer, and perform p times frequency multiplication processing on ② and ④ signals, after four electrical frequency multiplication and filtering units 6 obtain four frequency multiplications The signals are respectively ①×(p-1), ②×p, ③×(p-1) and ④×p; using two first-stage mixers 7, for the multiplied signal and the non-multiplied signal mix first frequency processing to obtain four first-stage mixed signals, respectively: ①-②=(nm)△f r1 , ②×p-①×(p-1)=2f a +△f 01 +(pgm -pgn+n)△f r1 +2△f a 、③-④=(nm)△f r2 、④×p-③×(p-1)=2f b +△f 02 +(pgm-pgn+n )△f r2 +2△f b ; use the second-stage mixing 8 to perform the second-stage mixing processing on the signal after the first-stage mixing, and obtain two signals after the second-stage mixing, respectively: (① -②)-③-④)=(nm)△f r 、[②×p-①×(p-1)]-[④×p-③×(p-1)]=2(f a -f b )+△f 0 +(pgm-pgn+n)△f r +2△(f a -f b );△f 0 +q△f r , the above two signals are substantially equal to (nm)△f r , △f 0 +q△f r , where n, m and q are positive integers, △fr r is the relative repetition frequency jitter of the two ultrashort pulse lasers 1, △f 0 is the jitter of the two ultrashort pulse lasers 1 Relative to the phase jitter of the carrier envelope, these two signals can be directly used as compensation signals for adaptive dual-comb spectroscopy. Wherein, the Michelson interferometer 2 includes a beam splitter or a fiber coupler, a first Faraday mirror, a second Faraday mirror, a time-delay crystal or a time-delay fiber, and a tape-driven acousto-optic modulator. The output light of the pulsed laser 1 passes through a beam splitter or a fiber coupler, and is divided into two beams of light at a power ratio of 1:1, and one beam of light is directly reflected back to the beam splitter or fiber coupler through the first Faraday reflector, The other beam of light first passes through the delay crystal or delay fiber, is reflected by the belt-driven acousto-optic modulator and the second Faraday mirror, and then passes through the belt-driven acousto-optic modulator and time-delay crystal or delay fiber in turn, and returns to the To the beam splitter or fiber coupler, the two beams of reflected light are synthesized on the beam splitter or fiber coupler and output from the Michelson interferometer 2.
下面以空间结构和光纤结构的自适应双光梳光谱补偿信号提取系统对本发明作进一步说明。The present invention will be further described below with an adaptive dual-comb spectral compensation signal extraction system with a spatial structure and an optical fiber structure.
实施例1Example 1
参阅附图2,本发明由两个脉冲激光器11、12、两个迈克尔逊干涉仪21、22、两个光学分束器31、32、四个光学滤波器41、42、43、44、四个光电探测器51、52、53、54、四个电学倍频及滤波单元61、62、63、64、两个一级混频器71、72、一个二级混频器8和双光梳光谱数据处理单元9组成空间结构的自适应双光梳光谱补偿信号的提取系统。Referring to accompanying drawing 2, the present invention is by two pulsed lasers 11,12, two Michelson interferometers 21,22, two optical beam splitters 31,32, four optical filters 41,42,43,44, four photodetectors 51, 52, 53, 54, four electrical frequency doubling and filtering units 61, 62, 63, 64, two primary mixers 71, 72, a secondary mixer 8 and dual optical combs The spectral data processing unit 9 constitutes a space-structured adaptive dual-comb spectral compensation signal extraction system.
所述超短脉冲激光器11的输出光入射到迈克尔逊干涉仪21,其中,迈克尔逊干涉仪21由分束片211、第一法拉第反射镜212、第二法拉第反射镜215、延时晶体213和带驱动的声光调制器214组成,调制频率为fa。超短脉冲激光器11的输出光经过分束片211按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜212反射回到分束片211,另一束光先经过延时晶体213和带驱动的声光调制器214,最后经过第二法拉第反射镜213反射,再依次经过带驱动的声光调制器214和延时晶体213后回到分束片211,两束反射光在分束片211上合成一束从迈克尔逊干涉仪21输出。The output light of the ultrashort pulse laser 11 is incident on the Michelson interferometer 21, wherein the Michelson interferometer 21 is composed of a beam splitter 211, a first Faraday reflector 212, a second Faraday reflector 215, a time delay crystal 213 and It is composed of an acousto-optic modulator 214 driven by a belt, and the modulation frequency is fa. The output light of the ultrashort pulse laser 11 passes through the beam splitter 211 and is divided into two beams of light at a power ratio of 1:1. One beam of light is directly reflected back to the beam splitter 211 by the first Faraday reflector 212, and the other beam passes through the The time-delay crystal 213 and the AOM 214 driven by the belt are finally reflected by the second Faraday reflector 213, and then return to the beam splitter 211 after passing through the AOM 214 and the time-delay crystal 213 driven by the belt in turn. The reflected light is combined into one beam on the beam splitter 211 and output from the Michelson interferometer 21 .
所述另一超短脉冲激光器12的输出光入射到另一迈克尔逊干涉仪22,其中,另一迈克尔逊干涉仪22由分束片221、第一法拉第反射镜222、第二法拉第反射镜225、延时晶体223和带驱动的声光调制器224组成,调制频率为fa。另一超短脉冲激光器12的输出光经过分束片221按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜222反射回到分束片221,另一束光先经过延时晶体223和带驱动的声光调制器224,最后经过第二法拉第反射镜225反射,再依次经过带驱动的声光调制器224和延时晶体223后回到分束片221,两束反射光在分束片221上合成一束从另一迈克尔逊干涉仪22输出。The output light of another ultrashort pulse laser 12 is incident on another Michelson interferometer 22, wherein, another Michelson interferometer 22 is composed of a beam splitter 221, a first Faraday reflector 222, a second Faraday reflector 225 , time-delay crystal 223 and band-driven acousto-optic modulator 224, and the modulation frequency is fa. The output light of another ultrashort pulse laser 12 is divided into two beams of light by a power ratio of 1:1 through the beam splitter 221, and one beam of light is directly reflected back to the beam splitter 221 through the first Faraday reflector 222, and the other beam of light First pass through the time-delay crystal 223 and the belt-driven AOM 224, and finally reflect through the second Faraday reflector 225, then pass through the belt-driven AOM 224 and the time-delay crystal 223, and then return to the beam splitter 221, The two beams of reflected light are synthesized on the beam splitter 221 and output from another Michelson interferometer 22 .
所述迈克尔逊干涉仪21的输出光经过光学分束器31分别耦合到光学滤波器41和另一光学滤波器43,采用两个光电探测器51、53分别测量滤波后的光学信号,探测到的信号分别记为①、②,信号①可以表示为2fa+△(nfr1+f01+2fa),n为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移△fr1、载波包络相位飘移△f01和声光调制器驱动频率的飘移△fa,信号②可以表示为2fa+△(mfr1+f01+2fa),m为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the Michelson interferometer 21 is respectively coupled to an optical filter 41 and another optical filter 43 through an optical beam splitter 31, and two photodetectors 51, 53 are used to measure the filtered optical signals respectively, and detect The signals are denoted as ① and ② respectively, the signal ① can be expressed as 2f a + △ (nf r1 + f 01 + 2f a ), n is a positive integer, its center frequency is 2fa, and the drift of the signal includes the ultrashort pulse laser 11 Repetition frequency drift △f r1 , carrier envelope phase drift △f 01 and AOM driving frequency drift △f a , signal ② can be expressed as 2f a + △(mf r1 +f 01 +2f a ), m is A positive integer whose center frequency is 2fa, and the signal drift includes the repetition frequency drift of the ultrashort pulse laser 11, the carrier envelope phase drift and the driving frequency drift of the acousto-optic modulator.
所述另一迈克尔逊干涉仪22的输出光经过光学分束器32分别耦合到光学滤波器42和另一光学滤波器44,采用两个光电探测器52、54分别测量滤波后的光学信号,探测到的信号分别记为③、④,信号③可以表示为2fb+△(nfr2+f02+2fb),n为正整数,其中心频率在2fb,信号的飘移包括超短脉冲激光器21的重复频率飘移△fr2、载波包络相位飘移△f02和声光调制器驱动频率的飘移△fb,信号④可以表示为2fb+△(mfr2+f02+2fb),m为正整数,其中心频率在2fb,信号的飘移包括另一超短脉冲激光器12的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the other Michelson interferometer 22 is respectively coupled to the optical filter 42 and another optical filter 44 through the optical beam splitter 32, and two photodetectors 52, 54 are used to measure the filtered optical signals respectively, The detected signals are denoted as ③ and ④ respectively, and the signal ③ can be expressed as 2f b + △ (nf r2 + f 02 + 2f b ), n is a positive integer, and its center frequency is at 2fb. The drift of the signal includes the ultrashort pulse laser 21 repetition frequency drift △f r2 , carrier envelope phase drift △f 02 and AOM drive frequency drift △f b , the signal ④ can be expressed as 2f b + △(mf r2 +f 02 +2f b ), m is a positive integer, its center frequency is 2fb, and the signal drift includes the repetition frequency drift of another ultrashort pulse laser 12, the carrier envelope phase drift and the drive frequency drift of the AOM.
将上述四个光电探测器51、52、53、54探测到的①、②、③和④信号按每个信号平均分为两路,一路不倍频,一路采用四个电学倍频及滤波单元61、62、63、64分别对上述①、②、③和④信号作倍频处理,对信号①和③作(p-1)倍的倍频处理,p为正整数,对信号②和④作p倍的倍频处理,得到四个倍频后的信号,分别为①×(p-1)、②×p、③×(p-1)、④×p;采用两个一级混频器71、72对倍频后的信号及不倍频的信号作第一级混频处理,得到四个第一级混频后的信号,分别为:①-②=(n-m)△fr1、②×p-①×(p-1)=2fa+△f01+(pgm-pgn+n)△fr1+2△fa、③-④=(n-m)△fr2、④×p-③×(p-1)=2fb+△f02+(pgm-pgn+n)△fr2+2△fb;采用一个二级混频器8,对第一级混频后的信号作第二级混频处理,得到两个第二级混频后的信号,分别为:(①-②)-(③-④)=(n-m)△fr、[②×p-①×(p-1)]-[④×p-③×(p-1)]=2(fa-fb)+△f0+(pgm-pgn+n)△fr+2△(fa-fb);△f0+q△fr,上述两个信号实质上分别等于(n-m)△fr,△f0+q△fr,其中n、m、q为正整数,△fr为两台超短脉冲激光器11、12的相对重复频率抖动,△f0为两台超短脉冲激光器11、12的相对载波包络相位抖动,这两个信号就可以直接作为自适应双光梳光谱的补偿信号接入双光梳光谱数据处理单元9。Divide the ①, ②, ③ and ④ signals detected by the above four photodetectors 51, 52, 53, 54 into two paths on average for each signal, one path without frequency multiplication, and one path with four electrical frequency multiplication and filtering units 61, 62, 63, and 64 perform frequency multiplication processing on the above-mentioned ①, ②, ③ and ④ signals respectively, and perform (p-1) times frequency multiplication processing on the signals ① and ③, p is a positive integer, and perform frequency multiplication processing on the signals ② and ④ Perform p-fold frequency multiplication processing to obtain four multiplied signals, which are ①×(p-1), ②×p, ③×(p-1), ④×p; using two first-level mixing Devices 71 and 72 perform first-stage mixing processing on the multiplied signal and the non-multiplied signal to obtain four first-stage mixed signals, which are respectively: ①-②=(nm)△f r1 , ②×p-①×(p-1)=2f a +△f 01 +(pgm-pgn+n)△f r1 +2△f a 、③-④=(nm)△f r2 、④×p- ③×(p-1)=2f b +△f 02 +(pgm-pgn+n)△f r2 +2△f b ; a second-stage mixer 8 is used to make the first-stage mixed signal The second-stage frequency mixing process obtains two second-stage mixed signals, which are: (①-②)-(③-④)=(nm)△f r , [②×p-①×(p -1)]-[④×p-③×(p-1)]=2(f a -f b )+△f 0 +(pgm-pgn+n)△f r +2△(f a -f b ); △f 0 +q△fr , the above two signals are substantially equal to (nm)△f r , △f 0 +q△f r , where n, m, and q are positive integers, and △f r is The relative repetition frequency jitter of the two ultrashort pulse lasers 11, 12, △ f 0 is the relative carrier envelope phase jitter of the two ultrashort pulse lasers 11, 12, these two signals can be directly used as adaptive dual optical comb spectrum The compensation signal is connected to the dual-comb spectral data processing unit 9.
实施例2Example 2
参阅附图3,本发明由两个脉冲激光器11、12、两个迈克尔逊干涉仪21、22、两个光学分束器31、32、四个光学滤波器41、42、43、45、四个光电探测器51、52、53、54、四个电学倍频及滤波单元61、62、63、64、两个一级混频器71、72和一个二级混频器8和双光梳光谱数据处理单元9组成光纤结构的自适应双光梳光谱补偿信号的提取系统。Referring to accompanying drawing 3, the present invention is by two pulsed lasers 11,12, two Michelson interferometers 21,22, two optical beam splitters 31,32, four optical filters 41,42,43,45, four A photodetector 51, 52, 53, 54, four electrical frequency doubling and filtering units 61, 62, 63, 64, two primary mixers 71, 72 and a secondary mixer 8 and dual optical combs The spectral data processing unit 9 constitutes an extraction system of an adaptive dual-comb spectral compensation signal with a fiber structure.
所述超短脉冲激光器11的输出光入射到迈克尔逊干涉仪21,其中,迈克尔逊干涉仪由光纤耦合器211、第一法拉第反射镜212、第二法拉第反射镜215、延时光纤213和带驱动的声光调制器214组成,调制频率为fa。脉冲激光器11的输出光经过光纤耦合器211按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜212反射回到光纤耦合器211,另一束光先经过延时光纤213,再经过带驱动的声光调制器214,最后经过第二法拉第反射镜215反射,再依次经过带驱动的声光调制器214和延时光纤213回到光纤耦合器211,两束反射光在光纤耦合器211上合成一束从迈克尔逊干涉仪21输出。The output light of the ultrashort pulse laser 11 is incident on the Michelson interferometer 21, wherein the Michelson interferometer consists of a fiber coupler 211, a first Faraday reflector 212, a second Faraday reflector 215, a time-delay fiber 213 and a ribbon The driven acousto-optic modulator 214 is composed of the modulation frequency fa. The output light of the pulsed laser 11 passes through the fiber coupler 211 and is divided into two beams of light according to a power ratio of 1:1. One beam of light is directly reflected back to the fiber coupler 211 by the first Faraday reflector 212, and the other beam of light is first delayed. The optical fiber 213 passes through the acousto-optic modulator 214 driven by the belt, and finally reflects through the second Faraday reflector 215, and then returns to the fiber coupler 211 through the acousto-optic modulator 214 driven by the belt and the time-delay optical fiber 213 in turn. The light is combined into one beam on the fiber coupler 211 and output from the Michelson interferometer 21 .
所述迈克尔逊干涉仪21的输出光经过光学分束器31分别耦合到光学滤波器41和光学滤波器43,采用两个光电探测器51、53分别测量滤波后的光学信号,探测到的信号分别记为①、②,信号①可以表示为2fa+△(nfr1+f01+2fa),n为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移△fr1、载波包络相位飘移△f01和声光调制器驱动频率的飘移△fa,信号②可以表示为2fa+△(mfr1+f01+2fa),m为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the Michelson interferometer 21 is respectively coupled to the optical filter 41 and the optical filter 43 through the optical beam splitter 31, and two photodetectors 51, 53 are used to measure the filtered optical signal respectively, and the detected signal Denoted as ① and ② respectively, the signal ① can be expressed as 2f a + △ (nf r1 + f 01 + 2f a ), n is a positive integer, its center frequency is 2fa, and the drift of the signal includes the repetition frequency of the ultrashort pulse laser 11 Drift △f r1 , carrier envelope phase drift △f 01 and AOM driving frequency drift △f a , signal ② can be expressed as 2f a + △(mf r1 +f 01 +2f a ), m is a positive integer , its center frequency is at 2fa, and the drift of the signal includes the drift of the repetition frequency of the ultrashort pulse laser 11, the drift of the carrier envelope phase and the drift of the driving frequency of the acousto-optic modulator.
所述另一超短脉冲激光器12的输出光入射到迈克尔逊干涉仪22,其中,迈克尔逊干涉仪22由光纤耦合器221、第一法拉第反射镜222、第二法拉第反射镜225、延时光纤223和带驱动的声光调制器224组成,调制频率为fb。另一超短脉冲激光器12的输出光经过光纤耦合器221按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜222反射回到光纤耦合器221,另一束光先经过延时光纤223,再经过带驱动的声光调制器224,最后经过第二法拉第反射镜225反射,再依次经过带驱动的声光调制器224和延时光纤223回到光纤耦合器221,两束反射光在光纤耦合器221上合成一束从迈克尔逊干涉仪22输出。The output light of another ultrashort pulse laser 12 is incident on the Michelson interferometer 22, wherein the Michelson interferometer 22 is composed of a fiber coupler 221, a first Faraday reflector 222, a second Faraday reflector 225, a time-delay optical fiber 223 and the band-driven acousto-optic modulator 224, the modulation frequency is fb. The output light of another ultrashort pulse laser 12 is divided into two beams of light by a power ratio of 1:1 through the fiber coupler 221, and one beam of light is directly reflected back to the fiber coupler 221 through the first Faraday reflector 222, and the other beam of light First pass through the time-delay fiber 223, then pass through the belt-driven acousto-optic modulator 224, and finally reflect through the second Faraday mirror 225, and then return to the fiber coupler 221 through the belt-driven acousto-optic modulator 224 and time-delay fiber 223 , the two beams of reflected light are synthesized on the fiber coupler 221 and output from the Michelson interferometer 22 .
所述另一迈克尔逊干涉仪22的输出光经过光学分束器32分别耦合到光学滤波器42和另一光学滤波器44,采用两个光电探测器52、54分别测量滤波后的光学信号,探测到的信号分别记为③、④,信号③可以表示为2fb+△(nfr2+f02+2fb),n为正整数,其中心频率在2fb,信号的飘移包括超短脉冲激光器12的重复频率飘移△fr2、载波包络相位飘移△f02和声光调制器驱动频率的飘移△fb,信号④可以表示为2fb+△(mfr2+f02+2fb),m为正整数,其中心频率在2fb,信号的飘移包括超短脉冲激光器12的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the other Michelson interferometer 22 is respectively coupled to the optical filter 42 and another optical filter 44 through the optical beam splitter 32, and two photodetectors 52, 54 are used to measure the filtered optical signals respectively, The detected signals are denoted as ③ and ④ respectively, and the signal ③ can be expressed as 2f b + △ (nf r2 + f 02 + 2f b ), n is a positive integer, and its center frequency is at 2fb. The drift of the signal includes the ultrashort pulse laser The repetition frequency drift △f r2 of 12, the carrier envelope phase drift △f 02 and the driving frequency drift △f b of the AOM, the signal ④ can be expressed as 2f b + △(mf r2 +f 02 +2f b ), m is a positive integer, and its center frequency is 2fb. The signal drift includes the repetition frequency drift of the ultrashort pulse laser 12, the carrier envelope phase drift and the driving frequency drift of the acousto-optic modulator.
将上述四个光电探测器51、52、53、54探测到的①、②、③和④信号按每个信号平均分为两路,一路不倍频,一路采用四个电学倍频及滤波单元61、2、63、64分别对上述①、②、③和④信号作倍频处理,对信号①和③作(p-1)倍的倍频处理,p为正整数,对信号②和④作p倍的倍频处理,得到四个倍频后的信号,分别为①×(p-1)、②×p、③×(p-1)、④×p;采用两个一级混频器71、72对倍频后的信号及不倍频的信号作第一级混频处理,得到四个第一级混频后的信号,分别为:①-②=(n-m)△fr1、②×p-①×(p-1)=2fa+△f01+(pgm-pgn+n)△fr1+2△fa、③-④=(n-m)△fr2、④×p-③×(p-1)=2fb+△f02+(pgm-pgn+n)△fr2+2△fb;采用一个二级混频器8对第一级混频后的信号作第二级混频处理,得到两个第二级混频后的信号,分别为:(①-②)-(③-④)=(n-m)△fr、[②×p-①×(p-1)]-[④×p-③×(p-1)]=2(fa-fb)+△f0+(pgm-pgn+n)△fr+2△(fa-fb);△f0+q△fr,上述两个信号实质上分别等于(n-m)△fr,△f0+q△fr,其中n、m、q为正整数,△fr为两台超短脉冲激光器11、21的相对重复频率抖动,△f0为两台超短脉冲激光器11、21的相对载波包络相位抖动,这两个信号就可以直接作为自适应双光梳光谱的补偿信号接入双光梳光谱数据处理单元9。Divide the ①, ②, ③ and ④ signals detected by the above four photodetectors 51, 52, 53, 54 into two paths on average for each signal, one path without frequency multiplication, and one path with four electrical frequency multiplication and filtering units 61, 2, 63, and 64 perform frequency multiplication processing on the above-mentioned ①, ②, ③ and ④ signals respectively, and perform (p-1) times frequency multiplication processing on the signals ① and ③, p is a positive integer, and perform frequency multiplication processing on the signals ② and ④ Perform p-fold frequency multiplication processing to obtain four multiplied signals, which are ①×(p-1), ②×p, ③×(p-1), ④×p; using two first-level mixing Devices 71 and 72 perform first-stage mixing processing on the multiplied signal and the non-multiplied signal to obtain four first-stage mixed signals, which are respectively: ①-②=(nm)△f r1 , ②×p-①×(p-1)=2f a +△f 01 +(pgm-pgn+n)△f r1 +2△f a 、③-④=(nm)△f r2 、④×p- ③×(p-1)=2f b +△f 02 +(pgm-pgn+n)△f r2 +2△f b ; adopt a second-stage mixer 8 to make the first-stage mixed signal Second-level mixing processing, to obtain two second-level mixing signals, respectively: (①-②)-(③-④)=(nm)△f r , [②×p-①×(p- 1)]-[④×p-③×(p-1)]=2(f a -f b )+△f 0 +(pgm-pgn+n)△f r +2△(f a -f b ); △f 0 +q△f r , the above two signals are substantially equal to (nm)△f r , △f 0 +q△f r , where n, m, q are positive integers, and △f r is two The relative repetition frequency jitter of two ultrashort pulse lasers 11, 21, Δf 0 is the relative carrier envelope phase jitter of two ultrashort pulse lasers 11, 21, these two signals can be directly used as the The compensation signal is connected to the dual-comb spectral data processing unit 9 .
实施例3Example 3
参阅附图4,本发明由两个脉冲激光器11、12、样品池10、干涉信号探测模块14、两个迈克尔逊干涉仪21、22、两个光学分束器31、32、四个光学滤波器41、42、43、44、四个光电探测器51、52、53、54、四个电学倍频及滤波单元61、62、63、64、两个一级混频器71、72、一个二级混频器8、自适应双光梳光谱系统数据处理单元9和光谱测量结果输出模块13组成。Referring to accompanying drawing 4, the present invention is made up of two pulsed lasers 11,12, sample pool 10, interference signal detection module 14, two Michelson interferometers 21,22, two optical beam splitters 31,32, four optical filter 41, 42, 43, 44, four photodetectors 51, 52, 53, 54, four electrical frequency multiplication and filtering units 61, 62, 63, 64, two primary mixers 71, 72, one It consists of a secondary mixer 8, an adaptive dual-comb spectral system data processing unit 9 and a spectral measurement result output module 13.
所述超短脉冲激光器11的输出分为两路,一路直接入射到样品池10,另一路入射到迈克尔逊干涉仪21,其中,迈克尔逊干涉仪由分束片211、第一法拉第反射镜212、第二法拉第反射镜215、延时晶体213和带驱动的声光调制器214组成,调制频率为fa。超短脉冲激光器11的输出光经过分束片211按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜212反射回到分束片211,另一束光先经过延时晶体213,再经过带驱动的声光调制器214,最后经过第二法拉第反射镜215反射,再依次经过带驱动的声光调制器214和延时晶体213回到分束片211,两束反射光在分束片211上合成一束,从迈克尔逊干涉仪21输出。The output of the ultrashort pulse laser 11 is divided into two paths, one path is directly incident to the sample cell 10, and the other path is incident to the Michelson interferometer 21, wherein the Michelson interferometer is composed of a beam splitter 211, a first Faraday reflector 212 , a second Faraday reflector 215, a delay crystal 213 and a band-driven acousto-optic modulator 214, and the modulation frequency is fa. The output light of the ultrashort pulse laser 11 passes through the beam splitter 211 and is divided into two beams of light at a power ratio of 1:1. One beam of light is directly reflected back to the beam splitter 211 by the first Faraday reflector 212, and the other beam passes through the The time-delay crystal 213 passes through the acousto-optic modulator 214 driven by the belt, and finally is reflected by the second Faraday reflector 215, and then returns to the beam splitter 211 through the acousto-optic modulator 214 and the time-delay crystal 213 driven by the belt in turn. The beams of reflected light are combined into one beam on the beam splitter 211 and output from the Michelson interferometer 21 .
所述另一超短脉冲激光器12的输出光入射到另一迈克尔逊干涉仪22,其中,另一迈克尔逊干涉仪22由分束片221、第一法拉第反射镜222、第二法拉第反射镜225、延时晶体223和带驱动的声光调制器224组成,调制频率为fb。另一超短脉冲激光器21的输出光经过分束片221按1:1的功率比分成两束光,一束光直接经过第一法拉第反射镜222反射回到分束片221,另一束光先经过延时晶体223,再经过带驱动的声光调制器224,最后经过第二法拉第反射镜225反射,再依次经过带驱动的声光调制器224和延时晶体223回到分束片221,两束反射光在分束片221上合成一束从另一迈克尔逊干涉仪22输出。The output light of another ultrashort pulse laser 12 is incident on another Michelson interferometer 22, wherein, another Michelson interferometer 22 is composed of a beam splitter 221, a first Faraday reflector 222, a second Faraday reflector 225 , time-delay crystal 223 and band-driven acousto-optic modulator 224, and the modulation frequency is fb. The output light of another ultrashort pulse laser 21 is divided into two beams of light by a power ratio of 1:1 through the beam splitter 221, and one beam of light is directly reflected back to the beam splitter 221 through the first Faraday reflector 222, and the other beam of light First pass through the delay crystal 223, then pass through the belt-driven acousto-optic modulator 224, and finally reflect through the second Faraday mirror 225, and then return to the beam splitter 221 through the belt-driven acousto-optic modulator 224 and time-delay crystal 223 , the two beams of reflected light are synthesized on the beam splitter 221 and output from another Michelson interferometer 22 .
所述迈克尔逊干涉仪21的输出光经过光学分束器31分别耦合到光学滤波器41和另一光学滤波器43,采用两个光电探测器51、53分别测量滤波后的光学信号,探测到的信号分别记为①、②,信号①可以表示为2fa+△(nfr1+f01+2fa),n为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移△fr1、载波包络相位飘移△f01和声光调制器驱动频率的飘移△fa,信号②可以表示为2fa+△(mfr1+f01+2fa),m为正整数,其中心频率在2fa,信号的飘移包括超短脉冲激光器11的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the Michelson interferometer 21 is respectively coupled to an optical filter 41 and another optical filter 43 through an optical beam splitter 31, and two photodetectors 51, 53 are used to measure the filtered optical signals respectively, and detect The signals are denoted as ① and ② respectively, the signal ① can be expressed as 2f a + △ (nf r1 + f 01 + 2f a ), n is a positive integer, its center frequency is 2fa, and the drift of the signal includes the ultrashort pulse laser 11 Repetition frequency drift △f r1 , carrier envelope phase drift △f 01 and AOM driving frequency drift △f a , signal ② can be expressed as 2f a + △(mf r1 +f 01 +2f a ), m is A positive integer whose center frequency is 2fa, and the signal drift includes the repetition frequency drift of the ultrashort pulse laser 11, the carrier envelope phase drift and the driving frequency drift of the acousto-optic modulator.
所述另一迈克尔逊干涉仪22的输出光经过光学分束器32分别耦合到光学滤波器42和另一光学滤波器44,采用两个光电探测器52、54分别测量滤波后的光学信号,探测到的信号分别记为③、④,信号③可以表示为2fb+△(nfr2+f02+2fb),n为正整数,其中心频率在2fb,信号的飘移包括超短脉冲激光器21的重复频率飘移△fr2、载波包络相位飘移△f02和声光调制器驱动频率的飘移△fb,信号④可以表示为2fb+△(mfr2+f02+2fb),m为正整数,其中心频率在2fb,信号的飘移包括另一超短脉冲激光器12的重复频率飘移、载波包络相位飘移和声光调制器驱动频率的飘移。The output light of the other Michelson interferometer 22 is respectively coupled to the optical filter 42 and another optical filter 44 through the optical beam splitter 32, and two photodetectors 52, 54 are used to measure the filtered optical signals respectively, The detected signals are denoted as ③ and ④ respectively, and the signal ③ can be expressed as 2f b + △ (nf r2 + f 02 + 2f b ), n is a positive integer, and its center frequency is at 2fb. The drift of the signal includes the ultrashort pulse laser 21 repetition frequency drift △f r2 , carrier envelope phase drift △f 02 and AOM drive frequency drift △f b , the signal ④ can be expressed as 2f b + △(mf r2 +f 02 +2f b ), m is a positive integer, its center frequency is 2fb, and the signal drift includes the repetition frequency drift of another ultrashort pulse laser 12, the carrier envelope phase drift and the drive frequency drift of the AOM.
将上述四个光电探测器51、52、53、54探测到的①、②、③和④信号按每个信号平均分为两路,一路不倍频,一路采用四个电学倍频及滤波单元61、2、63、64分别对上述①、②、③和④信号作倍频处理,对信号①和③作(p-1)倍的倍频处理,p为正整数,对信号②和④作p倍的倍频处理,得到四个倍频后的信号,分别为①×(p-1)、②×p、③×(p-1)、④×p;采用两个一级混频器71、72对倍频后的信号及不倍频的信号作第一级混频处理,得到四个第一级混频后的信号,分别为:①-②=(n-m)△fr1、②×p-①×(p-1)=2fa+△f01+(pgm-pgn+n)△fr1+2△fa、③-④=(n-m)△fr2、④×p-③×(p-1)=2fb+△f02+(pgm-pgn+n)△fr2+2△fb;采用二级混频器8对第一级混频后的信号作第二级混频处理,得到两个第二级混频后的信号,分别为:(①-②)-(③-④)=(n-m)△fr、[②×p-①×(p-1)]-[④×p-③×(p-1)]=2(fa-fb)+△f0+(pgm-pgn+n)△fr+2△(fa-fb);△f0+q△fr,上述两个信号实质上分别等于(n-m)△fr,△f0+q△fr,其中n、m、q为正整数,△fr为两台超短脉冲激光器11、12的相对重复频率抖动,△f0为两台超短脉冲激光器11、12的相对载波包络相位抖动,这两个信号就可以直接作为自适应双光梳光谱的补偿信号接入双光梳光谱数据处理单元9。Divide the ①, ②, ③ and ④ signals detected by the above four photodetectors 51, 52, 53, 54 into two paths on average for each signal, one path without frequency multiplication, and one path with four electrical frequency multiplication and filtering units 61, 2, 63, and 64 perform frequency multiplication processing on the above-mentioned ①, ②, ③ and ④ signals respectively, and perform (p-1) times frequency multiplication processing on the signals ① and ③, p is a positive integer, and perform frequency multiplication processing on the signals ② and ④ Perform p-fold frequency multiplication processing to obtain four multiplied signals, which are ①×(p-1), ②×p, ③×(p-1), ④×p; using two first-level mixing Devices 71 and 72 perform first-stage mixing processing on the multiplied signal and the non-multiplied signal to obtain four first-stage mixed signals, which are respectively: ①-②=(nm)△f r1 , ②×p-①×(p-1)=2f a +△f 01 +(pgm-pgn+n)△f r1 +2△f a 、③-④=(nm)△f r2 、④×p- ③×(p-1)=2f b +△f 02 +(pgm-pgn+n)△f r2 +2△f b ; use the second-stage mixer 8 to perform the second First-stage frequency mixing processing to obtain two second-stage mixed signals, which are: (①-②)-(③-④)=(nm)△f r , [②×p-①×(p-1 )]-[④×p-③×(p-1)]=2(f a -f b )+△f 0 +(pgm-pgn+n)△f r +2△(f a -f b ) ; △f 0 +q△f r , the above two signals are substantially equal to (nm)△f r , △f 0 +q△f r , where n, m, q are positive integers, and △f r is two The relative repetition frequency jitter of the ultrashort pulse lasers 11 and 12, △ f 0 is the relative carrier envelope phase jitter of the two ultrashort pulse lasers 11 and 12, these two signals can be directly used as compensation for the adaptive dual optical comb spectrum The signal is connected to the dual-comb spectral data processing unit 9.
本发明是这样对样品进行光谱测量的:将超短脉冲激光器11和另一超短脉冲激光器12两束光入射到样品池10中照射待测样品,然后两束入射光合为一束入射到干涉信号探测模块14,干涉信号探测模块14测量到两超短脉冲激光器11、12的干涉信号,将该干涉信号输入到自适应双光梳光谱系统数据处理单元9,自适应双光梳光谱系统数据处理单元9有三个输入信号,该三个输入信号分别表征为两台脉冲激光器11、12相对重复频率抖动的补偿信号、相对载波包络相位抖动的补偿信号和经过样品池10后的干涉信号,自适应双光梳光谱系统数据处理单元9通过采用表征两台脉冲激光器11、12相对重复频率抖动的补偿信号作为异步时钟采样信号,采用表征两台脉冲激光器11、12相对载波包络相位抖动的补偿信号与两台脉冲激光器11、12经过样品池10后的干涉信号混频的方法,分别消除两台脉冲激光器11、12相对重复频率和相对载波包络相位抖动对光谱测量的影响,获得高精度的光谱测量结果,最后通过光谱测量结果输出模块13将测量数据输出。The present invention measures the spectrum of the sample in this way: two beams of light from an ultrashort pulse laser 11 and another ultrashort pulse laser 12 are incident into the sample cell 10 to irradiate the sample to be measured, and then the two beams of incident light are combined into one beam and incident on the interferometer Signal detection module 14, the interference signal detection module 14 measures the interference signal of two ultrashort pulse lasers 11, 12, and the interference signal is input to the adaptive dual optical comb spectroscopy system data processing unit 9, and the adaptive dual optical comb spectroscopy system data The processing unit 9 has three input signals, and the three input signals are respectively characterized as the compensation signal of the two pulsed lasers 11 and 12 relative to the repetition frequency jitter, the compensation signal of the relative carrier envelope phase jitter, and the interference signal after passing through the sample cell 10, The data processing unit 9 of the adaptive dual-comb spectroscopy system adopts the compensation signal representing the relative repetition frequency jitter of the two pulsed lasers 11 and 12 as an asynchronous clock sampling signal, and adopts the signal representing the relative carrier envelope phase jitter of the two pulsed lasers 11 and 12 The compensation signal is mixed with the interference signal of the two pulsed lasers 11, 12 after passing through the sample cell 10, respectively eliminating the influence of the relative repetition frequency and relative carrier envelope phase jitter of the two pulsed lasers 11, 12 on the spectral measurement, and obtaining high The spectral measurement results with high precision are finally output by the spectral measurement result output module 13.
以上只是对本发明作进一步的说明,并非用以限制本专利的实施应用,凡为本发明等效实施,均应包含于本专利的权利要求范围之内。The above is only a further description of the present invention, and is not intended to limit the implementation and application of this patent. All equivalent implementations of the present invention should be included in the scope of claims of this patent.
Claims (4)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610308996.1A CN105823559B (en) | 2016-05-11 | 2016-05-11 | A kind of adaptive double light comb spectrally compensating method for extracting signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610308996.1A CN105823559B (en) | 2016-05-11 | 2016-05-11 | A kind of adaptive double light comb spectrally compensating method for extracting signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105823559A CN105823559A (en) | 2016-08-03 |
| CN105823559B true CN105823559B (en) | 2017-11-10 |
Family
ID=56529276
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610308996.1A Active CN105823559B (en) | 2016-05-11 | 2016-05-11 | A kind of adaptive double light comb spectrally compensating method for extracting signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105823559B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106197691B (en) * | 2016-07-05 | 2019-02-15 | 上海理工大学 | A detection device and detection method of ultrashort pulse carrier envelope phase |
| CN108613743B (en) * | 2018-06-14 | 2020-03-17 | 清华大学 | Noise compensation method of double-optical comb measurement system |
| CN109211414B (en) * | 2018-07-13 | 2020-10-16 | 华东师范大学 | Ultra-high precision optical frequency tester and its test method |
| CN111239075B (en) * | 2020-02-18 | 2022-10-11 | 华东师范大学重庆研究院 | A Measurement System of Gas Temperature and Multicomponent Concentration in Combustion Field Based on Adaptive Fiber Optic Comb |
| CN111786674B (en) * | 2020-07-09 | 2022-08-16 | 北京大学 | Analog bandwidth expansion method and system for analog-to-digital conversion system |
| CN111934806B (en) * | 2020-07-22 | 2022-06-21 | 上海交通大学 | Time-frequency synchronization device and synchronization method based on electro-optic modulation dual-optical comb |
| CN115507932B (en) * | 2022-08-23 | 2025-11-11 | 江苏科技大学 | Weak-reflection fiber bragg grating array weak vibration intensity and position demodulation system and method |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5647116B2 (en) * | 2008-07-25 | 2014-12-24 | サントル ナシオナル ドゥ ラ ルシェルシェサイアンティフィク(セエヌエールエス) | Synchronous interferometer with frequency comb |
| EP2939005A1 (en) * | 2012-12-28 | 2015-11-04 | CNR - Consiglio Nazionale Delle Ricerche | Frequency comb spectroscopy apparatus and method of frequency comb spectroscopy |
| US8541744B1 (en) * | 2013-03-09 | 2013-09-24 | Jian Liu | Method and apparatus for wavelength locking free optical frequency comb based differential absorption Lidar |
| CN104316180B (en) * | 2014-11-02 | 2016-06-01 | 华东师范大学 | Based on two optical frequency com optical imaging methods of continuous Frequency Stabilized Lasers |
| CN105548036B (en) * | 2015-12-08 | 2018-04-06 | 上海理工大学 | A kind of adaptive double light comb spectroscopic systems |
-
2016
- 2016-05-11 CN CN201610308996.1A patent/CN105823559B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105823559A (en) | 2016-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105823559B (en) | A kind of adaptive double light comb spectrally compensating method for extracting signal | |
| CN103712689B (en) | Continuous laser device spectral line width measurement device based on optical frequency comb | |
| US9164240B2 (en) | Optical buffering methods, apparatus, and systems for increasing the repetition rate of tunable light sources | |
| US11085831B2 (en) | High-resolution real-time time-frequency domain measuring device and method for ultra-short pulse | |
| CN106769952B (en) | Gas differential absorption lidar based on incoherent light source | |
| CN109357763B (en) | Atmospheric absorption spectrum measurement system and method based on time-resolved optical frequency comb | |
| CN106197691B (en) | A detection device and detection method of ultrashort pulse carrier envelope phase | |
| CN109346913B (en) | A dual femtosecond laser optical frequency comb locking device based on fiber delay line | |
| CN105548036B (en) | A kind of adaptive double light comb spectroscopic systems | |
| WO2020077735A1 (en) | Dual-frequency comb spectral focusing coherent anti-stokes raman spectroscopy detection system | |
| WO2013127370A1 (en) | Method and system for measuring optical asynchronous sampled signal | |
| US12444899B2 (en) | Method and apparatus for controlling a pulse repetition rate of a pulsed laser beam, and pulsed laser oscillator with stabilized pulse repetition rate | |
| CN109060150B (en) | Ultra-short pulse time width measuring device and method based on spectral interference | |
| CN106990089A (en) | The coherent anti-stokes raman scattering imaging system and imaging method of synchronous frequency reducing | |
| CN104914443B (en) | A kind of high-precision laser range-finding method of quick scanning | |
| CN205642621U (en) | Spectrum compensation signal extraction system is combed to two light of self -adaptation | |
| KR100996140B1 (en) | Cross correlator for measuring ultra-short laser pulses and method | |
| US20250172489A1 (en) | Time-resolved spectrum rapid measurement system and method | |
| US8995796B2 (en) | System for generating a beat signal | |
| US20240377323A1 (en) | Distributed time resolved fluorescence sensor using temporally correlated photons | |
| JP2012132704A (en) | Peak power monitoring device and peak power monitoring method | |
| JP4164599B2 (en) | Optical frequency measuring apparatus and measuring method using multicolor mode-locked laser | |
| CN201269780Y (en) | Dual optical path denoising carrier-envelope phase measurement device | |
| Cao et al. | Investigation on the pressure broadening of 85Rb 5S1/2− 5D3/2 monochromatic two− photon transition spectrum by multiple fluorescence detection | |
| Pan et al. | Modeling and Optimization of the OFC Spectral Filtering for High SNR Interference with CW Lasers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180830 Address after: 200237 District 2319, room 69, Lane 1985, Chunshen Road, Minhang District, Shanghai 1 district. Co-patentee after: East China Normal University Patentee after: Shanghai Langyan Optoelectronics Technology Co.,Ltd. Address before: 200237 District 2319, room 69, Lane 1985, Chunshen Road, Minhang District, Shanghai 1 district. Patentee before: Shanghai Langyan Optoelectronics Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211011 Address after: Room 4001, building 21, No. 123, Lane 1165, Jindu Road, Minhang District, Shanghai 201100 Patentee after: SHANGHAI LANGYAN OPTOELECTRONICS TECHNOLOGY Co.,Ltd. Patentee after: Chongqing Research Institute of East China Normal University Address before: 200237 zone 1, room 2319, No. 69, Lane 1985, Chunshen Road, Minhang District, Shanghai Patentee before: SHANGHAI LANGYAN OPTOELECTRONICS TECHNOLOGY Co.,Ltd. Patentee before: East China Normal University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230523 Address after: Building 1, No. 2 Huizhu Road, Kangmei Street, Liangjiang New District, Chongqing, 401120 (No. 1, 2nd floor) Patentee after: Chongqing Menghe Biotechnology Co.,Ltd. Address before: Room 4001, building 21, No. 123, Lane 1165, Jindu Road, Minhang District, Shanghai 201100 Patentee before: SHANGHAI LANGYAN OPTOELECTRONICS TECHNOLOGY Co.,Ltd. Patentee before: Chongqing Research Institute of East China Normal University |