CN108828132A

CN108828132A - The analysis method of Contents of Main Components in a kind of cigarette smoke

Info

Publication number: CN108828132A
Application number: CN201810653885.3A
Authority: CN
Inventors: 吴晓炯; 罗佳; 汪旭
Original assignee: Shanghai Tobacco Group Co Ltd
Current assignee: Shanghai Tobacco Group Co Ltd
Priority date: 2018-06-22
Filing date: 2018-06-22
Publication date: 2018-11-16
Anticipated expiration: 2038-06-22
Also published as: CN108828132B

Abstract

The present invention provides a kind of analysis method of Contents of Main Components in cigarette smoke, includes the following steps：Step S1 determines the underlying parameter of cigarette；Step S2 calculates the initial parameter of cigarette according to the underlying parameter of cigarette in step S1；Step S3 calculates the relative parameter of cigarette according to the initial parameter of cigarette in the underlying parameter of cigarette in step S1 and step S2；Step S4 obtains Contents of Main Components value in cigarette smoke to calculate according to the relative parameter of cigarette in the initial parameter of cigarette in step S2 and step S3.The analysis method of Contents of Main Components in a kind of cigarette smoke provided by the invention, it can apply in computer, calculate the influence of formula, auxiliary material and cigarette structure variation to main flume, can it is accurate calculate cigarette smoke in main component content, without being transformed to existing production of cigarettes process equipment, without carrying out trial test early period of a large amount of repeatability, there is easily and efficiently advantage.

Description

The analysis method of Contents of Main Components in a kind of cigarette smoke

Technical field

The invention belongs to technical field of cigarettes, it is related to a kind of analysis method of Contents of Main Components in cigarette smoke.

Background technique

Accelerate the desin speed of product from the fifties in last century in order to improve the quality of cigarette product, optimizes the body of product It tests, has carried out a large amount of correlative study both at home and abroad.External cigarette enterprise is started with from the physical chemistry rule of the primary mostly, thin Breathe freely rule, filter stick filter capacity, the combustion mechanism of cigarette etc. of point field such as cigarette paper has carried out correlative study, forms pair The whole understanding of cigarette burning filtration mechanism；Domestic research is generally basede on statistical test, analyzes some specific auxiliary materials, rolls The combination of parameter reaches the forecast function of certain main flume.But since the factor for influencing cigarette smoke is too many, both The combustion law of entire cigarette an entirety is not combined into, also without forming the analysis method based on cigarette population parameter.

Specifically, cigarette smoke composition and its content are by pipe tobacco component, cigarette structure, cigarette material parameter, burning item The joint effect of the factors such as part is related to leaf tobacco production, pipe tobacco processing, multiple links such as roll equipment, product design.For standard Really prediction cigarette mainstream flue gas, reaches the target of Cigarette design, a large amount of research work has been done in industry. CN105628646 discloses a kind of online tar prediction of cigarette and method for early warning, this method utilize infrared light collection pipe tobacco ingredient, The prediction model for constructing tar content calculates upper lower limit value by factorial analysis and forms threshold value of warning.But this method is by environment Influence is more significant, and measurement accuracy deviation is larger in data procedures.CN106108107 disclose it is a kind of based on cigarette structure with The modeled design method of the total ventilation rate of the cigarette of physical parameter, filter ventilation rates and resistance to suction with cigarette structure characteristic parameter and Cigarette physical property characteristic parameter is variable, using the total ventilation rate of cigarette, filter ventilation rates and resistance to suction as objective function, by adjusting Cigarette structure and physical parameter obtain a kind of optimization method for meeting Cigarette design requirement.But this method only obtains cigarette ventilation Rate and resistance to suction can not study the content of cigarette smoke main component deeper into ground.

Therefore index used in the Cigarette designs such as tobacco leaf formulation, cigarette auxiliary material and wrapped configuration how is integrated, and is formed Physical method based on the primary has substantial worth to Cigarette design, this also has realistic meaning to reduction entreprise cost.In addition, Using single factor design technique study cigarette smoke composition and content, cannot generally meet well cigarette smoke composition and The design requirement of content has codetermined cigarette production because the tobacco leaf formulation of cigarette and cigarette structure parameter are interactions Product.So being based on tobacco leaf formulation and wrapped configuration these two aspects, the analysis method of building cigarette smoke composition and content is very It is necessary.

Summary of the invention

In view of the foregoing deficiencies of prior art, the purpose of the present invention is to provide main components in a kind of cigarette smoke The analysis method of content, the algorithm model based on cigarette underlying parameter and relevant parameter can efficiently and accurately analyze cigarette cigarette The content of main component in gas improves Cigarette design efficiency to make further adjustment to tobacco leaf formulation and wrapped configuration.

In order to achieve the above objects and other related objects, the present invention provides a kind of point of Contents of Main Components in cigarette smoke Analysis method, includes the following steps：

Step S1 determines the underlying parameter of cigarette；

Preferably, in step S1, the underlying parameter of the cigarette includes but is not limited to pipe tobacco type, pipe tobacco weight, cigarette Segment length, filter lengths, filter stick pressure drop, tipping paper air permeability, tipping paper pore area overall width, tipping paper width, connects Cigarette circumference The outer pore area of paper is filled away from side, forming paper air permeability, air permeance for cigarette paper, the combustion-supporting agent content of cigarette paper, cigarette paper grammes per square metre, nicotine mistake Filter coefficient, tar filtration coefficient.The underlying parameter of the cigarette is the routine data of cigarette, can measure and obtains.Step S2, is adopted With smoldering time recurrence method, the initial parameter of cigarette is calculated according to the underlying parameter of cigarette in step S1；

Preferably, in step S2, the initial parameter of the cigarette include but is not limited to original smoldering time, original mouth number, Every mouthful of burning original length, every mouthful of burning original weight, every mouthful of burning original residual length.

Preferably, in step S2, the calculating of the initial parameter of the cigarette includes the following steps：

Step S21 after determining pipe tobacco type, calculates original smoldering time by formula (1), the formula (1) is：

T1=73.3 × l1+279 × c-33.1 × a-19.5 × mc+6.86 × ps+11.7 × pe+6.49 × pp+10.4 × pl,

Wherein, t1 is original smoldering time, s；L1 is cigarette segment length, cm；C is Cigarette circumference, cm；A helps for cigarette paper Fire agent content, %；Mc is cigarette paper grammes per square metre, g/m²；Ps is stem ratio, %；Pe is swollen silk ratio, %；Pp is thin slice ratio Example, %；Pl is cut tobacco ratio, %.

Stem ratio described above, swollen silk ratio, thin slice ratio, cut tobacco ratio are the setting numerical value of cigarette, are belonged to Fixed data parameter value obtained by designer.

Step S22 calculates original mouth number by formula (2) by original smoldering time, and the formula (2) is：

N=0.0121 × t1+0.487,

Wherein, n is original mouth number, mouth；T1 is original smoldering time, s.

Step S23, it is original by every mouthful of formula (3), (4), (5) calculating burning by original smoldering time and original mouth number Length, the formula (3) are：

V=l1/t1,

Wherein, v is smoldering velocity, cm/s；L1 is cigarette segment length, cm；T1 is original smoldering time, s；

The formula (4) is：

T2=58 × n,

Wherein, t2 is practical smoldering time, s；N is original mouth number, mouth；

The formula (5) is：

L2=v × (t1-t2)/n,

Wherein, l2 is every mouthful of burning original length, cm；V is smoldering velocity, cm/s；T1 is original smoldering time, s；T2 is Practical smoldering time, s；N is original mouth number, mouth.

Step S24 calculates every mouthful of burning original weight, the formula by formula (6) by every mouthful of burning original length (6) it is：M2=m1 × l2/l1,

Wherein, m2 is every mouthful of burning original weight, g；M1 is pipe tobacco weight, g；L2 is every mouthful of burning original length, cm；l1 For cigarette segment length, cm.

It is highly preferred that every mouthful of burning original weight is calculated separately the original calculation weight of nicotine by formula (a), (b) Amount, the original calculation weight of tar,

The formula (a) is：M3=m2 × nic,

Wherein, m3 is the original calculation weight of nicotine, g；M2 is every mouthful of burning original weight, g；Nic is the Buddhist nun of pipe tobacco Ancient fourth empirical coefficient, is set as 1.11.The nicotine empirical coefficient nic of the pipe tobacco is to adjust different leaf group nicotine deliveries Coefficient；

The formula (b) is：M4=m2 × tar,

Wherein, m4 is the original calculation weight of tar, g；M2 is every mouthful of burning original weight, g；Tar is the tar of pipe tobacco Empirical coefficient is set as 1.22.The tar empirical coefficient tar is the coefficient for adjusting different leaf group releasing content of coke tar.

Step S25 calculates every mouthful of burning original residual length, institute by formula (7), (8) by every mouthful of burning original length Stating formula (7) is：L3 [1]=l1-l2,

Wherein, every mouthful of burning original residual length that l3 [1] is the 1st mouthful, cm；L1 is cigarette segment length, cm；L2 is every mouthful Burn original length, cm；

The formula (8) is：L3 [i]=l3 [i-1] -58v-l2,

Wherein, every mouthful of burning original residual length that l3 [i] is i-th mouthful, i=2 ..., n, cm；L3 [i-1] is (i-1)-th mouthful Every mouthful of burning original residual length, i=2 ..., n, cm；V is smoldering velocity, cm/s；L2 is every mouthful of burning original length, cm.

The formula (7), (8) are iterative algorithm.

Preferably, in step S2, the smoldering time recurrence method refers to, (only and cigarette by the original smoldering time of cigarette Silk type, wrapped configuration, cigarette paper are related) the original mouth number of cigarette burning is calculated, and cigarette combustion is further calculated out with this Every mouthful of burning original length and every mouthful of burning original weight when burning, another step obtain the nicotine and the original meter of tar of every mouthful of cigarette It calculates weight (assuming that cigarette shreds are uniformly distributed), and calculates every mouthful of burning original residual length for obtaining every mouthful of cigarette.

Preferably, in step S2, the initial parameter of the cigarette is the parameter calculated under the conditions of ignoring tipping paper and punching. The tipping paper punching condition refers to the case where hole of tipping paper is blocked with adhesive plaster.

Step S3 matches regression algorithm using pressure drop, according to cigarette in the underlying parameter of cigarette in step S1 and step S2 Initial parameter calculate the relative parameter of cigarette；

Preferably, in step S3, the relative parameter of the cigarette includes but is not limited to that mouth number, every mouthful of burning relatively are relatively surplus Remaining length, every mouthful of burning relative weight.

Preferably, in step S3, the relative parameter of the cigarette is the parameter calculated under the conditions of tipping paper punching.

Preferably, in step S3, the calculating of the relative parameter of the cigarette includes the following steps：

Step S31 calculates opposite mouth number by formula (9), and the formula (9) is：

Rn=l1/ (l6avg+58v),

Wherein, rn is opposite mouth number, mouth；L1 is cigarette segment length, cm；L6avg is average relative combustion length, cm；V is Smoldering velocity, cm/s.

It is highly preferred that the average relative combustion length is calculated by formula (c), the formula (c) is：

L6avg=Simga ( l6 [ i ] ) /n,

Wherein, l6avg is average relative combustion length, cm；Simgais an add sign∑,is the arrayl6 [ i ] summation； L6 [i] is the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；N is original mouth number, mouth.

It is further preferred that the relative combustion length after i-th mouthful of ventilation is calculated by formula (d), the formula (d) For：L6 [i]=qi [i]/17.5 × l2,

Wherein, l6 [i] is the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；Qi [i] is the volume at i-th mouthful Cigarette combustion inlet flow velocity qi value, i=1 ..., n, ml/s；L2 is every mouthful of burning original length, cm.

It is further preferred that the cigarette burning inlet flow rate qi value at i-th mouthful is calculated by formula (e), the public affairs Formula (e) is：Qi [i]=0.9 × (17.5-q [i]),

Wherein, qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, i=1 ..., n, ml/s；Q [i] is i-th The ventilation flow rate q value of filter tip when length l3 [i] of mouth, i=1 ..., n, ml/s.

Still further preferably, the flow velocity identity in i-th mouthful of the cigarette burning inlet flow rate coincidence formula (f), The formula (f) is：Qi [i]+qc [i]+q [i]=17.5,

Wherein, qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, i=1 ..., n, ml/s；Qc [i] is i-th The flow velocity qc value for crossing cigarette paper of mouth, i=1 ..., n, ml/s；Q [i] is the ventilation of the filter tip in i-th mouthful of length l3 [i] Flow velocity q value, i=1 ..., n, ml/s.

Still further preferably, the flow speed of the filter tip is calculated by the recurrence pressure drop matching equation in formula (g), institute Stating formula (g) is：

Q=FindRoot [ptpw [q]=pr [q]+pre [q]+pf [q], q->0],

Wherein, q is the flow speed by filter tip, ml/s；FindRoot is to pass through initial value (q->0) equation is searched for Numerical solution, ml/s；q->0 refers to that q is intended to 0；Ptpw [q] is that tipping paper and forming paper combine the function of pressure drop at flow velocity q, cmH₂O；Pr [q] is that cigarette section considers function of pressure drop of the air permeance for cigarette paper at flow velocity q, cmH₂O；Pre [q] is flow velocity when being q Transposition section cigarette function of pressure drop, cmH₂O；The filter stick function of pressure drop that pf [q] is flow velocity when being q, cmH₂O。

It is particularly preferred that the function of pressure drop of the tipping paper and forming paper combination at flow velocity q calculates tipping by formula (h) Paper-forming paper pressure drop-flow function, the cigarette section consider that function of pressure drop of the air permeance for cigarette paper at flow velocity q presses formula (i) Pressure drop-flow function under ventilation condition containing cigarette paper is calculated, the transposition section cigarette function of pressure drop when flow velocity is q presses formula (j) pressure drop-flow function without cigarette under cigarette paper ventilation condition is calculated, the filter stick function of pressure drop when flow velocity is q is by public affairs Formula (k) calculates pressure drop-flow function of filter stick,

The formula (h) is：Ptpw [q]=(q/c)²×(2.56×10⁸/a2²+3.6×10⁴/a3²),

Wherein, ptpw [q] is that tipping paper and forming paper combine the function of pressure drop at flow velocity q, cmH₂O；Q is to pass through filter tip Flow speed, ml/s；C is Cigarette circumference, cm；A2 is forming paper air permeability, CU；A3 is tipping paper air permeability, CU；

The formula (i) is：

Pr [q]=pre0 × Tanh (Sqrt (a1 × c × l3 [i] × pre0 × (7 × 10^-5+2.9×10^-6×(17.5- q)))/Sqrt(a1×c×l3[i]×pre0×(7×10^-5+2.9×10^-6× (17.5-q))),

Wherein, pr [q] is that cigarette section considers function of pressure drop of the air permeance for cigarette paper at flow velocity q, cmH₂O；Pre0 is cigarette Branch section ignores the voltage drop value of air permeance for cigarette paper, cmH₂O；Tanh is hyperbolic tangent function；Sqrt is radical sign function；A1 is cigarette Paper air permeability, CU；C is Cigarette circumference, cm；Every mouthful of burning original residual length that l3 [i] is i-th mouthful, i=2 ..., n, cm；q To pass through the flow speed of filter tip, ml/s；

The formula (j) is：

Pre [q]=0.143 × (1-E)²×l5×(17.5-q)/E³/ s,

Wherein, the transposition section cigarette function of pressure drop that pre [q] is flow velocity when being q, cmH₂O；E is pipe tobacco porosity, is set as 0.76；L5 is the length of intersecting of tipping paper and cigarette paper, cm；S is cigarette sectional area, cm²；

The formula (k) is：Pf [q]=k × (17.5-q) × l4,

Wherein, the filter stick function of pressure drop that pf [q] is flow velocity when being q, cmH₂O；K is filter stick pressure-drop coefficient, cmH2O*s/cm³； Q is the flow speed by filter tip, ml/s；L4 is that tipping paper punches leading portion filter lengths, cm.The filter stick pressure-drop coefficient is root It calculates and obtains according to filter stick pressure drop and filter lengths, each filter stick unit all has corresponding filter stick pressure drop and filter lengths.

Most preferably, the voltage drop value that the cigarette section ignores air permeance for cigarette paper is calculated by formula (l), the formula (l) For：Pre0=2.5 × (1-E)²×l3[i]/E³/ s,

Wherein, pre0 is the voltage drop value that cigarette section ignores air permeance for cigarette paper, cmH₂O；E is pipe tobacco porosity, is set as 0.76；Every mouthful of burning original residual length that l3 [i] is i-th mouthful, i=2 ..., n, cm；S is cigarette sectional area, cm²。

Most preferably, the cigarette sectional area is calculated by formula (m), and the formula (m) is：

S=c × c/4/Pi,

Wherein, s is cigarette sectional area, cm²；C is Cigarette circumference, cm；Pi is pi.

Most preferably, the length of intersecting of the tipping paper and cigarette paper is calculated by formula (n), and the formula (n) is：

L5=lp-lf,

Wherein, l5 is the length of intersecting of tipping paper and cigarette paper, cm；Lp is tipping paper width, cm；Lf is filter lengths, cm。

Most preferably, the tipping paper punching leading portion filter lengths are calculated by formula (o), and the formula (o) is：

L4=lf-lt,

Wherein, l4 is that tipping paper punches leading portion filter lengths, cm；Lf is filter lengths, cm；Lt be the outer pore area of tipping paper away from Side, cm.

Step S32 calculates every mouthful of burning relative surplus length by formula (10), (11), and the formula (10) is：

L7 [1]=l1-l6 [1],

Wherein, every mouthful of burning relative surplus length that l7 [1] is the 1st mouthful, cm；L1 is cigarette segment length, cm；L6 [1] is Relative combustion length after 1st mouthful of ventilation, cm；

The formula (11) is：

L7 [i]=l7 [i-1] -58v-l6 [i],

Wherein, every mouthful of burning relative surplus length that l7 [i] is i-th mouthful, i=2 ..., n, cm；L7 [i-1] is (i-1)-th mouthful Every mouthful of burning relative surplus length, i=2 ..., n, cm；V is smoldering velocity, cm/s；L6 [i] is opposite after i-th mouthful of ventilation Burning length, i=1 ..., n, cm.

Step S33 calculates every mouthful of burning relative weight by formula (12), and the formula (12) is：

Rm2 [i]=m1 × l6 [i]/l1,

Wherein, rm2 [i] is every mouthful of burning relative weight, i.e., i-th mouthful burning relative weight, i=1 ..., n, g；M1 is cigarette Silk weight, g；L6 [i] is the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；L1 is cigarette segment length, cm.

It is highly preferred that every mouthful of burning relative weight is calculated separately every mouthful of burning phase of nicotine by formula (p), (q) To weight, every mouthful of burning relative weight of tar,

The formula (p) is：

Rm3 [i]=rm2 [i] × nic,

Wherein, rm3 [i] is every mouthful of burning relative weight of nicotine, i.e. the i-th of nicotine mouthful burning relative weight, i= 1 ..., n, g；Rm2 [i] is i-th mouthful of burning relative weight, i=1 ..., n, g；Nic is the nicotine empirical coefficient of pipe tobacco, setting It is 1.11；

The formula (q) is：

Rm4 [i]=rm2 [i] × tar,

Wherein, rm4 [i] is every mouthful of burning relative weight of tar, i.e. the i-th of tar mouthful burning relative weight, i= 1 ..., n, g；Rm2 [i] is i-th mouthful of burning relative weight, i=1 ..., n, g；Tar is the tar empirical coefficient of pipe tobacco, is set as 1.22。

Preferably, in step S3, the pressure drop matching regression algorithm refers to, every mouthful of the cigarette original burning obtained by S2 Based on length, pressure drop matching is carried out to cigarette pressure drop and tipping paper, forming paper and filter stick system by formula, is matched each time Corresponding i-th mouthful of filter tip ventilation flow rate and combustion cone inlet flow rate are obtained, combustion cone inlet flow rate determines opposite mouth number, every mouthful Burning relative surplus length, every mouthful of burning relative weight (weight including pipe tobacco relative combustion nicotine and tar).

Step S4, using pipe tobacco by mouth release-equation of filtration, according in step S2 in the initial parameter of cigarette and step S3 The relative parameter of cigarette obtains Contents of Main Components value in cigarette smoke to calculate.

Preferably, in step S4, main component includes but is not limited to tar, nicotine in the cigarette smoke.

Preferably, in step S4, the calculating of Contents of Main Components value (i.e. burst size) in the cigarette smoke, including it is following Step：

Step S41, the absolute burst size of whole branch Nicotine from Cigarette is calculated by formula (13), and the formula (13) is：f0tot =Sigma [f0 [i]],

Wherein, f0tot is the absolute burst size of whole branch Nicotine from Cigarette, mg；Simgais an add sign∑,is the array F0 [i] summation；F0 [i] is the absolute burst size of the nicotine of i-th mouthful of cigarette burning, i=1 ..., n, mg.

It is highly preferred that the absolute burst size of the nicotine of i-th mouthful of cigarette burning is calculated by formula (r), the formula (r) For：

F0 [i]=g × qi [i] × (m3+f1 [i-1,1]+f1 [i-2,2]+...+f1 [i, j]) × e^(-a4×l2), i= 1 ..., n,

Wherein, f0 [i] is the absolute burst size of the nicotine of i-th mouthful of cigarette burning, mg；G is empirical coefficient, and setting value is 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M3 is the original calculation weight of nicotine, g；a4 For empirical coefficient, setting value 0.22；E is the truth of a matter of natural logrithm；L2 is every mouthful of burning original length, cm；F1 [i, j] is the The absolute nicotine function that the cigarette of i mouthfuls of burnings is trapped in jth mouth burning zone, mg.

It is further preferred that the absolute nicotine function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone is by public affairs Formula (s) calculates, and the formula (s) is：

F1 [i, j]=g × qi [i] × (m3+f1 [i-1,1]+f1 [i-2,2]+...+f1 [1, j]) × e^{(-a4×l3[j+i-1])}× (1-e^(-a4×l2)),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, f1 [i, j] is the absolute nicotine function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is Empirical coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M3 is nicotine Original calculation weight, g；A4 is empirical coefficient, setting value 0.22；Every mouthful of burning that l3 [j+i-1] is jth+i-1 mouthfuls is original surplus Remaining length, cm；L2 is every mouthful of burning original length, cm；E is the truth of a matter of natural logrithm.

Step S42, the absolute burst size of tar in whole branch cigarette is calculated by formula (14), and the formula (14) is： F02tot=Sigma [f02 [i]],

Wherein, f02tot is the absolute burst size of tar in whole branch cigarette, mg；Simgais an add sign∑,is the array F02 [i] summation；F02 [i] is the absolute burst size of the tar of i-th mouthful of cigarette burning, i=1 ..., n, mg.

It is highly preferred that the absolute burst size of the tar of i-th mouthful of cigarette burning is calculated by formula (t), the formula (t) For：

F02 [i]=g × qi [i] × (m4+f2 [i-1,1]+f2 [i-2,2]+...+f2 [i, j]) × e^(-a4×l2), i= 1 ..., n,

Wherein, f02 [i] is the absolute burst size of the tar of i-th mouthful of cigarette burning, mg；G is empirical coefficient, and setting value is 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M4 is the original calculation weight of tar, g；A4 is Empirical coefficient, setting value 0.22；E is the truth of a matter of natural logrithm；L2 is every mouthful of burning original length, cm；F2 [i, j] is i-th The absolute tar function that the cigarette of mouth burning is trapped in jth mouth burning zone, mg.

It is further preferred that the absolute tar function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone is by public affairs Formula (u) calculates, and the formula (u) is：

F2 [i, j]=g × qi [i] × (m4+f2 [i-1,1]+f2 [i-2,2]+...+f2 [1, j]) × e^{(-a4×l3[j+i-1])}× (1-e^(-a4×l2)),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, f2 [i, j] is the absolute tar function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is Empirical coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M4 is the original of tar Beginning calculated weight, g；A4 is empirical coefficient, setting value 0.22；Every mouthful of burning original residual that l3 [j+i-1] is jth+i-1 mouthfuls Length, cm；L2 is every mouthful of burning original length, cm；E is the truth of a matter of natural logrithm.

Step S43, the opposite burst size of whole branch Nicotine from Cigarette is calculated by formula (15), and the formula (15) is： Rf0tot=Sigma [rf0 [i]] × (1-niceff),

Wherein, rf0tot is the opposite burst size of whole branch Nicotine from Cigarette, mg；Simgais an add sign∑,is the array Rf0 [i] summation；The opposite burst size for the nicotine that rf0 [i] is i-th mouthful, i=1 ..., n, mg；Niceff is the nicotine mistake of filter stick Filter coefficient.

It is highly preferred that the opposite burst size of i-th mouthful of the nicotine is calculated by formula (v), the formula (v) is：

Rf0 [i]=g × qi [i] × (rm3 [i]+rf1 [i-1,1]+rf1 [i-2,2]+...+rf1 [i, j]) × e^(-a4×l6[i]), i=1 ..., n,

Wherein, the opposite burst size for the nicotine that rf0 [i] is i-th mouthful, mg；G is empirical coefficient, setting value 0.297；qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm3 [i] is i-th mouthful of burning relative weight of nicotine, g； Rf1 [i, j] is the opposite nicotine function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone, mg；A4 is empirical coefficient, if Definite value is 0.22；E is the truth of a matter of natural logrithm；L6 [i] is the relative combustion length after i-th mouthful of ventilation, cm.

It is further preferred that the opposite nicotine function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone is by public affairs Formula (w) calculates, and the formula (w) is：

Rf1 [i, j]=g × qi [i] × (rm3 [i]+rf1 [i-1,1]+rf1 [i-2,2]+...+rf1 [1, j]) × e^{(-a4×l7[j+i-1])}×(1-e^{(-a4×l6[i+j-1])}),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, rf1 [i, j] is the opposite nicotine function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；g For empirical coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm3 [i] is Buddhist nun I-th mouthful of burning relative weight of Gu Ding, i=1 ..., n, g；A4 is empirical coefficient, setting value 0.22；L7 [j+i-1] be jth+ I-1 mouthfuls of every mouthful of burning relative surplus length, cm；L6 [j+i-1] is the relative combustion length after jth+i-1 mouthfuls of ventilation, cm；e For the truth of a matter of natural logrithm.

Step S44, the opposite burst size of tar in whole branch cigarette is calculated by formula (16), and the formula (16) is： Rf02tot=Sigma [rf02 [i]] × (1-tareff),

Wherein, rf02tot is the opposite burst size of tar in whole branch cigarette, mg；Simgais an add sign∑,is logarithm Group rf02 [i] summation；The opposite burst size for the tar that rf02 [i] is i-th mouthful, mg；Tareff is the tar filtration coefficient of filter stick.

It is highly preferred that the opposite burst size of i-th mouthful of the tar is calculated by formula (x), the formula (x) is：

Rf02 [i]=g × qi [i] × (rm4 [i]+rf2 [i-1,1]+rf2 [i-2,2]+...+rf2 [i, j]) × e^(-a4×l6[i]), i=1 ..., n,

Wherein, the opposite burst size for the tar that rf02 [i] is i-th mouthful, mg；G is empirical coefficient, setting value 0.297；qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm4 [i] is i-th mouthful of burning relative weight of tar, i= 1 ..., n, g；A4 is empirical coefficient, setting value 0.22；Rf2 [i, j] is that the cigarette of i-th mouthful of burning is cut in jth mouth burning zone The opposite tar function stayed, mg；E is the truth of a matter of natural logrithm；L6 [i] is the relative combustion length after i-th mouthful of ventilation, cm.

It is further preferred that the opposite tar function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone is by public affairs Formula (y) calculates, and the formula (y) is：

Rf2 [i, j]=g × qi [i] × (rm4 [i]+rf2 [i-1,1]+rf2 [i-2,2]+...+rf2 [1, j]) × e^{(-a4×l7[j+i-1])}×(1-e^{(-a4×l6[i+j-1])}),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, rf2 [i, j] is the opposite tar function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；g For empirical coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm4 [i] is coke I-th mouthful of burning relative weight of oil, i=1 ..., n, g；A4 is empirical coefficient, setting value 0.22；E is the bottom of natural logrithm Number；Every mouthful of burning relative surplus length that l7 [j+i-1] is jth+i-1 mouthfuls, cm；After l6 [j+i-1] is jth+i-1 mouthfuls of ventilation Relative combustion length, cm.

Preferably, in step S4, the pipe tobacco refers to by mouth release-equation of filtration, has obtained original yin by step S2 Time, original mouth number, every mouthful of burning original length, every mouthful of burning original weight, every mouthful of burning original residual length are fired, step is passed through Rapid S3 has obtained opposite mouth number, every mouthful of burning relative surplus length, every mouthful of burning relative weight, then auxiliary by the information of S2 and S3 It helps pipe tobacco by mouth release-equation of filtration, can be calculated absolute burst size, tar and the cigarette of tar and nicotine (reset condition) by mouth The opposite burst size of alkali, i.e. f0tot, f02tot, rf0tot, rf02tot.The absolute burst size of tar and the opposite release of tar The common content for determining tar in cigarette smoke of amount, the absolute burst size of nicotine and the opposite burst size of nicotine determine cigarette jointly The content of nicotine in flue gas.

The present invention also provides a kind of computer readable storage mediums, are stored thereon with computer program, described program quilt The step of processor realizes above-mentioned analysis method when executing.

Invention further provides a kind of computer-processing equipment, including processor and computer readable storage medium, The step of processor executes the computer program on the computer readable storage medium, realizes above-mentioned analysis method.

As described above, in a kind of cigarette smoke provided by the invention Contents of Main Components analysis method, having following has Beneficial effect：

(1) in a kind of cigarette smoke provided by the invention Contents of Main Components analysis method, this method be based on the primary Physics law, assist the semi-theoretical algorithm of semiempirical, can apply in computer, calculate formula, auxiliary material and cigarette structure become The dynamic influence to main flume, has certain universality.

(2) in a kind of cigarette smoke provided by the invention Contents of Main Components analysis method, accurate can predict The content of main component in cigarette smoke has easily and efficiently advantage without carrying out trial test early period of a large amount of repeatability.

(3) in a kind of cigarette smoke provided by the invention Contents of Main Components analysis method, only need to be by each parameter Between correlation reasonably adjusted, according to prediction result, suitably, reasonably adjust tobacco leaf formulation, auxiliary material and roll knot Structure parameter can meet practical Cigarette design target, can be fine without being transformed to existing production of cigarettes process equipment Achieve the effect that cost efficiency.

Detailed description of the invention

Fig. 1 is shown as the flow diagram of the analysis method of Contents of Main Components in a kind of cigarette smoke of the invention.

Specific embodiment

The present invention is further explained combined with specific embodiments below, it should be appreciated that these embodiments are merely to illustrate the present invention Rather than it limits the scope of the invention.

Illustrate embodiments of the present invention below by way of specific specific example, those skilled in the art can be by this specification Other advantages and efficacy of the present invention can be easily understood for disclosed content.The present invention can also pass through in addition different specific realities The mode of applying is embodied or practiced, the various details in this specification can also based on different viewpoints and application, without departing from Various modifications or alterations are carried out under spirit of the invention.

It should be clear that in the following example not specifically dated process equipment or device be all made of conventional equipment in the art or Device；All pressure values and range all refer to relative pressure.

In addition, it should also be understood that, one or more method and step mentioned in the present invention does not repel before and after the combination step It can also be inserted into other methods step there may also be other methods step or between these explicitly mentioned steps, unless separately It is described；It should also be understood that the combination connection relationship between one or more equipment/device mentioned in the present invention is not repelled The two equipment/devices specifically mentioned before and after the unit equipment/device there may also be other equipment/device or at these it Between can also be inserted into other equipment/device, unless otherwise indicated.Moreover, unless otherwise indicated, the number of various method steps is only Identify the convenient tool of various method steps, rather than for the arrangement order of limitation various method steps or limits the enforceable model of the present invention It encloses, relativeness is altered or modified, and without material changes in technical content, when being also considered as, the present invention is enforceable Scope.

As shown in Figure 1, the present invention provides a kind of analysis method of Contents of Main Components in cigarette smoke, including following step Suddenly：

Step S1 determines the underlying parameter of cigarette；

Step S2 calculates the original of cigarette according to the underlying parameter of cigarette in step S1 using smoldering time recurrence method Parameter；

In step sl, the underlying parameter of the cigarette includes but is not limited to pipe tobacco type, pipe tobacco weight, cigarette segment length Degree, Cigarette circumference, filter lengths, filter stick pressure drop, tipping paper air permeability, tipping paper pore area overall width, tipping paper width, tipping paper Outer pore area is away from side, forming paper air permeability, air permeance for cigarette paper, the combustion-supporting agent content of cigarette paper, cigarette paper grammes per square metre, nicotine filtering system Number, tar filtration coefficient.

In step s 2, the initial parameter of the cigarette includes but is not limited to original smoldering time, original mouth number, every mouthful of combustion Burn original length, every mouthful of burning original weight, every mouthful of burning original residual length.

In step s3, the relative parameter of the cigarette includes but is not limited to that mouth number, every mouthful of burning relative surplus relatively is long Degree, every mouthful of burning relative weight.

In step s 4, main component includes but is not limited to tar, nicotine in the cigarette smoke.

Embodiment 1

In a computer, according to Fig. 1, main component in cigarette smoke is analyzed, steps are as follows：

1, cigarette laboratory sample is chosen, determines the underlying parameter of cigarette, the specific data of the underlying parameter of cigarette are shown in Table 1.

The underlying parameter of 1 cigarette laboratory sample of table

Cigarette underlying parameter	Numerical value
		Pipe tobacco type	100% cut tobacco
Pipe tobacco weight	0.7g
		Cigarette segment length	5.1cm
Cigarette circumference	2.45cm
		Filter lengths	3cm
Filter stick pressure drop	8cmH2
		Nicotine filtration coefficient	0.4
Tar filtration coefficient	0.54
		Tipping paper air permeability	Electrostatic punches 150CU
Forming paper air permeability	10000CU
		The outer pore area of tipping paper is away from side	1.45cm
Tipping paper width	3cm
		Air permeance for cigarette paper	60CU
The combustion-supporting agent content of cigarette paper	1.8% pure sylvite
		Cigarette paper grammes per square metre	27g

2, using smoldering time recurrence method, the initial parameter of cigarette is calculated according to the underlying parameter of cigarette in step S1, The specific data of the initial parameter of cigarette are shown in Table 2.

The initial parameter of 2 cigarette laboratory sample of table

The initial parameter of cigarette	Numerical value
		Original smoldering time t1	570s
Original mouth number n	7.4 mouth
		Every mouthful of burning original weight m2 (g)	0.023、0.023、0.023、0.023、0.023、0.023、0.023、0.009
Every mouthful of burning original residual length l3 [i] (cm)	5.2、4.6、3.9、3.3、2.6、2.0、1.3、0.8

The specific calculating process of the initial parameter of cigarette is as follows：

Step 1, pipe tobacco type is determined to calculate original smoldering time t1, the formula by formula (1) after 100% cut tobacco (1) it is：T1=73.3 × l1+279 × c-33.1 × a-19.5 × mc+6.86 × ps+11.7 × pe+6.49 × pp+10.4 × Pl, acquisition t1 are 570s.

Step 2, original mouth number n is calculated by formula (2) by original smoldering time t1, the formula (2) is：N= 0.0121 × t1+0.487, obtaining n is 7.4 mouthfuls.

Step 3, original by every mouthful of formula (3), (4), (5) calculating burning by original smoldering time t1 and original mouth number n Length l2, the formula (3) are：V=l1/t1, the formula (4) are：T2=58 × n, the formula (5) are：L2=v × (t1-t2)/n, acquisition l2 are 0.18cm.

Step 4, by every mouthful of burning original length l2, every mouthful of burning original weight m2, the formula are calculated by formula (6) (6) it is：M2=m1 × l2/l1, obtaining m2 is 0.023.Wherein, every mouthful of burning original weight m2 by formula (a), (b) respectively based on The original calculation weight m3 of nicotine, the original calculation weight m4 of tar are calculated, the formula (a) is：M3=m2 × nic, it is described Formula (b) is：M4=m2 × tar, acquisition m3 are 0.026, m4 0.028.

Step 5, by every mouthful of burning original length l2, every mouthful of burning original residual length l3 is calculated by formula (7), (8) [i], the formula (7) are：L3 [1]=l1-l2, the formula (8) are：L3 [i]=l3 [i-1] -58v-l2 is obtained l3 [i] It is 5.2,4.6,3.9,3.3,2.6,2.0,1.3,0.8.

3, using pressure drop match regression algorithm, according in the underlying parameter of cigarette in step S1 and step S2 cigarette it is original Parameter

The relative parameter of cigarette is calculated, the specific data of the relative parameter of cigarette are shown in Table 3.

The relative parameter of 3 cigarette laboratory sample of table

The relative parameter of cigarette	Numerical value
		Opposite mouth number rn	7.7 mouth
Every mouthful of burning relative surplus length l7 [i] (mm)	5.2、4.6、4.0、3.3、2.7、2.0、1.4、0.8
		Every mouthful of burning relative weight rm2 [i] (g)	0.022、0.022、0.022、0.022、0.022、0.022、0.022、0.016
I-th mouthful of burning relative weight rm3 [i] (g) of nicotine	0.024,0.024,0.024,0.024,0.024,0.024,0.024,0.017
		I-th mouthful of burning relative weight rm4 [i] (g) of tar	0.027,0.027,0.027,0.027,0.027,0.027,0.027,0.019

The specific calculating process of the relative parameter of cigarette is as follows：

Step 1, pressure drop-flow function pf [q] of filter stick is calculated by formula (k), the formula (k) is：Pf [q]=k × (17.5-q) × l4 is obtained pf [q]=0.214 (17.5-q).

Step 2, the pressure drop without cigarette under cigarette paper ventilation condition-flow function pre [q], the public affairs are calculated by formula (j) Formula (j) is：Pre [q]=0.143 × (1-E)²×l5×(17.5-q)/E³/ s is obtained pre [q]=0.016 (17.5-q).Its In, cigarette sectional area s is calculated by formula (m), and the formula (m) is：S=c × c/4/Pi, obtaining s is 0.478.

Step 3, pressure drop-flow function pr [q] under ventilation condition containing cigarette paper, the formula (i) are calculated by formula (i) For：

Obtain pr [q]=0.07Tanh [52.7 × Sqrt [7 × 10^-5+2.9×10^-6(17.5-q)]]/Sqrt[7×10^-5 +2.9×10^-6(17.5-q)].Wherein, the voltage drop value pre0 that cigarette section ignores air permeance for cigarette paper is calculated by formula (l), described Formula (l) is：Pre0=2.5 × (1-E)²×l3[i]/E³/ s, obtain pre0 be (3.64,3.17,2.74,2.28,1.83, 1.38,0.93,0.56).

Step 4, tipping paper-forming paper pressure drop-flow function value ptpw [q] is calculated by formula (h), the formula (h) is： Ptpw [q]=(q/c)²×(2.56×10⁸/a2²+3.6×10⁴/a3²), obtaining ptpw [q] is 0.36q²。

Step 5, q value is calculated by the recurrence pressure drop matching equation in formula (g), the formula (g) is：Q=FindRoot [ptpw [q]=pr [q]+pre [q]+pf [q], q->0],

Obtaining q is (2.15,2.09,2.02,1.94,1.85,1.76,1.67,1.58).

Step 6, the flow speed value of cigarette burning entrance is calculated by formula (e), the formula (e) is：Qi [i]=0.9 × (17.5-q [i]), obtaining qi [i] is (13.8,13.87,13.94,14,14.1,14.2,14.25,14.3).

Step 7, by the flow speed value of cigarette burning entrance, it is long that the relative combustion after i-th mouthful of ventilation is calculated by formula (d) Degree, the formula (d) be l6 [i]=qi [i]/17.5 × l2, obtain l6 [i] be (0.168,0.169,0.17,0.171, 0.171,0.173,0.174,0.174).

Step 8, by the relative combustion length after i-th mouthful of ventilation, average relative combustion length, institute are calculated by formula (c) Stating formula (c) is：L6avg=Simga ( l6 [ i ] ) /n,to obtainl6avgis0.17.

Step 9, opposite mouth number rn is calculated by formula (9), the formula (9) is：Rn=l1/ (l6avg+58v) obtains rn It is 7.7.

Step 10, every mouthful of burning relative surplus length l7 [i] is calculated by formula (10), (11), the formula (10) is：l7 [1]=l1-l6 [1], the formula (11) are：L7 [i]=l7 [i-1] -58v-l6 [i], obtain l7 [i] be 5.2,4.6,4.0, 3.3、2.7、2.0、1.4、0.8。

Step S11 calculates i-th mouthful of burning relative weight rm2 [i] by formula (12), and the formula (12) is：Rm2 [i]= M1 × l6 [i]/l1, obtaining rm2 [i] is (0.022,0.022,0.022,0.022,0.022,0.022,0.022,0.016).So Afterwards, i-th mouthful of burning relative weight calculates separately i-th mouthful of nicotine relative weight of burning, tar by formula (p), (q) I-th mouthful of burning relative weight, wherein the formula (p) is：Rm3 [i]=rm2 [i] × nic, the formula (q) are：rm4[i] =rm2 [i] × tar, obtaining rm3 [i] is (0.024,0.024,0.024,0.024,0.024,0.024,0.024,0.017), Obtaining rm4 [i] is (0.027,0.027,0.027,0.027,0.027,0.027,0.027,0.019).

4, using pipe tobacco by mouth release-equation of filtration, according to cigarette in the initial parameter of cigarette in step S2 and step S3 Relative parameter, obtain Contents of Main Components value in cigarette smoke to calculate, specific data are shown in Table 4.

The flue gas Contents of Main Components of 4 cigarette laboratory sample of table

The main component of cigarette smoke	Numerical value
		The absolute burst size f0tot of original nicotine	1.7mg
Original nicotine experiment value	1.6mg
		The relative deviation of original nicotine	6.25%
The opposite burst size rf0tot of nicotine	0.92mg
		Nicotine experiment value	0.97mg
The relative deviation of nicotine	5.15%
		The absolute burst size f02tot of original tar	23.7mg
Original tar experiments value	23.8mg
		The relative deviation of original tar	0.4%
The opposite burst size rf02tot of tar	13mg
		Tar experiments value	12mg
The relative deviation of tar	8.33%

The specific calculating process of Contents of Main Components value is as follows in cigarette smoke：

Step 1, the absolute nicotine function being trapped by the cigarette that formula (s) calculates i-th mouthful of burning in jth mouth burning zone, The formula (s) is：F1 [i, j]=g × qi [i] × (m3+f1 [i-1,1]+f1 [i-2,2]+...+f1 [1, j]) × e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2)), i<J=1 ..., n；I+j-1=1 ..., n are obtained f1 [1,1]=0.0041, f1 [1,2] =0.0038 ... f [6,1]=0.0018.

