Electronic atomized matter component detection method and equipment

Abstract

The invention discloses a method and equipment for detecting components of an electronic atomized substance, and belongs to the technical field of component detection by using gas chromatography. According to the method, when the chromatographic peak corresponding to each known component is segmented, the length of the segmentation window is adjusted by taking the secondary peak as the increment, so that the method is more reasonable, and the accuracy of the segmentation result of the chromatographic peak of each component can be improved. And the segmented window is adjusted by taking the secondary peak as the increment, and the adjustment frequency is only the number of the secondary peaks, so that the frequency of adjusting the segmented window is reduced, and the data processing efficiency of determining the length of the segmented window is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of utilizing gas chromatography to detect components, and in particular to a method and device for detecting components of electronic atomized matter. Background Art

[0002] In the process of detecting the components of the smoke oil atomizer using a gas chromatograph, due to the fluctuation of the chromatographic column temperature and the carrier gas flow rate, the same component corresponds to a main peak and multiple secondary peaks at the same time (hereinafter referred to as a multi-peak state). Therefore, as described in Patent Document 1 (Publication No. CN118501325A), when dividing the corresponding retention time period for each component in the chromatogram, it is necessary to first identify a main peak and multiple secondary peaks corresponding to the same component, and then merge the peaks, and determine the length range of the window according to the merging results; then select a window of random length within the length range and perform parameter calculation with the retention time period of each known component (the parameters include evaluation factors, initial attribution, overlapping characteristics, and final attribution, etc.), and continuously adjust the window length according to the calculation results, and finally obtain a window with the highest degree of conformity with the retention time period of the known component and segment the chromatogram with this window.

[0003] However, the solution in Patent Document 1 still has the following shortcomings: 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, and finally obtain a window that has the highest degree of consistency 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 in the data processing process used to determine the final window length.

[0004] 2. In the process of adjusting the window length, the increment of the constraint window length adjustment is a manually preset unit time interval. If the unit time interval is set unreasonably, the accuracy of the segmentation results of the chromatographic peaks of each component will also be reduced. Summary of the invention

[0005] The purpose of the present invention is to solve at least one of the technical problems existing in the prior art and to provide a method and device for detecting components of electronic atomizers, which can improve the efficiency of data processing and enhance the accuracy of the chromatographic peak segmentation results of each component.

[0006] In order to achieve the above-mentioned object of the invention, the present invention provides the following technical solutions: A method for detecting components of electronic atomization, the method comprising: 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, the amount of substance of the known component is detected.

[0007] On the other hand, the present invention also provides an electronic atomizer component detection device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the steps of the above-mentioned electronic atomizer component detection method when executing the program.

[0008] Beneficial effect: Compared with the method of segmenting the retention time period corresponding to each component in Patent Document 1, the electronic atomization 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 secondary peaks, 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 secondary peaks, the number of adjustments is only the number of secondary peaks, which reduces the number of times the segmentation window is adjusted, and improves the data processing efficiency of determining the length of the segmentation window.

[0009] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The present invention is further described below in conjunction with the accompanying drawings and embodiments; Figure 1is a schematic diagram of a first window range corresponding to some known components on a chromatogram in an embodiment; Figure 2 A schematic diagram of a second window range corresponding to some known component groups on a chromatogram in one embodiment; Figure 3 It is a schematic diagram of the state when the maximum value and the minimum value in the acquisition window are in different position relationships in one embodiment; Figure 4 A schematic diagram of the principle of searching for maximum and minimum points within the current acquisition window in one embodiment; Figure 5 A schematic diagram of a flow chart of a method for detecting components of electronic atomized matter in an embodiment; Figure 6 A schematic diagram of a process for verifying a test result in one embodiment; Figure 7 A schematic diagram of a process for determining abnormal component information in one embodiment; Figure 8 A schematic diagram of a process for obtaining a maximum value point and a minimum value point in one embodiment; Fig. 9 The figure is a flow chart of determining the indication code corresponding to the current acquisition window in one embodiment. DETAILED DESCRIPTION

[0011] This section will describe in detail the specific embodiments of the present invention. The preferred embodiments of the present invention are shown in the accompanying drawings. The purpose of the accompanying drawings is to supplement the description of the text part of the specification with graphics, so that people can intuitively and vividly understand each technical feature and the overall technical solution of the present invention, but it cannot be understood as a limitation on the scope of protection of the present invention.

