亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        一種求解二次模型信賴域子問題的Admas4隱式算法

        2016-11-14 10:31:56王英慧王希云
        太原科技大學學報 2016年5期
        關(guān)鍵詞:李亮測試函數(shù)折線

        王英慧,王希云

        (太原科技大學 應用科學學院,太原 030024)

        ?

        一種求解二次模型信賴域子問題的Admas4隱式算法

        王英慧,王希云

        (太原科技大學 應用科學學院,太原 030024)

        基于信賴域子問題最優(yōu)曲線的微分方程模型,在Hessian矩陣正定及步長固定的前提下,采用求解微分方程的Admas4隱式公式構(gòu)造了一條折線,稱Admas4隱式折線,用其代替最優(yōu)曲線,提出求解子問題的新算法—Admas4隱式算法。數(shù)值結(jié)果表明Admas4隱式算法比R-K4算法效果好。

        微分方程模型;信賴域子問題;Adams4隱式算法

        在非線性優(yōu)化中,信賴域方法不僅具有很好的可靠性和強適性,而且具有較強的收斂性。因而,自其出現(xiàn)起就受到高度重視,成為熱點。信賴域方法實現(xiàn)的關(guān)鍵是每步迭代時要求解一個二次模型信賴域子問題,其形式如下:

        (1)

        (1)中的各個參數(shù)表示的含義分別是:g∈Rn:目標函數(shù)在當前點的梯度,B∈Rn×n:目標函數(shù)在當前點處的Hessian矩陣或者它的近似矩陣,△∈R:信賴域半徑,δ∈Rn:要求的變量。當△發(fā)生變化時,子問題(1)的解δ*在空間形成一條曲線——最優(yōu)曲線。

        針對最優(yōu)曲線的參數(shù)方程,在Hessian矩陣正定的前提下,有如下微分方程模型[1]:

        (2)

        對于子問題(1)的求解,目前提出的方法主要有單折線法、雙折線法、切線單折線法。后來李亮提出最優(yōu)曲線的微分方程模型之后,基于此微分方程模型,先后提出了解決信賴域子問題的隱式分段折線法、分段切線法、平均歐拉算法及休恩算法[1-4]。后來又提出了解二次模型子問題的R-K4算法。該算法主要是利用解微分方程模型的R-K4公式,構(gòu)造一條R-K4折線,進而用此折線代替最優(yōu)曲線來求解子問題(1)的解。為進一步提高數(shù)值計算結(jié)果的精度,本文結(jié)合求解微分方程的較高階的Admas4隱式公式[5],構(gòu)造了Admas4折線代替最優(yōu)曲線來求解子問題(1),并將數(shù)值結(jié)果與R-K4方法做了比較。

        1 Adams4隱式折線的構(gòu)造

        Admas4隱式公式如下:

        (n=2,3,4,…)

        (3)

        其中:

        μn=μ0+nh.具體構(gòu)造如下:

        第一步:從初始點P0(μ0,δ0)開始,其中μ0=0,δ0=-B-1g.選取步長,用公式:

        (4)

        第二步:由常用的四階龍格-庫塔公式:

        若記:f01=-(B+μ0I)-1δ0,

        則:

        (5)

        計算出δ1,得出下一個節(jié)點P1(μ1,δ1),其中μ1=μ0+h.同理計算出δ2,得到節(jié)點P2(μ2,δ2).

        (6)

        其中μ0=0,μn=μn-1+h,n=1,2,3,…,N.

        2 對步長h的限定

        為了使構(gòu)造的折線滿足引理6.4.1[6],則步長需要滿足下式:

        (7)

        3 Admas4隱式折線路徑的性質(zhì)定理

        定理1 設B對稱正定,且當n=0,1時有:

        gT[fn1+2fn2+2fn3+fn4]+

        (8)

        n=2,3,…時,有:

        (9)

        則δ(τ)滿足:

        (1)‖δ(τ)‖2關(guān)于τ為單調(diào)減函數(shù)。

        (2)q[δ(τ)]關(guān)于τ為單調(diào)增函數(shù)。

        證明:(1)當τ∈[μ0,μ1],即τ∈[0,h0]時:

        則:

        由式(7)得:

        因此,‖δ(τ)‖2在區(qū)間[μ0,μ1],[μ1,μ2]上為單調(diào)減函數(shù)。

        則:

        由式(7)得:

        因此‖δ(τ)‖2在區(qū)間[μi,μi+1],i=1,2,3,…,N-1上為單調(diào)減函數(shù)。

        (2)當τ∈[μ0,μ1],即τ∈[0,h0]時:

        則:

        由已知條件式(8):

        gT[f01+2f02+2f03+f04]+

        得:(q[δ(τ)])′≥0,τ∈[μ0,μ1].

        同上,可以證明(q[δ(τ)])′≥0,τ∈[μ1,μ2].

        所以q[δ(τ)]在區(qū)間[μ0,μ1],[μ1,μ2]上關(guān)于τ為單調(diào)增函數(shù)。

        對?τ∈[μi,μi+1],即:

        (τ-μi)∈(0,hi),i=1,2,3,…,N-1時:

        由已知條件(9):

        可得:(q[δ(τ)])′≥0,τ∈[μi,μi+1].

        所以:q[δ(τ)]在區(qū)間[μi,μi+1],i=0,1,2,3,…,N-1上關(guān)于τ為單調(diào)增函數(shù)。(證畢)。

        4 算法描述

        通過上面的討論,下面給出Admas4隱式折線算法的具體步驟:

        步0 給定梯度g,正定矩陣B,信賴域半徑△.令n∶=0.

        步1 令δ0=-B-1g.

        步2 如果△≥‖δ0‖2,則取δ*=δ0,停止計算。否則,令n∶=n+1,轉(zhuǎn)步3.

        步3 選取適當步長,令:

        其中μ0=0.令:

        如果△≥‖δ1‖2,則取:

        步4 令:

        如果△≥‖δ2‖2,則取:

        步5 根據(jù)公式:

        根據(jù)公式:

        計算δn+1,轉(zhuǎn)步6.

        步6 如果△≥‖δn+1‖2,則取:

        停止計算,否則令n∶=n+1,轉(zhuǎn)步5.

        5 Admas4隱式折線算法的適定性

        由定理1可知下列結(jié)論成立:

        定理2 設B對稱正定,在Admas4隱式算法中,對任意給定的信賴域半徑△<‖B-1g‖2,則存在自然數(shù)N使得‖δN‖2≤△.

        定理1和定理2說明對任意給定的信賴域半徑△,在Admas4隱式折線δ(τ)上的近似解存在且唯一。并且用Admas4隱式算法求信賴域子問題(2)的最優(yōu)解時,最優(yōu)解δ*在信賴域邊界上取得。

        6 數(shù)值實驗

        對于附錄中給定的測試函數(shù)1和測試函數(shù)2的信賴域子問題(1),取h=0.1,選取不同的信賴域半徑△,然后將Admas4隱式算法利用MATLAB進行了數(shù)值實驗。并且用該算法求得的測試函數(shù)在近似最優(yōu)解的函數(shù)值與R-K4方法進行比較,數(shù)值結(jié)果分別列在表1和表2中。

        表1 測試函數(shù)1的數(shù)值結(jié)果

        從表1和表2的數(shù)值結(jié)果可以看出,本論文所提出的Admas4隱式算法是有效且可行的。對于測試函數(shù)1和測試函數(shù)2,當信賴域半徑△<‖B-1g‖2時,除了信賴域半徑在‖B-1g‖2附近的情況,Admas4隱式算法求得的測試函數(shù)在近似最優(yōu)解的函數(shù)值比R-K4效果要好;當信賴域半徑△≥‖B-1g‖2時,則Admas4隱式算法與R-K4算法所求得的測試函數(shù)在近似最優(yōu)解的函數(shù)值相同。因此,本論文所構(gòu)造的Admas4隱式算法比R-K4算法更好地近似了最優(yōu)曲線,而且Admas4隱式算法也是一個很好的求解信賴域子問題(1)的折線法。

        對于測試函數(shù)1,‖B-1g‖2=15.34.對于測試函數(shù)2,‖B-1g‖2=10.01.

        附錄:測試函數(shù)

        Function1:

        s.t.‖δ‖2≤△.

        Function2:

        s.t.‖δ‖2≤△

        表2 測試函數(shù)2的數(shù)值結(jié)果

        [1] 王希云,李亮,于海波.解決信賴域子問題的隱式分段折線算法[J].應用數(shù)學和力學,2014,35(6):610-619.

        [2] 王希云,李亮,張雅琦,等.一種求解二次函數(shù)模型信賴域子問題的分段切線法[J].應用數(shù)學,2015,28(1):26-32.

        [3] 朱帥,李亮,王希云.一種求解二次模型信賴域子問題的新算法[J].西南民族大學學報,2014,40(1):91-96.

        [4] 李亮,王希云,張雅琦,等.一種求解二次模型信賴域子問題的休恩算法[J].太原科技大學學報,2014,35(2):151-155.

        [5] 林成森.數(shù)值計算方法[M].北京:科學出版社,2005.

        [6] 李董輝,童曉嬌,萬中.數(shù)值最優(yōu)化算法與理論[M].北京:科學出版社,2010.

        An implied Adams4′s Algorithm for Solving Trust-region Subproblems with Quadratic Model

        WANG Ying-hui,WANG Xi-yun

        (School of Applied Science,Taiyuan University of Science and Technology,Taiyuan 030024,China)

        In the premise of Hessian matrix as a positive definite matrix and fixed step,a broken line is constructed by implied Adams4′s method according to the differential equation model,which is defined as Admas4 implicit broken line instead of the optimal curve.Through the comparison with the R-K4 method,the results indicate that Adams4 implicit algorithm has obvious advantage over the R-K4 method.

        differential equation model,trust-region subproblems,implied Adams4′s method

        2015-09-23

        山西省高?!?31’項目基金

        王英慧(1987-)女,碩士研究生,主要研究方向為最優(yōu)化理論與應用。

        1673-2057(2016)05-0406-06

        O221

        A

        10.3969/j.issn.1673-2057.2016.04.014

        猜你喜歡
        李亮測試函數(shù)折線
        折線統(tǒng)計圖
        改天請你喝酒
        故事會(2022年3期)2022-02-10 21:13:35
        白磷燃燒實驗的改進
        折線的舞臺——談含絕對值的一次函數(shù)的圖象
        從唐詩看“一代有一代之文學”理論之失
        折線
        具有收縮因子的自適應鴿群算法用于函數(shù)優(yōu)化問題
        帶勢函數(shù)的雙調(diào)和不等式組的整體解的不存在性
        必修1復習測試題一
        約束二進制二次規(guī)劃測試函數(shù)的一個構(gòu)造方法
        婷婷开心五月亚洲综合| 久久精品无码一区二区三区蜜费 | 国产乱人视频在线播放| 欧美一级色图| 国产一区二区三区av香蕉| 国产流白浆视频在线观看| 国产成人无码精品久久二区三区 | 亚洲综合精品一区二区三区| 青青草视频在线观看色| 亚洲精品www久久久| 亚洲天堂手机在线| 极品少妇被后入内射视| av影片在线免费观看| 女邻居的大乳中文字幕| 一区二区无码中出| 国产在线观看免费不卡视频| 国产情侣一区二区| 国产99久久久久久免费看| 久久亚洲第一视频黄色| 久久精品亚洲一区二区三区画质| 免费网站内射红桃视频| 小12萝8禁在线喷水观看| 日韩av在线不卡一区二区三区| 蜜桃视频一区二区三区四| 男女高潮免费观看无遮挡 | 中文字幕亚洲精品专区| 国产xxx69麻豆国语对白| 亚洲香蕉视频| 人妻熟女中文字幕在线视频| 精品一区三区视频在线观看| 欧美国产一区二区三区激情无套| 免费一级国产大片| 国产精品亚洲一区二区三区在线| 国产成人a在线观看视频免费 | 色一情一乱一伦| 亚洲免费观看网站| 一级a免费高清免在线| 无码喷潮a片无码高潮| 熟妇人妻中文av无码| 国产高潮精品一区二区三区av| 国产一区国产二区亚洲精品|