西南交大2019數值分析上機實習報告要求

2021-09-28 22:31:54 字數 1477 閱讀 4772

下列資料為另外的對照記錄,它們可以作為近似函式的評價參考資料。

xi =

columns 1 through 7

1.5000 1.9000 2.3000 2.7000 3.1000 3.5000 3.9000

columns 8 through 14

4.3000 4.7000 5.1000 5.5000 5.9000 6.3000 6.7000

columns 15 through 17

7.1000 7.5000 7.9000

yi =

columns 1 through 7

42.1498 41.4620 35.1182 24.3852 11.2732 -1.7813 -12.3006

columns 8 through 14

-18.1566 -17.9069 -11.0226 2.0284 19.8549 40.3626 61.0840

columns 15 through 17

79.5688 93.7700 102.3677

3.用雅格比法與高斯-賽德爾迭代法解下列方程組ax=b,研究其收斂性,上機驗證理論分析是否正確,比較它們的收斂速度,觀察右端項對迭代收斂有無影響。

(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,

(3)a行分別為a1=[1,3],a2=[-7,1];b=[4,6]t,

4. 鬆弛因子對sor法收斂速度的影響。

用sor法求解方程組ax=b,其中

要求程式中不存係數矩陣a,分別對不同的階數取w=1.1, 1.2, ...

,1.9進行迭代,記錄近似解x(k)達到||x(k)-x(k-1)||<10-6時所用的迭代次數k,觀察鬆弛因子對收斂速度的影響,並觀察當w0或w2會有什麼影響?

5. 用runge-kutta 4階演算法對初值問題y/=-20*y,y(0)=1按不同步長求解,用於觀察穩定區間的作用,推薦兩種步長h=0.1,0.2。

注:此方程的精確解為:y=e-20x

6. 實驗內容

(1) 實際驗證梯形求積公式、simpson求積公式、newton-cotes求積公式的代數精度。

(2) 針對上述三個函式和積分區間[a,b],實驗觀察梯形求積公式、simpson求積公式和newton-cotes求積公式的復化求積公式的實際計算效果。

y=exp(-x.^2).*sin(10*x)+4; a=1; b=3;

y=sin(5*x)./x.^3;a=2*pi;b=4*pi;

y=sin(5*x)./x.^3;a=2*pi;b=9.4248;

數值分析上機報告

數學與計算科學學院 實驗報告 實驗專案名稱迭代法求線性方程組 所屬課程名稱數值方法a 實驗型別驗證型 實驗日期2011.11.03 班級學號 姓名成績 附錄1 源程式 附錄2 實驗報告填寫說明 1 實驗專案名稱 要求與實驗教學大綱一致。2 實驗目的 目的要明確,要抓住重點,符合實驗教學大綱要求。3 ...

數值分析上機實習報告

學號 姓名 專業 聯絡 2011 年 11 月 序言1.程式語言 我的數值分析上機實踐採用matlab語言作為程式語言,利用matlab集數值計算 符號計算和圖形視覺化三大功能與一體的特點程式設計實現,採用函式檔案編寫程式命令流,並將計算結果儲存為excel格式,並根據計算結果生成圖形,方便查閱 分...

《數值分析》上機實習報告

數值分析程式設計作業 一.housholder變換.求方程的解 二.求特徵值 一,程式要求 1 用household變換,把a化為三對角陣b 並列印b 2 用超鬆弛法求解bx b 取鬆弛因子 1.4,x 0 0,迭代9次 3 用列主元消去法求解bx b。4 用對分法求b 即a 的小於23且最接近23...