首先对41个食用油的原始光谱进行均值化处理，去除噪声。初始化窗口宽度为11个波长变量，窗口宽度增加的步长为10波长变量，依次建立了窗口宽度11～481的多个PLS模型。其中在窗口宽度为121个光谱数据点时，如图9-31所示：横轴为波数点值，纵轴为RMSECV的值，由RESECV随窗口位置变化的关系图可以计算得到最小的RMSECV。波段选取结果如图9-32所示。对应的波数点范围是5288～6213cm-1。
图9-31 RMSECV随窗口位置变化的关系图
图9-32 MWPLS法波段选择图
对之前按Kennard-Stone法划分好的样本进行建模，MWPLS法建模结果如图9-33所示。所建模型的主成分数为12时，模型最佳，决定系数R2为0.9514，RMSECV为0.6206，RMSEP为0.4668。
图9-33 MWPLS法棕榈酸模型分析结果
2.基于iPLS法的棕榈酸近红外光模型优化
首先对41个食用油的原始光谱进行均值化处理，去除噪声。将原始光谱分成2～50个子区间分别建模比较。其中将原始光谱分成19个区间得到最佳模型，图9-34所示为各个子区间模型的RMSECV的比较图。由图可知在第15个区间得到最小的RMSECV。波段选取结果如图9-35所示。该区间对应的波数范围是5789～6229cm-1。
图9-34 子区间模型的RMSECV的比较图
图9-35 iPLS法波段选择图
对之前按Kennard-Stone法划分好的样本集进行建模，iPLS法建模结果如图9-36所示。所建模型的主成分数为11时，模型最佳，决定系数R2为0.9432，RMSECV为0.6679，RMSEP为0.4437。
图9-36 iPLS法棕榈酸模型分析结果
3.基于BiPLS法的棕榈酸近红外光模型优化
首先对41个食用油的原始光谱进行均值化处理，去除噪声。将原始光谱分成2～55个子区间，选择每个划分区间下的最佳联合子区间分别建模比较。因为算法的计算量会随着划分区间数与联合子区间数的增大而剧烈增加，所以选择的子区间数应小于5。综合考虑RMSECV及RMSEP的值及联合子区间数，将原始光谱分成28个区间，4个子区间联合（14、16、22、23），波段选取结果如图9-37所示，区间对应的波数范围为5504～5797cm-1、5804～6090cm-1、7610～7902cm-1、8211～8504cm-1。(www.xing528.com)
对之前按Kennard-Stone法划分好的样本集进行建模，BiPLS法建模结果如图9-38所示。所建模型的主成分数为10时，模型最佳，决定系数R2为0.9628，RMSECV为0.5351，RMSEP为0.4441。
图9-37 BiPLS法波段选择图
图9-38 BiPLS法棕榈酸模型分析结果
4.基于SiPLS法的棕榈酸近红外光模型优化
首先对41个食用油的原始光谱进行均值化处理，去除噪声。将原始光谱分成4～50个子区间，分别选择联合子区间个数为2、3、4，进行联合建模得到相对应的最佳模型，见表9-7。
表9-7 基于SiPLS法的子区间组合建模
综合考虑RMSECV及RMSEP的值及联合子区间数，将原始光谱分成35个区间，4个子区间联合（3、6、15、23），波段选取结果如图9-39所示。区间对应的波数范围为4478～4709cm-1、5673～5904cm-1、8302～8541cm-1。建模结果如图9-40所示。
图9-39 SiPLS法波段选择图
图9-40 SiPLS法棕榈酸模型分析结果
5.棕榈酸定量模型比较
采用上述4种方法筛选波长后建立的最佳PLS模型见表9-8。从表中可以得出通过对划分区间数、联合区间数、区间选取等筛选得到的波长建立PLS回归模型，4种模型预测精度均明显优于全光谱建模。其中SiPLS（联合4个区间）所建模型指标最佳。决定系数R2为0.9446，RM-SECV为0.1812，RMSEP为0.1912。
表9-8 特征谱区筛选后建模结果比较
观察利用4种特征波长挑选方法优选的波数范围，4种方法建立的模型所对应波数范围的公共区域集中于5500～6100cm-1，而线性烷烃化合物的次甲基C-H的有两个主要吸收峰，在5680cm-1和5800cm-1处。其中5680cm-1处的吸收峰为亚甲基C-H反对称或者对称伸缩振动的一级倍频谱带，5800cm-1处的吸收峰为合频谱带。预测结果较好的SiPLS法和BiPLS法优选的波数范围在8302～8504cm-1也有公共区域，而链烷烃化合物中的端甲基C-H伸缩振动的二级倍频（8365～8400cm-1）是在该区域有主要吸收峰。本实验采用4种方法所挑选的特征波长与理论分析的特征峰相符。
免责声明：以上内容源自网络，版权归原作者所有，如有侵犯您的原创版权请告知，我们将尽快删除相关内容。
