Automatic peak searching method, system and device for multi-peak signal and storage medium

Abstract

The invention provides an automatic peak searching method, system and device for a multi-peak signal and a storage medium so as to realize automatic peak searching in the multi-peak signal and further calculate the main peak spectral line width. The method comprises the following steps: acquiring a fluorescence signal, performing microwave scanning of different frequencies on the fluorescence signal according to a preset microwave frequency scanning range and a preset step value to obtain the number of fluorescence signals under different frequencies, drawing a fluorescence signal spectral line, and further processing the fluorescence signal spectral line; according to the processed fluorescence signal spectral lines, obtaining the difference value of the number of the fluorescence signals of every two adjacent points and the main peak frequency; calculating a first-order derivative trend of the fluorescence signals according to a difference value and a stepping value of the number of the fluorescence signals of every two adjacent points; according to the trend, determining frequencies with derivative change point values closest to the main peak frequency on the left side and the right side of the main peak frequency approaching zero, namely a first frequency and a second frequency; and extracting a main peak in the fluorescence signal spectral line according to the main peak frequency, the first frequency and the second frequency.

Description

Technical Field

[0001] The present application relates to the field of microwave frequency standards, and in particular to an automatic peak-finding method, system, device and storage medium for multi-peak signals. Background Art

[0002] The mercury ion microwave frequency standard is a new type of frequency standard that uses a new working principle that is different from the traditional atomic frequency standards of hydrogen, rubidium, and cesium lamps. It has the characteristics of being basically unaffected by physical particles and external field disturbances, small motion effects, and long quantum state coherence time. The mercury ion microwave frequency standard uses a spectral lamp to pump the ultrafine energy levels of ions. The high-energy-level atoms transition to low-energy levels after spontaneous radiation, and generate fluorescence signals at the same time. This working mode is divided into continuous mode and Rabi mode. In the continuous mode, the fluorescence spectrum line shows a single peak and does not need to be searched.

[0003] However, when working in Rabi mode, the fluorescence spectrum line shape presents a combination of a single main peak and multiple short peaks. When processing data, it is necessary to calculate the main peak spectrum line width. Therefore, it is necessary to search for the main peak data in the acquired fluorescence signal and calculate the main peak spectrum line width based on this, which will serve as the basis for subsequent experiments.

[0004] In response to the above problems, how to automatically find peaks in multi-peak signals and then calculate the width of the main peak spectral line as a basis for subsequent experiments to achieve an accurate new atomic frequency standard is a problem that needs to be solved at present. Summary of the invention

[0005] The embodiments of the present application provide a method, system, device and storage medium for automatic peak finding of multi-peak signals, so as to realize automatic peak finding in multi-peak signals, and then calculate the main peak spectral line width as a basis for subsequent experiments, and finally realize an accurate new atomic frequency standard.

[0006] In order to achieve the above objectives, this application adopts the following technical solutions:

[0007] In a first aspect, the present application provides an automatic peak-finding method for a multi-peak signal, which is used to automatically extract a main peak from a multi-peak signal. The method includes: acquiring a fluorescence signal, and performing microwave scanning of the fluorescence signal at different frequencies according to a preset microwave frequency scanning range and step value to obtain the number of fluorescence signals at different frequencies, and drawing a fluorescence signal spectrum, performing data smoothing processing on the fluorescence signal spectrum to obtain a processed fluorescence signal spectrum; according to the processed fluorescence signal spectrum, obtaining a difference in the number of fluorescence signals at each two adjacent points and a main peak frequency; according to the difference in the number of fluorescence signals at each two adjacent points and the step value, calculating a first-order derivative trend of the fluorescence signal; according to the first-order derivative trend of the fluorescence signal, determining a frequency on each side of the main peak frequency whose derivative change point value is closest to the main peak frequency and approaches zero, which are respectively a first frequency and a second frequency; and extracting a main peak in the fluorescence signal spectrum according to the main peak frequency, the first frequency and the second frequency.

[0008] In a possible design scheme, the method of the first aspect also includes that the fluorescence signal spectrum is a two-dimensional data graph, the horizontal axis is the frequency, and the vertical axis is the number of fluorescence signals.

