題#12.docVIP

《數值分析》期中報告 袁惠茹 題目： 利用Gauss-Seidel法，以初始值[0, 0, 0]T ，求解聯立方程組 ， ， 。 並利用鬆馳法（over-relaxiation），求取最佳之w值。 相關概念： Gauss-Seidel Iteration:　 Proceed each x-value using always the most recent approximations of other variables. The rate of convergence for Gauss-Seidel is more rapid than that for Jacobi. 　　　　　　　　, 1, 2,…, n. Over-relaxing:　 The standard relationship for Gauss-Seidel iteration for the set of equations , for variable , can be written as , 1, 2,…, n, where the superscript (k+1) indicates the (k+1) iterate. An algebraically equivalent form is , 1, 2,…, n, because is both added to and subtracted from the right side. Relaxation can be applied to get , 1, 2,…, n. which is called over-relaxation (過度鬆弛法) if and called under-relaxation (不足鬆弛法) if . 解題思路： 求解聯立方程組： 變量定義： 首先定義int型變量icount作為計數器記錄迭代的次數，初始化為0， 定義二維數組A[3][3]來表示矩陣， 並定義數組b[3]來表示矩陣， 用x[3]來存放x1，x2，x3的值，初值均為0， 以xsave[3]來存放三個變量的臨時取值， 定義變量TOL表示允許的最大誤差，並賦初值1.0e-4， 定義errmax表示誤差值。 定義文件指針fPtr，指向用來存放輸出的文件。 b)運算過程： 建立while循環，每次循環使值icount值加1表示完成一次迭代，同時將x1,x2,x3的舊值分別存放到xsave[0],xsave[1],xsave[2]中，並以公式,i = 1,2,3計算x1,x2,x3的值，將計算得出的值寫入文件，並計算誤差errmax，一直進行上述循環直至誤差值小於允許的最大誤差TOL時，跳出循環並關閉文件。 c)運算結果： 代碼運行結束後可在創建的output.dat文件中查看疊代次數及每一次疊代過程中的三個未知數的值。 過度鬆馳法： 變量定義： 除了上面所定義的變量外再新增xx[3]和w。 運算過程： 建立for循環，循環變量為w，取值為0.8到1.7，每循環一次w值加0.1， 先將icount值賦為0，用來完成對於每一個w值所進行的迭代次數的統計， 在for循環中建立while循環，並以xx來代替x， 將xx的舊值存放在xsave中，並用公式，, i = 1,2,3來計算x1,x2,x3的值， 計算誤差errmax， 將當前w值，針對該w正在進行的迭代次數icount和本次迭代中的三個x值輸出， 一直進行該while循環直到誤差值errmax小於允許的最大誤差TOL時將此時的w和對應的最終icount寫入文件並跳出循環。 運算結果： 運算結束後可在生成的文件output.dat中查看對於每一個w值的迭代次數icount並找出可以令迭代次數最小的w值作為最佳w值。 代碼： 求解聯立方程組： #include stdio.h #include math.h int main() { int icount = 0; float A[3][3] = { 4.63 , -1.21 , 3.22 , -3.07 , 5.48 , 2.11 , 1.26 , 3.11 , 4.57 }; float b[3] = { 2.22 , -3.17 , 5.11 }; float x[3] = { 0.0 , 0.0 , 0.0 } , xsave[3]; float TOL = 1.0e-4 , errmax; FILE *fPtr; fPtr = fopen( output.dat , w ); fprintf( fPtr , %5s %8s %8s %8s\n, Iter. , x1 , x2 , x3 ); fp