Fuse powder burning and explosion delay time measuring system and method based on high-precision sensor

Abstract

The invention discloses a system and a method for measuring delay time of fuze explosive burning and explosion based on a high-precision sensor, and belongs to the technical field of industrial testing. The method specifically comprises the following steps: collecting light intensity change and pressure change signal data of fuze explosive after ignition in real time; collected light intensity change and pressure change signal data of the fuze powder after ignition are converted into time domain signals, the time domain signals are processed, the converted time domain signals are analyzed, an initial peak value of a burning and explosion signal is identified, delay time of the fuze powder is calculated according to change of signal intensity, and the delay time of the fuze powder is calculated according to the delay time. And the calculated delay time is optimized and improved by verifying the difference between the measured delay time and the actual burning and explosion time through experiments.

Description

Technical Field

[0001] The invention belongs to the technical field of industrial testing, and in particular to a system and method for measuring the delay time of a fuze explosion based on a high-precision sensor. Background Art

[0002] The fuze is the core device that controls the detonation time of the fuze. Through reasonable time delay, it can ensure that the fuze takes effect at the best time. Accurate delay control can not only effectively improve the strike efficiency of the fuze, but also improve safety. Traditional fuzes mostly use mechanical and electronic delay control, but their control accuracy and reliability are limited, especially in harsh environments such as high temperature and severe vibration, which are prone to failure.

[0003] The fuze explosion delay time measurement technology based on high-precision sensors is an emerging delay measurement and control method. By embedding high-precision sensors into the fuze, various data of the fuze explosion process can be monitored in real time, and the explosion trigger can be controlled according to the preset delay time requirements. However, in harsh environments such as high temperature and severe vibration, there will be certain errors, making it difficult to make accurate judgments based on the preset time.

[0004] For example, the Chinese patent with the authorization announcement number CN113515030B discloses a method for measuring the delay time of the direct combustion and explosion of a fuse, which belongs to the field of industrial testing technology. The purpose of this invention is to use laser irradiation to directly detonate the fuse to measure the delay time of the direct combustion and explosion of a fuse. The steps of this invention are: preparation; the time from the moment the laser beam is emitted and irradiated to the upper surface of the fuse to the combustion and explosion of the fuse is 50ms; the photoelectric sensor can capture the ejected flame instantly; the difference between the moment the fuse absorbs light energy until the moment of combustion and explosion and the flame signal at the lower end of the fuse combustion and explosion detection sensor capturing the sound signal in real time is the delay time of the fuse combustion and explosion; the measurement results are displayed on the touch screen. This invention will have the advantages of simple operation, accurate measurement, high degree of automation, and high reliability. It can be specifically applied to the ignition test and blasting of pyrotechnics.

[0005] The defect of the above patent is that in harsh environments such as foggy, windy and sandy weather, it is impossible to accurately capture the flame and light, resulting in inaccurate delay time. Summary of the invention

[0006] In view of the shortcomings of the prior art, the present invention proposes a system and method for measuring the delay time of the fuze explosion based on high-precision sensors. Through high-precision sensors, data processing and complex mathematical models, the accuracy of the delay time measurement of the fuze explosion is improved.

[0007] To achieve the above object, the present invention provides the following technical solutions:

[0008] The method for measuring the delay time of the fuze explosion based on a high-precision sensor includes the following specific steps:

[0009] Step S1: real-time acquisition of light intensity change and pressure change signal data of the fuse charge after ignition;

[0010] Step S2: converting the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and processing them;

[0011] Step S3: Analyze the converted time domain signal to identify the initial peak of the explosion signal;

[0012] Step S4: Calculating the delay time of the fuse according to the change of the signal strength;

[0013] Step S5: verifying the difference between the measured delay time and the actual explosion time through experiments, and optimizing and improving the calculated delay time.

[0014] Specifically, the specific steps of step S2 are:

[0015] Step S201: Set the frequency of real-time sampling of light intensity change and pressure change signals to f s , then the sampling period is T sample ,

[0016] Step S202: representing the collected light intensity change and pressure change signals as a discrete time series x[n], where n represents the time index of the sampling point;

[0017] Step S203: Use a high-pass filter to remove low-frequency noise. The specific formula of the transfer function of the high-pass filter is:

[0018]

[0019] Where H(f) represents the transfer function of the high-pass filter, f c represents the cutoff frequency, p represents the order of the filter, and f represents the frequency index;

[0020] Step S204: De-noising the light intensity change and pressure change signals from which low-frequency noise has been removed. The specific de-noising formula is:

[0021]

[0022] in, represents the denoised time domain signal of light intensity change and pressure change, and h[n] represents the impulse response of the high-pass filter.

[0023] Specifically, the specific steps of step S3 are:

[0024] Step S301: De-noising the light intensity change and pressure change time domain signals Perform Fourier transform to convert the denoised light intensity change and pressure change time domain signals into frequency domain signals. The specific formula is:

[0025]

[0026] Where X[f] represents the frequency domain signal of light intensity change and pressure change after denoising, N represents the signal length, j represents the imaginary unit, and e represents the exponential function;

[0027] Step S302: Calculate the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula is:

[0028]

[0029] Wherein, |X[f]| represents the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising, Re(X[f]) represents the real part of the frequency domain signal of the light intensity change and pressure change after denoising, and Im(X[f]) represents the imaginary part of the frequency domain signal of the light intensity change and pressure change after denoising.

[0030] Specifically, the specific steps of step S3 also include:

[0031] Step S303: Setting a dynamic threshold A thershold , the initial peak of the explosion signal is identified in the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula of the dynamic threshold is:

[0032] A threshold =α×max(|X[f]|),

[0033] Among them, α represents the weight coefficient, and max() represents the maximum value function

[0034] Step S304: traverse the amplitude spectrum |X[f]| of the frequency domain signal of the light intensity change and pressure change after denoising, and make a judgment on each frequency. The specific formula for judgment is:

[0035]

[0036] Among them, |X[f i ]| represents the amplitude of the ith frequency domain signal of the light intensity change and pressure change after denoising, |X[f i-1 ]| represents the amplitude of the i-1th frequency domain signal of the light intensity change and pressure change after denoising, |X[f i+1 ]| represents the amplitude of the i+1th frequency domain signal of the light intensity change and pressure change after denoising;

[0037] Step S305: When the judgment formula is satisfied, the peak frequency of the ith frequency domain signal of the denoised light intensity change and pressure change is set to P peak , and P peak =f i .

[0038] Specifically, the specific steps of step S4 are:

[0039] Step S401: Dynamic threshold A thershold Dynamic adjustment is performed, and the specific formula for adjustment is:

[0040] A′ threshold =β×S avg +ε,

[0041] Among them, A′ threshold represents the dynamic threshold after dynamic adjustment, β represents the sensitivity adjustment coefficient, S avg It represents the average value of the frequency domain signal of light intensity change and pressure change after denoising in the past certain period of time, and ε represents the adjustment coefficient;

[0042] Step S402: Analyze the changes of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal over time through linear regression, and calculate the change rate of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal per unit time. The specific formula is:

[0043]

[0044] Among them, k represents the rate of change of the intensity of the frequency domain signal of the light intensity change and pressure change after denoising per unit time, S init Indicates the strength of the initial signal, that is, the signal strength before ignition, T peak Indicates the timestamp when the peak value is detected, T init Indicates the timestamp of ignition;

[0045] Step S403: Calculate the delay time of the fuze explosion by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate. The specific formula is:

[0046]

[0047] Among them, T delay Indicates the delay time of the fuze explosion, T sample Indicates the sampling period.

[0048] Specifically, the specific steps of step S5 are:

[0049] Step S501: According to historical experimental data, the calculation result of the delay time of the fuze explosion is corrected, and the correction factor is set to C. The specific correction formula is:

[0050] T delay,correct =T delay ×(1+C),

[0051] Among them, T delay,correct Indicates the corrected delay time of the fuze explosion;

[0052] Step S502: Perform M experimental measurements on the measured delay time and take the average value. The final calculation formula for the delay time of the fuze explosion is:

[0053]

[0054] Among them, T delay,finalt It represents the final delay time of the fuze explosion, T delay,correct [q] represents the corrected fuze explosion delay time obtained from the qth experiment.

[0055] A delay time measurement system for the combustion and explosion of a fuze charge based on a high-precision sensor, which is used to implement the delay time measurement method for the combustion and explosion of a fuze charge based on a high-precision sensor, comprises: a data acquisition module, a signal conversion module, a signal analysis and recognition module, a delay time calculation module and a correction and optimization module;

[0056] The data acquisition module is used to collect the light intensity change and pressure change signal data of the fuse charge after ignition in real time;

[0057] The signal conversion module is used to convert the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and process them;

[0058] The signal analysis and identification module is used to analyze the converted time domain signal and identify the initial peak value of the explosion signal;

[0059] The delay time calculation module is used to calculate the delay time of the fuse according to the change of signal strength;

[0060] The correction and optimization module is used to verify the difference between the measured delay time and the actual explosion time through experiments, and to optimize and improve the calculated delay time. Specifically,

[0061] Specifically, the delay time calculation module includes: a dynamic threshold adjustment unit, a change rate calculation unit and a delay time calculation unit;

[0062] The dynamic threshold adjustment unit is used to dynamically adjust the dynamic threshold;

[0063] The change rate calculation unit is used to calculate the change rate of the denoised light intensity change and pressure change frequency domain signal intensity over time by linear regression analysis, and the change rate of the denoised light intensity change and pressure change frequency domain signal intensity per unit time;

[0064] The delay time calculation unit is used to calculate the delay time of the detonation of the fuse by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate.

[0065] An electronic device comprises a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of a method for measuring the delay time of a fuze explosion based on a high-precision sensor are implemented.

[0066] A computer-readable storage medium stores computer instructions, which, when executed, execute the steps of a method for measuring the delay time of a fuze explosion based on a high-precision sensor.

[0067] Compared with the prior art, the present invention has the following beneficial effects:

[0068] 1. The present invention proposes a method for measuring the delay time of the detonation of a fuse based on a high-precision sensor. The measurement method based on a high-precision sensor can achieve a measurement accuracy of nanoseconds or higher, and combined with signal conversion, ensures accurate control of the delay time.

[0069] 2. The present invention proposes a method for measuring the delay time of the explosion of a fuze based on a high-precision sensor. The method uses a high-precision sensor to monitor the explosion process of the medicine in real time, and can detect abnormal situations in time, such as incomplete combustion, too fast or too slow combustion, and respond quickly.

[0070] 3. The present invention proposes a method for measuring the delay time of the fuze explosion based on a high-precision sensor. Through multiple experiments, the delay time of the fuze explosion is measured, corrected and calculated, which greatly improves the measurement accuracy of the delay time of the fuze explosion. BRIEF DESCRIPTION OF THE DRAWINGS

[0071] Figure 1 A flow chart of a method for measuring the delay time of a fuze explosion based on a high-precision sensor provided by the present invention;

[0072] Figure 2 This is a diagram of the architecture of the delay time measurement system for the detonation of fuze based on a high-precision sensor provided by the present invention;

[0073] Figure 3 A schematic diagram of an electronic device provided by the present invention;

[0074] Figure 4 A schematic diagram of the computer-readable storage medium structure provided by the present invention. DETAILED DESCRIPTION

[0075] In order to make the technical means, creative features, objectives and effects achieved by the present invention easy to understand, in the description of the present invention, it should be noted that the terms "center", "up", "down", "left", "right", "vertical", "horizontal", "inside", "outside" and the like indicate the orientation or position relationship based on the orientation or position relationship shown in the accompanying drawings, which is only for the convenience of describing the present invention and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In addition, the terms "No. 1", "No. 2" and "No. 3" are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. The present invention is further explained below in conjunction with specific embodiments.

[0076] Example 1

[0077] See also Figure 1 The present invention provides an embodiment: a method for measuring the delay time of a fuze explosion based on a high-precision sensor, comprising the following specific steps:

[0078] Step S1: real-time acquisition of light intensity change and pressure change signal data of the fuse charge after ignition;

[0079] Highly sensitive photoelectric sensors and piezoelectric sensors are placed near the fuse to monitor light or pressure changes during the explosion process in real time.

[0080] Step S2: converting the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and processing them;

[0081] The specific steps of step S2 are:

[0082] Step S201: Set the frequency of real-time sampling of light intensity change and pressure change signals to f s , then the sampling period is T sample ,

[0083] Step S202: representing the collected light intensity change and pressure change signals as a discrete time series x[n], where n represents the time index of the sampling point;

