e******g
发帖数: 33 | 1 我是外行啊，实在不会，哭死了，拜托大家救我。
有这样一个程序
function [t,S,I,R] = Program_2_1(beta,gamma,S0,I0,MaxTime)
if nargin == 0
beta=1.4247;
gamma=0.14286;
S0=1-1e-6;
I0=1e-6;
MaxTime=70;
end
S=S0; I=I0; R=1-S-I;
% The main iteration
[t, pop]=ode45(@Diff_2_1,[0 MaxTime],[S I R],[],[beta gamma]);
S=pop(:,1); I=pop(:,2); R=pop(:,3);
% Calculates the differential rates used in the integration.
function dPop=Diff_2_1(t,pop, parameter)
beta=parameter(1); gamma=parameter(2);
S=pop(1); I=pop(2); R=pop(3); | b*****n
发帖数: 685 | | D******n
发帖数: 2836 | 3 不知道你要干啥
iteration在function里面，你要改啥？
【在 e******g 的大作中提到】
: 我是外行啊，实在不会，哭死了，拜托大家救我。
: 有这样一个程序
: function [t,S,I,R] = Program_2_1(beta,gamma,S0,I0,MaxTime)
: if nargin == 0
: beta=1.4247;
: gamma=0.14286;
: S0=1-1e-6;
: I0=1e-6;
: MaxTime=70;
: end
| e******g
发帖数: 33 | 4 我要改条件和iteration的次数。
I是在不停变大的，我希望当I到达0.1以后，接下去20个t是在beta等于0.4的情况下
iteration.
救命啊，真的很急，请指教一下。
【在 D******n 的大作中提到】
: 不知道你要干啥
: iteration在function里面，你要改啥？
|
|
