Data processing method for rapidly and accurately obtaining cigar temperature distribution based on temperature data of few detection points

Abstract

The invention discloses a data processing method for quickly and accurately obtaining cigar temperature distribution based on temperature data of a small number of detection points. The data processing method comprises the following steps: S1, collecting temperature data; s2, constructing a temperature curve of all grid points; and S3, obtaining temperature distribution. The method has the advantage that cigar temperature distribution can be quickly and accurately obtained by using temperature data of a small number of detection points acquired by repeated experiments. Suction time consistency is guaranteed by considering a suction interval translation temperature curve, and the efficient and fast processing process is guaranteed by combining the data processing process of translating and averaging original temperature data of multiple repeated experiments and supplementing the temperature curve of each depth interpolation point into one step.

Description

Technical Field

[0001] The invention relates to the field of data processing, and in particular to a data processing method for quickly and accurately obtaining the temperature distribution of a cigar based on temperature data of a small number of detection points. Background Art

[0002] In recent years, the sales volume of domestic cigars has increased year by year, the market scale has continued to expand, and the future development prospects are huge. Moreover, with the increase in cigar consumers, consumers' requirements for the quality of cigar smoke are constantly increasing. Therefore, improving the quality of cigar smoke and producing high-quality cigars that are popular with consumers is a future development trend.

[0003] Since cigars are made up of only about five raw materials, the raw materials have a great impact on product quality and style characteristics, especially the combustibility of the raw materials is an important factor affecting the smoke quality of cigars. At present, there are relatively few studies on the combustibility of cigar raw materials, especially the research on the temperature distribution of the burning cone of cigars is very weak.

[0004] The Chinese invention patent application with publication number CN117470400A discloses a method for rapid detection of the temperature distribution of a cigar during continuous smoking. The method uses a data processing method to perform a translational average of the original temperature data of multiple repeated experiments and fill in the temperature curve of each deep interpolation point on the temperature data of a small number of detection points collected by multiple repeated experiments to quickly obtain the temperature distribution of the burning cone of the cigar. Since the temperature of all grid points changes simultaneously when a cigar is smoked, the temperature curve should change simultaneously according to the period of the puff interval. However, this method does not take the value of the translational displacement of the temperature curve to an integer multiple of the puff interval when determining the translational displacement of the temperature curve. In addition, dividing the data processing process into two steps: translational average of the original temperature data of multiple repeated experiments and filling in the temperature curve of each deep interpolation point is relatively cumbersome and time-consuming. Summary of the invention

[0005] The purpose of the present invention is to provide a data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on the temperature data of a small number of detection points. The method can quickly and accurately obtain the temperature distribution of cigar smoke using the temperature data of a small number of detection points collected by repeated experiments. The consistency of the puffing time is ensured by considering the puffing interval to translate the temperature curve, and the data processing process of translating the average of the original temperature data of repeated experiments and filling the temperature curve of each deep interpolation point is combined into one step to ensure that the processing process is efficient and fast. To achieve the above purpose, the technical solution provided by the present invention is:

[0006] A data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data of a small number of detection points is provided, the method comprising the following steps:

[0007] S1: Collecting temperature data: Collecting the temperature data of each measuring point in the detection grid of the temperature distribution of the cigar smoke, performing repeated experiments for many times and ensuring that the puffing parameters are consistent in each experiment, and recording the puffing interval;

[0008] S2: Constructing the temperature curves of all grid points: constructing the temperature curves of all grid points in the cigar temperature distribution detection grid according to the temperature data collected by repeated experiments;

[0009] S3: Obtaining temperature distribution: The temperature data of the temperature curve of all grid points at a certain moment are interpolated by the Makima interpolation method in Matlab software to obtain the temperature distribution of the cigar burning cone at that moment.

[0010] Furthermore, the specific steps of collecting temperature data described in S1 include:

[0011] S11: Set up a detection grid for the temperature distribution of cigars: divide the cigars into several depths along the radial direction, set no less than 14 detection points in the axial direction at each depth, and the interval is no less than 2mm. The detection axial position is selected between 30 and 60mm from the burning end of the cigar where the density of the cigarettes is relatively uniform. In actual tests, in order to ensure that temperature data of all depths are obtained in one experiment, y (2≦y≦4 and y is an integer) actual measurement points are selected at each depth to insert thermocouples, and the detection points selected for inserting thermocouples at different depths cannot be in the same radial direction. The detection points where thermocouples are inserted in all the set detection points are called actual measurement points, and the detection points where thermocouples are not inserted are called interpolation points. In this way, the actual measurement points and interpolation points are evenly distributed on the longitudinal section of the cigar to form a number of grid nodes;

