반응형

전체 글 254

Gauss_Jordan 소거법

function x=gauss_jordan2(A)% A: 계수행렬% I: 단위행렬% x: inv(A)의 계산값 [m1,k1]=size(A);% [m2,k2]=size(b); x=zeros(m1,k2);I=eye(m1); % 상삼각 행력 형성for i=1:m1-1 m=-A(i+1:m1,i)/A(i,i); A(i+1:m1,:)=A(i+1:m1,:)+m*A(i,:); I(i+1:m1,:)=I(i+1:m1,:)+m*I(i,:);end%A%I % 하삼각 행렬 형성for i=1:m1-1 for j=i+1:k1 m=-A(i,j)/A(j,j); A(i,:)=A(i,:)+m*A(j,:); I(i,:)=I(i,:)+m*I(j,:); end %A %Iend%A%Ifor i=1:m1 x(i,:)=I(i,:)/A(i,i); A..

Gauss

Gauss33 , Gauss44, Gauss55를 아우르는 일반적인 가우스 일반해모든 크기의 행렬에 대해 사용 가능 function x=gauss(A,c)% A: 계수행렬% c: 상수행렬% x: 연립방정식의 해 [m1,k1]=size(A); % 상삼각행렬로 변환for i=1:m1-1m=-A(i+1:m1,i)/A(i,i);A(i+1:m1,:)=A(i+1:m1,:)+m*A(i,:);c(i+1:m1,1)=c(i+1:m1,1)+m*c(i,1);endx(m1,1)=c(m1,1)./A(m1,m1); % 후진대입법을 사용하여 근 x를 구함for i=m1-1:-1:1x(i,1)=((c(i,1))-A(i,i+1:m1)*x(i+1:m1,1))./A(i,i);end