阅读说明：本技术 一种应用于激光多普勒测速仪的运动起点自动定位方法 (Automatic movement starting point positioning method applied to laser Doppler velocimeter ) 是由 郝歌扬 吕沛 吴国俊 杨钰城 吕小鹏 于 2019-10-12 设计创作，主要内容包括：本发明涉及一种应用于激光多普勒测速仪的运动起点自动定位方法,其解决传统激光多普勒测速仪存在的两个弊端,1.起始运动阶段的速度测量精度较低；2.无法得到绝对时刻与被测物运动速度之间的关系。该方法包括以下步骤：1)数据获取；2)计算静止开始运动的时刻t<Sub>1</Sub>。该方法能够自动定位到被测物由静止状态开始运动的绝对时刻,便于得到绝对时刻与物体运动速度之间的关系,同时提高起始运动阶段的速度测量精度。(The invention relates to a movement starting point automatic positioning method applied to a laser Doppler velocimeter, which solves two defects of the traditional laser Doppler velocimeter, and 1. the speed measurement precision of the initial movement stage is lower; 2. the relationship between the absolute time and the velocity of the object to be measured cannot be obtained. The method comprises the following steps: 1) acquiring data; 2) calculating the time t of the stationary start movement 1 . The method can automatically position the absolute time when the measured object starts to move from a static state, is convenient to obtain the relation between the absolute time and the moving speed of the object, and simultaneously improves the speed measuring precision of the initial moving stage.)

1. A motion starting point automatic positioning method applied to a laser Doppler velocimeter is characterized by comprising the following steps:

1) data acquisition

1.1) setting the sampling frequency f of the laser Doppler velocimeter₀The sampling duration T, the laser Doppler velocimeter starts to collect data, and the motion state data of the measured object is obtained;

1.2) after the data acquisition of the laser Doppler velocimeter is finished, acquiring the time t when the data acquisition is started by an internal clock of the laser Doppler velocimeter₀；

1.3) wavelength lambda used for obtaining data collected by laser Doppler velocimeter₀And lower velocity limit v_min；

2) Calculating the time t of the stationary start movement₁

2.1) reading data collected by a laser Doppler velocimeter and parameter data obtained in the step 1);

2.2) calculating the total point number N of the acquired data and the minimum Doppler frequency difference omega according to the data obtained in the step 2.1)₀The calculation formula is as follows:

N＝f₀×T(4)

2.3) temporary counting bits i, j and temporary flag bits flag, flag required for initialization₁And flag₂Let i ═ j ═ 1, and flag ═ flag₁＝flag₂＝0；

2.4) putting the data read in the step 2.1) into an Array1, performing fast Fourier transform on the Array1 to obtain the frequency spectrum distribution of the data, finding the frequency point frequency with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency point frequency as f 1;

2.5) determining f1 and ω₀The size relationship of (1): if f 1-omega₀If the value is less than 0, the measured object is in a static state or the movement speed is less than the lower speed limit of the laser Doppler velocimeter, and the result that the measured object does not move is output after the calculation is finished; if f 1-omega₀If not less than 0, executing step 2.6);

2.6) dividing the Array1 into two groups, wherein the first group of data is put into an Array Temp1, and the second group of data is put into an Array Temp 2;

2.7) carrying out fast Fourier transform on the Temp1 logarithm group to obtain the frequency spectrum distribution of the Temp1 logarithm group, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 2;

2.8) determining f2 and ω₀The size relationship of (1): if f2 < omega₀Setting the Array Temp2 to the Array1, setting the counting position i to 2, and returning to step 2.6); if f2 is more than or equal to omega₀Then step 2.9) is executed;

2.9) calculating flag bit flag₁The calculation formula is as follows:

2.10) let Array2 ═ Array Temp 1;

2.11) equally dividing the Array2 into two parts, wherein the first part of data is placed in an Array Temp3, and the second part of data is placed in an Array Temp 4;

2.12) carrying out fast Fourier transform on the Temp4 logarithm group to obtain the frequency spectrum distribution of the Temp4 logarithm group, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 4;

2.13) judgment of f4- ω₀And ω₀The size relationship of (1): if f 4-omega₀＜ω₀Let the Array Temp3 be the Array2 and the count bit j be j +1, return to step 2.11); if f 4-omega₀≥ω₀Then step 2.14) is performed;

2.14) calculating flag bit flag₂The calculation formula is as follows:

2.15) calculating a flag bit flag, wherein the calculation formula is as follows:

flag＝flag₁+flag₂(8)

2.16) calculating the moment t of the measured object starting to move from the standstill according to the flag bit flag₁And output t₁The calculation formula of (2) is as follows:

2. a computer-readable storage medium having stored thereon a computer program, characterized in that: which computer program, when being executed by a processor, carries out the steps of the method as claimed in claim 1.

3. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein: the processor, when executing the program, implements the steps of the method of claim 1.

Technical Field

The invention relates to the field of laser Doppler velocity and acceleration measurement, in particular to an automatic positioning method for a motion starting point of a laser Doppler velocimeter.

Background

The laser Doppler velocity measurement technology is a high-precision non-contact velocity measurement technology based on the optical Doppler effect and is widely applied to various fields. Typical laserThe optical Doppler velocimeter has a structure shown in FIG. 1, wherein laser emitted from a laser passes through a front lens L₁The light beam is divided into two paths by a beam splitter, one path of reflected light irradiates an object, and the other path of transmitted light irradiates a reflector M₁A surface. Scattered light from the surface of the object to be measured and the beam M₁The two light beams reflected from the surface are reflected by the light beam M₂Reflected by the collimating mirror L₂And collimating and receiving by a detector, wherein the frequency of light scattered back by the measured object is changed due to the fact that the measured object is in a moving state, and the frequency of the light and the other light form a beat frequency on the surface of the detector.

If defined by M₁The reflected light is reference light with a frequency of the output frequency omega of the laser₁Defining the light scattered from the surface of the measured object as signal light, and changing the frequency of the signal light into omega due to Doppler effect₂＝ω₁+2u/λ, where u is the velocity of the object being measured and λ is the wavelength of the output light, the intensity of the two beams can be expressed as:

where E1 denotes the amplitude of the reference light, t denotes time,

indicating the phase of the reference light, E2 indicating the amplitude of the signal light,

indicating the phase of the signal light;

according to the interference formula of the light, the interference light intensity of the two beams of light can be expressed as:

due to the square law effect of the photodiode, only the frequency is ω₁-ω₂When the output frequency of the laser is fixed, the frequency of the interference fringe is only related to the speed of the motion of the object, so that the interference signal received by the photodiode is collected by the data acquisition card, the doppler frequency difference of the interference signal is calculated by the data processing system, and the motion speed of the object to be measured is:

when the laser Doppler velocimeter is adopted to measure the speed of a measured object which is converted from a static state to a moving state, the traditional laser Doppler velocimeter cannot automatically position the absolute moment when the moving state of the measured object changes, so that the following two disadvantages exist: 1. the speed measurement precision in the initial motion stage is low; 2. the relationship between the absolute time and the velocity of the object to be measured cannot be obtained.

Disclosure of Invention

In order to solve two defects of the traditional laser Doppler velocimeter, the invention provides an automatic movement starting point positioning method applied to the laser Doppler velocimeter.

In order to achieve the purpose, the invention is realized by the following technical scheme:

a motion starting point automatic positioning method applied to a laser Doppler velocimeter comprises the following steps:

1) data acquisition

1.1) setting the sampling frequency f of the laser Doppler velocimeter₀The sampling duration T, the laser Doppler velocimeter starts to collect data, and the motion state data of the measured object is obtained;

1.2) after the data acquisition of the laser Doppler velocimeter is finished, acquiring the time t when the data acquisition is started by an internal clock of the laser Doppler velocimeter₀；

1.3) wavelength lambda used for obtaining data collected by laser Doppler velocimeter₀And lower velocity limit v_min；

2) Calculating the time t of the stationary start movement₁

2.1) reading data collected by a laser Doppler velocimeter and parameter data obtained in the step 1);

2.2) calculating the total point number N of the acquired data and the minimum Doppler frequency difference omega according to the data obtained in the step 2.1)₀The calculation formula is as follows:

