資源描述:
《——數(shù)值分析上機(jī)題》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、課程名稱:數(shù)值分析上機(jī)實(shí)習(xí)報(bào)告姓名:學(xué)號(hào):專業(yè):聯(lián)系電話:第II頁(yè)數(shù)值分析上機(jī)題目目錄序言3第1章必做題41.1必做題第一題41.1.1題目41.1.2分析41.1.3計(jì)算結(jié)果41.1.3總結(jié)61.2必做題第二題61.2.1題目61.2.2分析61.2.3計(jì)算結(jié)果61.2.4結(jié)論81.1必做題第一題81.1.1題目8第2章選做題82.1選做題第一題82.1.1題目82.1.2分析82.1.3計(jì)算結(jié)果8附錄10附錄一:必做題第一題程序10附錄二:必做題第二題程序11附錄三:選做題第一題的程序13第15頁(yè)數(shù)值分析
2、上機(jī)作業(yè)序言本次數(shù)值分析上機(jī)實(shí)習(xí)采用數(shù)學(xué)軟件。是一種用于算法開(kāi)發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)技術(shù)計(jì)算語(yǔ)言和交互式環(huán)境。在數(shù)值分析應(yīng)用中可以直接調(diào)用軟件中已有的函數(shù),同時(shí)用戶也可以將自己編寫(xiě)的實(shí)用程序?qū)氲胶瘮?shù)庫(kù)中方便自己調(diào)用?;跀?shù)學(xué)軟件的各種實(shí)用性功能與優(yōu)點(diǎn),本次數(shù)值分析實(shí)習(xí)決定采用其作為分析計(jì)算工具。1.編程效率高M(jìn)ATLAB是一種面向科學(xué)與工程計(jì)算的高級(jí)語(yǔ)言,允許使用數(shù)學(xué)形式的語(yǔ)言編寫(xiě)程序,且比BASIC、FORTRAN和C等語(yǔ)言更加接近我們書(shū)寫(xiě)計(jì)算公式的思維方式,用MATLAB編寫(xiě)程序猶如
3、在演算紙上排列出公式與求解問(wèn)題。因此,MATLAB語(yǔ)言也可通俗地稱為演算紙式科學(xué)算法語(yǔ)言。由于它編寫(xiě)簡(jiǎn)單,所以編程效率高,易學(xué)易懂2.用戶使用方便MATLAB語(yǔ)言與其他語(yǔ)言相比,較好的解決了上述問(wèn)題,把編輯、編譯、鏈接和執(zhí)行融為一體。它能在同一畫(huà)面上進(jìn)行靈活操作,快速排除輸入程序中的書(shū)寫(xiě)錯(cuò)誤、語(yǔ)法錯(cuò)誤以至語(yǔ)義錯(cuò)誤,從而加快了用戶編寫(xiě)、修改和調(diào)試程序的速度,可以說(shuō)在編程和調(diào)試過(guò)程中它是一種比VB還要簡(jiǎn)單的語(yǔ)言。3.方便的繪圖功能MATLAB的繪圖是十分方便的,它有一系列繪圖函數(shù)(命令),例如線性坐標(biāo)、對(duì)數(shù)坐標(biāo)、
4、半對(duì)數(shù)坐標(biāo)及極坐標(biāo),均只需調(diào)用不同的繪圖函數(shù)(命令),在圖上標(biāo)出圖題、XY軸標(biāo)注,格(柵)繪制也只需調(diào)用相應(yīng)的命令,簡(jiǎn)單易行。另外,在調(diào)用繪圖函數(shù)時(shí)調(diào)整自變量可繪出不變顏色的點(diǎn)、線、復(fù)線或多重線。這種為科學(xué)研究著想的設(shè)計(jì)是通用的編程語(yǔ)言所不能及的。第15頁(yè)數(shù)值分析上機(jī)作業(yè)第1章必做題1.1必做題第一題1.1.1題目某過(guò)程涉及兩變量x和y,擬分別用插值多項(xiàng)式和多項(xiàng)式擬合給出其對(duì)應(yīng)規(guī)律的近似多項(xiàng)式,已知xi與yi之間的對(duì)應(yīng)數(shù)據(jù)如下,xi=1,2,…,10yi=34.658840.371914.6448-14.27
5、21-13.357024.823475.2795103.574397.484778.2392(1)請(qǐng)用次數(shù)分別為3,4,5,6的多項(xiàng)式擬合并給出最好近似結(jié)果f(x)。(2)請(qǐng)用插值多項(xiàng)式給出最好近似結(jié)果下列數(shù)據(jù)為另外的對(duì)照記錄,它們可以作為近似函數(shù)的評(píng)價(jià)參考數(shù)據(jù)。xi=Columns1through71.50001.90002.30002.70003.10003.50003.9000Columns8through144.30004.70005.10005.50005.90006.30006.7000Colum
6、ns15through177.10007.50007.9000yi=Columns1through742.149841.462035.118224.385211.2732-1.7813-12.3006Columns8through14-18.1566-17.9069-11.02262.028419.854940.362661.0840Columns15through1779.568893.7700102.36771.1.2分析(1)假定擬合函數(shù)的形式分別為:;利用matlab編程。1.1.3計(jì)算結(jié)果第15頁(yè)數(shù)
7、值分析上機(jī)作業(yè)擬合公式和擬合圖如下三次多項(xiàng)式擬合結(jié)果:四次多項(xiàng)式擬合結(jié)果:五次多項(xiàng)式擬合結(jié)果:六次多項(xiàng)式擬合結(jié)果:圖1三次多項(xiàng)式擬合圖2四次多項(xiàng)式擬合圖3五次多項(xiàng)式擬合圖4六次多項(xiàng)式擬合第15頁(yè)數(shù)值分析上機(jī)作業(yè)1.1.4總結(jié)不管是從書(shū)本中的理論還是實(shí)際的擬合分析中,我們知道在多項(xiàng)式擬合中,擬合的次數(shù)越高,則擬合函數(shù)越逼近原函數(shù),精確度越高,此外,多項(xiàng)式擬合不會(huì)出現(xiàn)龍格現(xiàn)象,數(shù)值較為穩(wěn)定,故該題中的最佳擬合為六次擬合。1.2必做題第二題1.2.1題目用雅格比法與高斯-賽德?tīng)柕ń庀铝蟹匠探MAx=b1或Ax=b
8、2,研究其收斂性。上機(jī)驗(yàn)證理論分析是否正確,比較它們的收斂速度,觀察右端項(xiàng)對(duì)迭代收斂有無(wú)影響。(1)A行分別為A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4];b1=[-3,2,4]T;b2=[100,-200,345]T。(2)A行分別為A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1];b1=[3,2,1]T;b2=[5,0,-10]T