clear
n=80;
x(:,1)=[0.3,0.5,0.2];m=10;
p=linspace(0,0.5,m);
for j=1:m
for i=1:n
for p=1:3
xc(1)=0.3*p(j);xd(1)=0.5*p(j);xe(1)=1-0.8*p(j);
xc(i+1)=(1-0.01*p)*xc(i)+0.08*xd(i)+0*xe(i);
xd(i+1)=0.01*p*xc(i)+0.87*xd(i)+0.01*xe(i);
xe(i+1)=0*xc(i)+0.05*xd(i)+0.99*xe(i);
x(i,j)=xc(i);
y(i,j)=xd(i);
z(i,j)=xe(i);
end
end
end
figure(1)
i=1:n;
j=1:m;
mesh(p(j),i, x(i,j));
xlabel('A城市的迁徙率p% ')
ylabel('年')
zlabel('A城人口')
ylim([1,n]);
figure(2)
i=1:n;
j=1:m;
mesh(p(j),i,y(i,j))
xlabel('B城市的迁徙率p% ');ylabel('年')
zlabel('B城市人口')
ylim([1,n]);
不知道你这个问题是否已经解决, 如果还没有解决的话: