clear
n=50;
X(:,1)=[0.5;0.2];%B,C城市的初始值,
for i=1:n;
A=[0.95 0.1;0.05 0.9];%我不会三维的,所以我打算先让BC迁移
X(:,i+1)=A*X(:,i);%第一次迁移完成
end
a=linspace(0.01,0.5, 20);%A到B的迁移率
x(:,1)=[0.3 0.5];%AB城市初始值
for j=1: m
for i=1: n
A=[(1-a(j)) 0.08;a(j) 0.92;];%AB城市迁移
X(:,i+1)=A*X(:,i);
x(i,j)=X(1, i);
y(i,j)=X(2,i);%完成第二次迁移
end
end
figure(1)
i=1:n;%a=[0.95 0.1;0.05 0.9]也可以直接初始值3个,我不会写3维
j=1:m;
mesh(a(j),i,x(i,j))
xlabel('城市到郊区的迁徙率'), ylabel('时间(年)')
zlabel('城市人口')
ylim([1,n]);
你好,我是有问必答小助手,非常抱歉,本次您提出的有问必答问题,目前超出我们的服务范围,暂时无法为您解答。
首次提问人员可免费体验一次有问必答服务。目前首次提问的问题服务范围为:编程语言、Java开发、python、数据库、前端开发 领域专业技术问题,为您提供问题的解决思路和指导。不提供源码代写、项目文档代写、论文代写、作业代写、安装包资源发送或安装、软件使用指导等服务。
我们后续会持续优化,扩大我们的服务范围,为您带来更好地服务。