[0009] A possible design scheme, the method of the first aspect also includes performing data smoothing processing on the fluorescence signal spectrum to obtain a processed fluorescence signal spectrum, including: based on a local weighted scatter point smoothing method, weighted smoothing is performed on the neighborhood of each fluorescence signal data point in the fluorescence signal spectrum to obtain the processed fluorescence signal spectrum, wherein the neighborhood is a preset range of adjacent data points of each fluorescence signal data point.

[0010] A possible design scheme, the method of the first aspect also includes calculating the first-order derivative trend of the fluorescence signal based on the difference and step value of the number of fluorescence signals at each two adjacent points, including: calculating the first-order derivatives of multiple fluorescence signals based on the difference and step value of the number of fluorescence signals at each two adjacent points in the fluorescence signal spectrum; obtaining the first-order derivative trend of the fluorescence signal based on the first-order derivatives of the multiple fluorescence signals.

[0011] A possible design scheme, the method of the first aspect also includes, after extracting the main peak in the fluorescence signal spectrum line, fitting the data of the main peak to obtain the function model of the main peak; based on the function model of the main peak, the first frequency, the second frequency and the main peak frequency, calculating the half-height value, and then based on the main peak function and the half-height value, calculating two frequency points, and the absolute value of the difference between the two frequency points is the main peak spectrum line width.

[0012] In a possible design scheme, the method of the first aspect further includes fitting the main peak data to obtain a function model of the main peak, including: fitting the main peak data according to the fluorescence signal, and the function model obtained by fitting is: Among them, W, S0, y0, 0 are all fitted coefficients, is the current working frequency of the microwave source, and f() represents the number of fluorescence signals.

[0013] A possible design scheme, the method of the first aspect also includes, based on the function model of the main peak, the first frequency, the second frequency and the main peak frequency, calculating the half-height value, and then calculating two frequency points based on the main peak function and the half-height value, the absolute value of the difference between the two frequency points is the main peak spectral line width, including: based on the function model of the main peak, the first frequency 1, the second frequency 2 and the main peak frequency v, calculating f(1), f(2) and f(v); by formula Calculate the half-height value, where f( h ) is the half-height value, f( max ) is f(v), f( min ) is the minimum value of f(1) and f(2); Substitute the half-height value into the function model of the main peak and calculate the two frequency points h1 , h2 , two frequency points h1 , h2 The absolute value of the difference is the spectral line width of the main peak.

[0014] In a second aspect, the present application provides an automatic peak-finding system for multi-peak signals, which is used to automatically extract a main peak from a multi-peak signal. The system includes: an acquisition module, a processing module, a calculation module and an extraction module; wherein the acquisition module is used to acquire a fluorescent signal; the processing module is used to perform microwave scanning of the fluorescent signal at different frequencies according to a preset microwave frequency scanning range and step value, obtain the number of fluorescent signals at different frequencies, and draw a fluorescent signal spectrum line, perform data smoothing processing on the fluorescent signal spectrum line to obtain a processed fluorescent signal spectrum line, and also to obtain the difference in the number of fluorescent signals at each two adjacent points and the main peak frequency based on the processed fluorescent signal spectrum line; the calculation module is used to calculate the first-order derivative trend of the fluorescent signal based on the difference in the number of fluorescent signals at each two adjacent points and the step value; and is also used to determine, based on the first-order derivative trend of the fluorescent signal, the frequencies on the left and right sides of the main peak frequency whose derivative change points closest to the main peak frequency are close to zero, which are the first frequency and the second frequency respectively; the extraction module is used to extract the main peak in the fluorescent signal spectrum line based on the main peak frequency, the first frequency and the second frequency.

[0015] In a third aspect, an automatic peak-finding device for a multi-peak signal is provided, wherein the automatic peak-finding device for a multi-peak signal comprises a module for executing the method of the first aspect.

[0016] In a possible design solution, the automatic peak-finding device for multi-peak signals of the third aspect may further include a transceiver. The transceiver may be a transceiver circuit or an interface circuit. The transceiver may be used for the automatic peak-finding device for multi-peak signals of the third aspect to communicate with other devices.

[0017] In a possible design solution, the automatic peak-finding device for multi-peak signals of the third aspect may further include a memory. The memory may be integrated with the processor or may be separately provided. The memory may be used to store instructions involved in the method of the first aspect.