[0012] S12: Collecting temperature: When testing the temperature distribution of cigar smoke, maintain the experimental temperature at 22±1°C and the humidity at 60±5% to achieve temperature and humidity balance, light the cigar with a handheld lighter, start the smoking machine to continuously smoke the cigar smoke and continuously collect the temperature data of each measuring point in real time, and record the smoking interval used in the experimental process as T, in seconds;

[0013] S13: Repeat the experiment: Repeat the steps of S12 above n times to collect the temperature data of each measuring point, and ensure that the experimental samples, experimental environment, suction parameters, experimental equipment, etc. are consistent in each repeated experiment.

[0014] Further, the specific steps of constructing the temperature curves of all grid points described in S2 include:

[0015] S21: Determine the translation displacement of two adjacent grid points:

[0016] According to the above S11, the detection grid of the cigar temperature distribution is set to perform n repeated experiments, and n groups of temperature curves of the actual measurement points at the i-th depth y can be obtained.

[0017] Taking the temperature curve of one of the measured points j at the i-th depth of the k-th repeated experiment as the reference, the temperature curves of the remaining measured points at the i-th depth of the k-th repeated experiment are translated to the measured point j according to the maximum correlation coefficient method, and the translation displacements p1, p2…, p y-1 , and the corresponding intervals of the remaining measured points at the i-th depth from the measured point j are h1, h2…, h y-1 , where the interval number h refers to the number of adjacent grid point intervals between two grid points at the same depth on the detection grid of cigar smoke temperature distribution.

[0018] According to the displacement p of the remaining measured points at the i-th depth in the k-th repeated experiment according to the maximum correlation coefficient method and the interval number h between the remaining measured points at the i-th depth and the measured point j, the displacement required for the remaining measured points at the i-th depth to translate one interval can be obtained. Denoted as l1, l2…, l y-1 , and then taking the average value of l can obtain the translation displacement of the interval between adjacent grid points at the i-th depth in the k-th repeated experiment:

[0019] For each repeated experiment, the translation displacement L of the adjacent grid points at the i-th depth is k Taking the average and taking it as an integer multiple of the suction interval T, the translation displacement of two adjacent grid points at the i-th depth can be obtained:

[0020] S22: Determine the translation reference curve:

[0021] The temperature curve N of a measured point j at the i-th depth in the k-th repeated experiment ij As a benchmark, the temperature curve of the yth measured point at the i-th depth of each repeated experiment is translated to the measured point j at the i-th depth of the k-th repeated experiment according to the maximum correlation coefficient method, and the translation displacement is recorded as f1, f2…, f β ,β=n*y-1. Then the translation displacement f is taken as an integer multiple of the suction interval T and recorded as F, then Then the temperature curve of the y measured points at the i-th depth in each repeated experiment is F1, F2…, F β ,β=n*y-1, press F1, F2…, F β The displacement of the temperature curve of the y measured points at the i-th depth of each repeated experiment is translated to the measured point j at the i-th depth of the k-th repeated experiment. In this way, there are n*y temperature curves at the measured point j at the i-th depth of the k-th repeated experiment. The n*y temperature data of these n*y temperature curves at the same time are averaged after eliminating abnormal values ​​through double standard deviation, and the translation reference curve N at the measured point j at the i-th depth can be obtained.ij ;

[0022] S23: Construct the temperature curve of all grid points:

[0023] The translation displacement L of two adjacent grid points at the i-th depth and the translation reference curve N at the measured point j at the i-th depth are known. ij , the distance between a grid point s at the i-th depth and the measured point j is recorded as h s , if the grid point s is to the left of the measured point j, then the reference curve N will be translated ij Move ΔX = L*h to the left along the horizontal axis; if the grid point s is to the right of the measured point j, the reference curve N will be translated ij Move rightward along the horizontal axis by ΔX=L*h. In this way, the temperature curve of all grid points at the i-th depth can be obtained.

