本文格式為Word版，下載可任意編輯——數(shù)值計算基礎(chǔ)試驗指導(dǎo)部分試驗源代碼復(fù)習(xí)指導(dǎo)三套試題及其答

數(shù)值計算基礎(chǔ)

試驗指導(dǎo)書

2023年

目錄

試驗一直接法解線性方程組的3試驗二插值方法12試驗三數(shù)值積分6試驗四常微分方程的數(shù)值解8試驗五迭代法解線性方程組與非線性方程10

試驗一直接法解線性方程組

一、試驗?zāi)康?/p>

把握全選主元消去法與高斯-塞德爾法解線性方程組。

二、試驗內(nèi)容

分別寫出Guass列選主元消去法與追趕法的算法，編寫程序上機調(diào)試出結(jié)果，要求所編程序適用于任何一解線性方程組問題，即能解決這一類問題，而不是某一個問題。試驗中以以下數(shù)據(jù)驗證程序的正確性。

1、用Guass列選主元消去法求解方程組

?2.5???5.3??8.1

2.3?5.1??x1???3.7??

??????9.61.53.8?2???????x??

?1.7?4.3?x3?????????5.5??

?

?

2、用追趕法求解方程組

??20000??x1????10??

?1?2000??x???0?????2?????01?200??x3???0??????????001?20??x4???0??0001?2x50

三、試驗儀器設(shè)備與材料

主流微型計算機

四、試驗原理

1、Guass列選主元消去法對于AX=B

~

AB是上三角矩陣。即：~~

1）、消元過程：將（A|B）進行變換為(|)，其中A

a?a1nb?ab?a???1a??11121121n1??????

?a21a22?a2nb2??01?a2nb2??

????????????????????????

??a????????n1an2?annbn????00?annbn

k從1到n-1a、列選主元

k?i?n

選取第k列中絕對值最大元素maxaik作為主元。b、換行

akj?aij,j?k?1,?,nbk?bi

c、歸一化

akj/akk?akj,j?k?1,?,n

bk/akk?bk

d、消元

aij?aikakj?aij,i?k?1,?,n;j?k?1,?,n

bi?aikbk?bi,i?k?1,?,n

2）、回代過程：由(~A|B~)解出x

n,xn?1,?,x1。

bn/ann?xn

n

bk???akjxj?xk,k?n?1,?,2,1

j?k?12、追趕法

線性方程組為：

??a1c1??x1????f1???a????f????b22c2??x2???2??b3a3c3??????x3???f??????????????????3??????????????

?????????????????

?b??????

n?1an?1cn?1?????

?fn?1???b???xn?1

na??????

nxn???

fn????

做LU分解為：

?

??

????1??1?1????2?2???1??L??

?3?3??2??,R????????????????????????????nn???

分解公式：

????????????

1???n?1??1???

?

?

???

??????a(i?2,3,?,n)i

?i

?bi??i?(i?2,3,?,n)??1?b1,?ii

?1

??ci???i??i?

則

(i?1,2,?,n?1)

?Ly?fAx?f?LUx?f???

?Ux?y

回代公式：

??f1y1???1????

??y?fi??iyi?1i??i?

?xn?yn??

ixi?xi?yi??

?1

(i?2,3,?,n)

(i?n?1,n?2,?,1)

五、試驗步驟

1、理解并把握全選主元消去法與高斯-塞德爾迭代法公式；2、畫出全選主元消去法與高斯-塞德爾迭代法的流程圖3、使用C語言編寫出相應(yīng)的程序并調(diào)試驗證通過

六、試驗報告要求

1、統(tǒng)一使用《武漢科技大學(xué)試驗報告》本書寫，試驗報告的內(nèi)容要求有：試驗?zāi)康摹⒃囼瀮?nèi)容、程序流程圖、源程序、運行結(jié)果及試驗小結(jié)六個部分。

2、源程序需打印后粘貼在試驗