[发明专利]一种基于遗传算法自适应选择分段点的红外光谱PLS建模方法

权利要求书

1.一种基于遗传算法自适应选择分段点的红外光谱PLS建模方法，其特征在于，包括以下步骤：

步骤一，样品的划分，每种类别准备多个不同的样品数据，每种类别的样品划分为三部分SET1，SET2，SET3，SET1用于分段模型的建立，SET2用于分段模型的测试，SET3用于最终分段模型的验证，由SE T1构成的矩阵为校正集，由SET2构成的矩阵为预测子集，由SET3构成的矩阵为独立预测集；

步骤二，设置Y变量的值，如果有n种类别，则这些类别的Y变量的值分别设置为1，2，3，…n，如有6种类别，则Y变量的值设置为1，2，3，4，5，6；

步骤三，对校正集样品做主成分分析，得到一个得分矩阵Sc、投影方向矩阵T和每列平均值构成的行矩阵M；

步骤四，创建初始种群：包括设置个体每位的进制数，代表最大的分段数p，如分p段则采用p进制，如p＝3，表示个体每位有0，1，2这3种取值可能；设置个体的长度n，代表类别数，即Crom矩阵的列数；设置种群的大小v，代表一个种群内个体的数量，即Crom矩阵的行数，如p＝3，n＝6，v＝4，随机产生的初始种群如下所示：

Crom=111111111222110022100000abcdef]]>

a，b，c，d，e，f表示类别，Crom矩阵的每一行对应一个个体，个体的每一位数字对应一种类别，分别对应a，b，c，d，e，f这6种类别，数字相同的类别归为一段，第一行[1 1 1 1 1 1]表示，a，b，c，d，e，f归为一段，即不分段的形式，第二行[1 1 1 2 2 2]表示a，b，c三类归为一段，d，e，f归为一段，一共只出现两种不同的数字，即全部类别分为两段，第三段[1 1 0 0 2 2]，表示a，b归为一段，c，d归为一段，e，f归为一段，一共出现三种不同的数字，即全部类别共分为三段，

Crom矩阵对应的Y变量矩阵如下所示，

Y=123456abcdef]]>

Y变量矩阵是建立PLS模型的因变量的值；

步骤五，对初始种群的每个个体代表的模型进行建模，方法如下：对矩阵的每一行数据，先看每个数字出现的次数，如果存在多种不同的数字而有个数字只出现一次，则这个模型的目标函数值直接设为一个较大的值，接着再判断是不是全部数字都一样的情况，如果是则将所有类别归为一段，建立相应的回归模型，回归方程的自变量矩阵就是全部的校正集样品的红外光谱数据，因变量矩阵就是全部的校正集样品的因变量值，每类样品的因变量的设置值如Y矩阵所示，用所求的回归方程模型求得预测子集样品的因变量的值，从而判断其类别，目标函数是预测子集样品的误判数，如果是其它的情况，则从0开始，用数字相同的类别建立相应的分段回归模型1：包括相应的PLS方程和类别信息矩阵，回归方程的自变量矩阵就是归为同一段的校正集样品的红外光谱数据，因变量矩阵就是这些样品对应的Y变量值，类别信息矩阵就是归为同一段的校正集样品的红外光谱数据做主成分分析后的得分，接着寻找标号为1的数字，重复建模的动作，建立分段回归模型2，依次类推，直到所有的样本都划分完毕，求这个模型的目标函数值，方法如下：先求测试样品数据的得分，再求得分和各个类别信息矩阵的马氏距离，然后把该测试样品数据带入马氏距离最小的分段回归模型，求出预测值，如果该测试样本不属于马氏距离最小的分段回归模型，则该样本定位错误，所有预测子集样品都预测完毕，通过和真实的Y变量的值进行比较，得到测试集样品的误判数，记录测试集样品定位错误数，目标函数值为误判数和样品定位错误数之和，所有个体的目标函数值构成一个列矩阵F；

步骤六，根据目标函数值矩阵F，种群，由遗传算法的相关操作得到下一代种群；

步骤七，判断遗传是否终止，是的话选出这代种群中目标函数值最小的个体，这个就是遗传算法选中的最终分段模型，否的话，回到步骤六；

步骤八，遗传过程结束，构造遗传算法找的的最优的模型；

步骤九，将独立预测集的样本代入最终的模型，算出误判数W1，和普通PLS模型的误判数W2进行比较，如果W1比W2小则接受这个新模型，否则还是采用普通PLS模型。