[0024] The invention has the advantages that the temperature distribution of cigar smoke can be obtained quickly and accurately by using the temperature data of a small number of detection points collected by repeated experiments. The consistency of the puffing time is ensured by considering the puffing interval to translate the temperature curve, and the data processing process of translating the average of the original temperature data of repeated experiments and filling the temperature curve of each deep interpolation point is combined into one step to ensure the processing process is efficient and fast. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 A schematic diagram of the overall flow of a data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data of a small number of detection points provided by the present invention;

[0026] Figure 2 A detection grid for cigar temperature distribution set in a specific implementation manner;

[0027] Figure 3 It is the temperature curve graph of the 8th depth 8 in the specific implementation mode;

[0028] Figure 4 is a temperature curve diagram of all grid points in a specific implementation manner;

[0029] Figure 5 1 is the temperature distribution inside the combustion cone before and after continuous puffing of the cigar in a specific embodiment. DETAILED DESCRIPTION

[0030] The present invention will further explain a data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data of a small number of detection points provided by the present invention in conjunction with the accompanying drawings below:

[0031] Example

[0032] Figure 1The overall flow chart of a data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data of a small number of detection points provided by the present invention is as follows: Figure 1 As shown, the overall flow diagram of a data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data of a small number of detection points provided by the present invention includes:

[0033] S1: Collecting temperature data: Collect the temperature data of each measuring point in the detection grid of the temperature distribution of the cigar smoke, repeat the experiment 4 times and ensure that the puffing parameters are consistent in each experiment, and record the puffing interval T = 40s;

[0034] S2: Constructing the temperature curves of all grid points: Constructing the temperature curves of all grid points in the cigar temperature distribution detection grid based on the temperature data collected from 4 repeated experiments;

[0035] S3: Obtain temperature distribution: Figure 5 As shown, the temperature data of 0s, 0.5s, 1s, 1.5s, and 2s of the temperature curve of all grid points are interpolated by the makima interpolation method in Matlab software to obtain the temperature distribution inside the combustion cone before and after continuous puffing of the cigar.

[0036] In this embodiment, refer to Figure 2 The specific steps of collecting temperature data described in S1 include:

[0037] S11: Set up a detection grid for the temperature distribution of cigars: Divide the cigars into 8 equal depths along the radial direction, set 16 detection points in the axial direction of each depth, with an interval of 2mm, and select the detection axial position between 30 and 60mm from the burning end of the cigar where the density of the cigarettes is relatively uniform. In actual testing, in order to ensure that temperature data of all depths are obtained in one experiment, two actual measurement points are selected at each depth to insert thermocouples, and the detection points selected for inserting thermocouples at different depths cannot be in the same radial direction. The detection points where thermocouples are inserted in all the set detection points are called actual measurement points, and the detection points where thermocouples are not inserted are called interpolation points. Figure 2 As shown, the measured points and interpolation points are evenly distributed on the longitudinal section of the cigar to form a number of grid nodes;

[0038] S12: Collecting temperature: When testing the temperature distribution of cigar smoke, maintain the experimental temperature at 22±1°C and the humidity at 60±5% to achieve temperature and humidity balance, light the cigar with a handheld lighter, start the smoking machine to continuously smoke the cigar smoke and continuously collect the temperature data of each measuring point in real time, and record the smoking interval used in the experimental process as 40 seconds;

[0039] S13: Repeat the experiment: Repeat the steps of S12 above 4 times to collect the temperature data of each measuring point, and ensure that the experimental samples, experimental environment, suction parameters, experimental equipment, etc. are consistent in each repeated experiment.

[0040] In this embodiment, see Figure 3 and Figure 4 The specific steps of constructing the temperature curves of all grid points described in S2 include:

[0041] S21: Determine the translation displacement of two adjacent grid points:

[0042] The specific processing steps of S21 are taken at the 8th depth as an example, and the processing methods of other depths are the same as that of the 8th depth.

[0043] According to the above S11, the detection grid of the cigar temperature distribution is set to repeat the experiment 4 times, and 4 groups of temperature curves of 2 measuring points at the 8th depth can be obtained.

[0044] Taking the temperature curve of the first measured point j at the 8th depth close to the burning end of the cigar in the first repeated experiment as a reference, the temperature curves of the remaining measured points at the 8th depth of the first repeated experiment were translated to the measured point j according to the maximum correlation coefficient method, and the translation displacement p1=718.3 of the remaining measured points at the 8th depth and the corresponding interval number h1=8 between the remaining measured points at the 8th depth and the measured point j were recorded.

[0045] According to the displacement p1 of the remaining measured points at the 8th depth in the first repeated experiment according to the maximum correlation coefficient method and the number of intervals h1=8 between the remaining measured points at the 8th depth and the measured point j, the displacement required for the remaining measured points at the 8th depth to translate one interval can be obtained. Let l1 = 89.8, then take the average value of l to get the translation displacement of the interval between adjacent grid points at the 8th depth in the first repeated experiment:

[0046] The translation displacement L of the adjacent grid points at the 8th depth in the 2nd, 3rd and 4th repetition experiments k The calculation method is the same as the first repeated experiment mentioned above. After calculation, it can be obtained that L2=84.9, L3=87.6, L4=94.3.

[0047] The translation displacement L of the interval between adjacent grid points at the 8th depth for each repeated experiment is k Taking the average and taking it as an integer multiple of the suction interval T, the translation displacement between two adjacent grid points at the 8th depth can be obtained:

[0048] S22: Determine the translation reference curve:

[0049] The specific processing steps of S22 are taken at the 8th depth as an example, and the processing methods of other depths are the same as that of the 8th depth.

[0050] The temperature curve N of one of the measured points j at the 8th depth in the first repeated experiment 8j As a benchmark, the temperature curves of the two measured points at the 8th depth of each repeated experiment were translated to the measured point j at the 8th depth of the first repeated experiment according to the maximum correlation coefficient method, and the translation displacements were recorded as f1=358.3, f2=255.3, f3=541.9, f4=718.3, f5=1118.1, f6=974.9, f7=1301.3.

[0051] The translation displacement f is taken as an integer multiple of the suction interval T and recorded as It can be obtained that F1=360.0, F2=240.0F3=560.0, F4=720.0, F5=1120.0, F6=960.0, F7=1320.0. The temperature curves of the two measured points at the 8th depth of each repeated experiment are translated to the measured point j at the 8th depth of the first repeated experiment according to the corresponding translation displacement F, as shown in Figure 3 As shown in the figure, there are 8 temperature curves at the measured point j at the 8th depth in the first repeated experiment. The 8 temperature data of these 8 temperature curves at the same time are averaged after eliminating abnormal values ​​through double standard deviation, and the translation reference curve N at the measured point j at the 8th depth can be obtained. 8j ;

[0052] S23: Construct the temperature curve of all grid points:

[0053] The processing step of S23 is taken at the 8th depth as an example, and the processing methods of other depths are the same as that of the 8th depth.

[0054] The translation displacement L = 80 between two adjacent grid points at the 8th depth and the translation reference curve N at the measured point j at the 8th depth are known. 8j ,like Figure 2 As shown, taking the grid point s at the 8th depth as an example, the interval number h between the grid point s and the measured point j is s = 3 and the grid point s is to the right of the measured point j, then the reference curve N will be translated 8j Move right along the horizontal axis by ΔX=80*3=240. The processing method for the remaining grid points at the 8th depth is the same as that for grid point s. In this way, the temperature curves of all grid points at the 8th depth can be obtained. The processing method for the remaining depths is the same as that for the 8th depth, such as Figure 4 As shown, the temperature curves of all grid points at each depth can be obtained.

Claims

1. A data processing method for quickly and accurately obtaining the temperature distribution of cigar smoke based on temperature data from a small number of detection points, characterized in that: The following steps are involved: S1: Collecting temperature data: Collecting the temperature data of each measuring point in the detection grid of the temperature distribution of the cigar smoke, performing repeated experiments for many times and ensuring that the puffing parameters are consistent in each experiment, and recording the puffing interval; S2: Constructing the temperature curves of all grid points: constructing the temperature curves of all grid points in the cigar temperature distribution detection grid according to the temperature data collected by repeated experiments; S3: Obtaining temperature distribution: The temperature data of the temperature curve of all grid points at a certain moment are interpolated by the Makima interpolation method in Matlab software to obtain the temperature distribution of the cigar burning cone at that moment.