N＝f₀×T (4)

2.3) temporary counting bits i, j and temporary flag bits flag, flag required in the initialization algorithm₁And flag₂Let i ═ j ═ 1, and flag ═ flag₁＝flag₂＝0；

2.4) putting the data read in the step 2.1) into an Array1, performing fast Fourier transform on the Array1 to obtain the frequency spectrum distribution of the data, finding the frequency point frequency with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency point frequency as f 1;

2.5) determining f1 and ω₀The size relationship of (1): if f 1-omega₀If the value is less than 0, the measured object is in a static state or the movement speed is less than the lower speed limit of the laser Doppler velocimeter, and the result that the measured object does not move is output after the calculation is finished; if f 1-omega₀If not less than 0, executing step 2.6);

2.6) dividing the Array1 into two groups, wherein the first group of data is put into an Array Temp1, and the second group of data is put into an Array Temp 2;

2.7) carrying out fast Fourier transform on the Temp1 logarithm group to obtain the frequency spectrum distribution of the Temp1 logarithm group, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 2;

2.8) determining f2 and ω₀The size relationship of (1): if f2 < omega₀Setting the Array Temp2 to the Array1, setting the counting position i to 2, and returning to step 2.6); if f2 is more than or equal to omega₀Then step 2.9) is executed;

2.9) calculating flag bit flag₁The calculation formula is as follows:

2.10) let Array2 ═ Array Temp 1;

2.11) equally dividing the Array2 into two parts, wherein the first part of data is placed in an Array Temp3, and the second part of data is placed in an Array Temp 4;

2.12) carrying out fast Fourier transform on the Temp4 logarithm group to obtain the frequency spectrum distribution of the Temp4 logarithm group, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 4;

2.13) judgment of f4- ω₀And ω₀The size relationship of (1): if f 4-omega₀＜ω₀Let the Array Temp3 be the Array2 and the count bit j be j +1, return to step 2.11); if f 4-omega₀≥ω₀Then step 2.14) is performed;

2.14) calculating flag bit flag₂The calculation formula is as follows:

2.15) calculating a flag bit flag, wherein the calculation formula is as follows:

flag＝flag₁+flag₂(8)

2.16) calculating the moment t of the measured object starting to move from the standstill according to the flag bit flag₁And output t₁The calculation formula of (2) is as follows:

meanwhile, the present invention also provides a computer-readable storage medium having stored thereon a computer program which, when being executed by a processor, carries out the steps of the above-mentioned method.

Furthermore, the present invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the above method when executing the program.

Compared with the prior art, the invention has the following beneficial effects:

1. the invention provides a method for automatically positioning a movement starting point of a measured object when a laser Doppler velocimeter measures the movement speed of the measured object converted from a static state to a movement state. The method can automatically position the absolute time when the measured object starts to move from a static state, is convenient to obtain the relation between the absolute time and the moving speed of the object, and simultaneously improves the speed measuring precision of the initial moving stage.

2. The method can accurately position the starting point moment of the movement of the object to be measured in a section of laser Doppler velocity measurement data containing the static state and the moving state of the object, thereby being convenient for obtaining the corresponding relation between the absolute time and the movement speed of the object.

3. When the object to be measured starts to move from a static state, the movement speed is generally low, the Doppler frequency formed by the laser Doppler velocimeter is low, and the method can accurately position the object to the movement starting point, so that the identification precision of the Doppler frequency can be improved, and the accuracy of speed measurement is further improved.

Drawings

Fig. 1 is a structural diagram of a typical conventional laser doppler velocimeter;

fig. 2 is a flowchart of the method for automatically positioning the movement starting point of the laser doppler velocimeter of the present invention.

Detailed Description

The invention is described in further detail below with reference to the figures and specific embodiments.

The invention provides a method for automatically positioning the change moment of a motion state with high precision when a laser Doppler velocimeter is used for measuring the speed of a measured object which is changed from a static state to a motion state.

As shown in fig. 2, the invention discloses a method for automatically positioning a motion starting point when a measured object starts to move from a standstill by using a laser doppler velocimeter, which specifically comprises the following steps:

1) parameter acquisition

