非線(xiàn)性數(shù)列變換在一些數(shù)學(xué)物理問(wèn)題中的應(yīng)用
發(fā)布時(shí)間:2021-07-28 07:04
非線(xiàn)性數(shù)列變換是一種加速收斂數(shù)列與級(jí)數(shù),或求發(fā)散級(jí)數(shù)和(summation of divergent series)的方法,該方法能有效地解決數(shù)值計(jì)算結(jié)果精度因舍入誤差積累而惡化的問(wèn)題。本文選用兩種不同的非線(xiàn)性數(shù)列變換,針對(duì)求歐拉常數(shù)γ與無(wú)窮耦合極限這兩類(lèi)實(shí)際問(wèn)題,進(jìn)行了詳細(xì)的研究和分析。歐拉常數(shù)γ的定義式是一個(gè)收斂速度極慢的數(shù)列,Sintamarian和Lu等人對(duì)其進(jìn)行了優(yōu)化修正并給出了明確的余項(xiàng)估計(jì)表達(dá)式。我們?cè)谛拚龤W拉常數(shù)數(shù)列基礎(chǔ)上,創(chuàng)新地采用Levin變換方法加速收斂修正歐拉常數(shù)數(shù)列,得到一種有效的新方法計(jì)算歐拉常數(shù)γ。非諧振子基態(tài)能量本征值的微擾解是一個(gè)迅速發(fā)散的級(jí)數(shù),我們采用Weniger變換求發(fā)散級(jí)數(shù)和。此外,我們借助計(jì)算機(jī)代數(shù)系統(tǒng)實(shí)現(xiàn)有理化的數(shù)值計(jì)算,解決了舍入誤差的問(wèn)題。隨著變換階的增加,微擾級(jí)數(shù)系數(shù)消耗的內(nèi)存迅速增加,極易導(dǎo)致內(nèi)存溢出的情況。針對(duì)這個(gè)問(wèn)題,我們?cè)赪eniger工作的基礎(chǔ)上,壓縮程序數(shù)組維數(shù)并將計(jì)算微擾級(jí)數(shù)系數(shù)從變換迭代過(guò)程的程序中分離出來(lái)從而克服了內(nèi)存的限制,得到精度極高的無(wú)窮耦合極限近似值。
【文章來(lái)源】:華東理工大學(xué)上海市 211工程院校 教育部直屬院校
【文章頁(yè)數(shù)】:69 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
圖2.1?Levin變換結(jié)果的測(cè)量精度Akm與一類(lèi)修正數(shù)列(p?=?4,6,???,16)之間的關(guān)系
第38頁(yè)?華東理工大學(xué)碩士學(xué)位論文??600?"1?*?"?}?|?v|?*?1?r???I?*?^??500?-?-??400?-?-??,X涔??■??<f?300?-?x’?-??,〇:?/?:??¥??米■-?.??岸’??100?-?-??ik??Q?,?j?j?I???l?,?|?,?|?1?l??0?50?100?150?200?250?300??k??圖2.?2待變序列{jUiG(K^i(K)i,‘",Mi3〇o}的Levin變換結(jié)果。余項(xiàng)估計(jì)6J??=?n_16,擬合函數(shù)Afc=??a?+?hfc?+?hfc2,?fc?=?10,20,.",300,其中a?=?64.337±?1.67835,?土??0.02496,?62=-0.0019?±7.81168E-5,RSS?=?221_19159,?COD?=?0.99967。??Fig.2.2?Numerical?results?of?Levin's?transformation?for?string?{/^i〇〇〇,/^i〇〇ii?"?,?Mi3〇〇}-?Remainder?estimate??〇)n?=?n—16,fitting?function?Afcm=?a?+?fcp,?intercept?a?=?64.337?±?1.67835,6!?=?2.31095?土??0.02496,?62?=-0.0019?土?7.81168E-5,RSS?=?221.19159,?COD?=?0.99967??實(shí)際上,還有另一種效率更高的方法。我們可以改變Levin變
?I?1?I?'?I?1?If??1400?-??1200?-?-?????M?-??1000?-?-?????Z??-?^?.??X??s;?800? ̄?-??J?,?.??600?-?-??400?-?,?-??/??'?^????200?-?^?-??Q?UZ?I?1?t?I?,?I?,?[?,?L??0?200?400?600?800?1000??km??圖2.3測(cè)量精度八)^與不同長(zhǎng)度待變序列{內(nèi),化+1,?,內(nèi)+?1}最高變換階/^的關(guān)系。擬合函數(shù)為??Akm=?a?+?+?b2k^,其中a?=?10S.8181?±?6.90436,心=1.69811?±?0.03217,=??-3.60658E-4?±3.02808E-5,RSS?=?951.44871,COD?=?0.99956!啊鳌北硎疚唇(jīng)??過(guò)Levin變換時(shí),修正歐拉常數(shù)數(shù)列jU2QQQ能給出52有效數(shù)字的y??Fig.2.3?Dependence?of?measured?accuracy?on?the?highest?transformation?order?km?of?various?lengths??of?strings?{?,??Fitting?function?a+b1/cm?+?ib2fc^l,a?=105.8181?土??6.90436?,?&?=?1.69811?±0.03217,?=?—3.60658E?-?4?±?3.02808E?-?5?,?RSS?=??951.44871,COD?=?0.99956.?“A”
本文編號(hào):3307470
【文章來(lái)源】:華東理工大學(xué)上海市 211工程院校 教育部直屬院校
【文章頁(yè)數(shù)】:69 頁(yè)
【學(xué)位級(jí)別】:碩士
【部分圖文】:
圖2.1?Levin變換結(jié)果的測(cè)量精度Akm與一類(lèi)修正數(shù)列(p?=?4,6,???,16)之間的關(guān)系
第38頁(yè)?華東理工大學(xué)碩士學(xué)位論文??600?"1?*?"?}?|?v|?*?1?r???I?*?^??500?-?-??400?-?-??,X涔??■??<f?300?-?x’?-??,〇:?/?:??¥??米■-?.??岸’??100?-?-??ik??Q?,?j?j?I???l?,?|?,?|?1?l??0?50?100?150?200?250?300??k??圖2.?2待變序列{jUiG(K^i(K)i,‘",Mi3〇o}的Levin變換結(jié)果。余項(xiàng)估計(jì)6J??=?n_16,擬合函數(shù)Afc=??a?+?hfc?+?hfc2,?fc?=?10,20,.",300,其中a?=?64.337±?1.67835,?土??0.02496,?62=-0.0019?±7.81168E-5,RSS?=?221_19159,?COD?=?0.99967。??Fig.2.2?Numerical?results?of?Levin's?transformation?for?string?{/^i〇〇〇,/^i〇〇ii?"?,?Mi3〇〇}-?Remainder?estimate??〇)n?=?n—16,fitting?function?Afcm=?a?+?fcp,?intercept?a?=?64.337?±?1.67835,6!?=?2.31095?土??0.02496,?62?=-0.0019?土?7.81168E-5,RSS?=?221.19159,?COD?=?0.99967??實(shí)際上,還有另一種效率更高的方法。我們可以改變Levin變
?I?1?I?'?I?1?If??1400?-??1200?-?-?????M?-??1000?-?-?????Z??-?^?.??X??s;?800? ̄?-??J?,?.??600?-?-??400?-?,?-??/??'?^????200?-?^?-??Q?UZ?I?1?t?I?,?I?,?[?,?L??0?200?400?600?800?1000??km??圖2.3測(cè)量精度八)^與不同長(zhǎng)度待變序列{內(nèi),化+1,?,內(nèi)+?1}最高變換階/^的關(guān)系。擬合函數(shù)為??Akm=?a?+?+?b2k^,其中a?=?10S.8181?±?6.90436,心=1.69811?±?0.03217,=??-3.60658E-4?±3.02808E-5,RSS?=?951.44871,COD?=?0.99956!啊鳌北硎疚唇(jīng)??過(guò)Levin變換時(shí),修正歐拉常數(shù)數(shù)列jU2QQQ能給出52有效數(shù)字的y??Fig.2.3?Dependence?of?measured?accuracy?on?the?highest?transformation?order?km?of?various?lengths??of?strings?{?,??Fitting?function?a+b1/cm?+?ib2fc^l,a?=105.8181?土??6.90436?,?&?=?1.69811?±0.03217,?=?—3.60658E?-?4?±?3.02808E?-?5?,?RSS?=??951.44871,COD?=?0.99956.?“A”
本文編號(hào):3307470
本文鏈接:http://sikaile.net/kejilunwen/wulilw/3307470.html
最近更新
教材專(zhuān)著