[0018] In a fourth aspect, an automatic peak-finding device for multi-peak signals is provided, comprising: a processor coupled to a memory, and configured to execute instructions stored in the memory, so that the automatic peak-finding device for multi-peak signals executes the method described in the first aspect.

[0019] In a possible design solution, the automatic peak-finding device for multi-peak signals described in the fourth aspect may further include a transceiver. The transceiver may be a transceiver circuit or an interface circuit. The transceiver may be used for the automatic peak-finding device for multi-peak signals described in the fourth aspect to communicate with other devices.

[0020] In a fifth aspect, an automatic peak-finding device for multi-peak signals is provided, comprising: a processor and a memory; the memory is used to store instructions, and when the processor executes the instructions, the automatic peak-finding device for multi-peak signals executes the method described in the first aspect.

[0021] In a possible design scheme, the automatic peak-finding device for multi-peak signals described in the fifth aspect may further include a transceiver. The transceiver may be a transceiver circuit or an interface circuit. The transceiver may be used for the automatic peak-finding device for multi-peak signals described in the fifth aspect to communicate with other devices.

[0022] In a sixth aspect, a computer-readable storage medium is provided, wherein the computer-readable storage medium includes a computer program or instruction stored therein, and when the computer program or instruction is executed, the automatic peak finding method for multi-peak signals described in the first aspect is executed.

[0023] In the embodiment of the present application, the main peak position is obtained according to the first-order derivative trend, and then the main peak is fitted. Finally, the accurate spectral line width of the main peak signal is obtained according to the fitting result, which provides an accurate basis for the subsequent microwave source setting and frequency hopping setting, and helps to ultimately realize the development of high-performance new atomic frequency standards.

[0024] Other features and advantages of the present application will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0026] Figure 1 A schematic diagram of a flow chart of an automatic peak finding method for a multi-peak signal provided in an embodiment of the present application;

[0027] Figure 2 A schematic diagram of the structure of the automatic peak-finding device for multi-peak signals provided in an embodiment of the present application Figure 1 ;

[0028] Figure 3 A schematic diagram of the structure of the automatic peak-finding device for multi-peak signals provided in an embodiment of the present application Figure 2 . DETAILED DESCRIPTION

[0029] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without making creative work are within the scope of protection of the present application. At the same time, in the description of the embodiments of the present application, the terms "first", "second", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance. Thus, the features defined as "first" and "second" may explicitly or implicitly include one or more features. In the description of the embodiments of the present application, the meaning of "multiple" is two or more, unless otherwise clearly and specifically defined.

[0030] Figure 1 A flow chart of the automatic peak finding method for multi-peak signals provided in the embodiment of the present application. The automatic peak finding method for multi-peak signals is applicable to mercury ion microwave frequency standard, when the Rabi mode is used for fluorescence signal detection, when it is necessary to find the main peak position of the signal, and then calculate the main peak spectral line width, or other possible scenarios, which are not limited here. To facilitate the understanding of this application, the following mainly introduces the mercury ion microwave frequency standard scenario as an example.

[0031] The process of the automatic peak finding method of the multi-peak signal is as follows:

[0032] Step S101, obtaining a fluorescence signal, and performing microwave scanning of different frequencies on the fluorescence signal according to a preset microwave frequency scanning range and step value, obtaining the number of fluorescence signals at different frequencies, and drawing a fluorescence signal spectrum, performing data smoothing processing on the fluorescence signal spectrum, and obtaining a processed fluorescence signal spectrum.

[0033] It can be understood that after the mercury ions are stably trapped in the ion trap, they are prepared by the mercury lamp to make them transition from the ground state to the excited state, and then they spontaneously radiate from the excited state to the ground state. 2 S 1 / 2 After the state preparation is completed, the ground state hyperfine structure F = 1 state is evacuated, and finally almost all ions fall into 2 S 1 / 2 F=0 state; at this time, 40.5GHz microwaves are applied (this frequency is the ground state 2 S 1 / 2 The ions undergo resonant absorption and transition from the F=0 state to the F=1 state, increasing the number of ions in the F=1 state and pumping them back to the excited state. At this time, the spontaneous fluorescence becomes stronger, and the fluorescence signal is obtained by using a photomultiplier tube (PMT) in conjunction with the collection optical path. The signal intensity here is very weak relative to the noise floor. If the voltage data of the photomultiplier tube is used as the basis for the fluorescence signal conversion, the signal will be submerged. Therefore, single photon counting is required for collection, and a certain integration time is set to perform photon counting to obtain the fluorescence signal.

[0034] After that, the microwave frequency scanning range is preset, the starting point of the range is, the stop point is, the attenuation Att, the step value s, and after starting the microwave source, the current frequency and the photon count value (that is, the number of fluorescence signals at the current frequency) are recorded through the self-compiled software. In addition, since the microwave frequency will be scanned and changed according to the step value, the number of ions transitioning from the F=0 state to the F=1 state is different, so the fluorescence signal value obtained is also different. When the microwave source frequency changes to the stop point, the microwave scanning stops.

[0035] After the scanning is completed, a fluorescence signal spectrum is drawn. The fluorescence signal spectrum is a two-dimensional data graph, where the horizontal axis is the frequency and the vertical axis is the number of fluorescence signals.

[0036] Since the fluorescence signal spectrum line has a lot of data fluctuations, the fluorescence signal spectrum line is smoothed to eliminate the fluctuations, and the short peaks beside the main peak can be further depressed, which is more conducive to extracting the main peak. For example: based on the local weighted scatter point smoothing method, the neighborhood of each fluorescence signal data point in the fluorescence signal spectrum line is weighted and smoothed to obtain the processed fluorescence signal spectrum line, where the neighborhood is the preset neighboring data point range of each fluorescence signal data point. Alternatively, the Savitzky-Golay filter, wavelet transform and other methods are used to reduce noise and smooth data, which are not limited here.

[0037] Step S102, obtaining the difference in the number of fluorescence signals between every two adjacent points and the main peak frequency according to the processed fluorescence signal spectrum.

[0038] For example, for 51 sets of fluorescence signal data, microwave frequencies and fluorescence signal intensities are in one-to-one correspondence, as shown below. The microwave frequencies are: 0.0, 0.04, 0.08, 0.12, 0.16, 0.20, 0.24, 0.28, 0.32, 0.36, 0.40, 0.44, 0.48, 0.52, 0.56, 0.6, 0.64, 0.68, 0.72, 0.76, 0.8, 0.84, 0.88, 0.92, 0.96, 1.0, 1.04, 1.08, 1.12, 1.16, 1.20 , 1.24, 1.28, 1.32, 1.36, 1.4, 1.44, 1.48, 1.52, 1.56, 1.6, 1.64, 1.68, 1.72, 1.76, 1.80, 1.84, 1.88, 1.92, 1.96, 2. The fluorescence signal intensities corresponding to the above microwave frequencies are: 952000, 951000, 946000, 948000, 950000, 950000 0, 947000, 949000, 948000, 952000, 953000, 951500, 951000, 950500, 950600, 946000, 946000, 946000, 952500, 955000, 962000, 965000, 970000, 976000, 980000, 979000, 981000, 977000, 97 6000, 970000, 965000, 960000, 955000, 953000, 950000, 948000, 949000, 951000, 953000, 951500, 952900, 953200, 950000, 950000, 948000, 950000, 949000, 951000, 950500, 952000, 951500.

[0039] It can be seen that the main peak appears at frequency 0 = 1.04, and the number of fluorescent signals is 981000. The difference in the number of fluorescent signals at each two adjacent points can be calculated one by one. For example, when the microwave frequencies are 0.0 and 0.04, the corresponding difference in the number of fluorescent signals is 951000 minus 952000, which is -1000, and so on.

[0040] Step S103, calculating the first-order derivative trend of the fluorescence signal according to the difference in the number of fluorescence signals between every two adjacent points and the step value.

[0041] The calculation method of the first-order derivative is: Among them, is the step value s, and is the difference in the number of fluorescence signals between every two adjacent points in the fluorescence signal spectrum.

[0042] That is to say, the first-order derivatives of multiple fluorescence signals are calculated based on the difference and step value of the fluorescence signals of each two adjacent points. The calculation formula is: Among them, y n+1 y n is the difference in the number of fluorescence signals between every two adjacent fluorescence signal data points in the fluorescence signal spectrum, n is a positive integer greater than 1; s is the step value.

[0043] According to the first-order derivatives of the multiple fluorescence signals, a first-order derivative trend of the fluorescence signal is obtained.

[0044] Corresponding to the example of step S102, according to the first-order derivative calculation formula and the preset step value s of 0.04, 50 first-order derivatives can be calculated, namely: -25000, -125000, etc., and so on.

[0045] Step S104, according to the first-order derivative trend of the fluorescence signal, determine the frequencies on the left and right sides of the main peak frequency whose derivative change points closest to the main peak frequency are close to zero, which are the first frequency and the second frequency respectively.

[0046] That is to say, when analyzing the fluorescence signal spectrum, we look for a derivative change point on each side of the main peak, which is closest to the main peak. The values ​​of these change points are close to 0 because these points represent the change in the slope of the spectral curve, that is, the position where the concavity of the curve changes. These points are usually the points where the first-order derivative of the function is zero, that is, the extreme value of the function, that is, the position of the peak or trough.

[0047] That is, when the first-order derivative of the fluorescence signal spectrum is calculated, at the peak of the spectrum, the slope of the curve (i.e., the first-order derivative) is positive on the left side of the peak and negative on the right side of the peak. At the peak point, the slope changes from positive to negative, so the first-order derivative passes through zero at the peak point. Similarly, at the trough, the slope changes from negative to positive, and the first-order derivative is also zero.

[0048] Corresponding to the example of step S103, through step S104, it can be determined that the first frequency and the second frequency are 0.6 and 1.4 respectively.

[0049] Step S105, extracting the main peak in the fluorescence signal spectrum according to the main peak frequency, the first frequency and the second frequency.

[0050] The main peak can be extracted according to three key points: the first frequency 1, the second frequency 2 and the main peak frequency v. That is, the main peak in the example of this application can be extracted by frequencies 0.6, 1.4 and 1.04. In practical applications, it may also be necessary to use curve fitting techniques, such as Gaussian fitting or Lorentz fitting, to extract the detailed shape of the main peak from the raw data.

[0051] In addition, after step S105, the automatic peak finding method for multi-peak signals may further include steps S106 and S107 to calculate the main peak spectral line width.

[0052] Step S106, fitting the data of the main peak to obtain a function model of the main peak.

[0053] That is to say, according to the fluorescence signal, the data of the main peak is fitted, and the function model obtained by fitting is:

[0054]

[0055] Among them, the data of the main peak are mainly multiple data points in the aforementioned fluorescence signal spectrum corresponding to the main peak, W, S0, y0, 0 are all fitted coefficients, is the current working frequency of the microwave source, and f() represents the number of fluorescence signals.

[0056] Corresponding to the example in the above steps, in the 51 sets of fluorescence signal data, the fitted coefficients are: W is -0.536087, S0 is -27598.3, y0 is 929570, and 0 is 1.01983.

[0057] Step S107, based on the function model of the main peak, the first frequency, the second frequency and the main peak frequency, the half-height value is calculated, and then based on the main peak function and the half-height value, two frequency points are calculated, and the absolute value of the difference between the two frequency points is the main peak spectral line width.

[0058] Based on the function model of the main peak, the first frequency 1, the second frequency 2 and the main peak frequency v, f(1), f(2) and f(v) are calculated; by the formula Calculate the half-height value, where f( h ) is the half-height value, f( max ) is f(v), f( min ) is the minimum value of f(1) and f(2).

[0059] Substitute the half-height value into the function model of the main peak and calculate two frequency points h1 , h2 .

[0060] Two frequency points h1 , h2 The absolute value of the difference is the line width of the main peak, that is, W =| 1h2 |.

[0061] Combining the above examples, we can calculate f( h ) is 962624, and substituting it into the function model of the main peak, we can get h1 is 0.800129, h2It is 1.19987, and then the spectral line width of the main peak is determined to be 0.399741.

[0062] In summary, the embodiment of the present application obtains the main peak position according to the first-order derivative trend, then fits the main peak, and finally obtains the accurate spectral line width of the main peak signal according to the fitting result, which provides an accurate basis for subsequent microwave source setting and frequency hopping setting, and helps to ultimately realize the development of high-performance new atomic frequency standards.

[0063] Combination of the above Figure 1 The automatic peak finding method for a multi-peak signal provided in an embodiment of the present application is described in detail. The automatic peak finding system for executing the multi-peak signal provided in an embodiment of the present application is described in detail below.

[0064] The system specifically includes: an acquisition module, a processing module, a calculation module and an extraction module, as shown below.

[0065] The acquisition module is used to acquire the fluorescence signal.

[0066] The processing module is used to perform microwave scanning of different frequencies on the fluorescent signal according to a preset microwave frequency scanning range and step value, obtain the number of fluorescent signals at different frequencies, draw fluorescent signal spectra, perform data smoothing on the fluorescent signal spectra, obtain processed fluorescent signal spectra, and obtain the difference in the number of fluorescent signals at each two adjacent points and the main peak frequency based on the processed fluorescent signal spectra.

[0067] The calculation module is used to calculate the first-order derivative trend of the fluorescence signal according to the difference and step value of the number of fluorescence signals at each two adjacent points; it is also used to determine the frequencies of the derivative change points closest to the main peak frequency on both sides of the main peak frequency, which are close to zero, based on the first-order derivative trend of the fluorescence signal, which are the first frequency and the second frequency, respectively.

[0068] The extraction module is used to extract the main peak in the fluorescence signal spectrum according to the main peak frequency, the first frequency and the second frequency.

[0069] In addition, for the above-mentioned specific implementation of the system, since it is basically similar to the method implementation, the description is relatively simple, and the relevant parts can be referred to the partial description of the method implementation. Moreover, it should be noted that in each module of the system of the present application, the components therein are logically divided according to the functions to be implemented, but the present application is not limited thereto, and the components can be re-divided or combined as needed.

[0070] The above describes the automatic peak finding method and system for multi-peak signals provided by the embodiments of the present application. Figure 2-Figure 3 The automatic peak-finding device for executing the multi-peak signal provided in the embodiment of the present application is described in detail.

[0071] Figure 2 The structure of the automatic peak-finding device for multi-peak signals provided in the embodiment of the present application is shown in FIG. Figure 1 For example, Figure 2 As shown, the automatic peak-finding device 200 for multi-peak signals includes: a transceiver module 201 and a processing module 202. For ease of description, Figure 2 Only the main components of the automatic peak finding device for the multi-peak signal are shown.

[0072] The transceiver module 201 is used to perform the transceiver function of the automatic peak finding method for multi-peak signals, and the processing module 202 is used to perform other functions of the automatic peak finding method for multi-peak signals except the transceiver function.

[0073] Optionally, the transceiver module 201 may include a sending module ( Figure 2 ) and a receiving module ( Figure 2 (not shown in the figure). The sending module is used to implement the sending function of the automatic peak-finding device 200 for multi-peak signals, and the receiving module is used to implement the receiving function of the automatic peak-finding device 200 for multi-peak signals.

[0074] Optionally, the automatic peak-finding device 200 for multi-peak signals may further include a storage module ( Figure 2 (not shown), the storage module stores a program or instruction. When the processing module 202 executes the program or instruction, the automatic peak-finding device 200 for multi-peak signals can execute the automatic peak-finding method for multi-peak signals in the embodiment of the present application.

[0075] Combine the following Figure 3 Each component of the automatic peak-finding device 300 for multi-peak signals is specifically introduced:

[0076] The processor 301 is the control center of the automatic peak-finding device 300 for multi-peak signals, and may be a processor or a general term for multiple processing elements. For example, the processor 301 is one or more central processing units (CPUs), or may be application specific integrated circuits (ASICs), or may be configured to implement one or more integrated circuits of the embodiments of the present application, such as one or more microprocessors (digital signal processors, DSPs), or one or more field programmable gate arrays (FPGAs).

[0077] Optionally, the processor 301 can execute various functions of the automatic peak finding device 300 for multi-peak signals, such as executing the automatic peak finding method for multi-peak signals in the embodiment of the present application, by running or executing a software program stored in the memory 302 and calling data stored in the memory 302.

[0078] In a specific implementation, as an embodiment, the processor 301 may include one or more CPUs, such as Figure 3 CPU0 and CPU1 are shown in FIG.

[0079] In a specific implementation, as an embodiment, the automatic peak-finding device 300 for multi-peak signals may also include multiple processors, such as Figure 3 The processor 301 and processor 304 shown in . Each of these processors can be a single-core processor (single-CPU) or a multi-core processor (multi-CPU). The processor here can refer to one or more devices, circuits, and / or processing cores for processing data (such as computer program instructions). Among them, the memory 302 is used to store the software program for executing the solution of the present application, and the execution is controlled by the processor 301. The specific implementation method can refer to the above method embodiment, which will not be repeated here.

[0080] Optionally, the memory 302 may be a read-only memory (ROM) or other types of static storage devices that can store static information and instructions, a random access memory (RAM) or other types of dynamic storage devices that can store information and instructions, or an electrically erasable programmable read-only memory (EEPROM), a compact disc read-only memory (CD-ROM) or other optical disc storage, optical disc storage (including compressed optical disc, laser disc, optical disc, digital versatile disc, Blu-ray disc, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto. The memory 302 may be integrated with the processor 301, or may exist independently, and may be connected to the processor 301 through the interface circuit ( Figure 3 (not shown) is coupled to the processor 301, which is not specifically limited in the embodiment of the present application.

[0081] The transceiver 303 is used for communication with other communication devices. For example, the automatic peak-finding device 300 for multi-peak signals is a first device, and the transceiver 303 can be used for communication with a second device or a third device.

[0082] Optionally, the transceiver 303 may include a receiver and a transmitter ( Figure 3 The receiver is used to implement a receiving function, and the transmitter is used to implement a sending function.

[0083] Optionally, the transceiver 303 may be integrated with the processor 301, or may exist independently, and communicate with the processor 301 through the interface circuit ( Figure 3 (not shown) is coupled to the processor 301, which is not specifically limited in the embodiment of the present application.

[0084] Understandably, Figure 3 The structure of the automatic peak-finding device 300 for multi-peak signals shown in the figure does not constitute a limitation of the automatic peak-finding device for multi-peak signals. An actual automatic peak-finding device for multi-peak signals may include more or fewer components than those shown in the figure, or combine certain components, or arrange the components differently.

[0085] In addition, the technical effects of the automatic peak-finding device 300 for multi-peak signals can refer to the technical effects of the methods described in the above method embodiments, which will not be described in detail here.

[0086] It should be understood that the processor in the embodiments of the present application may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0087] It should also be understood that the memory in the embodiments of the present application may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. Among them, the non-volatile memory may be 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), or a flash memory. The volatile memory may be a random access memory (RAM), which is used as an external cache. By way of example and not limitation, many forms of random access memory (RAM) are available, such as static RAM (SRAM), dynamic random access memory (DRAM), synchronous DRAM (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link DRAM (SLDRAM), and direct rambus RAM (DR RAM).

[0088] The above embodiments can be implemented in whole or in part by software, hardware (such as circuits), firmware or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions can be transmitted from one website site, computer, server or data center to another website site, computer, server or data center by wired (such as infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that contains one or more available media sets. The available medium can be a magnetic medium (for example, a floppy disk, a hard disk, a tape), an optical medium (for example, a DVD), or a semiconductor medium. The semiconductor medium can be a solid-state hard disk.

Claims

1. An automatic peak finding method for a multi-peak signal, used for automatically extracting a main peak from a multi-peak signal, characterized in that: The method comprises: Acquire a fluorescence signal, and perform microwave scanning of different frequencies on the fluorescence signal according to a preset microwave frequency scanning range and step value to obtain the number of fluorescence signals at different frequencies, draw a fluorescence signal spectrum, and perform data smoothing processing on the fluorescence signal spectrum to obtain a processed fluorescence signal spectrum; According to the processed fluorescence signal spectrum, the difference in the number of fluorescence signals at each two adjacent points and the main peak frequency are obtained; Calculating the first-order derivative trend of the fluorescence signal according to the difference in the number of fluorescence signals of each two adjacent points and the step value; According to the first-order derivative trend of the fluorescence signal, determine the frequencies on the left and right sides of the main peak frequency, the frequencies whose derivative change points closest to the main peak frequency are close to zero, which are the first frequency and the second frequency respectively; A main peak in the fluorescence signal spectrum is extracted according to the main peak frequency, the first frequency and the second frequency.

2. The automatic peak finding method for multi-peak signals according to claim 1, characterized in that: The fluorescence signal spectrum is a two-dimensional data graph, with the horizontal axis being the frequency and the vertical axis being the number of fluorescence signals.

3. The automatic peak finding method for multi-peak signals according to claim 1, characterized in that: The step of performing data smoothing on the fluorescence signal spectrum to obtain a processed fluorescence signal spectrum includes: Based on the local weighted scatter point smoothing method, the neighborhood of each fluorescence signal data point in the fluorescence signal spectrum is weighted and smoothed to obtain a processed fluorescence signal spectrum, wherein the neighborhood is a preset range of adjacent data points of each fluorescence signal data point.

4. The automatic peak finding method for multi-peak signals according to claim 1, characterized in that: The step of calculating the first-order derivative trend of the fluorescence signal according to the difference in the number of fluorescence signals of each two adjacent points and the step value comprises: Calculating a plurality of first-order derivatives of the fluorescence signal according to the difference in the number of fluorescence signals of each two adjacent points in the fluorescence signal spectrum and the step value; A first-order derivative trend of the fluorescence signal is obtained according to the first-order derivatives of the plurality of fluorescence signals.

5. The automatic peak finding method for multi-peak signals according to claim 1, characterized in that: After extracting the main peak in the fluorescence signal spectrum, the method further comprises: Fitting the data of the main peak to obtain a function model of the main peak; Based on the function model of the main peak, the first frequency, the second frequency and the main peak frequency, a half-height value is calculated, and then based on the main peak function and the half-height value, two frequency points are calculated, and the absolute value of the difference between the two frequency points is the main peak spectral line width.

6. The automatic peak finding method for multi-peak signals according to claim 5, characterized in that: The fitting of the main peak data to obtain a function model of the main peak includes: According to the fluorescence signal, the main peak data is fitted, and the function model obtained by fitting is: Among them, W, S0, y0, 0 are all fitted coefficients, is the current working frequency of the microwave source, and f() represents the number of fluorescence signals.

7. The automatic peak finding method for multi-peak signals according to claim 5, characterized in that: The function model based on the main peak, the first frequency, the second frequency and the main peak frequency is used to calculate the half-height value, and then based on the main peak function and the half-height value, two frequency points are calculated, and the absolute value of the difference between the two frequency points is the main peak spectral line width, including: Based on the function model of the main peak, the first frequency 1, the second frequency 2 and the main peak frequency v, f(1), f(2) and f(v) are calculated; By formula Calculate the half-height value, where f( h ) is the half-height value, f( max ) is f(v), f( min ) is the minimum value of f(1) and f(2); Substitute the half-height value into the function model of the main peak and calculate two frequency points h1 , h2 , The two frequency points h1 , h2 The absolute value of the difference is the line width of the main peak.

8. An automatic peak-finding system for multi-peak signals, used for automatically extracting a main peak from a multi-peak signal, characterized in that: The system comprises: an acquisition module, a processing module, a calculation module and an extraction module; The acquisition module is used to acquire the fluorescence signal; The processing module is used to perform microwave scanning of different frequencies on the fluorescent signal according to a preset microwave frequency scanning range and step value to obtain the number of fluorescent signals at different frequencies, draw fluorescent signal spectrum lines, perform data smoothing processing on the fluorescent signal spectrum lines to obtain processed fluorescent signal spectrum lines, and further used to obtain the difference in the number of fluorescent signals at each two adjacent points and the main peak frequency according to the processed fluorescent signal spectrum lines; The calculation module is used to calculate the first-order derivative trend of the fluorescence signal according to the difference in the number of fluorescence signals of each two adjacent points and the step value; and is also used to determine, according to the first-order derivative trend of the fluorescence signal, the frequencies of the derivative change points closest to the main peak frequency on both sides of the main peak frequency, which are close to zero, as the first frequency and the second frequency, respectively; The extraction module is used to extract the main peak in the fluorescence signal spectrum according to the main peak frequency, the first frequency and the second frequency.

9. An automatic peak-finding device for multi-peak signals, characterized in that: The apparatus comprises: a module for executing the method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a computer program or instruction stored therein, and when the computer program or instruction is executed, the method according to any one of claims 1 to 7 is performed.

Images