|
Gauss-seidel법으로 구한 I1,I3,I4 행렬
Itrue %% Cramer법으로 구한 참값... 비교하기 위해 만듬
Ea %% 백분율 상대오차 행렬
실행 결과
A =
55 0 -25
0 37 -4
-25 -4 29
a =
0 0 -0.4545
0 0 -0.1081
-0.8621 -0.1379 0
Inew =
-4.1081
-6.8690
-1.0406
Itrue =
-4.1103
-6.8695
-1.0426
Ea =
0.0762
0.
|
|
|
\\n\", k);
for(i=0;i<n;i++){
sum=0.;
for(j=0;j<n;j++){
if(i!=j)
sum=sum+A[i][j]*x[j];
}
next_x[i]=(A[i][n]-sum)/A[i][i];
printf(\"%lf \\t\", next_x[i]);
x[i]=next_x[i];
}
err = 0.;
sum = 0.;
for(i=0;i<n;i++){
for(j=0;j<n;j++)
sum = sum + pow(A[i][j],2);
}
err = pow(sum,0.5);
if(err <=
|
|
|
Gauss-Seidel
clear all
clc
format long
A=[10, 2, 6; 1, 5,3; 4, 12, 20];
b=[18; 9; 36];
xi=[0; 0; 0];
N=tril(A);
P=N-A;
n=1;
fprintf(\'n x1,x2,x3 \\n\\n\')
for k=1:16
xi=(N)\\(P*xi+b);
fprintf(\'%0.0f \',n);
fprintf(\'소수점 표기 \\n\')
disp(xi);
n=n+1;
end
%jacobi
clear all
clc
format long
A=[10, 2,
|
|
|
[i]);//자릿수 조정
String num2 = String.format(\"%.6f\", x[i]);
System.out.println(\"x[\"+i+\"]의 값은 =>\"+num2);
System.out.println(\"x[\"+i+\"]의 ea 값은 =>\"+num1);
if(ea<es)//ea<es일때 for문 종료
break
}
System.out.println(\"\");
}
System.out.println(\"------------------------------
|
|
|
Gauss-Seidel method와 같다. 만일 0≤〈1이라면 그 결과는 현재와 직전에 계산된 값들의 가중 평균을 취한 것이 되어 수렴하지 않는 시스템을 수렴하도록 만들거나, 진동하는 시스템에 감쇠를 부여하여 수렴을 촉진시킬 목적으로 사용된다. 또 1〈
|
|
메뉴펼치기