[0084] The collected signal is represented as a discrete time series, i.e., a time domain signal, according to the sampling period and the time index of the sampling point;

[0085] Step S203: Use a high-pass filter to remove low-frequency noise. The specific formula of the transfer function of the high-pass filter is:

[0086]

[0087] Where H(f) represents the transfer function of the high-pass filter, f c represents the cutoff frequency, p represents the order of the filter, and f represents the frequency index. The frequency here is the frequency of all signals. During acquisition, there are signal fluctuations, which causes the signal to not be completely the set frequency f. s ;

[0088] Step S204: De-noising the light intensity change and pressure change signals from which low-frequency noise has been removed. The specific de-noising formula is:

[0089]

[0090] in, represents the denoised time domain signal of light intensity change and pressure change, and h[n] represents the impulse response of the high-pass filter.

[0091] Step S3: Analyze the converted time domain signal to identify the initial peak of the explosion signal;

[0092] The specific steps of step S3 are:

[0093] Step S301: De-noising the light intensity change and pressure change time domain signals Perform Fourier transform to convert the denoised light intensity change and pressure change time domain signals into frequency domain signals. The specific formula is:

[0094]

[0095] Where X[f] represents the frequency domain signal of light intensity change and pressure change after denoising, N represents the signal length, j represents the imaginary unit, and e represents the exponential function;

[0096] Step S302: Calculate the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula is:

[0097]

[0098] Wherein, |X[f]| represents the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising, Re(X[f]) represents the real part of the frequency domain signal of the light intensity change and pressure change after denoising, and Im(X[f]) represents the imaginary part of the frequency domain signal of the light intensity change and pressure change after denoising;

[0099] Step S303: Setting a dynamic threshold A thershold , the initial peak of the explosion signal is identified in the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula of the dynamic threshold is:

[0100] A threshold =α×max(|X[f]|),

[0101] Among them, α represents the weight coefficient, which is a constant less than 1 and is used to adjust the sensitivity of the threshold. Max() represents the maximum value function.

[0102] Step S304: traverse the amplitude spectrum |X[f]| of the frequency domain signal of the light intensity change and pressure change after denoising, and make a judgment on each frequency. The specific formula for judgment is:

[0103]

[0104] Among them, |X[f i ]| represents the amplitude of the ith frequency domain signal of the light intensity change and pressure change after denoising, |X[f i-1 ]| represents the amplitude of the i-1th frequency domain signal of the light intensity change and pressure change after denoising, |X[f i+1 ]| represents the amplitude of the i+1th frequency domain signal of the light intensity change and pressure change after denoising;

[0105] Step S305: When the judgment formula is satisfied, the peak frequency of the ith frequency domain signal of the denoised light intensity change and pressure change is set to P peak , and P peak =f i .

[0106] Step S4: Calculating the delay time of the fuse according to the change of the signal strength;

[0107] The specific steps of step S4 are:

[0108] Step S401: Dynamic threshold A thershold Dynamic adjustment is performed, and the specific formula for adjustment is:

[0109] A′ threshold =β×S avg +ε,

[0110] Among them, A′ threshold It represents the dynamic threshold after dynamic adjustment, β represents the sensitivity adjustment coefficient, which is a coefficient less than 1 and is used to adjust the sensitivity, S avg It represents the average value of the frequency domain signal of the light intensity change and pressure change after denoising in the past certain period of time. ε represents the adjustment coefficient, which is a small constant used to ensure that the threshold is not zero.

[0111] The dynamic threshold is set based on the average value of the current signal strength S avg , and automatically adjust according to environmental conditions to ensure that the combustion and explosion moments can be stably detected under different experimental conditions, overcoming the errors caused by fixed thresholds in unstable experimental environments. By adaptively adjusting the threshold, the recognition accuracy of the combustion and explosion peaks is improved;

[0112] Step S402: Analyze the changes of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal over time through linear regression, and calculate the change rate of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal per unit time. The specific formula is:

[0113]

[0114] Among them, k represents the rate of change of the intensity of the frequency domain signal of the light intensity change and pressure change after denoising per unit time, S init Indicates the strength of the initial signal, that is, the signal strength before ignition, T peak Indicates the timestamp when the peak value is detected, T init Indicates the timestamp when the ignition is started;

