Brillouin frequency shift extraction method and device, equipment and storage medium

Abstract

The invention provides a Brillouin frequency shift extraction method and device, equipment and a storage medium, and relates to the technical field of optical fiber sensing, and the method comprises the steps: carrying out the algebraic conversion of a measured Brillouin scattering gain spectrum into a Brillouin center frequency curve; according to the Brillouin center frequency curve, the approximate value of the center frequency of the maximum value position of the Brillouin scattering gain spectrum is solved, then the center frequency of the Brillouin scattering gain spectrum is determined, that is, a frequency shift extraction result is given in a multi-center-frequency statistical mode, and the method has the advantages of being high in precision and stability. Meanwhile, a fitting tool does not need to be used in the calculation process, and the calculation complexity is extremely low; compared with the prior art, the method provided by the invention has the advantages that the method is not easily influenced by a multi-peak Brillouin form, and particularly, a high-precision result can be given for common bimodal form data, so that the technical problems of large error, high calculation complexity and great influence of an initial value on the fitting success rate and the stability are solved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of optical fiber sensing, in particular to a Brillouin frequency shift extraction method and device, equipment and a storage medium. BACKGROUND

[0002] In Brillouin scattering gain spectrum analysis, using a Lorentz function to fit the scattering peak can more accurately extract the frequency shift information and other parameters related to the scattering process. Due to the influence of phonon attenuation, the Brillouin scattering gain spectrum generated in the optical fiber is in the form of a Lorentz spectrum line. The Brillouin frequency shift can be extracted from the Brillouin scattering gain spectrum. The change in the frequency shift of the stimulated Brillouin scattering light or the self-Brillouin scattering light of the optical fiber is linearly related to the temperature and strain to realize sensing.

[0003] The accuracy of Brillouin frequency shift extraction directly affects the sensing measurement effect. At present, the commonly used single-peak fitting method based on the Lorentz function will cause a large error if there are multiple peak components in the data to be fitted, whether the least squares method or the LM method is used. The Lorentz fitting of multiple Brillouin components has the ability to fit multiple component Lorentz shape components. However, its defect is high computational complexity. Since the specific parameters of the fitting are unknown, the initial value cannot be accurately known during the fitting process using the LM method, and the fitting success rate and stability will be greatly affected by the initial value. SUMMARY

[0004] Therefore, the embodiments of the present application aim to provide a Brillouin frequency shift extraction method, device, equipment and storage medium. According to the Brillouin scattering gain spectrum obtained by measurement, the Brillouin center frequency curve is determined. The approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum is solved according to the Brillouin center frequency curve, and then the center frequency of the Brillouin scattering gain spectrum is determined, that is, the frequency shift extraction result is given in a statistical manner of multiple center frequencies, which has the characteristics of high precision and high stability. At the same time, no fitting tool is needed during the calculation process, and the computational complexity is extremely small. It is not easily affected by the multi-peak Brillouin shape, especially for the commonly used double-peak shape data, and can give a high-precision result, thereby solving the above technical problems of "large error, high computational complexity, and fitting success rate and stability greatly affected by the initial value".

[0005] In a first aspect, the embodiments of the present application provide a Brillouin frequency shift extraction method, which comprises:

[0006] determining a Brillouin center frequency curve according to the Brillouin scattering gain spectrum obtained by measurement;

[0007] finding an approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve;

[0008] determine the center frequency of the Brillouin scattering gain spectrum based on the degree of curve fluctuation of the approximate center frequency.

[0009] In the implementation process, the center frequency curve of Brillouin is directly given by algebraic conversion of the Brillouin scattering gain spectrum obtained by measurement, and the center frequency of the Brillouin scattering gain spectrum is determined according to the center frequency curve of Brillouin, which can give the frequency shift extraction result in a statistical manner of multiple center frequencies, and has the characteristics of high precision, high stability and extremely small calculation complexity.

[0010] Optionally, the Brillouin center frequency curve is determined according to the Brillouin scattering gain spectrum obtained by measurement, and the Brillouin center frequency curve is determined according to the Brillouin scattering gain spectrum obtained by measurement.

[0011] The Brillouin frequency scattering gain is collected according to a preset frequency scanning interval to obtain a Brillouin shape curve in a Lorentz function form;

[0012] The derivative of the Brillouin scattering gain curve is calculated to obtain a Brillouin scattering derivative curve;

[0013] The Brillouin center frequency curve is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

[0014] In the implementation process, the center frequency is obtained by simple mathematical algebraic operation on the Brillouin scattering gain spectrum obtained by measurement, without using a fitting tool, and the calculation complexity is extremely small.

[0015] Optionally, the Brillouin center frequency curve is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve, and the Brillouin center frequency curve is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

[0016] The relationship between any two points and the center frequency is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

[0017] Based on the relationship between the two points and the center frequency, the half-height, half-width and Lorentz amplitude in the Brillouin scattering gain curve or the Brillouin scattering derivative curve are eliminated to obtain the Brillouin center frequency curve.

[0018] In the implementation process, the Brillouin center frequency curve is determined based on the relationship curve of any two points and the center frequency, which is used to calculate the center frequency of the Brillouin data, without a fitting process, thereby saving the calculation resources.

[0019] Optionally, the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum is found according to the Brillouin center frequency curve, and the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum is found according to the Brillouin center frequency curve.

[0020] introducing a half-height half-width range super parameter to determine the abscissa of the maximum value position of the Brillouin scattering gain spectrum; wherein the half-height half-width range super parameter is: the integral of half of the number of value points greater than the half-height on the left and right sides of the maximum value position of the Brillouin scattering gain spectrum.

[0021] According to the Brillouin center frequency curve, a series of approximate center frequency values near the maximum value position of the Brillouin scattering gain spectrum are solved in sections to obtain the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum.