Step 2, the absolute burst size f0 [i] of the nicotine of i-th mouthful of cigarette burning, the formula (r) are calculated by formula (r) For：

Obtain f0 [i]=(0.16,0.18,0.2,0.22,0.25,0.27,0.29,0.09).

Step 3, the absolute burst size of nicotine (under reset condition) in whole branch cigarette, the formula are calculated by formula (13) (13) it is：F0tot=Sigma [f0 [i]], acquisition f0tot are 1.7mg.

Step 4, the absolute tar function being trapped by the cigarette that formula (u) calculates i-th mouthful of burning in jth mouth burning zone, The formula (u) is：F2 [i, j]=g × qi [i] × (m4+f2 [i-1,1]+f2 [i-2,2]+...+f2 [1, j]) × e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2)), i<J=1 ..., n；I+j-1=1 ..., n are obtained f2 [1,1]=0.285, f2 [1,2] =0.027 ..., f2 [6,1]=0.035.

Step 5, the absolute burst size of the tar of i-th mouthful of cigarette burning is calculated by formula (t), the formula (t) is：

Obtain f02 [i]=(1.24,1.74,2.4,3.3,4.3,5.7,4.3,0.73).

Step 6, the absolute burst size of tar (under reset condition) in whole branch cigarette, the formula are calculated by formula (14) (14) it is：F02tot=Sigma [f02 [i]], acquisition f02tot are 23.7mg.

