|
프로그램 하였다. 또한 Cramer법을 이용, 참값을 계산하여 결과 값과 비교를 하니 오차가 거의 나지 않았음을 알 수 있었다. << MATLAB REPORT >>
◉ 문제 11.29 회로도
◉ 프로그램 (Gauss-Seidel 이용)
◉ 실행 결과
◉ 결과 분석
|
|
|
matlab 검증
결론 b.
% 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 l
|
|
|
function Z(step,node)
A=[0 0.422/1.506 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0.211/1.422 0 0.211/1.422 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0.211/1.422 0 0.211/1.422 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0.211/1.422 0 0.211/1.422 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0.211/1.422 0 0.211/1.422 0 0 0 0 0 0
|
|
|
for n=0
for j=1:17
if (j==1)&(n==0)
t(n+1,j)=20;
end
if (j=2)&(n==0)
t(n+1,j)=20;
end
if (j=3)&(n==0)
t(n+1,j)=20;
end
if (j=4)&(n==0)
t(n+1,j)=20;
end
if (j=5)&(n==0)
t(n+1,j)=20;
end
|
|
|
Seidel method와 같다. 만일 0≤〈1이라면 그 결과는 현재와 직전에 계산된 값들의 가중 평균을 취한 것이 되어 수렴하지 않는 시스템을 수렴하도록 만들거나, 진동하는 시스템에 감쇠를 부여하여 수렴을 촉진시킬 목적으로 사용된다. 또 1〈≤2이
|
|
메뉴펼치기