[0022] In the above implementation process, a series of approximate center frequency values near the maximum value position of the Brillouin scattering gain spectrum are determined according to the Brillouin center frequency curve, which is fast and efficient, and can also give a high-precision calculation result for a multi-peak Brillouin morphology.

[0023] Optionally, the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum is obtained by solving a series of approximate center frequency values of the maximum value position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve, including:

[0024] According to the Brillouin center frequency curve, the Brillouin center frequency in the abscissa range on the left of the maximum value position determined according to the half-height half-width range super parameter is solved to obtain a series of approximate center frequency values on the left;

[0025] According to the Brillouin center frequency curve, the Brillouin center frequency in the abscissa range on the right of the maximum value position determined according to the half-height half-width range super parameter is solved to obtain a series of approximate center frequency values on the right;

[0026] The approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum is determined according to the series of approximate center frequency values on the left and the series of approximate center frequency values on the right.

[0027] In the above implementation process, the final Brillouin center frequency is determined by using the series of approximate center frequency values on the left and the right of the maximum value, which can also give a high-precision calculation result for a multi-peak Brillouin morphology, and improves the practicability.

[0028] Optionally, the approximate center frequency includes: a series of approximate center frequency values on the left and a series of approximate center frequency values on the right; and the center frequency of the Brillouin scattering gain spectrum is determined based on the curve fluctuation degree of the approximate center frequency, including:

[0029] The standard deviation, variance or entropy of the series of approximate center frequency values on the left is calculated to obtain the left curve fluctuation degree;

[0030] The standard deviation, variance or entropy of the series of approximate center frequency values on the right is calculated to obtain the right curve fluctuation degree;

[0031] determining the magnitude of the left curve fluctuation degree and the right curve fluctuation degree;

[0032] if it is determined that the left curve fluctuation degree is less than the right curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the maximum value or the average value of the left center frequency approximate series values;

[0033] if it is determined that the left curve fluctuation degree is greater than the right curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the maximum value or the average value of the right center frequency approximate series values;

[0034] if it is determined that the left curve fluctuation degree is equal to the right curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the average value of the maximum value of the left center frequency approximate series values and the maximum value of the right center frequency approximate series values.

[0035] In the above implementation process, the final Brillouin center frequency is determined according to the fluctuation degree of the center frequency approximate series values on both sides of the maximum value, avoiding the influence of the secondary peak, truly reflecting the center frequency position of the main peak, and improving the stability and accuracy of the calculation.

[0036] Optionally, the approximate center frequency includes left center frequency approximate series values and right center frequency approximate series values; and the curve fluctuation degree based on the approximate center frequency is used to determine the center frequency of the Brillouin scattering gain spectrum, including:

[0037] calculating the standard deviation or variance or entropy of the left center frequency approximate series values to obtain the left curve fluctuation degree;

[0038] calculating the standard deviation or variance or entropy of the right center frequency approximate series values to obtain the right curve fluctuation degree;

[0039] determining the magnitude of the left curve fluctuation degree and the right curve fluctuation degree;

[0040] if it is determined that the left curve fluctuation degree is less than the right curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the median of the left center frequency approximate series values;

[0041] if it is determined that the left curve fluctuation degree is greater than the right curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the median of the right center frequency approximate series values;

[0042] If it is determined that the left curve fluctuation degree is equal to the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the mean value of the median of the left center frequency approximate series value and the median of the right center frequency approximate series value.

[0043] In the implementation process, the data noise is large, and the use of the median is more stable.

[0044] In a second aspect, the embodiments of the present application provide a Brillouin frequency shift processing device, the device comprises:

[0045] A center frequency curve determination module is configured to determine a Brillouin center frequency curve according to a Brillouin scattering gain spectrum obtained by measurement.

[0046] A maximum gain value finding module is configured to find an approximate center frequency of a maximum value position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve.

[0047] A Brillouin frequency shift extraction module is configured to determine a center frequency of the Brillouin scattering gain spectrum based on a curve fluctuation degree of the approximate center frequency.

[0048] In a third aspect, the embodiments of the present application also provide an electronic device, which comprises a processor and a memory, the memory stores machine readable instructions executable by the processor, when the electronic device is running, the machine readable instructions are executed by the processor to perform the steps of the above method.

[0049] In a fourth aspect, the embodiments of the present application provide a storage medium, the storage medium stores a computer program, when the computer program is run by a processor, the steps of the above method are performed.

[0050] In order to make the above purposes, features and advantages of the present application more obvious and easy to understand, the following embodiments are described in detail below, and the accompanying drawings are described as follows. BRIEF DESCRIPTION OF DRAWINGS

[0051] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiments of the present application, it should be understood that the following drawings only show some embodiments of the present application, therefore should not be regarded as a limitation to the scope, for those skilled in the art, without paying creative labor, other related drawings can also be obtained according to these drawings.

[0052] Figure 1 A flow chart of a Brillouin frequency shift extraction method provided by the embodiments of the present application;

[0053] Figure 2 A Brillouin center frequency curve schematic diagram provided by the embodiments of the present application;

[0054] Figure 3 A Brillouin morphology curve provided for the embodiment of the present application;

[0055] Figure 4 A functional module schematic diagram of the Brillouin frequency shift extraction device provided for the embodiment of the present application;

[0056] Figure 5 A block schematic diagram of the electronic device of the Brillouin frequency shift extraction device provided for the embodiment of the present application.

[0057] Icon: 210 - central frequency curve determination module; 220 - maximum gain finding module; 230 - Brillouin frequency shift extraction module; 300 - electronic device; 311 - memory; 312 - storage controller; 313 - processor; 314 - peripheral interface; 315 - input / output unit; 316 - display unit. DETAILED DESCRIPTION

[0058] The technical solutions in the embodiments of the present application will be clearly and completely described with reference to the accompanying drawings of the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. The components of the embodiments of the present application described and shown in the accompanying drawings can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the present application.

[0059] It should be noted that: similar reference numerals and letters represent similar items in the following drawings, therefore, once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings. The term "comprise", "include" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, article or device including a series of elements not only includes those elements, but also includes other elements not explicitly listed, or further includes elements inherent to such process, method, article or device. Without more limitations, the element defined by the statement "comprising a" does not exclude the presence of additional identical elements in the process, method, article or device including the element. The terms "first", "second" and the like are only used to distinguish description, and cannot be understood as indicating or implying relative importance.

[0060] Before introducing the embodiments of the present application, first, briefly introduce the technical concepts involved in the present application.

[0061] Brillouin scattering: A phenomenon of inelastic scattering that occurs in solids, liquids, or gases, where an incident light wave, sound wave, or other form of wave interacts with thermal excitations or density fluctuations in the medium, resulting in a shift in the frequency of the scattered wave. This frequency shift is related to the speed of sound in the medium and the scattering angle, so by analyzing the frequency shift information of the Brillouin scattering gain spectrum, we can obtain important information about the internal structure and dynamic characteristics of the medium.

[0062] Spectral analysis and frequency shift extraction: Preprocess the collected scattering spectrum, such as smoothing, background subtraction, etc., to reduce noise and interference. Identify the Brillouin scattering peak, which is usually achieved by finding a specific peak value in the spectrum. Calculate the frequency shift of the Brillouin scattering peak, which can be achieved by comparing the difference between the scattering peak and the incident light frequency. The frequency shift (Δf) can be calculated by the following formula: Δf = f_scattered - f_incident; where f_scattered is the frequency of the scattered light, and f_incident is the frequency of the incident light.

[0063] Lorentzian function: A function commonly used to describe resonance phenomena with natural width, such as the resonance peak in Brillouin scattering. In the analysis of Brillouin scattering gain spectrum, using Lorentzian function to fit the scattering peak can more accurately extract the frequency shift information and other parameters related to the scattering process. The following are the general steps for using Lorentzian function fitting to calculate the frequency shift information: (1) Data preparation: Obtain the experimental data of the Brillouin scattering gain spectrum, including the correspondence between the scattering intensity and the frequency (or wavelength), ensure that the data points are sufficient and evenly distributed, so as to carry out accurate fitting.(2) Lorentzian function definition: The general form of Lorentzian function is:

[0064]

[0065] where A is the peak intensity, x0 is the center position of the resonance peak (i.e. the frequency shift), Γ 2 is the parameter related to the peak width (usually related to damping or natural width).(3) Fitting process: Choose appropriate fitting software or programming tools (such as SciPy library in Python), import experimental data into the fitting program. Define Lorentzian function as the fitting model, and set the initial parameter values (such as initial estimates of A, x0, Γ 2 ). Use nonlinear least squares method or other optimization algorithms for fitting to find the best fitting parameters.(4) Extract frequency shift information: After fitting, extract the center position x0 of the resonance peak from the fitting results, which is the frequency shift. At the same time, other parameters such as peak intensity A and peak width parameter Γ 2These parameters can provide more information about the scattering process. (4) Result analysis and verification: check the quality of the fitting results, such as the goodness of fit (R-squared value) and the residual distribution. If the fitting result is not ideal, you can try to adjust the parameters of the fitting model or consider other possible fitting functions. Combine the fitting results with experimental conditions, medium properties, etc. for further analysis and interpretation.

[0066] The inventors of the present application have noticed that the existing calculation method for the frequency shift information of the Brillouin scattering gain spectrum, the common way to calculate the frequency shift is the Lorentz function fitting, and the fitting process usually uses the least square method, LM method (Levenberg-Marquardt method, least square estimation method of regression parameters in nonlinear regression) and the like. (1) The single-peak fitting method based on the Lorentz function is commonly used, whether using the least square method or the LM method, due to the existence of multi-peak components in the data to be fitted, direct fitting will cause larger error. (2) The Lorentz fitting of multi-Brillouin components has the ability to fit multi-component Lorentz shape components, but the cost is high computational complexity, since the number of components to be fitted is unknown, the initial value cannot be accurately known during the fitting process using the LM method, and the fitting success rate and stability will be greatly affected by the initial value. (3) The LM method involves iteration, and in order to achieve the target accuracy, the computational complexity is very high.

[0067] Therefore, the Brillouin frequency shift extraction method, device, equipment and storage medium provided by the embodiments of the present application are introduced as follows, which extracts the Brillouin frequency shift from the Brillouin scattering gain spectrum, and then is applied to the application scenarios such as sensing, which utilizes the linear relationship between the change amount of the frequency shift of the stimulated Brillouin scattering light or the self-phase Brillouin scattering light of the optical fiber and the temperature and strain.

[0068] Please refer to Figure 1 , Figure 1 A flowchart of a Brillouin frequency shift extraction method provided by the embodiments of the present application is provided. The embodiments of the present application are described in detail below. The method comprises steps 100, 120 and 140.

[0069] Step 100: determining a Brillouin center frequency curve according to the Brillouin scattering gain spectrum obtained by measurement;

[0070] Step 120: finding an approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve;

[0071] Step 140: determining the center frequency of the Brillouin scattering gain spectrum based on the curve fluctuation degree of the approximate center frequency.

[0072] Exemplarily, the Brillouin center frequency curve can be: a Lorentz function form of the Brillouin scattering gain spectrum obtained by measurement, a derivative form of the Lorentz, or other associated relationship forms, a relationship of the Brillouin center frequency is directly calculated by using known data in the Brillouin scattering gain spectrum, and then the Brillouin center frequency curve corresponding to the relationship can be obtained, and only algebraic expression conversion is involved in the calculation process, without using fitting tools, and the calculation complexity is extremely small. The approximate center frequency of the maximum value position can be: a series of approximate center frequency values of the position coordinates of the peak value in the Brillouin scattering gain spectrum curve and its vicinity, which can be regarded as a plurality of approximate center frequency points corresponding to the position of the peak value point determined by the horizontal coordinate sampling frequency and the vertical coordinate gain or intensity on the Brillouin center frequency curve.

[0073] The Brillouin scattering gain spectrum corresponding to a certain point on the optical fiber is composed of the gain sum generated by each section of the optical fiber with the spatial resolution corresponding to the effective pulse length. When the optical fiber senses the existence of a local hot spot or deformation smaller than the spatial resolution, the Brillouin scattering gain spectrum in the transition region between the background optical fiber and the abnormal section is no longer a single Lorentz type spectrum, but a scattering spectrum with double peaks. If single-peak Lorentz fitting is still performed on the abnormal Brillouin scattering gain spectrum, the fitting effect cannot meet the accuracy requirement, and the Brillouin frequency shift has a large deviation, which cannot accurately reflect the actual value of the temperature or strain change. Since the double-peak situation in the multi-peak Brillouin scattering gain spectrum, the applicant finds that the influence of the secondary peak will cause the center frequency value on the secondary peak side of the main peak data to fluctuate strongly, while the data on the side not affected by the secondary peak is very stable and truly reflects the center frequency position of the main peak. Among them, the multi-peak can be double-peak, 3 peaks, 4 peaks, and more, and here the common double-peak is taken as an example for description.

[0074] Therefore, based on the above finding, by establishing a relationship through the Lorentz relationship of the Brillouin form, the Brillouin center frequency relationship curve can be directly obtained by using algebraic operation. According to the curve fluctuation degree of the approximate center frequency on the Brillouin center frequency relationship curve, the fluctuation degree of the numerical value on both sides of the maximum value position can be compared, and the center frequency on the side with smaller fluctuation is selected as the final center frequency of the Brillouin scattering gain spectrum, that is, the frequency shift extraction result is given in a statistical manner of a plurality of center frequencies, which can solve the defect that the main peak is affected by the secondary peak in the above double-peak. Generally, the double-peak is located near the fusion point of different optical fibers or near the temperature mutation. For long-distance measurement, it is generally considered that the Brillouin scattering gain spectrum with uniform numerical value has only one peak, and the mutation interval has two peaks.

[0075] In one embodiment, step 100 can include step 101, step 102, and step 103.

[0076] Step 101: Collect the gain data of the Brillouin frequency according to a preset scanning interval, and determine the Brillouin scattering gain curve;

[0077] Step 102: Calculate the derivative of the Brillouin scattering gain curve to obtain the Brillouin scattering derivative curve;

[0078] Step 103: Determine the Brillouin center frequency curve according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

[0079] Exemplarily, assuming that n is the number of frequency sweeps, X = [x1, x2, …, x n ] is the frequency sweep (Brillouin frequency), Y = [y1, y2, …, y n ] is the gain corresponding to the Brillouin frequency, the gain of the Brillouin frequency is scanned or collected according to a certain frequency interval (preset scanning interval), and the original Brillouin morphology curve can be measured, which is generally in the form of Lorentz, that is, the Brillouin scattering gain curve (that is, the Brillouin scattering gain spectrum). If the smoothness of the Brillouin scattering gain curve corresponding to Y is poor, a noise reduction method such as moving average or filtering can be used. Y' = [y1', y2', …, y n '] is the derivative curve of the gain corresponding to the Brillouin frequency, which can be approximated by using the central difference method to approximate the derivative of Y with respect to frequency, and the difference operation is performed by taking the values of the two sides. The derivative of the two ends has little effect on the calculation, and the left and right difference values can be replaced, or 0 can be replaced, so that the approximate value of the derivative of the Brillouin scattering gain spectrum can be obtained. Since the Brillouin scattering gain spectrum naturally has the form of the Lorentz function, the expression of the Lorentz function can be:

[0080]

[0081] The above function relationship is the Lorentz function form of the Brillouin scattering gain spectrum, that is, the Brillouin scattering gain curve, h can be the amplitude represented by the Lorentz function, w is the half-height half-width, f is the center frequency, and x is each sampling frequency. Then, the Brillouin scattering derivative curve, that is, the derivative of the Lorentz function, can be:

[0082]

[0083] By observation and calculation, it can be obtained that the relationship between the Lorentz function and its derivative is:

[0084]

[0085] According to the relationship between the Lorentz function and its derivative, the function relationship of the Brillouin center frequency can be further converted by algebra to determine the Brillouin center frequency curve.

[0086] In one embodiment, step 103 can include step 1031 and step 1032.

[0087] Step 1031: determining the relationship between any two points and the center frequency according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve;

[0088] Step 1032: eliminating the half-height half-width and the Lorentz amplitude in the Brillouin scattering gain curve or the Brillouin scattering derivative curve based on the relationship between any two points and the center frequency, to obtain the Brillouin center frequency curve.

[0089] Exemplarily, according to steps 101-103, the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve can be obtained, and for any point (x i , y i ), the center frequency f is moved to the left side of the equation, and through algebraic conversion, the relationship between any point and the center frequency can be obtained:

[0090]

[0091] Wherein, the amplitude h and the half-height half-width w represented by the Lorentz function are unknown terms, and the subscript i is an index, that is, x i , y i , y i ' respectively correspond to the i-th element of X, Y, Y', and by dividing the relationship between any two points and the center frequency, the solution is avoided, and for any two points i, j, the following can be obtained:

[0092]

[0093] It can be seen that the half-height half-width w and the Lorentz amplitude h in the Brillouin scattering gain curve or the Brillouin scattering derivative curve are eliminated in the above formula, and only the data of x, the Brillouin scattering gain curve y, and the Brillouin scattering derivative curve y' are contained, which can be regarded as known conditions. Therefore, the relationship formula of the final Brillouin center frequency curve can be:

[0094]