[0115] k can reflect the speed of signal strength increase. As the weight coefficient in the delay time calculation formula, when the signal strength increases rapidly, it means that the explosion is faster; when it increases slowly, the explosion delay is longer, making the delay time calculation more accurate.

[0116] Step S403: Calculate the delay time of the fuze explosion by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate. The specific formula is:

[0117]

[0118] Among them, T delay Indicates the delay time of the fuze explosion, T sample Indicates the sampling period.

[0119] In the formula, the difference between the signal strength peak and the dynamic threshold, combined with the signal strength change rate and the sampling period, can more accurately locate the delay time of the explosion;

[0120] Compared with the traditional direct measurement method of time points, this formula takes into account the signal rising rate by adding the signal change rate k, thereby achieving a more accurate estimation of the signal delay time.

[0121] Step S5: verifying the difference between the measured delay time and the actual explosion time through multiple experiments, and optimizing and improving the calculated delay time.

[0122] The specific steps of step S5 are:

[0123] Step S501: According to historical experimental data, the calculation result of the delay time of the fuze explosion is corrected, and the correction factor is set to C. The specific correction formula is:

[0124] T delay,correct =T delay ×(1+C),

[0125] Among them, T delay,correct Indicates the corrected delay time of the fuze explosion;

[0126] The correction factor C is based on historical experimental results and can more accurately reflect the system error. It is usually less than 1 to prevent over-correction of the delay time.

[0127] Step S502: Perform M experimental measurements on the measured delay time and take the average value. The final calculation formula for the delay time of the fuze explosion is:

[0128]

[0129] Among them, T delay,finalt It represents the final delay time of the fuze explosion, T delay,correct [q] represents the corrected fuze explosion delay time obtained from the qth experiment.

[0130] Example 2

[0131] See also Figure 2 , another embodiment provided by the present invention: a delay time measurement system for the detonation of a fuze based on a high-precision sensor, comprising: a data acquisition module, a signal conversion module, a signal analysis and recognition module, a delay time calculation module and a correction and optimization module;

[0132] The data acquisition module is used to collect the light intensity change and pressure change signal data of the fuse charge after ignition in real time;

[0133] The signal conversion module is used to convert the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and process them;

[0134] The signal analysis and identification module is used to analyze the converted time domain signal and identify the initial peak value of the explosion signal;

[0135] The delay time calculation module is used to calculate the delay time of the fuse according to the change of signal strength;

[0136] The correction and optimization module is used to verify the difference between the measured delay time and the actual explosion time through experiments, and to optimize and improve the calculated delay time.

[0137] The delay time calculation module includes: a dynamic threshold value adjustment unit, a change rate calculation unit and a delay time calculation unit;

[0138] The dynamic threshold adjustment unit is used to dynamically adjust the dynamic threshold;

[0139] The change rate calculation unit is used to calculate the change rate of the denoised light intensity change and pressure change frequency domain signal intensity over time by linear regression analysis, and the change rate of the denoised light intensity change and pressure change frequency domain signal intensity per unit time;

[0140] The delay time calculation unit is used to calculate the delay time of the detonation of the fuse by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate.

[0141] Example 3

[0142] Figure 3 1 is a schematic diagram of the structure of an electronic device provided by an embodiment of the present application. Figure 3 As shown, according to another aspect of the present application, an electronic device 500 is also provided. The electronic device 500 may include one or more processors and one or more memories. The memories store computer readable codes, and when the computer readable codes are executed by one or more processors, a method for measuring the delay time of the detonation of a fuze based on a high-precision sensor can be executed.

[0143] The method or system according to the embodiment of the present application can also be used by Figure 3 The electronic device architecture shown in FIG. Figure 3 As shown, the electronic device 500 may include a bus 501, one or more CPUs 502, a ROM 503, a RAM 504, a communication port 505 connected to a network, an input / output component 506, a hard disk 507, etc. The storage device in the electronic device 500, such as the ROM 503 or the hard disk 507, may store the delay time measurement method of the detonation of the fuze based on the high-precision sensor provided in the present application. Furthermore, the electronic device 500 may also include a user interface 508. Of course, Figure 3 The architecture shown is only exemplary and can be omitted according to actual needs when implementing different devices. Figure 3 One or more components of an electronic device are shown.

[0144] Example 4