1.1) setting the sampling frequency f of the laser Doppler velocimeter₀SamplingThe duration T is that the laser Doppler velocimeter starts to collect data, and the motion state data of the measured object is obtained;

1.2) after the data acquisition of the laser Doppler velocimeter is finished, the internal clock of the laser Doppler velocimeter acquires the accurate time t when the data acquisition is started₀；

1.3) wavelength lambda used for obtaining data collected by laser Doppler velocimeter₀And lower velocity limit v_min(ii) a (wavelength is a fixed parameter, and the lower speed measurement limit is set by a user);

2) calculating the time t of the stationary start movement₁

2.1) reading data collected by a laser Doppler velocimeter and parameter data obtained in the step 1);

2.2) calculating the total point number N of the acquired data and the minimum Doppler frequency difference omega which can be resolved by the algorithm according to the data acquired in the step 2.1)₀The calculation formula is as follows:

N＝f₀×T (4)

2.3) initializing temporary counting bits i, j and temporary flag bits flag, flag₁And flag₂Let i ═ j ═ 1, and flag ═ flag₁＝flag₂＝0；

2.4) putting the acquired data (the data read in the step 2.1)) into an Array1, performing Fast Fourier Transform (FFT) on the Array1 to obtain the frequency spectrum distribution of the Array, and finding the frequency point frequency with the maximum frequency spectrum amplitude in the frequency spectrum distribution, wherein the frequency point frequency is recorded as f 1;

2.5) determining f1 and ω₀The size relationship of (1): if f 1-omega₀If the value is less than 0, the measured object is in a static state or the movement speed is less than the lower speed limit of the laser Doppler velocimeter, the algorithm is ended, and a conclusion that the measured object does not move is output; if f 1-omega₀If the value is more than or equal to 0, the loop is continuously executed downwards, and the step 2.6) is executed;

2.6) dividing the Array1 into two groups, wherein the first group of data is put into an Array Temp1, and the second group of data is put into an Array Temp 2;

2.7) carrying out FFT on the array Temp1 to obtain the frequency spectrum distribution of the array Temp1, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 2;

2.8) determining f2 and ω₀The size relationship of (1): if f2 < omega₀Setting the Array Temp2 to the Array1, setting the counting position i to 2, and returning to step 2.6); if f2 is more than or equal to omega₀The loop continues to be executed downwards, and step 2.9) is executed;

2.9) calculating flag bit flag₁The calculation formula is as follows:

2.10) let Array2 ═ Array Temp 1;

2.11) equally dividing the Array2 into two parts, wherein the first part of data is placed in an Array Temp3, and the second part of data is placed in an Array Temp 4;

2.12) carrying out FFT on the array Temp4 to obtain the frequency spectrum distribution of the array Temp4, finding the frequency of the frequency point with the maximum frequency spectrum amplitude in the frequency spectrum distribution, and recording the frequency as f 4;

2.13) judgment of f4- ω₀And ω₀The size relationship of (1): if f 4-omega₀＜ω₀Let the Array Temp3 be the Array2 and the count bit j be j +1, return to step 2.11); if f 4-omega₀≥ω₀Then step 2.14) is performed;

2.14) calculating flag bit flag₂The calculation formula is as follows:

2.15) calculating a flag bit flag, wherein the calculation formula is as follows:

flag＝flag₁+flag₂(8)

2.16) calculating and outputting the time t1 when the measured object starts to move from the static state according to the flag bit flag, wherein the calculation formula of t1 is as follows:

the embodiment of the invention also provides a computer readable storage medium for storing a program, and the program realizes the steps of the automatic positioning method of the motion starting point applied to the laser Doppler velocimeter when being executed. In some possible embodiments, the various aspects of the invention may also be implemented in the form of a program product comprising program code means for causing a terminal device to carry out the steps according to the various exemplary embodiments of the invention described in the method part above of the description, when said program product is run on the terminal device.

A program product for implementing the above method, which may employ a portable compact disc read only memory (CD-ROM) and include program code, may be run on a terminal device, a computer device, such as a personal computer. However, the program product of the present invention is not limited in this regard and, in the present document, a readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The program product may employ any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the readable storage medium include: an electrical connection having one or more wires, a portable disk, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.