CN105181641A - 一种菜籽油品质的近红外检测方法及应用 - Google Patents

Abstract

本发明属于油脂分析技术领域，具体涉及一种菜籽油品质的近红外检测方法及应用。包括以下步骤：收集菜籽油样品、对菜籽油样品进行理化测试、对菜籽油样品进行近红外扫描、对所得的吸光值进行异常样本剔除和样品集划分、预处理、对得到的预处理吸光值分析，筛选得到近红外光谱特征波长；通过逐步回归法建立菜籽油品质的预测模型为等步骤。本发明获得的信息直观可靠，确定了菜籽油品质的特征波长，且特征波长数量少，运用特征波长建立模型的分析方法，提高模型的精度。在相同预测精度的情况下，本发明的预测速度快，同时本发明建立的菜籽油分析比较全面，易于推广使用。

Description

一种菜籽油品质的近红外检测方法及应用

技术领域

本发明属于油脂分析技术领域，具体涉及一种菜籽油品质的近红外检测方法及应用。

背景技术

菜籽油是世界第三大植物油，对其品质的科学评价及监控一直都是热门课题，但目前菜籽油品质指标的检测存在操作复杂、费时、劳动强度大等缺点。近红外光谱分析技术基于其具有无损、绿色、高效、便捷等优点，已应用于检测油脂的脂肪酸组成、酸价、过氧化值、碘价等指标。在国外，RaoYL等采用偏最小二乘法构建菜籽油酸价的近红外定量模型，校正集相关系数达到0.972，验证集相关系数达到0.938，均取得了满意的效果；LiH采用偏最小二乘法构建了高碘价油脂的碘价近红外预测模型，相关系数达到0.99；在国内，崔晓君等同样采用偏最小二乘法构建了花生油过氧化值的近红外定量模型，校正集及验证集的相关系数均达到0.991，模型拟合精度非常高。申请号200910071567.7文献公布了一种基于近红外光谱分析的食用油脂酸价检测方法，采用的是在780～2500nm近红外全谱区的吸光信息，用逐步回归算法找特征波长，再用多元线性回归算法构建校正模型；申请号200910071566文献公布了一种基于近红外光谱技术快速检测食用油脂中反式脂肪酸含量的方法，同样采用的是多元线性回归的方法建立模型。

但这些研究尚存在光谱沉余信息较多，针对每种指标的近红外特征波长不明确，模型不清晰，同时对食用油色泽、黏度的近红外定量模型报道较少等问题。

发明内容

本发明的目的是克服现有的技术缺陷，提供一种菜籽油品质的近红外检测方法及应用，为菜籽油多项品质指标的检测提供一种快速、无损、环保、简便的方法。

一种菜籽油品质的近红外检测方法，所述方法包括以下步骤：

1)收集菜籽油样品；

2)对所述步骤1)收集的菜籽油样品进行理化指标测试，得到理化指标测试值，记为y_mj，其中：m是第m个指标，m＝1,2,3，…，13；当m＝1时记为棕榈酸含量，m＝2时记为亚油酸含量，m＝3时记为油酸含量，m＝4时记为亚麻酸含量，m＝5时记为硬脂酸含量，m＝6时记为花生一烯酸含量，m＝7时记为芥酸含量，m＝8时记为酸价，m＝9时记为过氧化值，m＝10时记为碘价，m＝11时记为色泽黄值，m＝12时记为色泽红值，m＝13时记为黏度；j为第j个样品，共n个样品，n≥90；

3)对所述步骤1)收集的菜籽油样品进行近红外扫描，得到吸光值为X_ij，其中i表示波长数值；

4)异常样本剔除和样品集划分：对所述步骤2)得到的菜籽油理化指标测试值y_mj和所述步骤3)得到的吸光值x_ij，采用学生残差法检测异常样本，剔除异常样本后的样品集采用Kennard-Stone法将全部样品集划分为校正集和验证集；

5)吸光值预处理：对所述步骤4)校正集吸光值进行预处理，得到校正集预处理吸光值A_ij；

6)筛选近红外光谱特征波长：采用逐步回归法对步骤2)得到的理化指标y_mj和步骤5)得到的校正集预处理吸光值A_ij进行分析，筛选得到菜籽油品质的近红外光谱特征波长；

7)建立菜籽油品质预测模型：采用多元线性回归法对步骤2)得到的理化指标y_mj和步骤5)得到的近红外光谱特征波长的预处理吸光值A_ij进行分析，建立菜籽油品质的预测模型为Z_m＝b_m+∑a_miX_i，其中Z_m为菜籽油品质的近红外检测值，X_i为A_ij中第j个样品的吸光值，b_m、a_mi为回归系数；

8)预测模型的验证：将所述步骤4)划分出的样品验证集按照所述步骤5)方法进行预处理，得到验证集预处理吸光值，代入所述步骤7)的预测模型，得到菜籽油品质指标预测值，与所述步骤2)该样品理化指标测试值进行比对，得到模型预测的精确度；

9)近红外指纹波长的确定：

完成菜籽油品质的检测。

所述步骤3)的近红外扫描温度为20-40℃，光程为3-6mm。所述步骤3)的菜籽油样品进行近红外扫描波长为：波长i＝1000nm，1001nm，1002nm，…，1799nm。

所述步骤5)预处理方法为一阶导数或二阶导数或多元散射校正或或净分析信号或正交信号校正或标准正态变量变换或去趋势校正或标准化或均值中心化中的一种或其组合。

所述步骤6)菜籽油品质的近红外光谱特征波长为：

棕榈酸：1097nm、1102nm、1260nm、1266nm、1286nm、1299nm、1355nm、1372nm、1428nm、1438nm、1450nm、1521nm、1534nm、1536nm、1619nm、1628nm、1629nm、1633nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

亚油酸：1040nm、1072nm、1077nm、1098nm、1108nm、1110nm、1125nm、1176nm、1181nm、1187nm、1194nm、1236nm、1334nm、1344nm、1395nm、1401nm、1442nm、1451nm、1457nm、1497nm、1608nm、1641nm、1643nm、1646nm、1657nm、1662nm、1668nm、1682nm、1687nm、1740nm、1780nm、1791nm、1794nm；

油酸：1024nm、1027nm、1028nm、1039nm、1065nm、1076nm、1092nm、1107nm、1115nm、1131nm、1183nm、1200nm、1203nm、1215nm、1218nm、1226nm、1227nm、1304nm、1325nm、1343nm、1376nm、1401nm、1409nm、1436nm、1455nm、1505nm、1554nm、1570nm、1622nm、1629nm、1635nm、1662nm、1663nm、1665nm、1672nm、1685nm、1698nm、1743nm、1788nm；

亚麻酸：1000nm、1039nm、1078nm、1084nm、1094nm、1146nm、1168nm、1175nm、1188nm、1192nm、1198nm、1282nm、1289nm、1341nm、1381nm、1414nm、1426nm、1433nm、1434nm、1451nm、1517nm、1575nm、1593nm、1594nm、1595nm、1598nm、1605nm、1606nm、1614nm、1637nm、1673nm、1697nm、1711nm、1712nm、1714nm、1726nm、1727nm、1729nm、1737nm、1743nm、1758nm、1771nm、1773nm；

硬脂酸：1008nm、1014nm、1024nm、1063nm、1078nm、1083nm、1088nm、1096nm、1117nm、1150nm、1164nm、1193nm、1238nm、1260nm、1264nm、1286nm、1287nm、1302nm、1339nm、1356nm、1357nm、1359nm、1389nm、1399nm、1402nm、1421nm、1442nm、1465nm、1481nm、1501nm、1518nm、1531nm、1532nm、1533nm、1542nm、1543nm、1585nm、1631nm、1651nm、1656nm、1669nm；

花生一烯酸：1026nm、1036nm、1041nm、1053nm、1056nm、1070nm、1078nm、1085nm、1091nm、1141nm、1181nm、1197nm、1205nm、1217nm、1241nm、1263nm、1268nm、1279nm、1289nm、1321nm、1348nm、1361nm、1385nm、1390nm、1394nm、1413nm、1433nm、1447nm、1460nm、1523nm、1545nm、1577nm、1623nm、1625nm、1634nm、1640nm、1671nm、1675nm、1738nm；

芥酸：1073nm、1105nm、1110nm、1123nm、1140nm、1145nmn、1155nm、1190nm、1201nm、1205nm、1229nm、1258nm、1273nm、1291nm、1302nm、1351nm、1374nm、1403nm、1408nm、1409nm、1433nm、1444nm、1458nm、1511nm、1527nm、1528nm、1534nm、1568nm、1601nm、1630nm、1631nm、1637nm、1650nm、1660nm、1661nm、1679nm、1687nm、1688nm、1690nm、1693nm、1695nm、1698nm、1699nm、1701nm、1739nm、1743nm、1791nm、1793nm；

酸价：1004nm、1033nm、1035nm、1051nm、1063nm、1079nm、1118nm、1143nm、1170nm、1180nm、1204nm、1212nm、1230nm、1250nm、1264nm、1268nm、1289nm、1308nm、1309nm、1316nm、1329nm、1378nm、1396nm、1410nm、1414nm、1477nm、1478nm、1497nm、1577nm、1617nm；

过氧化值：1006nm、1018nm、1038nm、1053nm、1070nm、1088nm、1110nm、1113nm、1127nm、1132nm、1139nm、1222nm、1232nm、1240nm、1277nm、1285nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1494nm、1559nm、1571nm、1627nm、1660nm、1682nm、1765nm、1774nm、1784nm、1787nm；

碘价：1007nm、1008nm、1012nm、1017nm、1020nm、1041nm、1042nm、1051nm、1066nm、1078nm、1100nm、1185nm、1190nm、1228nm、1243nm、1266nm、1283nm、1305nm、1311nm、1312nm、1346nm、1354nm、1368nm、1369nm、1374nm、1378nm、1382nm、1394nm、1398nm、1399nm、1423nm、1444nm、1445nm、1491nm、1510nm、1653nm、1654nm、1655nm、1671nm、1688nm、1698nm、1704nm、1705nm、1707nm、1708nm、1717nm、1721nm、1729nm、1733nm、1743nm、1746nm、1749nm、1763nm、1771nm、1782nm、1785nm、1794nm；

色泽黄值：1003nm、1005nm、1051nm、1098nm、1149nm、1157nm、1174nm、1179nm、1185nm、1199nm、1218nm、1295nm、1299nm、1357nm、1364nm、1372nm、1398nm、1419nm、1430nm、1479nm、1503nm、1517nm、1526nm、1527nm、1531nm、1532nm、1533nm、1537nm、1556nm、1608nm、1610nm、1643nm、1662nm、1680nm、1697nm、1737nm、1739nm、1767nm、1772nm、1788nm、1790nm、1799nm；

色泽红值：1031nm、1051nm、1071nm、1100nm、1105nm、1123nm、1125nm、1150nm、1161nm、1169nm、1179nm、1184nm、1218nm、1232nm、1294nm、1295nm、1300nm、1312nm、1357nm、1370nm、1389nm、1470nm、1492nm、1528nm、1529nm、1538nm、1551nm、1565nm、1590nm、1602nm、1639nm、1644nm、1657nm、1690nm、1694nm、1702nm、1703nm、1743nm、1744nm、1778nm、1790nm；

黏度：1015nm、1023nm、1025nm、1117nm、1129nm、1140nm、1163nm、1179nm、1210nm、1249nm、1261nm、1359nm、1365nm、1370nm、1379nm、1382nm、1414nm、1426nm、1468nm、1493nm、1508nm、1560nm、1572nm、1582nm、1590nm、1591nm、1604nm、1628nm、1696nm、1792nm；

上述特征波长允许有±2nm的偏差。

所述步骤7)中预测模型的评价指标为决定系数R²、校正集标准差RMSEC和验证集标准差RMSEP。

所述步骤7)中建立菜籽油品质预测模型：

棕榈酸：

Z₁＝7.122-554.923X₁₀₉₇+559.861X₁₁₀₂+1154.024X₁₂₆₀-1717.392X₁₂₆₆+1001.188X₁₂₈₆

-520.795X₁₂₉₉-400.896X₁₃₅₅+485.139X₁₃₇₂-540.567X₁₄₂₈+1052.793X₁₄₃₈-374.591X₁₄₅₀

+418.760X₁₅₂₁+501.261X₁₅₃₄-1063.741X₁₅₃₆+386.163X₁₆₁₉-2644.180X₁₆₂₈+2391.044X₁₆₂₉

-666.401X₁₆₃₃+729.704X₁₆₃₉-765.265X₁₆₇₀+1022.987X₁₆₇₉-969.315X₁₆₈₇+482.195X₁₆₉₃-

316.015X₁₇₇₇+277.649X₁₇₉₁；

亚油酸：

Z₂＝12.881+1494.756X₁₀₄₀+1188.261X₁₀₇₂-2106.146X₁₀₇₇+2004.052X₁₀₉₈-3264.258X₁₁₀₈-

2164.257X₁₁₁₀+2926.439X₁₁₂₅-1919.937X₁₁₇₆+3415.763X₁₁₈₁-1835.691X₁₁₈₇

+2643.582X₁₁₉₄-1986.239X₁₂₃₆+2298.319X₁₃₃₄-1712.230X₁₃₄₄-2152.860X₁₃₉₅

+1301.989X₁₄₀₁+669.145X₁₄₄₂-2266.494X₁₄₅₁+1861.268X₁₄₅₇-1366.906X₁₄₉₇

-2751.327X₁₆₀₈+1042.784X₁₆₄₁+7273.727X₁₆₄₃-4483.278X₁₆₄₆-3243.621X₁₆₅₇

+6191.148X₁₆₆₂-3773.298X₁₆₆₈+2006.371X₁₆₈₂-1839.867X₁₆₈₇+883.481X₁₇₄₀

-1004.192X₁₇₈₀+1035.786X₁₇₉₁-1120.310X₁₇₉₄；

油酸：

Z₃＝131.518+54.648X₁₀₂₄-145.218X₁₀₂₇+238.569X₁₀₂₈-124.154X₁₀₃₉-38.688X₁₀₆₅

-187.688X₁₀₇₆-139.253X₁₀₉₂+222.844X₁₁₀₇+215.092X₁₁₁₅+63.101X₁₁₃₁+86.268X₁₁₈₃

-134.256X₁₂₀₀-78.406X₁₂₀₃-74.438X₁₂₁₅+224.804X₁₂₁₈+267.659X₁₂₂₆-134.32748X₁₂₂₇

-212.718X₁₃₀₄-88.665X₁₃₂₅-97.920X₁₃₄₃+100.588X₁₃₇₆-120.263X₁₄₀₁+137.713X₁₄₀₉

-219.328X₁₄₃₆+145.883X₁₄₅₅-189.689X₁₅₀₅+187.495X₁₅₅₄-299.930X₁₅₇₀+201.052X₁₆₂₂

+168.477X₁₆₂₉-186.007X₁₆₃₅+446.389X₁₆₆₂-323.377X₁₆₆₃+367.357X₁₆₆₅-135.115X₁₆₇₂

+72.825X₁₆₈₅-50.836X₁₆₉₈+46.310X₁₇₄₃-53.128X₁₇₈₈；

亚麻酸：

Z₄＝7506.006+1591.699X₁₀₀₀-4004.831X₁₀₃₉-2633.555X₁₀₇₈+3094.150X₁₀₈₄

-2992.399X₁₀₉₄-3152.559X₁₁₄₆+1920.268X₁₁₆₈-2117.618X₁₁₇₅+3385.247X₁₁₈₈

-4188.517X₁₁₉₂-1635.062X₁₁₉₈-3958.504X₁₂₈₂+2138.468X₁₂₈₉-1634.357X₁₃₄₁-927.100X₁₃₈₁

-2385.798X₁₄₁₄-1935.714X₁₄₂₆+2986.193X₁₄₃₃+366.652X₁₄₃₄-1159.740X₁₄₅₁

-3490.094X₁₅₁₇-1700.527X₁₅₇₅-9384.497X₁₅₉₃+14487X₁₅₉₄-10710X₁₅₉₅-5118.240X₁₅₉₈

+12228X₁₆₀₅-5748.613X₁₆₀₆+4074.247X₁₆₁₄-2355.427X₁₆₃₇-2148.543X₁₆₇₃-959.665X₁₆₉₇

-921.807X₁₇₁₁+1177.955X₁₇₁₂-417.392X₁₇₁₄-882.869X₁₇₂₆+1267.076X₁₇₂₇-877.131X₁₇₂₉

-423.845X₁₇₃₇+822.571X₁₇₄₃-910.270X₁₇₅₈+952.185X₁₇₇₁-1287.709X₁₇₇₃；

硬脂酸：

Z₅＝-11.140+3255.650X₁₀₀₈+3240.414X₁₀₁₄-2422.871X₁₀₂₄+2690.775X₁₀₆₃

-2097.645X₁₀₇₈-2325.170X₁₀₈₃-6746.814X₁₀₈₈-827.644X₁₀₉₆-2125.511X₁₁₁₇-5597.104X₁₁₅₀

+4606.841X₁₁₆₄+4574.228X₁₁₉₃+8693.499X₁₂₃₈-8975.580X₁₂₆₀+2426.146X₁₂₆₄-3632.476X

₁₂₈₆-53.705X₁₂₈₇+1334.176X₁₃₀₂-5409.880X₁₃₃₉+7075.527X₁₃₅₆-9719.556X₁₃₅₇+11419X₁₃₅₉

+3817.161X₁₃₈₉+2471.477X₁₃₉₉+5158.646X₁₄₀₂+5644.206X₁₄₂₁+1723.091X₁₄₄₂

+6683.316X₁₄₆₅+2456.902X₁₄₈₁+5891.291X₁₅₀₁-757.199X₁₅₁₈-2085.148X₁₅₃₁

-5314.264X₁₅₃₂-6929.488X₁₅₃₃+13.519X₁₅₄₂+6958.481X₁₅₄₃-2976.579X₁₅₈₅+4905.371X₁₆₃₁

+1954.005X₁₆₅₁+1529.944X₁₆₅₆-1580.192X₁₆₆₉；

花生一烯酸：

Z₆＝2.815+1284.950X₁₀₂₆-1461.413X₁₀₃₆-404.856X₁₀₄₁-260.283X₁₀₅₃+805.799X₁₀₅₆

-1159.894X₁₀₇₀+1577.242X₁₀₇₈-808.352X₁₀₈₅-1020.663X₁₀₉₁+389.510X₁₁₄₁

-1156.822X₁₁₈₁-1502.081X₁₁₉₇+985.943X₁₂₀₅+342.352X₁₂₁₇+1878.121X₁₂₄₁-1334.758X₁₂₆₃

+475.320X₁₂₆₈-1755.228X₁₂₇₉+1593.834X₁₂₈₉-672.329X₁₃₂₁+631.634X₁₃₄₈-644.578X₁₃₆₁

+1496.912X₁₃₈₅+694.086X₁₃₉₀-812.705X₁₃₉₄-446.055X₁₄₁₃+501.298X₁₄₃₃+458.211X₁₄₄₇

-530.111X₁₄₆₀+866.661X₁₅₂₃+597.164X₁₅₄₅-1612.002X₁₅₇₇+1330.369X₁₆₂₃+1429.233X₁₆₂₅

-2031.183X₁₆₃₄-1173.266X₁₆₄₀+445.451X₁₆₇₁+895.088X₁₆₇₅-298.929X₁₇₃₈；

芥酸：

Z₇＝10.284+4659.521X₁₀₇₃-2923.995X₁₁₀₅+1345.165X₁₁₁₀-6375.893X₁₁₂₃+4111.558X₁₁₄₀

+5199.399X₁₁₄₅+4806.482X₁₁₅₅-6702.304X₁₁₉₀+7739.103X₁₂₀₁+4164.891X₁₂₀₅

-1952.258X₁₂₂₉-1801.109X₁₂₅₈-10555X₁₂₇₃-4023.318X₁₂₉₁-6897.226X₁₃₀₂-3135.779X₁₃₅₁

-2720.412X₁₃₇₄+7058.360X₁₄₀₃-4950.972X₁₄₀₈-3924.462X₁₄₀₉-7708.661X₁₄₃₃

-2204.383X₁₄₄₄-3509.622X₁₄₅₈-5880.686X₁₅₁₁-4326.789X₁₅₂₇+6809.207X₁₅₂₈+12682X₁₅₃₄

+4803.327X₁₅₆₈+7116.305X₁₆₀₁+3908.016X₁₆₃₀+5744.782X₁₆₃₁-3704.439X₁₆₃₇

+6322.864X₁₆₅₀-5132.705X₁₆₆₀-6511.577X₁₆₆₁+4127.098X₁₆₇₉+2712.062X₁₆₈₇

-5831.724X₁₆₈₈-1914.506X₁₆₉₀+1515.624X₁₆₉₃+1262.943X₁₆₉₅-2057.092X₁₆₉₈

+3295.167X₁₆₉₉-1365.754X₁₇₀₁-919.912X₁₇₃₉-1699.429X₁₇₄₃+2978.033X₁₇₉₁

+2350.872X₁₇₉₃；

酸价：

Z₈＝0.119-2538.959X₁₀₀₄-3100.461X₁₀₃₃+4886.971X₁₀₃₅+2897.828X₁₀₅₁-2006.278X₁₀₆₃

+2760.098X₁₀₇₉-2640.380X₁₁₁₈+2648.164X₁₁₄₃-2753.983X₁₁₇₀+5003.504X₁₁₈₀

-1000.648X₁₂₀₄+1400.125X₁₂₁₂-3286.331X₁₂₃₀+4502.311X₁₂₅₀-2293.994X₁₂₆₄

-5782.938X₁₂₆₈+2145.841X₁₂₈₉+10122X₁₃₀₈-11554X₁₃₀₉+2774.508X₁₃₁₆+567.939X₁₃₂₉

-1925.473X₁₃₇₈+1430.852X₁₃₉₆-3064.104X₁₄₁₀+2059.191X₁₄₁₄+21246X₁₄₇₇-20903X₁₄₇₈

-1061.150X₁₄₉₇+3780.992X₁₅₇₇-4940.442X₁₆₁₇；

过氧化值：

Z₉＝-1.095+7893.961X₁₀₀₆+8674.625X₁₀₁₈-2189.742X₁₀₃₈+2747.823X₁₀₅₃-4545.457X₁₀₇₀

-5891.021X₁₀₈₈-11945X₁₁₁₀+7041.339X₁₁₁₃-6603.093X₁₁₂₇+8646.126X₁₁₃₂+5028.386X₁₁₃₉

+7597.665X₁₂₂₂-11089X₁₂₃₂-1108.292X₁₂₄₀+9275.776X₁₂₇₇+8261.865X₁₂₈₅-12310X₁₃₆₉

+13064X₁₃₈₂-5946.595X₁₄₀₅+7664.969X₁₄₁₄-11912X₁₄₃₅+9507.563X₁₄₆₀+7903.013X₁₄₉₄

-13248X₁₅₅₉-7108.786X₁₅₇₁-2923.063X₁₆₂₇-3778.807X₁₆₆₀-4053.265X₁₆₈₂-3191.812X₁₇₆₅

+4108.125X₁₇₇₄-5542.854X₁₇₈₄+7032.594X₁₇₈₇；

碘价：

Z₁₀＝70.379+5655.836X₁₀₀₇+3835.644X₁₀₀₈+2255.332X₁₀₁₂+8315.045X₁₀₁₇+5967.677X₁₀₂₀

+5946.963X₁₀₄₁+4961.757X₁₀₄₂+8824.426X₁₀₅₁+9546.300X₁₀₆₆+2556.441X₁₀₇₈

+10214X₁₁₀₀-3767.757X₁₁₈₅-6114.234X₁₁₉₀-7551.222X₁₂₂₈-12332X₁₂₄₃-2271.972X₁₂₆₆

-8074.840X₁₂₈₃-1175.525X₁₃₀₅+3644.068X₁₃₁₁-12947X₁₃₁₂-6439.407X₁₃₄₆-7045.824X₁₃₅₄

-14294X₁₃₆₈+8944.760X₁₃₆₉-8323.329X₁₃₇₄-6197.887X₁₃₇₈-8325.039X₁₃₈₂-9125.409X₁₃₉₄

-3157.877X₁₃₉₈-2646.257X₁₃₉₉+5074.726X₁₄₂₃-4409.695X₁₄₄₄+16935X₁₄₄₅-7959.432X₁₄₉₁

-9135.176X₁₅₁₀+3952.970X₁₆₅₃+9651.345X₁₆₅₄+2212.403X₁₆₅₅+4346.486X₁₆₇₁

-4980.917X₁₆₈₈+5301.666X₁₆₉₈+4013.108X₁₇₀₄+2637.914X₁₇₀₅-3167.676X₁₇₀₇

+5617.745X₁₇₀₈-3577.202X₁₇₁₇-4002.850X₁₇₂₁+5140.303X₁₇₂₉+4737.340X₁₇₃₃

-2469.897X₁₇₄₃-3276.167X₁₇₄₆-959.039X₁₇₄₉+2539.913X₁₇₆₃-2006.402X₁₇₇₁-3899.264X₁₇₈₂

-5321.969X₁₇₈₅+3019.113X₁₇₉₄；

色泽黄值：

Z₁₁＝7.200+15691X₁₀₀₃+14907X₁₀₀₅-12573X₁₀₅₁+50407X₁₀₉₈+38618X₁₁₄₉-22311X₁₁₅₇

+21980X₁₁₇₄-46337X₁₁₇₉+15773X₁₁₈₅+21827X₁₁₉₉+20574X₁₂₁₈+49588X₁₂₉₅-52794X₁₂₉₉

+46525X₁₃₅₇-23933X₁₃₆₄+11818X₁₃₇₂-10222X₁₃₉₈+25684X₁₄₁₉-25402X₁₄₃₀+32483X₁₄₇₉

-42643X₁₅₀₃-68673X₁₅₁₇+74279X₁₅₂₆-30897X₁₅₂₇+120133X₁₅₃₁-149057X₁₅₃₂+167820X₁₅₃₃

-88396X₁₅₃₇-28409X₁₅₅₆+68851X₁₆₀₈-30206X₁₆₁₀-46112X₁₆₄₃+11539X₁₆₆₂+10951X₁₆₈₀

-10465X₁₆₉₇+7317.372X₁₇₃₇-14058X₁₇₃₉-7434.215X₁₇₆₇+13105X₁₇₇₂+12795X₁₇₈₈

-27958X₁₇₉₀-6001.537X₁₇₉₉；

色泽红值：

Z₁₂＝0.891-1868.697X₁₀₃₁-2572.488X₁₀₅₁-288.888X₁₀₇₁+5057.192X₁₁₀₀-3072.222X₁₁₀₅

+1257.980X₁₁₂₃-3508.014X₁₁₂₅+3855.554X₁₁₅₀-3758.369X₁₁₆₁+1483.714X₁₁₆₉

-2301.038X₁₁₇₉+2555.615X₁₁₈₄+2017.671X₁₂₁₈-2454.693X₁₂₃₂+4498.395X₁₂₉₄

+2850.912X₁₂₉₅-4340.720X₁₃₀₀-2012.020X₁₃₁₂+3303.660X₁₃₅₇-1920.040X₁₃₇₀

+707.359X₁₃₈₉+3523.153X₁₄₇₀-3226.896X₁₄₉₂+1275.462X₁₅₂₈+5537.216X₁₅₂₉

-4646.654X₁₅₃₈+1219.908X₁₅₅₁+4206.180X₁₅₆₅-6255.102X₁₅₉₀+5286.438X₁₆₀₂

+5131.157X₁₆₃₉-5381.917X₁₆₄₄+3517.545X₁₆₅₇+2421.794X₁₆₉₀-2134.497X₁₆₉₄

-1276.800X₁₇₀₂+1863.362X₁₇₀₃-2434.757X₁₇₄₃+2167.080X₁₇₄₄+1611.073X₁₇₇₈

-851.937X₁₇₉₀；

黏度：

Z₁₃＝53.980-21386X₁₀₁₅-26783X₁₀₂₃+37751X₁₀₂₅-26813X₁₁₁₇-35111X₁₁₂₉+25465X₁₁₄₀

-32237X₁₁₆₃-21558X₁₁₇₉-21465X₁₂₁₀-12785X₁₂₄₉-61270X₁₂₆₁+57634X₁₃₅₉-60468X₁₃₆₅

-65030X₁₃₇₀-9524.085X₁₃₇₉+26949X₁₃₈₂-21009X₁₄₁₄-11847X₁₄₂₆-34364X₁₄₆₈+28553X₁₄₉₃

-59846X₁₅₀₈-66740X₁₅₆₀+34126X₁₅₇₂+51567X₁₅₈₂-69656X₁₅₉₀-33573X₁₅₉₁+74660X₁₆₀₄

-41775X₁₆₂₈-15134X₁₆₉₆-10211X₁₇₉₂。

所述步骤9)指纹波长确定方法为对所述步骤7)预测模型的每个特征波长对应的回归系数进行显著性分析，极显著的特征波长即为该菜籽油理化指标的近红外光谱指纹波长。

所述步骤9)菜籽油品质的近红外光谱指纹波长为：

亚油酸：1125nm、1194nm、1236nm、1657nm、1662nm、1668nm、1687nm、1740nm；

油酸：1107nm、1218nm、1409nm、1436nm、1698nm、1743nm；

亚麻酸：1039nm、1673nm、1697nm、1727nm、1729nm、1758nm；

棕榈酸：1260nm、1266nm、1286nm、1372nm、1428nm、1438nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

硬脂酸：1150nm、1164nm、1193nm、1238nm、1260nm、1359nm、1389nm、1402nm、1421nm、1465nm；

花生一烯酸：1197nm、1241nm、1634nm、1738nm；

芥酸：1190nm、1201nm、1205nm、1403nm、1433nm、1534nm、1650nm、1679nm、1743nm、1791nm；

酸价：1035nm、1170nm、1180nm、1414nm；

过氧化值：1018nm、1222nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1682nm、1765nm、1787nm；

碘价：1100nm、1688nm、1698nm、1708nm、1717nm、1721nm、1729nm、1733nm、1746nm、1763nm；

色泽黄值：1098nm、1179nm、1479nm、1643nm；

色泽红值：1100nm、1357nm、1657nm、1690nm、1694nm；

黏度：1163nm、1210nm、1604nm、1628nm、1696nm；

上述波长允许有±2nm的偏差。

所述的菜籽油品质的近红外检测方法在玉米油或橄榄油或山茶油或棕榈油或芥花子油或葵花子油或大豆油或芝麻油检测中的应用。

本发明提供的一种菜籽油品质的近红外检测方法及应用，有益效果如下：

1、本发明确定了检测最适宜的条件，适宜的光程和扫描温度，保证检测结果的准确性和稳定性；

2、确定了菜籽油品质的近红外特征波长，且特征波长数量少，运用特征波长来建立模型的分析方法可提高模型的精度；基于特征波长建立的菜籽油品质指标的近红外预测模型的方法，在相同预测精度的情况下，本发明预测速度快；

3、实现了菜籽油品质的快速检测，常规的菜籽油品质检测，耗时长，而近红外检测只需要几分钟，且可以一次性测出多种指标，分析更加全面；

4、实现了菜籽油品质的无损检测，常规的菜籽油品质检测需要破坏油脂本身的性质，检测后不可回收利用，而近红外检测后，油脂仍然可以正常使用；

5、对菜籽油品质分析较全面，建立了菜籽油的7种脂肪酸组成、酸价、过氧化值、碘价、色泽、黏度的预测模型和测试方法。

附图说明

图1：是本发明不同光程条件下采集同一样品的原始光谱图；

图2：是本发明不同光程条件下采集同一样品的标准偏差光谱图；

图3：是本发明不同扫描温度条件下采集同一样品的原始光谱图；

图4：是本发明不同扫描温度条件下采集同一样品的标准偏差光谱图；

图5：是本发明实施例1-2中菜籽油样品原始近红外光谱图；

图6：是本发明实施例1中异常样本剔除图；

其中：1为光程10mm的原始光谱数据，2为光程5mm的原始光谱数据，3为光程2mm的原始光谱数据，4为光程10mm的标准偏差光谱数据，5为光程5mm的标准偏差光谱数据，6为光程2mm的标准偏差光谱数据，7为扫描温度10℃的标准偏差光谱数据，8为扫描温度20℃的标准偏差光谱数据，9为扫描温度30℃的标准偏差光谱数据，10为扫描温度20℃的标准偏差光谱数据，11为扫描温度50℃的标准偏差光谱数据。

具体实施方式

实验材料与方法

近红外光谱采集温度的确定

菜籽油样品摇匀，置于光程为5mm的石英石样品池内，扫描次数10次，固定波长范围为1000～1799nm，仪器带宽为1nm，光谱重复性优于0.2nm，置于不同光程下(2mm，5mm，10mm)条件下，采集同一菜籽油样的近红外光谱图5次，得到不同光程条件下采集同一样品的原始光谱图如图1所示，光程为10mm时菜籽油吸光值最大，光程为2mm，吸光值最小，但从图1看不出来，对原始光谱数据进行标准差处理，得到不同光程条件下采集同一样品的标准偏差光谱图2，可见光程为10mm时，各波长点下吸光值的相对标准偏差最大，说明菜籽油光谱最不稳定。而光程为5mm时，各波长点下吸光值的相对标准偏差最小，菜籽油光谱最稳定。因此选择菜籽油近红外扫描光程为5mm。

在不同温度(10℃、20℃、30℃、40℃、50℃)条件下，采集同一菜籽油样的近红外原始光谱图如图3所示。由图3可知，不同扫描温度条件下采集的近红外光谱图重叠严重，看不出来哪种条件下光谱图最为稳定。因此经过同样的处理，得到图4不同扫描温度条件下同一菜籽油样品的标准偏差光谱图。由图4可见扫描温度在20℃、30℃、40℃时光谱图稳定性的差别不大，因此选择扫描温度为20-40℃。

菜籽油理化指标检测方法

1)脂肪酸组成：GB/T17376-2008、GB/T17377-2008测定。

2)酸价、过氧化值：GB/T5009.37-2003《食用植物油卫生标准的分析方法》测定。

3)碘值：GB/T5532-2008《动植物油脂碘值的测定》测定菜籽油样品碘值。

4)色泽：GB/T22460-2008《动植物油脂罗维朋色泽的测定方法》测定。

5)黏度：NDJ-5S型数字旋转黏度计测定。

模型评价指标

以决定系数R²、定标标准差RMSEC、校验标准差RMSEP来评价模型的定标效果和预测能力。

$R^2=1-\frac{\underset{i=1}{\overset{n}{\Σ}}{(y_i-{\hat{y}}_i)}^2}{\underset{i=1}{\overset{n}{\Σ}}{(y_i-{\stackrel{\‾}{y}}_i)}^2}$

式中，n、m分别为定标集和验证集样品数，y_i为第i个样品的某一品质指标的理化检测值，为第i个样品的品质指标近红外预测值，为样品检测值的平均值。R²越接近于1，说明回归效果显著，RMSEC和RMSEP越接近于0，说明模型具有很好的稳定性和预测能力。

本发明利用软件SAS8.1进数据分析。

实施例1：菜籽油棕榈酸的近红外检测

1)菜籽油样品：共125个菜籽油样品为不同产地、不同级别的菜籽纯油样品，另外还有20个菜籽调和油样品，共计145个；

2)按照GB/T17377-2008测定菜籽油棕榈酸含量，菜籽油样品棕榈酸含量的基本统计数据见表1；

表1菜籽油棕榈酸含量的基本统计数据(％)

3)将菜籽油样品摇匀，置于光程为5mm的石英石样品池内，盖好盖子，用近红外光谱仪扫描。扫描前需开机预热30min，预热后调试扫描次数为10次，检测室温度为40℃，固定波长范围为1000～1799nm，扫描间隔为1nm，仪器带宽为1nm，信噪比优于2000：1，光谱重复性优于0.2nm等参数后，将整个样品池置于检测室内并关闭检测室，再等3min左右直至检测室温度达到40℃时，开始采集菜籽油样的近红外光谱图，原始光谱见附图5；

4)异常样品剔除和样品集划分：将145个菜籽油样品品质检测值及近红外吸光值采用学生残差法检测异常样本，设定学生残差临界值为2.5，学生残差大于2.5的样品被认定为异常样品，棕榈酸异常样品剔除2个，见附图6。异常样本剔除的影响见实施例4。剩下的样品用Kennard-Stone(K-S)法划分为定标集样品114个及验证集样品29个；规定样品集中的80％为校正集，用于构建预测模型，20％为验证集，用来检验所建模型的准确定与可靠性；

5)近红外光谱预处理：对所述步骤4)校正集吸光值x_ij先用S-G卷积平滑进行预处理，设置窗口参数为9，再用正交信号校正的方法进行预处理，得到预处理吸光值A_ij，预处理方法的影响见实施例5；

6)筛选近红外光谱特征波长：采用逐步回归法对步骤2)得到的理化指标y_mj和步骤5)得到的校正集预处理吸光值A_ij进行分析，筛选出的棕榈酸特征波长为：1097nm、1102nm、1260nm、1266nm、1286nm、1299nm、1355nm、1372nm、1428nm、1438nm、1450nm、1521nm、1534nm、1536nm、1619nm、1628nm、1629nm、1633nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

7)预测模型的构建：采用多元线性回归法对步骤2)得到的理化指标y_mj和步骤5)得到的近红外光谱特征波长的预处理吸光值A_ij进行分析，建立菜籽油品质的预测模型为Z₁＝b₁+∑a_1iXi，回归系数及显著性结果见表2；Z₁为菜籽油棕榈酸近红外预测值，b₁为回归常数项，此时b₁＝7.122，a_1i为各特征波长的回归系数，x_i为特征波长的吸光值经过消噪和预处理后的数值，i为特征波长，拟合出的模型R²为0.9995，RMSEC为0.0751，可看出模型的拟合精度非常高。

表2菜籽油棕榈酸的近红外模型显著性分析

8)预测模型的验证：将步骤4)中划分出的验证集29个样品验证集按照所述步骤5)方法进行预处理，得到验证集预处理吸光值，代入到所述步骤7)的预测模型，得到29个样品的棕榈酸近红外预测值，与所述步骤2)该样品棕榈酸理化测试值进行比对，结果见表3。从表3的数据得到RMSEP为0.5980，对菜籽油棕榈酸验证集样品的预测值与理化测量值进行显著性分析，给定显著水平0.05，棕榈酸预测值和理化测量值之间不存在显著性差异，表明模型预测的精确度较高。

表3菜籽油棕榈酸近红外模型的检测结果

9)指纹波长确定：对所述步骤7)预测模型的每个特征波长对应的回归系数进行显著性分析结果见表2，取回归系数极显著的特征波长为棕榈酸的近红外光谱指纹波长，为1260nm、1266nm、1286nm、1372nm、1428nm、1438nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

完成菜籽油棕榈酸的检测。

实施例2：菜籽油油酸、亚油酸、亚麻酸、硬脂酸、花生一烯酸、芥酸、酸价、过氧化值、碘值、色泽黄值、红值、黏度的近红外检测

1)用实施例1步骤1)的方法检测菜籽油亚油酸、亚麻酸、油酸、硬脂酸、花生一烯酸、芥酸的含量，测定菜籽油酸价、过氧化值、碘值、色泽黄值、红值、黏度指标，见表4。

表4菜籽油脂肪酸基本统计数据(mg/g)

3)按照实施例1步骤3)的方法对样品进行近红外光谱扫描；

4)按照实施例1步骤4)的方法对样品进行异常样品剔除和样品集划分；

5)近红外光谱预处理：对所述步骤4)校正集吸光值x_ij进行预处理得到预处理吸光值A_ij，

其中预处理方法：硬脂酸和芥酸采用的是一阶导数加正交信号校正，油酸采用的是一阶导数加标准正态变量变换，亚麻酸采用的是多元散射校正；酸价、过氧化值和碘价采用正交信号校正；色泽黄值、色泽红值、黏度采用净分析信号(NAS)；

6)采用实施例1步骤6)的方法分别确定了各指标的特征波长：

亚油酸：1040nm、1072nm、1077nm、1098nm、1108nm、1110nm、1125nm、1176nm、1181nm、1187nm、1194nm、1236nm、1334nm、1344nm、1395nm、1401nm、1442nm、1451nm、1457nm、1497nm、1608nm、1641nm、1643nm、1646nm、1657nm、1662nm、1668nm、1682nm、1687nm、1740nm、1780nm、1791nm、1794nm；

油酸：1024nm、1027nm、1028nm、1039nm、1065nm、1076nm、1092nm、1107nm、1115nm、1131nm、1183nm、1200nm、1203nm、1215nm、1218nm、1226nm、1227nm、1304nm、1325nm、1343nm、1376nm、1401nm、1409nm、1436nm、1455nm、1505nm、1554nm、1570nm、1622nm、1629nm、1635nm、1662nm、1663nm、1665nm、1672nm、1685nm、1698nm、1743nm、1788nm；

亚麻酸：1000nm、1039nm、1078nm、1084nm、1094nm、1146nm、1168nm、1175nm、1188nm、1192nm、1198nm、1282nm、1289nm、1341nm、1381nm、1414nm、1426nm、1433nm、1434nm、1451nm、1517nm、1575nm、1593nm、1594nm、1595nm、1598nm、1605nm、1606nm、1614nm、1637nm、1673nm、1697nm、1711nm、1712nm、1714nm、1726nm、1727nm、1729nm、1737nm、1743nm、1758nm、1771nm、1773nm；

硬脂酸：1008nm、1014nm、1024nm、1063nm、1078nm、1083nm、1088nm、1096nm、1117nm、1150nm、1164nm、1193nm、1238nm、1260nm、1264nm、1286nm、1287nm、1302nm、1339nm、1356nm、1357nm、1359nm、1389nm、1399nm、1402nm、1421nm、1442nm、1465nm、1481nm、1501nm、1518nm、1531nm、1532nm、1533nm、1542nm、1543nm、1585nm、1631nm、1651nm、1656nm、1669nm；

花生一烯酸：1026nm、1036nm、1041nm、1053nm、1056nm、1070nm、1078nm、1085nm、1091nm、1141nm、1181nm、1197nm、1205nm、1217nm、1241nm、1263nm、1268nm、1279nm、1289nm、1321nm、1348nm、1361nm、1385nm、1390nm、1394nm、1413nm、1433nm、1447nm、1460nm、1523nm、1545nm、1577nm、1623nm、1625nm、1634nm、1640nm、1671nm、1675nm、1738nm；

芥酸：1073nm、1105nm、1110nm、1123nm、1140nm、1145nmn、1155nm、1190nm、1201nm、1205nm、1229nm、1258nm、1273nm、1291nm、1302nm、1351nm、1374nm、1403nm、1408nm、1409nm、1433nm、1444nm、1458nm、1511nm、1527nm、1528nm、1534nm、1568nm、1601nm、1630nm、1631nm、1637nm、1650nm、1660nm、1661nm、1679nm、1687nm、1688nm、1690nm、1693nm、1695nm、1698nm、1699nm、1701nm、1739nm、1743nm、1791nm、1793nm；

酸价：1004nm、1033nm、1035nm、1051nm、1063nm、1079nm、1118nm、1143nm、1170nm、1180nm、1204nm、1212nm、1230nm、1250nm、1264nm、1268nm、1289nm、1308nm、1309nm、1316nm、1329nm、1378nm、1396nm、1410nm、1414nm、1477nm、1478nm、1497nm、1577nm、1617nm；

过氧化值：1006nm、1018nm、1038nm、1053nm、1070nm、1088nm、1110nm、1113nm、1127nm、1132nm、1139nm、1222nm、1232nm、1240nm、1277nm、1285nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1494nm、1559nm、1571nm、1627nm、1660nm、1682nm、1765nm、1774nm、1784nm、1787nm；

碘价：1007nm、1008nm、1012nm、1017nm、1020nm、1041nm、1042nm、1051nm、1066nm、1078nm、1100nm、1185nm、1190nm、1228nm、1243nm、1266nm、1283nm、1305nm、1311nm、1312nm、1346nm、1354nm、1368nm、1369nm、1374nm、1378nm、1382nm、1394nm、1398nm、1399nm、1423nm、1444nm、1445nm、1491nm、1510nm、1653nm、1654nm、1655nm、1671nm、1688nm、1698nm、1704nm、1705nm、1707nm、1708nm、1717nm、1721nm、1729nm、1733nm、1743nm、1746nm、1749nm、1763nm、1771nm、1782nm、1785nm、1794nm；

色泽黄值：1003nm、1005nm、1051nm、1098nm、1149nm、1157nm、1174nm、1179nm、1185nm、1199nm、1218nm、1295nm、1299nm、1357nm、1364nm、1372nm、1398nm、1419nm、1430nm、1479nm、1503nm、1517nm、1526nm、1527nm、1531nm、1532nm、1533nm、1537nm、1556nm、1608nm、1610nm、1643nm、1662nm、1680nm、1697nm、1737nm、1739nm、1767nm、1772nm、1788nm、1790nm、1799nm；

色泽红值：1031nm、1051nm、1071nm、1100nm、1105nm、1123nm、1125nm、1150nm、1161nm、1169nm、1179nm、1184nm、1218nm、1232nm、1294nm、1295nm、1300nm、1312nm、1357nm、1370nm、1389nm、1470nm、1492nm、1528nm、1529nm、1538nm、1551nm、1565nm、1590nm、1602nm、1639nm、1644nm、1657nm、1690nm、1694nm、1702nm、1703nm、1743nm、1744nm、1778nm、1790nm；

黏度：1015nm、1023nm、1025nm、1117nm、1129nm、1140nm、1163nm、1179nm、1210nm、1249nm、1261nm、1359nm、1365nm、1370nm、1379nm、1382nm、1414nm、1426nm、1468nm、1493nm、1508nm、1560nm、1572nm、1582nm、1590nm、1591nm、1604nm、1628nm、1696nm、1792nm；

7)采用实施例1步骤7)的方法分别建立了各指标的预测模型如下所示，从各模型R²及RMSEC可看出各模型的拟合精度均较高；

亚油酸(模型R²为0.9997，RMSEC为0.2107)：

Z₂＝12.881+1494.756X₁₀₄₀+1188.261X₁₀₇₂-2106.146X₁₀₇₇+2004.052X₁₀₉₈-3264.258X₁₁₀₈-

2164.257X₁₁₁₀+2926.439X₁₁₂₅-1919.937X₁₁₇₆+3415.763X₁₁₈₁-1835.691X₁₁₈₇

+2643.582X₁₁₉₄-1986.239X₁₂₃₆+2298.319X₁₃₃₄-1712.230X₁₃₄₄-2152.860X₁₃₉₅

+1301.989X₁₄₀₁+669.145X₁₄₄₂-2266.494X₁₄₅₁+1861.268X₁₄₅₇-1366.906X₁₄₉₇

-2751.327X₁₆₀₈+1042.784X₁₆₄₁+7273.727X₁₆₄₃-4483.278X₁₆₄₆-3243.621X₁₆₅₇

+6191.148X₁₆₆₂-3773.298X₁₆₆₈+2006.371X₁₆₈₂-1839.867X₁₆₈₇+883.481X₁₇₄₀

-1004.192X₁₇₈₀+1035.786X₁₇₉₁-1120.310X₁₇₉₄；

油酸(模型R²为0.9993，RMSEC为0.2516)：

Z₃＝131.518+54.648X₁₀₂₄-145.218X₁₀₂₇+238.569X₁₀₂₈-124.154X₁₀₃₉-38.688X₁₀₆₅

