一种卫星电源主母线电流区间预测方法
【专利摘要】本发明公开了一种卫星电源主母线电流区间预测方法，该方法基于优化的核极限学习机训练的预测模型，利用比例系数法确定预测区间，其中比例系数法的参数由差分进化算法进行优化，具体过程包括：卫星电源主母线电流数据预处理，去除噪声数据并获取归一化数据；采用差分进化算法优化核极限学习机参数；使用优化的核极限学习机建立初始预测模型；给定评估预测区间质量的综合指标，采用比例系数法确定预测区间，评估预测区间满意度；最后，利用差分进化算法优化区间预测比例系数，获得最优的卫星电源主母线电流预测区间。本发明的优点是：针对复杂的卫星电源主母线电流数据，给出了一种卫星电源主母线电流区间预测方法，该方法预测精度更高，效果更好。
【专利说明】
一种卫星电源主母线电流区间预测方法
技术领域
[0001] 本发明涉及一种卫星电源主母线电流区间预测方法，该方法基于优化的核极限学 习机训练的预测模型，利用差分进化算法优化的比例系数法确定预测区间，本发明属于航 天科学与计算机科学的交叉领域。
【背景技术】
[0002] 卫星是一种结合遥感、通信等技术研制出来的大型多功能复杂系统，它是人类探 索宇宙世界的主要途径。卫星电源系统是系统的功能子系统，它为卫星所有设备正常工作 提供能量，主要完成产生电能、储存能量、变换和调节能量、传输分配电能等任务。它的工作 原理是通过物理或化学变化，将太阳能、核能等转化成电能，并完成对电能的存储、变换和 调节，然后将电能传输分配给卫星其他系统。
[0003] 目前常用的卫星电源系统结构为太阳能/蓄电池供电系统：它由供能设备和控制 设备组成，其中太阳电池阵、蓄电池组是供能设备；电池充电控制器、电池放电控制器、分流 调节器、升压调节器等是控制设备。太阳电池阵是卫星电源系统的发电装置，当卫星处于光 照期时，太阳电池阵经过光伏效应把太阳能转换成电能，经过分流调节器进行分流，经过升 压调节进行稳压调节，通过电池充电控制器对蓄电池进行充电，并向卫星各负载设备供电。 蓄电池是卫星电源系统的供电和储能设备，当卫星处于阴影期时，蓄电池经过电池放电调 节器向卫星各负载设备供电。
[0004]控制设备是电源系统的核心部件，其中主母线的优劣直接影响卫星电源系统的供 电品质。主母线的功能是在光照期用于调节太阳电池阵输出功率，在阴影期用于控制和保 护蓄电池放电过程，为卫星电源系统和其他设备提供与地面连接的接口，接受并执行遥控 指令及地面有限指令，采样并预处理遥测信号等。
[0005] 虽然卫星在轨异常发生概率很小，但是一旦发生，极有可能中断卫星工作任务甚 至导致卫星陨落，造成的损失是致命的。例如通信卫星发生故障，将导致信号中断，通信设 备无法继续使用;又如气象卫星出现工作异常，将导致气象工作人员不能及时、有效的获得 气象信息；导航卫星发生故障将会导致地面接收端无法获得地理位置信号，进而无法对地 理位置做出准确的判断。由此可见，卫星故障及异常状况给国家和普通民众带来的影响和 损失是不可估量的。调查发现，在种种卫星故障和异常中，电源系统故障是主要因素，电源 系统故障对于整个卫星的工作有严重的影响，因此，有必要对卫星电源系统的关键参数如 主母线电流进行准确预测，及时发现和避免卫星故障。
[0006] 通常情况下，及时发现和避免卫星电源系统故障的有效做法是对主母线电流的未 来趋势进行预测，通过主母线电流的预测值和其正常范围的对比来判断卫星电源主母线是 否可能出现异常，本发明的目的是通过对卫星电源主母线电流数据进行区间预测模型的训 练，利用该区间预测模型获得卫星电源主母线电流的预测值及预测区间，从而实现对卫星 电源主母线电流未来趋势的预测，以此来及时发现和避免卫星电源主母线的异常和故障。
【发明内容】
[0007] 本发明的目的是基于卫星电源系统的主母线电流数据，发明一个对卫星电源主母 线电流区间预测的方法，基于该方法，可以预测卫星电源主母线电流的预测值以及预测区 间，从而解决卫星电源主母线电流数据异常监测等问题，实现对卫星电源系统主母线可疑 异常的及时发现和避免，为卫星在轨运行的辅助控制管理打下坚实基础。
[0008] 本发明的具体技术方案包括以下几个步骤：
[0009] 步骤一:针对卫星电源主母线电流数据，对原始数据进行清洗、压缩、转换和归一 化处理。
[0010] 步骤二:借助差分进化算法优化核极限学习机参数，确定最优的核极限学习机的 核参数以及惩罚系数。
[0011] 步骤三:借助步骤二中得到的核极限学习机的参数，结合比例系数法建立核极限 学习机预测模型。
[0012] 步骤四：借助步骤三中得到的原始预测模型，结合集成方法和差分进化算法优化 比例系数法建立卫星电源主母线电流预测模型。
[0013] 本发明的有益效果是：
[0014] 本发明所提的方法从数据挖掘的角度出发，针对卫星电源主母线电流数据本身进 行特征提取与预测分析，提出了一种新颖的预测方法，有效地解决了卫星电源主母线电流 数据噪声数据多、维度高、分析难度高的问题。本发明经过扩展后，还可以用于航天器其他 遥测数据的预测，为在轨卫星异常检测与分析提供了一种新方法。
【附图说明】
[0015] 图1是本发明方法的总体流程图。
【具体实施方式】
[0016] 下面结合附图和相关算法，对本发明做进一步的说明。
[0017] 本发明的总体流程如图1所示。
[0018] 本发明针对卫星电源主母线电流数据设计了一种主母线电流区间预测方法，借助 数据预处理方法去除数据中的噪声数据，提取归一化数据。借助差分进化算法优化核极限 学习机参数，确定最优的核极限学习机的核参数以及惩罚系数，降低预测误差，提高预测模 型精度。针对优化后的核极限学习机参数，结合比例系数法建立原始预测模型，针对原始预 测模型，借助集成方法和差分进化算法优化比例系数法建立最终预测模型。具体实施步骤 如下：
[0019] 1.数据预处理
[0020] 针对卫星电源遥测数据，首先将数据按照一分钟进行分段，将每分钟超出正常范 围30倍且出现次数小于3次的数据作为噪声数据并将其去除；随后对数据进行统计分析，将 去噪后的主母线电流数据进行压缩;再将单维的时间序列转换成矩阵的形式，获得适合核 极限学习机的样本集;最后对样本集矩阵按列分别进行归一化处理。其过程可用算法1描 述：
[0021] 算法1:数据预处理
[0022] 输入:数据集D，嵌入维度m
[0023]输出：新数据集ND
[0024] 1 :temp_D =临时数据集
[0025] 2: count =主母线电流数据数目
[0026] 3:minute_number =主母线电流数据按分钟分段数
[0027] 4:将主母线电流数据temp_D按照分钟进行分段，分段结果存放在Seg_D*
[0028] 5：for(i = 0； i<minute_number ； i++)
[0029] 6:统计seg_D[ i ]中超出正常范围30倍的数据出现次数
[0030] 7:end for
[0031] 8:将每分钟超出正常范围30倍且出现次数小于3次的数据去除
[0032] 9:对数据进行统计分析，压缩主母线电流数据
[0033] 10 ：for( j = 0 ； j<count-m； j++)
[0034] 11:将单维的时间序列转换成矩阵matrix的形式
[0035] 12:end for
[0036] 13 ：for(k = 0 ；k<m；k++)
[0037] 14:ND[k]=样本集矩阵按列进行归一化处理后的数据
[0038] 15:end for
[0039] 2.借助差分进化算法优化核极限学习机参数
[0040] 针对经过数据预处理的数据集ND，本发明借助差分进化算法优化核极限学习机参 数，确定最优的核极限学习机的核参数以及惩罚系数，降低预测误差，提高预测模型精度。
[0041] 首先，初始化差分进化算法参数和个体基因序列;然后计算目标函数适应值，并选 择出最佳适应值和最差适应值。随后经过变异、交叉、选择操作，保留最优基因到下一代， 检查是否满足终止条件，若满足则结束，否则返回继续计算新种群的目标函数适应值、最佳 适应值、最差适应值，继续经过变异、交叉、选择操作，直到满足终止条件为止。其具体过程 如算法2描述：
[0042]算法2:差分进化算法优化核极限学习机参数
[0043]输入:种群个数NP，搜索空间维数D，最大迭代次数tmax，最大缩放因子Fmax，最小缩 放因子Fmin，最大交叉概率CRmax，最小交叉概率CRmin
[0044] 输出：全局最优解〇best,G，Cbest,G
[0045] l：G = 0,Fmax=l,Fmin = 0,CRmax=l,CRmin = 0.5
[0046] 2：for(i = 0：i<NP；i++)
[0047] 3: crmin …d
[0048] Cmm = {c\，c2 ,…，},Cmax = {c】，c2，…，cD }
[0049] 4:根据公式(4)初始化种群
[0050] 5： end for
[0051] 6:while(t<tmax| | 目标函数值收敛）
[0052] 7：for(i = 0； i<NP； i++)
[0053] 8:计算目标函数适应值
[0054] 9： end for
[0055] 10:选择最佳适应值、最差适应值
[0056] ll：for(i = 0；i<NP； i++)
[0057] 12:根据公式(7)更新缩放因子
[0058] 13:然后选择公式(8)_(13)中的一个公式根据原始个体产生一个新的变异个体
[0059] 14:判断变异个体是否满足边界条件，如果不满足边界条件，则通过公式(14)修改
[0060] 15:end for
[0061 ] 16:根据公式(15)更新交叉概率CR，
[0062] 17：for(i = 0；i<NP；i++)
[0063] 18:九-二卜酬机】)* D」
[0064] 19：for( j = 0； j<D； j++)
[0065] 20:根据交叉概率和公式（16)从原始个体和变异个体中选择一个作为新的交叉个 体
[0066] 21:end for
[0067] 22：end for
[0068] 23：for(i = 0；i<NP；i++)
[0069] 24:采用贪心算法选择出最优的个体进入下一代
[0070] 25：end for
[0071] 26:t-1+1
[0072] 27:end while
[0073] 28:return 〇best,G，Gbest,G
[0074] 3.建立核极限学习机预测模型
[0075]利用步骤2中得到的核极限学习机参数，本发明结合比例系数法建立核极限学习 机预测模型。首先，利用优化后的核极限学习机参数和经过预处理的数据建立预测模型，将 得到的预测值使用比例系数法确定预测区间。其具体实现过程如算法3、4所描述：
[0076]算法3:建立核极限学习机模型
[0077]输入:优化后的核极限学习机参数〇和(：，主母线电流数据ND，数据维数M，N [0078]输出：主母线电流预测值Y
[0079] 1 :K(x，y) =exp(_ | | x_y | 12/2〇2)
[0080] 2：for(i = 0； i<M； i++)
[0081] 3：for( j = 0； j<N； j++)
[0082] 4: Q i，j = h(xi) ? h(Xj)=K(Xi，Xj)
[0083] 5：end for
[0084] 6：end for
[0085] 7： ^ = /fr(Q￡LM+//〇-1r AT
[0086] 8： /(x) = h(x)p = h(x)HT(HHT + //C)~lT = ??? +1/CylT K{x,xn)
[0087] 9：Y = f(x)
[0088]算法4:确定预测区间
[0089] 输入:主母线电流预测值Y，预测值数目count，预测值变化范围R，比例参数a和0
[0090] 输出：主母线电流预测区间[L，U]，预测区间覆盖率PICP，预测区间平均宽度百分 比PIMWP，平均带宽误差MWD
[0091] l:u = Y*(l+a)
[0092] 2:L = Y*(l-0)
[0093] 3：SPICP = 0,SMWP = 0,SMWD = 0
[0094] 4：for(i = l ； i<count； i++)
[0095] 5：if(Y[i]^L[i]&&Y[i]^U[i])
[0096] 6：CP[i] = l
[0097] 7：else
[0098] 8：CP[i]=0
[0099] 9：end if
[0105] 15:end if
[0106] 16:end for
[0107] 17：for(i = l； i<count； i++)
[0108] 18：SPISP = SPICP+CP[i]
[0109] 19：SMWP = SMffP+MWP[i]
[0110] 20：SMWD = SMffD+WD[i]
[0111] 21：end for
[0115] 4.建立卫星电源主母线电流预测模型
[0116]本发明选取卫星电源主母线电流数据特征向量集，训练得到核极限学习机预测模 型，并结合集成方法和差分进化算法优化比例系数法对原始模型进行优化得到最终预测模 型，利用该模型对卫星电源主母线电流进行区间预测，判断是否出现数据异常。
[0117]其本质是将卫星电源主母线电流数据的特征向量按照5: 1:1的比例分割为训练 集、验证集和测试集，在训练集上训练k个核极限学习机模型，在验证集上使用比例系数法 进行事后预测与评估，根据评估结果使用差分进化算法优化比例系数，得到最优的比例系 数，同时应用集成技术集成k个核极限学习机的优化结果，将k组比例系数的平均数作为最 终的比例系数，建立稳定的预测模型。其具体实现过程如算法5、6所描述：
[0118]算法5:差分进化算法优化比例系数
[0119]输入:种群个数NP，搜索空间维数D，最大迭代次数tmax，最大缩放因子Fmax，最小缩 放因子Fmin，最大交叉概率CRmax，最小交叉概率CRmin
[0120] 输出：全局最优解abest,G，0best,G
[0121] l：G = 0,Fmax=l,Fmin = 0,CRmax=l,CRmin = 0.5
[0122] 2：for(i = 0；i<NP；i++)
[0123] 3: amm ={?…，d〇rmax ={?...尤}，
[0124] 1={心心...,〇 1={心心???,<}
[0125] 4：for(j = 0；j<D；j++)
[0126] 5： < = aJmm + ra?J(0,l).(<axj = \2,…,D
[0127] 6： end for
[0128] 7:end for
[0129] 8:while(t<tmax| 旧标函数值收敛）
[0130] 9：for(i = 0；i<NP； i++)
[0131] 10:F(a[i]，机i]) = l-[l-(l+2 ? MWD) ? PBWP ? (1+e-50(PICP-90))]
[0132] 11：end for
[0133] 12:选择最佳适应值、最差适应值
[0134] 13：for(i = 0；i<NP；i++)
[0136] 15:选择公式(8)_(13)中的一个公式根据原始个体产生一个新的变异个体
[0137] 16:判断变异个体是否满足边界条件，如果不满足边界条件，则通过公式(14)修改
[0138] 17：end for
[0140] 19：for(i = 0；i<NP；i++)
[0141] 20:九w二卜邱0,】)*乃」
[0142] 21：for(j = 0； j<D； j++)
[0143] 22:根据交叉概率和公式（16)从原始个体和变异个体中选择一个作为新的交叉个 体
[0144] 23：end for
[0145] 24：end for
[0146] 25：for(i = 0；i<NP；i++)
[0147] 26:采用贪心算法选择出最优的个体进入下一代
[0148] 27：end for
[0149] 28:t-1+1
[0150] 29：end while
[0151] 30:return abest,G，0best,G
[0152] 算法6:利用集成技术建立最终预测模型
[0153]输入:k个核极限学习机模型的比例系数法的优化比例系数，主母线电流测试集预 测值TY，测试集数目count，预测值变化范围R
[0154] 输出：预测区间[TL，TU]，预测区间覆盖率PICP，预测区间平均宽度百分比PMWP， 平均带宽误差MWD
[0155] 1 :sum_a = 0，sumJ3 = 0
[0156] 2：for(i = l；i<k； i++)
[0157] 3:sum_a=sum_a+a[i]
[0158] 4:sum_0 = sum_0+0[i]
[0159] 5： end for
[0162] 8:TU = TY*(l+a)
[0163] 9:TL = TY*(l-0)
[0164] i〇:sPICP = 0，SMWP = 0，SMWD = 0
[0165] ll：for(i = l； i<count； i++)
[0166] 12:if(TY[i]彡TL[i]&&TY[i]彡TU[i])
[0167] 13：CP[i] = l
[0168] 14：else
[0169] 15：CP[i]=0
[0170] 16:end if
[0174] 20：else if(TY[i]>TU[i])
[0176] 22：end if
[0177] 23：end for
[0178] 24：for(i = l ； i<count ； i++)
[0179] 25：SPISP = SPICP+CP[i]
[0180] 26：SMWP = SMWP+MffP[i]
[0181] 27：SMWD = SMWD+ffD[i]
[0182] 28：end for
【主权项】
1. 一种卫星电源主母线电流区间预测方法，其特征包括如下步骤： (1) 针对卫星电源主母线电流数据的特点，对原数据进行预处理； (2) 借助差分进化算法优化核极限学习机参数，确定最优的核极限学习机的核参数以 及惩ill系数； (3) 借助经过差分进化算法优化的核极限学习机的参数，结合比例系数法建立核极限 学习机的预测模型； (4) 借助集成技术和差分进化算法优化比例系数法建立卫星电源主母线电流预测模 型。2. 根据权利要求1所述的一种卫星电源数据区间预测方法，其特征在于，步骤（1)针对 卫星电源主母线电流数据的预处理，其实现方法包括： (21) 将每分钟超出正常范围30倍且出现次数小于三次的数据认定为噪声数据，将其删 除； (22) 提取原始数据中的特征数据，并进行压缩； (23) 根据给定的嵌入维度，将单维的时间序列转换成矩阵的形式，获得适合核极限学 习机的样本集矩阵； (24) 对样本集矩阵分别按列进行归一化处理。3. 根据权利要求2所述的方法，其特征在于，根据给定的嵌入维度，将原始数据的单维 数据转换成矩阵形式，并对矩阵分别按列进行归一化处理，其实现步骤包括： (31) 确定嵌入维度 针对卫星电源主母线电流数据某时刻输入值xt和其输出值yt存在某种函数关系： yt = f (xt)= f (yt-1，yt-2，…，yt-m) (1) 其中，m是嵌入维度，其值决定时间窗口大小即输入向量的维度。 借助自定义的参数设置方法选择一个合适的嵌入维度，分别进行10组不同的实验，设 置嵌入维度分别为1~10,核函数参数为10,惩罚系数为10,以均方根误差为评价指标，得到 可靠性和准确性最高的嵌入维度。 均方根误差计算公式为：(32) 数据转换 根据已确定的嵌入维度，将经过去噪和压缩的数据xt转换为如下矩阵的数据格式： T={(Xl,Yl),---,(Xi,Yi),---,(Xn-m+l,Yn-m+l)}e(XXY) n-m+1 (3) (33) 归一化处理 针对卫星电源主母线电流数据复杂性，对数据转换得到的样本集矩阵T进行归一化处 理。通过公式(2 )，将样本集矩阵归一化到[0，1 ]范围内，其中，Xupper和Xl_r分别代表归一化上界和下界，Xmax和Xmin分别代表主母线电流的最大 值和最小值。 4 .根据权利要求1所述方法，其特征在于，借助差分进化算法优化核极限学习机参数， 确定最优的核极限学习机的核参数以及惩罚系数，其实现步骤包括： (41) 初始化 初始化差分进化算法参数和个体基因序列：其中，A? = Uim，4…4}和Aax = d，4ax . ?}为参数变量的边界，r and (0，1)产 生0到1之间的均匀随机数。 (42) 计算目标函数适应值，并选择出最佳适应值和最差适应值。 采用k折交叉检验法来评估每个个体的目标函数适应值。首先将训练数据分成大致相 等且互不重叠的k个子集。随机选择其中k-Ι个子集作为训练集，剩下的一个子集作为测试 集，采用均方根误差作为衡量模型预测效果的指标，循环该过程k次，直到每个子集都轮流 作为过测试集一次且仅一次。每个子代的适应值取这k次均方根误差的平均值：选择最小目标函数适应值作为最佳适应值，最大目标函数适应值作为最差适应值。 (43) 变异 首先，使用自适应的方法来调节缩放因子F，其值随目标函数值改变而变化的趋势为：其中，fitneSS(Xi,G)表示第G代个体i的目标函数值，fitnesS(X best,G)表示第G代最优目 标函数值，行加6^仏。^,(；)表示第6代最差目标函数值(^_和?_分别代表?的最大值和最 小值，由于F通常限定在[0，1]之间，所以Fmax=l，Fmin = 0.1。 从以下几种常用的变异策略中选择一种产生一个新的变异个体。其中，F表示缩放因子，Xu表示第G代个体i当前的基因序列，Xb(3St,c表示第G代种群当前 最优个体的基因序列，心是在[1，NP]之间不等于i且互斥的随机整数， 表示第G代种群个体i的第个基因片段。 最后，判断新生成的变异个体是否满足边界条件，如果不满足边界条件，则通过公式 (14)修改其值，使其合法有效。 Vi,G = Xmin+rand(0，l) · (Xmax-Xmin) (14) (44) 交叉 首先，使用一种基于对数函数的方法调节交叉因子，计算公式为：其中，G是当前的迭代次数;max_iterations表示最大的迭代次数;a，b是调节参数，参 数a确保CR大于0,参数b用于放大当前迭代次数和最大迭代次数的比值;CRmax和CR_分别代 表CR的最大值和最小值，通常取值CR max = 1，CRmin = 0.5。 然后，根据交叉概率，从原始个体和变异个体中选择一个作为新的交叉个体：其中，CR表示交叉概率，用于控制基因复制变异基因的概率;ran山[0，1)表示第j个基因 片段随机差生的交叉概率，数值在[〇，1)之间；jrand代表取值为[1，2，一，0]的随机整数。 (45) 选择 采用贪心算法选择出最优的个体进入下一代，如果交叉个体仏^的目标函数适应值小 于或等于相应的原种群个体Xi c的目标函数适应值，则交叉个体C将取代原种群个体c， 进入下一代;否则，原种群个体L c将继续保留到下一代中。相应的选择公式为：(46) 检查是否满足终止条件 检查是否满足终止条件，若满足则结束，否则返回步骤(42)。此迭代过程的结束条件是 达到最大迭代次数或者得到收敛的目标函数值，即满足如下条件： f itness(Xi,c+i)-f itness(Xi,c) | ^Ξ1 .Oe-5 (18) 〇5.根据权利要求1所述方法，其特征在于借助经过差分进化算法优化的核极限学习机 参数，结合比例系数法建立原始的预测模型。其实现方法包括： (51)本发明定义由核极限学习机算法建立预测模型，其实现过程为： 根据差分进化算法优化得到的参数，设置核函数参数σ和惩罚系数C。 应用Mercer条件为核极限学习机算法定义一个核矩阵Ω并计算其值：其中，核矩阵Ω替代ELM算法中的随机矩阵HHT，h(x)为从样本空间映射到隐藏层空间的 函数，核函数K(x，y)将所有的输入样本从η维输入空间映射到高维隐藏层空间，通过差分进 化算法确定核参数之后，核矩阵Ω的映射值为定值。 应用矩阵运算，计算核极限学习机的输出权值：其中，I是单位矩阵;C是惩罚系数，用于权衡结构风险和经验风险之间的比例，解决矩 阵HHT不可逆的问题，对角矩阵HHT可以通过核函数映射产生。 根据已经计算出的核矩阵和输出权值，计算样本输出值：则输出权值为(52)比例系数法 比例系数法是将预测模型得到的预测值分别放大α倍或缩小β倍作为预测区间的上界 和下界。其中，是第i个输入量对应的预测值，U4PU分别是第i个输入量对应的预测上界和下 界，α和β是未知的比例系数。 为了获得最优的预测区间，将卫星电源主母线电流数据分割为训练集、验证集和测试 集，通过优化的核极限学习算法在训练集上建立预测模型，结合比例系数法构建验证集的 预测区间，应用评估指标对这一预测区间进行综合评估，在验证的过程中不断调整参数α和 β的值，直到获得最好的评估结果为止，此时α和β的值即为最优值，这一过程可以转换为求 解如下最优问题：PICP是预测区间覆盖率，用于评价预测区间的准确度和可靠性：其中，Ν表示预测样本数，0&是布尔值，表示第i个样本是否落入预测区间内。其中，^是样本真实值，U和仏是预测区间的上下限。 PIMWP是预测区间平均宽度百分比，反映预测区间的清晰度：其中，N表示预测样本数，R是预测目标值的变化范围，"是样本真实值，U和仏是预测区 间的上下界。 MWD是预测区间平均带宽误差，用于衡量预测区间的误差和精确度：其中，N表示预测样本数，"是样本真实值，^和仏是预测区间的上下界，WDi是样本的带 宽误差。 这一优化问题可以采用步骤(61)差分进化算法来求解，目标函数最小值对应的α和β的 值即为最优值，再根据公式(23)确定预测区间。6.根据权利要求1所述方法，其特征在于采用集成技术和差分进化算法优化比例系数 法建立最终的预测模型。其实现方法包括： (61) 采用差分进化算法优化比例系数法的参数 借助差分进化算法优化比例系数法的参数，确定最优的比例系数，其实现步骤如下： a) 初始化 初始化差分进化算法参数并根据公式(5)初始化种群个体基因序列。 b) 计算目标适应值 根据公式(24)计算目标函数适应值，选择最小目标函数适应值作为最佳适应值，最大 目标函数适应值作为最差适应值。 c) 变异 首先根据公式(7)更新缩放因子，然后选择公式(8)-(13)中的一个公式根据原始个体 产生一个新的变异个体，再判断这个变异个体是否满足边界条件，如果不满足边界条件，则 通过公式(14)修改其值，使其合法有效。 d) 交叉 首先根据公式（15)更新交叉概率CR，然后根据交叉概率通过公式（16)从原始个体和变 异个体中选择一个作为新的交叉个体。 e) 选择 采用贪心算法选择出最优的个体进入下一代，如果交叉个体仏^的目标函数适应值小 于或等于相应的原种群个体Xw，则交叉个体Uw将取代原种群个体Xw，进入下一代;否则， 原种群个体Xw将留在下一代中。 f) 检查是否满足终止条件 检查是否满足终止条件，若满足则结束，否则返回b)。本迭代过程的结束条件是达到最 大迭代次数或者得到收敛的目标函数值。 (62) 借助集成技术建立最终预测模型 为增加预测模型的稳定性并且提高预测精度，将实验数据分割为训练集、验证集和测 试集，通过集成技术使用训练集分别训练k个核极限学习机模型，在验证集上采用差分进化 优化的比例系数法分别确定k组预测区间参数，取k组预测区间参数的平均值作为预测模型 的最终比例系数，建立最终预测模型。
【文档编号】G01R19/00GK105929216SQ201610223421
【公开日】2016年9月7日
【申请日】2016年4月8日
【发明人】皮德常, 康旭
【申请人】南京航空航天大学