Step 7, the opposite nicotine function being trapped by the cigarette that formula (w) calculates i-th mouthful of burning in jth mouth burning zone, The formula (w) is：

i<J=1 ..., n；I+j-1=1 ..., n obtain rf1 [1,1]=0.004, rf1 [1,2]=0.0038 ..., rf1 [6,1]=0.0018.

Step 8, the opposite burst size of i-th mouthful of nicotine is calculated by formula (v), the formula (v) is：Rf0 [i]=g × qi[i]×(rm3[i]+rf1[i-1,1]+rf1[i-2,2]+…+rf1[i,j])×e^(-a4×l6[i]), i=1 ..., n are obtained Rf0 [i]=(0.08,0.095,0.11,0.12,0.13,0.15,0.14,0.09).

Step 9, the opposite burst size rf0tot of whole branch Nicotine from Cigarette is calculated by formula (15), the formula (15) is： Rf0tot=Sigma [rf0 [i]] × (1-niceff), acquisition rf0tot are 0.92mg.

Step S10, the opposite tar letter being trapped by the cigarette that formula (y) calculates i-th mouthful of burning in jth mouth burning zone Number, the formula (y) are：

i<J=1 ..., n；I+j-1=1 ..., n obtain rf2 [1,1]=0.028, rf2 [1,2]=0.027 ..., rf2 [6,1]=0.035.

