利用matlab軟件求解常數(shù)e和歐拉常數(shù)γ.docx

1/5數(shù)學(xué)實驗報告利用MATLAB軟件求解常數(shù)E和歐拉常數(shù)Γ實驗?zāi)康睦肕ATLAB軟件計算常數(shù)E和Γ，并嘗試?yán)貌煌乃惴ㄓ嬎?#xff0c;比較計算精度和時間，找到較好的算法。掌握MATLAB程序求和、求極限的方法，學(xué)會尋找更優(yōu)算法。實驗內(nèi)容1、求EE可以來源于兩個數(shù)列的極限和，即ENLIM11/XX,X∞(1式)SN1/01/11/21/31/41/5(2式)，根據(jù)1式，可在MATLAB上設(shè)計如下代碼FORN115N10NE11/NN求常數(shù)E的循環(huán)語句ENDATLONG使結(jié)果顯示16位雙精度數(shù)結(jié)果N算出的EN算出的E125937424601000029271828205201156022704813829421529102718282053234788327169239322355941127182820533571102/54271814592682492612271852349603723852718268237192298132716110034086901627182804690957531427161100340870237271828169413208215303503520654926282718281798347358E的標(biāo)準(zhǔn)值約為271828182845904523536由上述結(jié)果可知，使用1式，有很大的缺陷，不僅精度連107都沒有，而且當(dāng)N109誤差開始變大。根據(jù)2式，可得如下代碼SUM0T1FORN118TNTSUMSUM1/TEND求常數(shù)E的循環(huán)語句E1SUMATLONG使結(jié)果顯示16位雙精度數(shù)結(jié)果N算出的EN算出的E121027182818011463852250000000000000011271828182619849332666666666666667122718281828286169427083333333333341327182818284467593/55271666666666666714271828182845823062718055555555556152718281828458995727182539682539681627182818284590428271827876984127017271828182845904692718281525573192182718281828459046E的標(biāo)準(zhǔn)值約為271828182845904523536如上所示，隨著N的增大，E的計算值越來越接近E的真實值但是，當(dāng)N的值大于17后，計算的精度不再提高，原因是雙精度型數(shù)只能精確到16位，所以結(jié)果只有個位以及小數(shù)點后15位(最后一位是近似取的)，而1/181561016，所以N超過18再往下計算不會更精確。。在1式代碼中，11/N和N都只能精確到小數(shù)點后16位，兩者相乘，結(jié)果精度將只能精確到8位。在2式中，,每一項都能精確到小數(shù)點后16位，而E是所有項的和，求和后仍然能夠精確到小數(shù)點后16位。所以，對于某些使用數(shù)學(xué)軟件求解的問題，如果對精度有要求，應(yīng)該盡量使用加、減運算，少用其他的運算例如乘、除、乘方、對數(shù)等，這樣可以提高運算精度。2、求Γ如此歐拉常數(shù)Γ也可以使用MATLAB求出較為精確的值。可由公式ΓLIMN→∞11/21/31/NLNN得出。FORN110S0FORI110NSS1/IENDYSLOG10N求常數(shù)Γ的循環(huán)語句ENDATLONG使結(jié)果顯示16位雙精度數(shù)結(jié)果4/5N算出的ΓN算出的Γ106263831609742086057721616490071520582207331651529705772157148989513057771558156820680577215669900188405772656640681659057721566540213950577220664893106100577215665057043Γ的權(quán)威數(shù)值約等于057721566490153286060651209。上述結(jié)果精度約有108，雖然精度還有提高空間，然而MATLAB上運行時，結(jié)果表明，增加N值，精度提高，運算時間也將大大加長(通常N超過10就需要好幾分鐘)使用級數(shù)來計算FORN19S0FORI110NSS1/ILOG11/IENDYS求常數(shù)Γ的循環(huán)語句ENDATLONG使結(jié)果顯示16位雙精度數(shù)結(jié)果N算出的ΓN算出的Γ10531072981169884605772151649020985/520572257000798361705772156149011663057671608123512580577215660395299405771656690678709057721566039536650577210664943251Γ的權(quán)威數(shù)值約等于057721566490153286060651209。上述結(jié)果精度也約有108，N再大就計算時間大大增加了，且可能會出錯。上面兩個算法都不能算出更為精確的歐拉常數(shù)的值，因為調(diào)和級數(shù)收斂較慢，因此MATLAB計算很緩慢。這也表明，這個算法還有待改進，需要更好的算法才能計算出Γ更精確的值。實驗總結(jié)本文探索了使用數(shù)學(xué)軟件求解常見常數(shù)的方法，并比較了算法的優(yōu)劣，而算法不同，計算精度、時間相差很大，這表明，數(shù)學(xué)不能完全依靠計算機，人所編寫的算法也非常重要，計算機只是人的工具，人的思維能力才是最重要的。另外，我們也需要計算機強大計算能力的幫助，學(xué)會使用數(shù)學(xué)軟件，才能更加有效地發(fā)展數(shù)學(xué)。

總結(jié)

以上是生活随笔為你收集整理的利用matlab求解常数e,利用matlab软件求解常数e和欧拉常数γ.docx的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。