[0012] It should be noted that when performing chromatographic separation, each component has a different retention time depending on its interaction 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 balance between the stationary phase and the mobile phase. Changes in the carrier gas flow rate will affect the migration time of the substance between the stationary phase and the mobile phase. Therefore, in the process of detecting the components of the cigarette oil atomizer using a gas chromatograph, due to fluctuations in the chromatographic column temperature and the carrier gas flow rate, the same component corresponds to a main peak and multiple secondary peaks at the same time (hereinafter referred to as a multi-peak state). Specifically, Figure 1 As shown in the figure, there are four chromatographic peaks corresponding to the known component B in the interval from t3 to t4 on the horizontal axis. The peak with the largest peak height is defined as the main peak, and the other three peaks with lower heights are defined as secondary peaks. Under ideal conditions, the chromatographic peak of a component is a single peak. Figure 1The chromatographic peak corresponding to the known component A shown in the interval from t1 to t2 on the horizontal axis is a single peak. However, due to the limited control accuracy of the chromatographic column temperature and the carrier gas flow rate, the corresponding chromatographic peak of a certain component may be a single peak during this test due to the ideal control effect of the chromatographic column temperature and the carrier gas flow rate. In the next inspection, although the same set of chromatographs and the same chromatograph control parameters are used, the chromatographic column temperature and the carrier gas flow rate may fluctuate due to the limited control accuracy of the chromatographic column temperature and the carrier gas flow rate, which may cause the chromatographic peak of the component to be multi-peaked. For known components, under the same chromatograph control parameters, their retention time can generally be determined by measuring the standard substance. For example, if the retention time of a chromatographic peak detected in the sample is consistent with the retention time of the chromatographic peak in the standard solution of nicotine, for example, within plus or minus 10 seconds of 32 minutes, it is considered that the component corresponding to the chromatographic peak is nicotine. The area is obtained by integrating the chromatographic peak, and then the area of ​​nicotine is determined based on the relationship between the area and the amount of substance that has been calibrated in advance. For example, the relationship between area and amount of substance can be determined by the area internal standard method, specifically, take the standard measured component, and add each to the quantitative internal standard substance specified by each monomer according to the known stage amount that increases or decreases in sequence, and modulate the standard solution. Take a certain amount of this standard solution and inject it into the gas chromatograph column respectively. According to the chromatogram on the gas chromatograph, take the ratio of the peak area and peak height of the standard measured component to the peak area and peak height of the internal standard substance as the vertical coordinate, take the ratio of the amount of the standard measured component to the mass of the internal standard substance, or the amount of the standard measured component as the horizontal coordinate, and make a standard curve. Then modulate the sample solution according to the method specified in the monomer. When modulating the sample solution, add the same amount of internal standard substance as when modulating the standard solution in advance. Then, according to the chromatogram obtained under the same conditions when making the standard curve, find the ratio of the peak area or peak height of the measured component to the peak product or peak height of the internal standard substance, and then find the content of the measured component according to the standard curve.

[0013] See Patent Document 2 (Publication No. CN111595992A). For a 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, and then the starting point and end point of the chromatographic peak (i.e., the retention time period) can be determined by the two left and right minimum points. The retention time period is integrated to obtain the amount of substance of the component. For a multi-peak state, the retention time period corresponding to the chromatographic peak of the multi-peak state of each known component can be determined by an expert through manual calibration of 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 divided into corresponding components, and the difference between the area of ​​the component under different division methods and the amount of substance of the component in the sample of the same specification in a single peak state is calculated, and the division method is selected when the difference in the amount of substance meets the allowable amount error of the component in production to determine the retention time period corresponding to the component in the multi-peak state. If the chromatographic peak of a known component B of a sample of this specification is in a single-peak state and the calculated amount of substance is 50 units when the chromatographic instrument detects it, the allowable amount error of the component in production is 3%, and in another test, the chromatographic peak of the known component B in the sample of the same specification forms a multi-peak state, and the two secondary peaks on the left side and the single secondary peak on the right side of the main peak of the known component B are manually divided into the chromatographic peak corresponding to the known component B, and the left minimum of the leftmost secondary peak is used as the starting point, and the right minimum of the right secondary peak is used as the end point to determine the retention time period of the known component B in the multi-peak state, and the amount of substance of the known component B obtained under this division method is calculated to be 52 units, and the error is within 3%, while other division methods do not meet the 3% error, so the retention time period corresponding to the component B can be determined, and multiple retention time periods corresponding to the known component B are obtained by multiple samples of the same specification, and the starting points corresponding to the multiple retention time periods corresponding to the known component B obtained by the manual method are averaged as the left boundary t3 of the first window corresponding to the known component B, and the end points are averaged to obtain the right boundary t4. Since the mass fluctuation of the components of the same batch of samples with the same specifications is low, the retention time period corresponding to each known component in the multi-peak state is manually determined as a reference, so as to determine the mass of the corresponding substance of each known component in the multi-peak state. For the scheme of detecting each known component in the multi-peak state, refer to Patent Document 1.

[0014] As described in Patent Document 1, when dividing the corresponding retention time period for each component in the chromatogram, it is necessary to first identify a main peak and multiple secondary peaks corresponding to the same component and then merge the peaks, and determine the length range of the window based on the merging results; then select a window of random length within the length range and perform parameter calculation with the retention time period of each known component (the parameters include evaluation factors, initial attribution, overlapping characteristics, and final attribution, etc.), and continuously adjust the window length based on the calculation results, and finally obtain a window with the highest degree of conformity with the retention time period of the known component and use the window to segment the chromatogram.

[0015] However, the solution in Patent Document 1 still has the following shortcomings: 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, and finally obtain a window that has the highest degree of consistency 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 in the data processing process used to determine the final window length.

[0016] 2. In the process of adjusting the window length, the increment of the constraint window length adjustment is a manually preset unit time interval. If the unit time interval is set unreasonably, the accuracy of the segmentation results of the chromatographic peaks of each component will also be reduced.

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

[0018] Based on this, the present invention provides a method for detecting components of electronic atomizers, which can improve the efficiency of data processing and enhance the accuracy of the chromatographic peak segmentation results of each component.

[0019] A method for detecting components of electronic atomized matter provided by the present invention will be described in detail below in conjunction with the accompanying drawings and embodiments.

[0020] like Figure 5 As shown, in one embodiment, a method for detecting components of electronic atomization is provided, the method comprising: Step S202, obtaining each maximum point and each minimum point within a preset range on the chromatogram of the atomized sample to be detected.

[0021] Specifically, each chromatogram is a voltage-time graph generated by the different voltage signals exhibited in the detector when a portion of electronic cigarette oil (i.e., atomizer 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, the chromatogram corresponding to each atomizer sample has a total of 8,000 data points, arranged in sequence starting from 0, and every 1,000 data points represent one minute. The chromatogram corresponding to each atomizer sample is the data for all 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, such as Figure 1 As shown, the preset range on the chromatogram can be three intervals, namely, t1 to t2, t3 to t4, and t5 to t6.

[0022] Step S204, the coordinates of each maximum point, the left minimum point adjacent to each maximum point, and the right minimum point within the 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.

[0023] 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 the various information that constitute a peak. The maximum point records the peak height, and the distance between the left minimum point and the right minimum point on the horizontal axis is the peak width. The horizontal coordinate value of the maximum value in the peak descriptor is used as the position of the peak descriptor. For example, if the horizontal coordinate of the maximum value is t7, then the position of the peak descriptor is t7.

[0024] 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 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). The peak descriptor is a six-tuple, which 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 on the horizontal coordinate can be determined, that is, H_n[0]=t8; by accessing the fourth element of the array H_n[], the value of the maximum on the vertical coordinate can be determined, that is, H_n[3]=U9.

[0025] Step S206 , 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.

[0026] Step S208, for each known component, based on the first window preset for the known component, the sequence fragment corresponding to the peak descriptor in the first sequence that is located within the first window is used as the second sequence corresponding to the known component. The range of the first window is determined based on 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.

[0027] In one example, the preset range on the chromatogram is Figure 1 The range shown, such as Figure 1As shown, a total of 9 chromatographic peaks corresponding to the four components appearing on the chromatogram are shown, among which the known component B has 1 main peak and 3 secondary peaks, and the known component C has 1 main peak and 2 secondary peaks. At this time, the first sequence is represented by {H_1[], H_2[], H_3[], H_4[], H_5[], H_6[], H_7[], H_8[], H_9[]}, among which the sequence fragment corresponding to the known component A is H_1[], the sequence fragment corresponding to the known component B is H_2[], H_3[], H_4[], H_5[], the sequence fragment corresponding to the abnormal component R (i.e., the unknown component) is H_6[], and the sequence fragment corresponding to the known component C is H_7[], H_8[], H_9[].

[0028] Step S210 : using the peak descriptor with the smallest ranking in the second sequence as the left peak descriptor, and using the peak descriptor with the largest ranking in the second sequence as the right peak descriptor.

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

[0030] Step S212, for each known component, determine the starting point of the retention time period of the known component based on the left minimum point of the left peak descriptor and determine the end point of the retention time period of the known component based on the right minimum point of the right peak descriptor, to obtain the retention time period corresponding to each known component.

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

[0032] Step S214, detecting the amount of substance of each known component based on the retention time period corresponding to each known component.

[0033] After determining the retention time period corresponding to the known component B, the amount of the known component B can be detected by integration.

[0034] In summary, compared with the method of segmenting the retention time period corresponding to each component in Patent Document 1, when the chromatographic peaks corresponding to each known component are segmented in this embodiment, the length of the segmentation window is adjusted in increments of secondary peaks, 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 secondary peaks, the number of adjustments is only the number of secondary peaks, which reduces the number of times the segmentation window is adjusted and improves 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 secondary peaks, there is no need to first perform the process of merging the main peak and the secondary peak, which further improves the data processing efficiency.

[0035] like Figure 6 As shown, in one embodiment, the method further includes a process step of verifying the component test results of the aerosol sample of the same specification that has been tested: Step S302, sorting the first windows corresponding to the known components according to the time sequence of the chromatographic peaks of the known components appearing on the chromatogram to obtain a third sequence.

[0036] In one example, if Figure 1 As 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, according to the time sequence of the chromatographic peaks of each known component appearing on the chromatogram, the third sequence obtained after sorting the three first windows is: CK_A, CK_B, CK_C.

[0037] Step S304: 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. The first windows that are not merged after grouping and the merged second windows together constitute a fourth sequence.

[0038] 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. The two adjacent first windows CK_B and CK_C that are grouped into the same group are merged to obtain Figure 2 The second window is shown in FIG. 1 , and the second window is named CK_D. At this time, the unmerged first window CK_A and the merged second window CK_D together form a fourth sequence.

[0039] In one example, merging a plurality of adjacent first windows that are grouped into the same group to obtain a second window specifically includes: Step S3041: 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.

[0040] Specifically, the second window CK_D is obtained by merging two adjacent first windows CK_B and CK_C that are grouped into the same group. The left edge of the first window CK_B in the group is t3, and the left edge of the first window CK_C is t5. Figure 2 As shown, t3<t5, therefore, the left border of the first window CK_B is used as the left border of the second window CK_D, that is, the left border of the second window CK_D is t3.

[0041] Step S3042: 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.

[0042] Specifically, the second window CK_D is obtained by merging two adjacent first windows CK_B and CK_C that are grouped into the same group. The right edge of the first window CK_B in the group is t4, and the left edge of the first window CK_C is t6. Figure 2 As shown, t6>t4, therefore, the right border of the first window CK_C is used as the right border of the second window CK_D, that is, the right border of the second window CK_D is t6.

[0043] Step S3043: determine the range of the second window based on the left border of the second window and the right border of the second window.

[0044] Therefore, if Figure 2 As shown, the range of the second window is from t3 to t6.

[0045] Step S306, determining the retention time period corresponding to the corresponding single known component based on the first window in the fourth sequence, and detecting the substance amount of the known component based on the retention time period corresponding to the single known component, and using the substance amount of the known component obtained by the detection to verify the component detection result of the same known component in the same atomized sample that has been detected.

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

[0047] Step S308, based on the second window in the fourth sequence, the retention time period corresponding to the corresponding known component cluster is determined, and the substance amount of the known component cluster is detected based on the retention time period corresponding to the known component cluster, and the substance amount of the known component cluster obtained by the detection is used to verify the sum of the substance amounts 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.

[0048] In one example, the known component cluster is composed of known component B and known component C. The known component cluster corresponds to the second window CK_D. When the known component cluster is regarded as one component, the method corresponding to steps S202 to S214 and the first window CK_A are used to determine the corresponding retention time period. The area corresponding to the known component cluster is obtained by integrating the retention time period, and then the amount of substance corresponding to the known component cluster is obtained. When a batch of atomized samples of the same specification is tested for the first time, the corresponding amounts of known components B and known components C in the batch of samples are obtained. Specifically, assuming that there are 50 samples in the first batch of tests, the average amount of the known component B in the 50 samples is used as the amount of the known component B in the first test. Similarly, the amount of the known component C in the first test is obtained. The sum of the amounts of known components B and known components C in the first test is compared and verified with the known component cluster to determine whether the difference between the two exceeds the allowable deviation range. It should be noted that the allowable deviation range is the allowable mass error range of the component in the production process. For example, if the mass of the standard substance in the production of known component B is 50 units, the allowable error range is 3%; if the mass of the standard substance in the production of known component C is 100 units, the allowable error range is 3%. In this case, it is normal for the mass calculated after integration of the known component mass and the sum of the known components B and C to have an error within 3%.