Step S11, the opposite burst size of i-th mouthful of tar is calculated by formula (x), and the formula (x) is：

Obtain rf02 [i]=(0.63,0.9,1.3,1.7,2.3,3.1,2.3,0.73).

Step 12, the opposite burst size of tar in whole branch cigarette is calculated by formula (16), the formula (16) is： Rf02tot=Sigma [rf02 [i]] × (1-tareff), acquisition rf02tot are 13mg.

In above-mentioned 2 and 4, is provided according in standard GB/T/T 19609-2004, tar in cigarette smoke is examined It surveys, obtains original tar experiments value and tar experiments value；It is provided according in standard GB/T/T 23355-2009, to cigarette cigarette Nicotine is detected in gas, obtains original nicotine experiment value and nicotine experiment value.Specific detection device uses standard GB/T/T Regulation sets experimental facilities parameter in 16450-2004：Pumping volume is 35mL, 60s/ mouthfuls of suction frequencies, suction time For 2s.The relative deviation of Contents of Main Components is calculated by formula (z) in cigarette smoke, and the formula (z) is w=Abs ((w1-w2)/w2), in formula, w is relative deviation, and w1 is the predicted value of a certain ingredient in flue gas, and w2 is tie element in flue gas Experiment value.

As shown in Table 4, provide through the invention method analysis cigarette smoke in tar and nicotine content, experiment value and Predicted value is very close, and the relative deviation of the content of tar and nicotine is 8.33% and 5.15%, obtained cigarette cigarette respectively The experiment value relative deviation of the predicted value of Contents of Main Components and corresponding ingredient fully meets Cigarette design and wants less than 10% in gas It asks.

In conclusion in a kind of cigarette smoke provided by the invention Contents of Main Components analysis method, can apply to count In calculation machine, calculate formula, auxiliary material and cigarette structure change influence to main flume, can be in accurate prediction cigarette smoke The content of main component, without being transformed to existing production of cigarettes process equipment, the early period without carrying out a large amount of repeatability is pre- Test has easily and efficiently advantage.So the present invention effectively overcomes various shortcoming in the prior art and has height and produce Industry utility value.

The above-described embodiments merely illustrate the principles and effects of the present invention, and is not intended to limit the present invention.It is any ripe The personage for knowing this technology all without departing from the spirit and scope of the present invention, carries out modifications and changes to above-described embodiment.Cause This, institute is complete without departing from the spirit and technical ideas disclosed in the present invention by those of ordinary skill in the art such as At all equivalent modifications or change, should be covered by the claims of the present invention.

Claims

1. the analysis method of Contents of Main Components, includes the following steps in a kind of cigarette smoke：

Step S1 determines the underlying parameter of cigarette；

Step S2 calculates the initial parameter of cigarette according to the underlying parameter of cigarette in step S1 using smoldering time recurrence method；

Step S3 matches regression algorithm using pressure drop, according to the original of cigarette in the underlying parameter of cigarette in step S1 and step S2 Beginning parameter calculates the relative parameter of cigarette；

Step S4, using pipe tobacco by mouth release-equation of filtration, according to cigarette in the initial parameter of cigarette in step S2 and step S3 Relative parameter, thus calculate obtain cigarette smoke in Contents of Main Components value.

2. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 1, which is characterized in that described The underlying parameter of cigarette includes but is not limited to pipe tobacco type, pipe tobacco weight, cigarette segment length, Cigarette circumference, filter lengths, filter stick Pressure drop, tipping paper air permeability, tipping paper pore area overall width, tipping paper width, the outer pore area of tipping paper away from side, forming paper air permeability, The combustion-supporting agent content of air permeance for cigarette paper, cigarette paper, cigarette paper grammes per square metre, nicotine filtration coefficient, tar filtration coefficient.

3. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 1, which is characterized in that step In S2, the initial parameter of the cigarette includes but is not limited to original smoldering time, original mouth number, every mouthful of burning original length, every Mouth burning original weight, every mouthful of burning original residual length；The calculating of the initial parameter of the cigarette, includes the following steps：Step Rapid S21 after determining pipe tobacco type, calculates original smoldering time by formula (1), the formula (1) is：

Wherein, t1 is original smoldering time, s；L1 is cigarette segment length, cm；C is Cigarette circumference, cm；A is cigarette paper combustion adjuvant Content, %；Mc is cigarette paper grammes per square metre, g/m²；Ps is stem ratio, %；Pe is swollen silk ratio, %；Pp is thin slice ratio, %； Pl is cut tobacco ratio, %；

N=0.0121 × t1+0.487, wherein n is original mouth number, mouth；T1 is original smoldering time, s；

Step S23 calculates every mouthful of burning original length by formula (3), (4), (5) by original smoldering time and original mouth number, The formula (3) is：V=l1/t1, wherein v is smoldering velocity, cm/s；L1 is cigarette segment length, cm；When t1 is original glows Between, s；

The formula (4) is：T2=58 × n, wherein t2 is practical smoldering time, s；N is original mouth number, mouth；

The formula (5) is：L2=v × (t1-t2)/n, wherein l2 is every mouthful of burning original length, cm；V is smoldering velocity, cm/s；T1 is original smoldering time, s；T2 is practical smoldering time, s；N is original mouth number, mouth；

Step S24 calculates every mouthful of burning original weight by formula (6), the formula (6) is by every mouthful of burning original length： M2=m1 × l2/l1, wherein m2 is every mouthful of burning original weight, g；M1 is pipe tobacco weight, g；L2 is every mouthful of original length of burning Degree, cm；L1 is cigarette segment length, cm；

Step S25 calculates every mouthful of burning original residual length, the public affairs by formula (7), (8) by every mouthful of burning original length Formula (7) is：L3 [1]=l1-l2, wherein every mouthful of burning original residual length that l3 [1] is the 1st mouthful, cm；L1 is cigarette segment length Degree, cm；L2 is every mouthful of burning original length, cm；

The formula (8) is：L3 [i]=l3 [i-1] -58v-l2, wherein l3 [i] is that i-th mouthful of every mouthful of burning original residual is long Degree, i=2 ..., n, cm；Every mouthful of burning original residual length that l3 [i-1] is (i-1)-th mouthful, i=2 ..., n, cm；V is speed of glowing Degree, cm/s；L2 is every mouthful of burning original length, cm.

4. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 3, which is characterized in that step In S24, every mouthful of burning original weight is calculated separately original calculation weight, the original of tar of nicotine by formula (a), (b) Beginning calculated weight,

The formula (a) is：M3=m2 × nic, wherein m3 is the original calculation weight of nicotine, g；M2 is that every mouthful of burning is former Starting weight amount, g；Nic is the nicotine empirical coefficient of pipe tobacco, is set as 1.11.The nicotine empirical coefficient nic of the pipe tobacco is to adjust Save the coefficient of different leaf group nicotine deliveries；

The formula (b) is：M4=m2 × tar, wherein m4 is the original calculation weight of tar, g；M2 is that every mouthful of burning is original Weight, g；Tar is the tar empirical coefficient of pipe tobacco, is set as 1.22.The tar empirical coefficient tar is to adjust different leaf groups cokes The coefficient of oily burst size.

5. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 1, which is characterized in that step In S3, the relative parameter of the cigarette includes but is not limited to that mouth number, every mouthful of burning relative surplus length, every mouthful of burning relatively are opposite Weight；The calculating of the relative parameter of the cigarette, includes the following steps：

Step S31 calculates opposite mouth number by formula (9), and the formula (9) is：Rn=l1/ (l6avg+58v),

Wherein, rn is opposite mouth number, mouth；L1 is cigarette segment length, cm；L6avg is average relative combustion length, cm；V is to glow Speed, cm/s；

Step S32 calculates every mouthful of burning relative surplus length by formula (10), (11), and the formula (10) is：L7 [1]=l1- L6 [1],

Wherein, every mouthful of burning relative surplus length that l7 [1] is the 1st mouthful, cm；L1 is cigarette segment length, cm；L6 [1] is the 1st mouthful Relative combustion length after ventilation, cm；

The formula (11) is：L7 [i]=l7 [i-1] -58v-l6 [i],

Wherein, every mouthful of burning relative surplus length that l7 [i] is i-th mouthful, i=2 ..., n, cm；L7 [i-1] is (i-1)-th mouthful every Mouth burning relative surplus length, i=2 ..., n, cm；V is smoldering velocity, cm/s；L6 [i] is the relative combustion after i-th mouthful of ventilation Length, i=1 ..., n, cm；

Step S33 calculates every mouthful of burning relative weight by formula (12), and the formula (12) is：Rm2 [i]=m1 × l6 [i]/ L1, wherein rm2 [i] is every mouthful of burning relative weight, i.e., i-th mouthful burning relative weight, i=1 ..., n, g；M1 is pipe tobacco weight Amount, g；L6 [i] is the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；L1 is cigarette segment length, cm.

6. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 5, which is characterized in that step In S31,

The average relative combustion length is calculated by formula (c), and the formula (c) is：L6avg=Simga ( l6 [ i ] ) /n,

Wherein, l6avg is average relative combustion length, cm；Simgais an add sign∑,is the arrayl6 [ i ] summation；l6[i] For the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；N is original mouth number, mouth；

Relative combustion length after i-th mouthful of ventilation is calculated by formula (d), and the formula (d) is：L6 [i]=qi [i]/ 17.5 × l2,

Wherein, l6 [i] is the relative combustion length after i-th mouthful of ventilation, i=1 ..., n, cm；Qi [i] is the cigarette combustion at i-th mouthful Burn inlet flow rate qi value, i=1 ..., n, ml/s；L2 is every mouthful of burning original length, cm；

The cigarette burning inlet flow rate qi value at i-th mouthful is calculated by formula (e), and the formula (e) is：Qi [i]=0.9 × (17.5-q [i]),

Wherein, qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, i=1 ..., n, ml/s；Q [i] is at i-th mouthful The ventilation flow rate q value of filter tip when length l3 [i], i=1 ..., n, ml/s；

Flow velocity identity in i-th mouthful of the cigarette burning inlet flow rate coincidence formula (f), the formula (f) are：qi[i]+ Qc [i]+q [i]=17.5,

Wherein, qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, i=1 ..., n, ml/s；Qc [i] is at i-th mouthful Cross the flow velocity qc value of cigarette paper, i=1 ..., n, ml/s；Q [i] is the ventilation flow rate q of the filter tip in i-th mouthful of length l3 [i] Value, i=1 ..., n, ml/s；

The flow speed of the filter tip is calculated by the recurrence pressure drop matching equation in formula (g), and the formula (g) is：Q= FindRoot [ptpw [q]=pr [q]+pre [q]+pf [q], q->0],

Wherein, q is the flow speed by filter tip, ml/s；FindRoot is to pass through initial value (q->0) numerical value of equation is searched for Solution, ml/s；q->0 refers to that q is intended to 0；Ptpw [q] is that tipping paper and forming paper combine the function of pressure drop at flow velocity q, cmH₂O；Pr [q] is that cigarette section considers function of pressure drop of the air permeance for cigarette paper at flow velocity q, cmH₂O；Pre [q] is flow velocity when being q Transposition section cigarette function of pressure drop, cmH₂O；The filter stick function of pressure drop that pf [q] is flow velocity when being q, cmH₂O。

7. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 6, which is characterized in that described The function of pressure drop of tipping paper and forming paper combination at flow velocity q calculates tipping paper-forming paper pressure drop-flow function by formula (h), The cigarette section considers that air permeance for cigarette paper is calculated under ventilation condition containing cigarette paper in the function of pressure drop under flow velocity q by formula (i) Pressure drop-flow function, the flow velocity be q when transposition section cigarette function of pressure drop by formula (j) calculate without cigarette paper air duct slats Pressure drop-flow function of cigarette under part, the filter stick function of pressure drop when flow velocity is q are calculated pressure drop-stream of filter stick by formula (k) Flow function,

The formula (h) is：Ptpw [q]=(q/c)²×(2.56×10⁸/a2²+3.6×10⁴/a3²),

Wherein, ptpw [q] is that tipping paper and forming paper combine the function of pressure drop at flow velocity q, cmH₂O；Q is the stream by filter tip Measure speed, ml/s；C is Cigarette circumference, cm；A2 is forming paper air permeability, CU；A3 is tipping paper air permeability, CU；

The formula (i) is：

Pr [q]=pre0 × Tanh (Sqrt (a1 × c × l3 [i] × pre0 × (7 × 10^-5+2.9×10^-6×(17.5-q)))/ Sqrt(a1×c×l3[i]×pre0×(7×10^-5+2.9×10^-6× (17.5-q))),

Wherein, pr [q] is that cigarette section considers function of pressure drop of the air permeance for cigarette paper at flow velocity q, cmH₂O；Pre0 is that cigarette section is neglected The slightly voltage drop value of air permeance for cigarette paper, cmH₂O；Tanh is hyperbolic tangent function；Sqrt is radical sign function；A1 is ventilative for cigarette paper Degree, CU；C is Cigarette circumference, cm；Every mouthful of burning original residual length that l3 [i] is i-th mouthful, i=2 ..., n, cm；Q is to pass through The flow speed of filter tip, ml/s；

The formula (j) is：Pre [q]=0.143 × (1-E)²×l5×(17.5-q)/E³/ s,

Wherein, the transposition section cigarette function of pressure drop that pre [q] is flow velocity when being q, cmH₂O；E is pipe tobacco porosity, is set as 0.76； L5 is the length of intersecting of tipping paper and cigarette paper, cm；S is cigarette sectional area, cm²；

The formula (k) is：Pf [q]=k × (17.5-q) × l4,

Wherein, the filter stick function of pressure drop that pf [q] is flow velocity when being q, cmH₂O；K is filter stick pressure-drop coefficient, cmH2O*s/cm³；Q is Pass through the flow speed of filter tip, ml/s；L4 is that tipping paper punches leading portion filter lengths, cm.The filter stick pressure-drop coefficient is basis Filter stick pressure drop and filter lengths, which calculate, to be obtained, and each filter stick unit all has corresponding filter stick pressure drop and filter lengths.

8. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 7, which is characterized in that

The voltage drop value that the cigarette section ignores air permeance for cigarette paper is calculated by formula (l), and the formula (l) is：Pre0=2.5 × (1-E)²×l3[i]/E³/ s,

Wherein, pre0 is the voltage drop value that cigarette section ignores air permeance for cigarette paper, cmH₂O；E is pipe tobacco porosity, is set as 0.76； Every mouthful of burning original residual length that l3 [i] is i-th mouthful, i=2 ..., n, cm；S is cigarette sectional area, cm²；

The cigarette sectional area is calculated by formula (m), and the formula (m) is：S=c × c/4/Pi,

Wherein, s is cigarette sectional area, cm²；C is Cigarette circumference, cm；Pi is pi；

The length of intersecting of the tipping paper and cigarette paper is calculated by formula (n), and the formula (n) is：L5=lp-lf,

Wherein, l5 is the length of intersecting of tipping paper and cigarette paper, cm；Lp is tipping paper width, cm；Lf is filter lengths, cm；

The tipping paper punching leading portion filter lengths are calculated by formula (o), and the formula (o) is：L4=lf-lt,

Wherein, l4 is that tipping paper punches leading portion filter lengths, cm；Lf is filter lengths, cm；Lt is the outer pore area of tipping paper away from side, cm。

9. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 5, which is characterized in that step In S33, every mouthful of burning relative weight is calculated separately every mouthful of burning relative weight, tar of nicotine by formula (p), (q) Every mouthful of burning relative weight,

The formula (p) is：Rm3 [i]=rm2 [i] × nic,

The formula (q) is：Rm4 [i]=rm2 [i] × tar,

Wherein, rm4 [i] is every mouthful of burning relative weight of tar, i.e. the i-th of tar mouthful burning relative weight, i=1 ..., n, g；Rm2 [i] is i-th mouthful of burning relative weight, i=1 ..., n, g；Tar is the tar empirical coefficient of pipe tobacco, is set as 1.22.

10. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 1, which is characterized in that step In rapid S4, main component includes but is not limited to tar, nicotine in the cigarette smoke；Contents of Main Components in the cigarette smoke Value is the calculating of burst size, is included the following steps：

Step S41, the absolute burst size of whole branch Nicotine from Cigarette is calculated by formula (13), and the formula (13) is：F0tot= Sigma [f0 [i]],

Wherein, f0tot is the absolute burst size of whole branch Nicotine from Cigarette, mg；Simgais an add sign∑,is the arrayf0 [ i ] Summation；F0 [i] is the absolute burst size of the nicotine of i-th mouthful of cigarette burning, i=1 ..., n, mg；

Step S42, the absolute burst size of tar in whole branch cigarette is calculated by formula (14), and the formula (14) is：F02tot= Sigma [f02 [i]],

Wherein, f02tot is the absolute burst size of tar in whole branch cigarette, mg；Simgais an add sign∑,is the arrayf02 [i] summation；F02 [i] is the absolute burst size of the tar of i-th mouthful of cigarette burning, i=1 ..., n, mg；

Step S43, the opposite burst size of whole branch Nicotine from Cigarette is calculated by formula (15), and the formula (15) is：Rf0tot= Sigma [rf0 [i]] × (1-niceff),

Wherein, rf0tot is the opposite burst size of whole branch Nicotine from Cigarette, mg；Simgais an add sign∑,is the arrayrf0 [i] summation；The opposite burst size for the nicotine that rf0 [i] is i-th mouthful, i=1 ..., n, mg；Niceff is that the nicotine of filter stick filters system Number；

Step S44, the opposite burst size of tar in whole branch cigarette is calculated by formula (16), and the formula (16) is：Rf02tot= Sigma [rf02 [i]] × (1-tareff),

Wherein, rf02tot is the opposite burst size of tar in whole branch cigarette, mg；Simgais an add sign∑,is the array Rf02 [i] summation；The opposite burst size for the tar that rf02 [i] is i-th mouthful, mg；Tareff is the tar filtration coefficient of filter stick.

11. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 10, which is characterized in that step In rapid S41, the absolute burst size of the nicotine of i-th mouthful of cigarette burning is calculated by formula (r), and the formula (r) is：

F0 [i]=g × qi [i] × (m3+f1 [i-1,1]+f1 [i-2,2]+...+f1 [i, j]) × e^(-a4×l2), i=1 ..., n,

Wherein, f0 [i] is the absolute burst size of the nicotine of i-th mouthful of cigarette burning, mg；G is empirical coefficient, setting value 0.297； Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M3 is the original calculation weight of nicotine, g；A4 is experience Coefficient, setting value 0.22；E is the truth of a matter of natural logrithm；L2 is every mouthful of burning original length, cm；F1 [i, j] is i-th mouthful of combustion The absolute nicotine function that the cigarette of burning is trapped in jth mouth burning zone, mg；

The cigarette of i-th mouthful of burning is calculated in the absolute nicotine function that jth mouth burning zone is trapped by formula (s), the public affairs Formula (s) is：

F1 [i, j]=g × qi [i] × (m3+f1 [i-1,1]+f1 [i-2,2]+...+f1 [1, j]) × e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2)),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, f1 [i, j] is the absolute nicotine function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is experience Coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M3 is the original of nicotine Calculated weight, g；A4 is empirical coefficient, setting value 0.22；L3 [j+i-1] is that jth+i-1 mouthfuls of every mouthful of burning original residual is long Degree, cm；L2 is every mouthful of burning original length, cm；E is the truth of a matter of natural logrithm.

12. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 10, which is characterized in that step In rapid S42, the absolute burst size of the tar of i-th mouthful of cigarette burning is calculated by formula (t), and the formula (t) is：f02[i] =g × qi [i] × (m4+f2 [i-1,1]+f2 [i-2,2]+...+f2 [i, j]) × e^(-a4×l2), i=1 ..., n,

Wherein, f02 [i] is the absolute burst size of the tar of i-th mouthful of cigarette burning, mg；G is empirical coefficient, and setting value is 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M4 is the original calculation weight of tar, g；A4 is Empirical coefficient, setting value 0.22；E is the truth of a matter of natural logrithm；L2 is every mouthful of burning original length, cm；F2 [i, j] is i-th The absolute tar function that the cigarette of mouth burning is trapped in jth mouth burning zone, mg；

The cigarette of i-th mouthful of burning is calculated in the absolute tar function that jth mouth burning zone is trapped by formula (u), the public affairs Formula (u) is：

F2 [i, j]=g × qi [i] × (m4+f2 [i-1,1]+f2 [i-2,2]+...+f2 [1, j]) × e^{(-a4×l3[j+i-1])}×(1-e^(-a4×l2)),

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, f2 [i, j] is the absolute tar function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is experience Coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；M4 is the original meter of tar Calculate weight, g；A4 is empirical coefficient, setting value 0.22；L3 [j+i-1] is that jth+i-1 mouthfuls of every mouthful of burning original residual is long Degree, cm；L2 is every mouthful of burning original length, cm；E is the truth of a matter of natural logrithm.

13. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 10, which is characterized in that step In rapid S43, the opposite burst size of i-th mouthful of the nicotine is calculated by formula (v), and the formula (v) is：

Rf0 [i]=g × qi [i] × (rm3 [i]+rf1 [i-1,1]+rf1 [i-2,2]+...+rf1 [i, j]) × e^(-a4×l6[i]), i =1 ..., n,

Wherein, the opposite burst size for the nicotine that rf0 [i] is i-th mouthful, mg；G is empirical coefficient, setting value 0.297；Qi [i] is In i-th mouthful of cigarette burning inlet flow rate qi value, ml/s；Rm3 [i] is i-th mouthful of burning relative weight of nicotine, g；rf1[i, J] it is the opposite nicotine function that the cigarette of i-th mouthful of burning is trapped in jth mouth burning zone, mg；A4 is empirical coefficient, and setting value is 0.22；E is the truth of a matter of natural logrithm；L6 [i] is the relative combustion length after i-th mouthful of ventilation, cm；

The cigarette of i-th mouthful of burning is calculated in the opposite nicotine function that jth mouth burning zone is trapped by formula (w), the public affairs Formula (w) is：

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, rf1 [i, j] is the opposite nicotine function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is warp Test coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm3 [i] is nicotine I-th mouthful of burning relative weight, i=1 ..., n, g；A4 is empirical coefficient, setting value 0.22；L7 [j+i-1] is jth+i-1 Every mouthful of burning relative surplus length of mouth, cm；L6 [j+i-1] is the relative combustion length after jth+i-1 mouthfuls of ventilation, cm；E is certainly The truth of a matter of right logarithm.

14. the analysis method of Contents of Main Components in a kind of cigarette smoke according to claim 10, which is characterized in that step In rapid S44, the opposite burst size of i-th mouthful of the tar is calculated by formula (x), and the formula (x) is：

Rf02 [i]=g × qi [i] × (rm4 [i]+rf2 [i-1,1]+rf2 [i-2,2]+...+rf2 [i, j]) × e^(-a4×l6[i]), i =1 ..., n,

Wherein, the opposite burst size for the tar that rf02 [i] is i-th mouthful, mg；G is empirical coefficient, setting value 0.297；qi[i] For in i-th mouthful of cigarette burning inlet flow rate qi value, ml/s；I-th mouthful burning relative weight of the rm4 [i] for tar, i=1 ..., N, g；A4 is empirical coefficient, setting value 0.22；Rf2 [i, j] is what the cigarette of i-th mouthful of burning was trapped in jth mouth burning zone Opposite tar function, mg；E is the truth of a matter of natural logrithm；L6 [i] is the relative combustion length after i-th mouthful of ventilation, cm；

The cigarette of i-th mouthful of burning is calculated in the opposite tar function that jth mouth burning zone is trapped by formula (y), the public affairs Formula (y) is：

i<J=1 ..., n；I+j-1=1 ..., n,

Wherein, rf2 [i, j] is the opposite tar function that i-th mouthful of cigarette to burn is trapped in jth mouth burning zone, mg；G is warp Test coefficient, setting value 0.297；Qi [i] is the cigarette burning inlet flow rate qi value at i-th mouthful, ml/s；Rm4 [i] is tar I-th mouthful of burning relative weight, i=1 ..., n, g；A4 is empirical coefficient, setting value 0.22；E is the truth of a matter of natural logrithm；l7 Every mouthful of burning relative surplus length that [j+i-1] is jth+i-1 mouthfuls, cm；L6 [j+i-1] is opposite after jth+i-1 mouthfuls of ventilation Burning length, cm.

15. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is by processor The step of claim 1-14 any analysis method is realized when execution.

16. a kind of computer-processing equipment, special including computer readable storage medium described in processor and claim 15 Sign is that the processor executes the computer program on the computer readable storage medium, realizes that claim 1-14 appoints The step of one analysis method.