t를 어떻게 선언해야하는지 모르겠습니다ㅠㅠ

function [dy]=SILAR(t,y)
t=dt;
a=0.000000000184;
b=0.8;
c=0.7;
d=0.002;
L=1235/47041434;
M=1235/47041434;
dydt=zeors(4,1);
dydt(1)={(L-M)*y(1)-a*y(1)*y(3)};
dydt(2)={a*y(1)*y(3)-b*y(2)-M*y(2)};
dydt(3)={b*y(2)-c*y(3)-d*y(3)-M*y(3)};
dydt(4)={c*y(3)-M*y(4)};
end
clear all;
clc;
[t,y]=ode(@seir,[0,300],[47041334,0,100,0])
plot(t,y)
legend('S','E','I','R')
SIER을 모델링 하기위해 다음과 같은 코드를 짰는데 계속 t가 선언되지 않았다며 실행이 되지 않습니다ㅠㅠ 도와주십시오

 채택된 답변

madhan ravi
madhan ravi 2020년 6월 27일
편집: madhan ravi 2020년 6월 27일

2 개 추천

[t,y]=ode45(@seir,[0,300],[47041334,0,100,0])
plot(t,y) % :(
legend('S','E','I','R')
function dydt = seir(t,y)
a=0.000000000184;
b=0.8;
c=0.7;
d=0.002;
L=1235/47041434;
M=1235/47041434;
dydt=zeros(4,1);
dydt(1)=(L-M)*y(1)-a*y(1)*y(3);
dydt(2)= a*y(1)*y(3)-b*y(2)-M*y(2);
dydt(3)=b*y(2)-c*y(3)-d*y(3)-M*y(3);
dydt(4)=c*y(3)-M*y(4);
end

댓글 수: 1

조현 김
조현 김 2020년 6월 27일
감사합니다ㅠㅠ

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 상미분 방정식에 대해 자세히 알아보기

태그

질문:

2020년 6월 27일

댓글:

2020년 6월 27일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!