|
Gauss elimination에서 구한 해 ( c1=55.91566086, c2=55.74110349, c3=93.11626404, c4=100.2373485 ) 와 비교한 결과 거의 같은 값을 얻을 수 있었다.
Truss문제와 전기회로, 마지막 농도에 관련된 문제를 통해 Gauss-Seidel Method의 장단점을 생각해 보았다. 장점으로는 변
|
|
|
method converged
10.0000
4.0003
8.0008
-1.9995
반복횟수 10회만에 수렴을 하였는데 값이 정확한 값이 아닌 약간의 오차가 발생하였다. 이 역시 컴퓨터 계산상에서 나오는 round-off 오차이다.
이완법에서 =1이면 이 방법은 Gauss-Seidel method와 같다. 만일 0≤〈1
|
|
|
C++ 코드 입니다..
한학기 내용 정리하고 C++ 코드 올려놨어요.. 1.Choleski
2.Bairstow Method
3.Bisection Method
4.Crout
5.Doolittle
6.False Position Method
7.Gauss-Elimination Method
8.Gauss-Jordan Method
9.Gauss-Seidel Method
10.Inverse Matrix Method
11.LU-Decomposition Method
12
|
|
|
GAUSS ELIMINATION ... 231
CHPATER 10 LU DECOMPOSITION AND MATRIX INVERSION ... 264
CHPATER 11 SPECIAL MATRICES AND GAUSS-SEIDEL ... 285
CHPATER 12 CASE STUDIES : LINEAR ALGEBRAIC EQUATIONS ... 305
EPILOGUE : PART THREE ... 327
PART FOUR OPTIMIZATION ... 331
CHPATER 13 ONE- DIMENSIONAL UNCONS
|
|
|
Gauss-seidel Method
double a,b,c,a0,b0,c0,ea;
b=0,c=0;
while (iter<=1000)
{
iter++;
a=(sumz-sumx*b-sumy*c)/n;
ea=fabs((a-a0)/a*100);
a0=a;
b=(sumxz-sumx*a0-sumxy*c)/sumx2;
b0=b;
c=(sumyz-sumy*a0-sumxy*b0)/sumy2;
c0=c;
if (ea<0.001)
{
//cout << " "<< iter << "\t" << a &
|
|
메뉴펼치기
