CN107490635A - 茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 - Google Patents
茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 Download PDFInfo
- Publication number
- CN107490635A CN107490635A CN201710675588.4A CN201710675588A CN107490635A CN 107490635 A CN107490635 A CN 107490635A CN 201710675588 A CN201710675588 A CN 201710675588A CN 107490635 A CN107490635 A CN 107490635A
- Authority
- CN
- China
- Prior art keywords
- concentration
- regression equation
- equation
- coefficient
- tea
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- A—HUMAN NECESSITIES
- A23—FOODS OR FOODSTUFFS; TREATMENT THEREOF, NOT COVERED BY OTHER CLASSES
- A23F—COFFEE; TEA; THEIR SUBSTITUTES; MANUFACTURE, PREPARATION, OR INFUSION THEREOF
- A23F3/00—Tea; Tea substitutes; Preparations thereof
- A23F3/16—Tea extraction; Tea extracts; Treating tea extract; Making instant tea
- A23F3/18—Extraction of water soluble tea constituents
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/251—Colorimeters; Construction thereof
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/75—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
- G01N21/77—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
- G01N21/78—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator producing a change of colour
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Plasma & Fusion (AREA)
- Food Science & Technology (AREA)
- Polymers & Plastics (AREA)
- Non-Alcoholic Beverages (AREA)
- Tea And Coffee (AREA)
Abstract
本发明提供了一种茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法:对茶样进行多次浸提,对每次浸提得到的茶汤中含有的7种滋味特征成分进行定量分析;针对7种滋味特征成分,分别建立相应的浓度变化线性方程,从线性方程中提取修正函数,建立相应的多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型。本发明对于随机抽取茶叶样本,可通过定量检测样品第1泡茶汤中各滋味特征成分浓度,预测第x(x为大于1的自然数)泡茶汤中各滋味特征成分浓度。该方法具有成本低、效率高和预测准确等特点,其平均预测准确率高达92.0%,进一步可预测随机样本的耐泡特性,同时对于原叶茶多次冲泡所得茶汤的味觉可视化分析等方面具有重要应用价值和指导意义。
Description
技术领域
本发明涉及一种茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法。
背景技术
茶是我国的传统饮料,适合采用冲泡的方式饮用。传统的饮茶习惯是多次反复冲泡,随着冲泡次数的增加,茶汤中滋味特征成分的浓度和比例都在不断变化。而茶汤的滋味味型和强度是由几种主要特征成分决定的,例如,绿茶的滋味特征成分主要是呈苦味的咖啡碱和呈涩味的酯型儿茶素。
市场上茶叶品类繁多。受品种、采摘嫩度、加工工艺和外形的影响,茶叶的耐冲泡特性(可冲泡次数)存在较大差异。市场上某些茶叶产品在宣传中往往以耐冲泡作为其产品特点,但是茶叶到底能泡几泡,众说纷纭,缺乏科学的判定依据。现有的国家和行业标准,均未对茶叶的耐泡性作出规范性要求,因此,在判定茶叶的耐泡特性时尚无标准可依。
有研究表明,西湖龙井茶在多次冲泡时,含量较高的茶多酚类和氨基酸类化合物的浸出具有各自不同的特点,且茶汤中茶多酚与氨基酸浓度的比值并非固定不变,而是随冲泡次数的增加不断变化(邵晓林等,2006),这为基于滋味特征成分来预测茶汤的呈味特性提供了可能性。
本发明依据实测数据,通过分析变化规律,建立数学模型,实现了随机抽取茶叶样本任意泡次茶汤滋味特征成分浓度的预测并在耐泡特性预判和原叶茶饮料浸提工艺的质量控制等方面具有重要应用价值和指导意义。
发明内容
针对现有技术存在的问题,本发明的目的在于提供一种茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法。
本发明采用的技术方案是:
茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法,所述方法包括以下步骤:
(1)对茶叶样品A1进行多次浸提,收集每次浸提得到的茶汤、定容并标记浸提次数;
(2)对步骤⑴所得不同浸提次数的茶汤进行定量分析,分析其中含有的7种滋味特征成分浓度;
(3)针对步骤⑵中的7种滋味特征成分,分别建立特征成分浓度变化与茶汤泡数之间的回归方程,
(4)选取与步骤(1)同类的多种茶叶样品A2~Am,重复步骤(1)~(3),建立茶叶样品A2~Am的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程;m为3~10的自然数;
(5)将m种茶叶样品A1、A2、~Am的各自7种滋味成分的浓度变化与茶汤泡数之间的回归方程的系数与第一泡茶汤的对应滋味成分的测量值之间建立修正函数,建立相应的多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型。
所述预测模型包括:
7种特征成分浓度变化与茶汤泡数之间的回归方程共七个,所述回归方程为指数回归方程、线性回归方程或幂函数回归方程,每个回归方程各有2个待定系数,每个待定系数根据该系数的修正函数确定,修正函数是每种特征成分的回归方程的系数与该特征成分的第一泡茶汤的测量值之间的回归方程,每个待定系数对应一个修正函数,因此7个回归方程共有14个待定系数和14个对应的修正函数。所述修正函数为指数函数、对数函数、线性函数或幂函数函数。
对于待测未知茶叶样品,按照步骤(1)的浸提方法对待测茶叶样品进行一次浸提,得到的第一泡茶汤按照步骤(2)进行定量分析,得到第一泡茶汤中含有的7种滋味特征成分浓度的测量值,将第一泡茶汤的测量值分别代入修正函数,获得各特征成分的回归方程的系数,再将系数分别代入7个回归方程,即得到该样品的7种特征成分浓度变化与茶汤泡数之间的回归方程,根据回归方程预测第二泡茶汤及以后泡数茶汤的7种特征成分的浓度,并可根据Dot值判定茶叶可冲泡次数,所述Dot值=第x泡下该滋味特征成分的浓度(mg/ml)/该滋味特征成分的味觉阈值,Dot值大于1表示人能感知到该成分的滋味,数值越大,味觉强度越大,小于1表示该化合物不呈味。
所述的步骤⑴中,茶叶样品A1可以为各类茶叶,如绿茶、乌龙茶等,优选为绿茶,优选为炒青绿茶。
所述步骤⑴中,茶叶样品A1的浸提步骤优选为:
(a)取3g茶叶样品A1加入100℃蒸馏水150ml后摇匀,浸提5min后迅速用脱脂棉过滤,滤液冷却后用蒸馏水定容至200ml,标记为第1泡茶汤,
(b)过滤时得到的叶底重复步骤(a)进行二次浸提,得到第2泡茶汤,多次重复步骤(a),分别得到第3至第n泡茶汤,n为大于3的自然数,优选n为4~10的自然数,更优选为4~8的自然数。
所述的步骤⑵中,所述的7种滋味特征成分为咖啡碱(Caf)、没食子酸(GA)、儿茶素(C)、表没食子儿茶素没食子酸酯(EGCG)、表儿茶素没食子酸酯(ECG)、茶多酚(TP)和游离氨基酸(AA)。其中咖啡碱的浓度采用高效液相色谱法(GB/T 8312-2013)测定,没食子酸、儿茶素、表没食子儿茶素没食子酸酯和表儿茶素没食子酸酯的浓度采用高效液相色谱法(GB/T8313-2013)测定,茶多酚总量采用福林酚法(GB/T 8313-2013)测定,游离氨基酸(AA)总量采用茚三酮法(GB/T 8314-2013)测定。
所述的步骤⑶中,针对7种滋味特征成分,采用SPSS统计软件分别建立茶汤泡数x(x为大于1的自然数)和特征成分浓度y之间的回归方程。特征成分浓度y包括咖啡碱的浓度y1、没食子酸的浓度y2、儿茶素的浓度y3、表没食子儿茶素没食子酸酯的浓度y4、表儿茶素没食子酸酯的浓度y5、茶多酚总量的浓度y6、游离氨基酸总量的浓度y7,由此分别建立咖啡碱的浓度y1和茶汤泡数x的回归方程Y1、没食子酸的浓度y2和茶汤泡数x的回归方程Y2、儿茶素的浓度y3和茶汤泡数x的回归方程Y3、表没食子儿茶素没食子酸酯的浓度y4和茶汤泡数x的回归方程Y4、表儿茶素没食子酸酯的浓度y5和茶汤泡数x的回归方程Y5、茶多酚总量的浓度y6和茶汤泡数x的回归方程Y6、游离氨基酸总量的浓度y7和茶汤泡数x的回归方程Y7。
对于七个回归方程,Y1有两个待定系数a1、b1;Y2有两个待定系数a2、b2;Y3有两个待定系数a3、b3;Y4有两个待定系数a4、b4;Y5有两个待定系数a5、b5;Y6有两个待定系数a6、b6;Y7有两个待定系数a7、b7;
对于茶叶样品A1~Am,根据咖啡碱浓度的测量值,分别得到m个咖啡碱的回归方程Y1(A1)、Y1(A2)、……、Y1(Am),每个方程有2个系数,方程Y1(A1)的系数为a1(A1)、b1(A1),方程Y1(A2)的系数为a1(A2)、b1(A2),……,方程Y1(Am)的系数为a1(Am)、b1(Am);
对于茶叶样品A1~Am,根据没食子酸浓度的测量值,分别得到m个没食子酸的回归方程Y2(A1)、Y2(A2)、……、Y2(Am),每个方程有2个系数,方程Y2(A1)的系数为a2(A1)、b2(A1),方程Y2(A2)的系数为a2(A2)、b2(A2),……,方程Y2(Am)的系数为a2(Am)、b2(Am);
对于茶叶样品A1~Am,根据儿茶素浓度的测量值,分别得到m个儿茶素回归方程Y3(A1)、Y3(A2)、……、Y3(Am),每个方程有2个系数,方程Y3(A1)的系数为a3(A1)、b3(A1),方程Y3(A2)的系数为a3(A2)、b3(A2),……,方程Y3(Am)的系数为a3(Am)、b3(Am);
对于茶叶样品A1~Am,根据表没食子儿茶素没食子酸酯浓度的测量值,分别得到m个表没食子儿茶素没食子酸酯的回归方程Y4(A1)、Y4(A2)、……、Y4(Am),每个方程有2个系数,方程Y4(A1)的系数为a4(A1)、b4(A1),方程Y4(A2)的系数为a4(A2)、b4(A2),……,方程Y4(Am)的系数为a4(Am)、b4(Am);
对于茶叶样品A1~Am,根据表儿茶素没食子酸酯浓度的测量值,分别得到m个表儿茶素没食子酸酯的回归方程Y5(A1)、Y5(A2)、……、Y5(Am),每个方程有2个系数,方程Y5(A1)的系数为a5(A1)、b5(A1),方程Y5(A2)的系数为a5(A2)、b5(A2),……,方程Y5(Am)的系数为a5(Am)、b5(Am);
对于茶叶样品A1~Am,根据茶多酚总量的测量值,分别得到m个茶多酚总量的回归方程Y6(A1)、Y6(A2)、……、Y6(Am),每个方程有2个系数,方程Y6(A1)的系数为a6(A1)、b6(A1),方程Y6(A2)的系数为a6(A2)、b6(A2),……,方程Y6(Am)的系数为a6(Am)、b6(Am);
对于茶叶样品A1~Am,根据游离氨基酸总量浓度的测量值,分别得到m个游离氨基酸总量的回归方程Y7(A1)、Y7(A2)、……、Y7(Am),每个方程有2个系数,方程Y7(A1)的系数为a7(A1)、b7(A1),方程Y7(A2)的系数为a7(A2)、b7(A2),……,方程Y7(Am)的系数为a7(Am)、b7(Am);
进一步,对于炒青茶叶,回归方程的形式如下:
咖啡碱的浓度变化回归方程Y1=a1eb1x;没食子酸的浓度变化回归方程Y2=a2eb2x;儿茶素的浓度变化回归方程Y3=a3eb3x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4=a4eb4x;表儿茶素没食子酸酯的浓度变化回归方程Y5=a5eb5x;茶多酚总量的浓度变化回归方程Y6=a6eb6x;游离氨基酸总量的浓度变化回归方程Y7=a7xb7,a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7均为回归方程的待定系数。
对于其他种类的茶叶,回归方程的形式有多种形式,如回归方程、幂函数方程、指数方程等,每个回归方程都有2个待定系数。
进一步,对于炒青茶的茶叶样品A1,建立得到咖啡碱的浓度变化回归方程Y1(A1)=a1(A1)eb1(A1)x,没食子酸的浓度变化回归方程Y2(A1)=a2(A1)eb2(A1)x;儿茶素的浓度变化回归方程Y3(A1)=a3(A1)eb3(A1)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(A1)=a4(A1)eb4(A1 )x;表儿茶素没食子酸酯的浓度变化回归方程Y5(A1)=a5(A1)eb5(A1)x;茶多酚总量的浓度变化回归方程Y6(A1)=a6(A1)eb6(A1)x;游离氨基酸总量的浓度变化回归方程Y7(A1)=a7(A1)xb7(A1),a1(A1)、b1(A1)、a2(A1)、b2(A1)、a3(A1)、b3(A1)、a4(A1)、b4(A1)、a5(A1)、b5(A1)、a6(A1)、b6(A1)、a7(A1)、b7(A1)均为茶叶样品A1的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数。对于茶叶样品A2,咖啡碱的浓度变化回归方程Y1(A2)=a1(A2)eb1(A2)x;没食子酸的浓度变化回归方程Y2(A2)=a2(A2)eb2(A2)x;儿茶素的浓度变化回归方程Y3(A2)=a3(A2)eb3(A2)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(A2)=a4(A2)eb4(A2)x;表儿茶素没食子酸酯的浓度变化回归方程Y5(A2)=a5(A2)eb5(A2)x;茶多酚总量的浓度变化回归方程Y6(A2)=a6(A2)eb6(A2)x;游离氨基酸总量的浓度变化回归方程Y7(A2)=a7(A2)xb7(A2),a1(A2)、b1(A2)、a2(A2)、b2(A2)、a3(A2)、b3(A2)、a4(A2)、b4(A2)、a5(A2)、b5(A2)、a6(A2)、b6(A2)、a7(A2)、b7(A2)均为茶叶样品A2的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数。
以此类推,对于炒青茶叶样品Am,咖啡碱的浓度变化回归方程Y1(Am)=a1(Am)eb1(Am)x;没食子酸的浓度变化回归方程Y2(Am)=a2(Am)eb2(Am)x;儿茶素的浓度变化回归方程Y3(Am)=a3(Am)eb3(Am)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(Am)=a4(Am)eb4(Am)x;表儿茶素没食子酸酯的浓度变化回归方程Y5(Am)=a5(Am)eb5(Am)x;茶多酚总量的浓度变化回归方程Y6(Am)=a6(Am)eb6(Am)x;游离氨基酸总量的浓度变化回归方程Y7(Am)=a7(Am)xb7(Am),a1(Am)、b1(Am)、a2(Am)、b2(Am)、a3(Am)、b3(Am)、a4(Am)、b4(Am)、a5(Am)、b5(Am)、a6(Am)、b6(Am)、a7(Am)、b7(Am)均为茶叶样品Am的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数。
所述步骤(4)中,同类的多种茶叶中,分类主要是以加工工艺为分类依据,根据加工工艺将茶叶分为绿茶、红茶、黑茶、白茶、黄茶和乌龙茶六大类,例如其中绿茶的特征工艺为杀青。而根据加工工艺绿茶又分为炒青绿茶、烘青绿茶和蒸青绿茶等,例如其中炒青绿茶的特征工艺为以炒干的形式进行干燥。炒青绿茶是一类,属于这类的茶叶很多,如龙井、香茶等等。
所述步骤(5)中,将茶叶样品A1、A2、~Am的7种滋味成分的浓度变化与茶汤泡数之间的回归方程的系数与第一泡茶汤的对应滋味成分的测量值建立修正函数,优选按以下方法进行:
(i)将茶叶样品A1、A2、~Am分别得到的m个关于咖啡碱的浓度变化与茶汤泡数之间的回归方程Y1(A1)、Y1(A2)、……、Y1(Am)中的系数a1(A1)、a1(A2)、……、a1(Am)为纵坐标Fa1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、……、Y1(Am,x=1)为横坐标z,建立Fa1与z之间的回归方程,即为系数a1的修正函数;修正函数也有多种形式,根据回归拟合结果而定,修正函数为指数函数、对数函数、线性函数或幂函数函数。
同样的,以系数b1(A1)、b1(A2)、……~b1(Am)为纵坐标Fb1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、~Y1(Am,x=1)为横坐标z,建立Fb1与z之间的回归方程,即为系数b1的修正函数;
(ii)参照上述步骤(i),分别得到没食子酸的浓度变化与茶汤泡数之间的回归方程的系数a2、b2与第一泡茶汤的没食子酸的测量值之间的修正函数Fa2、Fb2;儿茶素的浓度变化与茶汤泡数之间的回归方程的系数a3、b3与第一泡茶汤的儿茶素的测量值之间的修正函数Fa3、Fb3;表没食子儿茶素没食子酸酯的浓度变化与茶汤泡数之间的回归方程的系数a4、b4与第一泡茶汤的表没食子儿茶素没食子酸酯的测量值之间的修正函数Fa4、Fb4;表儿茶素没食子酸酯的浓度变化与茶汤泡数之间的回归方程的系数a5、b5与第一泡茶汤的表儿茶素没食子酸酯的测量值之间的修正函数Fa5、Fb5;茶多酚的浓度变化与茶汤泡数之间的回归方程的系数a6、b6与第一泡茶汤的茶多酚的测量值之间的修正函数Fa6、Fb6;游离氨基酸的浓度变化与茶汤泡数之间的回归方程的系数a7、b7与第一泡茶汤的游离氨基酸的测量值之间的修正函数Fa7、Fb7;
进一步,对于炒青绿茶的茶叶样品,其步骤(i)优选按以下步骤进行:
将炒青绿茶的茶叶样品A1、A2、~Am分别得到的七个关于咖啡碱的浓度变化与茶汤泡数之间的回归方程Y1(A1)=a1(A1)eb1(A1)x、Y1(A2)=a1(A2)eb1(A2)x、~、Y1(Am)=a1(Am)eb1(Am)x中的系数a1(A1)、a1(A2)、……~a1(Am)为纵坐标Fa1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、~Y1(Am,x=1)为横坐标z,建立Fa1与z之间的回归方程,即为系数a1的修正函数;同样的,以系数b1(A1)、b1(A2)、……~b1(Am)为纵坐标Fb1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、~Y1(Am,x=1)为横坐标z,建立Fb1与z之间的回归方程,即为系数b1的修正函数;
进一步,步骤(5)中,所建立的预测模型进一步具体为:7种特征成分浓度变化与茶汤泡数之间的回归方程共七个:咖啡碱的浓度y1和茶汤泡数x的回归方程Y1、没食子酸的浓度y2和茶汤泡数x的回归方程Y2、儿茶素的浓度y3和茶汤泡数x的回归方程Y3、表没食子儿茶素没食子酸酯的浓度y4和茶汤泡数x的回归方程Y4、表儿茶素没食子酸酯的浓度y5和茶汤泡数x的回归方程Y5、茶多酚总量的浓度y6和茶汤泡数x的回归方程Y6、游离氨基酸总量的浓度y7和茶汤泡数x的回归方程Y7;所述回归方程为指数回归方程、线性回归方程或幂函数回归方程;
对于七个回归方程,Y1有两个待定系数a1、b1;Y2有两个待定系数a2、b2;Y3有两个待定系数a3、b3;Y4有两个待定系数a4、b4;Y5有两个待定系数a5、b5;Y6有两个待定系数a6、b6;Y7有两个待定系数a7、b7;
待定系数a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7分别由该系数的修正函数确定;修正函数是每种特征成分的回归方程的系数与该特征成分的第一泡茶汤的测量值之间的回归方程,a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7各自与第一泡茶汤的对应滋味成分的测量值之间的修正函数为Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7。
对于待测未知茶叶样品,按照步骤(1)的浸提方法对待测茶叶样品进行一次浸提,得到的第一泡茶汤按照步骤(2)进行定量分析,得到第一泡茶汤中含有的7种滋味特征成分浓度的测量值,将第一泡茶汤的测量值分别代入修正函数方程Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7,得到a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7,再将所得系数分别代入回归方程Y1、Y2、Y3、Y4、Y5、Y6、Y7,即得到该样品的7种特征成分浓度变化与茶汤泡数之间的回归方程,即为预测模型。
根据回归方程预测第二泡茶汤及以后泡数茶汤的特征成分的浓度,并可根据Dot值判定茶叶可冲泡次数,所述Dot值=第x泡下该滋味特征成分的浓度(mg/ml)/该滋味特征成分的味觉阈值,Dot值大于1表示人能感知到该成分的滋味,数值越大,味觉强度越大,小于1表示该化合物不呈味。
对于炒青绿茶,所建立的多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型如下:
针对7种滋味特征成分,七个回归方程如下:咖啡碱的浓度变化回归方程Y1=a1eb1x;没食子酸的浓度变化回归方程Y2=a2eb2x;儿茶素的浓度变化回归方程Y3=a3eb3x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4=a4eb4x;表儿茶素没食子酸酯的浓度变化回归方程Y5=a5eb5x;茶多酚总量的浓度变化回归方程Y6=a6eb6x;游离氨基酸总量的浓度变化回归方程Y7=a7xb7,a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7均为回归方程的待定系数。
待定系数a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7分别由该系数的修正函数确定:a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7各自与第一泡茶汤的对应滋味成分的测量值之间的修正函数为Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7,通过检测第一泡茶汤的滋味成分的测量值,代入修正函数得到回归方程的待定系数。
对于未知待测样品,按照步骤(1)的浸提方法对待测茶叶样品进行一次浸提,得到的第一泡茶汤按照步骤(2)进行定量分析,得到第一泡茶汤中含有的7种滋味特征成分浓度的测量值,然后将所得其中测量值分别代入修正函数方程Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7,得到a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7,再将所得系数代入回归方程,即得到待测茶叶样品的7种特征成分浓度变化与茶汤泡数之间的回归方程。
针对现有的技术空白,本发明创建了一种针对多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测方法,根据不同茶叶种类,可建立不同种类的茶叶如绿茶、乌龙茶的预测模型,以绿茶模型为例,对于随机抽取绿茶样本,可通过定量检测样品第1泡茶汤中各滋味特征成分浓度,即可预测第x泡茶汤中各滋味特征成分浓度。该方法具有成本低、效率高和预测准确等特点,其平均预测准确率高达92.0%,在随机样本(盲样)的茶汤第x泡滋味特征成分浓度预测,随机抽取样本(盲样)的耐泡特性预测,以及原叶茶多次冲泡所得茶汤的味觉可视化分析等方面具有重要应用价值和指导意义。
具体实施方式
下面结合实施例对本发明的技术方案作进一步描述,实施例是在本发明技术方案前提下实施的,给出了详细的实施方式和步骤,但本发明的保护范围不限于下述的实施例。
实施例1
建立炒青茶叶预测模型:
炒青样品A1选取为香茶
(1)浸提茶汤
(a)取3g茶叶样品加入100℃蒸馏水150ml后摇匀,浸提5min后迅速用脱脂棉过滤,滤液冷却后用蒸馏水定容至200ml,标记为第1泡茶汤,
(b)过滤时得到的叶底重复步骤(a)进行二次浸提,得到第2泡茶汤,多次重复步骤(a),分别得到第3至第4泡茶汤。
(2)对步骤⑴所得不同浸提次数的茶汤进行定量分析,分析其中含有的7种滋味特征成分浓度;7种滋味特征成分为咖啡碱(Caf)、没食子酸(GA)、儿茶素(C)、表没食子儿茶素没食子酸酯(EGCG)、表儿茶素没食子酸酯(ECG)、茶多酚(TP)和游离氨基酸(AA)。其中咖啡碱的浓度采用高效液相色谱法(GB/T 8312-2013)测定,没食子酸、儿茶素、表没食子儿茶素没食子酸酯和表儿茶素没食子酸酯的浓度采用高效液相色谱法(GB/T 8313-2013)测定,茶多酚总量采用福林酚法(GB/T 8313-2013)测定,游离氨基酸(AA)总量采用茚三酮法(GB/T8314-2013)测定。结果如表1所示:
表1.炒青绿茶样本A1的第1-4泡茶汤特征成分浓度(mg/ml)
(3)针对步骤⑵中的7种滋味特征成分,采用SPSS统计软件分别建立特征成分浓度变化与茶汤泡数x(x为大于1的自然数)之间的线性方程,所得线性方程如下:
咖啡碱的回归方程为y=1.4727e-0.974x(R2=0.9956);没食子酸的回归方程为y=0.2417e-1.113x(R2=0.9967);儿茶素的回归方程为y=0.0362e-0.584x(R2=0.9858);表没食子儿茶素没食子酸酯的回归方程为y=1.6176e-0.558x(R2=0.9845);表儿茶素没食子酸酯的回归方程为y=0.3191e-0.457x(R2=0.9782);茶多酚的回归方程为y=1.9713e-0.462x(R2=0.9569);游离氨基酸回归方程为y=0.1384x-1.913(R2=0.9993);
(4)选取同为炒青茶叶的瓶炒青、其他炒青1、其他炒青2、其他炒青3、其他炒青4的茶叶样品,重复步骤(1)~(3),建立上述茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程;
瓶炒青茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程如下:
咖啡碱的回归方程为y=2.1949e-1.133x(R2=0.9945);没食子酸的回归方程为y=0.0879e-1.049x(R2=0.9929);儿茶素的回归方程为y=0.0471e-0.656x(R2=0.9831);表没食子儿茶素没食子酸酯的回归方程为y=2.0556e-0.644x(R2=0.982);表儿茶素没食子酸酯的回归方程为y=0.4647e-0.588x(R2=0.9738);茶多酚的回归方程为y=1.6384e-0.384x(R2=0.9795);游离氨基酸回归方程为y=0.1288x-1.86(R2=0.9941);
其他炒青1茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程如下:
咖啡碱的回归方程为y=1.919e-1.034x(R2=0.9993);没食子酸的回归方程为y=0.2674e-1.137x(R2=0.9984);儿茶素的回归方程为y=0.0373e-0.502x(R2=0.9201);表没食子儿茶素没食子酸酯的回归方程为y=2.085e-0.608x(R2=0.999);表儿茶素没食子酸酯的回归方程为y=0.4157e-0.494x(R2=0.994);茶多酚的回归方程为y=2.2075e-0.483x(R2=0.9812);游离氨基酸回归方程为y=0.1235x-1.82(R2=0.9948);
其他炒青2茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程如下:
咖啡碱的回归方程为y=1.9862e-1.027x(R2=0.9996);没食子酸的回归方程为y=0.2817e-1.142x(R2=0.9999);儿茶素的回归方程为y=0.0519e-0.675x(R2=0.996);表没食子儿茶素没食子酸酯的回归方程为y=2.1558e-0.6x(R2=0.9983);表儿茶素没食子酸酯的回归方程为y=0.4586e-0.518x(R2=0.9985);茶多酚的回归方程为y=2.4357e-0.522x(R2=0.9837);游离氨基酸回归方程为y=0.1581x-2.013(R2=0.9966);
其他炒青3茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程如下:
咖啡碱的回归方程为y=2.4719e-1.198x(R2=0.9954);没食子酸的回归方程为y=0.0926e-1.092x(R2=0.9952);儿茶素的回归方程为y=0.0457e-0.656x(R2=0.9868);表没食子儿茶素没食子酸酯的回归方程为y=2.0823e-0.645x(R2=0.9914);表儿茶素没食子酸酯的回归方程为y=0.4573e-0.546x(R2=0.9879);茶多酚的回归方程为y=2.0278e-0.439x(R2=0.9608);游离氨基酸回归方程为y=0.1401x-1.925(R2=0.9936);
其他炒青4茶叶样品的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程如下:
咖啡碱的回归方程为y=2.2023e-1.129x(R2=0.9951);没食子酸的回归方程为y=0.0826e-1.035x(R2=0.9832);儿茶素的回归方程为y=0.0525e-0.681x(R2=0.9677);表没食子儿茶素没食子酸酯的回归方程为y=1.9478e-0.615x(R2=0.9536);表儿茶素没食子酸酯的回归方程为y=0.4186e-0.514x(R2=0.9465);茶多酚的回归方程为y=2.7885e-0.577x(R2=0.9996);游离氨基酸回归方程为y=0.135x-1.896(R2=0.9893);
(5)建立修正函数:
以咖啡碱为例,以咖啡碱的浓度变化与茶汤泡数之间的回归方程
Y1(A1)=1.4727e-0.974x;Y1(A2)=2.1949e-1.133x;Y1(A3)=1.919e-1.034x;Y1(A4)=1.9862e-1.027x;Y1(A5)=2.4719e-1.198x;Y1(A6)=2.2023e-1.129x中的系数1.4727(A1)、2.1949(A2)、1.919(A3)、1.9862(A4)、2.4719(A5)、2.2023(A6)为纵坐标Fa1,以茶叶样品A1、A2、A3、A4、A5、A6第一泡茶汤的咖啡碱的测量值0.520(A1,x=1)、0.731(A2,x=1)、0.661(A3,x=1)、0.699(A4,x=1)、0.696(A5,x=1)、0.739(A6,x=1)为横坐标z,采用SPSS统计软件建立F a1与z之间的回归方程,所得修正函数为Fa1=0.536e1.9648z(R2=0.7905),同样的,系数b1的修正函数为Fb1=-0.408ln(z)-1.2262(R2=0.7164)。
按照上述方法,最后得到修正函数如下:
没食子酸的修正函数Fa2=5.1993z1.2063(R2=0.9933),Fb2=-0.081ln(z)-1.3313(R2=0.8145);儿茶素的修正函数Fa3=1.331z0.9211(R2=0.9367),Fb3=-0.24ln(z)-1.525(R2=0.9489);表没食子儿茶素没食子酸酯的修正函数Fa4=1.8195z0.7339(R2=0.8114),Fb4=-0.199ln(z)-0.5884(R2=0.6334);表儿茶素没食子酸酯的修正函数Fa5=1.3782z0.891(R2=0.9157),Fb5=-0.227ln(z)-0.8217(R2=0.8463);茶多酚总量的修正函数Fa6=2.8688ln(z)+1.5709(R2=0.9727),Fb6=-0.466ln(z)-0.3793(R2=0.9124);游离氨基酸总量的修正函数Fa7=0.051e7.0175z(R2=0.7186),Fb7=-0.761ln(z)-3.3987(R2=0.7425)。
因此,得到炒青品种多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型如下:
由茶汤泡数x(x为大于1的自然数)与特征成分浓度y之间的浓度变化线性方程和修正函数组合而成。其中咖啡碱的预测模型Y1=a1eb1x,Fa1=0.536e1.9648z,Fb1=-0.408ln(z)-1.2262;没食子酸的预测模型Y2=a2eb2x,Fa2=5.1993z1.2063,Fb2=-0.081ln(z)-1.3313;儿茶素的预测模型Y3=a3eb3x,Fa3=1.331z0.9211,Fb3=-0.24ln(z)-1.525;表没食子儿茶素没食子酸酯的预测模型Y4=a4eb4x,Fa4=1.8195z0.7339,Fb4=-0.199ln(z)-0.5884;表儿茶素没食子酸酯的预测模型Y5=a5eb5x,Fa5=1.3782z0.891,Fb5=-0.227ln(z)-0.8217;茶多酚总量的预测模型Y6=a6eb6x,Fa6=2.8688ln(z)+1.5709,Fb6=-0.466ln(z)-0.3793;游离氨基酸总量的预测模型Y7=a7xb7,Fa7=0.051e7.0175z,Fb7=-0.761ln(z)-3.3987。
实施例2
随机炒青绿茶样本(盲样)的茶汤第x泡特征成分浓度预测,具体步骤如下:
⑴随机抽取某炒青绿茶样本3g置于三角瓶中,加入100℃蒸馏水150ml后摇匀,浸提5min后迅速用脱脂棉过滤,冷却后用蒸馏水定容至200ml,茶汤标记为第1泡。
⑵对第1泡茶汤中的滋味特征成分进行定量分析。采用高效液相色谱法(GB/T8312-2013)测定咖啡碱的浓度,采用高效液相色谱法(GB/T 8313-2013)测定没食子酸、儿茶素、表没食子儿茶素没食子酸酯和表儿茶素没食子酸酯的浓度,采用福林酚法(GB/T8313-2013)测定茶多酚总量,采用茚三酮法(GB/T 8314-2013)测定游离氨基酸(AA)总量。所得第一泡茶汤中咖啡碱、没食子酸、儿茶素、表没食子儿茶素没食子酸酯、表儿茶素没食子酸酯、茶多酚和游离氨基酸的测量值分别为0.520mg/ml、0.074mg/ml、0.020mg/ml、0.869mg/ml、0.195mg/ml、1.110mg/ml、0.137mg/ml。
⑶将滋味特征成分的浓度代入实施例1得到的多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型,首先代入修正模型,分别得到a1=1.490、b1=-0.960、a2=0.226、b2=-1.121、a3=0.036、b3=-0.588、a4=1.641、b4=-0.560、a5=0.321、b5=-0.451、a6=1.871、b6=-0.428、a7=0.133、b7=-1.884。再将系数代入回归方程,得到咖啡碱的预测模型为Y1=1.490e-0.960x;没食子酸的浓度变化回归方程Y2=0.226e-1.121x;儿茶素的浓度变化回归方程Y3=0.036e-0.588x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4=1.641e-0.560x;表儿茶素没食子酸酯的浓度变化回归方程Y5=0.321e-0.451x;茶多酚总量的浓度变化回归方程Y6=1.871e-0.428x;游离氨基酸总量的浓度变化回归方程Y7=0.133x-1.884。
根据预测模型预测第x(x=2,3,4,5,6,7,8)泡茶汤中滋味特征成分的浓度。如表2。
并按照步骤⑵对第2-4泡茶汤中的滋味特征成分进行定量分析,实测值结果见表2,结果表明,预测值与实测值误差较小。
表2.某随机炒青绿茶样本的茶汤第n泡特征成分浓度预测(mg/ml)
实施例3
随机抽取样本(盲样)的耐泡特性预测,具体步骤如下:
⑴随机抽取某炒青绿茶样本3g置于三角瓶中,加入100℃蒸馏水150ml后摇匀,浸提5min后迅速用脱脂棉过滤,冷却后用蒸馏水定容至200ml,茶汤标记为第1泡。
⑵对第1泡茶汤中的滋味特征成分进行定量分析。采用高效液相色谱法(GB/T8312-2013)测定咖啡碱的浓度,采用高效液相色谱法(GB/T 8313-2013)测定没食子酸、儿茶素、表没食子儿茶素没食子酸酯和表儿茶素没食子酸酯的浓度,采用福林酚法(GB/T8313-2013)测定茶多酚总量,采用茚三酮法(GB/T 8314-2013)测定游离氨基酸(AA)总量。所得第一泡茶汤中咖啡碱、没食子酸、儿茶素、表没食子儿茶素没食子酸酯、表儿茶素没食子酸酯、茶多酚和游离氨基酸的测量值分别为0.520mg/ml、0.074mg/ml、0.020mg/ml、0.869mg/ml、0.195mg/ml、1.110mg/ml、0.137mg/ml。
⑶将滋味特征成分的浓度代入多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型,先代入修正函数得到系数,再将系数代入回归方程,得到预测模型如下:
咖啡碱的预测模型为Y1=1.490e-0.960x;没食子酸的浓度变化回归方程Y2=0.226e-1.121x;儿茶素的浓度变化回归方程Y3=0.036e-0.588x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4=1.641e-0.560x;表儿茶素没食子酸酯的浓度变化回归方程Y5=0.321e-0.451x;茶多酚总量的浓度变化回归方程Y6=1.871e-0.428x;游离氨基酸总量的浓度变化回归方程Y7=0.665x-1.81。
根据预测模型预测第x(x=2,3,4,5,6)泡茶汤中滋味特征成分的浓度,所得结果如表3。
⑷通过Dot值(Dot值=第x泡下该滋味特征成分的浓度(mg/ml)/该滋味特征成分的味觉阈值,Dot值大于1表示人能感知到该成分的滋味,数值越大,味觉强度越大,小于1表示该化合物不呈味。)判断第x泡时7种滋味特征成分的呈味情况。其中咖啡碱的阈值为0.097mg/ml(Scharbert and Hofmann,2005),没食子酸的阈值为0.034mg/ml(Kaneko etal,2006),儿茶素、表没食子儿茶素没食子酸酯和表儿茶素没食子酸酯的阈值分别为0.011mg/ml、0.087mg/ml和0.115mg/ml(Scharbert et al,2004)。由如表3的预测结果可知,该炒青绿茶在冲泡5次以后,仍有EGCG的DOT值大于1,而在第六泡时,已经没有任何滋味成分的浓度能够超过阈值。
表3.随机抽取样本(盲样)的耐泡特性预测
Claims (9)
1.茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法,其特征在于所述方法包括以下步骤:
(1)对茶叶样品A1进行多次浸提,收集每次浸提得到的茶汤、定容并标记浸提次数;
(2)对步骤⑴所得不同浸提次数的茶汤进行定量分析,分析其中含有的7种滋味特征成分浓度;
(3)针对步骤⑵中的7种滋味特征成分,分别建立特征成分浓度变化与茶汤泡数之间的回归方程,
(4)选取与步骤(1)同类的多种茶叶样品A2~Am,重复步骤(1)~(3),建立茶叶样品A2~Am的7种滋味特征成分浓度变化与茶汤泡数之间的回归方程;m为3~10的自然数;
(5)将m种茶叶样品A1、A2、~Am的各自7种滋味成分的浓度变化与茶汤泡数之间的回归方程的系数与第一泡茶汤的对应滋味成分的测量值之间建立修正函数,建立相应的多次冲泡茶叶所得茶汤中滋味特征成分浓度的预测模型。
2.如权利要求1所述的方法,其特征在于所述预测模型包括:
7种特征成分浓度变化与茶汤泡数之间的回归方程共七个,所述回归方程为指数回归方程、线性回归方程或幂函数回归方程,每个回归方程各有2个待定系数,每个待定系数根据该系数的修正函数确定,修正函数是每种特征成分的回归方程的系数与该特征成分的第一泡茶汤的测量值之间的回归方程,每个待定系数对应一个修正函数,7个回归方程共有14个待定系数和14个对应的修正函数。
3.如权利要求2所述的方法,其特征在于所述预测模型按以下方法进行预测:
对于待测未知茶叶样品,按照步骤(1)的浸提方法对待测茶叶样品进行一次浸提,得到的第一泡茶汤按照步骤(2)进行定量分析,得到第一泡茶汤中含有的7种滋味特征成分浓度的测量值,将第一泡茶汤的测量值分别代入修正函数,获得各特征成分的回归方程的系数,再将系数分别代入7个回归方程,即得到该样品的7种特征成分浓度变化与茶汤泡数之间的回归方程,根据回归方程可预测第二泡茶汤及以后泡数茶汤的7种特征成分的浓度,并可根据Dot值判定茶叶可冲泡次数,所述Dot值=第x泡下该滋味特征成分的浓度/该滋味特征成分的味觉阈值,Dot值大于1表示人能感知到该成分的滋味,数值越大,味觉强度越大,小于1表示该化合物不呈味。
4.如权利要求1所述的方法,其特征在于所述步骤⑴中,茶叶样品A1的浸提步骤为:
(a)取3g茶叶样品A1加入100℃蒸馏水150ml后摇匀,浸提5min后迅速用脱脂棉过滤,滤液冷却后用蒸馏水定容至200ml,标记为第1泡茶汤,
(b)过滤时得到的叶底重复步骤(a)进行二次浸提,得到第2泡茶汤,多次重复步骤(a),分别得到第3至第n泡茶汤,n为大于3的自然数。
5.如权利要求1所述的方法,其特征在于所述的步骤⑵中,所述的7种滋味特征成分为咖啡碱、没食子酸、儿茶素、表没食子儿茶素没食子酸酯、表儿茶素没食子酸酯、茶多酚和游离氨基酸。
6.如权利要求1所述的方法,其特征在于所述的步骤⑶中,针对7种滋味特征成分,采用SPSS统计软件分别建立茶汤泡数x和特征成分浓度y之间的回归方程,x为大于1的自然数,特征成分浓度y包括咖啡碱的浓度y1、没食子酸的浓度y2、儿茶素的浓度y3、表没食子儿茶素没食子酸酯的浓度y4、表儿茶素没食子酸酯的浓度y5、茶多酚总量的浓度y6、游离氨基酸总量的浓度y7,由此分别建立咖啡碱的浓度y1和茶汤泡数x的回归方程Y1、没食子酸的浓度y2和茶汤泡数x的回归方程Y2、儿茶素的浓度y3和茶汤泡数x的回归方程Y3、表没食子儿茶素没食子酸酯的浓度y4和茶汤泡数x的回归方程Y4、表儿茶素没食子酸酯的浓度y5和茶汤泡数x的回归方程Y5、茶多酚总量的浓度y6和茶汤泡数x的回归方程Y6、游离氨基酸总量的浓度y7和茶汤泡数x的回归方程Y7;
对于七个回归方程,Y1有两个待定系数a1、b1;Y2有两个待定系数a2、b2;Y3有两个待定系数a3、b3;Y4有两个待定系数a4、b4;Y5有两个待定系数a5、b5;Y6有两个待定系数a6、b6;Y7有两个待定系数a7、b7;
对于茶叶样品A1~Am,根据咖啡碱浓度的测量值,分别得到m个咖啡碱的回归方程Y1(A1)、Y1(A2)、……、Y1(Am),每个方程有2个系数,方程Y1(A1)的系数为a1(A1)、b1(A1),方程Y1(A2)的系数为a1(A2)、b1(A2),……,方程Y1(Am)的系数为a1(Am)、b1(Am);
对于茶叶样品A1~Am,根据没食子酸浓度的测量值,分别得到m个没食子酸的回归方程Y2(A1)、Y2(A2)、……、Y2(Am),每个方程有2个系数,方程Y2(A1)的系数为a2(A1)、b2(A1),方程Y2(A2)的系数为a2(A2)、b2(A2),……,方程Y2(Am)的系数为a2(Am)、b2(Am);
对于茶叶样品A1~Am,根据儿茶素浓度的测量值,分别得到m个儿茶素回归方程Y3(A1)、Y3(A2)、……、Y3(Am),每个方程有2个系数,方程Y3(A1)的系数为a3(A1)、b3(A1),方程Y3(A2)的系数为a3(A2)、b3(A2),……,方程Y3(Am)的系数为a3(Am)、b3(Am);
对于茶叶样品A1~Am,根据表没食子儿茶素没食子酸酯浓度的测量值,分别得到m个表没食子儿茶素没食子酸酯的回归方程Y4(A1)、Y4(A2)、……、Y4(Am),每个方程有2个系数,方程Y4(A1)的系数为a4(A1)、b4(A1),方程Y4(A2)的系数为a4(A2)、b4(A2),……,方程Y4(Am)的系数为a4(Am)、b4(Am);
对于茶叶样品A1~Am,根据表儿茶素没食子酸酯浓度的测量值,分别得到m个表儿茶素没食子酸酯的回归方程Y5(A1)、Y5(A2)、……、Y5(Am),每个方程有2个系数,方程Y5(A1)的系数为a5(A1)、b5(A1),方程Y5(A2)的系数为a5(A2)、b5(A2),……,方程Y5(Am)的系数为a5(Am)、b5(Am);
对于茶叶样品A1~Am,根据茶多酚总量的测量值,分别得到m个茶多酚总量的回归方程Y6(A1)、Y6(A2)、……、Y6(Am),每个方程有2个系数,方程Y6(A1)的系数为a6(A1)、b6(A1),方程Y6(A2)的系数为a6(A2)、b6(A2),……,方程Y6(Am)的系数为a6(Am)、b6(Am);
对于茶叶样品A1~Am,根据游离氨基酸总量浓度的测量值,分别得到m个游离氨基酸总量的回归方程Y7(A1)、Y7(A2)、……、Y7(Am),每个方程有2个系数,方程Y7(A1)的系数为a7(A1)、b7(A1),方程Y7(A2)的系数为a7(A2)、b7(A2),……,方程Y7(Am)的系数为a7(Am)、b7(Am)。
7.如权利要求6所述的方法,其特征在于所述步骤(5)中,将茶叶样品A1、A2、~Am的7种滋味成分的浓度变化与茶汤泡数之间的回归方程的系数与第一泡茶汤的对应滋味成分的测量值建立修正函数,按以下方法进行:
(i)将茶叶样品A1、A2、~Am分别得到的七个关于咖啡碱的浓度变化与茶汤泡数之间的回归方程Y1(A1)、Y1(A2)、……、Y1(Am)中的系数a1(A1)、a1(A2)、……、a1(Am)为纵坐标Fa1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、……、Y1(Am,x=1)为横坐标z,建立Fa1与z之间的回归方程,即为系数a1的修正函数;
同样的,以系数b1(A1)、b1(A2)、……~b1(Am)为纵坐标Fb1,以茶叶样品A1、A2、~Am第一泡茶汤的咖啡碱的测量值Y1(A1,x=1)、Y1(A2,x=1)、~Y1(Am,x=1)为横坐标z,建立Fb1与z之间的回归方程,即为系数b1的修正函数;
(ii)参照上述步骤(i),分别得到没食子酸的浓度变化与茶汤泡数之间的回归方程的系数a2、b2与第一泡茶汤的没食子酸的测量值之间的修正函数Fa2、Fb2;儿茶素的浓度变化与茶汤泡数之间的回归方程的系数a3、b3与第一泡茶汤的儿茶素的测量值之间的修正函数Fa3、Fb3;表没食子儿茶素没食子酸酯的浓度变化与茶汤泡数之间的回归方程的系数a4、b4与第一泡茶汤的表没食子儿茶素没食子酸酯的测量值之间的修正函数Fa4、Fb4;表儿茶素没食子酸酯的浓度变化与茶汤泡数之间的回归方程的系数a5、b5与第一泡茶汤的表儿茶素没食子酸酯的测量值之间的修正函数Fa5、Fb5;茶多酚的浓度变化与茶汤泡数之间的回归方程的系数a6、b6与第一泡茶汤的茶多酚的测量值之间的修正函数Fa6、Fb6;游离氨基酸的浓度变化与茶汤泡数之间的回归方程的系数a7、b7与第一泡茶汤的游离氨基酸的测量值之间的修正函数Fa7、Fb7。
8.如权利要求6所述的方法,其特征在于,对于炒青茶叶,回归方程的形式如下:
咖啡碱的浓度变化回归方程Y1=a1eb1x;没食子酸的浓度变化回归方程Y2=a2eb2x;儿茶素的浓度变化回归方程Y3=a3eb3x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4=a4eb4x;表儿茶素没食子酸酯的浓度变化回归方程Y5=a5eb5x;茶多酚总量的浓度变化回归方程Y6=a6eb6x;游离氨基酸总量的浓度变化回归方程Y7=a7xb7,a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7均为回归方程的待定系数;
对于炒青茶的茶叶样品A1,建立得到咖啡碱的浓度变化回归方程Y1(A1)=a1(A1)eb1(A1)x,没食子酸的浓度变化回归方程Y2(A1)=a2(A1)eb2(A1)x;儿茶素的浓度变化回归方程Y3(A1)=a3(A1)eb3(A1)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(A1)=a4(A1)eb4(A1)x;表儿茶素没食子酸酯的浓度变化回归方程Y5(A1)=a5(A1)eb5(A1)x;茶多酚总量的浓度变化回归方程Y6(A1)=a6(A1)eb6(A1)x;游离氨基酸总量的浓度变化回归方程Y7(A1)=a7(A1)xb7(A1),a1(A1)、b1(A1)、a2(A1)、b2(A1)、a3(A1)、b3(A1)、a4(A1)、b4(A1)、a5(A1)、b5(A1)、a6(A1)、b6(A1)、a7(A1)、b7(A1)均为茶叶样品A1的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数;对于茶叶样品A2,咖啡碱的浓度变化回归方程Y1(A2)=a1(A2)eb1(A2)x;没食子酸的浓度变化回归方程Y2(A2)=a2(A2)eb2(A2 )x;儿茶素的浓度变化回归方程Y3(A2)=a3(A2)eb3(A2)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(A2)=a4(A2)eb4(A2)x;表儿茶素没食子酸酯的浓度变化回归方程Y5(A2)=a5(A2)eb5 (A2)x;茶多酚总量的浓度变化回归方程Y6(A2)=a6(A2)eb6(A2)x;游离氨基酸总量的浓度变化回归方程Y7(A2)=a7(A2)xb7(A2),a1(A2)、b1(A2)、a2(A2)、b2(A2)、a3(A2)、b3(A2)、a4(A2)、b4(A2)、a5(A2)、b5(A2)、a6(A2)、b6(A2)、a7(A2)、b7(A2)均为茶叶样品A2的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数;
以此类推,对于茶叶样品Am,咖啡碱的浓度变化回归方程Y1(Am)=a1(Am)eb1(Am)x;没食子酸的浓度变化回归方程Y2(Am)=a2(Am)eb2(Am)x;儿茶素的浓度变化回归方程Y3(Am)=a3(Am)eb3(Am)x;表没食子儿茶素没食子酸酯的浓度变化回归方程Y4(Am)=a4(Am)eb4(Am)x;表儿茶素没食子酸酯的浓度变化回归方程Y5(Am)=a5(Am)eb5(Am)x;茶多酚总量的浓度变化回归方程Y6(Am)=a6(Am)eb6 (Am)x;游离氨基酸总量的浓度变化回归方程Y7(Am)=a7(Am)xb7(Am),a1(Am)、b1(Am)、a2(Am)、b2(Am)、a3(Am)、b3(Am)、a4(Am)、b4(Am)、a5(Am)、b5(Am)、a6(Am)、b6(Am)、a7(Am)、b7(Am)均为茶叶样品Am的7种滋味成分的浓度变化与茶汤泡数x之间的回归方程的系数。
9.如权利要求7所述的方法,其特征在于所述步骤(5)中,所建立的预测模型具体为:7种特征成分浓度变化与茶汤泡数之间的回归方程共七个:咖啡碱的浓度y1和茶汤泡数x的回归方程Y1、没食子酸的浓度y2和茶汤泡数x的回归方程Y2、儿茶素的浓度y3和茶汤泡数x的回归方程Y3、表没食子儿茶素没食子酸酯的浓度y4和茶汤泡数x的回归方程Y4、表儿茶素没食子酸酯的浓度y5和茶汤泡数x的回归方程Y5、茶多酚总量的浓度y6和茶汤泡数x的回归方程Y6、游离氨基酸总量的浓度y7和茶汤泡数x的回归方程Y7;所述回归方程为指数回归方程、线性回归方程或幂函数回归方程;
对于七个回归方程,Y1有两个待定系数a1、b1;Y2有两个待定系数a2、b2;Y3有两个待定系数a3、b3;Y4有两个待定系数a4、b4;Y5有两个待定系数a5、b5;Y6有两个待定系数a6、b6;Y7有两个待定系数a7、b7;
待定系数a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7分别由该系数的修正函数确定;修正函数是每种特征成分的回归方程的系数与该特征成分的第一泡茶汤的测量值之间的回归方程,a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7各自与第一泡茶汤的对应滋味成分的测量值之间的修正函数为Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7;
对于待测未知茶叶样品,按照步骤(1)的浸提方法对待测茶叶样品进行一次浸提,得到的第一泡茶汤按照步骤(2)进行定量分析,得到第一泡茶汤中含有的7种滋味特征成分浓度的测量值,将第一泡茶汤的测量值分别代入修正函数方程Fa1、Fb1、Fa2、Fb2、Fa3、Fb3、Fa4、Fb4、Fa5、Fb5、Fa6、Fb6、Fa7、Fb7,得到a1、b1、a2、b2、a3、b3、a4、b4、a5、b5、a6、b6、a7、b7,再将所得系数分别代入回归方程Y1、Y2、Y3、Y4、Y5、Y6、Y7,即得到该样品的7种特征成分浓度变化与茶汤泡数之间的回归方程,即为预测模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710675588.4A CN107490635B (zh) | 2017-08-09 | 2017-08-09 | 茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710675588.4A CN107490635B (zh) | 2017-08-09 | 2017-08-09 | 茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107490635A true CN107490635A (zh) | 2017-12-19 |
CN107490635B CN107490635B (zh) | 2020-03-03 |
Family
ID=60644067
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710675588.4A Active CN107490635B (zh) | 2017-08-09 | 2017-08-09 | 茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107490635B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161166A (zh) * | 2019-05-23 | 2019-08-23 | 广西轻工业科学技术研究院有限公司 | 一种茶多酚中egcg含量的快速检测方法 |
CN111596007A (zh) * | 2020-05-13 | 2020-08-28 | 广东工业大学 | 一种基于口感预认识的寻找冲泡茶叶条件最佳参数的算法 |
CN114343030A (zh) * | 2022-02-16 | 2022-04-15 | 恩施市紫竹茶业有限公司 | 一种一叶七味的制茶工艺 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1763543A (zh) * | 2005-09-27 | 2006-04-26 | 浙江大学 | 红茶品质鉴定方法 |
CN104133047A (zh) * | 2014-07-10 | 2014-11-05 | 中国农业科学院茶叶研究所 | 一种分析茶汤滋味强度的方法 |
CN104569328A (zh) * | 2013-10-28 | 2015-04-29 | 云南天士力帝泊洱生物茶集团有限公司 | 普洱晒青茶滋味品质和香气品质的测定方法 |
CN105651712A (zh) * | 2015-12-29 | 2016-06-08 | 安徽农业大学 | 一种绿茶涩味强度的定量判别方法 |
-
2017
- 2017-08-09 CN CN201710675588.4A patent/CN107490635B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1763543A (zh) * | 2005-09-27 | 2006-04-26 | 浙江大学 | 红茶品质鉴定方法 |
CN104569328A (zh) * | 2013-10-28 | 2015-04-29 | 云南天士力帝泊洱生物茶集团有限公司 | 普洱晒青茶滋味品质和香气品质的测定方法 |
CN104133047A (zh) * | 2014-07-10 | 2014-11-05 | 中国农业科学院茶叶研究所 | 一种分析茶汤滋味强度的方法 |
CN105651712A (zh) * | 2015-12-29 | 2016-06-08 | 安徽农业大学 | 一种绿茶涩味强度的定量判别方法 |
Non-Patent Citations (3)
Title |
---|
凌彩金等: "茶叶审评技术研究进展", 《广东农业科》 * |
沙济琴等: "茶叶矿质元素浸出量与茶叶冲泡时间及冲泡次数的关系", 《福建师范大学学报》 * |
邵晓林等: "西湖龙井茶主要呈味物质浸出浓度与速率的研究", 《茶叶》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161166A (zh) * | 2019-05-23 | 2019-08-23 | 广西轻工业科学技术研究院有限公司 | 一种茶多酚中egcg含量的快速检测方法 |
CN111596007A (zh) * | 2020-05-13 | 2020-08-28 | 广东工业大学 | 一种基于口感预认识的寻找冲泡茶叶条件最佳参数的算法 |
CN114343030A (zh) * | 2022-02-16 | 2022-04-15 | 恩施市紫竹茶业有限公司 | 一种一叶七味的制茶工艺 |
Also Published As
Publication number | Publication date |
---|---|
CN107490635B (zh) | 2020-03-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Adnan et al. | Chemical composition and sensory evaluation of tea (Camellia sinensis) commercialized in Pakistan | |
CN104106696B (zh) | 一种茶饮料及其制造方法 | |
TWI386165B (zh) | Containers mixed with tea drinks | |
CN107490635A (zh) | 茶汤中滋味特征成分浓度及茶叶可冲泡次数的预测方法 | |
CN110477172A (zh) | 一种快速冷萃咖啡的制备方法 | |
KR101314132B1 (ko) | 헛개나무열매를 포함하는 라떼 음료 및 이의 제조방법 | |
Severini et al. | How the variance of some extraction variables may affect the quality of espresso coffees served in coffee shops | |
Caracostea et al. | Determination of caffeine content in Arabica and Robusta green coffee of indian origin | |
CN102870925A (zh) | 茉莉花茶的加工方法 | |
JP5681367B2 (ja) | 飲料用原料茶の製造方法 | |
TWI568363B (zh) | Green tea beverage and its manufacturing method | |
CN108157551A (zh) | 一种采用冰萃技术制得的茶萃取液、茶饮料及其制备方法 | |
JP2018033366A (ja) | コーヒー飲料の製造方法 | |
Garbowska et al. | The content of minerals, bioactive compounds and anti-nutritional factors in tea infusions | |
Ma et al. | Effect of brewing conditions using a single‐serve coffee maker on black tea (Lapsang Souchong) quality | |
Fadhil et al. | Sensory Assessment of Gayo arabica coffee taste based on various varieties and manual brewing devices | |
CN109329504A (zh) | 一种白茶冷泡茶 | |
CN102711501A (zh) | 容器装焙煎茶饮料 | |
Gebretsadik et al. | Levels of metals in commercially available Ethiopian black teas and their infusions | |
SA00210428B1 (ar) | منتج شاي | |
CN110646267A (zh) | 一种武夷肉桂茶等级及烘焙程度判别方法 | |
CN107198401B (zh) | 数字化冲泡普洱茶的方法及辅助泡茶用具 | |
Derun et al. | The determination of the boron amounts of teas that are sold in Turkey by using the ICP-OES technique | |
EP2800474A1 (en) | Coffee formulation | |
Fibrianto et al. | Brewing characterization for optimum functional properties of Dampit Robusta (Coffea canephora) and Liberica (Coffea Liberica) coffee leaves tea |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |