全 文 :第31卷第6期
2013年12月
上 海 交 通 大 学 学 报 (农 业 科 学 版)
JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY(AGRICULTURAL SCIENCE)
Vol.31No.6
Dec.2013
文章编号:1671-9964(2013)06-0041-07 DOI:10.3969/J.ISSN.1671-9964.2013.06.008
收稿日期:2013-01-15
基金项目:上海市科学技术委员会资助(09dz1205003);国家科技支撑项目(2008BAJ10B04-5)
作者简介:张宝智(1987-),男,硕士生,研究方向:园林植物与观赏园艺,E-mail:baozhizhang@gmail.com;
刘群录(1970-)为本文通讯作者,男,博士,副教授,研究方向:园林植物应用,E-mail:liuql@sjtu.edu.cn
紫叶李叶片花色苷超声辅助有机溶剂提取工艺
张宝智1,刘群录1,2,3,4,张慧婷1
(1.上海交通大学 农业与生物学院,上海200240;2.农业部都市农业(南方)重点开放实验室,上海200240;
3.上海城市植物资源开发与应用工程技术研究中心,上海200231;
4.上海交通大学 陆伯勋食品安全研究中心,上海200240)
摘 要:在单因素试验的基础上,运用Box-Behnken中心组合试验和响应面法考察了液料比、超声
提取时间、温度3个因素对花色苷提取率的影响,并优化了提取工艺。结果表明,紫叶李叶片花色
苷最佳工艺条件为:提取剂0.1%盐酸甲醇,超声(功率300W)提取时间5min,温度78℃,液料比
51∶1(V/m)。在此条件下,花色苷的提取率预测值为2 291.96±72.37mg/kg,验证值为2 467.42
±69.58mg/kg,与预测值的相对误差为7.42%,说明利用响应面法优化超声辅助有机溶剂提取紫
叶李叶片花色苷工艺可行。
关键词:紫叶李;超声辅助提取;花色苷;响应面法
中图分类号:S687.9 文献标识码:A
Optimization of the Ultrasonic Extraction of Anthocyanins
from Prunus cerasifera Leares
ZHANG Bao-zhi1,LIU Qun-lu1,2,3,ZHANF Hui-ting
(1.School of Agriculture and Biology,Shanghai Jiaotong University,Shanghai 200240,China;
2.Key Lab of Urban Agriculture(South),Ministry of Agriculture,Shanghai Jiaotong University,Shanghai 200240,China;
3.Shanghai Engineering Research Center of Sustainable Plant Innovation,Shanghai 200231,China;
4.Bor.S.Luh Food Safety Research Center,Shanghai Jiaotong University,Shanghai 200240,China)
Abstract:Single-factor extraction tests and response surface methodology (RSM)were employed to
optimize the extraction conditions for anthocyanins from leaves of Prunus cerasifera.Extraction time,
solvent-solid ratio and ultrasound temperature were investigated.The results demonstrated that the
optimum extraction conditions for anthocyanins were as folows:extraction time 5 min,ultrasound
temperature 78℃ and liquid-solid ratio 51∶1 (V/m),ultrasound power 300 W.Under the optimum
conditions,the theoretical extraction rate was 2 291.96±72.37 mg/kg,while the extraction rate of
anthocyanins actualy measured was 2 467.42±69.58mg/kg,with relative error of 7.42%,which revealed
response surface methodology was feasible to optimize the ultrasonic assisted extraction of anthocyandins
from leaves of P.cerasifera using organic solvent.
Key words:Prunus cerasifera var.atropurpurea;ultrasonic assisted extraction;anthocyanins;response
surface methodology
上 海 交 通 大 学 学 报 (农 业 科 学 版) 第31卷
紫 叶 李 (Prunus cerasifera var.
atropurpurea),蔷 薇 科 (Rosaceae)李 亚 科
(Prunoideae)李属(Prunus)中的著名观叶树种,原
产亚洲西部,叶片全年呈紫红色或红色。紫叶李叶
片中含有丰富的花色苷,主要为矢车菊色素[1-2]。紫
叶李叶片转色期间中部功能叶花色苷含量最高,增
幅最大,其次为基部老叶、稍部新叶[3],叶片中花色
苷主要在上表皮层,栅栏组织和海绵组织中成片分
布,分布范围与花色苷含量呈正相关[4]。花色苷是
一种生物活性物质等,具有抗氧化[5]、清除自由
基[6]、抗突变活性[7]、减轻肝功能障碍[8]、以及抗肿
瘤[9]的功能,因此,紫叶李是极具开发前景的天然植
物色素资源。
超声辅助提取技术因其提取效率高、成本低廉,
在食品和植物化学领域得到广泛的应用[10-12]。但
紫叶李叶片花色苷超声辅助提取工艺的研究未见报
道。本文以紫叶李叶片为试验材料,采用超声波辅
助提取紫叶李叶片中的花色苷,并用响应面法优化
了提取条件,以期为紫叶李叶片中天然色素的进一
步开发利用提供参考。
1 材料与方法
1.1 试验材料
紫叶李(Prunus cerasifera var.atropurea)叶
片采自上海交通大学西校区,选择其中的4棵,于
08:00~10:00采样,每棵树在相同高度上从东、
南、西、北4个方向采集枝条中部的功能叶片(带小
枝),放入塑料封口袋,采后立即带回实验室,用自来
水冲洗3~4次,然后用蒸馏水冲洗l~2次,晾干,
去叶柄及叶脉后剪碎混匀,备用。
1.2 仪器与试剂
GeneQuant pro型紫外―可见光分光光度计
(Biochrom有限公司);PHB-4型便携式pH 计(上
海精密科学仪器有限公司);JY92-ⅡD型超声波提
取仪(宁波新芝);DK-S22型水浴锅(上海精宏实验
设备有限公司);TDL-40B型高速离心机(上海安
亭);其他试剂均为国产分析纯。
1.3 试验设计
1.3.1 测定指标与方法
称取一定质量的样品,按比例加入提取剂超
声处理,3 000r/min离心5min后取上清。花色
苷含量的测定采用pH 示差法[13],吸取花色苷提
取液1mL,分别用pH 1.0和pH 4.5缓冲液稀释
定容至3mL。室温静置15min,以蒸馏水作参比,
在最大吸收波长520nm和700nm处测定二者的
吸光值,计算花色苷含量。以矢车菊素-3-葡萄糖
苷作为花色苷标准,总花色苷含量以矢车菊素-3-
葡萄糖苷含量表示,用以下公式计算花色苷提
取率:
紫叶李叶片总花色苷提取率R (mg·kg-1)=
A
εL×MW×DF×
V
m×1000
;A=(A520pH1.0-A700
pH1.0)-(A520pH4.5-A700pH4.5),式中:A,吸光
值;ε,矢车菊-3-葡萄糖苷摩尔消光系数(29 600);
MW,花色苷分子量(449.2);DF,稀释倍数;V,体积
(mL);m,样品质量(g);L,光程(1cm)。
1.3.2 单因素试验
以剪碎的紫叶李叶片为原料,在恒温水浴条件
下,进行花色苷提取产量的单因素试验。提取剂的
种类分别为0.1%盐酸甲醇,乙醇,丙酮;液料比是
提取液体积(mL)比剪碎的紫叶李叶片质量(g),为
10∶1~70∶1,提取温度30~80℃,提取时间为1~
7min。
1.3.3 响应面法优化总花色苷提取条件
根据单因素实验结果,以紫叶李花色苷提取时
间(A)、料液比(B)、提取温度(C)3个因素,采用
Box-Behnken设计3因素3水平试验并与花色苷提
取产量进行响应面试验设计,结合响应面分析法对
紫叶李叶片中总花色苷的提取工艺进行优化,因素
水平设计见表1。
表1 紫叶李花色苷提取工艺优化响应面
试验因素及水平表
Tab.1 Coded values and corresponding actual values of
the optimization parameters used in
Box-Behnken design
因素Factors
水平Level
-1 0 1
A时间Time/min 4 5 6
B液料比Liquid solid ratio/(mL·g-1) 40 50 60
C温度Temperature/℃ 60 70 80
1.4 数据处理
通过Design Expert 8.0软件对实验数据进行
回归分析,研究各自变量交互作用并预测紫叶李花
色苷提取的最优工艺参数。所有试验均重复3次,
取平均值用于分析。采用Origin 8.0软件进行数据
整理和作图,用SAS 9.0软件进行统计分析,统计
24
第6期 张宝智,等:紫叶李叶片花色苷超声辅助有机溶剂提取工艺
方法采用One-Way ANOVA。
2 结果与分析
2.1 单因素试验
2.1.1 不同提取剂对紫叶李花色苷提取率的影响
不同有机溶剂对紫叶李叶片提取率的影响如
图1所示。由图1可知,提取剂0.1%的盐酸甲醇
提取效率显著高于乙醇和丙酮(P<0.01),酸化的
甲醇提取剂适用的基质范围较为宽泛[14-15],提取
效率高,并且有利于提高花色苷的稳定性,防止非
酰基化的花色苷降解[16];乙醇作为一种食品和医
药工业常用的绿色提取剂,环保无污染,符合人体
健康标准的要求,但是极性较弱,花色苷溶出效率
低;丙酮虽然能够溶解细胞壁中的果胶类物质[17],
加速色素溶出,防止因有机酸的使用而导致的花
色苷分子的甲酰化和乙酰化[18],但也有研究指出
在以丙酮为媒介的提取过程中,花色苷发生缩合
反应形成甲基化的原花青素[19],从而导致丙酮的
提取效率显著低于另外2种提取剂。因此,选择
0.1%的盐酸甲醇作为紫叶李叶片花色苷超声辅
助提取剂。
图1 提取溶剂对紫叶李叶片花色苷提取率的影响
注:值为3次重复的平均值±标准误(n=3),“*”表示不同提取
剂之间存在显著差异(P<0.05);“**”表示不同提取剂之间存在
极显著差异(P<0.01)。图2同。
Fig.1 Effect of extractant on extraction yield of
anthocyandins fromPrunus cerasiferaleaves
Notes:Values were means±SE(n=3).Significant differences
(P<0.05)between different extractant according to One-Way
ANOVA was shown by an asterisk.Significant differences(P<0.
01)between different extractant according to One-Way ANOVA was
shown by two asterisks.The same as Fig.2.
2.1.2 单因素处理对紫叶李花色苷提取率的影响
不同单因素(超声时间、温度、液料比)处理对紫
叶李叶片花色苷提取率的影响如图2所示。在提取
温度为40℃,超声的功率为300W,液料比为20∶1
的条件下,不同超声时间对紫叶李叶片总花色苷得
率的影响见图2A。由图2A可知,提取时间在1~
6min范围内,紫叶李花色苷提取产量呈上升趋势;
当提取时间超过5min时,原料中紫叶李花色苷的
溶出量达到平衡,花色苷提取产量无显著变化。花
色苷提取时间超过6min时,花色苷的提取率下降,
可能是由于环境中的光、氧和微生物等作用致使其
稳定性下降所致[20]。因此,紫叶李花色苷最佳提取
时间为5min。
在超声提取时间为3min,超声功率为300W,
提取温度为40℃的条件下,不同液料比对紫叶李叶
片总花色苷得率的影响见图2B。由图2B可知,在
液料比为10∶1~30∶1范围内,随料液比的增加,提
取率显著提高,液料比为30∶1~50∶1时,提取液中
花色苷含量上升的幅度减小,且相邻两处理间差异
不显著(P>0.01),但当液料比大于50l时,提取
率下降;前期溶剂量越大,整个溶媒的传质过程越
快,有效成分浸出越完全,提取率也就越大,但当溶
剂过大时,探头式超声能量对沉滞低层的紫叶李叶
片空化效应和机械效应等的强度减弱,影响了超声
提取的效果,同时会造成溶剂和能源的浪费。因此,
选用50∶1作为提取液料比。
在超声提取时间为3min,超声功率为300W,
液料比为20∶1的条件下,不同超声温度对紫叶李叶
片总花色苷得率的影响见图2C。温度在30~60
℃范围内,紫叶李花色苷提取产量显著升高(P<
0.01)。温度超过60℃时,提取产量增加不显著(P
>0.01)。其原因可能是较高的提取温度加速花色
苷的溶出,也可能使紫叶李矢车菊等色素结构发生
改变,生成无色的查尔酮式结构[21]。从生产成本考
虑,提取温度越高,能耗越大,同时温度过高会造成
部分花色苷分解而降低花色苷提取产量,故提取温
度以70℃为宜。
2.2 响应面法优化紫叶李叶片花色苷超声辅助提
取工艺
取3个自变量,分别为提取时间(A)、料液比
(B)、提取温度(C),以紫叶李叶片中总花色苷提取
率为响应值(R)进行3因素3水平响应面分析试
验,试验设计与结果见表2,方差分析结果见表3。
回归模型达到极显著水平(P<0.01),误差项
不显著,失拟项(P>0.05)不显著,回归模型的决定
系数为0.962,说明该模型能够解释96.2%的变化。
34
上 海 交 通 大 学 学 报 (农 业 科 学 版) 第31卷
图2 不同单因素(A:时间;B:液料比;C:温度)处理对紫叶李叶片花色苷提取率的影响
Fig.2 Effect of single factor treatment(A:time;B:liquid solid ratio;C:temperature)on
anthocyandins extraction yield of anthocyandins fromPrunus cerasifer leaves
表2 响应面试验方案及试验结果
Tab.2 The experiments design and results of
response surface methodology
试验
序号
Experiment
No.
A:时间
/min
Time
B:液料比
/(V∶m)
Liquid-
solid ratio
C:温度/℃
Temperature
提取率
/(mg·kg-1)
Extraction ratio
1 4 40 70 1 693.72
2 6 40 70 1 845.50
3 4 60 70 1 981.54
4 6 60 70 1 929.09
5 4 50 60 1 754.06
6 6 50 60 1 761.17
7 4 50 80 1 976.28
8 6 50 80 2 130.80
9 5 40 60 1 770.99
10 5 60 60 1 986.76
11 5 40 80 2 180.05
12 5 60 80 2 099.32
13 5 50 70 2 230.29
14 5 50 70 2 227.19
15 5 50 70 2 272.02
16 5 50 70 2 246.70
17 5 50 70 2 205.36
回归方程具有较高的拟合度和可信度,说明回归方
程与实际情况吻合较好,试验误差小,因此可用该回
归方程代替真实试验对试验结果进行分析,分析和
预测紫叶李叶片花色苷的超声辅助提取工艺。
利用Design Expert 8.0软件对表2数据进行
二次多元回归拟合,得到紫叶李花色苷提取产量预
测值对自变量A、B和C的二次多项回归方程:
提取率R=2 236.31+31.97×A+61.98×B
+141.16×C-49.74×A×B+34.23×A×C-
72.84×B×C-237.48×A2-137.68×B2-90.63
×C2;模型中一次项温度(C)、二次项A2、B2、C2,交
互项BC极显著(P<0.01),由此说明一次项、二次
项都对紫叶李叶片花色苷提取率有较显著影响,响
应值的变化相当复杂,各个具体试验因素对响应值
的影响不是简单的线性关系,简单的分析方法难以
解析,响应面法能较好地描述它们之间的关系,并能
较好地优化花色苷的提取工艺。
对回归方程中一次项系数的绝对值进行比较,
可判断因子影响的主次性,结合方差分析结果,确定
对紫叶李叶片花色苷提取率的影响从大到小的顺序
依次为提取温度、超声时间、料液比。
图3直观地给出了各个因素交互作用的响应面
3D图和等高线图。等高线的形状可反映出交互效
应的强弱,椭圆形表示两因素交互作用显著,而圆形
44
第6期 张宝智,等:紫叶李叶片花色苷超声辅助有机溶剂提取工艺
表3 二次响应面回归模型方差分析
Tab.3 Variance analysis of quadratic response surface regression model
方差来源
Source
平方和
Sum of squares
自由度
df
均方
Mean square
F值
Fvalue
Pvalue
Prob>F
显著性
Significance
Model 6.19E5 9 68 828.04 46.2 <0.000 1 **
A-时间 8 175.37 1 8 175.37 5.49 0.051 7
B-液料比 30 735.88 1 30 735.37 20.63 0.002 7 *
C-温度 1.59E5 1 1.59E5 106.99 <0.000 1 **
AB 9 897.27 1 9 897.27 6.64 0.036 6 *
AC 4 687.46 1 4 687.46 3.15 0.119 4
BC 21 225.58 1 2 122.58 14.25 0.006 9 **
A2 2.38E5 1 2.37E5 159.39 <0.000 1 **
B2 79 814.98 1 79 814.98 53.57 0.000 2 **
C2 34 585.17 1 34 585.17 23.21 0.001 9 **
Residual 10 429.30 7 1 489.90
Lack of Fit 7 968.83 3 2 656.28 4.32 0.095 8
Pure Error 2 460.47 4 615.12
Cor Total 6.29E5 16
则与之相反[22]。提取温度为70℃时,提取时间和
料液比对花色苷提取率的交互作用不显著[图3
(Ⅰ)]。适当加大料液比能提高花色苷的提取率,缩
短色素的提取时间。由响应曲面图可知,在一定的
料液比范围内,随着料液比的增大,紫叶李花色苷的
提取率增加;料液比超过一定范围后,花色苷的提取
率反而下降。原因可能是料液比过大时,色素浓度
下降,色素分子间作用力变小,其稳定性下降而容易
分解,导致色素提取率下降[23]。
从图3(Ⅱ)等高线图可以直观地看出,在液料
比为50∶l时,提取时间与温度对紫叶李花色苷提取
率的交互作用不显著。在提取温度60℃时,得到较
高色素提取产量需要5min;温度提高到70℃时,
提取时间也只需5min。这表明,在本试验水平范
围内,适当增加提取温度有利于紫叶李花色苷的溶
出。但提取温度过高,花色素变得不稳定,提取产率
反而下降。
提取温度和料液比对紫叶李花色苷提取率的交
互作用较为显著[图3(Ⅲ)]。提取温度较低时,要
得到较高的提取率必须加大料液比,但是液料比过
大,提取率反而下降;而提取温度高时,较低的料液
比就能达到较高的提取率。由此可见,适当提高
温度,降低提取液的用量,可提高花色苷的提取率。
从响应面的最高点和等高线图(图3)还可以看
出,在所选的范围内存在极值,即是响应面的最高
点。由Design Expert 8.0软件分析得到最佳提取
条件为提取时间5.13min、提取温度77.57℃、料
液比50.63∶1,理论总花色苷提取率为2 291.96±
72.37mg/kg。为操作方便,总花色苷提取优化工
艺参数选为提取时间5min、提取温度78℃、料液
比51∶1。在此条件下进行3次平行试验,总花色苷
提取率为2 467.42±69.58mg/kg,相对标准偏差
为7.42%。这表明响应面法优化紫叶李叶片中总
花色苷的提取工艺可靠,具有一定的实用价值。
3 结论
该试验确定了紫叶李叶片总花色苷的最佳提取
工艺条件为:提取剂0.1%的盐酸甲醇,提取(功率
300W)时间5min,提取温度78℃,料液比51∶1。
在此条件下,花色苷的得率预测值为2 291.96±
72.37mg/kg,验证值为2 467.42±69.58mg/kg,
与预测值的相对误差为7.42%,超声辅助有机溶剂
提取紫叶李叶片花色苷工艺可行。
54
上 海 交 通 大 学 学 报 (农 业 科 学 版) 第31卷
图3 不同自变量交互作用对紫叶李叶片总花色苷提取量影响的响应面和等高线图
(Ⅰ)提取时间与液料比;(Ⅱ)提取时间与温度;(Ⅲ)液料比与温度
Fig.3 Response surface and contour plots for the interaction effects of three independent variables on the extraction yield of
total anthocyandins fromPrunus cerasiferaleaves
(Ⅰ)Response surface method of time and liquid-solid ratio;(Ⅱ)Response surface method of time and temperature;(Ⅲ)Response surface
method of liquid solid ratio and temperature
参考文献:
[1] 朱书香.李属彩叶植物叶片成色机理的研究[D].保
定:河北农业大学,2010.
[2] 鲁仪增,亓文英,王晓云.紫叶李叶片花色素及花色苷
HPLC定量分析 [J].山东农业科学,2011(11):
100-102.
[3] 王庆菊.李属红叶树种叶片花色苷合成规律及影响因
素研究[D].泰安:山东农业大学,2007.
[4] 李亚蒙,王庆菊,沈向.四种李属彩叶树种枝叶花色苷
分布的解剖研究[J].山东农业大学学报(自然科学
版),2006,37(4):489-494.
[5] Chang Y C,Huang K X,Huang A C,et al.Hibiscus
anthocyanins-rich extract inhibited LDL oxidation and
oxLDL-mediated macrophages apoptosis [J].Food
and Chemical Toxicology,2006,44(7):1015-1023.
[6] Geng M J,Ren M X.Free radical scavenging activities
of pigment extract from Hibiscus syriacus L.petals in
vitro[J].African Journal of Biotechnology,2012,11
(2):429-435.
[7] Rosa R M,Melecchi M S,Abad F C,et al.Antioxi-
dant and antimutagenic properties of Hibiscus tilia-
ceus L.methanolic extract[J].Journal of Agricultural
and Food Chemistry,2006,54(19):7324-7330.
[8] Wang C J,Wang J M,Lin W L,et al.Protective effect
of Hibiscus anthocyanins against tert-butyl hydroper-
oxide-induced hepatic toxicity in rats[J].Food and
Chemical Toxicology,2000,38(5):411-416.
[9] Hou D X,Tong X H,Terahara N,et al.Delphinidin
3-sambubioside,aHibiscus anthocyanin,induces ap-
64
第6期 张宝智,等:紫叶李叶片花色苷超声辅助有机溶剂提取工艺
optosis in human leukemia cels through reactive oxy-
gen species-mediated mitochondrial pathway[J].Ar-
chives of Biochemistry and Biophysics,2005,440(1):
101-109.
[10] Huang W,Xue A.,Niu H,et al.Optimised ultrason-
ic-assisted extraction of flavonoids from Folium eu-
commiae and evaluation of antioxidant activity in
multi-test systems in vitro[J].Food Chemistry,2009,
114(3):1147-1154.
[11] Chen F,Sun Y,Zhao G,et al.Optimization of ultra-
sound-assisted extraction of anthocyanins in red rasp-
berries and identification of anthocyanins in extract u-
sing high-performance liquid chromatography-mass
spectrometry[J].Ultrasonics Sonochemistry,2007,14
(6):767-778.
[12] Ghafoor K,Choi Y H,Jeon J Y,et al.Optimization of
ultrasound-assisted extraction of phenolic com-
pounds,antioxidants,and anthocyanins from grape
(Vitis vinifera)seeds [J].Journal of Agricultural
and Food Chemistry,2009,57(11):4988-4994.
[13] Lee J,Durst R W,Wrolstad R E.Determination of to-
tal monomeric anthocyanin pigment content of fruit
juices,beverages,natural colorants,and wines by the
pH differential method:Colaborative study[J].Jour-
nal of AOAC International,2005,88(5):1269-1278.
[14] Longo L,Vasapolo G.Determination of anthocyanins
in Ruscus aculeatus L.berries[J].Journal of Agricul-
tural and Food Chemistry,2005,53(2):475-479.
[15] Bhandari P,Kumar N,Gupta A P,et al.A rapid RP
‐ HPTLC densitometry method for simultaneous de-
termination of major flavonoids in important medicinal
plants[J].Journal of Separation Science,2007,30
(13):2092-2096.
[16] 庞志申.花色苷研究概况[J].北京农业科学,2000
(5):37-42.
[17] Garcia-Viguera C,Zafrila P,Francisco A,et al.The
use of acetone as an extraction solvent for anthocya-
nins from strawberry fruit[J].Phytochemical Analy-
sis,1998,9(6):274-277.
[18] Timberlake C F,Bridle P.The anthocyanins of apples
and pears:The occurrence of acyl derivatives [J].
Journal of the Science of Food and Agriculture,1971,
22(10):509-513.
[19] MyjavcováR,Marhol P,Kren V,et al.Analysis of
anthocyanin pigments in Lonicera(Caerulea)extracts
using chromatographic fractionation folowed by mi-
crocolumn liquid chromatography-mass spectrometry
[J].Journal of Chromatography A,2010,1217(51):
7932-7941.
[20] 范龚健,韩永斌,顾振新.用响应面法优化红甘蓝色素
提取工艺参数[J].南京农业大学学报,2006,29(1):
103-107.
[21] Chigurupati N,Saiki L,Gayser C,et al.Evaluation of
red cabbage dye as a potential natural color for phar-
maceutical use[J].International Journal of Pharma-
ceutics,2002,241(2):293-299.
[22] 王允祥,吕风霞,陆兆新.杯伞发酵培养基的响应曲面
法优化研究[J].南京农业大学学报,2004,27(3):
89-94.
[23] Dyrby M,Westergaard N,Stapelfeldt H.Light and
heat sensitivity of red cabbage extract in soft drink
model system[J].Food Chemistry,2001,72:431-437.
74