1、精品資料實驗報告課程名稱：數值分析實驗項目：解方程/解方程組專業班級：機械工程姓名：李奇學號： 201820065實驗室號：實驗組號：實驗時間：批閱時間：指導教師：曲紹波成績：可編輯修改精品資料沈陽工業大學實驗報告（適用計算機程序設計類）專業班級：機械工程學號：201820065姓名：李奇實驗名稱：解方程可編輯修改精品資料1. 實驗目的：掌握牛頓法和割線法求解方程能夠使用 MATLAB編程求解2. 實驗內容：(1) 學習掌握牛頓法，割線法原理。(2) 學習 MATLAB 語言的語法格式。(3) 編寫 MATLAB程序 ,求解方程 cos x xex0 的最小正根，取初值x0 0，當xk 1xk

2、106時，迭代終止。3. 實驗方案（程序設計說明）牛頓法是一種線性化的近似方法，其基本思想是將非線性方程轉化為線性方程來進行迭代求解，當計算結果滿足精度時，計算中止。割線法具有超線性的收斂速度，求解速度快。根據所掌握知識將牛頓法和割線法用計算機語言表達出來，并運行求解，比較兩種求解方法的優劣。4. 實驗步驟或程序（經調試后正確的源程序）見附件 A5 程序運行結果見附件 A6 出現的問題及解決方法出現的問題：在編輯器內編寫程序時出現“字符數無效”的警告。解決方法：將輸入法切換至英文模式，重新輸入。可編輯修改精品資料可編輯修改精品資料附件 A沈陽工業大學實驗報告（適用計算機程序設計類）專業班級：機

3、械工程學號：201820065姓名：李奇實驗步驟或程序 ：1、 編寫牛頓法 M 文件2 、調用 newton.m文件，輸入求解可編輯修改精品資料3 、運行結果如下4 、編寫割線法 M 文件（ 1）子程序（2 ）主程序可編輯修改精品資料5 、調用 secant.m文件，輸入方程6 、運行結果如下可編輯修改精品資料沈陽工業大學實驗報告（適用計算機程序設計類）專業班級：機械工程學號：201820065姓名：李奇可編輯修改1. 實驗目的：精品資料掌握高斯消去法，列主元高斯消去法和LU 分解法求解方程組能夠使用 MATLAB編程2. 實驗內容：(1) 學習掌握高斯消元，列主元， LU 分解原理。(2)

4、學習 MATLAB 語言的語法格式。(3) 編寫 MATLAB 程序求解下列方程組4 x5.3y5.6z3m3.4n100.165x2.1y3.2z4m8n75.722x4 y7.2z5m2.4n98.25x3y8z2.3m3n57.14.2x3y2n3.723. 實驗方案（程序設計說明）高斯消去法的基本思想是通過消元把線性方程組化為等價的上三角方程組，再進行求解。高斯消去一般由“消元過程”和“回代過程“兩部分組成。消元過程就是按確定的計算過程對方程組的增廣矩陣進行初等變換，將原方程化為與之等價的上三角方程組；回代過程就是對得到的上三角方程組求解的過程。列主元高斯消去法的算法步驟：(1)對 k

5、 =1,2, ,n-1, 選主元素，選出 |a (k) kk | ，|a (k) k+1 ， k | , ,|a (k) nk |中絕對值最大的元素 |a (k) mk | ，將第 k 行和第 m 行交換。 (2) 消元過程。 (3) 回代過程。 (4) 輸出方程組的解。LU 分解法的算法步驟： (1)首先進行 LU 分解，計算 U 的第一行和 L 的第一列。(2)對 k=2 ，3, ，n, 計算 U 的第 k 行和 L 的第 k 列。 (4) 計算 y, 求解 Ly=b 。(4)求解 x, 求解 Ux=y 。4. 實驗步驟或程序（經調試后正確的源程序）見附件 B5 程序運行結果可編輯修改見附件 B6 出現的問題及解決方法精品資料實驗名稱：解方程組可編輯修改精品資料附件 B沈陽工業大學實驗報告（適用計算機程序設計類）專業班級：機械工程學號：201820065姓名：李奇實驗步驟或程序 ：1、 編寫高斯消元法M 文件2、 調用 Gaussxiaoqu.m文件，輸入方程組可編輯修改精品資料3 、運行結果如下可編輯修改精品資料4、 編寫高斯列主元法M 文件可編輯修改精品資料5、 調用 Gauss.m文件，輸入