Gauss-Seidel 방법에 의한 연립 방정식 풀이

소개글

% 다중 회기식 구하기
% Gauss-Seidel 방법이용

% y(x) = a*exp(b*x)의 식 양변에 로그를 취하면
% ln(y) = ln(a) + (b*x)가 되어
% Y = A + (b*x)의 형태로 식을 변형 시킬 수 있다.
% Gauss-Seidel 방법사용

목차

없음

본문내용

% 다중 회기식 구하기
% Gauss-Seidel 방법이용

xi=[1 1 2 2 3 3 4 4];
yi=[1 2 1 2 1 2 1 2];
A=[1 1 1];
for n=2:length(xi)
A=[A; 1 xi(n) yi(n)];
end
b=[18 12.8 25.7 20.6 35.0 29.8 45.5 40.3]';
xj=[0 0 0]';
N=tril(A,1);
P=N-A;
fprintf(' k a b c \n')
for k=1:10
xj=N\(P*xj+b);
fprintf('%2.0f %6.6f %6.6f %2.6e \n',k,xj)
k=k+1;
end

참고 자료

없음