數(shù)值分析上機實習(xí)報告-wenkub
2023-04-07 08:39:13 本頁面
【正文】 B與CX = D是同解方程組。在新版本中也加入了對C、FORTRAN、c++、JAVA的支持，使用時可以直接調(diào)用，也可將編寫的實用程序?qū)氲絤atlab函數(shù)庫中方便以后使用時調(diào)用。在當(dāng)今世界流行的30多個數(shù)學(xué)類軟件中，MATLAB語言處于數(shù)值計算型軟件的主導(dǎo)地位，適用范圍涵蓋了工程數(shù)學(xué)的各個方面。Python 具有腳本語言中最豐富和強大的類庫，足以支持絕大多數(shù)日常應(yīng)用。而在利用計算機解決實際問題時，要根據(jù)具體情況作出可靠的理論分析，才能夠?qū)懗霰容^可靠的程序。現(xiàn)在面向數(shù)值分析問題的計算機軟件有：C、C++、MATLAB、Python、Fortran等。Fortran為“公式翻譯器”，它是世界上最早出現(xiàn)的計算機高級程序設(shè)計語言，廣泛應(yīng)用于科學(xué)和工程計算領(lǐng)域。它的有點主要有：matlab是以矩陣為基礎(chǔ)的工具，若是編一些對速度沒有要求的，進行數(shù)值計算或者信號處理的小程序，可以用matlab，且簡單。本次編程所用的軟件為MATLAB，希望通過這次作業(yè)，能夠?qū)λ辛顺醪降恼J(rèn)識，為以后的學(xué)習(xí)和工作奠定一定的基礎(chǔ)。所以我們可以用初等行變換把增廣矩陣轉(zhuǎn)換為行階梯陣，然后回代求出方程的解。3 1 4]。這里只給出運行結(jié)果：Gauss消元法GaussSeidel法理論值3) 其輸入方法同上，不再贅述。假定第k步之后的計算不再有舍入誤差，只是由δk引起的擾動δm，都有︳δm︳︳δk︳，則稱此方法是絕對穩(wěn)定的。表6 h= 表7 h=xkyky(xk)εk 表8 h=xk（*103）yk（*103）y(xk) （*103）εk（*103）（1）用標(biāo)準(zhǔn)的四階 RungeKutta 算法計算常微分方程時，不同的步長算法的穩(wěn)定性有影響。. 基本理論將積分區(qū)間[a,b]劃分n等分，步長，求積節(jié)點，在每個小區(qū)間上應(yīng)用梯形公式：然后將它們累加求和，作為所求積分I的近似值.記 將積分區(qū)間[a,b]劃分n等分，記子區(qū)間的中點為在每個小區(qū)間上應(yīng)用辛普森公式，則有其中 記 對給定的精度，若則終止計算，并取作為所求結(jié)果；否則，重復(fù)2~4步，直到滿足精度為止。運行結(jié)果：T = 2) 復(fù)化Simpson公式程序如下： fun=inline(39。運行結(jié)果：Q = 3) Romberg算法，在命令行進行調(diào)用，調(diào)用命令如下：[q,s]=Rg(39。. 基本理論首先，選擇一個接近函數(shù)零點的，計算相應(yīng)的和切線斜率（這里表示函數(shù)的導(dǎo)數(shù)）。 并且，如果不為0, 那么牛頓法將具有平方收斂的性能. 粗略的說，這意味著每迭代一次，牛頓法結(jié)果的有效數(shù)字將增加一倍。圖42 收斂情況比較表. 方法總結(jié)當(dāng)初始值為1時，計算結(jié)果不收斂。也讓我意識到了學(xué)習(xí)編程語言的重要性。本學(xué)期學(xué)習(xí)的數(shù)值分析方法，讓我在對數(shù)學(xué)的理解上又有了新的認(rèn)識。最后，感謝老師給我這樣的機會去接觸這門語言，雖然只了解了皮毛，可是仍然收獲頗多。 RB=rank(B)。 for p= 1:n1 for k=p+1:n
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1