CN115184945B - Signal restoration method for measuring pulse laser based on array detection method - Google Patents
Signal restoration method for measuring pulse laser based on array detection method Download PDFInfo
- Publication number
- CN115184945B CN115184945B CN202210945190.9A CN202210945190A CN115184945B CN 115184945 B CN115184945 B CN 115184945B CN 202210945190 A CN202210945190 A CN 202210945190A CN 115184945 B CN115184945 B CN 115184945B
- Authority
- CN
- China
- Prior art keywords
- pulse
- photodetector
- signal
- sampling
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 55
- 230000010354 integration Effects 0.000 claims abstract description 50
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 238000011084 recovery Methods 0.000 claims abstract description 16
- 230000004044 response Effects 0.000 claims description 14
- 239000003990 capacitor Substances 0.000 claims description 13
- 238000013507 mapping Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 6
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 2
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 1
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/4861—Circuits for detection, sampling, integration or read-out
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及脉冲激光信号复原,具体涉及一种基于阵列探测法测量脉冲激光的信号复原方法。The invention relates to pulse laser signal restoration, and in particular to a signal restoration method for measuring pulse laser based on an array detection method.
背景技术Background technique
在激光斜程大气传输试验中,基于光电探测器的阵列探测法是测量激光远场光斑强度时空分布的有效手段。现有技术中,光电探测器阵列由若干光电探测器按点阵分布排列组成,用于对光斑进行空间取样,该光电探测器阵列将光信号转换为电信号后进行信号采集处理,并计算得到各个光电探测器在阵列分布点处的功率值,最后按照功率的空间分布值进行图像复原,进而得到光斑图像及远场光斑参数。也就是说,只有获得每个脉冲的波形,才能获得其积分值来反映激光单脉冲能量,对采集电路的帧频要求较高,尤其是在测量窄脉宽、低重频脉冲激光应用场合很难满足要求。In the laser oblique atmospheric transmission test, the array detection method based on photodetectors is an effective means to measure the spatiotemporal distribution of the laser far-field spot intensity. In the prior art, the photodetector array is composed of a number of photodetectors arranged in a dot matrix distribution, which is used to perform spatial sampling of the spot. The photodetector array converts the optical signal into an electrical signal and then performs signal acquisition processing, and calculates the power value of each photodetector at the array distribution point. Finally, the image is restored according to the spatial distribution value of the power, and then the spot image and far-field spot parameters are obtained. In other words, only by obtaining the waveform of each pulse can its integral value be obtained to reflect the energy of a single laser pulse, and the frame rate requirements of the acquisition circuit are relatively high, especially in the application of measuring narrow pulse width and low repetition rate pulse laser, it is difficult to meet the requirements.
发明内容Summary of the invention
本发明的目的是解决现有的复原方法不能满足窄脉宽、低重频脉冲激光测量需求的技术问题,提供一种基于阵列探测法测量脉冲激光的信号复原方法。The purpose of the present invention is to solve the technical problem that the existing restoration method cannot meet the requirements of narrow pulse width and low repetition rate pulse laser measurement, and to provide a signal restoration method for measuring pulse laser based on array detection method.
为解决上述技术问题,本发明提供的技术解决方案如下:To solve the above technical problems, the technical solutions provided by the present invention are as follows:
一种基于阵列探测法测量脉冲激光的信号复原方法,其特殊之处在于:包括以下步骤:A signal recovery method for measuring pulsed laser based on array detection method, which is special in that it includes the following steps:
步骤1:将脉冲激光信号辐照在阵列探测器上,所述阵列探测器包括按照一定映射关系排布的N个光电探测器,用于将所述脉冲激光信号转换为N路脉冲电信号,其中N为正整数;Step 1: irradiate a pulsed laser signal onto an array detector, wherein the array detector includes N photoelectric detectors arranged according to a certain mapping relationship, and is used to convert the pulsed laser signal into N pulsed electrical signals, where N is a positive integer;
步骤2:所述阵列探测器使用N个电荷积分电路分别对N路所述脉冲电信号进行展宽,得到N路脉冲展宽电信号,并在展宽过程中,分别生成N个所述电荷积分电路对应的N条电荷放电曲线;Step 2: The array detector uses N charge integration circuits to respectively widen the N pulse electrical signals to obtain N pulse widened electrical signals, and in the widening process, generates N charge discharge curves corresponding to the N charge integration circuits respectively;
步骤3,对N路所述脉冲展宽电信号进行异步采样,得到N组采样数据;Step 3, asynchronously sampling the N pulse-stretched electrical signals to obtain N groups of sampled data;
步骤4,在所得的每组采样数据中分别选取M个点作为零点测量点,得到M*N个零点测量点,并根据M*N个零点测量点确定唯一的时间零点t0,其中M为正整数;Step 4, select M points from each set of sampled data as zero point measurement points to obtain M*N zero point measurement points, and determine a unique time zero point t 0 according to the M*N zero point measurement points, where M is a positive integer;
步骤5,在N条所述电荷放电曲线上分别选取一点作为采样点,得到N个采样点;Step 5, selecting one point on each of the N charge-discharge curves as a sampling point, to obtain N sampling points;
步骤6,通过N个所述采样点和唯一的时间零点t0,分别计算得到每个光电探测器对脉冲激光信号的响应电流i0;Step 6, calculating the response current i 0 of each photodetector to the pulse laser signal through the N sampling points and the unique time zero point t 0 ;
步骤7,根据每个所述光电探测器对脉冲激光信号的响应电流i0,分别计算得到脉冲激光信号照射在每个光电探测器上的光强I;Step 7, calculating the light intensity I of the pulsed laser signal irradiating each photodetector according to the response current i 0 of each photodetector to the pulsed laser signal;
步骤8,将所得N个光电探测器上的光强I,按照与N个光电探测器的映射关系一一排布,得到所述脉冲激光信号的功率密度分布,进而复原所述脉冲激光信号。Step 8, arranging the light intensities I on the N photodetectors one by one according to the mapping relationship with the N photodetectors, obtaining the power density distribution of the pulse laser signal, and then restoring the pulse laser signal.
进一步地,步骤2中的N个所述电荷积分电路串联。Furthermore, the N charge integration circuits in step 2 are connected in series.
进一步地,步骤2中,每个所述电荷积分电路均为一阶RC电路,包括负载R和积分电容C,定义所述电荷积分电路的时间常数为τ,则时间常数τ与负载R、积分电容C满足以下关系:Furthermore, in step 2, each of the charge integration circuits is a first-order RC circuit, including a load R and an integration capacitor C, and the time constant of the charge integration circuit is defined as τ, then the time constant τ and the load R and the integration capacitor C satisfy the following relationship:
τ=RCτ=RC
其中,每个所述电荷积分电路中负载R和积分电容C的值均根据以下步骤确定:The values of the load R and the integral capacitor C in each of the charge integration circuits are determined according to the following steps:
步骤2.1,定义Is为光电探测器的饱和阈值,Umax为光电探测器的最高输出电压,A为光电探测器光敏面的面积,Re表示光电探测器的响应率,则根据下述公式计算得到与之对应的电荷积分电路中负载R的值:Step 2.1, define Is as the saturation threshold of the photodetector, Umax as the maximum output voltage of the photodetector, A as the area of the photosensitive surface of the photodetector, and Re as the response rate of the photodetector. Then, the value of the load R in the corresponding charge integration circuit is calculated according to the following formula:
步骤2.2,将负载R的值代入τ=RC,计算得到所述电荷积分电路中积分电容C的值,进而确定与之对应的所述电荷积分电路。Step 2.2, substituting the value of the load R into τ=RC, calculating the value of the integral capacitor C in the charge integration circuit, and then determining the corresponding charge integration circuit.
进一步地,步骤4中确定唯一的时间零点t0,包括以下步骤:Furthermore, determining the unique time zero point t 0 in step 4 includes the following steps:
步骤4.1,将M*N个等效测量点按照时序进行排列;Step 4.1, arrange M*N equivalent measurement points in time sequence;
步骤4.2,依次比较M*N个等效测量点,将最先出现由零信号转变为脉冲信号的两个等效测量点作为零点判据;Step 4.2, compare M*N equivalent measurement points in sequence, and take the two equivalent measurement points where the zero signal turns into the pulse signal first as the zero point judgment criteria;
步骤4.3,选取零点判据中时序在前的等效测量点作为时间零点t0。Step 4.3, select the equivalent measurement point that is earlier in time sequence in the zero point criterion as the time zero point t 0 .
进一步地,步骤6中,所述光电探测器对脉冲激光信号响应电流i0的计算公式为:Furthermore, in step 6, the calculation formula of the photodetector's response current i0 to the pulsed laser signal is:
其中,t1表示采样点的采样时刻,U1表示t1时刻光电探测器的输出电压,(t1,U1)表示电荷放电曲线上所选取的采样点,w为脉冲激光信号经光电探测器转换后的脉冲电信号的脉宽。Wherein, t 1 represents the sampling time of the sampling point, U 1 represents the output voltage of the photodetector at time t 1 , (t 1 , U 1 ) represents the sampling point selected on the charge discharge curve, and w is the pulse width of the pulse electrical signal after the pulse laser signal is converted by the photodetector.
进一步地,步骤7中,脉冲激光信号照射在对应光电探测器上的光强I的计算公式为:Furthermore, in step 7, the calculation formula for the light intensity I of the pulsed laser signal irradiated on the corresponding photodetector is:
其中,α为光强衰减倍数。Among them, α is the light intensity attenuation multiple.
进一步地,步骤3中,通过采用模拟多路开关对多个光电探测器通道进行选通时序,以此实现多个光电探测器通道的异步采样,其中,每一帧采样内所有光电探测器切换选通一次。Furthermore, in step 3, multiple photodetector channels are gated in time by using an analog multiplexer, thereby achieving asynchronous sampling of multiple photodetector channels, wherein all photodetectors are switched on once in each frame of sampling.
进一步地,步骤5中,多个所述电荷放电曲线上的采样点均选取所述电荷放电曲线上与所述脉冲展宽电信号下降沿对应的曲线部分。Furthermore, in step 5, the sampling points on the plurality of charge discharge curves are all selected from the curve portion on the charge discharge curve corresponding to the falling edge of the pulse stretching electrical signal.
本发明相比现有技术具有的有益效果:Compared with the prior art, the present invention has the following beneficial effects:
1、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,通过阵列探测器将辐照在其上的脉冲激光信号转换为多路脉冲电信号,使用电荷积分电路对多路脉冲电信号分别进行展宽,得到多路脉冲展宽信号。也就是说,本发明在无需改变采集电路的情况下,通过电荷积分电路对脉冲信号的展宽,使其能够适用于测量窄脉宽、低重频信号。1. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method, which converts the pulsed laser signal irradiated thereon into multiple pulse electrical signals through an array detector, and uses a charge integration circuit to respectively widen the multiple pulse electrical signals to obtain multiple pulse widened signals. In other words, the present invention widens the pulse signal through a charge integration circuit without changing the acquisition circuit, so that it can be applied to measuring narrow pulse width and low repetition rate signals.
2、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,通过采用模拟多路开关对多个探测器通道进行选通时序,以实现多通道的异步采样,通过异步采样得到的多组采样数据确定时间零点,利用多帧频数据对比优化,最大限度降低时间零点的不确定度。2. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method, which realizes asynchronous sampling of multiple channels by adopting analog multi-way switch to select timing of multiple detector channels, determines time zero point by multiple groups of sampling data obtained by asynchronous sampling, and utilizes multi-frame rate data comparison optimization to minimize the uncertainty of time zero point.
3、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,本发明中信号采样只需在电荷放电曲线上选取一个采样点,相较于现有技术中的基于奈奎斯特采样定律的信号采样至少需要两个采样点来说,有效降低对采样频率的要求,使得采样频率的需求至少降低两倍。3. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method. In the present invention, signal sampling only requires selecting one sampling point on the charge discharge curve. Compared with the signal sampling based on Nyquist sampling theorem in the prior art which requires at least two sampling points, the requirement for sampling frequency is effectively reduced, so that the sampling frequency requirement is reduced by at least two times.
4、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,相较于现有技术只有采样频率高时才能获得靶面处激光信号功率随时间变化的曲线来说,本发明在采样频率低的情况下,只需通过计算得到每个探测器上所照射的脉冲激光信号的光强数据,即可获得该变化曲线,使用范围广。4. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method. Compared with the prior art, which can only obtain the curve of laser signal power changing with time at the target surface when the sampling frequency is high, the present invention can obtain the change curve by simply calculating the light intensity data of the pulsed laser signal irradiated on each detector when the sampling frequency is low, and has a wide range of applications.
5、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,相较于现有技术中窄脉宽和低重频脉冲信号需要更高的采样频率,本发明具有采样频率低的特点,因此针对窄脉宽和低重频脉冲信号可体现更大的优势。5. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method. Compared with the prior art that narrow pulse width and low repetition rate pulse signals require a higher sampling frequency, the present invention has the characteristic of low sampling frequency, so it can show greater advantages for narrow pulse width and low repetition rate pulse signals.
6、本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法,可同时适用于连续体制激光信号的测量,其脉冲周期可根据实际需求选择,具有广泛的适用性。6. The present invention provides a signal recovery method for measuring pulsed laser based on array detection method, which can also be applied to the measurement of continuous laser signals. Its pulse period can be selected according to actual needs and has wide applicability.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明提供的一种基于阵列探测法测量脉冲激光的信号复原方法实施例的流程框图;FIG1 is a flowchart of an embodiment of a signal recovery method for measuring pulsed laser based on an array detection method provided by the present invention;
图2为本发明实施例中电荷积分电路对应脉冲电信号展开生成的电荷放电曲线示意图;其中(a)为探测器对应脉冲激光信号转换后的脉冲电信号的波形示意图;(b)为电荷放电曲线的示意图;FIG2 is a schematic diagram of a charge discharge curve generated by the charge integration circuit corresponding to the pulse electrical signal in an embodiment of the present invention; wherein (a) is a schematic diagram of the waveform of the pulse electrical signal after the detector corresponds to the pulse laser signal conversion; (b) is a schematic diagram of the charge discharge curve;
图3为本发明实施例中电荷积分电路的电路示意图;FIG3 is a circuit diagram of a charge integration circuit according to an embodiment of the present invention;
图4为本发明实施例中N个探测器通道的时序切换图;FIG4 is a timing switching diagram of N detector channels in an embodiment of the present invention;
图5为本发明实施例中确定时间零点的原理图。FIG. 5 is a schematic diagram showing a principle of determining the time zero point in an embodiment of the present invention.
具体实施方式Detailed ways
下面结合附图和具体实施例对本发明作进一步地说明。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.
本发明提供一种基于阵列探测法测量脉冲激光的信号复原方法,参照图1,该方法包括以下步骤:The present invention provides a signal recovery method for measuring pulsed laser based on an array detection method. Referring to FIG. 1 , the method comprises the following steps:
步骤1:将脉冲激光信号辐照在阵列探测器上,阵列探测器包括按照一定映射关系排布的N个光电探测器,用于将脉冲激光信号转换为N路脉冲电信号,其中N为正整数;Step 1: irradiate the pulsed laser signal onto an array detector, the array detector comprising N photoelectric detectors arranged according to a certain mapping relationship, for converting the pulsed laser signal into N pulsed electrical signals, where N is a positive integer;
在本实施例中,当脉冲激光信号辐照在阵列光电探测器上,脉冲激光信号的功率密度分布在空间上呈现差异,而本实施例中的N个光电探测器通道可表示为[d1,d2,…,dn]。如图2中(a)所示,每路所述脉冲电信号的波形表达式均为:In this embodiment, when the pulse laser signal is irradiated on the array photodetector, the power density distribution of the pulse laser signal presents differences in space, and the N photodetector channels in this embodiment can be expressed as [d 1 , d 2 , …, d n ]. As shown in FIG. 2 (a), the waveform expression of each pulse electrical signal is:
其中,i0为光电探测器对脉冲激光信号的响应电流,u(t)为阶跃函数,w为脉冲电信号的脉宽,T为脉冲电信号的周期。Where i0 is the response current of the photodetector to the pulse laser signal, u(t) is the step function, w is the pulse width of the pulse electrical signal, and T is the period of the pulse electrical signal.
步骤2:阵列探测器使用N个电荷积分电路分别对N路所述脉冲电信号进行展宽,得到N路脉冲展宽电信号,并在展宽过程中,分别生成N个电荷积分电路对应的N条电荷放电曲线,如图2中(b)所示;Step 2: The array detector uses N charge integration circuits to respectively widen the N pulse electrical signals to obtain N pulse widened electrical signals, and in the widening process, generates N charge discharge curves corresponding to the N charge integration circuits, as shown in FIG2(b);
参照图3,在本实施例中,N个电荷积分电路串联,并且N个电荷积分电路均为一阶RC电路,包括负载R和积分电容C,定义电荷积分电路的时间常数为τ,则时间常数τ与负载R、积分电容C满足以下关系:3 , in this embodiment, N charge integration circuits are connected in series, and the N charge integration circuits are all first-order RC circuits, including a load R and an integration capacitor C. The time constant of the charge integration circuit is defined as τ, and the time constant τ, the load R, and the integration capacitor C satisfy the following relationship:
τ=RCτ=RC
为了能够实现多路脉冲电信号的展宽,需要先确定每个光电探测器相应的电荷积分电路,即确定每个电荷积分电路中负载R和积分电容C的值。在本实施例中,每个电荷积分电路中负载R和积分电容C的值均根据以下步骤确定:In order to achieve the broadening of multi-channel pulse electrical signals, it is necessary to first determine the charge integration circuit corresponding to each photodetector, that is, to determine the value of the load R and the integration capacitor C in each charge integration circuit. In this embodiment, the value of the load R and the integration capacitor C in each charge integration circuit is determined according to the following steps:
步骤2.1,定义Is为光电探测器的饱和阈值,Umax为光电探测器的最高输出电压,A为光电探测器光敏面的面积,Re表示光电探测器的响应率,则根据下述公式计算得到与该光电探测器相应的电荷积分电路中负载R的值:Step 2.1, define Is as the saturation threshold of the photodetector, Umax as the maximum output voltage of the photodetector, A as the area of the photosensitive surface of the photodetector, and Re as the response rate of the photodetector. Then, the value of the load R in the charge integration circuit corresponding to the photodetector is calculated according to the following formula:
步骤2.2,根据该电荷积分电路中负载R的值,以及τ=RC,计算得到该电荷积分电路中积分电容C的值,进而得到该光电探测器相应的电荷积分电路。Step 2.2, according to the value of the load R in the charge integration circuit and τ=RC, the value of the integration capacitor C in the charge integration circuit is calculated, and then the charge integration circuit corresponding to the photodetector is obtained.
并且,在完成多路脉冲电信号的展宽后,根据上述波形表达式和电路基础知识,计算得到每个电荷积分电路中负载R两端的电压:Furthermore, after completing the broadening of the multi-channel pulse electrical signals, the voltage across the load R in each charge integration circuit is calculated based on the above waveform expression and basic circuit knowledge:
在本实施例中,设定电荷积分电路的放电过程至其峰值的电压的1%时截止,则根据上述负载R两端的电压表达式计算得到下述关系式:In this embodiment, the discharge process of the charge integrator circuit is set to be terminated when the voltage reaches 1% of its peak value. Then, the following relationship is calculated based on the voltage expression across the load R:
其中,U(T)为自脉冲激光输出零时刻后,T时刻光电探测器的输出电压,U(w)为自脉冲激光输出零时刻后,W时刻光电探测器的输出电压。根据上述关系式计算得到时间常数τ和脉冲周期T(T>>w,计算时可忽略)的关系式:Among them, U(T) is the output voltage of the photodetector at time T after the pulse laser outputs zero time, and U(w) is the output voltage of the photodetector at time W after the pulse laser outputs zero time. According to the above relationship, the relationship between the time constant τ and the pulse period T (T>>w, which can be ignored in the calculation) is calculated:
τ≈T/4.6τ≈T/4.6
当然,在其它实施例中,也可将电荷积分电路的放电截止电压设置为其峰值电压的2%,或者根据需求设定为其它数值,则根据上述计算过程,仍可计算出时间常数τ和周期T(T>>w,计算时可忽略)的关系式。Of course, in other embodiments, the discharge cut-off voltage of the charge integration circuit can also be set to 2% of its peak voltage, or set to other values according to requirements. Then, according to the above calculation process, the relationship between the time constant τ and the period T (T>>w, which can be ignored during calculation) can still be calculated.
步骤3,对N路所述脉冲展宽电信号进行异步采样,得到N组采样数据;Step 3, asynchronously sampling the N pulse-stretched electrical signals to obtain N groups of sampled data;
参照图4,在本实施例中,通过采用模拟多路开关对多个光电探测器通道进行选通时序,以此实现多个光电探测器通道的异步采样,其中,每一帧采样内所有光电探测器切换选通一次,保证在实现异步采样的同时,能够完成对N路所述脉冲展宽信号的全部采样。4 , in this embodiment, an analog multi-way switch is used to perform a selection timing on multiple photodetector channels, thereby realizing asynchronous sampling of multiple photodetector channels, wherein all photodetectors are switched and selected once in each frame sampling, thereby ensuring that all sampling of the N-channel pulse stretching signals can be completed while realizing asynchronous sampling.
步骤4,在所得的每组采样数据中分别选取M个点作为零点测量点,得到M*N个零点测量点,并根据M*N个零点测量点确定唯一的时间零点t0,其中M为正整数;Step 4, select M points from each set of sampled data as zero point measurement points to obtain M*N zero point measurement points, and determine a unique time zero point t 0 according to the M*N zero point measurement points, where M is a positive integer;
在本实施例中,确定唯一的时间零点t0,包括以下步骤:In this embodiment, determining a unique time zero point t 0 includes the following steps:
步骤4.1,将M*N个等效测量点按照时序进行排列;Step 4.1, arrange M*N equivalent measurement points in time sequence;
步骤4.2,依次比较M*N个等效测量点,将最先出现由零信号(忽略噪声影响)转变为脉冲信号的两个等效测量点作为零点判据;Step 4.2, compare M*N equivalent measurement points in sequence, and take the two equivalent measurement points that first appear to transform from zero signal (ignoring the influence of noise) to pulse signal as zero point judgment criteria;
步骤4.3,选取零点判据中的时序在前的等效测量点作为时间零点t0。Step 4.3, select the equivalent measurement point that comes first in the time sequence in the zero point criterion as the time zero point t 0 .
在本实施例中,光电探测器通道数为N,如图4所示,定义Ts为光电探测器的采样周期,每个光电探测器通道上的采样点数均为M,则等效测量点的数量为M*N,该等效测量点用于时间零点的确定,则相邻光电探测器通道间的切换时间间隔为:In this embodiment, the number of photodetector channels is N, as shown in FIG4 , Ts is defined as the sampling period of the photodetector, and the number of sampling points on each photodetector channel is M, then the number of equivalent measurement points is M*N, and the equivalent measurement points are used to determine the time zero point, then the switching time interval between adjacent photodetector channels is:
如图5所示,根据上述计算结果,则由零点位置误差引起的峰值电压误差可表示为:As shown in FIG5 , according to the above calculation results, the peak voltage error caused by the zero position error can be expressed as:
在本实施例中,假设上述峰值电压误差不超过其峰值电压的1%,则可计算得到在不超出该误差的情况下,所选取的等效测量点数M*N的最小值,即当ΔU<1%时,计算得到M*N的值大于458,即最小值为458。基于该计算结果,根据实际误差需求确定M*N的实际值,此时通过计算能够得到相邻光电探测器间的切换时间间隔。当然,在其它实施例中,该峰值电压误差也可以设定为其峰值电压的3%,或者根据需求设定为其它数值。需要注意的是,尽管在本实施例中,设定电荷积分电路中的放电截止电压为其峰值电压的1%,其峰值电压误差不超过其峰值电压的1%,但上述设定仅为本发明的一种优选实施例,在其它实施例中,也可以设定为其它数值,而且两者的设定并无关联性。In this embodiment, assuming that the peak voltage error does not exceed 1% of its peak voltage, the minimum value of the selected equivalent measurement points M*N can be calculated without exceeding the error, that is, when ΔU<1%, the value of M*N calculated is greater than 458, that is, the minimum value is 458. Based on the calculation result, the actual value of M*N is determined according to the actual error requirement, and the switching time interval between adjacent photodetectors can be obtained by calculation. Of course, in other embodiments, the peak voltage error can also be set to 3% of its peak voltage, or set to other values according to requirements. It should be noted that although in this embodiment, the discharge cut-off voltage in the charge integration circuit is set to 1% of its peak voltage, and its peak voltage error does not exceed 1% of its peak voltage, the above setting is only a preferred embodiment of the present invention. In other embodiments, it can also be set to other values, and the two settings are not related.
步骤5,在N条所述电荷放电曲线上分别选取一点作为采样点,得到N个采样点;Step 5, selecting one point on each of the N charge-discharge curves as a sampling point, to obtain N sampling points;
在本实施例中,为了得到稳定的采样数据,该采样点通常选取所述电荷放电曲线上与所述脉冲展宽电信号下降沿对应的曲线部分。In this embodiment, in order to obtain stable sampling data, the sampling point is usually selected from the curve portion of the charge discharge curve corresponding to the falling edge of the pulse stretching electrical signal.
步骤6,通过N个所述采样点和唯一的时间零点t0,分别计算得到每个光电探测器对脉冲激光信号的响应电流i0;Step 6, calculating the response current i 0 of each photodetector to the pulse laser signal through the N sampling points and the unique time zero point t 0 ;
在本实施例中,光电探测器对脉冲激光信号响应电流i0的计算公式为:In this embodiment, the calculation formula of the photodetector's response current i0 to the pulse laser signal is:
其中,t1表示采样点的采样时刻,U1表示t1时刻光电探测器的输出电压,(t1,U1)表示电荷放电曲线上所选取的采样点,根据上述公式可依次计算得到每个光电探测器对脉冲激光信号的响应电流i0。Wherein, t 1 represents the sampling time of the sampling point, U 1 represents the output voltage of the photodetector at time t 1 , (t 1 , U 1 ) represents the sampling point selected on the charge discharge curve, and the response current i 0 of each photodetector to the pulse laser signal can be calculated in turn according to the above formula.
步骤7,根据每个所述光电探测器对脉冲激光信号的响应电流i0,分别计算得到脉冲激光信号照射在每个光电探测器上的光强I;Step 7, calculating the light intensity I of the pulsed laser signal irradiating each photodetector according to the response current i 0 of each photodetector to the pulsed laser signal;
在本实施例中,脉冲激光信号照射在对应光电探测器上的光强I的计算公式为:In this embodiment, the calculation formula of the light intensity I of the pulse laser signal irradiated on the corresponding photodetector is:
其中,α为光强衰减倍数,根据上述公式可依次计算得到脉冲激光信号照射在每个光电探测器上的光强I。Among them, α is the light intensity attenuation multiple. According to the above formula, the light intensity I of the pulsed laser signal irradiated on each photodetector can be calculated in turn.
步骤8,将所得N个光电探测器上的光强I,按照与N个光电探测器的映射关系一一排布,得到所述脉冲激光信号的功率密度分布,进而复原所述脉冲激光信号。Step 8, arranging the light intensities I on the N photodetectors one by one according to the mapping relationship with the N photodetectors, obtaining the power density distribution of the pulse laser signal, and then restoring the pulse laser signal.
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制,对于本领域的普通专业技术人员来说,可以对前述各实施例所记载的具体技术方案进行修改,或者对其中部分技术特征进行等同替换,而这些修改或者替换,并不使相应技术方案的本质脱离本发明所保护技术方案的范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit them. For ordinary professional and technical personnel in this field, the specific technical solutions recorded in the aforementioned embodiments can be modified, or some of the technical features therein can be replaced by equivalents, and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions protected by the present invention.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210945190.9A CN115184945B (en) | 2022-08-08 | 2022-08-08 | Signal restoration method for measuring pulse laser based on array detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210945190.9A CN115184945B (en) | 2022-08-08 | 2022-08-08 | Signal restoration method for measuring pulse laser based on array detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115184945A CN115184945A (en) | 2022-10-14 |
CN115184945B true CN115184945B (en) | 2024-06-21 |
Family
ID=83523473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210945190.9A Active CN115184945B (en) | 2022-08-08 | 2022-08-08 | Signal restoration method for measuring pulse laser based on array detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115184945B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108007906A (en) * | 2017-11-30 | 2018-05-08 | 百色学院 | A kind of long-afterglow material phosphorescence excitation spectrum measuring system and method |
CN109238133A (en) * | 2018-07-23 | 2019-01-18 | 中国科学院长春光学精密机械与物理研究所 | The collecting method and device of the facula position measurement device of Gao Zhongying laser |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5283426A (en) * | 1993-05-17 | 1994-02-01 | X-Ray Scanner Corporation | Methods and apparatus for obtaining nonlinear responses from photodetector arrays |
CN102112054B (en) * | 2008-08-04 | 2014-07-02 | 皇家飞利浦电子股份有限公司 | Data acquisition |
CN108955872B (en) * | 2018-04-26 | 2020-12-01 | 中国科学院光电技术研究所 | A pulsed laser beacon wavefront detection device and method for eliminating the influence of sky light background |
CN113219477A (en) * | 2020-01-21 | 2021-08-06 | 华为技术有限公司 | Laser detector and method for manufacturing the same |
CN114567772B (en) * | 2022-02-23 | 2024-06-21 | 西北核技术研究所 | Measuring method and measuring system for optical gating time characteristics of image intensifier |
-
2022
- 2022-08-08 CN CN202210945190.9A patent/CN115184945B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108007906A (en) * | 2017-11-30 | 2018-05-08 | 百色学院 | A kind of long-afterglow material phosphorescence excitation spectrum measuring system and method |
CN109238133A (en) * | 2018-07-23 | 2019-01-18 | 中国科学院长春光学精密机械与物理研究所 | The collecting method and device of the facula position measurement device of Gao Zhongying laser |
Also Published As
Publication number | Publication date |
---|---|
CN115184945A (en) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5815532B2 (en) | Device for processing signals supplied by radiation detectors | |
JP6096408B2 (en) | Gamma ray detection system and gamma ray detection method | |
Park et al. | A wide dynamic-range CMOS image sensor using self-reset technique | |
CN114503545B (en) | Dual ramp pixel readout | |
KR20160010347A (en) | Sampling circuitry and sampling method for a plurality of electrodes | |
KR20190044261A (en) | Single-slope comparator with low-noise, and analog to digital converting apparatus and cmos image sensor thereof | |
TWI856038B (en) | Method for measuring amplitude of an incident particle flux | |
CN115184945B (en) | Signal restoration method for measuring pulse laser based on array detection method | |
JP2015049148A (en) | Distance measuring device | |
CN109283569B (en) | Apparatus for measuring photon information and photon measuring device | |
CN105572721A (en) | Device, equipment and method for measuring sensor gain | |
CN105866557B (en) | A kind of realization GHZThe time of pulse percent of pass, energy bispectrum synchronized measurement system | |
Broesch et al. | A digital long pulse integrator | |
Simpson et al. | Discriminator design considerations for time-interval measurement circuits in collider detector systems | |
Menezes et al. | A preamplifier with digital output for semiconductor detectors | |
CN115389017B (en) | A pulse laser power measurement method based on integration circuit | |
CN112084731A (en) | FPGA digital circuit and method for improving peak power measurement trigger dynamic range | |
Peña et al. | A simulation model of front-end electronics for high-precision timing measurements with low-gain avalanche detectors | |
CN113252965B (en) | Detection circuit, device and method | |
CN218239068U (en) | Optical power monitoring circuit | |
Clemêncio et al. | An all-digital coincidence-selection and coincidence-trigger generation for a small animal RPC-PET camera | |
Stevenson et al. | Comparative analysis of circuit performance for photocurrent estimation | |
KR20220129711A (en) | Time digital converter and lidar including same | |
Zhong et al. | A 6.25 M-pixel 270-fps Floating Gate Image Sensor Using High-speed and Low-power Readout Technology | |
Adamczewski-Musch et al. | Tests of the CBM Rich readout and Daq prototype |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |