发明名称 基于Lagrange函数的最小二乘准则多目标优化方法
摘要 本发明公开了一种基于Lagrange函数的最小二乘准则多目标优化方法,包括:S1、将多目标规划转化为单目标规划;S2、对建立的模型拟合并得到拟合误差;S3、得到最优解模型;S4、将加权因子β介入到目标函数Q(n)中;S5、得到扩展函数;S6、在最优解模型中去掉非负权重;S7、采用最小二乘准则计算最优加权系数。本发明的方法不管系统初始状态如何,系统都可以演化到平衡态;亦即不管优化变量初始值如何,神经网络都将使变量收敛到全局唯一最优解。
申请公布号 CN102662916A 申请公布日期 2012.09.12
申请号 CN201210087506.1 申请日期 2012.03.28
申请人 高俊文;刘建成 发明人 高俊文;刘建成
分类号 G06F17/10(2006.01)I;G06N3/02(2006.01)I 主分类号 G06F17/10(2006.01)I
代理机构 广州市华学知识产权代理有限公司 44245 代理人 杨晓松
主权项 1.一种基于Lagrange函数的最小二乘准则多目标优化方法,其特征在于,包括下述步骤:S1、在凸目标函数和凸约束的情况下,将多目标规划转化为单目标规划;S2、假定|x<sub>i</sub>|为观测序列,共有J个数学模型拟合描述,拟合误差e<sub>i</sub>为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>e</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow></munder><msub><mi>w</mi><mi>j</mi></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><msub><mi>e</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中<img file="FDA0000147888060000012.GIF" wi="139" he="72" />为第j个模型对x<sub>i</sub>的拟合值,e<sub>i</sub>(j)为拟合误差,i=1,2,…,N;S3、取e<sub>i</sub>为误差统计量,由矩阵理论,此时最小二乘准则下综合模型权系数约束最小的最优解由以下最优解模型获得:<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>min</mi></mtd><mtd><mi>Q</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>e</mi><mi>i</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi></mtd><mtd><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo><msub><mi>w</mi><mi>j</mi></msub><mo>></mo><mn>0</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>J</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>S4、将加权因子β介入到目标函数Q(n)中,其中观测数据n是可变的,变化之后的目标函数为:<maths num="0003"><![CDATA[<math><mrow><mi>Q</mi><mrow><mo>(</mo><mi>n</mi><mo>;</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>&beta;</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow><msup><mrow><mo>|</mo><mi>e</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>上式中采用指数加权函数β(n,i)=λ<sup>n-i</sup>(i=1,…,N),且满足0<β(n,i)≤1约束;S5、根据约束最小二乘方法估计是一个约束最佳问题,采用扩展Lagrange函数方法的解决0<β(n,i)≤1约束问题,其扩展Lagrange函数如下:<maths num="0004"><![CDATA[<math><mrow><msub><mi>L</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>,</mo><mi>l</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mi>&lambda;</mi><mrow><mi>n</mi><mo>-</mo><mi>i</mi></mrow></msup><msup><mrow><mo>|</mo><mi>e</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mi>&Sigma;</mi><mi>m</mi><mn>2</mn></munderover><msub><mi>z</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>g</mi><mi>m</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>l</mi><mi>m</mi><mn>2</mn></msubsup></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>c</mi><munderover><mi>&Sigma;</mi><mi>m</mi><mn>2</mn></munderover><msup><mrow><mrow><mo>(</mo><msub><mi>g</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><msubsup><mi>l</mi><mi>m</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>式中l是辅助变量,z是Lagrange乘子,常数c是补偿参数,λ是遗忘因子;<maths num="0006"><![CDATA[<math><mrow><msub><mi>g</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>&lambda;</mi><mo>&DoubleRightArrow;</mo><msub><mi>g</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>l</mi><mn>1</mn><mn>2</mn></msubsup><mo>=</mo><mn>0</mn></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>&lambda;</mi><mo>-</mo><mn>1</mn><mo>&DoubleRightArrow;</mo><msub><mi>g</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&lambda;</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>l</mi><mn>2</mn><mn>2</mn></msubsup><mo>=</mo><mn>0</mn></mrow></math>]]></maths>其中g<sub>1</sub>和g<sub>2</sub>是约束函数,对g<sub>1</sub>求导为-1,对g<sub>2</sub>求导为1;S6、把Q(n;β)作为新的目标函数,在步骤S3的最优解模型中去掉非负权重w<sub>j</sub>≥0,j=1,2,…,J,利用Lagrange乘子法求得解析解,最优权向量W<sup>*</sup>=(P<sup>T</sup>R<sup>-1</sup>P)<sup>-1</sup>R<sup>-1</sup>P,最小Q值Q<sup>*</sup>=(P<sup>T</sup>R<sup>-1</sup>P)<sup>-1</sup>;S7、采用最小二乘准则计算最优加权系数,并利用神经网络对步骤S3中最优解模型的优化问题进行优化计算。
地址 510507 广东省广州市天河区粤垦路198号