나는 이것을 research paper에 따르려고합니다. 20 페이지의 그림 7에있는 솔루션 그래프를 복제하려고합니다. 그림 7의 스크린 샷이 있습니다. MATLAB에서 비선형 미분 방정식의 시스템 해법을 그리는 방법은 무엇입니까?
먼저 왼쪽 그림을 다시 만들고 싶습니다. 문제의 시스템은 내가 가지고있는 것입니다 dX
. 여기가 m-파일에있는 것입니다 :
function dX = CompetitionModel(t,X)
bs = 8*10^(-3);
bl = 4*10^(-3);
bh = 6.4*10^(-3);
N = bs + bl + bh;
K = 10^8;
m1 = 2*10^(-5);
m2 = 9*10^(-9);
p = 5*10^(-13);
I = 10^(-3);
T = 0;
a = 0;
dX = [X(1) * (bs * (1 - N/K) - I - T - m1) - p * X(1) * (X(2) + X(3));
X(2) * (bl * (1 - N/K) - I - a*T - m2) + m1 * X(1) + p * X(2) * (X(1) - X(3));
X(3) * (bh * (1 - N/K) - I - a*T) + m2 * X(2) + p * X(3) * (X(1) + X(2))];
end
ode45
구문이 있습니다 [T,Y] = solver(odefun,tspan,y0)
. 나는 내가 게시 한 그림에서 tspan을 얻는다. 내 초기 조건은 : S0 = 10^4; Rl0 = 0; Rh0 = 0
, 그래서 이것은 내가 가지고있는 것입니다 y0
. 나는 명령 창에 다음을 입력 :
>>[t,X1] = ode45('CompetitionModel', [0,45000], [10^4, 0, 0]);
>>[t,X2] = ode45('CompetitionModel', [0,45000], [10^4, 0, 0]);
>>[t,X3] = ode45('CompetitionModel', [0,45000], [10^4, 0, 0]);
MATLAB 지난 30 분 동안 바빴다 내 노트북이 정말 섹시지기 시작한다. 따라서 작업이 끝날 때까지 플로팅을 할 수 없으며 코드에 오류가 있는지 알 수 없습니다. 시스템의 솔루션을 얻을 수있는 더 좋은 방법이 있는지 궁금합니다. dX
.
실수를 시정 해 주셔서 감사합니다. DE의 시스템 솔루션의 플롯으로 "위상 초상화"를 이해합니다. 너는 내 혼란을 해결했다. 시간 내 주셔서 감사합니다. –