2. The data processing method according to claim 1, characterized in that: The specific steps of collecting temperature data described in S1 include: S11: Setting a detection grid for the temperature distribution of cigars: Divide the cigars into several depths along the radial direction, set no less than 14 detection points in the axial direction of each depth, and the interval is no less than 2mm. The detection axial position is selected between 30 and 60mm from the burning end of the cigar where the density of the cigarettes is relatively uniform; in the actual test, in order to ensure that the temperature data of all depths are obtained in one experiment, y actual measurement points are selected at each depth to insert thermocouples, and the detection points selected for inserting thermocouples at different depths cannot be in the same radial direction, and y is an integer 2≦y≦4; the detection points where thermocouples are inserted in all the set detection points are called actual measurement points, and the detection points where thermocouples are not inserted are called interpolation points; in this way, the actual measurement points and interpolation points are evenly distributed on the longitudinal section of the cigar to form a number of grid nodes; S12: Collecting temperature: When testing the temperature distribution of cigar smoke, maintain the experimental temperature at 22±1°C and the humidity at 60±5% to achieve temperature and humidity balance, light the cigar with a handheld lighter, start the smoking machine to continuously smoke the cigar smoke and continuously collect the temperature data of each measuring point in real time, and record the smoking interval used in the experimental process as T, in seconds; S13: Repeat the experiment: Repeat the steps of S12 above n times to collect the temperature data of each measuring point, and ensure that the experimental samples, experimental environment, suction parameters and experimental equipment are consistent in each repeated experiment.

3. The data processing method according to claim 1, characterized in that: The specific steps of constructing the temperature curves of all grid points in S2 include: S21: Determine the translation displacement of two adjacent grid points: According to the above S11, the detection grid of the cigar temperature distribution is set to repeat the experiment n times, and n groups of temperature curves of the measured points at the i-th depth y can be obtained; Taking the temperature curve of one of the measured points j at the i-th depth of the k-th repeated experiment as the reference, the temperature curves of the remaining measured points at the i-th depth of the k-th repeated experiment are translated to the measured point j according to the maximum correlation coefficient method, and the translation displacements p1, p2…, p y-1 , and the corresponding intervals of the remaining measured points at the i-th depth from the measured point j are h1, h2…, h y-1 , where the interval number h refers to the number of adjacent grid point intervals between two grid points at the same depth on the detection grid of the cigar smoke temperature distribution; According to the displacement p of the remaining measured points at the i-th depth in the k-th repeated experiment according to the maximum correlation coefficient method and the interval number h between the remaining measured points at the i-th depth and the measured point j, the displacement required for the remaining measured points at the i-th depth to translate one interval can be obtained. Denoted as l1, l2…, l y-1 , and then taking the average value of l can obtain the translation displacement of the interval between adjacent grid points at the i-th depth in the k-th repeated experiment: For each repeated experiment, the translation displacement L of the adjacent grid points at the i-th depth is k Taking the average and taking it as an integer multiple of the suction interval T, the translation displacement of two adjacent grid points at the i-th depth can be obtained: S22: Determine the translation reference curve: The temperature curve N of a measured point j at the i-th depth in the k-th repeated experiment ij As a benchmark, the temperature curve of the yth measured point at the i-th depth of each repeated experiment is translated to the measured point j at the i-th depth of the k-th repeated experiment according to the maximum correlation coefficient method, and the translation displacement is recorded as f1, f2…, f β ,β=n*y-1; and then the translation displacement f is taken as an integer multiple of the suction interval T and recorded as F, then Then the temperature curve of the y measured points at the i-th depth in each repeated experiment is F1, F2…, F β ,β=n*y-1, press F1, F2…, F β The displacement of the temperature curve of the y measured points at the i-th depth of each repeated experiment is translated to the measured point j at the i-th depth of the k-th repeated experiment. In this way, there are n*y temperature curves at the measured point j at the i-th depth of the k-th repeated experiment. The n*y temperature data of these n*y temperature curves at the same time are averaged after eliminating abnormal values ​​through double standard deviation, and the translation reference curve N at the measured point j at the i-th depth can be obtained. ij ; S23: Construct the temperature curve of all grid points: The translation displacement L of two adjacent grid points at the i-th depth and the translation reference curve N at the measured point j at the i-th depth are known. ij , the distance between a grid point s at the i-th depth and the measured point j is recorded as h s , if the grid point s is to the left of the measured point j, then the reference curve N will be translated ij Move ΔX = L*h to the left along the horizontal axis; if the grid point s is to the right of the measured point j, the reference curve N will be translated ij Move rightward along the horizontal axis by a displacement of ΔX=L*h; in this way, the temperature curve of all grid points at the i-th depth can be obtained.

Images