[0145] Figure 4 Schematic diagram of a computer-readable storage medium structure provided by an embodiment of the present application. Figure 4As shown, a storage medium 600 according to an embodiment of the present application is shown. Computer-readable instructions are stored on the storage medium 600. When the computer-readable instructions are executed by the processor, the method for measuring the delay time of the fuze explosion based on a high-precision sensor according to the embodiment of the present application described with reference to the above figures can be executed. The storage medium 600 includes, but is not limited to, for example, volatile memory and / or non-volatile memory. Volatile memory may include, for example, random access memory (RAM) and cache memory (cache). Non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0146] It should be understood that the methods, apparatuses, and devices of the present application can be implemented in many ways. For example, the methods, apparatuses, and devices of the present application can be implemented by software, hardware, firmware, or any combination of software, hardware, and firmware. The above order of steps for the method is only for illustration, and the steps of the method of the present application are not limited to the order specifically described above unless otherwise specified. In addition, in some embodiments, the present application can also be implemented as a program recorded in a recording medium, which includes machine-readable instructions for implementing the method according to the present application. Therefore, the present application also covers a recording medium storing a program for executing the method according to the present application.

[0147] In addition, the parts of the above-mentioned technical solutions provided in the embodiments of the present application that are consistent with the implementation principles of the corresponding technical solutions in the prior art are not described in detail to avoid excessive redundancy.

[0148] The specific implementation modes as described above further describe the purpose, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation mode of the present invention and is not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for measuring the delay time of fuze explosion based on a high-precision sensor, characterized in that: The specific steps include: Step S1: real-time acquisition of light intensity change and pressure change signal data of the fuse charge after ignition; Step S2: converting the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and processing them; Step S3: Analyze the converted time domain signal to identify the initial peak of the explosion signal; Step S4: Calculating the delay time of the fuse according to the change of the signal strength; Step S5: verifying the difference between the measured delay time and the actual explosion time through experiments, and optimizing and improving the calculated delay time.

2. The method for measuring the delay time of fuze explosion based on a high-precision sensor as claimed in claim 1 is characterized in that: The specific steps of step S2 are: Step S201: Set the frequency of real-time sampling of light intensity change and pressure change signals to f s , then the sampling period is T sample , Step S202: representing the collected light intensity change and pressure change signals as a discrete time series x[n], where n represents the time index of the sampling point; Step S203: Use a high-pass filter to remove low-frequency noise. The specific formula of the transfer function of the high-pass filter is: Where H(f) represents the transfer function of the high-pass filter, f c represents the cutoff frequency, p represents the order of the filter, and f represents the frequency index; Step S204: De-noising the light intensity change and pressure change signals from which low-frequency noise has been removed. The specific de-noising formula is: in, represents the denoised time domain signal of light intensity change and pressure change, and h[n] represents the impulse response of the high-pass filter.

3. The method for measuring the delay time of fuze explosion based on a high-precision sensor as claimed in claim 2, characterized in that: The specific steps of step S3 are: Step S301: De-noising the light intensity change and pressure change time domain signals Perform Fourier transform to convert the denoised light intensity change and pressure change time domain signals into frequency domain signals. The specific formula is: Where X[f] represents the frequency domain signal of light intensity change and pressure change after denoising, N represents the signal length, j represents the imaginary unit, and e represents the exponential function; Step S302: Calculate the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula is: Wherein, |X[f]| represents the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising, Re(X[f]) represents the real part of the frequency domain signal of the light intensity change and pressure change after denoising, and Im(X[f]) represents the imaginary part of the frequency domain signal of the light intensity change and pressure change after denoising.

4. The method for measuring the delay time of fuze explosion based on a high-precision sensor as claimed in claim 3 is characterized in that: The specific steps of step S3 also include: Step S303: Setting a dynamic threshold A thershold , the initial peak of the explosion signal is identified in the amplitude spectrum of the frequency domain signal of the light intensity change and pressure change after denoising. The specific formula of the dynamic threshold is: A threshold =α×max(|X[f]|), Among them, α represents the weight coefficient, and max() represents the maximum value function Step S304: traverse the amplitude spectrum |X[f]| of the frequency domain signal of the light intensity change and pressure change after denoising, and make a judgment on each frequency. The specific formula for judgment is: Among them, |X[f i ]· represents the amplitude of the ith frequency domain signal of the light intensity change and pressure change after denoising, |X[f i-1 ]| represents the amplitude of the i-1th frequency domain signal of the light intensity change and pressure change after denoising, |X[f i+1 ]| represents the amplitude of the i+1th frequency domain signal of the light intensity change and pressure change after denoising; Step S305: When the judgment formula is satisfied, the peak frequency of the ith frequency domain signal of the denoised light intensity change and pressure change is set to P peak , and P peak =f i .

5. The method for measuring the delay time of fuze explosion based on a high-precision sensor as claimed in claim 4 is characterized in that: The specific steps of step S4 are: Step S401: Dynamic threshold A thershold Dynamic adjustment is performed, and the specific formula for adjustment is: A' threshold =β×S avg +e, Among them, A' threshold represents the dynamic threshold after dynamic adjustment, β represents the sensitivity adjustment coefficient, S avg It represents the average value of the frequency domain signal of light intensity change and pressure change after denoising in the past certain period of time, and ε represents the adjustment coefficient; Step S402: Analyze the changes of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal over time through linear regression, and calculate the change rate of the intensity of the denoised light intensity change and the intensity of the pressure change frequency domain signal per unit time. The specific formula is: Among them, k represents the rate of change of the intensity of the frequency domain signal of the light intensity change and pressure change after denoising per unit time, S init Indicates the strength of the initial signal, that is, the signal strength before ignition, T peak Indicates the timestamp when the peak value is detected, T init Indicates the timestamp when the ignition is started; Step S403: Calculate the delay time of the fuze explosion by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate. The specific formula is: Among them, T delay Indicates the delay time of the fuze explosion, T sample Indicates the sampling period.

6. The method for measuring the delay time of fuze explosion based on a high-precision sensor as claimed in claim 5, characterized in that: The specific steps of step S5 are: Step S501: According to historical experimental data, the calculation result of the delay time of the fuze explosion is corrected, and the correction factor is set to C. The specific correction formula is: T delay,correct =T delay ×(1+C), Among them, T delay,correct Indicates the corrected delay time of the fuze explosion; Step S502: Perform M experimental measurements on the measured delay time and take the average value. The final calculation formula for the delay time of the fuze explosion is: Among them, T delay,finalt It represents the final delay time of the fuze explosion, T delay,correct [q] represents the corrected fuze explosion delay time obtained from the qth experiment.

7. A system for measuring the delay time of a fuze explosion based on a high-precision sensor, used to implement the method for measuring the delay time of a fuze explosion based on a high-precision sensor as claimed in any one of claims 1 to 6, characterized in that: include: Data acquisition module, signal conversion module, signal analysis and identification module, delay time calculation module and correction and optimization module; The data acquisition module is used to collect the light intensity change and pressure change signal data of the fuse charge after ignition in real time; The signal conversion module is used to convert the collected signal data of light intensity change and pressure change of the fuse charge after ignition into time domain signals and process them; The signal analysis and identification module is used to analyze the converted time domain signal and identify the initial peak value of the explosion signal; The delay time calculation module is used to calculate the delay time of the fuse according to the change of signal strength; The correction and optimization module is used to verify the difference between the measured delay time and the actual explosion time through experiments, and to optimize and improve the calculated delay time.

8. The delay time measurement system for fuze explosion based on high-precision sensor as claimed in claim 7 is characterized in that: The delay time calculation module includes: a dynamic threshold value adjustment unit, a change rate calculation unit and a delay time calculation unit; The dynamic threshold adjustment unit is used to dynamically adjust the dynamic threshold; The change rate calculation unit is used to calculate the change rate of the denoised light intensity change and pressure change frequency domain signal intensity over time by linear regression analysis, and the change rate of the denoised light intensity change and pressure change frequency domain signal intensity per unit time; The delay time calculation unit is used to calculate the delay time of the detonation of the fuse by combining the peak signal strength, the set dynamic threshold, and the signal strength change rate.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for measuring the delay time of the detonation of a fuse based on a high-precision sensor as described in any one of claims 1-6 are implemented.

10. A computer-readable storage medium, characterized in that: Computer instructions are stored thereon, and when the computer instructions are executed by the processor, the steps of the method for measuring the delay time of the detonation of a fuse based on a high-precision sensor as described in any one of claims 1-6 are executed.

Images