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

        ?

        絕對值方程的光滑牛頓算法

        2011-12-23 00:51:48鄧永坤
        關(guān)鍵詞:線性方程組方程組牛頓

        鄧永坤, 張 萍

        (中國礦業(yè)大學(xué) 理學(xué)院,江蘇 徐州 221116)

        絕對值方程的光滑牛頓算法

        鄧永坤, 張 萍

        (中國礦業(yè)大學(xué) 理學(xué)院,江蘇 徐州 221116)

        針對絕對值方程Ax+=b的求解問題,給出了光滑牛頓法。通過引進(jìn)極大熵函數(shù)將絕對值方程進(jìn)行光滑化處理,進(jìn)而轉(zhuǎn)化為非線性光滑方程組,利用光滑牛頓算法對其進(jìn)行求解,并對算法的收斂性和收斂速度進(jìn)行了驗(yàn)證。數(shù)值實(shí)驗(yàn)結(jié)果表明該算法是有效的。

        絕對值方程;極大熵方法;光滑牛頓算法

        0 引言

        考慮絕對值方程:其中A、B∈n×n,b∈n,絕對值依分量形式而取。絕對值方程首次由Rohn在文獻(xiàn)[1]中提出,此后許多學(xué)者對其進(jìn)行了廣泛研究。文獻(xiàn)[2]給出了該方程存在唯一解或者存在非平凡解的擇一定理及幾個等價形式,文獻(xiàn)[3]證明了該方程的求解是NP-h(huán)ard問題,文獻(xiàn)[4]給出了該方程唯一可解性的條件,文獻(xiàn)[5]給出了求解該方程的廣義牛頓法,并證明了該方程與二階錐互補(bǔ)問題的等價性,進(jìn)而延伸到求解二階錐互補(bǔ)問題。

        1 光滑化處理

        引理1[4]若矩陣A的最小奇異值大于矩陣B的最大奇異值,則對任意的b∈n,絕對值方程(1)存在唯一解。

        引理1給出了方程(1)存在唯一解的條件,為保證存在唯一解,文中總是假設(shè)該條件成立,下面考慮如何求解該類絕對值方程。

        然后對絕對值方程(1)的每一行都進(jìn)行相應(yīng)的光滑化處理,得到非線性光滑方程組:

        將方程(2)簡記為

        定理1 當(dāng)p→0時,非線性光滑方程組(3)的解即為絕對值方程(1)的解。令f(x)=Ax+-b,則誤差限:

        當(dāng)p→0時,‖F(xiàn)(x)-f(x)‖→0,所以p→0時方程(3)的解即為方程(1)的解。得證。

        2 算法及收斂性分析

        算法1 (1)給定一個初始向量x0,容許誤差ε,k:=0,構(gòu)造非線性方程組F(x)。

        (2)計算xk+1,xk+1=xk-[F'(xk)]-1F(xk),k=k+1。

        (3)若‖xk+1-xk‖≤ε,則停止,得到解x*= xk;否則轉(zhuǎn)入步驟(2)。

        記非線性光滑方程組(3)中矩陣A、B及相關(guān)等式的具體表達(dá)式為

        則F'(x)可以化為

        引理2[11]若F:D?n→n在x*處F可導(dǎo),x*是非線性方程組F(x)=0的解,S0?D是x*的鄰域,又設(shè)M:S0?D→L(n)在x*處連續(xù),M(x*)非奇異,則存在閉球S=ˉS(x*,δ)?S0,使映象G(x)=x-[M(x)]-1F(x)在x∈S上有定義,且在x*處有F導(dǎo)數(shù):

        若還有ρ(I-[M(x*)]-1F'(x*))<1,則x*是迭代序列xk+1=xk-[M(xk)]-1F(xk)的吸引點(diǎn)。

        引理3[11]假定映象F:D?n→n及D0?D,若存在常數(shù)c≥0及p∈[0,1],對所有的x、y∈D0有‖F(xiàn)'(x)-F'(y)‖≤c‖x-y‖p成立,則有以下不等式成立:

        引理2和3的證明在文獻(xiàn)[11]中有詳細(xì)描述,文中從略。

        定理2 假定x*是非線性光滑方程組(3)的解,則存在閉球S=ˉS(x*,δ)?S0(δ>0),使映象G(x)=x-[F'(x)]-1F(x)對所有x∈S有定義,且算法1生成的迭代序列{xk}超線性收斂到x*,如果還假定:

        其中α>0為常數(shù),則迭代序列{xk}至少2階收斂。

        假定式(4)成立,則由引理3[11]可知:

        3 數(shù)值實(shí)驗(yàn)

        運(yùn)用MATLAB7.0進(jìn)行編程計算,參數(shù)設(shè)為:計算精度ε=1.0e-4;p為極大熵因子;n為迭代次數(shù);t為算法運(yùn)行的CPU時間,s。

        為保證方程有唯一解,利用矩陣的奇異值分解構(gòu)造矩陣

        (此矩陣A的最小奇異值大于矩陣B的最大奇異值),

        計算結(jié)果見表1。

        表1 計算結(jié)果Table 1 Results of calculation

        算例2 求解絕對值方程:Ax+B|x|=b,其中

        為保證方程有唯一解,同上例構(gòu)造矩陣

        計算結(jié)果見表2。

        表2 計算結(jié)果Table 2 Results of calculation

        通過表1和2可以看出,光滑牛頓法在求解絕對值問題(1)時,不僅得出了方程的解,而且具有較快的收斂速度,從而進(jìn)一步驗(yàn)證了算法的有效性。

        4 結(jié)束語

        [1]ROHN J.Systems of linear interval equations[J].Linear Algebra and its Application,1989,126:39-78.

        [2]ROHN J.A theorem of the alternatives for the equation Ax+= b[J].Linear and Multilinear Algebra,2004,52(6):421-426.

        [3]MANGASARIAN O L.Absolute value programming[J].Computational Optimization and Applications,2007,36(1):43-53.

        [4]ROHN J.On unique solvability of the absolute value equation[J].Optimization Letters,2009,3(4):603-606.

        [5]HU S L,HUANG Z H,ZHANG Q.A generalized Newton method for absolute value equations associated with second order cones[J].Computational Optimization and Applications,2011,235: 1490-1501.

        [6]MANGASARIAN O L,MEYER R R.Absolute value equations[J].Linear Algebra and Its Application,2006,419(2/3):359-367.

        [7]MANGASARIAN O L.A generalized Newton method for absolute value equations[J].Optimization Letters,2009,3(1):101-108.

        [8]王愛祥,王海軍.絕對值方程的區(qū)間算法[J].貴州大學(xué)學(xué)報:自然科學(xué)版,2010,27(2):7-10.

        [9]雍龍泉,孫培民,高 凱.極大熵自適應(yīng)微粒子群混合算法求解絕對值方程[J].計算機(jī)應(yīng)用研究,2011,28(7):2479-2481.

        [10]李興斯.一類不可微優(yōu)化問題的有效解法[J].中國科學(xué):A輯,1994,24(4):371-377.

        [11]李慶揚(yáng),莫孜中,祁力群.非線性方程組的數(shù)值解法[M].北京:科學(xué)出版社,1987:8-50.

        Smooth Newton method for absolute value equations

        DENG Yongkun, ZHANG Ping
        (College of Sciences,China University of Mining&Technology,Xuzhou 221116,China)

        Aimed at the solution to the absolute value equation Ax+B x =b,this paper presents a smooth Newton method.By using the maximum entropy function,absolute value equations problem could be transformed into the nonlinear smooth equations.Then the smooth Newton method is proposed for solving the absolute value equations and its convergence and convergent rate are proved.Numerical results indicate that the method is effective.

        absolute value equation;maximum entropy method;smooth Newton method

        O242.2

        A

        1671-0118(2011)06-0499-04

        2011-10-26

        鄧永坤(1988-),男,山東省青州人,碩士,研究方向:計算數(shù)學(xué),E-mail:dengyongkun@126.com。

        (編輯王 冬)

        猜你喜歡
        線性方程組方程組牛頓
        深入學(xué)習(xí)“二元一次方程組”
        求解非線性方程組的Newton迭代與Newton-Kazcmarz迭代的吸引域
        《二元一次方程組》鞏固練習(xí)
        牛頓忘食
        一類次臨界Bose-Einstein凝聚型方程組的漸近收斂行為和相位分離
        風(fēng)中的牛頓
        失信的牛頓
        勇于探索的牛頓
        線性方程組解的判別
        非自治耗散Schr?dinger-Boussinesq方程組緊致核截面的存在性
        国产精品久久中文字幕第一页 | 免费a级毛片无码无遮挡| 国产午夜在线观看视频播放| 国产三级在线观看性色av| 羞羞色院99精品全部免| 情人伊人久久综合亚洲| 激情偷乱人伦小说视频在线| www.亚洲天堂.com| 国产专区亚洲专区久久| 一本加勒比hezyo无码专区| 欧美性猛交xxxx黑人猛交| 野外三级国产在线观看| 亚洲五月天中文字幕第一页| 亚洲av成人片色在线观看| 67194熟妇在线永久免费观看| 性做久久久久久久| 日本一道本加勒比东京热| 国产日韩精品suv| 蜜桃av噜噜一区二区三区| 国产亚洲欧美另类久久久| 亚洲av午夜福利精品一区不卡| 男人的天堂av网站| 亚洲黄色电影| 亚洲av毛片成人精品| 日本男人精品一区二区| 国产精品久久久久久婷婷| 久久亚洲高清观看| 一本到亚洲av日韩av在线天堂 | 亚洲中文字幕乱码免费看| 无码熟妇人妻av在线网站 | 97精品人人妻人人| 国产在线精品一区二区不卡| 一区二区三区不卡免费av| 白白色发布免费手机在线视频观看| 久久久久久亚洲av无码蜜芽| 成年女人免费v片| 漂亮人妻被强中文字幕乱码| 激情综合色综合啪啪开心| 少妇人妻真实偷人精品视频| 欧美日一本| 日韩一区二区三区精品视频|