부산대 기계공학 수치해석 과제2(중간고사 범위)

소개글

"부산대 기계공학 수치해석 과제2(중간고사 범위)"에 대한 내용입니다.

목차

없음

본문내용

2. 교재의 그림 5.7을 참조하여 이분법에 대한 M-파일을 개발하라. 위의 문제를 풀어 개발한 함수를 검증하라. 근사 상대오차는 2%, 최대 반복 횟수는 100으로 설정하라.
(MATLAB 코드)

func.m파일로 func함수를 다음과같은코드로 내장함수로 먼저 만들었다.
function [y] = func(x)

m= 80 ;
v=36;
t=4 ;
g= 9.81;
y = sqrt(g.*m./x).*tanh(sqrt(g.*x./m).*t)-v;

이 다음 bisect.m파일을 다음과 같이 코딩하고 ,bisect를 실행해보았다.

function [root,fx,ea,iter] = bisect(x1,xu,es,maxit)
m= 80 ;
v=36;
t=4 ;
g= 9.81;

test = func(x1)*func(xu);
if test >0, error('no sign change'), end
iter=0; xr=x1; ea = 1;
while(1)
xrold = xr;
xr = (x1 +xu)/2;
iter = iter +1;
if xr~=0, ea = abs((xr -xrold)/xr)*100; end
test = func(x1)*func(xr);
if test<0
xu = xr;
elseif test>0
x1 = xr;
else

참고 자료

없음