/
SIRPSpatialModelDriver.m
52 lines (39 loc) · 1.03 KB
/
SIRPSpatialModelDriver.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
runs=100;
tmax=4000;
tFinal=0:1:tmax;
SFinal=zeros(length(tFinal),runs);
IFinal=zeros(length(tFinal),runs);
RFinal=zeros(length(tFinal),runs);
PFinal=zeros(length(tFinal),runs);
tic
for i=1:runs
i
[Scount,Icount,Rcount,Pcount,Tcount]=SIRPSpatialModel;
j=0;
for t=tFinal
j=j+1;
tPos=(Tcount<=t);
temp=Tcount.*tPos;
[~,x]=max(temp);
SFinal(j,i)=Scount(x);
IFinal(j,i)=Icount(x);
RFinal(j,i)=Rcount(x);
PFinal(j,i)=Pcount(x);
end
end
toc
% save(['Spatialruns',num2str(runs),'S.mat'],'SFinal')
% save(['Spatialruns',num2str(runs),'I.mat'],'IFinal')
% save(['Spatialruns',num2str(runs),'P.mat'],'PFinal')
% save(['Spatialruns',num2str(runs),'Cavg.mat'],'Cavg')
% save(['Spatialruns',num2str(runs),'Pavg.mat'],'Pavg')
figure(1)
hold on
plot(tFinal,mean(SFinal'),'b','linewidth',2)
plot(tFinal,mean(IFinal'),'r','linewidth',2)
plot(tFinal,mean(RFinal'),'g','linewidth',2)
PlotFont
figure(2)
hold on
plot(tFinal,mean(PFinal'),'m','linewidth',2)
PlotFont