[0095] Through the above formula, the data in the three lists of the sweep frequency X, the gain Y of the Brillouin frequency, and the center difference Y' of the gain of the Brillouin frequency can be used, and through algebraic conversion operation, the Brillouin center frequency can be directly obtained.

[0096] In one embodiment, step 120 can include step 121 and step 122.

[0097] Step 121: introduce a half-height half-width range hyperparameter to determine the abscissa of the maximum value position of the Brillouin scattering gain spectrum; wherein the half-height half-width range hyperparameter is: the integer of half of the number of value points greater than the half-height on both sides of the maximum value position of the Brillouin scattering gain spectrum.

[0098] Step 122: according to the Brillouin center frequency curve, piecewise solve a series of center frequency approximate values near the maximum value position abscissa of the Brillouin scattering gain spectrum to obtain the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum.

[0099] Exemplarily, find the Y maximum value h in the Brillouin scattering gain spectrum, and the Brillouin scattering gain spectrum data index value c corresponding to the maximum value h, the index value c corresponds to x c Corresponding. Define a hyperparameter v, which generally represents the number of half-height half-width points. When the system parameters are specified, this value can be considered as a constant, or it can be found by numerical methods. In the Brillouin scattering gain spectrum, the full width at half maximum, also known as the full width at half maximum, refers to the full width of the spectrum band when the height of the spectrum band is half the maximum height, that is, the scattering peak width when the peak height is half. It can also represent the energy resolution. Explanatorily, in the above Lorentz function expression of the Brillouin scattering gain spectrum, when f-w

[0100] For the hyperparameter v representing the number of half-height half-width points, if the total number of points greater than h / 2 on both sides of c is Q, then generally v=int(Q / 2). In order to prevent the index from exceeding the limit, adjust the width of the hyperparameter v: if c-v<1, then v=c-1; if c+v>n, then v=n-c, n is the n in Y=[y1,y2,...,y n ] corresponding to the index value c of the Y maximum value h in the Brillouin scattering gain spectrum, there are Q points greater than the half-height, and the number of points in the half-height half-width range is int(Q / 2). Explanatorily, the index value is the serial number corresponding to the frequency from small to large according to a certain frequency interval (preset scanning interval) scanning or collecting the Brillouin frequency of the Brillouin scattering gain spectrum, for example, the index value corresponding to x1 can be 1; the index value corresponding to x2 can be 2.

[0101] According to the Brillouin center frequency curve determined in steps 1031-1032, a series of center frequency approximate values near the maximum index c of the Brillouin scattering gain spectrum can be piecewise solved according to subsequent steps 1221-1223 to obtain the approximate center frequency of the maximum value position of the Brillouin scattering gain spectrum.

[0102] In one embodiment, step 122 may include steps 1221, 1222, and 1223.

[0103] Step 1221: Based on the Brillouin center frequency curve, solve for the Brillouin center frequency within the range of the horizontal axis to the left of the maximum value position determined by the hyperparameters of the half-height and half-width range, and obtain an approximate series of values ​​for the left center frequency.

[0104] Step 1222: Based on the Brillouin center frequency curve, solve for the Brillouin center frequency within the range of the horizontal axis to the right of the maximum value position determined by the hyperparameters of the half-height and half-width range, and obtain an approximate series of center frequency values ​​on the right side.

[0105] Step 1223: Determine the approximate center frequency of the Brillouin scattering gain spectrum based on the approximate series values ​​of the center frequency on the left and right sides.

[0106] For example, such as Figure 2 As shown, the horizontal axis represents the index corresponding to the sampling frequency, and the vertical axis represents the center frequency. The curve is precisely the Brillouin center frequency curve. The index is the index of the data setting X and Y list, for example, X = [x1, x2, ..., x...]. n If X is a given index i, then the index i of X represents X[i] = x. i If i = 1, then take the first element x1 in the list or range X.

[0107] Set the interval hyperparameter d, where d is an integer. For example, d = 1:

[0108] Calculate the center frequency f(i, id) corresponding to index i within the left range [cv, c]. Substitute this into the Brillouin center frequency curve formula for mathematical solution to obtain an approximate left-hand value of the list's center frequency: f left ={f(i,i-1)|i=[cv,c-v+1,…,c]}.

[0109] Calculate the center frequency corresponding to index i within the range [c, c+v] on the right: f(i, i+d). Substitute this into the Brillouin center frequency curve formula for mathematical solution to obtain an approximate value of the right-hand side of the list center frequency: f right ={f(i,i+1)|i=[c,c+1,…,c+v]}, then we have the expression f left =[f(cv,cvd),f(c-v+1,c-v+1-d),…,f(c,cd)]; f right =[f(c,c+d),f(c+1,c+1+d),…,f(c+v,c+v+d)].

[0110] Figure 2The center frequency approximation with v=3, d=1 is shown in the figure, which includes 4 left center frequency approximation series values and 4 right center frequency approximation series values:

[0111] f left =[f(c-3,c-4),f(c-2,c-3),…,f(c,c-1)];f right =[f(c,c+1),f(c+1,c+2),…,f(c+3,c+4)], the yellow triangle in the figure is f left , and the green inverted triangle is f left .

[0112] Based on the above approximation, the final center frequency of the Brillouin scattering gain spectrum is further obtained according to the fluctuation degree.

[0113] In one embodiment, the approximate center frequency of the maximum value position includes: left center frequency approximation series values and right center frequency approximation series values; step 140 can include: step 141, step 142, step 143, step 144, step 145, and can further include step 146.

[0114] Step 141: calculating the standard deviation or variance or entropy of the left center frequency approximation series values to obtain the left curve fluctuation degree;

[0115] Step 142: calculating the standard deviation or variance or entropy of the right center frequency approximation series values to obtain the right curve fluctuation degree;

[0116] Step 143: judging the size of the left curve fluctuation degree and the right curve fluctuation degree;

[0117] Step 144: if it is determined that the left curve fluctuation degree is less than the right curve fluctuation degree, then the center frequency of the Brillouin scattering gain spectrum is determined as the maximum value or the mean value of the left center frequency approximation series values;

