목차
1. ode solver의 신속성
2. ode solver의 정확성
2. ode solver의 정확성
본문내용
과제 : ode solver를 모두 사용해서 ode problem을 풀어보고
걸리는 시간과 정확성을 비교하라.
Example 1.7
function [ dx ] = exam7( t,x )
ep=1;
dx(1,1)=x(2);
dx(2,1)=-x(1)+ep*(1-(x(1))^2)*x(2);
end
t_f=500;
x0=[1 0];
fprintf('ode113\n');
tic;
[t_113 x_113]=ode113(@exam7,[0 t_f],x0);
toc;
fprintf('ode15s\n');
tic;
[t_15s x_15s]=ode15s(@exam7,[0 t_f],x0);
toc;
fprintf('ode23\n');
tic;
[t_23 x_23]=ode23(@exam7,[0 t_f],x0);
toc;
fprintf('ode23s\n');
tic;
[t_23s x_23s]=ode23s(@exam7,[0 t_f],x0);
toc;
fprintf('ode23t\n');
tic;
[t_23t x_23t]=ode23t(@exam7,[0 t_f],x0);
toc;
fprintf('ode23tb\n');
tic;
[t_23tb x_23tb]=ode23tb(@exam7,[0 t_f],x0);
toc
걸리는 시간과 정확성을 비교하라.
Example 1.7
function [ dx ] = exam7( t,x )
ep=1;
dx(1,1)=x(2);
dx(2,1)=-x(1)+ep*(1-(x(1))^2)*x(2);
end
t_f=500;
x0=[1 0];
fprintf('ode113\n');
tic;
[t_113 x_113]=ode113(@exam7,[0 t_f],x0);
toc;
fprintf('ode15s\n');
tic;
[t_15s x_15s]=ode15s(@exam7,[0 t_f],x0);
toc;
fprintf('ode23\n');
tic;
[t_23 x_23]=ode23(@exam7,[0 t_f],x0);
toc;
fprintf('ode23s\n');
tic;
[t_23s x_23s]=ode23s(@exam7,[0 t_f],x0);
toc;
fprintf('ode23t\n');
tic;
[t_23t x_23t]=ode23t(@exam7,[0 t_f],x0);
toc;
fprintf('ode23tb\n');
tic;
[t_23tb x_23tb]=ode23tb(@exam7,[0 t_f],x0);
toc
소개글