An electronic aerosol composition detection method and device

Abstract

The present invention discloses a method and device for detecting components of an electronic aerosol, belonging to the technical field of detecting components by gas chromatography. When segmenting the chromatographic peaks corresponding to each known component, the length of the segmentation window is adjusted in increments of the secondary peak, which is more reasonable and can improve the accuracy of the segmentation results of the chromatographic peaks of each component. And because the segmentation window is adjusted in increments of the secondary peak, the number of adjustments is only the number of secondary peaks, reducing the number of times of adjusting the segmentation window and improving the data processing efficiency of determining the length of the segmentation window.

Description

Technical Field

[0001] The present invention relates to the technical field of detecting components by gas chromatography, and particularly relates to a method and device for detecting components of an electronic aerosol. Background Art

[0002] In the process of detecting the components of a vape juice aerosol using a gas chromatograph, due to fluctuations in the column temperature and the carrier gas flow rate of the chromatograph, a state where the same component corresponds to one main peak and multiple sub-peaks simultaneously (hereinafter referred to as the multi-peak state) is caused. Therefore, as described in Patent Document 1 (publication number: CN118501325A), when dividing the corresponding retention time periods for each component in the chromatogram, it is necessary to first identify one main peak and multiple sub-peaks corresponding to the same component, and then merge each peak. Based on the merging result, the length range of the window is determined; then a window with a random length within this length range is selected to calculate parameters (parameters including evaluation factors, initial attribution degrees, overlapping characteristics, and final attribution degrees, etc.) with the retention time periods of each known component. According to the calculation results, the window length is continuously adjusted until a window with the highest degree of conformity with the retention time period of the known component is obtained, and the chromatogram is segmented using this window.

[0003] However, the solution in Patent Document 1 still has the following deficiencies:

[0004] 1. When segmenting, the window length is randomly selected. Therefore, it is necessary to continuously adjust windows of different lengths to perform multiple parameter calculations before finally obtaining a window with the highest degree of conformity with the retention time period of the known component for segmentation. That is, the parameter calculation process needs to be repeated for each window length, resulting in a low efficiency of the data processing flow for determining the final window length.

[0005] 2. Moreover, during the process of adjusting the window length, the increment that restricts the window length adjustment is a unit time interval preset manually. If this unit time interval is set unreasonably, it will also reduce the accuracy of the segmentation results of the chromatographic peaks of each component. Summary of the Invention

[0006] An object of the present invention is to at least solve one of the technical problems existing in the prior art, and provide a method and device for detecting components of an electronic aerosol, which can improve the efficiency of data processing and enhance the accuracy of the segmentation results of the chromatographic peaks of each component.

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

[0008] A method for detecting components of an electronic aerosol, the method comprising:

[0009] Obtaining each maximum point and each minimum point within a preset range on the chromatogram of the aerosol sample to be detected;

[0010] The coordinates of each maximum point, the left minimum point adjacent to each maximum point, and the right minimum point within a preset range on the chromatogram are jointly used as the peak descriptor corresponding to a maximum point; wherein, the abscissa of the maximum point in the peak descriptor is used as the position of the peak descriptor.

[0011] According to the forward order of the abscissas of each maximum point within a preset range on the chromatogram, the peak descriptors corresponding to each maximum point are sorted to obtain a first sequence.

[0012] For each known component, based on the preset first window of the known component, the sequence segment corresponding to the peak descriptor whose position is within the first window in the first sequence is used as the second sequence corresponding to the known component; wherein, the range of the first window is determined according to the retention time periods corresponding to each known component in the component detection results of other detected aerosol samples of the same specification as the aerosol sample to be detected.

[0013] The peak descriptor with the smallest sort order in the second sequence is used as the left peak descriptor, and the peak descriptor with the largest sort order in the second sequence is used as the right peak descriptor.

[0014] For each known component, the starting point of the retention time period of the known component is determined based on the left minimum point of the left peak descriptor, and the end point of the retention time period of the known component is determined based on the right minimum point of the right peak descriptor, so as to obtain the retention time period corresponding to each known component.

[0015] Based on the retention time periods corresponding to each known component, the amount of the known component is detected.

[0016] On the other hand, the present invention also provides an aerosol component detection device, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, it executes the steps of the above-mentioned aerosol component detection method.

[0017] Advantageous effects: Compared with the method of segmenting the retention time periods corresponding to each component in Patent Document 1, in the aerosol component detection method and device provided by the present invention, when segmenting the chromatographic peaks corresponding to each known component, the length of the segmentation window is adjusted in increments of the sub-peak, which is more reasonable and can improve the accuracy of the segmentation results of the chromatographic peaks of each component. And because the segmentation window is adjusted in increments of the sub-peak, the number of adjustments is only the number of sub-peaks, reducing the number of times of adjusting the segmentation window and improving the data processing efficiency of determining the length of the segmentation window.

[0018] The additional aspects and advantages of the present invention will be partially given in the following description, partially become obvious from the following description, or be understood through the practice of the present invention. Brief Description of the Drawings

[0019] The present invention will be further described below in conjunction with the drawings and embodiments;

[0020] Figure 1 Schematic diagram of the first window range corresponding to some known components on the chromatogram in one embodiment;

[0021] Figure 2 Schematic diagram of the second window range corresponding to some known component groups on the chromatogram in one embodiment;

[0022] Figure 3 Schematic diagram of the state when the maximum and minimum values in the acquisition window are in different positional relationships in one embodiment;

[0023] Figure 4 Schematic diagram of the principle of searching for the maximum value point and the minimum value point within the current acquisition window range in one embodiment;

[0024] Figure 5 Schematic diagram of the flow of a method for detecting the components of an electronic aerosol in one embodiment;

[0025] Figure 6 Schematic diagram of the flow for verifying the detection result in one embodiment;

[0026] Figure 7 Schematic diagram of the flow for determining abnormal component information in one embodiment;

[0027] Figure 8 Schematic diagram of the flow for obtaining the maximum value point and the minimum value point in one embodiment;

[0028] Figure 9 Schematic diagram of the flow for determining the indication code corresponding to the current acquisition window in one embodiment. Detailed Description of the Embodiments

[0029] This part will describe in detail the specific embodiments of the present invention. The preferred embodiments of the present invention are shown in the drawings. The role of the drawings is to supplement the description in the text part of the specification, enabling people to intuitively and vividly understand each technical feature and the overall technical solution of the present invention, but it should not be construed as a limitation on the protection scope of the present invention.

[0030] It should be noted that during chromatographic separation, each component has a different retention time according to its different interactions with the stationary phase in the chromatographic column. The efficiency of chromatographic separation strongly depends on the uniformity of the column temperature. Changes in temperature will change the interaction force between the sample and the stationary phase, thereby affecting the retention time and separation behavior. At the same time, chromatographic separation also depends on the distribution equilibrium between the stationary phase and the mobile phase. Changes in the carrier gas flow rate will affect the migration time of substances between the stationary phase and the mobile phase. Therefore, during the detection of the components of e-liquid aerosols using a gas chromatograph, due to fluctuations in the chromatographic column temperature and the carrier gas flow rate, a state where the same component corresponds to a main peak and multiple sub-peaks simultaneously (hereinafter referred to as the multi-peak state) occurs. Specifically, as Figure 1 shown, in the interval of the horizontal axis from t3 to t4, there are four chromatographic peaks corresponding to the known component B at the same time. Among them, the peak with the largest peak height is defined as the main peak, and the other three peaks with lower peak heights are defined as sub-peaks. In an ideal state, the chromatographic peak of a component is in a single-peak state, Figure 1 and the chromatographic peak corresponding to the known component A shown in the interval of the horizontal axis from t1 to t2 in is a single peak. However, due to the limited control accuracy of the chromatographic column temperature and the carrier gas flow rate, it may lead to a situation where for a certain component during this detection, due to the ideal control effect of the chromatographic column temperature and the carrier gas flow rate, its corresponding chromatographic peak is a single peak. However, in the next inspection, although the same set of chromatographs and the same chromatograph control parameters are used, due to the limited control accuracy of the chromatographic column temperature and the carrier gas flow rate, fluctuations in the chromatographic column temperature and the carrier gas flow rate may occur, which may lead to the chromatographic peak of this component being in a multi-peak state. For known components, under the same chromatograph control parameters, their retention times can generally be determined by measuring standard substances. For example, if the retention time of a certain chromatographic peak detected in the sample is the same as the retention time of the chromatographic peak in the nicotine standard solution, for example, within plus or minus 10 seconds of 32 minutes, it is considered that the component corresponding to this chromatographic peak is nicotine. The area is obtained by integrating the chromatographic peak, and then the amount of nicotine is determined according to the pre-calibrated relationship between the area and the amount of substance. For example, the relationship between the area and the amount of substance can be determined by the area internal standard method. Specifically, take the standard component to be measured, and add it to the specified quantitative internal standard substances of each monomer according to the known stepwise amounts that increase or decrease in turn to prepare a standard solution. Take a certain amount of this standard solution and inject it into the chromatographic column of the gas chromatograph. Using the ratio of the peak area and peak height of the standard component to be measured and the peak area and peak height of the internal standard substance on the chromatogram of the gas chromatograph as the ordinate, and taking the ratio of the amount of the standard component to be measured and the amount of the internal standard substance, or the amount of the standard component to be measured as the abscissa, a standard curve is made. Then prepare the sample solution according to the method specified in the monomer. When preparing the sample solution, add the same amount of internal standard substance as when preparing the standard solution in advance. Then, according to the chromatogram obtained under the same conditions as when making the standard curve, find the ratio of the peak area or peak height of the component to be measured and the peak area or peak height of the internal standard substance, and then find the content of the component to be measured according to the standard curve.

[0031] Referring to Patent Document 2 (publication number: CN111595992A), for the single-peak state, two minimum points can be searched from the maximum value of the chromatographic peak to the left and right sides of the horizontal axis respectively, and then the starting point and ending point (i.e., the retention time period) of the chromatographic peak are determined by the left and right minimum points, and the retention time period is integrated to obtain the amount of substance of this component. For the multi-peak state, the retention time periods corresponding to the chromatographic peaks in the multi-peak state of each known component can be determined by an expert through the method of manually calibrating a small number of samples in the same batch of samples of the same specification; for example, the secondary peaks on both sides of the main peak are manually assigned to the corresponding components, and the difference between the area of this component under different assignment methods and the amount of substance of this component in the single-peak state in the same-specification samples is calculated, and the assignment method when the difference in the amount of substance meets the allowable amount-of-substance error of this component in production is selected to determine the retention time period corresponding to this component in the multi-peak state. If the chromatographic peak of known component B in a sample of this specification is in the single-peak state during chromatograph detection and the calculated amount of substance is 50 units, and the allowable amount-of-substance error of this component in production is 3%, and in another detection, the chromatographic peak of known component B in the same-specification sample forms a multi-peak state, when the two secondary peaks on the left side of the main peak and the single secondary peak on the right side of known component B are manually assigned to the chromatographic peak corresponding to known component B, and the left minimum value of the leftmost secondary peak is used as the starting point and the right minimum value of the rightmost secondary peak is used as the ending point to determine the retention time period of known component B in the multi-peak state, and the calculated amount of substance of known component B obtained under this assignment method is 52 units, with the error within 3%, and other assignment methods do not meet the 3% error, so the retention time period corresponding to this component B can be determined. By obtaining the retention time periods corresponding to multiple known component Bs through multiple samples of the same specification, the average value of the starting points corresponding to the multiple retention time periods of known component B obtained through the manual method is used as the left boundary t3 of the first window corresponding to known component B, and similarly, the average value of the ending points is obtained as the right boundary t4. Since the fluctuation degree of the amount of substance of the components in the same batch of samples of the same specification is low, the retention time periods corresponding to each known component in the multi-peak state determined manually are used as a reference subsequently, so as to determine the amount of substance of each known component corresponding in the multi-peak state. For the scheme of detecting each known component in the multi-peak state, Patent Document 1 can be referred to.

[0032] As described in Patent Document 1, when dividing the corresponding retention time periods for each component in a chromatogram, it is necessary to first identify that the same component corresponds to a main peak and multiple secondary peaks, and then merge each peak. According to the merging result, determine the length range of the window; then select a window with a random length within this length range to calculate parameters (parameters include evaluation factors, initial attribution degrees, overlapping characteristics, and final attribution degrees, etc.) with the retention time periods of each known component. According to the calculation results, continuously adjust the window length until finally obtaining a window with the highest degree of compliance with the retention time periods of the known components, and segment the chromatogram with this window.

[0033] However, the solution in Patent Document 1 still has the following deficiencies:

[0034] 1. When segmenting, the window length is randomly selected. Therefore, it is necessary to continuously adjust windows of different lengths to perform multiple parameter calculations. Finally, a window with the highest degree of compliance with the retention time periods of the known components can be obtained for segmentation. That is, the parameter calculation process needs to be repeated for each length of the window, resulting in a low efficiency of the data processing flow for determining the final length of the window.

[0035] 2. And during the process of adjusting the window length, the increment that restricts the window length adjustment is a manually preset unit time interval. If this unit time interval is set unreasonably, it will also reduce the accuracy of the segmentation results of each component chromatographic peak.

[0036] 3. In the multi-peak state, in order to determine the length range of the window, it is also necessary to first perform the process of merging the main peak and multiple secondary peaks, and the data processing efficiency is limited.

[0037] Based on this, the present invention provides an electronic aerosol component detection method, which can improve the efficiency of data processing and enhance the accuracy of the segmentation results of each component chromatographic peak.

[0038] The following will combine the accompanying drawings and embodiments to detail an electronic aerosol component detection method provided by the present invention.

[0039] As Figure 5 shown, in one embodiment, an electronic aerosol component detection method is provided, and the method includes:

[0040] Step S202, obtain each maximum point and each minimum point within a preset range on the chromatogram of the aerosol sample to be detected.

[0041] Specifically, each chromatogram is a voltage-time graph generated by different voltage signals shown in the detector when an e-cigarette oil (i.e., an aerosol sample) passes through the end of the chromatographic column. That is, the vertical axis is voltage and the horizontal axis is time. In one example, each chromatogram corresponding to an aerosol sample has a total of 8,000 data points, arranged in sequence starting from 0. Every 1,000 data points represent one minute, and each chromatogram corresponding to an aerosol sample is data for the entire 8 minutes. In one example, the preset range on the chromatogram can be the entire range from 0 to 8,000 on the horizontal axis. In one example, the preset range on the chromatogram can be the range corresponding to the left boundary to the right boundary of each first window. As Figure 1 shown, the preset range on the chromatogram can be three intervals such as from t1 to t2, from t3 to t4, and from t5 to t6.

[0042] Step S204: Use the coordinates of each maximum point, the left minimum point adjacent to each maximum point, and the right minimum point on the chromatogram within the preset range as the peak descriptor corresponding to one maximum point. Among them, use the abscissa of the maximum point in the peak descriptor as the position of the peak descriptor.

[0043] It should be noted that the left minimum point adjacent to the maximum point is the first minimum point encountered by the maximum point along the negative direction of the horizontal axis (i.e., the direction in which the value of the horizontal axis decreases), and the right minimum point adjacent to the maximum point is the first minimum point encountered by the maximum point along the positive direction of the horizontal axis (i.e., the direction in which the value of the horizontal axis increases). The peak descriptor records various information due to constituting a peak. The maximum point records the peak height of the peak, and the distance between the left minimum point and the right minimum point on the horizontal axis is the peak width. Use the abscissa value of the maximum value in the peak descriptor as the position of the peak descriptor. For example, if the abscissa of the maximum value is t7, then the position of the peak descriptor is t7.

[0044] It can be seen that each data point is recorded in the form of a coordinate. For example, in one example, the maximum point of a certain peak description is recorded as (t9, U9), the left minimum point is recorded as (t8, U8), and the right minimum point is recorded as (t10, U10). This peak descriptor is a six-tuple and is recorded as a one-dimensional array H_n[] = {t8, U8, t9, U9, t10, U10}. By accessing the first element of the array H_n[], the value of the left minimum point on the abscissa can be determined, that is, H_n[0] = t8; by accessing the fourth element of the array H_n[], the value of the maximum value on the vertical axis can be determined, that is, H_n[3] = U9.

[0045] Step S206: Sort the peak descriptors corresponding to each maximum point in the forward order of the abscissa on the chromatogram within the preset range to obtain the first sequence.

[0046] Step S208: For each known component, based on the preset first window of the known component, the sequence fragment corresponding to the peak descriptor whose position in the first sequence is within the first window is used as the second sequence corresponding to the known component. Among them, the range of the first window is determined according to the retention time periods corresponding to each known component in the component detection results of other detected aerosol samples of the same specification as the aerosol sample to be detected.

[0047] In one example, the preset range on the chromatogram is Figure 1 the range shown, such as Figure 1 shown, which shows a total of 9 chromatographic peaks corresponding to four components that appear on the chromatogram. Among them, known component B has 1 main peak and 3 sub-peaks, and known component C has 1 main peak and 2 sub-peaks. At this time, the first sequence is expressed as {H_1[], H_2[], H_3[], H_4[], H_5[], H_6[], H_7[], H_8[], H_9[]}, where the sequence fragment corresponding to known component A is H_1[], the sequence fragment corresponding to known component B is H_2[], H_3[], H_4[], H_5[], the sequence fragment corresponding to abnormal component R (i.e., the unknown component) is H_6[], and the sequence fragment corresponding to known component C is H_7[], H_8[], H_9[].

[0048] Step S210: For each known component, the peak descriptor with the smallest sort order in the second sequence is used as the left peak descriptor, and the peak descriptor with the largest sort order in the second sequence is used as the right peak descriptor.

[0049] Taking known component B as an example, its second sequence is: H_2[], H_3[], H_4[], H_5[], the peak descriptor with the smallest sort order is H_2[], and the peak descriptor with the largest sort order is H_5[]. Therefore, H_2[] is the left peak descriptor and H_5[] is the right peak descriptor.

[0050] Step S212: For each known component, based on the left minimum point of the left peak descriptor, determine the starting point of the retention time period of the known component, and based on the right minimum point of the right peak descriptor, determine the ending point of the retention time period of the known component, so as to obtain the retention time periods corresponding to each known component.

[0051] In one example, the starting point of the retention time period corresponding to known component B is determined by accessing H_2[1], and the ending point of the retention time period corresponding to known component B is determined by accessing H_5[1]. The range between the ending point and the starting point is the retention time period corresponding to known component B.

[0052] Step S214: Based on the retention time periods corresponding to each known component, detect the amount of the known component.

[0053] After determining the retention time period corresponding to the known component B, the integration can be performed to detect the amount of the known component B.

[0054] In summary, compared with the method of segmenting the retention time periods corresponding to each component in Patent Document 1, when segmenting the chromatographic peaks corresponding to each known component in this embodiment, the length of the segmentation window is adjusted in increments of the sub-peak, which is more reasonable and can improve the accuracy of the segmentation results of the chromatographic peaks of each component. And because the segmentation window is adjusted in increments of the sub-peak, the number of adjustments is only the number of sub-peaks, reducing the number of times of adjusting the segmentation window and improving the data processing efficiency of determining the length of the segmentation window. In addition, because the length of the segmentation window is adjusted in increments of the sub-peak, there is no need to first perform the process of merging the main peak and the sub-peak, further improving the data processing efficiency.

[0055] As Figure 6 shown, in one embodiment, the method further includes a process step of verifying the component detection results of the atomized sample of the same specification that has been detected:

[0056] Step S302, sort the first windows corresponding to each known component in the order of the appearance time of the chromatographic peaks of each known component in the chromatogram to obtain a third sequence.

[0057] In an example, as Figure 1 shown, the range of the first window corresponding to the known component A is from t1 to t2, named CK_A, the range of the first window corresponding to the known component B is from t3 to t4, named CK_B, and the range of the first window corresponding to the known component C is from t5 to t6, named CK_C. Therefore, in the order of the appearance time of the chromatographic peaks of each known component in the chromatogram, the third sequence obtained after sorting the three first windows is: CK_A, CK_B, CK_C.

[0058] Step S304, according to the preset grouping strategy, divide the third sequence into multiple sequence segments to group each first window, and merge multiple adjacent first windows assigned to the same group to obtain a second window. The fourth sequence is composed of the unmerged first windows after grouping and the merged second windows.

[0059] Specifically, the third sequence is divided into two sequence segments. The first segment includes CK_A, and the second segment includes: CK_B, CK_C. Merge the two adjacent first windows CK_B and CK_C assigned to the same group to obtain the second window as Figure 2 shown. At this time, the second window is named CK_D. At this time, the unmerged first window CK_A and the merged second window CK_D together constitute the fourth sequence.

[0060] In one example, multiple adjacent first windows that are assigned to the same group are merged to obtain a second window, specifically including:

[0061] Step S3041: Use the smallest left boundary among the multiple first windows that are assigned to the same group as the left boundary of the second window.

[0062] Specifically, the second window CK_D is obtained by merging two adjacent first windows CK_B and CK_C that are assigned to the same group. The left boundary of the first window CK_B in this group is t3, and the left boundary of the first window CK_C is t5. As Figure 2 shown, t3 < t5. Therefore, use the left boundary of the first window CK_B as the left boundary of the second window CK_D, that is, the left boundary of the second window CK_D is t3.

[0063] Step S3042: Use the largest right boundary among the multiple first windows that are assigned to the same group as the right boundary of the second window.

[0064] Specifically, the second window CK_D is obtained by merging two adjacent first windows CK_B and CK_C that are assigned to the same group. The right boundary of the first window CK_B in this group is t4, and the left boundary of the first window CK_C is t6. As Figure 2 shown, t6 > t4. Therefore, use the right boundary of the first window CK_C as the right boundary of the second window CK_D, that is, the right boundary of the second window CK_D is t6.

[0065] Step S3043: Determine the range of the second window based on the left boundary and the right boundary of the second window.

[0066] Therefore, as Figure 2 shown, the range of the second window is from t3 to t6.

[0067] Step S306: Determine the retention period corresponding to a single known component based on the first window in the fourth sequence, detect the amount of the known component based on the retention period corresponding to the single known component, and verify the component detection result of the same known component in the same atomized sample that has been detected using the detected amount of the known component.

[0068] As Figure 2 shown, for the known component A, which is in a single-peak state, the corresponding retention period can be determined using the methods corresponding to steps S202 to S214 and the first window CK_A, or other methods in the prior art can also be used to determine the corresponding retention period, which will not be elaborated here.

[0069] Step S308: Determine the retention period corresponding to the corresponding known component group based on the second window in the fourth sequence, detect the amount of substance of the known component group based on the retention period corresponding to the known component group, and verify the sum of the amounts of substances of the known components that make up the known component group in the same atomized sample that has been detected using the detected amount of substance of the known component group. Among them, the known component group is an overall composed of multiple known components, and the second window corresponding to the known component group is obtained by merging the first windows corresponding to the multiple known components that make up the known component group.

[0070] In an example, the known component group is composed of known component B and known component C. This known component group corresponds to the second window CK_D. When considering the known component group as a single component, use the methods corresponding to steps S202 to S214 and the first window CK_A to determine its corresponding retention period. By integrating the retention period, obtain the area corresponding to the known component group, and then obtain the amount of substance corresponding to the known component group. When detecting a batch of atomized samples of the same specification for the first time, obtain the amounts of substances corresponding to known component B and known component C in this batch of samples. Specifically, assume there are 50 samples in the first batch of detections. Then, use the average value of the amounts of substance of known component B in the 50 samples as the amount of substance of known component B during the first detection. Similarly, obtain the amount of substance of known component C during the first detection. Compare and verify the sum of the amounts of substance of known component B and known component C during the first detection with the known component group, and determine whether the difference between the two exceeds the allowable deviation range. It should be noted that this allowable deviation range is the allowable error range of the amount of substance of this component during the production process. For example, the standard amount of substance of known component B during production is 50 units, and the allowable error range is 3%. The standard amount of substance of known component C during production is 100 units, and the allowable error range is 3%. Then, an error within 3% between the amount of substance calculated by integrating the known component group and the sum of known component B and C is considered normal.

[0071] In an example, determining the retention period corresponding to the corresponding known component group based on the second window in the fourth sequence specifically includes:

[0072] Step S3081: For each known component group, based on the second window corresponding to the known component group, use the sequence segment corresponding to the peak descriptor whose position in the first sequence is within the second window as the fifth sequence corresponding to the known component group.

[0073] Step S3082: Use the peak descriptor with the smallest sorting in the fifth sequence as the left peak descriptor, and use the peak descriptor with the largest sorting in the fifth sequence as the right peak descriptor.

[0074] Step S3083: For each known component group, determine the start point of the retention period of the known component group based on the left minimum point of the left peak descriptor and determine the end point of the retention period of the known component group based on the right minimum point of the right peak descriptor, so as to obtain the retention period corresponding to each known component group.

[0075] As Figure 2 shown, there are a total of 9 chromatographic peaks corresponding to four components appearing on the chromatogram. Among them, known component B has 1 main peak and 3 sub-peaks, abnormal component R has 1 chromatographic peak, and known component C has 1 main peak and 2 sub-peaks. At this time, the first sequence is expressed as {H_1[], H_2[], H_3[], H_4[], H_5[], H_6[], H_7[], H_8[], H_9[]}, where the sequence segment corresponding to known component A is H_1[], the sequence segments corresponding to known component B are H_2[], H_3[], H_4[], H_5[], the sequence segment corresponding to abnormal component R (i.e., the unknown component) is H_6[], and the sequence segments corresponding to known component C are H_7[], H_8[], H_9[]}. Then the fifth sequence is {H_2[], H_3[], H_4[], H_5[], H_6[], H_7[], H_8[], H_9[]}, where the peak descriptor with the smallest sorting is H_2[], and the peak descriptor with the largest sorting is H_9[]. Therefore, H_2[] is the left peak descriptor and H_9[] is the right peak descriptor. Determine the start point of the retention period corresponding to known component group w by accessing H_2[1], and determine the end point of the retention period corresponding to known component group w by accessing H_9[1]. The range between the end point and the start point is the retention period corresponding to known component group w.

[0076] In this embodiment, when verifying the detection result under the condition of knowing the amounts of substances of each known component in the first detection, since two adjacent components are combined into a known component group to jointly calculate the amount of substance by integration, when two adjacent known components are in a multi-peak state, the sub-peaks located between the two adjacent known components do not need to be merged first as in Patent Document 1, which improves the data processing efficiency.

[0077] In one embodiment, the method further includes:

[0078] Step S309: Subtract the sum of the amounts of substances of the known components that make up the known component group in the same atomized sample that has been detected from the amount of substance of the known component group in the verification results of each atomized sample, to obtain a plurality of verification differences corresponding to the known component group.

[0079] Step S310: Calculate the average value of the plurality of verification differences corresponding to the known component group.

[0080] Step S311, if the average value is greater than a preset difference threshold, it is prompted that there is an abnormal component.

[0081] In one embodiment, as Figure 2 shown, there is an abnormal component R between the known component B and the known component C. The abnormal component R is located within the second window. If the standard mass of the known component B in production is 50 units and the allowable error range is 3%, and the standard mass of the known component C in production is 100 units and the allowable error range is 3%, then the error between the mass calculated after integrating the known component group W and the sum of the known components B and C within 3% is normal. In this embodiment, due to the existence of the abnormal component R, the error between the mass calculated after integrating multiple known component groups W statistically obtained and the sum of the known components B and C is about 8%, which is greater than the preset difference threshold of 6%, then it is prompted that there is an abnormal component.

[0082] It can be understood that if the number of aerosol samples is 10, there are 10 verification results of the aerosol samples, and the masses of 10 known component groups W are obtained. When these 10 aerosol samples are detected for the first time, the masses of the known components B and C corresponding to the known component group W are obtained. Therefore, subtract the sum of the masses of the known components that make up the known component group in the same aerosol sample that has been detected (that is, the masses of the known components B and C corresponding to the known component group) from the mass of the known component group W in the verification result of each aerosol sample to obtain the verification difference corresponding to the known component group W in this aerosol sample. By analogy, 10 verification differences corresponding to the known component group W can be obtained, and the average value of the 10 verification differences is calculated. If the average value is greater than the preset difference threshold, it is prompted that there is an abnormal component. In this embodiment, using the average value obtained from multiple aerosol samples to judge the abnormal component can eliminate the interference caused by equipment problems during the detection of a single sample, and the detection result is more accurate.

[0083] As Figure 7 shown, in one embodiment, after it is prompted that there is an abnormal component, the method further includes:

[0084] Step S400, if the average value of the multiple verification differences corresponding to the known component group is greater than the preset difference threshold, then take this known component group as the target known component group.

[0085] Step S401, take the second window corresponding to the target known component group as the target second window.

[0086] Step S402, take each known component in the third sequence that is located within the target second window as the target known component.

[0087] Step S403, take the first window corresponding to each target known component as the target first window.

[0088] In this example, according to the third sequence, it can be seen that the target second window CK_D is obtained by merging two adjacent target first windows CK_B and CK_C that are assigned to the same group. Therefore, the target first windows include: the first windows CK_B and CK_C.

[0089] Step S404, determine the information of the abnormal component based on the peak descriptors between the target first windows with adjacent positions in the first sequence.

[0090] As Figure 1 shown, after obtaining the peak descriptors between the target first windows with adjacent positions in the first sequence, the retention time corresponding to the abnormal component R and the area of the chromatographic peak can be determined, so as to compare with the retention times corresponding to the components in the database, determine the component corresponding to the abnormal component R, and determine the amount of substance of the abnormal component R through the area.

[0091] In the process of verifying the detection results of known components in this embodiment, the existence of abnormal components can also be discovered, and the positions of the abnormal components can be quickly locked, improving the data processing efficiency.

[0092] In one embodiment, as Figure 8 shown, the obtaining of each maximum point and each minimum point within a preset range on the chromatogram of the to-be-detected aerosol sample specifically includes:

[0093] Step S502, starting from the horizontal axis starting coordinate of the preset range on the chromatogram, slide the acquisition window in the positive direction of the horizontal axis with a step size of one less than the width of the acquisition window.

[0094] In one example, the width of the acquisition window can be set to an initial preset value, such as 20 or 30. In one example, if the preset range on the chromatogram is the first window corresponding to a certain known component, the width of the acquisition window can be adjusted according to the characteristics of the chromatographic peak. For example, if the length of the first window is L and the number of peaks (including the main peak and sub-peaks) when the component is in a multi-peak state is N, then the width of the acquisition window is determined by the formula (L - (L % N)) / N. For example, for the known component B, the length L of the first window is 103 and the number of peaks N when in a multi-peak state is 4. Substituting into the formula, the width of the acquisition window is (103 - (103 % 4)) / 4 = 25. Through this formula, the possibility of having two maximum or minimum values in one acquisition window can be reduced as much as possible, avoiding missing the acquisition of a certain maximum or minimum value. It should be noted that the number of peaks of the known component B when in a multi-peak state can be obtained by manual observation of samples of the same specification in advance. When determining the first windows corresponding to each known component, the average value of the number of peaks corresponding to each known component when in a multi-peak state is used as N in the formula.

[0095] Step S503, when the acquisition window is at the starting position, perform the first extreme value acquisition operation, and record the monotonic characteristics of each part of the chromatogram picture segment within the current acquisition window in the temporary storage area.

[0096] The first extreme value acquisition operation includes:

[0097] Search for the maximum point and the minimum point within the current acquisition window;

[0098] As Figure 3 shown in a of, if the maximum point is located at the left boundary of the current acquisition window and the minimum point is located at a non-boundary position of the current acquisition window, add the minimum point to the extreme value point sequence and end the current first extreme value acquisition operation.

[0099] As Figure 3 shown in b of, if the maximum point is located at the left boundary of the current acquisition window and the minimum point is located at the right boundary of the current acquisition window, end the current first extreme value acquisition operation.

[0100] As Figure 3 shown in c of, if the maximum point is located at a non-boundary position of the current acquisition window and the minimum point is located at the left boundary of the current acquisition window, add the maximum point to the extreme value point sequence and end the current first extreme value acquisition operation.

[0101] As Figure 3 shown in d of, if both the maximum point and the minimum point are located at non-boundary positions of the current acquisition window and the position of the minimum point on the horizontal axis is less than that of the maximum point, add the minimum point and the maximum point to the extreme value point sequence in sequence and end the current first extreme value acquisition operation.

[0102] As Figure 3 shown in e of, if both the maximum point and the minimum point are located at non-boundary positions of the current acquisition window and the position of the maximum point on the horizontal axis is less than that of the minimum point, add the maximum point and the minimum point to the extreme value point sequence in sequence and end the current first extreme value acquisition operation.

[0103] As Figure 3 shown in f of, if the maximum point is located at a non-boundary position of the current acquisition window and the minimum point is located at the right boundary of the current acquisition window, add the maximum point to the extreme value point sequence and end the current first extreme value acquisition operation.

[0104] As Figure 3 shown in g of, if the maximum point is located at the right boundary of the current acquisition window and the minimum point is located at the left boundary of the current acquisition window, end the current first extreme value acquisition operation.

[0105] As Figure 3As shown in h, if the maximum point is located at the right boundary of the current acquisition window and the minimum point is located at a non-boundary position within the current acquisition window, add the minimum point to the extreme point sequence and end the current first extreme value acquisition operation.

[0106] Step S504: Each time the acquisition window slides one step to reach a new position, perform a second extreme value acquisition operation until the acquisition window covers the horizontal axis end coordinate of the preset range on the chromatogram. Among them, each time the acquisition window slides one step to reach a new position, record the monotonic characteristics of each part of the chromatogram segment within the current acquisition window in the temporary storage area; the monotonic characteristics of each part of the chromatogram segment along the horizontal axis are represented by an indication code, where 1 in the indication code represents the monotonically increasing part of the chromatogram segment, and 0 in the indication code represents the monotonically decreasing part of the chromatogram segment.

[0107] It can be understood that when the acquisition window slides to the horizontal axis end coordinate of the preset range on the chromatogram, the acquisition window may cover data points outside the preset range. At this time, use the horizontal axis end of the preset range on the chromatogram as the right boundary of the acquisition window, re-determine the range of the acquisition window, and then perform the second extreme value acquisition operation within this range.

[0108] It can be understood that the starting position of the acquisition window is regarded as the position when the acquisition window slides to the 0th step, the new position reached after sliding 1 step from the starting position is regarded as the position when the acquisition window slides to the 1st step, and the new position reached after sliding 2 steps from the starting position is regarded as the position when the acquisition window slides to the 2nd step, and so on for the positions corresponding to other steps, which will not be elaborated here.

[0109] The second extreme value acquisition operation includes:

[0110] Step A: Search for the maximum point and minimum point within the current acquisition window; among them, record the position where the acquisition window currently slides to the i-th step; i is a positive integer variable.

[0111] Step B: Based on the position relationship between the maximum point and the minimum point within the current acquisition window, obtain each maximum point and each minimum point on the chromatogram of the to-be-detected aerosol sample.

[0112] Specifically, as Figure 3 shown, the obtaining of each maximum point and each minimum point on the chromatogram of the to-be-detected aerosol sample based on the position relationship between the maximum point and the minimum point within the current acquisition window specifically includes:

[0113] As Figure 3As shown in a of [Figure ID], if the maximum point is located at the left boundary of the current acquisition window and the minimum point is located at a non-boundary position of the current acquisition window, compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add this minimum point to the extreme point sequence and end the current second extreme value acquisition operation; if not, add this minimum point and the maximum point to the extreme point sequence in sequence and end the current second extreme value acquisition operation.

[0114] As Figure 3 As shown in b of [Figure ID], if the maximum point is located at the left boundary of the current acquisition window and the minimum point is located at the right boundary of the current acquisition window, compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, end the current second extreme value acquisition operation; if not, add this maximum point to the extreme point sequence and end the current second extreme value acquisition operation.

[0115] As Figure 3 As shown in c of [Figure ID], if the maximum point is located at a non-boundary position of the current acquisition window and the minimum point is located at the left boundary of the current acquisition window, compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add this maximum point to the extreme point sequence and end the current second extreme value acquisition operation; if not, add this minimum point and the maximum point to the extreme point sequence in sequence and end the current second extreme value acquisition operation.

[0116] As Figure 3 As shown in d of [Figure ID], if both the maximum point and the minimum point are located at non-boundary positions of the current acquisition window and the position of the minimum point on the horizontal axis is less than that of the maximum point, compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add this minimum point and the maximum point to the extreme point sequence in sequence and end the current second extreme value acquisition operation; if not, add the points located at the left boundary of the current acquisition window, this minimum point, and the maximum point to the extreme point sequence in sequence and end the current second extreme value acquisition operation.

[0117] As Figure 3As shown in e, if both the maximum point and the minimum point are located at non-boundary positions within the current acquisition window and the position of the maximum point on the horizontal axis is less than that of the minimum point, then compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add the maximum point and the minimum point to the extreme point sequence in sequence, and end the current second extreme value acquisition operation; if not, add the points located on the left boundary of the current acquisition window, the maximum point, and the minimum point to the extreme point sequence in sequence, and end the current second extreme value acquisition operation.

[0118] As Figure 3 As shown in f, if the maximum point is located at a non-boundary position within the current acquisition window and the minimum point is located at the right boundary of the current acquisition window, then compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add the maximum point to the extreme point sequence and end the current second extreme value acquisition operation; if not, add the points located on the left boundary of the current acquisition window and the maximum point to the extreme point sequence in sequence, and end the current second extreme value acquisition operation.

[0119] As Figure 3 As shown in g, if the maximum point is located at the right boundary of the current acquisition window and the minimum point is located at the left boundary of the current acquisition window, then compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, end the current second extreme value acquisition operation; if not, add the minimum point to the extreme point sequence and end the current second extreme value acquisition operation.

[0120] As Figure 3 As shown in h, if the maximum point is located at the right boundary of the current acquisition window and the minimum point is located at a non-boundary position within the current acquisition window, then compare whether the last digit of the indication code recorded when the acquisition window slides to the (i - 1)-th step is the same as the first digit of the indication code recorded in the current acquisition window. If so, add the minimum point to the extreme point sequence and end the current second extreme value acquisition operation; if not, add the points located on the left boundary of the current acquisition window and the minimum point to the extreme point sequence in sequence, and end the current second extreme value acquisition operation.

[0121] In this embodiment, based on the position relationship between the maximum point and the minimum point within the current acquisition window, the maximum value and the minimum value within the acquisition window can be quickly determined. Compared with the prior art solution of using derivative calculation, this embodiment provides a simpler and more efficient extreme value search algorithm based on the characteristics of the chromatogram.

[0122] In one embodiment, as Figure 9As shown, the indication code corresponding to the current acquisition window is determined based on the following steps:

[0123] Step S11: Add the points on the left boundary, right boundary, minimum value point, and maximum value point of the current acquisition window to the indication set as elements in the indication set, and remove duplicate elements in the indication set to eliminate duplicate elements in the indication set;

[0124] Step S12: Sort the elements in the de-duplicated indication set in ascending order according to their positions on the horizontal axis, and let j = 2; where j is a variable;

[0125] Step S13: Subtract the value of the j-th element in the sorted indication set on the vertical axis from the value of the (j - 1)-th element on the vertical axis to obtain the coding result corresponding to the (j - 1)-th bit in the indication code; where if the value of the j-th element in the sorted indication set on the vertical axis minus the value of the (j - 1)-th element on the vertical axis is greater than zero, the coding result is 1, and if the value of the j-th element in the sorted indication set on the vertical axis minus the value of the (j - 1)-th element on the vertical axis is less than zero, the output result is 0;

[0126] Step S14: Determine whether j is equal to M. If so, jump to step S16 to end the current process; if not, adjust to step S15, let j + 1, and then jump to step S13. Where M is the number of elements in the de-duplicated indication set.

[0127] Taking Figure 3 a in

[0128] as an example, after sorting the elements in the de-duplicated indication set in ascending order according to their positions on the horizontal axis, they are successively: the point on the left boundary of the current acquisition window, the minimum value point, and the point on the right boundary, a total of three data points. Then, subtracting the value of the second point on the vertical axis from the value of the first point results in a value less than 0, and subtracting the value of the third point on the vertical axis from the value of the second point results in a value greater than 0. Therefore, the obtained indication code is 01.

[0129] In one embodiment, the method further includes:

[0130] If the maximum value point added to the extreme value point sequence when the acquisition window slides to the i-th step is adjacent to the maximum value point added to the extreme value point sequence when the acquisition window slides to the (i - 1)-th step, or the minimum value point added to the extreme value point sequence when the acquisition window slides to the i-th step is adjacent to the minimum value point added to the extreme value point sequence when the acquisition window slides to the (i - 1)-th step, then trigger a leak-check operation.

[0131] The leak detection operation includes:

[0132] If the maximum point added to the extreme point sequence when the acquisition window slides to the i-th step is adjacent to the maximum point added to the extreme point sequence when the acquisition window slides to the (i - 1)-th step, then the range between the two adjacent maximum points is used as the range of the leak detection window, and the minimum point is searched for in the leak detection window and inserted between the two adjacent maximum points in the extreme point sequence;

[0133] If the minimum point added to the extreme point sequence when the acquisition window slides to the i-th step is adjacent to the minimum point added to the extreme point sequence when the acquisition window slides to the (i - 1)-th step, then the range between the two adjacent minimum points is used as the range of the leak detection window, and the maximum point is searched for in the leak detection window and inserted between the two adjacent minimum points in the extreme point sequence.

[0134] It should be noted that when an acquisition window contains two maximum or minimum values at the same time, it may cause one of the extreme values to be missed, and by using the solution of this embodiment, the situation where the extreme value is missed can be detected, improving the accuracy of extreme value acquisition.

[0135] Such as Figure 4 As shown, in one embodiment, the searching for the maximum point and the minimum point within the current acquisition window specifically includes:

[0136] The first comparison process:

[0137] Step S21, obtain the first two points within the current acquisition window range in ascending order on the horizontal axis, compare them, store the point with the smaller vertical axis value in the first minimum value area, and store the point with the larger vertical axis value in the first maximum value area; let k = 3; where k is a variable; and if the vertical axis value of the second point is greater than the vertical axis value of the first point, then point the value address of the first reference point to the first maximum value area and jump to step S22; if the vertical axis value of the second point is less than the vertical axis value of the first point, then point the value address of the first reference point to the first minimum value area and jump to step S24;

[0138] Step S22, obtain the k-th point within the current acquisition window range in ascending order on the horizontal axis as the first target point, and compare the vertical axis value p of the first target point with the vertical axis value q of the first reference point. If q ≥ p, then update the first maximum value area with the first target point and jump to step S23; if q < p, then end the first process and jump to the second comparison process;

[0139] Step S23: Determine whether k is greater than the number of data points within the current acquisition window. If not, increment k by 1 and jump to step S22; if so, end the first comparison process and jump to the final comparison process;

[0140] Step S24: Obtain the k-th point within the current acquisition window range in ascending order along the horizontal axis as the first target point, and compare the value p of the vertical axis of the first target point with the value q of the vertical axis of the first reference point. If q ≤ p, update the first minimum value area using the first target point and jump to step S25; if q < p, end the first comparison process and jump to the second comparison process;

[0141] Step S25: Determine whether k is greater than the number of data points within the current acquisition window. If not, increment k by 1 and jump to step S24; if so, end the first process and jump to the final comparison process;

[0142] Second comparison process:

[0143] Step S31: Obtain the last two points within the current acquisition window range in ascending order along the horizontal axis for comparison, store the point with the smaller vertical axis value in the first minimum value area, and store the point with the larger vertical axis value in the first maximum value area; let m = 3, where m is a variable; and if the value of the vertical axis of the second-to-last point is greater than the value of the vertical axis of the last point, point the value address of the second reference point to the first maximum value area and jump to step S32; if the value of the vertical axis of the second-to-last point is less than the value of the vertical axis of the last point, point the value address of the second reference point to the first minimum value area and jump to step S34;

[0144] Step S32: Obtain the m-th point from the end within the current acquisition window range in ascending order along the horizontal axis as the second target point, and compare the value p of the vertical axis of the second target point with the value q of the vertical axis of the second reference point. If q ≥ p, update the first maximum value area using the second target point and jump to step S33; if q < p, end the second process and jump to the final comparison process;

[0145] Step S33: Increment m by 1 and jump to step S32;

[0146] Step S34: Obtain the m-th point from the end within the current acquisition window range in ascending order along the horizontal axis as the second target point, and compare the value p of the vertical axis of the second target point with the value q of the vertical axis of the second reference point. If q ≤ p, update the first minimum value area using the second target point and jump to step S35; if q > p, end the second process and jump to the final comparison process;

[0147] Step S35: Increment m by 1 and jump to step S34;

[0148] Final comparison process:

[0149] If both the second smallest value area and the second largest value area are the initial default values, the point in the first smallest value area is used as the minimum value point within the current acquisition window range, and the point in the first largest value area is used as the maximum value point within the current acquisition window range;

[0150] If neither the second smallest value area nor the second largest value area is the initial default value, compare the values of the points in the first smallest value area and the second smallest value area on the vertical axis, and use the smaller point as the minimum value point within the current acquisition window range. If the values of the points in the first smallest value area and the second smallest value area on the vertical axis are equal, use the point in the first smallest value area as the minimum value point within the current acquisition window range; compare the values of the points in the first largest value area and the second largest value area on the vertical axis, and use the larger point as the maximum value point within the current acquisition window range. If the values of the points in the first largest value area and the second largest value area on the vertical axis are equal, use the point in the first largest value area as the maximum value point within the current acquisition window range.

[0151] The maximum and minimum search algorithm provided in this embodiment utilizes the characteristics of the data points in the chromatogram within the acquisition window. Compared with the maximum and minimum search algorithms in the prior art that require 2n or 1.5n comparisons to search for the maximum and minimum values, where n is the number of elements in the dataset for which the maximum and minimum values need to be searched, this embodiment only requires n or n + 2 comparisons to obtain the maximum and minimum values among n data. Its working principle is as follows, as Figure 4 shown, for the n data points in the chromatogram segment that first increases and then decreases within the acquisition window, the data points in the increasing part are only compared with the data points in the first largest value area, and the data points in the decreasing part are only compared with the data points in the second largest value area. Therefore, the total number of comparisons for size is n times. Finally, if both the second smallest value area and the second largest value area are the initial default values, there is no need to compare anymore, so the total number of comparisons is n times. If neither the second smallest value area nor the second largest value area is the initial default value, then compare the first smallest value area with the second smallest value area, and compare the first largest value area with the second largest value area, for a total of n + 2 comparisons.

[0152] In one embodiment, there is also provided an electronic aerosol component detection device, including: a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it performs the steps of the above-mentioned electronic aerosol component detection method. The steps of the electronic aerosol component detection method here can be the steps in the above-mentioned various embodiments.

[0153] In one embodiment, a computer-readable storage medium is further provided. The computer-readable storage medium stores computer-executable instructions for causing a computer to execute the steps of the above-described electronic aerosol composition detection method. The steps of the electronic aerosol composition detection method here may be the steps in the electronic aerosol composition detection methods of the above various embodiments.

[0154] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it may include the processes of the embodiments of the above various methods. Among them, any reference to a memory, storage, database, or other medium used in the various embodiments provided in the present application may include non-volatile and / or volatile memories. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or an external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0155] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered to be within the scope described in this specification.

Claims

1. A method for detecting components of electronic atomized matter, characterized in that: The method comprises: Obtaining each maximum value point and each minimum value point within a preset range on the chromatogram of the atomized sample to be detected; The coordinates of each maximum point, the left minimum point adjacent to each maximum point, and the right minimum point within a preset range on the chromatogram are collectively used as a peak descriptor corresponding to a maximum point; wherein the abscissa of the maximum point in the peak descriptor is used as the position of the peak descriptor; According to the positive order of each maximum point on the horizontal axis within a preset range on the chromatogram, the peak descriptors corresponding to each maximum point are sorted to obtain a first sequence; For each known component, based on a first window preset for the known component, a sequence fragment corresponding to a peak descriptor in the first sequence that is located within the first window is used as a second sequence corresponding to the known component; wherein the range of the first window is determined according to the retention time period corresponding to each known component in the component detection results of other detected aerosol samples that belong to the same specification as the aerosol sample to be detected; The peak descriptor with the smallest ranking in the second sequence is used as the left peak descriptor, and the peak descriptor with the largest ranking in the second sequence is used as the right peak descriptor; For each known component, the starting point of the retention time period of the known component is determined based on the left minimum point of the left peak descriptor, and the end point of the retention time period of the known component is determined based on the right minimum point of the right peak descriptor, to obtain the retention time period corresponding to each known component; Based on the retention time period corresponding to each known component, detecting the amount of the known component; The method also includes the process step of verifying the composition test results of the aerosol samples of the same specification that have been tested: According to the time sequence of the appearance of the chromatographic peaks of the known components on the chromatogram, the first windows corresponding to the known components are sorted to obtain a third sequence; According to a preset grouping strategy, the third sequence is divided into a plurality of sequence segments to group the first windows, and a plurality of adjacent first windows grouped into the same group are merged to obtain second windows, and the first windows that are not merged after the grouping and the merged second windows together constitute a fourth sequence; Determining a retention time period corresponding to a corresponding single known component based on the first window in the fourth sequence, detecting the substance amount of the known component based on the retention time period corresponding to the single known component, and verifying the component detection result of the same known component in the same atomized sample that has been detected using the substance amount of the known component obtained by the detection; Based on the second window in the fourth sequence, a retention time period corresponding to the corresponding known component cluster is determined, and based on the retention time period corresponding to the known component cluster, the mass of the known component cluster is detected, and the mass of the known component cluster obtained by the detection is used to verify the sum of the masses of each known component constituting the known component cluster in the same atomized sample that has been detected; wherein the known component cluster is a whole composed of multiple known components, and the second window corresponding to the known component cluster is obtained by merging the first windows corresponding to the multiple known components constituting the known component cluster.

2. The method for detecting components of electronic atomized matter according to claim 1, characterized in that: The method further comprises: The mass of the known component cluster in the verification result of each aerosol sample is subtracted from the sum of the mass of each known component constituting the known component cluster in the same aerosol sample that has been tested, to obtain a plurality of verification differences corresponding to the known component cluster; Calculate the average of multiple verification differences corresponding to known component groups; If the average value is greater than the preset difference threshold, it indicates that there is an abnormal component.