[0118] Step 145: if it is determined that the left curve fluctuation degree is greater than the right curve fluctuation degree, then the center frequency of the Brillouin scattering gain spectrum is determined as the maximum value or the mean value of the right center frequency approximation series values;

[0119] In some embodiments, step 146 is: if it is determined that the left curve fluctuation degree is equal to the right curve fluctuation degree, then the center frequency of the Brillouin scattering gain spectrum is determined as the mean value of the maximum value of the left center frequency approximation series values and the maximum value of the right center frequency approximation series values.

[0120] For example, by comparing the fluctuation degrees of the left and right center frequency approximation values of the Brillouin center frequency curve, the maximum value of the center frequency approximation values on the side with smaller fluctuation can be selected as the final result, as shown inFigure 3 As shown, the abscissa is the index corresponding to the sampling frequency, and the ordinate is the gain or intensity. The curve data is the original measured Brillouin profile curve, and the value corresponding to the vertical axis result is the final determined Brillouin center frequency. The function std measuring the fluctuation degree can be the standard deviation, and the function can also be replaced by variance, or entropy, etc. The standard deviation or variance or entropy of the left center frequency approximate series value and the standard deviation or variance or entropy of the right center frequency approximate series value are calculated to obtain the left curve fluctuation degree and the right curve fluctuation degree: std(f left ), std(f right ).

[0121] If std(f left ) < std(f right ), the center frequency of the Brillouin scattering gain spectrum is max(f left ); otherwise, the center frequency of the Brillouin scattering gain spectrum is max(f right ).

[0122] The maximum value max and the mean value mean can both provide relatively high accuracy in the case of small data noise, and therefore, the max can also be replaced by the mean.

[0123] If std(f left ) = std(f right ), because the fluctuation degrees of the right center frequency approximate series value f right and the left center frequency approximate series value f left are equal, 1 / 2 × [max(f left ) + max(f right )] can be taken as the final center frequency at this time.

[0124] In particular, if v = 1, std(f left ) = std(f right ) = 0, because the right center frequency approximate series value f right and the left center frequency approximate series value f left each only have one element, 1 / 2 × (f left + f right ) can be taken as the final center frequency at this time.

[0125] In one embodiment, after step 143, the method can further include steps 1431, 1432 and 1433.

[0126] Step 1431: If it is determined that the left curve fluctuation degree is less than the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the median of the left center frequency approximate series value.

[0127] Step 1432: If it is determined that the fluctuation degree of the left curve is greater than the fluctuation degree of the right curve, the center frequency of the Brillouin scattering gain spectrum is determined as the median of the series of approximate values of the right center frequency;

[0128] Step 1433: If it is determined that the fluctuation degree of the left curve is equal to the fluctuation degree of the right curve, the center frequency of the Brillouin scattering gain spectrum is determined as the average of the median of the series of approximate values of the left center frequency and the median of the series of approximate values of the right center frequency.

[0129] Exemplarily, the median is a statistical term used to describe the value in the middle of a sorted set of data. When the number of data is odd, the median is the number in the middle of the sorted set; when the number of data is even, the median is the average of the two numbers in the middle of the sorted set. One of the main advantages of the median is that it is not sensitive to extreme values (i.e. very large or very small values), which makes it particularly useful in describing the central tendency of the data distribution, especially when the data distribution is skewed or contains outliers. In contrast, the mean (arithmetic mean) is more sensitive to extreme values. If the Brillouin scattering gain spectrum data has large noise, using the median is more stable, that is, after determining the fluctuation degree of the approximate center frequency on the Brillouin center frequency curve, the median of the series of approximate values of the center frequency on the side with smaller fluctuation (left or right) is determined as the final Brillouin center frequency.

[0130] If std(f left )=std(f right ), because the fluctuation degrees of the series of approximate values of the right center frequency f right and the series of approximate values of the left center frequency f left are equal, at this time, 1 / 2×[median(f left )+median(f right )] can be taken as the final center frequency.

[0131] In particular, if v=1, std(f left )=std(f right )=0, because the series of approximate values of the right center frequency f right and the series of approximate values of the left center frequency f left each have only one element, at this time, 1 / 2×(f left +f right ) can be taken as the final center frequency.

[0132] Please refer to Figure 4 , Figure 4 for a functional block diagram of a Brillouin frequency shift processing device provided by an embodiment of the present application. The device comprises:

[0133] a central frequency curve determination module 210, configured to determine a Brillouin central frequency curve according to the Brillouin scattering gain spectrum obtained by measurement;

[0134] a maximum gain position searching module 220, configured to search for an approximate central frequency of a maximum value position of the Brillouin scattering gain spectrum according to the Brillouin central frequency curve;

[0135] a Brillouin frequency shift extraction module 230, configured to determine the central frequency of the Brillouin scattering gain spectrum based on a curve fluctuation degree of the approximate central frequency.

[0136] Optionally, the Brillouin central frequency curve is determined according to the Brillouin scattering gain spectrum obtained by measurement, and the Brillouin central frequency curve includes:

[0137] The gain data of the Brillouin frequency is collected at a preset scanning interval, and a Brillouin scattering gain curve is determined;

[0138] A derivative of the Brillouin scattering gain curve is calculated to obtain a Brillouin scattering derivative curve;

[0139] The Brillouin central frequency curve is determined according to a relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

[0140] Optionally, the Brillouin central frequency curve is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve, and the Brillouin central frequency curve includes:

[0141] The relationship between any two points and the central frequency is determined according to the relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve;

[0142] Based on the relationship between the any two points and the central frequency, the half-height half-width and the Lorentz amplitude in the Brillouin scattering gain curve or the Brillouin scattering derivative curve are eliminated to obtain the Brillouin central frequency curve.

