매틀랩 과제 ( how to solve ode?)-(A+)
- 최초 등록일
- 2014.06.07
- 최종 저작일
- 2013.12
- 15페이지/ 한컴오피스
- 가격 3,000원
목차
2. 프로그램과 설명
3. 프로그램 실행결과
4. Discussion
본문내용
2. 프로그램과 설명
<각 함수에 대한 설명>
-----------Euler Method Mfile----------
function [t,w]=Euler_M(yprime,tvar,wi,h)
tvari=tvar(1);
tvarf=tvar(2);
t = (tvari:h:tvarf)';
n=length(t);
w=wi*ones(n,1);
for ind= 1:n-1
w(ind+1) = w(ind) + feval(yprime,t(ind))*h;
end
end
☞Euler Method 에 대한 알고리즘은 수업시간에 작성했던 소스를 활용하여 했습니다. 그리고 w에 대해 쓴 for문의 경우는 처음에 w = wi*ones(n, 1); 이 이해가 안되서 많은 시간을 할애해서 이해하려 한 결과 t가 열방향으로 정렬 되있기에 w도 열방향 정의를 위해 ones(n, 1)을 이용함을 알 수 있었습니다. 이렇게 되면서 w의 2번째 행부터는 그 전 행의 w값을 이용하여 누적 시킨 결과가 저장됨을 알 수 있었습니다.
----------Taylor 2th Method Mfile---------
function [t, w] = Taylor2M(Y1prime, tvar, wi, h)
tvari = tvar(1);
tvarf = tvar(2);
t = (tvari : h : tvarf)';
n = length(t);
w = wi * ones(n, 1);
Y2prime = inline('-1.2*exp(2*t)+0.2*cos(t)+0.4*sin(t)','t');
for ind = 1:n-1
yprimeSol = feval(Y1prime, t(ind));
y2primeSol = feval(Y2prime, t(ind));
w(ind+1) = w(ind) + h*(yprimeSol + (h/2)*y2primeSol);
end
참고 자료
없음