3. The method for detecting components of electronic atomization according to claim 2, characterized in that: After indicating the presence of abnormal components, the method further includes: If the average of multiple verification differences corresponding to a known component cluster is greater than a preset difference threshold, the known component cluster is taken as a target known component cluster; The second window corresponding to the target known component cluster is used as the target second window; taking each known component in the third sequence located in the target second window as the target known component; The first window corresponding to each target known component is used as the target first window; The information of abnormal components is determined based on peak descriptors located between adjacent target first windows in the first sequence.

4. The method for detecting components of electronic atomized matter according to claim 1, characterized in that: Determining a retention time period corresponding to a corresponding known component group based on the second window in the fourth sequence specifically includes: For each known component cluster, based on the second window corresponding to the known component cluster, the sequence fragment corresponding to the peak descriptor in the first sequence whose position is within the second window is used as the fifth sequence corresponding to the known component cluster; The peak descriptor with the smallest ranking in the fifth sequence is used as the left peak descriptor, and the peak descriptor with the largest ranking in the fifth sequence is used as the right peak descriptor; For each known component cluster, the starting point of the retention time period of the known component cluster is determined based on the left minimum point of the left peak descriptor, and the end point of the retention time period of the known component cluster is determined based on the right minimum point of the right peak descriptor, to obtain the retention time period corresponding to each known component cluster.

5. The method for detecting components of electronic atomized matter according to claim 1, characterized in that: Merging multiple adjacent first windows that are grouped into the same group to obtain a second window includes: The smallest left border among the plurality of first windows grouped into the same group is used as the left border of the second window; The largest right border among the plurality of first windows grouped into the same group is used as the right border of the second window; Based on the left border of the second window and the right border of the second window, a range of the second window is determined.

6. The method for detecting components of electronic atomized matter according to claim 1, characterized in that: The step of obtaining each maximum point and each minimum point within a preset range on the chromatogram of the atomized sample to be detected specifically includes: Starting from the horizontal axis starting point coordinate of the preset range on the chromatogram, slide the acquisition window in the positive direction of the horizontal axis with the width of the acquisition window minus one as the step length; When the acquisition window is at the starting point, the first extreme value acquisition operation is performed, and the monotonic characteristics of each part of the chromatogram segment within the current acquisition window range are recorded in the temporary storage area; The acquisition window performs a second extreme value acquisition operation each time it slides one step to a new position, until the acquisition window covers the horizontal axis end point coordinate of a preset range on the chromatogram; wherein, the monotonic characteristics of each part of the chromatogram segment within the current acquisition window range are recorded in a temporary storage area each time the acquisition window slides one step to a new position; the monotonic characteristics of each part of the chromatogram segment along the horizontal axis are represented by an indicator code, wherein 1 in the indicator code represents a monotonically increasing part of the chromatogram segment, and 0 in the indicator code represents a monotonically decreasing part of the chromatogram segment; The second extreme value acquisition operation includes: Search for the maximum and minimum points within the current acquisition window; where the acquisition window is currently sliding to the position of the i-th step; Based on the positional relationship between the maximum value point and the minimum value point in the current acquisition window, each maximum value point and each minimum value point on the chromatogram of the atomized sample to be detected is obtained.

7. The method for detecting components of electronic atomized matter according to claim 6, characterized in that: The step of obtaining each maximum point and each minimum point on the chromatogram of the atomized sample to be detected based on the positional relationship between the maximum point and the minimum point in the current acquisition window specifically includes: If the maximum value point is located at the left boundary of the current acquisition window and the minimum value point is located at the non-boundary position of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, add the minimum value point to the extreme value point sequence and end the current extreme value acquisition operation; if not, add the minimum value point and the maximum value point to the extreme value point sequence in turn and end the current extreme value acquisition operation; If the maximum value point is located at the left boundary of the current acquisition window and the minimum value point is located at the right boundary of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, end the current extreme value acquisition operation; if not, add the maximum value point to the extreme value point sequence and end the current extreme value acquisition operation; If the maximum value point is located at the non-boundary position of the current acquisition window and the minimum value point is located at the left boundary of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, add the maximum value point to the extreme value point sequence and end the current extreme value acquisition operation; if not, add the minimum value point and the maximum value point to the extreme value point sequence in turn and end the current extreme value acquisition operation; If both the maximum point and the minimum point are located at the non-boundary position of the current acquisition window and the position of the minimum point on the horizontal axis is smaller than the maximum point, then compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, then add the minimum point and the maximum point to the extreme value point sequence in turn, and end the current extreme value acquisition operation; if not, then add the points located on the left boundary of the current acquisition window and the minimum point and the maximum point to the extreme value point sequence in turn, and end the current extreme value acquisition operation; If both the maximum point and the minimum point are located at the non-boundary position of the current acquisition window and the position of the maximum point on the horizontal axis is smaller than the minimum point, then compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, then add the maximum point and the minimum point to the extreme value point sequence in turn, and end the current extreme value acquisition operation; if not, then add the points located on the left boundary of the current acquisition window and the maximum point and the minimum point to the extreme value point sequence in turn, and end the current extreme value acquisition operation; If the maximum value point is located at a non-boundary position of the current acquisition window and the minimum value point is located at the right boundary of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, add the maximum value point to the extreme value point sequence and end the current extreme value acquisition operation; if not, add the point located on the left boundary of the current acquisition window and the maximum value point to the extreme value point sequence in sequence, and end the current extreme value acquisition operation; If the maximum value point is located at the right boundary of the current acquisition window and the minimum value point is located at the left boundary of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, end the current extreme value acquisition operation; if not, add the minimum value point to the extreme value point sequence and end the current extreme value acquisition operation; If the maximum value point is located at the right boundary of the current acquisition window and the minimum value point is located at the non-boundary position of the current acquisition window, compare whether the last bit of the indicator code recorded when the acquisition window slides to the i-1th step is the same as the first bit of the indicator code recorded in the current acquisition window. If so, add the minimum value point to the extreme value point sequence and end the current extreme value acquisition operation; if not, add the points on the left boundary of the current acquisition window and the minimum value point to the extreme value point sequence in turn, and end the current extreme value acquisition operation.

8. The method for detecting components of electronic atomized matter according to claim 6, characterized in that: Determine the indicator code corresponding to the current acquisition window based on the following steps: Step S11, adding the points on the left boundary, the points on the right boundary, the minimum point and the maximum point of the current acquisition window to the indication set as elements in the indication set, and deduplicating the indication set to remove duplicate elements in the indication set; Step S12, sorting the elements in the indication set after deduplication from small to large according to their positions on the horizontal axis, and setting j=2; where j is a variable; Step S13, subtract the value of the j-1th element on the vertical axis from the value of the j-th element in the sorted indicator set, to obtain the encoding result corresponding to the j-1th bit in the indicator code; wherein, if the value of the j-th element on the vertical axis in the sorted indicator set minus the value of the j-1th element on the vertical axis is greater than zero, the encoding result is 1; if the value of the j-th element on the vertical axis in the sorted indicator set minus the value of the j-1th element on the vertical axis is less than zero, the output result is 0; Step S14, determine whether j is equal to M. If so, end the current process; if not, set j+1 and jump to step S13; where M is the number of elements in the indicator set after deduplication.

9. An electronic atomization component detection device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the method for detecting components of electronic atomizers as described in any one of claims 1 to 8 are performed.

Images