[0143] Optionally, the approximate central frequency of the maximum value position of the Brillouin scattering gain spectrum is searched for according to the Brillouin central frequency curve, and the Brillouin central frequency curve includes:

[0144] A half-height half-width range hyperparameter is introduced to determine the abscissa of the maximum value position of the Brillouin scattering gain spectrum, and the half-height half-width range hyperparameter is an integer of half of the number of value points greater than the half-height on the left and right sides of the maximum value position of the Brillouin scattering gain spectrum.

[0145] According to the Brillouin center frequency curve, a series of approximate center frequency values of the maximum position of the Brillouin scattering gain spectrum are solved in segments to obtain an approximate center frequency of the maximum position of the Brillouin scattering gain spectrum.

[0146] Optionally, the solving of the series of approximate center frequency values of the maximum position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve comprises:

[0147] According to the Brillouin center frequency curve, a series of left center frequency values are solved in a range of the Brillouin center frequency on the left side of the maximum position determined according to the half-height half-width range hyperparameter to obtain a series of left center frequency approximate values;

[0148] According to the Brillouin center frequency curve, a series of right center frequency values are solved in a range of the Brillouin center frequency on the right side of the maximum position determined according to the half-height half-width range hyperparameter to obtain a series of right center frequency approximate values;

[0149] According to the series of left center frequency approximate values and the series of right center frequency approximate values, an approximate center frequency of the maximum position of the Brillouin scattering gain spectrum is determined.

[0150] Optionally, the approximate center frequency comprises the series of left center frequency approximate values and the series of right center frequency approximate values, and the determination of the center frequency of the Brillouin scattering gain spectrum based on the degree of fluctuation of the curve on both sides of the approximate center frequency of the maximum position comprises:

[0151] A standard deviation, a variance or an entropy of the series of left center frequency approximate values is calculated to obtain a left curve fluctuation degree;

[0152] A standard deviation, a variance or an entropy of the series of right center frequency approximate values is calculated to obtain a right curve fluctuation degree;

[0153] It is judged whether the left curve fluctuation degree is greater than the right curve fluctuation degree;

[0154] If it is judged that the left curve fluctuation degree is less than the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as a maximum value or a mean value of the series of left center frequency approximate values;

[0155] If it is judged that the left curve fluctuation degree is greater than the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as a maximum value or a mean value of the series of right center frequency approximate values;

[0156] If it is determined that the left curve fluctuation degree is equal to the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the mean value of the maximum of the left center frequency approximate series value and the maximum of the right center frequency approximate series value.

[0157] Optionally, after the determination of the left curve fluctuation degree and the right curve fluctuation degree, the method further comprises:

[0158] If it is determined that the left curve fluctuation degree is less than the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the median of the left center frequency approximate series value.

[0159] If it is determined that the left curve fluctuation degree is greater than the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the median of the right center frequency approximate series value.

[0160] If it is determined that the left curve fluctuation degree is equal to the right curve fluctuation degree, the center frequency of the Brillouin scattering gain spectrum is determined as the mean value of the median of the left center frequency approximate series value and the median of the right center frequency approximate series value. Please refer to Figure 5 , Figure 5 is a block diagram of an electronic device. The electronic device 300 can include a memory 311, a storage controller 312, a processor 313, a peripheral interface 314, an input / output unit 315, and a display unit 316. Those skilled in the art can understand that Figure 5 The structure shown is only a schematic, which does not limit the structure of the electronic device 300. For example, the electronic device 300 can include more or fewer components than those shown in Figure 5 , or have a different configuration from that shown in Figure 5 .

[0161] The above-mentioned memory 311, storage controller 312, processor 313, peripheral interface 314, input / output unit 315, and display unit 316 are directly or indirectly electrically connected to each other to realize data transmission or interaction. For example, these elements can be electrically connected to each other through one or more communication buses or signal lines. The above-mentioned processor 313 is used to execute the executable modules stored in the memory.

[0162] The memory 311 can be, but is not limited to, a random access memory (RAM), a read only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), and the like. The memory 311 is configured to store a program. After receiving an execution instruction, the processor 313 executes the program. The method performed by the electronic device 300 defined in the process disclosed in any of the embodiments of the present application can be applied to the processor 313 or implemented by the processor 313.

[0163] The processor 313 can be an integrated circuit chip having a signal processing capability. The processor 313 can be a general purpose processor, including a central processing unit (CPU), a network processor (NP), and the like. The processor 313 can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The processor 313 can implement or execute the methods, steps, and logical block diagrams disclosed in the embodiments of the present application. The general purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0164] The peripheral interface 314 is configured to couple various input / output devices to the processor 313 and the memory 311. In some embodiments, the peripheral interface 314, the processor 313, and the memory controller 312 can be implemented in a single chip. In other embodiments, they can be implemented by independent chips respectively.

[0165] The input / output unit 315 is configured to provide input data for a user. The input / output unit 315 can be, but is not limited to, a mouse, a keyboard, and the like.

[0166] The display unit 316 described above provides an interaction interface (for example, a user operation interface) between the electronic device 300 and the user for the user to refer. In this embodiment, the display unit 316 can be a liquid crystal display or a touch display. The liquid crystal display or the touch display can display the process of executing the program by the processor.

[0167] The electronic device 300 in this embodiment can be used to execute each step in each method provided in the embodiments of the present application.

[0168] In addition, the embodiments of the present application further provide a storage medium, and the storage medium stores a computer program. When the computer program is run by a processor, the steps in the above method embodiments are executed.

[0169] The computer program product of the above method provided in the embodiments of the present application includes a storage medium storing program codes. The program codes include instructions for executing the steps in the above method embodiments. For details, refer to the above method embodiments, which will not be described here.

[0170] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are only schematic. For example, the division of the modules is only a logical function division. In actual implementation, another division mode can be used. For example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some communication interfaces. The coupling or communication connection can be electrical, mechanical or in other forms. The various functional modules in the embodiments of the present application can be integrated together to form a single independent part, or each module can exist independently, or two or more modules can be integrated to form a single independent part.