[0049] In one example, determining the retention time period corresponding to the corresponding known component cluster based on the second window in the fourth sequence specifically includes: Step S3081: 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.

[0050] Step S3082: 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.

[0051] Step S3083, for each known component cluster, determine the starting point of the retention time period of the known component cluster based on the left minimum point of the left peak descriptor and determine the end point of the retention time period of the known component cluster based on the right minimum point of the right peak descriptor, to obtain the retention time period corresponding to each known component cluster.

[0052] like Figure 2 As shown, a total of 9 chromatographic peaks corresponding to the four components appearing on the chromatogram are shown, among which the known component B has 1 main peak and 3 secondary peaks, the abnormal component R has 1 chromatographic peak, and the known component C has 1 main peak and 2 secondary peaks. At this time, the first sequence is represented by {H_1[], H_2[], H_3[], H_4[], H_5[], H_6[], H_7[], H_8[], H_9[]}, among which the sequence fragment corresponding to the known component A is H_1[], the sequence fragment corresponding to the known component B is H_2[], H_3[], H_4[], H_5[], and the sequence fragment corresponding to the abnormal component R (i.e., the unknown component) is H_6[], The sequence fragments corresponding to the 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[]}, among which the peak descriptor with the smallest ranking is H_2[], and the peak descriptor with the largest ranking is H_9[], therefore, H_2[] is the left peak descriptor, and H_9[] is the right peak descriptor. The starting point of the retention time period corresponding to the known component group w is determined by accessing H_2[1], and the end point of the retention time period corresponding to the known component group w is determined by accessing H_9[1]. The range between the end point and the starting point is the retention time period corresponding to the known component group w.

[0053] In this embodiment, when the mass of each known component in the first detection is known, when verifying the detection result, the two adjacent components are merged into one known component group and the mass is integrated and calculated together. Therefore, when the two adjacent known components are in a multi-peak state, the secondary peak between the two adjacent known components does not need to be merged first as in Patent Document 1, thereby improving data processing efficiency.

[0054] In one embodiment, the method further comprises: Step S309 , subtract the sum of the masses of the known components constituting the known component cluster in the same atomized sample that has been tested from the mass of the known component cluster in the verification result of each atomized sample, to obtain a plurality of verification differences corresponding to the known component cluster.

[0055] Step S310, calculating the average of multiple verification differences corresponding to the known component clusters.

[0056] Step S311: If the average value is greater than a preset difference threshold, it indicates that an abnormal component exists.

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

[0058] It can be understood that if the number of aerosol samples is 10, there are 10 aerosol sample verification results, and 10 known component mass W masses are obtained. When these 10 aerosol samples are tested for the first time, the mass corresponding to the known component B and the known component C constituting the known component mass W is obtained. Therefore, the mass of the known component mass W in each aerosol sample verification result is subtracted from the sum of the mass of each known component constituting the known component mass in the same aerosol sample that has been tested (i.e., the mass corresponding to the known component B and the known component C) to obtain the verification difference corresponding to the known component mass W in the aerosol sample. By analogy, 10 verification differences corresponding to the known component mass W can be obtained, and the average value corresponding to the 10 verification differences is calculated. If the average value is greater than the preset difference threshold, it indicates that an abnormal component exists. In this embodiment, the average value obtained from multiple aerosol samples is used to judge the abnormal component, which can eliminate the interference caused by equipment problems when testing a single sample, and the test result is more accurate.

[0059] like Figure 7 As shown, in one embodiment, after indicating the presence of abnormal components, the method further includes: Step S400: If the average of multiple verification differences corresponding to the known component cluster is greater than a preset difference threshold, the known component cluster is used as a target known component cluster.

[0060] Step S401: taking the second window corresponding to the target known component cluster as the target second window.

[0061] Step S402: taking each known component in the third sequence located in the target second window as the target known component.

[0062] Step S403: taking the first window corresponding to each target known component as the target first window.

[0063] In this example, according to the third sequence, the target second window CK_D is obtained by merging two adjacent target first windows CK_B and CK_C that are grouped into the same group. Therefore, the target first window includes: the first windows CK_B and CK_C.

[0064] Step S404: determining information of abnormal components based on peak descriptors located between adjacent target first windows in the first sequence.

