一般非线性预测模型参数优化的数值方法
更新时间:2023-04-22 05:12:01 阅读量: 实用文档 文档下载
- 非线性模型参数估计推荐度:
- 相关推荐
数值计算或是参数优化
《预测》1997年第1期
理论与方法研究
Ξ
(河海大学 210098) (摘 要 。
关键词 1 引言
;任意实常数q为1时即为最小乘准则,为2时为最小二乘准则,等等,视实际预测要求而定。
加速基因算法包括如下8步。
步1:模型参数变化空间的离散和二进制编码。设编码长度为l,把模型每个参数的变化区间等分成2l-1个子区间,于是模型参数变化空间被离散成(2l)p个
非线性预测模型参数优化(估计)的传统方法存在两方面的不足,一是对模型结构(如可线性化、可微)作了限定[1,2],二是方法受参数优化准则形式制约
(如目标规划法、预松驰算法)[3],而对实际预测环境、
参数格网点。GA中称每个格网点为个体,它对应模型
p个参数的一种可能取值状态,并用p个l二进制数表
示。这样,模型p个参数的取值状态、格网点、p个二进制数建立了一一对应关系。算法的模型参数优化工作将主要对这些二进制数进行操作运算。
步2:初始父代个体群的随机生成。从上述(2l)p个格网点中均匀随机选取n个点作为初始父代个体群。
步3:父代个体的适应能力评价。把第i个个体代入式(1)得相应的优化准则值fi,fi越小则该个体的适应能力越强。
步4:父代个体的概率选择。把已有父代个体按优化准则值fi从小到大排序。称排序后最前面的几个个体为优秀个体(群)。构造与fi值成反比的函数pi且满以概率pi选择第i个个体,共选择两组各n个个体。
步5:父代个体的杂交。由上得到的两组个体两两配对成为n对双亲。将每对双亲的二进制数的任意一段值互换,得到两组子代个体。
步6:子代个体的变异。任取上步中的一组子代个体,将它们的二进制数的任意一段值依某概率(即变异率)进行翻转(原值为0的变为1,反之变为0)。
步7:进化迭代。由上步得到的n个子代个体作为新的父代,算法转入步3进入下一次进化过程,重新评价、选择、杂交和变异,如此循环往复,使优秀个体逼近最优点。
测目标缺乏足够的考虑,致使其适用性不强,影响了实际预测效果。基于自然选择和自然基因机制的基因算法(GeneticAlgorithm,简称GA)突破了对模型结构是否线性、连续、可微等各种限制,同时可不受优化准则形式、优化参数数目、约束条件等的束缚,直接在优化准则值的引导下进行全局自适应参数优化,简便实用,是目前处理参数优化问题的一种通用的数值方法。这种方法已引起工程界的广泛应用和研究
[4,5]
。
我们在应用GA时发现,GA对各种实际预测模型参数变化区间的大小变化的适应能力较差,计算量但利用优秀个体逐步调整模型参数变化区间,可形成一种称为加速基因算法(AccelerationGenetic
Algorithm,简称AGA)的新方法。经大量数值试验和
大,算法自身参数设置技术目前尚无明确准则指导。足pi>0和p1+p2+…pn=1。从这些父代个体中
实际应用,表明AGA对GA有明显改进。下面给出
AGA的原理、步骤、算法参数设置准则。
2 加速基因算法的原理、步骤及其参数设置
不失一般性,设一般非线性预测模型为p-1阶多项式,模型参数优化准则为极小化下式
f=
6
m
i=1
-1
c1+c2xi+…+cpxp-yi qi
(1)
式中,{cj}为p个待优化参数;{(xi,yi)}为m对自变
Ξ
收稿日期:1996-11-11
49
数值计算或是参数优化
以上7步构成传统的基因算法(GA)。
步8:加速循环。大量研究、应用表明,GA中的杂交、变异的寻优功能随进化迭代次数的增加而迅速减弱,在实际应用中常出现在远离全局最优点的地方
GA即停滞寻优工作;GA的计算量大,全局优化速度
这些在加速基因算法中却是不必要的。实际上,由于
AGA利用进化迭代过程中产生的优秀个体所包含的
优化准则函数在最优点附近各参数方向的变化特性的信息来调整模型参数变化空间的大小,使AGA能同时在p个参数方向寻优且收敛,nsnAGA对实。经n、ns的配置应满(l 2l)ns n>n
(2)
慢,GA是从(2l)p种模型参数取值状态中找最优,优化结果的精度受编模长度l的控制;GA自身参数设置技术复杂,尚无好的准则指导。特别是当实际预测模型参数的变化区间很大时,上述问题就十分突出GA就极为困难。
受正交设计法[6],GA一、,算法转入步1,如此循环往复,优秀个体参数变化空间将逐步收缩,与最优点的距离越来越近,直至最优个体的优化准则值小于某一设定值或达到预定加速(循环)次数,结束整个算法的运行。
以上8步构成加速基因算法(AGA)。其中的算法参数有编码长度l、父代个体数目n、优秀个体数目ns和变异率ca。
从AGA角度看,ca反映了个体向其它参数格网点随机变迁的概率。因AGA的全局优化能力比GA强,
ca越大,优秀个体包围,接近最优点的机会越大,ca一
在实际应用时n一般取300或以上,ns相应取10或以上。
笔者对文献〔1〕、〔3〕、〔7〕中的指数模型、双曲线模型、威市尔分布模型、逻辑斯蒂曲线等模型的参数估计应用AGA法来处理,经加速循环5~10次后拟合效果均好于原实例使用的常规方法(具体过程、结果略),这说明,AGA法的确具有优化参数选择的实用功能。
(本文得到河海大学水环院朱元生生教授、刘权授
教授和四川联合大学水电学院丁晶教授的热情指点,特此深表谢意!)参 考 文 献
[1]唐五湘.指数曲线模型参数估计的两步最小二乘
法.预测,1994(6)
[2]陈性生.一类非线性回归模型LS估计的两步逼近
般取定1.0。
l越大参数格网越细,格网点就越多,要求父代个
求法.预测,1990(2)
[3]袁修贵.“残差绝对值和最小”准则的目标规划法.
体数目也须相应增大,导致计算量猛增。因AGA在算法步8中能逐步缩小模型参数的变化空间,随着AGA的运行参数格网自动分细,因此l常取定10。
由AGA原理可知,n个父代个体是按优化准则从
(2)个参数格网点中的一个随机抽样(相当于一组广
lp
预测,1994(1)
[4]GoldbergDE.geneticalgorithmsinsearch,opti2mizationandmachinelearning.1989
[5]WangQJ.thegeneticalgorithmanditsapplica2
tiontocalibratingconceputalrainfall_runoff.waterresourcesresearch,1991,27(9):models2467-2471
[6]金良超.正交设计多指标分析.北京:中国铁道
addison2wesly,
义的正交试验),n太小就不能充分探测优化准则函数在最优点附近的足够信息;当n增大而优秀个体数目
ns仍维持较小值不变时,优秀个体包围、接近最优点
的机会也会减少;AGA的计算量与n、ns都成正比,
n、ns又不能取得很大。
22
以往的基因算法认为,当模型参数变化空间很大时,GA的编码长度、父代个体数目也须很大,导致计算量大增,而且算法参数的设置技术也趋复杂化。而
出版社,1988
[7]司马锡生.
(1)
皮尔曲线的参数估计.预测,1994
50
正在阅读:
一般非线性预测模型参数优化的数值方法04-22
五年级英语下册Unit2《SpecialDays》(Lesson2)教学设计人教新起04-12
律师工作心得体会03-31
论语上说:师者,所以传道、授业、解惑者也09-12
2014年二级建造师《建筑工程管理与实务》考前冲刺试卷(第一部分)09-08
提高水泥混凝土路面平整度的修筑技术04-19
一年级造句09-18
温暖的中秋作文500字07-01
小组合作学习制在班级管理中的应用05-12
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 非线性
- 数值
- 模型
- 优化
- 参数
- 一般
- 预测
- 方法
- 2022MPACC会计硕士考研必须知道的内容,会计硕士mpacc论坛精选
- 高校共青团工作的定位
- 形象设计形式美法则
- 油漆涂料颜色的基础知识
- The role of phytochelatins in constitutive and adaptive heav
- 德国旅行注意事项,德国自助游旅游攻略
- 2.2线性代数第二章张第2节
- 人教版七年级上《一元一次方程》期末复习知识点+检测试卷3
- 常微分方程的数值解法
- 六年级奥数题及答案周汇总
- _天人合一_思想及其当代启示
- 2015年4月入党积极分子思想汇报
- 动物体的结构层次PPT
- 电路课件_第10章(第五版_邱关源_高等教育出版社)
- 2014辽宁省驾校考试科目一C2考试答题技巧
- 防火、防盗应急预案
- 2013年东营市初中学生学业考试英语试题
- 2015年4月1日出生的宝宝生辰八字五行缺什么喜用神起名分析
- 人力资源部部门职责规定
- 固定高清雷达测速电子警察