[发明专利]蚁群算法优化最小二乘支持向量机的溶解氧预测方法

摘要

本发明公开了一种蚁群算法优化最小二乘支持向量机的溶解氧预测方法，其先对原始溶解氧数据序列进行预处理以丢弃或修复坏点数据；然后将预处理后的溶解氧数据序列分成训练集和测试集；接着根据训练集中的所有溶解氧数据，利用蚁群算法优化最小二乘支持向量机的惩罚因子和核函数宽度参数，得到各自的最佳优化值；最后根据惩罚因子和核函数宽度参数的最佳优化值，计算测试集中的每个溶解氧数据的预测值；优点是利用蚁群算法优化最小二乘支持向量机的惩罚因子和核函数宽度参数，并在迭代搜索过程中，对头蚁的位置进行精细搜索获得头蚁的新位置，且保留每只蚂蚁搜索过程中的最佳位置，将蚂蚁的最佳位置引入搜索过程，提高了溶解氧的预测精度。

主权项

1.一种蚁群算法优化最小二乘支持向量机的溶解氧预测方法，其特征在于包括以下步骤：①在需要预测溶解氧的水产品养殖水体中布置一个溶解氧传感器，溶解氧传感器实时等时间间隔连续采集溶解氧数据；②将溶解氧传感器等时间间隔连续采集的T个溶解氧数据按采集时间先后顺序排列构成原始溶解氧数据序列，记为Y，Y＝{y₁,y₂,…,y_t,…,y_T}；再对Y进行预处理，以丢弃或修复Y中的坏点数据而保留Y中的有效数据，得到预处理后的溶解氧数据序列，记为X，X＝{x₁,x₂,…,x_s,…,x_S}；其中，T＞Z，Z表示设定的溶解氧数据最小需求量，1≤t≤T，y₁,y₂,…,y_t,…,y_T对应表示Y中的第1个溶解氧数据、第2个溶解氧数据、…、第t个溶解氧数据、…、第T个溶解氧数据，S≤T，1≤s≤S，x₁,x₂,…,x_s,…,x_S对应表示X中的第1个溶解氧数据、第2个溶解氧数据、…、第s个溶解氧数据、…、第S个溶解氧数据，坏点数据定义为：若Y中的任一个溶解氧数据大于Y中的所有溶解氧数据的平均值的60％或者小于Y中的所有溶解氧数据的平均值的40％，则确定该溶解氧数据为坏点数据；有效数据定义为：若Y中的任一个溶解氧数据大于或等于Y中的所有溶解氧数据的平均值的40％且小于或等于Y中的所有溶解氧数据的平均值的60％，则确定该溶解氧数据为有效数据；所述的步骤②中的X的获取过程为：②_1、找出Y中的所有有效数据和所有坏点数据，并将有效数据的标记设为1；②_2、按序对Y中的每个坏点数据进行处理，将Y中当前待处理的坏点数据定义为当前数据；②_3、如果当前数据为Y中的第1个溶解氧数据或当前数据为Y中的第T个溶解氧数据或时间上先于当前数据的所有溶解氧数据均为坏点数据或时间上后于当前数据的所有溶解氧数据均为坏点数据，则将当前数据的标记设为0，再执行步骤②_4；如果当前数据不为Y中的第1个溶解氧数据且当前数据不为Y中的第T个溶解氧数据且时间上先于当前数据的所有溶解氧数据不全为坏点数据且时间上后于当前数据的所有溶解氧数据不全为坏点数据，则对当前数据进行线性插值以修复当前数据，使当前数据转变成有效数据，假设当前数据为Y中的第t'个溶解氧数据y_t'，那么y_t'经线性插值修复后等于然后将经线性插值修复后的当前数据的标记设为2，再执行步骤②_4；其中，表示时间上先于y_t'且与y_t'最邻近的标记为1的有效数据，亦为Y中的第t₁个溶解氧数据，表示时间上后于y_t'且与y_t'最邻近的标记为1的有效数据，亦为Y中的第t₂个溶解氧数据，1≤t₁≤T，1≤t₂≤T，t₁＜t'＜t₂；②_4、将Y中下一个待处理的坏点数据作为当前数据，然后返回步骤②_3继续执行，直至Y中的所有坏点数据处理完毕，再执行步骤②_5；②_5、丢弃Y中所有标记为0的坏点数据，而保留其余溶解氧数据，得到新的溶解氧数据序列，作为预处理后的溶解氧数据序列X；③判断S≥Z是否成立，如果成立，则执行步骤④；否则，返回步骤②重新获取原始溶解氧数据序列，再进行预处理；④选取X中的前ST个溶解氧数据构成训练集，记为Train，Train＝{x₁,x₂,…,x_st,…,x_ST}；并将剩余的S‑ST个溶解氧数据构成测试集，记为Test，Test＝{x_ST+1,x_ST+2,…,x_S}；其中，S×60％≤ST≤S×80％，x₁,x₂,…,x_st,…,x_ST对应表示Train中的第1个溶解氧数据、第2个溶解氧数据、…、第st个溶解氧数据、…、第ST个溶解氧数据，x₁,x₂,…,x_st,…,x_ST亦对应表示X中的第1个溶解氧数据、第2个溶解氧数据、…、第st个溶解氧数据、…、第ST个溶解氧数据，x_ST+1,x_ST+2,…,x_S对应表示Test中的第1个溶解氧数据、第2个溶解氧数据、…、第S‑ST个溶解氧数据，x_ST+1,x_ST+2,…,x_S亦对应表示X中的第ST+1个溶解氧数据、第ST+2个溶解氧数据、…、第S个溶解氧数据；⑤根据Train中的所有溶解氧数据，利用蚁群算法优化最小二乘支持向量机的惩罚因子和核函数宽度参数，得到最小二乘支持向量机的惩罚因子和核函数宽度参数各自的最佳优化值，具体过程为：⑤_1、选择径向基函数作为最小二乘支持向量机的核函数；然后选择最小二乘支持向量机的惩罚因子c和核函数宽度参数σ作为优化选择参数；再将区间[0,150]作为惩罚因子的值选择范围，将区间[0.1,10]作为核函数宽度参数的值选择范围；其中，k()为径向基函数的表示形式，x为最小二乘支持向量机的核函数的输入变量，D为最小二乘支持向量机的核函数中心，exp()表示以自然基数e为底的指数函数，σ为最小二乘支持向量机的核函数宽度参数，x、D和σ的值均为实数；⑤_2、设定蚁群算法中所需的蚂蚁数量为N；然后从惩罚因子的值选择范围中随机选择N个值分别作为惩罚因子c的值，对应记为c₁,c₂,…,c_p,…,c_N；并从核函数宽度参数的值选择范围中随机选择N个值分别作为核函数宽度参数σ的值，对应记为σ₁,σ₂,…,σ_p,…,σ_N；其中，1＜N≤ST，1≤p≤N，c₁、c₂、c_p和c_N对应为从惩罚因子的值选择范围中随机选择的第1个值、第2个值、第p个值和第N个值，σ₁、σ₂、σ_p和σ_N对应为从核函数宽度参数的值选择范围中随机选择的第1个值、第2个值、第p个值和第N个值；⑤_3、获取蚁群算法第1次搜索时每只蚂蚁的位置，将蚁群算法第1次搜索时第p只蚂蚁的位置记为XA_p,1，XA_p,1＝[c_p,σ_p]；其中，[c_p,σ_p]表示由第p个惩罚因子c的值c_p和第p个核函数宽度参数σ的值σ_p构成的第p只蚂蚁的初始位置矢量；⑤_4、确定蚁群算法第1次搜索时的头蚁及其位置，将蚁群算法第1次搜索时的头蚁的位置记为并确定蚁群算法第1次搜索时每只蚂蚁的最佳位置，将蚁群算法第1次搜索时第p只蚂蚁的最佳位置记为和的获取过程为：a1、根据Train中的所有溶解氧数据和蚁群算法第1次搜索时每只蚂蚁的位置，获取蚁群算法第1次搜索时每只蚂蚁对应的ST个拉格朗日乘子和偏差，将蚁群算法第1次搜索时第p只蚂蚁对应的ST个拉格朗日乘子对应记为a_1,p,a_2,p,…,a_ST,p，将蚁群算法第1次搜索时第p只蚂蚁对应的偏差记为b_p，a_1,p,a_2,p,…,a_ST,p和b_p是根据Train中的所有溶解氧数据和XA_p,1，计算得到的；其中，a2、根据Train中的所有溶解氧数据和蚁群算法第1次搜索时每只蚂蚁的位置，计算蚁群算法第1次搜索时每只蚂蚁对应的Train中的每个溶解氧数据的估计值，将蚁群算法第1次搜索时第p只蚂蚁对应的Train中的第st个溶解氧数据x_st的估计值记为其中，a_i,p表示蚁群算法第1次搜索时第p只蚂蚁对应的第i个拉格朗日乘子；a3、计算蚁群算法第1次搜索时每只蚂蚁的位置对应的预测误差，将蚁群算法第1次搜索时第p只蚂蚁的位置对应的预测误差记为MSE(XA_p,1)，a4、计算蚁群算法第1次搜索时每只蚂蚁的信息素，将蚁群算法第1次搜索时第p只蚂蚁的信息素记为τ(p,1)，τ(p,1)＝exp(‑MSE(XA_p,1))；a5、将蚁群算法第1次搜索时N只蚂蚁的信息素按序构成一个维数为1×N的信息素矩阵，记为In₁，In₁＝[τ(1,1)τ(2,1)…τ(N,1)]；其中，τ(1,1)表示蚁群算法第1次搜索时第1只蚂蚁的信息素，τ(2,1)表示蚁群算法第1次搜索时第2只蚂蚁的信息素，τ(N,1)表示蚁群算法第1次搜索时第N只蚂蚁的信息素；a6、将In₁中值最大的信息素对应的蚂蚁确定为蚁群算法第1次搜索时的头蚁，将In₁中值最大的信息素对应的蚂蚁的位置作为蚁群算法第1次搜索时的头蚁的位置a7、确定蚁群算法第1次搜索时第p只蚂蚁的最佳位置为⑤_5、令g表示蚁群算法的搜索次数，并令g_max表示设定的最大搜索次数；其中，g的初始值为2，2≤g≤g_max，g_max＞2；⑤_6、获取蚁群算法第g次搜索时每只蚂蚁的位置，将蚁群算法第g次搜索时第p只蚂蚁的位置记为XA_p,g，其中，λ表示启发因子，λ∈[0,1]，XA_p,g‑1表示蚁群算法第g‑1次搜索时第p只蚂蚁的位置，表示蚁群算法第g‑1次搜索时的头蚁的位置，表示蚁群算法第g‑1次搜索时第p蚂蚁的最佳位置；⑤_7、确定蚁群算法第g次搜索时的头蚁及其位置，将蚁群算法第g次搜索时的头蚁的位置记为并确定蚁群算法第g次搜索时每只蚂蚁的最佳位置，将蚁群算法第g次搜索时第p只蚂蚁的最佳位置记为和的获取过程为：b1、计算蚁群算法第g次搜索时每只蚂蚁的信息素，将蚁群算法第g次搜索时第p只蚂蚁的信息素记为τ(p,g)，τ(p,g)＝(1‑ρ)×τ(p,g‑1)+τ(p,1)；其中，ρ表示信息素挥发数，ρ∈[0,1]，τ(p,g‑1)表示蚁群算法第g‑1次搜索时第p只蚂蚁的信息素，τ(p,1)表示蚁群算法第1次搜索时第p只蚂蚁的信息素；b2、构建一个维数为g×N的信息素矩阵，记为In_g，In_g中的第h行中的N个信息素对应为蚁群算法第h次搜索时N只蚂蚁的信息素；其中，1≤h≤g；b3、将In_g中的第p列中的g个信息素中值最大的信息素对应的蚂蚁的位置确定为蚁群算法第g次搜索时第p只蚂蚁的最佳位置b4、将In_g中值最大的信息素对应的蚂蚁确定为蚁群算法第g次搜索时的头蚁，将In_g中值最大的信息素对应的蚂蚁的位置作为蚁群算法第g次搜索时的头蚁的位置b5、对进行精细搜索，得到蚁群算法第g次搜索时的头蚁的新位置，记为其中，Cauchy(0,1)为用于产生0～1之间分布的随机数的函数；b6、更新将更新后的位置记为再令其中，表示蚁群算法第g次搜索时的头蚁的新位置对应的预测误差，表示蚁群算法第g次搜索时的头蚁的位置对应的预测误差，中的“＝”为赋值符号；⑤_8、判断g＜g_max是否成立，如果成立，则令g＝g+1，然后返回步骤⑤_6继续执行；否则，确定蚁群算法搜索结束，得到最终的头蚁的位置，记为XA^obj，然后将XA^obj对应的位置矢量[c^obj,σ^obj]中的c^obj作为最小二乘支持向量机的惩罚因子c的最佳优化值，并将XA^obj对应的位置矢量[c^obj,σ^obj]中的σ^obj作为最小二乘支持向量机的核函数宽度参数σ的最佳优化值；其中，g＝g+1中的“＝”为赋值符号，表示蚁群算法第g_max次搜索时的头蚁的位置；⑥根据最小二乘支持向量机的惩罚因子c的最佳优化值c^obj和核函数宽度参数σ的最佳优化值σ^obj，计算Test中的每个溶解氧数据的预测值，将Test中的第r个溶解氧数据的预测值记为其中，1≤r≤S‑ST，x_ST+j表示Test中的第j个溶解氧数据，x_ST+r表示Test中的第r个溶解氧数据，a_j表示Test对应的第j个拉格朗日乘子，b表示Test对应的偏差，a_j和b是根据Test中的所有溶解氧数据及c^obj和σ^obj，计算得到的，