-187.688X₁₀₇₆-139.253X₁₀₉₂+222.844X₁₁₀₇+215.092X₁₁₁₅+63.101X₁₁₃₁+86.268X₁₁₈₃

-134.256X₁₂₀₀-78.406X₁₂₀₃-74.438X₁₂₁₅+224.804X₁₂₁₈+267.659X₁₂₂₆-134.32748X₁₂₂₇

-212.718X₁₃₀₄-88.665X₁₃₂₅-97.920X₁₃₄₃+100.588X₁₃₇₆-120.263X₁₄₀₁+137.713X₁₄₀₉

-219.328X₁₄₃₆+145.883X₁₄₅₅-189.689X₁₅₀₅+187.495X₁₅₅₄-299.930X₁₅₇₀+201.052X₁₆₂₂

+168.477X₁₆₂₉-186.007X₁₆₃₅+446.389X₁₆₆₂-323.377X₁₆₆₃+367.357X₁₆₆₅-135.115X₁₆₇₂

+72.825X₁₆₈₅-50.836X₁₆₉₈+46.310X₁₇₄₃-53.128X₁₇₈₈；

亚麻酸(模型R²为0.9939，RMSEC为0.2631)：

Z₄＝7506.006+1591.699X₁₀₀₀-4004.831X₁₀₃₉-2633.555X₁₀₇₈+3094.150X₁₀₈₄

-2992.399X₁₀₉₄-3152.559X₁₁₄₆+1920.268X₁₁₆₈-2117.618X₁₁₇₅+3385.247X₁₁₈₈

-4188.517X₁₁₉₂-1635.062X₁₁₉₈-3958.504X₁₂₈₂+2138.468X₁₂₈₉-1634.357X₁₃₄₁-927.100X₁₃₈₁

-2385.798X₁₄₁₄-1935.714X₁₄₂₆+2986.193X₁₄₃₃+366.652X₁₄₃₄-1159.740X₁₄₅₁

-3490.094X₁₅₁₇-1700.527X₁₅₇₅-9384.497X₁₅₉₃+14487X₁₅₉₄-10710X₁₅₉₅-5118.240X₁₅₉₈

+12228X₁₆₀₅-5748.613X₁₆₀₆+4074.247X₁₆₁₄-2355.427X₁₆₃₇-2148.543X₁₆₇₃-959.665X₁₆₉₇

-921.807X₁₇₁₁+1177.955X₁₇₁₂-417.392X₁₇₁₄-882.869X₁₇₂₆+1267.076X₁₇₂₇-877.131X₁₇₂₉

-423.845X₁₇₃₇+822.571X₁₇₄₃-910.270X₁₇₅₈+952.185X₁₇₇₁-1287.709X₁₇₇₃；

硬脂酸(模型R²为0.9943，RMSEC为0.0999)：

Z₅＝-11.140+3255.650X₁₀₀₈+3240.414X₁₀₁₄-2422.871X₁₀₂₄+2690.775X₁₀₆₃

-2097.645X₁₀₇₈-2325.170X₁₀₈₃-6746.814X₁₀₈₈-827.644X₁₀₉₆-2125.511X₁₁₁₇-5597.104X₁₁₅₀

+4606.841X₁₁₆₄+4574.228X₁₁₉₃+8693.499X₁₂₃₈-8975.580X₁₂₆₀+2426.146X₁₂₆₄-3632.476X

₁₂₈₆-53.705X₁₂₈₇+1334.176X₁₃₀₂-5409.880X₁₃₃₉+7075.527X₁₃₅₆-9719.556X₁₃₅₇+11419X₁₃₅₉

+3817.161X₁₃₈₉+2471.477X₁₃₉₉+5158.646X₁₄₀₂+5644.206X₁₄₂₁+1723.091X₁₄₄₂

+6683.316X₁₄₆₅+2456.902X₁₄₈₁+5891.291X₁₅₀₁-757.199X₁₅₁₈-2085.148X₁₅₃₁

-5314.264X₁₅₃₂-6929.488X₁₅₃₃+13.519X₁₅₄₂+6958.481X₁₅₄₃-2976.579X₁₅₈₅+4905.371X₁₆₃₁

+1954.005X₁₆₅₁+1529.944X₁₆₅₆-1580.192X₁₆₆₉；

花生一烯酸(模型R²为0.9990，RMSEC为0.0899)：

Z₆＝2.815+1284.950X₁₀₂₆-1461.413X₁₀₃₆-404.856X₁₀₄₁-260.283X₁₀₅₃+805.799X₁₀₅₆

-1159.894X₁₀₇₀+1577.242X₁₀₇₈-808.352X₁₀₈₅-1020.663X₁₀₉₁+389.510X₁₁₄₁

-1156.822X₁₁₈₁-1502.081X₁₁₉₇+985.943X₁₂₀₅+342.352X₁₂₁₇+1878.121X₁₂₄₁-1334.758X₁₂₆₃

+475.320X₁₂₆₈-1755.228X₁₂₇₉+1593.834X₁₂₈₉-672.329X₁₃₂₁+631.634X₁₃₄₈-644.578X₁₃₆₁

+1496.912X₁₃₈₅+694.086X₁₃₉₀-812.705X₁₃₉₄-446.055X₁₄₁₃+501.298X₁₄₃₃+458.211X₁₄₄₇

-530.111X₁₄₆₀+866.661X₁₅₂₃+597.164X₁₅₄₅-1612.002X₁₅₇₇+1330.369X₁₆₂₃+1429.233X₁₆₂₅

-2031.183X₁₆₃₄-1173.266X₁₆₄₀+445.451X₁₆₇₁+895.088X₁₆₇₅-298.929X₁₇₃₈；

芥酸(模型R²为0.9997，RMSEC为0.0972)：

Z₇＝10.284+4659.521X₁₀₇₃-2923.995X₁₁₀₅+1345.165X₁₁₁₀-6375.893X₁₁₂₃+4111.558X₁₁₄₀

+5199.399X₁₁₄₅+4806.482X₁₁₅₅-6702.304X₁₁₉₀+7739.103X₁₂₀₁+4164.891X₁₂₀₅

-1952.258X₁₂₂₉-1801.109X₁₂₅₈-10555X₁₂₇₃-4023.318X₁₂₉₁-6897.226X₁₃₀₂-3135.779X₁₃₅₁

-2720.412X₁₃₇₄+7058.360X₁₄₀₃-4950.972X₁₄₀₈-3924.462X₁₄₀₉-7708.661X₁₄₃₃

-2204.383X₁₄₄₄-3509.622X₁₄₅₈-5880.686X₁₅₁₁-4326.789X₁₅₂₇+6809.207X₁₅₂₈+12682X₁₅₃₄

+4803.327X₁₅₆₈+7116.305X₁₆₀₁+3908.016X₁₆₃₀+5744.782X₁₆₃₁-3704.439X₁₆₃₇

+6322.864X₁₆₅₀-5132.705X₁₆₆₀-6511.577X₁₆₆₁+4127.098X₁₆₇₉+2712.062X₁₆₈₇

-5831.724X₁₆₈₈-1914.506X₁₆₉₀+1515.624X₁₆₉₃+1262.943X₁₆₉₅-2057.092X₁₆₉₈

+3295.167X₁₆₉₉-1365.754X₁₇₀₁-919.912X₁₇₃₉-1699.429X₁₇₄₃+2978.033X₁₇₉₁

+2350.872X₁₇₉₃；

酸价(模型R²为0.9684，RMSEC为0.2216)：

Z₈＝0.119-2538.959X₁₀₀₄-3100.461X₁₀₃₃+4886.971X₁₀₃₅+2897.828X₁₀₅₁-2006.278X₁₀₆₃

+2760.098X₁₀₇₉-2640.380X₁₁₁₈+2648.164X₁₁₄₃-2753.983X₁₁₇₀+5003.504X₁₁₈₀

-1000.648X₁₂₀₄+1400.125X₁₂₁₂-3286.331X₁₂₃₀+4502.311X₁₂₅₀-2293.994X₁₂₆₄

-5782.938X₁₂₆₈+2145.841X₁₂₈₉+10122X₁₃₀₈-11554X₁₃₀₉+2774.508X₁₃₁₆+567.939X₁₃₂₉

-1925.473X₁₃₇₈+1430.852X₁₃₉₆-3064.104X₁₄₁₀+2059.191X₁₄₁₄+21246X₁₄₇₇-20903X₁₄₇₈

-1061.150X₁₄₉₇+3780.992X₁₅₇₇-4940.442X₁₆₁₇；

过氧化值(模型R²为0.9972，RMSEC为0.4195)：

Z₉＝-1.095+7893.961X₁₀₀₆+8674.625X₁₀₁₈-2189.742X₁₀₃₈+2747.823X₁₀₅₃-4545.457X₁₀₇₀

-5891.021X₁₀₈₈-11945X₁₁₁₀+7041.339X₁₁₁₃-6603.093X₁₁₂₇+8646.126X₁₁₃₂+5028.386X₁₁₃₉

+7597.665X₁₂₂₂-11089X₁₂₃₂-1108.292X₁₂₄₀+9275.776X₁₂₇₇+8261.865X₁₂₈₅-12310X₁₃₆₉

+13064X₁₃₈₂-5946.595X₁₄₀₅+7664.969X₁₄₁₄-11912X₁₄₃₅+9507.563X₁₄₆₀+7903.013X₁₄₉₄

-13248X₁₅₅₉-7108.786X₁₅₇₁-2923.063X₁₆₂₇-3778.807X₁₆₆₀-4053.265X₁₆₈₂-3191.812X₁₇₆₅

+4108.125X₁₇₇₄-5542.854X₁₇₈₄+7032.594X₁₇₈₇；

碘价(模型R²为0.9935，RMSEC为0.5154)：

Z₁₀＝70.379+5655.836X₁₀₀₇+3835.644X₁₀₀₈+2255.332X₁₀₁₂+8315.045X₁₀₁₇+5967.677X₁₀₂₀

+5946.963X₁₀₄₁+4961.757X₁₀₄₂+8824.426X₁₀₅₁+9546.300X₁₀₆₆+2556.441X₁₀₇₈

+10214X₁₁₀₀-3767.757X₁₁₈₅-6114.234X₁₁₉₀-7551.222X₁₂₂₈-12332X₁₂₄₃-2271.972X₁₂₆₆

-8074.840X₁₂₈₃-1175.525X₁₃₀₅+3644.068X₁₃₁₁-12947X₁₃₁₂-6439.407X₁₃₄₆-7045.824X₁₃₅₄

-14294X₁₃₆₈+8944.760X₁₃₆₉-8323.329X₁₃₇₄-6197.887X₁₃₇₈-8325.039X₁₃₈₂-9125.409X₁₃₉₄

-3157.877X₁₃₉₈-2646.257X₁₃₉₉+5074.726X₁₄₂₃-4409.695X₁₄₄₄+16935X₁₄₄₅-7959.432X₁₄₉₁

-9135.176X₁₅₁₀+3952.970X₁₆₅₃+9651.345X₁₆₅₄+2212.403X₁₆₅₅+4346.486X₁₆₇₁

-4980.917X₁₆₈₈+5301.666X₁₆₉₈+4013.108X₁₇₀₄+2637.914X₁₇₀₅-3167.676X₁₇₀₇

+5617.745X₁₇₀₈-3577.202X₁₇₁₇-4002.850X₁₇₂₁+5140.303X₁₇₂₉+4737.340X₁₇₃₃

-2469.897X₁₇₄₃-3276.167X₁₇₄₆-959.039X₁₇₄₉+2539.913X₁₇₆₃-2006.402X₁₇₇₁-3899.264X₁₇₈₂

-5321.969X₁₇₈₅+3019.113X₁₇₉₄；

色泽黄值(模型R²为0.9873，RMSEC为0.2335)：

Z₁₁＝7.200+15691X₁₀₀₃+14907X₁₀₀₅-12573X₁₀₅₁+50407X₁₀₉₈+38618X₁₁₄₉-22311X₁₁₅₇

+21980X₁₁₇₄-46337X₁₁₇₉+15773X₁₁₈₅+21827X₁₁₉₉+20574X₁₂₁₈+49588X₁₂₉₅-52794X₁₂₉₉

+46525X₁₃₅₇-23933X₁₃₆₄+11818X₁₃₇₂-10222X₁₃₉₈+25684X₁₄₁₉-25402X₁₄₃₀+32483X₁₄₇₉

-42643X₁₅₀₃-68673X₁₅₁₇+74279X₁₅₂₆-30897X₁₅₂₇+120133X₁₅₃₁-149057X₁₅₃₂+167820X₁₅₃₃

-88396X₁₅₃₇-28409X₁₅₅₆+68851X₁₆₀₈-30206X₁₆₁₀-46112X₁₆₄₃+11539X₁₆₆₂+10951X₁₆₈₀

-10465X₁₆₉₇+7317.372X₁₇₃₇-14058X₁₇₃₉-7434.215X₁₇₆₇+13105X₁₇₇₂+12795X₁₇₈₈

-27958X₁₇₉₀-6001.537X₁₇₉₉；

色泽红值(模型R²为0.9838，RMSEC为0.5154)：

Z₁₂＝0.891-1868.697X₁₀₃₁-2572.488X₁₀₅₁-288.888X₁₀₇₁+5057.192X₁₁₀₀-3072.222X₁₁₀₅

+1257.980X₁₁₂₃-3508.014X₁₁₂₅+3855.554X₁₁₅₀-3758.369X₁₁₆₁+1483.714X₁₁₆₉

-2301.038X₁₁₇₉+2555.615X₁₁₈₄+2017.671X₁₂₁₈-2454.693X₁₂₃₂+4498.395X₁₂₉₄

+2850.912X₁₂₉₅-4340.720X₁₃₀₀-2012.020X₁₃₁₂+3303.660X₁₃₅₇-1920.040X₁₃₇₀

+707.359X₁₃₈₉+3523.153X₁₄₇₀-3226.896X₁₄₉₂+1275.462X₁₅₂₈+5537.216X₁₅₂₉

-4646.654X₁₅₃₈+1219.908X₁₅₅₁+4206.180X₁₅₆₅-6255.102X₁₅₉₀+5286.438X₁₆₀₂

+5131.157X₁₆₃₉-5381.917X₁₆₄₄+3517.545X₁₆₅₇+2421.794X₁₆₉₀-2134.497X₁₆₉₄

-1276.800X₁₇₀₂+1863.362X₁₇₀₃-2434.757X₁₇₄₃+2167.080X₁₇₄₄+1611.073X₁₇₇₈

-851.937X₁₇₉₀；

黏度(模型R²为0.8395，RMSEC为2.9731)：

Z₁₃＝53.980-21386X₁₀₁₅-26783X₁₀₂₃+37751X₁₀₂₅-26813X₁₁₁₇-35111X₁₁₂₉+25465X₁₁₄₀

-32237X₁₁₆₃-21558X₁₁₇₉-21465X₁₂₁₀-12785X₁₂₄₉-61270X₁₂₆₁+57634X₁₃₅₉-60468X₁₃₆₅

-65030X₁₃₇₀-9524.085X₁₃₇₉+26949X₁₃₈₂-21009X₁₄₁₄-11847X₁₄₂₆-34364X₁₄₆₈+28553X₁₄₉₃

-59846X₁₅₀₈-66740X₁₅₆₀+34126X₁₅₇₂+51567X₁₅₈₂-69656X₁₅₉₀-33573X₁₅₉₁+74660X₁₆₀₄

-41775X₁₆₂₈-15134X₁₆₉₆-10211X₁₇₉₂；

8)预测模型的验证：按照实施例1步骤8)计算出上述9个品质指标验证集样品的近红外预测值Z_m，与该样品理化指标测试值进行对比，计算出各指标验证集的决定系数及校验标准差，结果见表5。由表5可知，上述模型近红外检测的精度较高，模型的决定系数R²均大于0.8997，理化检测值和近红外测量值的偏差小，具有较好的实用价值。

表5菜籽油12项品质指标近红外预测模型的检验效果

9)近红外指纹波长的确定：按照实施例1步骤9)方法确定各理化指标的指纹波长，结果如下：

亚油酸：1125nm、1194nm、1236nm、1657nm、1662nm、1668nm、1687nm、1740nm；

油酸：1107nm、1218nm、1409nm、1436nm、1698nm、1743nm；

亚麻酸：1039nm、1673nm、1697nm、1727nm、1729nm、1758nm；

硬脂酸：1150nm、1164nm、1193nm、1238nm、1260nm、1359nm、1389nm、1402nm、1421nm、1465nm；

花生一烯酸：1197nm、1241nm、1634nm、1738nm；

芥酸：1190nm、1201nm、1205nm、1403nm、1433nm、1534nm、1650nm、1679nm、1743nm、1791nm；

酸价：1035nm、1170nm、1180nm、1414nm；

过氧化值：1018nm、1222nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1682nm、1765nm、1787nm；

碘价：1100nm、1688nm、1698nm、1708nm、1717nm、1721nm、1729nm、1733nm、1746nm、1763nm；

色泽黄值：1098nm、1179nm、1479nm、1643nm；

色泽红值：1100nm、1357nm、1657nm、1690nm、1694nm；

黏度：1163nm、1210nm、1604nm、1628nm、1696nm；

上述波长允许有±2nm的偏差。

实施例4：异常样本剔除对模型精度的影响

以菜籽油亚麻酸指标为例，对异常样本剔除作进一步说明(其他指标异常样本剔除同本实施例，本发明不一一穷举，本领域的技术人员利用本发明的方法可以重复本发明的效果)。

菜籽油的近红外光谱扫描方法、样品亚麻酸检测方法同实施例1，采用偏最小二乘法分别对未经处理的原始样本集以及经过异常样本提除后的样本集，建立亚麻酸定标模型，两种模型的定标集决定系数及标准差，和验证集标准差的计算结果见表6。由表6可知，经过异常样本剔除后所建立的模型的决定系数R²为0.9794，明显高于未经处理的0.8967，定标集、验证集的标准差也均明显小于未经处理的，说明异常样本剔除对建模影响非常大。

表6剔除异常样本处理对模型精度的影响

实施例5：预处理对模型精度的影响

以菜籽油棕榈酸指标为例，对预处理方法作进一步说明(其他指标的预处理方法同本实施例，本发明不一一穷举，本领域的技术人员利用本发明的方法可以重复本发明的效果)。

菜籽油脂肪酸组成的则定及近红外光谱扫描的方法同实施例1，运用S-G卷积平滑法进行消噪，分别运用一阶导数、二阶导数、多元散射校正、净分析信号、正交信号校正、标准正态变量变换、标准化、均值中心化、一阶导数加正交信号校正、一阶导数加正交信号校正加标准化、一阶导数加标准正态变量变换、一阶导数加标准正态变量变换加均值中心化13种方法对消噪后的菜籽油近红外光谱进行预处理，得到预处理后的吸光值，每种预处理后的数据集均采用偏最小二乘法构建棕榈酸定标模型，模型的定标集、验证集的决定系数及标准差的计算结果见表7。

表7不同的预处理方法对模型精度的影响

由表7可知，S-G平滑后模型精度明显提高，另外经过正交信号校正法处理的预测模型的定标集决定系数R²(0.9999)最高，标准差(0.0430)最小，是最适宜的菜籽油棕榈酸近红外光谱的预处理方法。

上述的实施例仅为本发明的优选技术方案，而不应视为对于本发明的限制，本申请中的实施例及实施例中的特征在不冲突的情况下，可以相互任意组合。本发明的保护范围应以权利要求记载的技术方案，包括权利要求记载的技术方案中技术特征的等同替换方案为保护范围。即在此范围内的等同替换改进，也在本发明的保护范围之内。

Claims (10)

1.一种菜籽油品质的近红外检测方法，其特征在于，所述方法包括以下步骤：

1)收集菜籽油样品；

2)对所述步骤1)收集的菜籽油样品进行理化指标测试，得到理化指标测试值，记为y_mj，其中：m是第m个指标，m＝1,2,3，…，13；当m＝1时记为棕榈酸含量，m＝2时记为亚油酸含量，m＝3时记为油酸含量，m＝4时记为亚麻酸含量，m＝5时记为硬脂酸含量，m＝6时记为花生一烯酸含量，m＝7时记为芥酸含量，m＝8时记为酸价，m＝9时记为过氧化值，m＝10时记为碘价，m＝11时记为色泽黄值，m＝12时记为色泽红值，m＝13时记为黏度；j为第j个样品，共n个样品，n≥90；

3)对所述步骤1)收集的菜籽油样品进行近红外扫描，得到吸光值为Xij，其中i表示波长数值；

4)异常样本剔除和样品集划分：对所述步骤2)得到的菜籽油理化指标测试值y_mj和所述步骤3)得到的吸光值Xij采用学生残差法将检测出来的异常样本剔除掉后，采用Kennard-Stone法将样品集划分为校正集和验证集；

5)吸光值预处理：对所述步骤4)校正集吸光值进行预处理，得到校正集预处理吸光值A_ij；

6)筛选近红外光谱特征波长：采用逐步回归法对步骤2)得到的理化指标y_mj和步骤5)得到的校正集预处理吸光值A_ij进行分析，筛选得到菜籽油品质的近红外光谱特征波长；

7)建立菜籽油品质预测模型：采用多元线性回归法对步骤2)得到的理化指标y_mj和步骤5)得到的近红外光谱特征波长的预处理吸光值A_ij进行分析，建立菜籽油品质的预测模型为Z_m＝b_m+∑a_miXi，其中Z_m为菜籽油品质的近红外检测值，Xi为A_ij中第j个样品的吸光值，b_m、a_mi为回归系数；

8)预测模型的验证：将所述步骤4)划分出的样品验证集按照所述步骤5)方法进行预处理，得到验证集预处理吸光值，代入所述步骤7)的预测模型，得到菜籽油品质指标预测值，与所述步骤2)该样品理化指标测试值进行比对，得到模型预测的精确度；

9)近红外指纹波长的确定：

完成菜籽油品质的检测。

2.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤3)的近红外扫描温度为20-40℃，光程为5mm。

3.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤3)的菜籽油样品进行近红外扫描波长为：波长i＝1000nm，1001nm，1002nm，…，1799nm。

4.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤5)预处理方法为一阶导数或二阶导数或多元散射校正或或净分析信号或正交信号校正或标准正态变量变换或去趋势校正或标准化或均值中心化中的一种或其组合。

5.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤6)菜籽油品质的近红外光谱特征波长为：

棕榈酸：1097nm、1102nm、1260nm、1266nm、1286nm、1299nm、1355nm、1372nm、1428nm、1438nm、1450nm、1521nm、1534nm、1536nm、1619nm、1628nm、1629nm、1633nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

亚油酸：1040nm、1072nm、1077nm、1098nm、1108nm、1110nm、1125nm、1176nm、1181nm、1187nm、1194nm、1236nm、1334nm、1344nm、1395nm、1401nm、1442nm、1451nm、1457nm、1497nm、1608nm、1641nm、1643nm、1646nm、1657nm、1662nm、1668nm、1682nm、1687nm、1740nm、1780nm、1791nm、1794nm；

油酸：1024nm、1027nm、1028nm、1039nm、1065nm、1076nm、1092nm、1107nm、1115nm、1131nm、1183nm、1200nm、1203nm、1215nm、1218nm、1226nm、1227nm、1304nm、1325nm、1343nm、1376nm、1401nm、1409nm、1436nm、1455nm、1505nm、1554nm、1570nm、1622nm、1629nm、1635nm、1662nm、1663nm、1665nm、1672nm、1685nm、1698nm、1743nm、1788nm；

亚麻酸：1000nm、1039nm、1078nm、1084nm、1094nm、1146nm、1168nm、1175nm、1188nm、1192nm、1198nm、1282nm、1289nm、1341nm、1381nm、1414nm、1426nm、1433nm、1434nm、1451nm、1517nm、1575nm、1593nm、1594nm、1595nm、1598nm、1605nm、1606nm、1614nm、1637nm、1673nm、1697nm、1711nm、1712nm、1714nm、1726nm、1727nm、1729nm、1737nm、1743nm、1758nm、1771nm、1773nm；

硬脂酸：1008nm、1014nm、1024nm、1063nm、1078nm、1083nm、1088nm、1096nm、1117nm、1150nm、1164nm、1193nm、1238nm、1260nm、1264nm、1286nm、1287nm、1302nm、1339nm、1356nm、1357nm、1359nm、1389nm、1399nm、1402nm、1421nm、1442nm、1465nm、1481nm、1501nm、1518nm、1531nm、1532nm、1533nm、1542nm、1543nm、1585nm、1631nm、1651nm、1656nm、1669nm；

花生一烯酸：1026nm、1036nm、1041nm、1053nm、1056nm、1070nm、1078nm、1085nm、1091nm、1141nm、1181nm、1197nm、1205nm、1217nm、1241nm、1263nm、1268nm、1279nm、1289nm、1321nm、1348nm、1361nm、1385nm、1390nm、1394nm、1413nm、1433nm、1447nm、1460nm、1523nm、1545nm、1577nm、1623nm、1625nm、1634nm、1640nm、1671nm、1675nm、1738nm；

芥酸：1073nm、1105nm、1110nm、1123nm、1140nm、1145nmn、1155nm、1190nm、1201nm、1205nm、1229nm、1258nm、1273nm、1291nm、1302nm、1351nm、1374nm、1403nm、1408nm、1409nm、1433nm、1444nm、1458nm、1511nm、1527nm、1528nm、1534nm、1568nm、1601nm、1630nm、1631nm、1637nm、1650nm、1660nm、1661nm、1679nm、1687nm、1688nm、1690nm、1693nm、1695nm、1698nm、1699nm、1701nm、1739nm、1743nm、1791nm、1793nm；

酸价：1004nm、1033nm、1035nm、1051nm、1063nm、1079nm、1118nm、1143nm、1170nm、1180nm、1204nm、1212nm、1230nm、1250nm、1264nm、1268nm、1289nm、1308nm、1309nm、1316nm、1329nm、1378nm、1396nm、1410nm、1414nm、1477nm、1478nm、1497nm、1577nm、1617nm；

过氧化值：1006nm、1018nm、1038nm、1053nm、1070nm、1088nm、1110nm、1113nm、1127nm、1132nm、1139nm、1222nm、1232nm、1240nm、1277nm、1285nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1494nm、1559nm、1571nm、1627nm、1660nm、1682nm、1765nm、1774nm、1784nm、1787nm；

碘价：1007nm、1008nm、1012nm、1017nm、1020nm、1041nm、1042nm、1051nm、1066nm、1078nm、1100nm、1185nm、1190nm、1228nm、1243nm、1266nm、1283nm、1305nm、1311nm、1312nm、1346nm、1354nm、1368nm、1369nm、1374nm、1378nm、1382nm、1394nm、1398nm、1399nm、1423nm、1444nm、1445nm、1491nm、1510nm、1653nm、1654nm、1655nm、1671nm、1688nm、1698nm、1704nm、1705nm、1707nm、1708nm、1717nm、1721nm、1729nm、1733nm、1743nm、1746nm、1749nm、1763nm、1771nm、1782nm、1785nm、1794nm；

色泽黄值：1003nm、1005nm、1051nm、1098nm、1149nm、1157nm、1174nm、1179nm、1185nm、1199nm、1218nm、1295nm、1299nm、1357nm、1364nm、1372nm、1398nm、1419nm、1430nm、1479nm、1503nm、1517nm、1526nm、1527nm、1531nm、1532nm、1533nm、1537nm、1556nm、1608nm、1610nm、1643nm、1662nm、1680nm、1697nm、1737nm、1739nm、1767nm、1772nm、1788nm、1790nm、1799nm；

色泽红值：1031nm、1051nm、1071nm、1100nm、1105nm、1123nm、1125nm、1150nm、1161nm、1169nm、1179nm、1184nm、1218nm、1232nm、1294nm、1295nm、1300nm、1312nm、1357nm、1370nm、1389nm、1470nm、1492nm、1528nm、1529nm、1538nm、1551nm、1565nm、1590nm、1602nm、1639nm、1644nm、1657nm、1690nm、1694nm、1702nm、1703nm、1743nm、1744nm、1778nm、1790nm；

黏度：1015nm、1023nm、1025nm、1117nm、1129nm、1140nm、1163nm、1179nm、1210nm、1249nm、1261nm、1359nm、1365nm、1370nm、1379nm、1382nm、1414nm、1426nm、1468nm、1493nm、1508nm、1560nm、1572nm、1582nm、1590nm、1591nm、1604nm、1628nm、1696nm、1792nm；

上述特征波长允许有±2nm的偏差。

6.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤7)中预测模型的评价指标为决定系数R²、校正集标准差RMSEC和验证集标准差RMSEP。

7.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤7)中建立菜籽油品质预测模型：

棕榈酸：

Z₁＝7.122-554.923X₁₀₉₇+559.861X₁₁₀₂+1154.024X₁₂₆₀-1717.392X₁₂₆₆+1001.188X₁₂₈₆

-520.795X₁₂₉₉-400.896X₁₃₅₅+485.139X₁₃₇₂-540.567X₁₄₂₈+1052.793X₁₄₃₈-374.591X₁₄₅₀

+418.760X₁₅₂₁+501.261X₁₅₃₄-1063.741X₁₅₃₆+386.163X₁₆₁₉-2644.180X₁₆₂₈+2391.044X₁₆₂₉

-666.401X₁₆₃₃+729.704X₁₆₃₉-765.265X₁₆₇₀+1022.987X₁₆₇₉-969.315X₁₆₈₇+482.195X₁₆₉₃-

316.015X₁₇₇₇+277.649X₁₇₉₁；

亚油酸：

Z₂＝12.881+1494.756X₁₀₄₀+1188.261X₁₀₇₂-2106.146X₁₀₇₇+2004.052X₁₀₉₈-3264.258X₁₁₀₈-

2164.257X₁₁₁₀+2926.439X₁₁₂₅-1919.937X₁₁₇₆+3415.763X₁₁₈₁-1835.691X₁₁₈₇

+2643.582X₁₁₉₄-1986.239X₁₂₃₆+2298.319X₁₃₃₄-1712.230X₁₃₄₄-2152.860X₁₃₉₅

+1301.989X₁₄₀₁+669.145X₁₄₄₂-2266.494X₁₄₅₁+1861.268X₁₄₅₇-1366.906X₁₄₉₇

-2751.327X₁₆₀₈+1042.784X₁₆₄₁+7273.727X₁₆₄₃-4483.278X₁₆₄₆-3243.621X₁₆₅₇

+6191.148X₁₆₆₂-3773.298X₁₆₆₈+2006.371X₁₆₈₂-1839.867X₁₆₈₇+883.481X₁₇₄₀

-1004.192X₁₇₈₀+1035.786X₁₇₉₁-1120.310X₁₇₉₄；

油酸：

Z₃＝131.518+54.648X₁₀₂₄-145.218X₁₀₂₇+238.569X₁₀₂₈-124.154X₁₀₃₉-38.688X₁₀₆₅

-187.688X₁₀₇₆-139.253X₁₀₉₂+222.844X₁₁₀₇+215.092X₁₁₁₅+63.101X₁₁₃₁+86.268X₁₁₈₃

-134.256X₁₂₀₀-78.406X₁₂₀₃-74.438X₁₂₁₅+224.804X₁₂₁₈+267.659X₁₂₂₆-134.32748X₁₂₂₇

-212.718X₁₃₀₄-88.665X₁₃₂₅-97.920X₁₃₄₃+100.588X₁₃₇₆-120.263X₁₄₀₁+137.713X₁₄₀₉

-219.328X₁₄₃₆+145.883X₁₄₅₅-189.689X₁₅₀₅+187.495X₁₅₅₄-299.930X₁₅₇₀+201.052X₁₆₂₂

+168.477X₁₆₂₉-186.007X₁₆₃₅+446.389X₁₆₆₂-323.377X₁₆₆₃+367.357X₁₆₆₅-135.115X₁₆₇₂

+72.825X₁₆₈₅-50.836X₁₆₉₈+46.310X₁₇₄₃-53.128X₁₇₈₈；

亚麻酸：

Z₄＝7506.006+1591.699X₁₀₀₀-4004.831X₁₀₃₉-2633.555X₁₀₇₈+3094.150X₁₀₈₄

-2992.399X₁₀₉₄-3152.559X₁₁₄₆+1920.268X₁₁₆₈-2117.618X₁₁₇₅+3385.247X₁₁₈₈

-4188.517X₁₁₉₂-1635.062X₁₁₉₈-3958.504X₁₂₈₂+2138.468X₁₂₈₉-1634.357X₁₃₄₁-927.100X₁₃₈₁

-2385.798X₁₄₁₄-1935.714X₁₄₂₆+2986.193X₁₄₃₃+366.652X₁₄₃₄-1159.740X₁₄₅₁

-3490.094X₁₅₁₇-1700.527X₁₅₇₅-9384.497X₁₅₉₃+14487X₁₅₉₄-10710X₁₅₉₅-5118.240X₁₅₉₈

+12228X₁₆₀₅-5748.613X₁₆₀₆+4074.247X₁₆₁₄-2355.427X₁₆₃₇-2148.543X₁₆₇₃-959.665X₁₆₉₇

-921.807X₁₇₁₁+1177.955X₁₇₁₂-417.392X₁₇₁₄-882.869X₁₇₂₆+1267.076X₁₇₂₇-877.131X₁₇₂₉

-423.845X₁₇₃₇+822.571X₁₇₄₃-910.270X₁₇₅₈+952.185X₁₇₇₁-1287.709X₁₇₇₃；

硬脂酸：

Z₅＝-11.140+3255.650X₁₀₀₈+3240.414X₁₀₁₄-2422.871X₁₀₂₄+2690.775X₁₀₆₃

-2097.645X₁₀₇₈-2325.170X₁₀₈₃-6746.814X₁₀₈₈-827.644X₁₀₉₆-2125.511X₁₁₁₇-5597.104X₁₁₅₀

+4606.841X₁₁₆₄+4574.228X₁₁₉₃+8693.499X₁₂₃₈-8975.580X₁₂₆₀+2426.146X₁₂₆₄-3632.476X

₁₂₈₆-53.705X₁₂₈₇+1334.176X₁₃₀₂-5409.880X₁₃₃₉+7075.527X₁₃₅₆-9719.556X₁₃₅₇+11419X₁₃₅₉

+3817.161X₁₃₈₉+2471.477X₁₃₉₉+5158.646X₁₄₀₂+5644.206X₁₄₂₁+1723.091X₁₄₄₂

+6683.316X₁₄₆₅+2456.902X₁₄₈₁+5891.291X₁₅₀₁-757.199X₁₅₁₈-2085.148X₁₅₃₁

-5314.264X₁₅₃₂-6929.488X₁₅₃₃+13.519X₁₅₄₂+6958.481X₁₅₄₃-2976.579X₁₅₈₅+4905.371X₁₆₃₁

+1954.005X₁₆₅₁+1529.944X₁₆₅₆-1580.192X₁₆₆₉；

花生一烯酸：

Z₆＝2.815+1284.950X₁₀₂₆-1461.413X₁₀₃₆-404.856X₁₀₄₁-260.283X₁₀₅₃+805.799X₁₀₅₆

-1159.894X₁₀₇₀+1577.242X₁₀₇₈-808.352X₁₀₈₅-1020.663X₁₀₉₁+389.510X₁₁₄₁

-1156.822X₁₁₈₁-1502.081X₁₁₉₇+985.943X₁₂₀₅+342.352X₁₂₁₇+1878.121X₁₂₄₁-1334.758X₁₂₆₃

+475.320X₁₂₆₈-1755.228X₁₂₇₉+1593.834X₁₂₈₉-672.329X₁₃₂₁+631.634X₁₃₄₈-644.578X₁₃₆₁

+1496.912X₁₃₈₅+694.086X₁₃₉₀-812.705X₁₃₉₄-446.055X₁₄₁₃+501.298X₁₄₃₃+458.211X₁₄₄₇

-530.111X₁₄₆₀+866.661X₁₅₂₃+597.164X₁₅₄₅-1612.002X₁₅₇₇+1330.369X₁₆₂₃+1429.233X₁₆₂₅

-2031.183X₁₆₃₄-1173.266X₁₆₄₀+445.451X₁₆₇₁+895.088X₁₆₇₅-298.929X₁₇₃₈；

芥酸：

Z₇＝10.284+4659.521X₁₀₇₃-2923.995X₁₁₀₅+1345.165X₁₁₁₀-6375.893X₁₁₂₃+4111.558X₁₁₄₀

+5199.399X₁₁₄₅+4806.482X₁₁₅₅-6702.304X₁₁₉₀+7739.103X₁₂₀₁+4164.891X₁₂₀₅

-1952.258X₁₂₂₉-1801.109X₁₂₅₈-10555X₁₂₇₃-4023.318X₁₂₉₁-6897.226X₁₃₀₂-3135.779X₁₃₅₁

-2720.412X₁₃₇₄+7058.360X₁₄₀₃-4950.972X₁₄₀₈-3924.462X₁₄₀₉-7708.661X₁₄₃₃

-2204.383X₁₄₄₄-3509.622X₁₄₅₈-5880.686X₁₅₁₁-4326.789X₁₅₂₇+6809.207X₁₅₂₈+12682X₁₅₃₄

+4803.327X₁₅₆₈+7116.305X₁₆₀₁+3908.016X₁₆₃₀+5744.782X₁₆₃₁-3704.439X₁₆₃₇

+6322.864X₁₆₅₀-5132.705X₁₆₆₀-6511.577X₁₆₆₁+4127.098X₁₆₇₉+2712.062X₁₆₈₇

-5831.724X₁₆₈₈-1914.506X₁₆₉₀+1515.624X₁₆₉₃+1262.943X₁₆₉₅-2057.092X₁₆₉₈

+3295.167X₁₆₉₉-1365.754X₁₇₀₁-919.912X₁₇₃₉-1699.429X₁₇₄₃+2978.033X₁₇₉₁

+2350.872X₁₇₉₃；

酸价：

Z₈＝0.119-2538.959X₁₀₀₄-3100.461X₁₀₃₃+4886.971X₁₀₃₅+2897.828X₁₀₅₁-2006.278X₁₀₆₃

+2760.098X₁₀₇₉-2640.380X₁₁₁₈+2648.164X₁₁₄₃-2753.983X₁₁₇₀+5003.504X₁₁₈₀

-1000.648X₁₂₀₄+1400.125X₁₂₁₂-3286.331X₁₂₃₀+4502.311X₁₂₅₀-2293.994X₁₂₆₄

-5782.938X₁₂₆₈+2145.841X₁₂₈₉+10122X₁₃₀₈-11554X₁₃₀₉+2774.508X₁₃₁₆+567.939X₁₃₂₉

-1925.473X₁₃₇₈+1430.852X₁₃₉₆-3064.104X₁₄₁₀+2059.191X₁₄₁₄+21246X₁₄₇₇-20903X₁₄₇₈

-1061.150X₁₄₉₇+3780.992X₁₅₇₇-4940.442X₁₆₁₇；

过氧化值：

Z₉＝-1.095+7893.961X₁₀₀₆+8674.625X₁₀₁₈-2189.742X₁₀₃₈+2747.823X₁₀₅₃-4545.457X₁₀₇₀

-5891.021X₁₀₈₈-11945X₁₁₁₀+7041.339X₁₁₁₃-6603.093X₁₁₂₇+8646.126X₁₁₃₂+5028.386X₁₁₃₉

+7597.665X₁₂₂₂-11089X₁₂₃₂-1108.292X₁₂₄₀+9275.776X₁₂₇₇+8261.865X₁₂₈₅-12310X₁₃₆₉

+13064X₁₃₈₂-5946.595X₁₄₀₅+7664.969X₁₄₁₄-11912X₁₄₃₅+9507.563X₁₄₆₀+7903.013X₁₄₉₄

-13248X₁₅₅₉-7108.786X₁₅₇₁-2923.063X₁₆₂₇-3778.807X₁₆₆₀-4053.265X₁₆₈₂-3191.812X₁₇₆₅

+4108.125X₁₇₇₄-5542.854X₁₇₈₄+7032.594X₁₇₈₇；

碘价：

Z₁₀＝70.379+5655.836X₁₀₀₇+3835.644X₁₀₀₈+2255.332X₁₀₁₂+8315.045X₁₀₁₇+5967.677X₁₀₂₀

+5946.963X₁₀₄₁+4961.757X₁₀₄₂+8824.426X₁₀₅₁+9546.300X₁₀₆₆+2556.441X₁₀₇₈

+10214X₁₁₀₀-3767.757X₁₁₈₅-6114.234X₁₁₉₀-7551.222X₁₂₂₈-12332X₁₂₄₃-2271.972X₁₂₆₆

-8074.840X₁₂₈₃-1175.525X₁₃₀₅+3644.068X₁₃₁₁-12947X₁₃₁₂-6439.407X₁₃₄₆-7045.824X₁₃₅₄

-14294X₁₃₆₈+8944.760X₁₃₆₉-8323.329X₁₃₇₄-6197.887X₁₃₇₈-8325.039X₁₃₈₂-9125.409X₁₃₉₄

-3157.877X₁₃₉₈-2646.257X₁₃₉₉+5074.726X₁₄₂₃-4409.695X₁₄₄₄+16935X₁₄₄₅-7959.432X₁₄₉₁

-9135.176X₁₅₁₀+3952.970X₁₆₅₃+9651.345X₁₆₅₄+2212.403X₁₆₅₅+4346.486X₁₆₇₁

-4980.917X₁₆₈₈+5301.666X₁₆₉₈+4013.108X₁₇₀₄+2637.914X₁₇₀₅-3167.676X₁₇₀₇

+5617.745X₁₇₀₈-3577.202X₁₇₁₇-4002.850X₁₇₂₁+5140.303X₁₇₂₉+4737.340X₁₇₃₃

-2469.897X₁₇₄₃-3276.167X₁₇₄₆-959.039X₁₇₄₉+2539.913X₁₇₆₃-2006.402X₁₇₇₁-3899.264X₁₇₈₂

-5321.969X₁₇₈₅+3019.113X₁₇₉₄；

色泽黄值：

Z₁₁＝7.200+15691X₁₀₀₃+14907X₁₀₀₅-12573X₁₀₅₁+50407X₁₀₉₈+38618X₁₁₄₉-22311X₁₁₅₇

+21980X₁₁₇₄-46337X₁₁₇₉+15773X₁₁₈₅+21827X₁₁₉₉+20574X₁₂₁₈+49588X₁₂₉₅-52794X₁₂₉₉

+46525X₁₃₅₇-23933X₁₃₆₄+11818X₁₃₇₂-10222X₁₃₉₈+25684X₁₄₁₉-25402X₁₄₃₀+32483X₁₄₇₉

-42643X₁₅₀₃-68673X₁₅₁₇+74279X₁₅₂₆-30897X₁₅₂₇+120133X₁₅₃₁-149057X₁₅₃₂+167820X₁₅₃₃

-88396X₁₅₃₇-28409X₁₅₅₆+68851X₁₆₀₈-30206X₁₆₁₀-46112X₁₆₄₃+11539X₁₆₆₂+10951X₁₆₈₀

-10465X₁₆₉₇+7317.372X₁₇₃₇-14058X₁₇₃₉-7434.215X₁₇₆₇+13105X₁₇₇₂+12795X₁₇₈₈

-27958X₁₇₉₀-6001.537X₁₇₉₉；

色泽红值：

Z₁₂＝0.891-1868.697X₁₀₃₁-2572.488X₁₀₅₁-288.888X₁₀₇₁+5057.192X₁₁₀₀-3072.222X₁₁₀₅

+1257.980X₁₁₂₃-3508.014X₁₁₂₅+3855.554X₁₁₅₀-3758.369X₁₁₆₁+1483.714X₁₁₆₉

-2301.038X₁₁₇₉+2555.615X₁₁₈₄+2017.671X₁₂₁₈-2454.693X₁₂₃₂+4498.395X₁₂₉₄

+2850.912X₁₂₉₅-4340.720X₁₃₀₀-2012.020X₁₃₁₂+3303.660X₁₃₅₇-1920.040X₁₃₇₀

+707.359X₁₃₈₉+3523.153X₁₄₇₀-3226.896X₁₄₉₂+1275.462X₁₅₂₈+5537.216X₁₅₂₉

-4646.654X₁₅₃₈+1219.908X₁₅₅₁+4206.180X₁₅₆₅-6255.102X₁₅₉₀+5286.438X₁₆₀₂

+5131.157X₁₆₃₉-5381.917X₁₆₄₄+3517.545X₁₆₅₇+2421.794X₁₆₉₀-2134.497X₁₆₉₄

-1276.800X₁₇₀₂+1863.362X₁₇₀₃-2434.757X₁₇₄₃+2167.080X₁₇₄₄+1611.073X₁₇₇₈

-851.937X₁₇₉₀；

黏度：

Z₁₃＝53.980-21386X₁₀₁₅-26783X₁₀₂₃+37751X₁₀₂₅-26813X₁₁₁₇-35111X₁₁₂₉+25465X₁₁₄₀

-32237X₁₁₆₃-21558X₁₁₇₉-21465X₁₂₁₀-12785X₁₂₄₉-61270X₁₂₆₁+57634X₁₃₅₉-60468X₁₃₆₅

-65030X₁₃₇₀-9524.085X₁₃₇₉+26949X₁₃₈₂-21009X₁₄₁₄-11847X₁₄₂₆-34364X₁₄₆₈+28553X₁₄₉₃

-59846X₁₅₀₈-66740X₁₅₆₀+34126X₁₅₇₂+51567X₁₅₈₂-69656X₁₅₉₀-33573X₁₅₉₁+74660X₁₆₀₄

-41775X₁₆₂₈-15134X₁₆₉₆-10211X₁₇₉₂。

8.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤9)指纹波长确定方法为对所述步骤7)预测模型的每个特征波长对应的回归系数进行显著性分析，极显著的特征波长即为该菜籽油理化指标的近红外光谱指纹波长。

9.根据权利要求1所述的菜籽油品质的近红外检测方法，其特征在于：所述步骤9)菜籽油品质的近红外光谱指纹波长为：

亚油酸：1125nm、1194nm、1236nm、1657nm、1662nm、1668nm、1687nm、1740nm；

油酸：1107nm、1218nm、1409nm、1436nm、1698nm、1743nm；

亚麻酸：1039nm、1673nm、1697nm、1727nm、1729nm、1758nm；

棕榈酸：1260nm、1266nm、1286nm、1372nm、1428nm、1438nm、1639nm、1670nm、1679nm、1687nm、1693nm、1777nm、1791nm；

硬脂酸：1150nm、1164nm、1193nm、1238nm、1260nm、1359nm、1389nm、1402nm、1421nm、1465nm；

花生一烯酸：1197nm、1241nm、1634nm、1738nm；

芥酸：1190nm、1201nm、1205nm、1403nm、1433nm、1534nm、1650nm、1679nm、1743nm、1791nm；

酸价：1035nm、1170nm、1180nm、1414nm；

过氧化值：1018nm、1222nm、1369nm、1382nm、1405nm、1414nm、1435nm、1460nm、1682nm、1765nm、1787nm；

碘价：1100nm、1688nm、1698nm、1708nm、1717nm、1721nm、1729nm、1733nm、1746nm、1763nm；

色泽黄值：1098nm、1179nm、1479nm、1643nm；

色泽红值：1100nm、1357nm、1657nm、1690nm、1694nm；

黏度：1163nm、1210nm、1604nm、1628nm、1696nm；

上述波长允许有±2nm的偏差。

10.权利要求1-9任意一项所述的菜籽油品质的近红外检测方法在玉米油或橄榄油或山茶油或棕榈油或芥花子油或葵花子油或大豆油或芝麻油检测中的应用。