线性l1问题的一个迭代法.pdfVIP

第9卷 第 l0期 2009年5月 科 学 技 术 与 工 程 Vo1．9 No．10 May2009 1671．1819(2009)10—2551·04 ScienceTechnologyandEngineering @ 2009 Sci．Tech．Engng． 线性 问题的一个迭代法 封b健康 (河海大学理学院，南京210098) 摘 要 运用极大熵技术，提出一个求解线性z。问题的迭代法，-q传统的直接解法相比，提出的极大熵方法具有简单易行、占 用计算机资源少的优点，最后的两个算例表明此算法是可行的。 关键词 线性z。问题 迭代法 极大熵 数值解 中图法分类号 o224； 文献标志码 A 在科技工程的许多领域中需要求解像超定线 ()=sgn(arx—b)，i岳E()。 性方程组和曲线拟合之类的问题，这时需要求解线 ()=aTx—b，i=I，2，…，rn， 性 Z问题。设 r()=(r。()，r2()，…，r())。 A=(口l，口2，…，口)= 于是， b1 厂()=II —bII。=∑ l0 一bI-∑ lr() b2 ∈R ～ ． b= i = ∑max{Fi()，一／i()}。 i=1 b 所以 )的极大熵函数为 线性Z问题可以表述为 F()=1妻ln( ) )。 i=1 m )= 一6IIJ=∑ IT。一bI(1) 其中，P0为一个可调参数，由于p +∞时，(， En i=1 问题 (1)是一个不可微优化问题，目前问题 (1) P)一致收敛于 )，所以可以用P0充分大时min 已有了一些解法，都是直接解法，就是运用适当的 F(，P)的最优解来作为问题 (1)的近似解。 手法将其转化为一个线性规划问题来求解，进而得 到其精确解。这些方法缺点是所转化的等价问题 1 极大熵函数F(x，p)的性质 规模都比原来大，当问题(1)的规模较大时，解法比 较复杂，占用计算机的资源较多。利用极大熵技术 根据文献[1]和文献 [2]可得到下列结论。 将问题(1)转化为无约束可微优化问题，再利用最 引理 l 当p一++∞时，F(x，P) )，且有 优化的解法来得到问题(1)的满足一定精度的数值 ／．()F(，P)≤乏厂()+ ln2。 口 解。与直接法相比，此方法简单易于操作，占用计 p 引理2 函数F(，p)关于p严格单调递减。 口 算机的资源较少。 引理3 函数F(，p)是关于 的凸函数。 为叙述方便，使用以下记号：