Optimization of Full-spectrum Color Matching Model of Wood Veneer Dyeing Based on ADE Algorithm
-
摘要: 采用自适应的交叉因子与变异因子在增强全局搜索能力的同时提高收敛速度,通过循环迭代取值的方法确定Stearns−Noechel模型中参数M的最优值。采用新模型对8组标准样进行染色配方预测,以CIEDE 2000色差评价标准对实验结果进行评价。结果表明:基于ADE优化模型预测染色配方得到的拟合样与标准样间ΔE00均小于3;迭代次数是遗传算法的13.04%,差分进化算法的50.00%,寻优速度更快;对于标准样1,新模型预测配方得到拟合样与标准样间ΔE00约为最小二乘法、遗传算法和差分进化算法优化模型的1/5、1/3与1/2,预测精度更高,说明了新模型在染色木材单板全光谱配色方面具有较高的使用价值。Abstract: The adaptive crossover factor and mutation factor are used to enhance the global search ability and improve the convergence speed. The optimal value of the parameter M in the Stearns−Noechel model is determined by the method of loop iteration. The new model was used to predict the dyeing formula of 8 groups of standard samples, and the experimental results were evaluated by the CIEDE2000 color difference evaluation standard. The experimental results show that the ΔE00 between the fitted sample and the standard sample obtained by predicting the dyeing formula based on the ADE algorithm optimization model is less than 3. The number of iterations is 13.04% of the genetic algorithm and 50.00% of the differential evolution, and the optimization speed is faster. For standard sample 1, the new formula predicts that the ΔE00 between the fitted sample and the standard sample is about 1/5, 1/3, and 1/2 of the least squares, GA and DE algorithm optimization models, the prediction accuracy is higher, which shows that the new model has higher value in the full spectrum color matching of dyed wood veneer.
-
木材材色是决定木制品价值的重要因素,为了提高木制品的装饰作用和产品价值,对木材及木质材料进行适当的材色改进十分重要,通过染色技术改良劣质材、仿制优质材,从而全面提高木材的价值,满足人们对色彩多样性的需求。木材染色中一个十分重要环节就是配色,其对染色后木材的材色质量至关重要[1-2]。采用人工配色,其对配色人员的素质要求较高,既费时又难以适应现代工业生产的要求,且成本高、准确性差。有学者提出将计算机配色的方法用于木材染色中加快染料配方生成的速度,提高工作效率,节约成本[3]。然而,计算机配色技术为纺织业的配色而研制,利用库贝尔卡−芒克理论对测得的颜色相关量进行计算,采用了大量的假设,自适应能力差,难于精确获得木材染色这类影响因素较多,过程变化复杂的染色配方。近年来,有学者将神经元网络等智能算法应用于木材染色计算机配色模型的建立中,取得了一定的效果[4]。但仅从标准观察者与施照态条件下进行了配方的预测,染色木材单板存在的同色异谱现象亟待解决。
物体表面的光谱反射率是表征物体表面颜色的固有属性,根据该原则,Stearns−Noechel光学模型在计算机配色应用中具有较高的配色精度,有效地解决了配色过程中存在的同色异谱现象。但由于在全光谱波段进行预测,方程组规模特别巨大,但在配色过程中染料组分的数量通常不超过5,导致染料配方预测模型是病态且超定的,无法获得预测配方精确解。马崇启等[5]中针对混色纱线配色问题采用传统的LS法对超定方程组进行求解,但随着染料种类的增加以及光谱分辨率的提高,方程规模非常巨大,常规的解法会消耗大量的时间与内存,而且LS方法获得的模糊度浮点解与真值相差较大,严重影响了预测精度[6]。宋丛威等[7]提出了一种基于主成分分析的快速降维算法求解印染行业中的线性超定方程组,虽然缩短了计算时间,但在主成分增多时会导致算法结果不稳定。随着智能算法的兴起,将差分进化算法(DE)应用于优化求解过程逐渐成为研究热点。DE是一种基于优胜劣汰的自然选择原理对种群中的个体进行淘汰和保留的算法,算法本身具有一定的自组织、自学习与自适应的特点,但也不可避免的存在陷入局部最优解和搜索停滞的问题[8]。本研究通过对DE的变异因子和交叉因子进行改进,提出了一种自适应差分进化算法,使其随着迭代过程自适应变化,避免陷入局部最优解以及搜索停止问题,提高病态超定方程组解算精度,进而提高染色配方预测精度。
1. 实验材料和仪器
1.1 材料
选取水曲柳(Fraxinus mandshurica)单板为染色基材,购于黑龙江省哈尔滨市,气干旋切制成规格为60 mm×40 mm×0.5 mm的单板,干燥至含水率8%。
木材染色染料:活性红X−3B染料、活性黄X−RG染料、活性艳蓝染料,购于黑龙江省哈尔滨市。染色渗透剂JFC、固色剂(无水碳酸钠)、促染剂(无水乙酸)均购于北京化工厂,均为分析纯。
根据郭明辉等[9]中的染色工艺选取单位浓度的活性红X−3B染料、活性艳蓝染料、活性黄X−RG染料制备单一染料组分染色材,水曲柳单板及单一染料组分染色单板如图1所示。
标准样的染色基材选用水曲柳单板,将活性红X−3B染料、活性艳蓝染料、活性黄X−RG染料按比例进行组合,共得到36组染色方案,去掉其中相同比例的28组染色方案后,采用余下的8组染色方案用来制备标准样,染料比例如表1所示。
表 1 制备标准样染料比例Table 1. Preparation of standard sample dye ratio试样编号 R−Y−B染料比例 试样编号 R−Y−B染料比例 1 1∶1∶8 5 2∶2∶6 2 1∶2∶7 6 2∶3∶5 3 1∶3∶6 7 2∶4∶4 4 1∶4∶5 8 3∶3∶4 注∶R为活性红,Y为活性黄,B为活性艳蓝。 1.2 仪器
在颜色测量仪器选取中,通常采用RGB原理的彩色相机或分光光度计进行测量,但是彩色相机颜色感知精度较低,分光光度计虽具有较高的感知精度但仅能用于单点测量。高光谱成像技术具有高光谱分辨率,可实现成像与测谱同步进行。故本研究选用SPECIM高光谱成像工作站(SPECIM公司,芬兰)进行颜色测量,其搭载VNIR−FX10相机,配备双卤素灯线性光源(400~2 500 nm)主要工作参数见表2。
表 2 SPECIM高光谱成像工作站工作参数Table 2. SPECIM hyperspectral imaging workstation operating parameters工作参数 参数值 光谱范围 400~1 000 nm 扫描速度 65 mm/s 成像速度 330~9 900 Hz 视场 38° 空间分辨率 1 024 SPECIM高光谱成像工作站将可见光近红外光谱技术与高分辨率成像相结合,采用推扫式成像技术对实验样品进行逐线全波段光谱采集并同步生成图像,图像中每一像素点都记录了其对应实验样品颜色信息的光谱特征,适用于对实验样品颜色特征进行定性、定量分析。SPECIM高光谱成像工作站工作方式如图2所示。实验中标准样本对单点的光谱反射率进行3次测量取其平均值,单一组分染料染色材及拟合样本在其正反面各取6点测量取其平均值,在400~700 nm波段内,步长5 nm进行取值,得到共61个波长下的光谱反射率。
2. 实验模型与研究方法
2.1 Stearns−Noechel模型
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$f[{R_{\text{混}}}(\lambda )]{\rm{ = }}\sum\limits_i {{x_i}f} [{R_i}(\lambda )]$
(1) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
Stearns等[11]基于前人研究基础,通过实验验证建立了混合组分反射率计算经验公式,被称为Stearns−Noechel模型,如式(2)所示
$f[R(\lambda )] = \frac{{1 - R(\lambda )}}{{M[R(\lambda ) - 0.01] + 0.01}}$
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$R(\lambda ) = \frac{{0.01(M - 1)f[R(\lambda )] + 1}}{{Mf[R(\lambda )] + 1}}$
(3) 木材染色单板反射光谱是唯一决定其表面颜色的因素。当前常用的计算机配色方法为三刺激值匹配法,但染色得到的拟合样的反射光谱和标准样并不一致,而是利用同色异谱理论进行配色的,所染单板在测色阶段施照态和观察者两个角度与配色阶段环境条件不符时,等色即被破坏,从而出现色差。为达到更高的染色品质,实现木材染色单板与标样间的无条件匹配,根据参考文献[2]中所提出的全光谱匹配方法,使标样与拟合样的光谱反射率完全相同,如式(4)所示。
$R_\lambda ^{{\rm{Standard}}} \approx R_\lambda ^{{\rm{Fitting}}}$
(4) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${F^{({\rm{Standand}})}} \approx {F^{({\rm{Fitting}})}}$
(5) ${F^{{\rm{(Standand)}}}} \approx F \times X$
(6) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ {F^{({\rm{Standand}})}} = \left[ {\begin{array}{*{20}{c}} {f(R_{400}^{{\rm{Standand}}})}\\ {f(R_{405}^{{\rm{Standand}}})}\\ \vdots \\ {f(R_{700}^{{\rm{Standand}}})} \end{array}} \right] $
(7) $ {F^{({\rm{Fitting}})}} = \left[ {\begin{array}{*{20}{c}} {f(R_{400}^{{\rm{Fitting}}})}\\ {f(R_{405}^{{\rm{Fitting}}})}\\ \vdots \\ {f(R_{700}^{{\rm{Fitting}}})} \end{array}} \right] $
(8) $ F = \left[ {\begin{array}{*{20}{c}} {f\left( {R_{400}^1} \right)}&{f\left( {R_{400}^2} \right)}&{f\left( {R_{400}^3} \right)}\\ {f\left( {R_{405}^1} \right)}&{f\left( {R_{405}^2} \right)}&{f\left( {R_{405}^3} \right)}\\ \vdots & \vdots & \vdots \\ {f\left( {R_{700}^1} \right)}&{f\left( {R_{700}^2} \right)}&{f\left( {R_{700}^3} \right)} \end{array}} \right] $
(9) $ X = \left[ \begin{array}{l} {x_1}\\ {x_2}\\ {x_3} \end{array} \right] $
(10) $ f(R_\lambda ^{{\rm{Standand}}}) \approx f(R_\lambda ^{{\rm{Fitting}}}) = {x_1}f(R_\lambda ^1) + {x_2}f(R_\lambda ^2) + {x_3}f(R_\lambda ^3) $
(11) $ \begin{array}{l} {{x'}_1} = \dfrac{{{x_1}}}{{{x_1} + {x_2} + {x_3}}}\\ {{x'}_2} = \dfrac{{{x_2}}}{{{x_1} + {x_2} + {x_3}}}\\ {{x'}_3} = \dfrac{{{x_3}}}{{{x_1} + {x_2} + {x_3}}} \end{array} $
(12) $ f({R'_\lambda }) = {x'_1}f\left( {R_\lambda ^1} \right) + {x'_2}f\left( {R_\lambda ^2} \right) + {x'_3}f\left( {R_\lambda ^3} \right) $
(13) 2.2 自适应差分进化算法优化全光谱配色模型
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
DE的基本思路是,取2个个体之间的差分向量,将这个差分向量的1个加权形式加到第3个个体上从而产生1个新的候选解。DE由Storn等[12]提出,全局寻优能力强、收敛速度快并具有良好的鲁棒性。DE的主要步骤如下:
1)设置种群规模N、最大进化代数(MaxIter)、变异因子(F)和交叉因子(CR)等参数,初始化种群。
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${v_i} = {x_{{r_1}}} + F\left( {{x_{{r_2}}} - {x_{{r_3}}}} \right)$
(14) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
通过对DE的分析与介绍可以得到算法流程图如图3所示。
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$\begin{array}{l} {F^{\rm{*}}} = \left\{ \begin{array}{l} 2\left| {\cos \left( {\dfrac{3}{4}\pi \dfrac{{{N_i}}}{{{N_{\max }}}}} \right)} \right|,{N_i} \leqslant \dfrac{2}{3}{N_{\max }} \\ 16\left[{{\rm{e}}^{{{\left( {\frac{{{N_{\max }} - {N_i}}}{{2{N_{\max }}}}} \right)}^3}}}\right],{N_i} > \dfrac{2}{3}{N_{\max }} \\ \end{array} \right. \\ {\rm{CR}}{^{\rm{*}}} = \dfrac{4}{5} - \dfrac{1}{2}\cos \left( {\dfrac{{\pi {N_i}}}{{2{N_{\max }}}}} \right) \\ \end{array} $
(15) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
ADE具体步骤如下:
1)设置种群规模N、最大进化代数(MaxIter),进行种群初始化。
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${v_i} = {x_{{r_1}}} + {F^{\rm{*}}}\left( {{x_{{r_2}}} - {x_{{r_3}}}} \right)$
(16) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
6)判断全局最优解精度是否满足要求,是则输出结果,否则转入步骤2)继续执行。
通过对ADE算法的分析与介绍可以得到算法流程图如图4所示。
2.3 参数M最优值选取
M值作为Stearns−Noechel模型中的唯一可变参量,合理选取M值对达到标准样与拟合样颜色最优匹配具有重要意义。传统选取M值的方法是已知混色样品配方,将单组分样品反射率以及质量分数代入Stearns−Noechel模型中反推M值,通过匹配实验,选取所有试样平均色差最小时的M值为最优值[14]。韩瑞叶等[15]和沈加加等[16]研究发现模型参数M与波长在部分波段存在线性相关,采用数据拟合的方法建立最优M值与波长的线性相关方程,增加了模型的预测能力,但采取这种方式的最优M值并不适用于所有的实验样本,依然对拟合样的匹配造成了一定的误差。本研究采用王泉等[17]提出的循环迭代取值方法求取最优M值,经过针对木材单板配色最优M值采用反推法进行多次实验的基础上,确定M通常取值范围为[0, 1],设置初始值M0=0,步长为0.01循环迭代计算M值,选取色差最小时的M值作为最优M值,求得8组混合比例标准样最小色差时的M值为0.21。
综上,本研究提出的基于ADE算法优化Stearns−Noechel模型的单板染色全光谱配色过程见图5。
3. 实验结果与分析
为验证本研究提出方法的有效性,根据以上流程设计实验。采用SPECIM FX10高光谱成像工作站分别采集单组分染料染色单板及标准样的反射光谱曲线如图6所示。通过软件ENVI 4.8在400~700 nm波段每隔5 nm进行取值,进而得到单色染色单板光谱反射率矩阵(F)与标准样光谱反射率矩阵(F (Standand)),以标准样1为例,采用ADE算法对式(6)进行求解。
设置种群规模=50,最大迭代次数=100,传统DE算法中,变异因子=0.5和交叉因子=0.9,ADE算法中变异因子及交叉因子按式(15)自适应选择。
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
根据表3,利用ADE求解染色配方得到的拟合样与标准样1的色差最小,与传统的LS法、GA和基本DE相比,色差分别约为1/5、1/3、1/2。其中,利用LS法求解染色配方得到的拟合样与标准样1间的色差较大,产生的色差人眼可以分辨。
通过表4可以得出,8组实验标准样与拟合样的明度差绝对值均值为1.15,拟合样与标准样明度相近,其中第1、4、6组明度差为负值,表明拟合样颜色略深于标准样,第2、3、5、7、8组明度差为正值,表明拟合样颜色略浅于标准样;8组实验标准样与拟合样的彩度差绝对值均值为1.6875,拟合样与标准样颜色纯度相近;8组实验标准样与拟合样的色差均值为2.09,其中第4组拟合样与标准样的色差最大为2.9,第6组拟合样与标准样的色差最小为1.3,均满足木材染色单板生产颜色质量要求。
表 3 标准样1解算结果对比Table 3. Standard sample 1 solution comparison计算方法 拟合样染色配方 $\Delta {E_{00}}$ R Y B LS 0.062 7 0.190 7 0.746 6 11.7 GA 0.059 3 0.135 4 0.805 3 6.3 DE 0.062 2 0.184 3 0.753 5 3.9 ADE 0.051 1 0.166 4 0.782 5 1.9 注:R表示活性红,Y表示活性黄,B表示活性艳蓝。 表 4 ADE优化全光谱配色模型预测结果Table 4. ADE optimized full-spectrum color matching model prediction results试样编号 比例 拟合样 标准样 $\Delta {E_{00}}$ R Y B ${L^{'}}$ $C_{ab}^{'}$ $h_{ab}^{'}$ ${L^{'}}$ $C_{ab}^{'}$ $h_{ab}^{'}$ 1 0.051 1 0.166 4 0.782 5 28.4 22.8 293.0 30.5 21.2 295.3 1.9 2 0.112 8 0.196 3 0.690 9 31.2 20.5 300.3 29.8 22.0 300.7 1.4 3 0.106 7 0.331 8 0.561 5 30.3 8.5 320.6 28.6 10.4 310.6 2.4 4 0.097 2 0.362 7 0.540 1 30.6 17.8 307.2 31.7 15.0 316.1 2.9 5 0.224 1 0.187 4 0.588 5 32.4 22.7 303.6 31.4 19.2 305.2 2.0 6 0.235 8 0.320 3 0.443 9 31.6 18.9 306.5 31.9 19.3 304.1 1.3 7 0.197 6 0.386 4 0.416 0 29.4 23.1 334.4 28.5 24.7 329.5 2.0 8 0.314 7 0.287 6 0.397 7 30.4 21.8 327.8 29.7 21.6 319.5 2.8 基于Matlab 2016a平台,分别使GA、DE、ADE迭代次数满足搜索到最优解的要求,得到不同进行解算优化过程曲线如图7所示。其中,GA迭代次数约为230,DE迭代次数约为60,ADE迭代次数约为30。可见基于ADE优化Stearns−Noechel模型的木材单板染色全光谱配色可以高效解算出精准的染色配方。
4. 结论
通过高光谱成像技术对染色木材单板进行反射光谱采集,采用ADE优化Stearns−Noechel全光谱配色模型,并通过循环迭代法选取模型中M参数的最优值。将ADE的计算结果与其他优化算法的计算结果进行对比,可以看出,通过循环迭代法在[0, 1]区间内确定染色木材单板配方预测模型中参数M的最优值为0.21。通过采用不同优化算法对标准样1进行配方预测,LS、GA、DE和ADE预测配方得到的拟合样与标准样间的色差分别为11.7、6.3、3.9和1.9,结果表明,ADE优化模型得到的配方具有更高的预测精度。采用ADE优化模型对8组标准样进行配方预测,得到拟合样与标准样色差分别为1.9、1.4、2.4、2.9、2.0、1.3、2.0、2.8,色差均小于3,颜色差异人眼不可辨识。GA与DE分别在迭代230次与60次时搜索到最优解,而ADE搜索到最优解的迭代次数仅为30次,结果表明,相同条件下,ADE的寻优效率更高。
综上,基于ADE优化的染色木材单板全光谱配色模型可以有效地预测染色配方,且具有较高的计算效率与预测精度,满足染色木材单板颜色质量要求。
-
表 1 制备标准样染料比例
Table 1 Preparation of standard sample dye ratio
试样编号 R−Y−B染料比例 试样编号 R−Y−B染料比例 1 1∶1∶8 5 2∶2∶6 2 1∶2∶7 6 2∶3∶5 3 1∶3∶6 7 2∶4∶4 4 1∶4∶5 8 3∶3∶4 注∶R为活性红,Y为活性黄,B为活性艳蓝。 表 2 SPECIM高光谱成像工作站工作参数
Table 2 SPECIM hyperspectral imaging workstation operating parameters
工作参数 参数值 光谱范围 400~1 000 nm 扫描速度 65 mm/s 成像速度 330~9 900 Hz 视场 38° 空间分辨率 1 024 表 3 标准样1解算结果对比
Table 3 Standard sample 1 solution comparison
计算方法 拟合样染色配方 $\Delta {E_{00}}$ R Y B LS 0.062 7 0.190 7 0.746 6 11.7 GA 0.059 3 0.135 4 0.805 3 6.3 DE 0.062 2 0.184 3 0.753 5 3.9 ADE 0.051 1 0.166 4 0.782 5 1.9 注:R表示活性红,Y表示活性黄,B表示活性艳蓝。 表 4 ADE优化全光谱配色模型预测结果
Table 4 ADE optimized full-spectrum color matching model prediction results
试样编号 比例 拟合样 标准样 $\Delta {E_{00}}$ R Y B ${L^{'}}$ $C_{ab}^{'}$ $h_{ab}^{'}$ ${L^{'}}$ $C_{ab}^{'}$ $h_{ab}^{'}$ 1 0.051 1 0.166 4 0.782 5 28.4 22.8 293.0 30.5 21.2 295.3 1.9 2 0.112 8 0.196 3 0.690 9 31.2 20.5 300.3 29.8 22.0 300.7 1.4 3 0.106 7 0.331 8 0.561 5 30.3 8.5 320.6 28.6 10.4 310.6 2.4 4 0.097 2 0.362 7 0.540 1 30.6 17.8 307.2 31.7 15.0 316.1 2.9 5 0.224 1 0.187 4 0.588 5 32.4 22.7 303.6 31.4 19.2 305.2 2.0 6 0.235 8 0.320 3 0.443 9 31.6 18.9 306.5 31.9 19.3 304.1 1.3 7 0.197 6 0.386 4 0.416 0 29.4 23.1 334.4 28.5 24.7 329.5 2.0 8 0.314 7 0.287 6 0.397 7 30.4 21.8 327.8 29.7 21.6 319.5 2.8 -
蔡绍祥, 王新洲, 李延军. 高温水热处理对马尾松木材尺寸稳定性和材色的影响 [J]. 西南林业大学学报(自然科学), 2019, 39(1): 160−165. 秦磊, 卢丹钰, 何海珊, 等. PEG加固处理木质遗存变色的原因探析 [J]. 西南林业大学学报(自然科学), 2018, 38(5): 173−179. Liu W, Yu Z M, Li W Z. Preliminary study on application of computer color matching technique for wood dyeing [J]. Key Engineering Materials, 2011(476/474): 1228−1232.
Guan X M, Fan F X, Zhu Y R, et al. Application of RBF neural network optimized globally by genetic algorithm in intelligent color matching of wood dyeing [J]. Journal of Intelligent & Fuzzy Systems, 2017, 33(5): 2895−2901.
马崇启, 程璐, 金晓, 等. Stearns−Noechel模型的全光谱纱线配色算法 [J]. 光谱学与光谱分析, 2018, 38(8): 2488−2492. 纪元法, 朱亮亮, 孙希延, 等. 差分进化算法结合正则化解算病态方程 [J]. 系统工程与电子技术, 2018, 40(7): 1573−1577. DOI: 10.3969/j.issn.1001-506X.2018.07.23 宋丛威, 张晓明. 基于PCA的解大型超定线性方程组快速算法及应用 [J]. 智能计算机与应用, 2019, 9(4): 91−95. DOI: 10.3969/j.issn.2095-2163.2019.04.020 Zhang H F, Zhou J Z, Zhang Y C, et al. Short term hydrothermal scheduling using multi-objective differential evolution with three chaotic sequences [J]. International Journal of Electrical Power & Energy Systems, 2013, 47(1): 85−99.
郭明辉, 管雪梅. 桦木单板仿珍贵材三刺激值配色技术 [J]. 林业科学, 2010, 46(9): 130−135. DOI: 10.11707/j.1001-7488.20100921 章斐燕, 李启正, 张声诚, 等. 基于Stearns−Noechel优化模型的交织混色织物配色设计系统 [J]. 丝绸, 2015, 52(1): 26−30. DOI: 10.3969/j.issn.1001-7003.2015.01.005 Steams E I, Noechel F. Spectrophotometric prediction of color of wool blends [J]. American Dyestuff Reporter, 1994, 33(9): 177−180.
Storn R, Price K. Differential evolution: a simple and efficient heuristic for global optimization over continuous spaces [J]. Journal of Global Optimization, 1997, 11(4): 341−359. DOI: 10.1023/A:1008202821328
丁青锋, 尹晓宇. 差分进化算法综述 [J]. 智能系统学报, 2017, 12(4): 431−442. Burlone D A. Theoretical and practical aspects of selected fiber-blend color formulation functions [J]. Color Research & Application, 1984, 9(4): 213−219.
韩瑞叶, 杨瑞华, 薛元, 等. 数码转杯纺的Stearns−Noechel配色模型 [J]. 纺织学报, 2017, 38(12): 27−32, 42. 沈加加, 程凤侠, 陈维国, 等. Stearns−Noechel模型优化及毛条混色配色系统开发 [J]. 纺织学报, 2009, 30(3): 121−125. DOI: 10.3321/j.issn:0253-9721.2009.03.025 王泉, 姚佳, 李艳清, 等. Stearns−Noechel模型在天然彩色棉混色中的应用 [J]. 纺织学报, 2016, 37(1): 93−97.