理學(xué)院,上海200237)1 引 言歐拉常數(shù)γ是在數(shù)學(xué)及物理領(lǐng)域十分重要的常數(shù),應(yīng)用也很廣泛[1],如計(jì)"/>

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

        ?

        Levin變換應(yīng)用于歐拉常數(shù)數(shù)列及連分式數(shù)列改進(jìn)算法的加速收斂

        2021-01-12 02:17:42姚清照賀黎明
        大學(xué)數(shù)學(xué) 2021年1期
        關(guān)鍵詞:歐拉分式表達(dá)式

        姚清照, 賀黎明

        (華東理工大學(xué) >理學(xué)院,上海200237)

        1 引 言

        歐拉常數(shù)γ是在數(shù)學(xué)及物理領(lǐng)域十分重要的常數(shù),應(yīng)用也很廣泛[1],如計(jì)算黎曼ζ函數(shù)[2]、Somos二次遞歸常數(shù)[3],以及Ising鏈模型中的關(guān)聯(lián)長度[4]等等.對歐拉常數(shù)的計(jì)算十分重要,然而要得到高精度的近似值,僅僅對歐拉常數(shù)的原始定義

        進(jìn)行直接計(jì)算是比較困難的,這就需要有更高效的理論方法.有許多針對原定義數(shù)列Dn表達(dá)式進(jìn)行修正實(shí)現(xiàn)加速收斂的工作[5-7].Lu等人[7]使用了連分式修正(modification)的方法提高收斂速度,并給出明確的余項(xiàng)表達(dá)式.本文在此修正數(shù)列的基礎(chǔ)上,首次采用通用Levin變換的方法進(jìn)行二次加速收斂修正歐拉常數(shù)數(shù)列,數(shù)值計(jì)算表明由此可以得到更為高效、精確的計(jì)算結(jié)果.

        2 Lu等人的連分式修正[7]

        用來定義歐拉常數(shù)的數(shù)列Dn收斂緩慢,收斂速度和n-1相當(dāng).Lu等人在Dn基礎(chǔ)上提出一種運(yùn)用連分式修正原定義后更快收斂到歐拉常數(shù)的數(shù)列

        (1)

        其中p,q是正整數(shù).可以證明,如果取p=1時(shí)有

        3 非線性數(shù)列變換[8]

        非線性數(shù)列變換T是一種將已知部分和數(shù)列集合{[sm]},變換成一個(gè)有望更快收斂到原數(shù)列極限s的新數(shù)列集合{[s′m]}的運(yùn)算規(guī)則

        T∶{[sm]}→{[s′m]},

        其中sm表示原數(shù)列的1到m項(xiàng)之和,m是非負(fù)整數(shù).在實(shí)際計(jì)算中,只考慮對有限個(gè)元素構(gòu)成的序列

        {sl,sl+1,…,sNm}

        進(jìn)行變換.

        Weniger對常用的加速收斂非線性數(shù)列變換方法做了很好的歸納和總結(jié).在眾多算法中,尤以Levin變換具有收斂速度快,計(jì)算效率高的突出優(yōu)點(diǎn).Levin變換假設(shè)數(shù)列部分和集合{[sm]}與它們的極限s滿足關(guān)系式

        (2)

        其中ωm是余項(xiàng)估計(jì),β為偏移常數(shù).將上式移項(xiàng)化簡并用k階線性的差分算符Δk消去關(guān)于m的k-1階多項(xiàng)式,經(jīng)過整理就可以得到通用Levin變換

        (3)

        再用Sidi[9]的方法能得到簡潔且適用于計(jì)算的通用Levin變換3項(xiàng)遞推公式:

        (4)

        本文采用Weniger的建議只考慮取β=1.

        余項(xiàng)估計(jì)值對于Levin變換是至關(guān)重要的.基于不同的余項(xiàng)模型的Levin變換稱之為u變換、v變換等.對于有明確余項(xiàng)表達(dá)式的部分和序列,可采用(2)式的模型表達(dá)式,由此得到的變換式稱通用Levin變換.此余項(xiàng)表達(dá)式的近似程度越好,Levin變換的收斂速度越快,這是Levin變換的顯著優(yōu)點(diǎn).反之,Levin變換就不能很好地收斂,這同時(shí)也是它的顯著弱點(diǎn).

        4 計(jì)算結(jié)果及討論

        Dn的余項(xiàng)大致為n-1,所以歐拉常數(shù)的直接計(jì)算是非常困難的.Lu等人[7]給出了改進(jìn)的歐拉常數(shù)的系列表達(dá)式.(1)式中p越大,連分式的形式越復(fù)雜,其余項(xiàng)越小,收斂越快.如取p=3,可有

        其余項(xiàng)為n-8,顯然比Dn收斂快得多.

        表1 連分式修正前后數(shù)列收斂的數(shù)值結(jié)果比較

        表2 通用Levin變換對連分式修正前后數(shù)列加速收斂的數(shù)值結(jié)果比較

        為了更細(xì)致地研究這個(gè)問題,利用計(jì)算機(jī)代數(shù)系統(tǒng)Maple把這里的所有結(jié)果重復(fù)進(jìn)行了計(jì)算,結(jié)果列于表3右側(cè).使用Maple的所有計(jì)算結(jié)果都采用有理分式的形式,所以就完全克服了舍入誤差的影響.與Fortran采用有限字長的浮點(diǎn)計(jì)算結(jié)果形成鮮明對比的是,由Maple得到的Levin變換結(jié)果的精度呈現(xiàn)出不斷增加的趨勢.

        表3 通用Levin變換對加速收斂的Fortran和Maple計(jì)算結(jié)果比較

        5 結(jié) 論

        歐拉常數(shù)是一個(gè)具有重要應(yīng)用價(jià)值的數(shù)學(xué)常數(shù),有針對歐拉常數(shù)的專門算法研究,但是近幾十年一直沒有很好的進(jìn)展.本文在已有的歐拉常數(shù)修正算法基礎(chǔ)上再利用更具廣泛性的非線性變換算法進(jìn)行二次加速,極大的提高了計(jì)算效率.此外針對舍入誤差對非線性變換的負(fù)面影響,可以利用Maple的有理式計(jì)算功能予以克服.通過與Fortran的浮點(diǎn)運(yùn)算結(jié)果比較,克服了舍入誤差的Levin變換可以使歐拉常數(shù)的計(jì)算精度不斷增長的趨勢保持不變.

        致謝本文工作的主要結(jié)果是建立在Lu等人的對歐拉常數(shù)數(shù)列Dn連分式修正基礎(chǔ)之上,有關(guān)非線性變換的算法主要參閱了Weniger的綜述文章.在這里向這些作者一并表示感謝.

        猜你喜歡
        歐拉分式表達(dá)式
        歐拉閃電貓
        汽車觀察(2022年12期)2023-01-17 02:20:42
        歐拉魔盒
        精致背后的野性 歐拉好貓GT
        車迷(2022年1期)2022-03-29 00:50:26
        一個(gè)混合核Hilbert型積分不等式及其算子范數(shù)表達(dá)式
        表達(dá)式轉(zhuǎn)換及求值探析
        淺析C語言運(yùn)算符及表達(dá)式的教學(xué)誤區(qū)
        如何認(rèn)識分式
        1.3 分式
        歐拉的疑惑
        拆分在分式題中的應(yīng)用
        免费人成在线观看视频播放| 人妖在线一区二区三区| 免费午夜爽爽爽www视频十八禁| 日韩吃奶摸下aa片免费观看| 日韩免费一区二区三区在线| 极品美女销魂一区二区三| 国产免费观看久久黄av麻豆| 又色又爽又高潮免费视频观看| 肉体裸交丰满丰满少妇在线观看| 亚洲国产不卡av一区二区三区| av网站免费在线浏览| 内射合集对白在线| chinese国产乱在线观看 | 手机免费在线观看av网址| 性色av闺蜜一区二区三区| 亚洲成人观看| av网站一区二区三区| 国产亚洲精品久久久久5区| 亚洲欧美激情在线一区| av深夜福利在线| 射进去av一区二区三区| 亚洲国产美女精品久久久久∴| 欧美精品久久久久久久自慰| AV无码人妻一区二区三区牛牛| 国产成人国产三级国产精品| 小sao货水好多真紧h无码视频| 国产一区a| 日本视频一区二区这里只有精品 | 9999毛片免费看| 在线看片免费人成视久网不卡| 亚洲视频一区二区三区视频| 国产精品igao视频网| 白白色免费视频一区二区| 亚洲女同免费在线观看| 国产啪亚洲国产精品无码| 97超在线视频免费| 永久免费看黄在线观看| 美女mm131爽爽爽| 亚洲av无码乱码国产麻豆穿越| 国内自拍视频在线观看| 少妇激情av一区二区三区|