1 简介
基于元胞自动机模拟生态池塘养殖模型。
2 部分代码
function [area,fish]=MyChange(area,fish,range,alpha)
n=length(area);
for i =1:n
for j=1:n
%周围9个方块共同impact
if MyRound(area,i,j)<range(1)%15
area(i,j)=MyRound(area,j,j)*(1-rand*mean(fish)*alpha(1));
%(1-rand*mean(fish)*alpha(1))鱼虾蟹等影响系数
elseif MyRound(area,i,j)<range(2)%50
area(i,j)=MyRound(area,j,j)*(1-rand*mean(fish)*alpha(2));
elseif MyRound(area,i,j)<range(3)%150
area(i,j)=MyRound(area,j,j)*(1-rand*mean(fish)*alpha(3));
elseif MyRound(area,i,j)<range(4)%500
area(i,j)=MyRound(area,j,j)*(1-rand*mean(fish)*alpha(4));
else
area(i,j)=MyRound(area,j,j)*(1-rand*mean(fish)*alpha(5));
end
end
end
avg=mean(mean(area));
%成长系数设置的很迷惑
%蟹类主要以鱼虾等动物尸体为食物,所以其数目会受到鱼虾的影响
if avg<range(1)
fish=fish*(1+(rand-0.5)*0.01);
elseif avg<range(2)
fish=fish.*[1.006 1.005 1.007];
elseif avg<range(3)
fish=fish.*[1.013 1.011 1.013];
elseif avg<range(4)
fish=fish.*[1.03 1.03 1.03];
else
fish=fish.*[1.05 1.04 1.05];
end
fish=fish+fish.*unifrnd(-0.12,0.1,1,3);%有随机性,更有说服力
end
- 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.
3 仿真结果
4 参考文献
[1]谢花林, and 李秀彬. "基于元胞自动机模型的区域生态用地调控情景模拟研究——以北京市为例." 地理学核心问题与主线——中国地理学会2011年学术年会暨中国科学院新疆生态与地理研究所建所五十年庆典论文摘要集 2011.