[0065] like Figure 1 As shown, after obtaining the peak descriptor located between adjacent target first windows in the first sequence, the retention time corresponding to the abnormal component R and the area of ​​the chromatographic peak can be determined, thereby comparing with the retention time corresponding to each component in the database to determine the component corresponding to the abnormal component R, and determining the mass of the abnormal component R through the area.

[0066] In the process of verifying the detection results of known components, this embodiment can also discover the existence of abnormal components and quickly locate the positions of the abnormal components, thereby improving data processing efficiency.

[0067] In one embodiment, Figure 8 As shown, 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: Step S502 , starting from the horizontal axis starting point coordinate of the preset range on the chromatogram, the acquisition window is slid in the positive direction of the horizontal axis with a step length of the acquisition window width minus one.

[0068] 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 known component, the width of the acquisition window can be adjusted according to the characteristics of the chromatographic peak. For example, the length of the first window is L, and the number of peaks (including main peaks and secondary 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, the length L of the first window of the known component B is 103, and the number of peaks N in the multi-peak state is 4. Substituting into the formula, the acquisition window width is (103-(103%4)) / 4=25. Through this formula, the possibility of two maximum or minimum values ​​appearing simultaneously in a collection window can be reduced as much as possible to avoid missing the collection of a certain maximum or minimum value. It should be noted that the number of peaks of the known component B in the multi-peak state can be obtained by artificially observing the same specification sample in advance. When determining the first window corresponding to each known component, the average value of the number of peaks corresponding to each known component in the multi-peak state is counted as N in the formula.

[0069] Step S503: when the acquisition window is at the starting point, a first extreme value acquisition operation is performed, and the monotonic characteristics of each part of the chromatogram segment within the current acquisition window are recorded in a temporary storage area.

[0070] The first extreme value acquisition operation includes: Search for the maximum and minimum points within the current acquisition window; like Figure 3 As shown in a, 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, the minimum value point is added to the extreme value point sequence and the current first extreme value acquisition operation is ended.

[0071] like Figure 3 As shown in b, 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, the current first extreme value acquisition operation is terminated.

[0072] like Figure 3 As shown in c, 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, the maximum value point is added to the extreme value point sequence and the current first extreme value acquisition operation is ended.

[0073] like Figure 3 As shown in d, if the maximum point and the minimum point are both 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, the minimum point and the maximum point are added to the extreme point sequence in sequence, and the current first extreme value acquisition operation is terminated.

[0074] like Figure 3 As shown in e, if the maximum point and the minimum point are both 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, the maximum point and the minimum point are added to the extreme point sequence in sequence, and the current first extreme value acquisition operation is terminated.

[0075] like Figure 3 As shown in f, 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 right boundary of the current acquisition window, the maximum value point is added to the extreme value point sequence and the current first extreme value acquisition operation is ended.

[0076] like Figure 3 As shown in g, 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, the current first extreme value acquisition operation is terminated.

[0077] like Figure 3As shown in h, 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, the minimum value point is added to the extreme value point sequence and the current first extreme value acquisition operation is ended.

[0078] Step S504, the acquisition window performs the second extreme value acquisition operation once each time it slides one step to a new position, until the acquisition window covers the horizontal axis end coordinate of the 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 the 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, 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.

[0079] It is understandable that when the acquisition window slides to the horizontal axis end point coordinate of the preset range on the chromatogram, the acquisition window may cover data points outside the preset range. At this time, the horizontal axis end point of the preset range on the chromatogram is used as the right boundary of the acquisition window, and the range of the acquisition window is re-determined, and then the second extreme value acquisition operation is performed within this range.

[0080] It can be understood that the starting position of the acquisition window is taken as the position when the acquisition window slides to step 0, and the new position reached after sliding 1 step from the starting position is taken as the position when the acquisition window slides to step 1, and the new position reached after sliding 2 steps from the starting position is taken as the position when the acquisition window slides to step 2, and the positions corresponding to other steps are similar and so on, which will not be repeated here.

[0081] The second extreme value acquisition operation includes: Step A, searching for the maximum value point and the minimum value point within the current acquisition window; wherein, the acquisition window is currently slid to the position of the i-th step; i is a positive integer variable.

[0082] Step B, based on the positional relationship between the maximum value point and the minimum value point in the current acquisition window, obtaining each maximum value point and each minimum value point on the chromatogram of the atomized sample to be detected.

[0083] Specifically, Figure 3 As shown, the method 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: like Figure 3As shown in a, 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, 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, add the minimum value point to the extreme value point sequence and end the current second 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 second extreme value acquisition operation.

[0084] like Figure 3 As shown in b, 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, 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, end the current second extreme value acquisition operation; if not, add the maximum value point to the extreme value point sequence and end the current second extreme value acquisition operation.

[0085] like Figure 3 As shown in c, 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, 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, add the maximum value point to the extreme value point sequence and end the current second 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 second extreme value acquisition operation.

[0086] like Figure 3 As shown in d, if the maximum point and the minimum point are both 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, add the minimum point and the maximum point to the extreme value point sequence in turn, 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 and the maximum point to the extreme value point sequence in turn, and end the current second extreme value acquisition operation.

[0087] like Figure 3As shown in e, if the maximum point and the minimum point are both 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, add the maximum point and the minimum point to the extreme value point sequence in turn, 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 and the minimum point to the extreme value point sequence in turn, and end the current second extreme value acquisition operation.

[0088] like Figure 3 As shown in f, 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 right boundary of the current acquisition window, 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, add the maximum value point to the extreme value point sequence and end the current second 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 turn, and end the current second extreme value acquisition operation.

[0089] like Figure 3 As shown in g, 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, 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, end the current second extreme value acquisition operation; if not, add the minimum value point to the extreme value point sequence and end the current second extreme value acquisition operation.

[0090] like Figure 3 As shown in h, 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, 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, add the minimum value point to the extreme value 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 value point to the extreme value point sequence in turn, and end the current second extreme value acquisition operation.

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

[0092] In one embodiment, Fig. 9As shown, the indicator code corresponding to the current acquisition window is determined 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, jump to step S16 and end the current process; if not, adjust to step S15, set j+1, and then jump to step S13. Among them, M is the number of elements in the indication set after deduplication.

[0093] by Figure 3 For example, the elements in the indicator set after deduplication are sorted from small to large according to their positions on the horizontal axis. They are: the point on the left boundary of the current acquisition window, the minimum point, and the point on the right boundary, a total of three data points. If the value of the second point on the vertical axis is subtracted from the value of the first point on the vertical axis, the result is less than 0, and if the value of the third point on the vertical axis is subtracted from the value of the second point on the vertical axis, the result is greater than 0, so the indicator code is 01.

[0094] In this embodiment, combined with the characteristics of the chromatogram, it is only necessary to calculate whether the difference of several points such as the points on the left boundary, the points on the right boundary, the minimum point and the maximum point of the current acquisition window is greater than zero to determine the monotonic characteristics of each part of the chromatogram segment. There is no need to perform derivative calculations, and the operation is simple and efficient.

[0095] In one embodiment, the method further comprises: 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, the leak detection operation is triggered.

[0096] The leak detection operation includes: 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, then the range between the two adjacent maximum values ​​is used as the range of the leak detection window, and the minimum value point is found in the leak detection window, and the minimum value point is inserted between the two adjacent maximum values ​​in the extreme value point sequence; If 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 the range between the two adjacent minimum value points is used as the range of the leak detection window, and the maximum value point is found in the leak detection window, and the maximum value point is inserted between the above two adjacent minimum value points in the extreme value point sequence.

[0097] It should be noted that when a collection window contains two maximum values ​​or minimum values ​​at the same time, one of the extreme values ​​may be missed. By adopting the solution of this embodiment, it is possible to detect the situation where the extreme value is missed, thereby improving the accuracy of extreme value collection.

[0098] like Figure 4 As shown, in one embodiment, the searching for the maximum value point and the minimum value point within the current acquisition window range specifically includes: First comparison process: Step S21, in order from small to large on the horizontal axis, obtain the first two points within the current acquisition window for comparison, store the point with the smaller value on the vertical axis into the first small value area, and store the point with the larger value on the vertical axis into the first large value area; let k=3; where k is a variable; and if the value of the second point on the vertical axis is greater than the value of the first point on the vertical axis, point the value address of the first reference point to the first large value area, and jump to step S22; if the value of the second point on the vertical axis is less than the value of the first point on the vertical axis, point the value address of the first reference point to the first small value area, and jump to step S24; Step S22, in order from small to large on the horizontal axis, obtain the kth point within the current acquisition window 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 large value area using the first target point, and jump to step S23; if q<p, end the first process and jump to the second comparison process; Step S23, determine whether k is greater than the number of data points in the current acquisition window. If not, add k+1 and jump to step S22; if yes, end the first comparison process and jump to the final comparison process; Step S24, in order from small to large on the horizontal axis, obtain the kth point within the current acquisition window 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 small 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; Step S25, determine whether k is greater than the number of data points in the current acquisition window. If not, add k+1 and jump to step S24; if yes, end the first process and jump to the final comparison process; Second comparison process: Step S31, in order from small to large on the horizontal axis, obtain the last two points in the current acquisition window range for comparison, store the point with the smaller value on the vertical axis into the first minimum value area, and store the point with the larger value on the vertical axis into the first maximum value area; let m=3; where m is a variable; and if the value of the second-to-last point on the vertical axis is greater than the value of the first-to-last point on the vertical axis, point the value address of the second reference point to the second maximum value area, and jump to step S32; if the value of the second-to-last point on the vertical axis is less than the value of the first-to-last point on the vertical axis, point the value address of the second reference point to the second minimum value area, and jump to step S34; Step S32, in order from small to large on the horizontal axis, obtain the mth point from the end within the current acquisition window as the second target point, and compare the value p of the second target point on the vertical axis with the value q of the second reference point on the vertical axis. If q≥p, update the second 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; Step S33, after setting m+1, jump to step S32; Step S34, in order from small to large on the horizontal axis, obtain the mth point from the end within the current acquisition window as the second target point, and compare the value p of the second target point on the vertical axis with the value q of the second reference point on the vertical axis. If q≤p, update the second 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; Step S35, after setting m+1, jump to step S34; Final comparison process: If both the second minimum value area and the second maximum value area are initial default values, the point in the first minimum value area is used as the minimum value point in the current acquisition window range, and the point in the first maximum value area is used as the maximum value point in the current acquisition window range; If neither the second minimum value area nor the second maximum value area is the initial default value, compare the values ​​of the points in the first minimum value area and the second minimum value area on the vertical axis, and take the smaller point as the minimum value point within the current acquisition window. If the values ​​of the points in the first minimum value area and the second minimum value area on the vertical axis are equal, take the point in the first minimum value area as the minimum value point within the current acquisition window. Compare the values ​​of the points in the first maximum value area and the second maximum value area on the vertical axis, and take the larger point as the maximum value point within the current acquisition window. If the values ​​of the points in the first maximum value area and the second maximum value area on the vertical axis are equal, take the point in the first maximum value area as the maximum value point within the current acquisition window.

[0099] The maximum and minimum value search algorithm provided in this embodiment utilizes the characteristics of the data points in the chromatogram in the acquisition window. Compared with the maximum and minimum value search algorithm in the prior art, which needs to compare 2n or 1.5n times to search for the maximum and minimum values, where n is the number of elements in the data set for which the maximum and minimum values ​​need to be searched, this embodiment only needs to compare n or n+2 times to obtain the maximum and minimum values ​​in n data. Its working principle is as follows: Figure 4 As shown, for n data points in the chromatogram segment that increases and then decreases in the acquisition window, the data points in the increasing part are only compared with the data points in the first large value area, and the data points in the decreasing part are only compared with the data points in the second large value area, so the total number of comparisons is n times. Finally, if both the second small value area and the second large value area are the initial default values, the comparison is invalid, so the total number of comparisons is n times. If neither the second small value area nor the second large value area is the initial default value, the first small value area is compared with the second small value area, and the first large value area is compared with the second large value area, for a total of n+2 comparisons.

[0100] In one embodiment, an electronic atomization component detection device is also provided, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the steps of the electronic atomization component detection method when executing the program. The steps of the electronic atomization component detection method here may be the steps in the above-mentioned embodiments.

[0101] In one embodiment, a computer-readable storage medium is further provided, wherein the computer-readable storage medium stores computer-executable instructions, wherein the computer-executable instructions are used to enable a computer to execute the steps of the above-mentioned electronic atomization component detection method. The steps of the electronic atomization component detection method here may be the steps of the electronic atomization component detection method in each of the above-mentioned embodiments.

[0102] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and / or volatile memory. 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 external cache memory. As an 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 (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRA), memory bus (Rambus) direct RAM (RDRA), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0103] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, 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, they should be considered to be within the scope of 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, the amount of substance of the known component is detected.

2. The method for detecting components of electronic atomized matter according to claim 1, characterized in that: 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.

3. The method for detecting components of electronic atomization according to claim 2, 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.

4. The method for detecting components of electronic atomized matter according to claim 3, 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.

5. The method for detecting components of electronic atomized matter according to claim 2, 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.

6. The method for detecting components of electronic atomized matter according to claim 2, 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.

7. 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.

8. The method for detecting components of electronic atomized matter according to claim 7, 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.

9. The method for detecting components of electronic atomized matter according to claim 7, 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.

10. 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 9 are performed.

Images