DOI:10.3969/j.issn.2095-6002.2020.06.010
中图分类号:TS264.4
刘亚1, 雷声1, 朱大洲2, 高莉1, 刘国荣3, 王成涛3, 刘娟1, 郭青1
| 【作者机构】 | 1云南中烟工业有限责任公司技术中心; 2农业农村部食物与营养发展研究所; 3北京工商大学北京市食品添加剂工程技术研究中心 |
| 【分 类 号】 | TS264.4 |
| 【基 金】 | 云南中烟工业有限责任公司科技项目(2018CP03) 云南省科技厅青年项目(2018FD164) |
食品着色剂(食用色素)是食品添加剂的重要门类,广泛应用于食品、医药、化妆品、印染、烟草等工业领域[1-2]。在食品及包装中应用食用色素,可赋予其鲜艳的颜色,吸引消费者,刺激其视觉、听觉、嗅觉与味觉,在生理和心理方面影响消费行为[3-7]。近年来,发现多种天然食用色素还具有一定生理活性、保健功能,可作为食品营养强化剂[8-9]。
食用色素呈现出不同的颜色,是其对不同波长的光选择性吸收的结果。光谱曲线则是以物质的吸光度为依据,反映物质颜色特性的曲线。通常色素的光谱曲线也被称为色素颜色的“指纹”[3,10]。目前,天然食用色素的调配仍大多依靠人工配色,其高度依赖配色者的配色经验,受其心理因素、生理因素和其他条件的影响,生产效率低,产品质量稳定性差[11-12]。除人工配色外,计算机配色技术已有一些报道。近年来,基于光谱反射率重建模型建立RGB颜色空间到L*、a*、b*颜色空间的映射方法,用于颜色测量,并通过多元线性回归模型、神经网络模型等进行色素调配、织物染色的研究已见报道[13-16]。相关研究表明,采用光谱曲线对天然食用色素建立配色模型,可有效提高配色模型的精确度,减少预测颜色与目标颜色之间的色差[17]。
本文基于天然食用色素的质量浓度与其在特征最大吸收波长下吸光度呈正比关系的原理,利用光谱曲线对颜色描述的唯一性特性[17-18],构建了天然食用色素的2种配色模型。通过分析预测多元线性模型和神经网络模型对天然食用色素的复配效果,以及对比2种模型的预测性能,选出较优模型,并完成了对较优模型的复配验证。希望研究结果可以为天然食用色素的配色模型构建、智能配色等提供新的理论依据和技术参数。
红曲红、栀子黄、栀子蓝色素,河南中大恒源生物科技股份有限公司。
UV- 3600 Plus型紫外- 可见分光光度计,岛津仪器有限公司;CM- 3610A型色差仪,日本柯尼卡美年达公司;BSG- 300型光照培养箱,上海博讯实业有限公司医疗设备厂。
选择红曲红、栀子黄和栀子蓝三原色色素按照一定比例进行复配,得到复配色素溶液,共计156组样品,其制备配比见表1,其中样品1-31为单一色素溶液,样品32-156为复配色素溶液。
1.3.1 色素溶液的特征吸收光谱分析
紫外- 可见分光光度计于350~720 nm波长内,测定表1样品溶液的可见光吸收光谱,每样品测量3次,取其平均值进行训练和优化模型[18]。
1.3.2 单色及复配色素溶液的颜色参数测定
选用孔径9 mm、100% UV、D65光源、10°视角时,全自动色差仪检测表1中156组样品溶液的亮度CIEL*(L*)、红值(a*)和黄值(b*),以超纯水为空白,每个样品检测3次,取其平均值[17-18]。
1.4.1 多元线性回归模型
多元线性回归模型[14]可用于描述因变量y与多个自变量x1,x2,……xn之间的线性相关关系,模型见式(1)。
(1)
式(1)中,x1, x2,…,xk(k≥2)为回归因子;y为观测值;b0, b1, b2,…,bk为回归系数;ε~N(0, σ2)为随机误差。
1.4.2 神经网络模型
神经网络在近红外光谱分析中得到广泛应用,其主要由信号的正向传播和误差的反向传播组成[15-16]。正向传播时,光谱数据从输入层传入,经隐含层处理后传向输出层,若实际输出与期望输出的误差较大,此时将误差反向传播,以修改各个神经单元的权值和阈值,直到网络的输出误差减小到规定范围,或达到设定的学习次数为止。
在GB/T 5525—2008[19]中2.2节重铬酸钾溶液比色法基础上修改。准确称取0.1 g试样,于100 mL容量瓶中用适当溶剂稀释,用溶剂作参比液于1 cm比色皿中,并用紫外可见分光光度计于最大波长处测定其吸光度,然后计算单一色素的色价,见式 (2)。
色价![]()
(2)
式(2)中,Aλ为实测试样在样品的最大吸收波长λ下的吸光度;f为稀释倍数;m为试样质量。
根据CIELAB色差公式,计算色素溶液的色差,见式(3)。
(3)
式(3)中,ΔL*是指实际配方和预测配方溶液亮度的差值;Δa*是指实际配方和预测配方溶液红值的差值;Δb*是指实际配方和预测配方溶液黄值的差值。
为更加直观评价预测模型的精确度,采用平均绝对误差(MAE)、平均相对误差(MAPE)、均方误差(MSE)和均方根误差(RMSE)等4项指标,对模型的性能进行评估,见式(4)~式(7)。
(4)
(5)
(6)
(7)
式(4)~式(7)中,Yi表示实际值;Y′表示模型预测值;i=1, 2, 3,……n;n为数据量个数。
表1 不同质量浓度三原色色素的样品制备配比
Tab.1 Sample preparation proportion of trichromatic pigments with different concentration g·L-1
试样编号ρ(复配)红黄蓝试样编号ρ(复配)红黄蓝试样编号ρ(复配)红黄蓝10.25400530.05300.13410500.1400.19120.20800540.10700.06710600.2790.09630.11000550.02300.04610700.1590.21840.16200560.03700.26710800.3180.10950.14200570.07300.13410900.2120.2960.12200580.03700.21811000.4240.14570.08100590.07300.10911100.5450.39880.04100600.23000.39811200.5450.309900.3580610.23000.30911300.5450.2571000.4190620.23000.25711400.4120.3981100.2930630.30900.28011500.4120.3091200.4770640.30900.22111600.4120.2571300.5850650.30900.18611700.3320.3981400.6360660.0330.092011800.3320.3091500.1960670.0690.279011900.3320.2571600.1060680.1390.140012000.4390.2801700.2820690.0530.196012100.4390.2211800.1410700.1070.098012200.4390.1861900.2120710.0370.390012300.2370.19920000.215720.3880.545012400.1620.25121000.287730.3880.412012500.2690.37822000.201740.3880.332012600.5270.45723000.327750.2890.545012700.6550.18924000.401760.2890.412012800.3900.13425000.436770.2890.332012900.1950.26726000.196780.2300.54501300.6662.3841.42427000.106790.2300.41201310.8980.8030.95128000.282800.2300.33201320.2570.3680.28029000.141810.3090.36801330.2570.3680.22130000.212820.3090.27901340.2570.3680.18631000.245830.3090.22501350.2570.2790.280320.55200.613840.2710.20501360.2570.2790.221330.89800.162850.6310.35001370.2570.2790.186340.21200.160860.8650.23701380.2570.2250.280350.28500.276870.3980.52701390.2570.2250.221360.43100.160880.7390.10801400.2570.2250.186370.29600.160890.8011.43501410.1520.2250.280380.29600.207900.6561.19801420.1520.2250.186390.80100.856910.5521.02901430.1520.2790.280400.38800.398920.3661.43501440.8981.5930.951410.38800.309930.8010.72301451.2360.7230.856420.38800.257940.0910.10701460.3661.4350.345430.28900.398950.0460.24401470.0370.0980.134440.28900.309960.0790.04601480.0370.1190.145450.28900.257970.0690.16001490.0711.1250.561461.01800.537980.0730.19501500.0950.3750.374470.03100.092990.0370.31801510.0850.1190.145480.04600.0461000.0730.15901520.0850.0980.067490.03200.04610101.4351.7081530.0850.1950.145500.03200.06410201.7910.3591540.0530.1190.134510.06900.19110301.4351.1401550.0370.2120.145520.13900.09610400.9011.0691560.0370.2120.067
于350~720 nm波长内绘制分析红曲红、栀子黄和栀子蓝色素的吸收光谱曲线(见图1),确认单一色素的最大特征吸收峰,发现其最大特征吸收波长分别为442、485、596 nm。根据单一色素溶液样品1-31的光谱曲线数据,采用式(2)计算得到红曲红、栀子黄和栀子蓝色素的色价分别为95.33、20.90和30.13。
图1 不同质量浓度红曲红、栀子黄、栀子蓝色素的吸收光谱
Fig.1 Absorption spectral analysis of different concentrations of monascus red pigment,
gardenia yellow pigment and gardenia blue pigment
由图1可知,红曲红、栀子黄和栀子蓝色素样品溶液的吸光度与其质量浓度呈正比,随着色素质量浓度的增加,样品的吸光度也随之增加。因此,可通过色素质量浓度与色价相乘得积,计算得到3种单一色素的吸光浓度,用于构建预测模型。
按照不同质量比,对红曲红、栀子黄和栀子蓝3种色素进行复配,于350~720 nm测定复配色素溶液的吸收光谱,并绘制其吸收光谱曲线(见图2),发现复配色素溶液分别于442、485、596 nm呈现特征吸收峰,与单一色素的最大特征吸收峰基本一致,且3个特征吸收峰下的色素质量浓度与其吸光度呈正相关性。此结果进一步验证了光谱曲线对颜色描述的唯一性[12,17]。
图2 不同质量浓度三原色色素复配溶液的吸收光谱
Fig.2 Absorption spectral analysis of different compound concentrations of trichromatic pigments
根据图1(单一色素)和图2(复配色素)的吸收光谱,为适应色素复配的不同色价需求,扩大预测复配模型的适应性,以混合色素溶液在3个特征吸收波长的吸光度为自变量,混合色素溶液内单一色素溶液的色价质量浓度为因变量,构建多元线性和神经网络预测模型,并分析比较2个预测模型的适应性。
2.3.1 多元线性预测模型构建
由于不能明确各影响因素对混合色素溶液特征吸收值的影响是否显著,为建立最优回归模型,采用逐步回归方法,对各变量进行F检验,使用F值判定回归模型的变量参数,设置步进方法的F值为:0.05<F<0.1。运用SPSS 23.0软件作为多元线性拟合工具,为提高数据的准确性和模型的稳定性,每个数据为3次试验平行数据的平均值。
根据得到的多元线性拟合结果,得到变量移入移除表和模型拟合度检验表(见表2),根据R2 判定回归方程的拟合度,优化得到最优拟合方程。由表2中红曲红、栀子黄和栀子蓝色价质量浓度为因变量,得到优化拟合方程的R2 分别为0.992、0.971、0.977,调整后的R2为0.984、0.942、0.954,R2满足要求,且较接近1,说明构建的回归方程拟合度良好。这与丁海霞[14]的多元线性回归模型分析织物染色配色结果相一致。
根据前述优化后拟合方程,通过方差分析得到F值和显著性结果(见表3)。以红曲红、栀子黄和栀子蓝色价质量浓度为目标值的回归模型F值分别为482.44、415.88、802.94,P值均为0.000,说明建立的回归拟合模型具有统计学上的显著性,回归效果有效,回归方程拟合良好。SPSS软件模型处理得到的回归系数(见表4),系数的共线性统计量中,模型的容忍度均大于0.1,VIF(方差膨胀系数)均小于10,表明该模型不存在多重共线性。
2.3.2 多层前馈BP神经网络模型构建
采用IBM SPSS Modeler构建神经网络模型(见图3),将442、485、596 nm的复配色素溶液吸光度设置为神经网络模型输入层,将红曲红、栀子黄、栀子蓝色价质量浓度为神经网络模型输出层,创建标准模型,选择建立多层感知器的神经网络模型,隐藏层数为1个,隐藏层节点为软件自动计算,以不能进一步降低误差为中止规则,将数据以7∶3分为训练区、测试区。隐含层和输出层的传递函数分别为tanh函数,学习速率取0.05,训练目标为0.001,最大迭代次数为1 000,训练的误差精度为0.05。通过软件的多次迭代计算,最终确定当神经网络模型的隐藏层神经元为5个时,准确度和运行时间较佳,R2为0.957。
表2 多元回归分析模型拟合度检验和变量移入表
Tab.2 Fitting degree inspection and variable value of multiple regression analysis model
因变量类型因素水平R2调整的R2标准估计误差输入变量方法10.7460.5560.188Aλ=485步进红曲红色价质量浓度20.9620.9250.078Aλ=442步进30.9920.9840.036Aλ=596步进10.8950.8020.215Aλ=442步进栀子黄色价质量浓度20.9670.9340.124Aλ=485步进30.9710.9420.118Aλ=596步进栀子蓝色价质量浓度10.9740.9490.074Aλ=596步进20.9770.9540.071Aλ=442步进
表3 方差分析
Tab.3 Analysis of variance
因变量类型差异源平方和df均方FSig回归6.21532.072482.440.000红曲红色价质量浓度1残差0.0992310.001总计6.314233回归17.29825.766415.880.000栀子黄色价质量浓度2残差1.0682310.014总计18.365233回归8.06924.034802.940.000栀子蓝色价质量浓度3残差0.3922310.005总计8.46233
表4 多元回归系数
Tab.4 Coefficients of multiple regression
因变量类型因变量非标准化系数B标准误差标准系数tSig.共线性统计量容差VIF常量0.0520.0077.580.00红曲红色价质量浓度Aλ=4851.3160.0222.1160.810.000.1695.909Aλ=442-0.6740.018-1.35-38.070.000.1626.174Aλ=596-0.2060.012-0.28-17.070.000.7421.348常量0.0940.0234.160.00栀子黄色价质量浓度Aλ=4421.4870.0581.7525.590.000.1626.174Aλ=485-0.9340.071-0.88-13.150.000.1695.909Aλ=596-0.1260.040-0.10-3.170.000.7421.348常量-0.0600.013-4.760.00栀子蓝色价质量浓度Aλ=5960.7860.0240.9332.950.000.7431.347Aλ=4420.0480.0160.082.9230.000.7431.347
图3 BP神经网络模型
Fig.3 BP neural network model
构建的多元线性回归模型和神经网络模型,对预测验证区的156组数据进行预测红曲红、栀子黄、栀子蓝色素色价质量浓度。根据模型结果,通过计算多元线性模型和神经网络模型的平均绝对误差(MAE)、平均相对误差(MAPE)、均方误差(MSE)和均方根误差(RMSE)等4项指标评价模型的精确度(见表5),表明神经网络模型的各个误差指标均小于多元线性回归拟合模型,说明神经网络模型精度较高,且相比于多元线性回归模型较稳定。因此,采用BP神经网络模型作为天然色素的复配模型,该结果与对橄榄油中二元、三元混合物的拟合结果相一致[20-21]。
表5 模型预测性能对比
Tab.5 Comparison of model predictive ability
模型类型均方误差均方根误差平均绝对误差平均相对误差多元线性模型0.0130.3000.0890.437神经网络模型0.0030.2390.0570.347
采用神经网络模型进行复配模型的验证,从样品制备表中随机选择5个色素复配配方,根据获得样品的光谱吸收曲线;通过模型预测计算得到预测配方,使用全自动色差仪分别测量原配方和拟合配方获得的颜色(色调)参数;根据CIELAB色差公式[见式(3)],计算两者的色差(见表6)。由表6可知,根据神经网络模型计算得到预测配方与原配方溶液的平均色差为1.604,肉眼无法分辨两者的差别。
表6 原配方与拟合配方的色差比较
Tab.6 Comparison of color aberration between sample
formula and fitting formula
试样编号原始配方预测配方红曲红栀子黄栀子蓝红曲红栀子黄栀子蓝色差10.0970.7930.3600.0970.7700.3712.1120.0970.3970.3600.0980.4190.3581.3230.0070.1620.1040.0080.1490.1051.5840.0460.1620.0000.0470.1600.0000.6250.0430.6750.1350.0470.6690.1232.39
本文分别测定了红曲红、栀子黄和栀子蓝单一色素溶液及3种色素复配溶液的吸收光谱,通过分析吸收光谱发现,复配色素溶液的特征吸收峰与参与复配的单一色素的最大特征吸收峰基本一致,且色素溶液的质量浓度与其吸光度呈正比。根据色素复配溶液的特征吸收峰吸光度及其复配比例,构建了多元线性回归模型和神经网络模型,通过模型误差分析,表明神经网络模型相较于多元线性模型具有更高的精度和稳定性,神经网络模型适用于天然色素复配模型构建。通过试验验证得出,预测配方和原配方溶液的色差在3以内,肉眼无法区分两者的差别,表明使用该配色模型对天然食用色素溶液配色精确度、重复性良好,可以较好地满足实际生产及应用的需要。
[1] 孙宝国. 食品添加剂[M]. 北京:化学工业出版社,2008:25-58.
SUN B G. Food additives[M]. Beijing:Chemical Industry Press,2008:25-58.
[2] 孙宝国,曹雁平,李健,等. 食品科学研究前沿动态[J]. 食品科学技术学报,2014,2:1-8.
SUN B G,CAO Y P,LI J,et al. Frontier dynamic of food science research[J]. Journal of Food Science and Technology,2014,2:1-8.
[3] 蔡雨杭,陶荣静,郭荣辉. 天然染料的应用及发展[J]. 纺织科学与工程学报,2018,35(3):137-142.
CAI Y H,TAO R J,GUO R H. Application and development of natural dyes[J]. Journal of Textile Science and Engineering,2018,35(3):137-142.
[4] 曲兴卫,王自珍. 基于色彩心理学分析商品包装设计的色彩选择[J]. 智库时代,2018,39:139.
[5] CAO S,LI Y,MENG X,et al. Dietary natural products and lung cancer: effects and mechanisms of action[J]. Journal of Functional Foods,2019,52:316-331.
[6] RODRIGUEZ-AMAYA D B. Update on natural food pigments: a mini-review on carotenoids,anthocyanins,and betalains[J]. Food Research International,2019,124:200-205.
[7] PARMAR R S, SINGH C A. Comprehensive study of eco-friendly natural pigment and its applications[J]. Biochemistry and Biophysics Reports,2018,13:22-26.
[8] 袁超,鲁晶鸽. 天然食用色素的功能及应用研究进展[J]. 粮食与油脂,2015,2:5-8.
YUAN C,LU J G. Research progress on functions and application of natural edible pigments[J]. Cereals & Oils,2015,2:5-8.
[9] 蒋新龙,蒋益花. 黑豆皮花色苷酯化修饰及其降解与抗氧化特性[J]. 中国粮油学报,2018,33(9):34-41.
JIANG X L,JIANG Y H. The black soybean peel anthocyanin process of esterification modified and its degradation and antioxidant characteristics[J]. Journal of the Chinese Cereals and Oils Association,2018,33(9):34-41.
[10] 宋杰. 生态印染技术在流行服装图案设计中的应用[J]. 染整技术,2018,40(10):53-56.
SONG J. Application of ecological printing and dyeing technologyin fashion pattern design[J]. Textile Dyeing and Finishing Journal,2018,40(10):53-56.
[11] 叶清珠,陈东生,罗奋涛,等. 天然栀子黄色素的提取及其对棉织物的染色特征值研究[J]. 毛纺科技,2018,46(7):40-43.
YE Q Z,CHEN D S,LUO F T,et al. Study on natural gardenia yellow dyestuff extracting and dyeing characteris-tic value on cotton fabric[J]. Wool Textile Journal,2018,46(7):40-43.
[12] 朱莉娜,鲁丹丹,鲍正壮. 天然染料橄榄叶色素提取及其真丝绸染色[J]. 印染助剂,2018,35(6):45-51.
ZHU L N,LU D D,BAO Z Z. Extraction of pigments of canarium album leaves and its dying performance on silk fabric[J]. Textile Auxiliaries,2018,35(6):45-51.
[13] 马崇启,程璐,金晓,等. Stearns-Noechel模型的全光谱纱线配色算法[J]. 光谱学与光谱分析,2018,38(8):2488-2492.
MA C Q,CHENG L,JIN X,et al. All spectral yarn color matching algorithm based on Stearns-Noechel model[J]. Spectroscopy and Spectral Analysis,2018,38(8):2488-2492.
[14] 丁海霞. 回归分析在织物染色计算机配色中的应用研究[D]. 青岛:青岛大学,2011.
DING H X. Research on application of regression analysis in computer color matching for textile dyeing[D]. Qingdao: Qingdao University,2011.
[15] 文孟俊. 基于神经网络算法的胶印专色油墨配色研究[D]. 济南:齐鲁工业大学,2013.
WEN M J. Study on the offset spot ink color matching based on neural network algorithm[D]. Jinan: Qilu University of Technology,2013.
[16] 朱道路. 基于卷积神经网络的颜色测量模型及对中小规模样本处理策略的研究[D]. 北京:北京交通大学,2018.
ZHU D L. Research on the color measurement model based on convolutional neural network and processing strategy of small and medium scale samples[D]. Beijing:Beijing Jiaotong University,2018.
[17] 方恩印. 基于光谱曲线油墨配色方法的研究[J]. 出版与印刷,2016,3:43-45.
[18] 米琳. 烟用天然色素的筛选及其在再造烟叶的应用[D]. 北京:北京工商大学,2019.
[19] 国家质量监督检验检疫总局,国家标准化管理委员会. 植物油脂 透明度、气味、滋味鉴定法:GB/T 5525—2008[S].北京:中国标准出版社,2008:1-3.
[20] AROCA-SANTOS R,LASTRA-MEJIAS M,CANCILLA J C,et al. Linear and non-linear quantification of extra virgin olive oil, soybean oil, and sweet almond oil in blends to assess their commercial labels[J]. Journal of Food Composition and Analysis,2019,75:70-74.
[21] AROCA-SANTOS R,CANCILLA J C,PEREZ-PEREZ A,et al. Quantifying binary and ternary mixtures of monovarietal extra virgin olive oils with UV-vis absorption and chemometrics[J]. Sensors and Actuators B:Chemical,2016,234:115-121.
LIU Ya, LEI Sheng, ZHU Dazhou, et al. Comparison of color matching between multivariate linear model and neural network model of natural food pigments[J]. Journal of Food Science and Technology, 2020,38(6):76-83.
X