[0171] It should be noted that, if the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application can be embodied in the form of a software product in essence or in the form of a part or the technical solutions of the present application that make contributions to the prior art. The computer software product is stored in a storage medium, includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (Read-Only Memory, ROM), a random access memory (Random Access Memory, RAM), a magnetic disk or an optical disk, and various media that can store program codes.

[0172] In this document, relational terms such as first and second and the like can be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.

[0173] The above only describes the embodiments of the present application and is not used to limit the protection scope of the present application. For those skilled in the art, the present application can have various modifications and changes. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A Brillouin frequency shift extraction method, characterized by, The method comprises: determining a Brillouin central frequency curve according to a Brillouin scattering gain spectrum obtained by measurement; finding an approximate central frequency of a maximum value position of the Brillouin scattering gain spectrum according to the Brillouin central frequency curve; determining the central frequency of the Brillouin scattering gain spectrum based on a curve fluctuation degree of the approximate central frequency.

2. The method of claim 1, wherein, The method comprises: measuring a Brillouin scattering gain at a preset frequency scanning interval to obtain a Brillouin scattering gain curve; calculating a derivative of the Brillouin scattering gain curve to obtain a Brillouin scattering derivative curve; determining a Brillouin central frequency curve according to a relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve.

3. The method of claim 2, wherein, The method comprises: determining a relationship between any two points and a central frequency according to a relationship between the Brillouin scattering gain curve and the Brillouin scattering derivative curve; eliminating a half-height half-width and a Lorentz amplitude in the Brillouin scattering gain curve or the Brillouin scattering derivative curve based on the relationship between the any two points and the central frequency to obtain the Brillouin central frequency curve.

4. The method of claim 1, wherein, The method comprises: introducing a half-height half-width range hyperparameter to determine an abscissa of the maximum value position of the Brillouin scattering gain spectrum; wherein the half-height half-width range hyperparameter is an integer of a half of a number of value points greater than a half-height on left and right sides of the maximum value position of the Brillouin scattering gain spectrum. segmentally solving approximate series values of the central frequency near the abscissa of the maximum value position of the Brillouin scattering gain spectrum according to the Brillouin central frequency curve to obtain the approximate central frequency of the maximum value position of the Brillouin scattering gain spectrum.

5. The method of claim 4, wherein, The method comprises: solving the Brillouin central frequency in a left abscissa range of the maximum value position according to the half-height half-width range hyperparameter to obtain left central frequency approximate series values according to the Brillouin central frequency curve; solving the Brillouin central frequency in a right abscissa range of the maximum value position according to the half-height half-width range hyperparameter to obtain right central frequency approximate series values according to the Brillouin central frequency curve; determining the approximate central frequency of the maximum value position of the Brillouin scattering gain spectrum according to the left central frequency approximate series values and the right central frequency approximate series values.

6. The method according to claim 1 or 5, characterized in that, The approximate central frequency comprises the left central frequency approximate series values and the right central frequency approximate series values; and the method comprises: determining the central frequency of the Brillouin scattering gain spectrum based on a curve fluctuation degree of the approximate central frequency. calculating the standard deviation or variance or entropy of the left side center frequency approximate series values to obtain the left side curve fluctuation degree; calculating the standard deviation or variance or entropy of the right side center frequency approximate series values to obtain the right side curve fluctuation degree; judging the size of the left side curve fluctuation degree and the right side curve fluctuation degree; if it is judged that the left side curve fluctuation degree is less than the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the maximum value or the mean value of the left side center frequency approximate series values; if it is judged that the left side curve fluctuation degree is greater than the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the maximum value or the mean value of the right side center frequency approximate series values; if it is judged that the left side curve fluctuation degree is equal to the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the mean value of the maximum value of the left side center frequency approximate series values and the maximum value of the right side center frequency approximate series values.

7. The method according to claim 1 or 5, characterized in that, wherein, the approximate center frequency comprises left side center frequency approximate series values and right side center frequency approximate series values; and the center frequency of the Brillouin scattering gain spectrum is determined based on the curve fluctuation degree of the approximate center frequency, comprising: calculating the standard deviation or variance or entropy of the left side center frequency approximate series values to obtain the left side curve fluctuation degree; calculating the standard deviation or variance or entropy of the right side center frequency approximate series values to obtain the right side curve fluctuation degree; judging the size of the left side curve fluctuation degree and the right side curve fluctuation degree; if it is judged that the left side curve fluctuation degree is less than the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the median of the left side center frequency approximate series values; if it is judged that the left side curve fluctuation degree is greater than the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the median of the right side center frequency approximate series values; if it is judged that the left side curve fluctuation degree is equal to the right side curve fluctuation degree, determining the center frequency of the Brillouin scattering gain spectrum as the mean value of the median of the left side center frequency approximate series values and the median of the right side center frequency approximate series values.

8. A Brillouin shift processing apparatus characterized by comprising: The device comprises: a center frequency curve determination module configured to determine a Brillouin center frequency curve according to a Brillouin scattering gain spectrum obtained by measurement; a maximum gain value searching module configured to search for an approximate center frequency of a maximum value position of the Brillouin scattering gain spectrum according to the Brillouin center frequency curve; a Brillouin frequency shift extraction module configured to determine the center frequency of the Brillouin scattering gain spectrum based on the curve fluctuation degree of the approximate center frequency.

9. An electronic device, comprising: comprising: a processor and a memory, wherein the memory stores machine readable instructions executable by the processor, and when the electronic device is running, the machine readable instructions are executed by the processor to perform the steps of the method according to any one of claims 1 to 7.

10. A storage medium, characterized by a storage medium having a computer program stored thereon, wherein the computer program is executed by a processor to perform the steps of the method according to any one of claims 1 to 7.

Images