clear;clc %清空所有變量，清除命令欄 B=5; %設定B的值 D=0.2; %設定D的值 Dc=0.01; %設定Dc的值 dt=0.1; %設定dt的值 x(1)=13; %設定人口初值為13 for n=1:100;x(n+1)=x(n)+dt*(B*x(n)-D*x(n)-Dc*x(n)*x(n)); %使用for循環，計算equation of differences，得到人口數量 endt=dt*(0:100); %計算time plot(t,x) %繪制time-Population曲線 xlabel('time') %設定圖的橫坐標名稱 ylabel('Population size') %設定圖的縱坐標名稱 title('Population') %設定圖的標題xn=x(1:end-1); %x的第1個數至倒數第2個數為x(n) xn1=x(2:end); %x的第2個數至最后一個數為x(n+1) figure plot(xn,xn1) %繪制steady states圖，橫坐標為x(n)，縱坐標為x(n+1) hold on %保持當前坐標 plot([0 600],[0 600],'r') %繪制y=x直線 scatter(x(end-1),x(end),'k') %標出兩條線的交點，用黑色圓圈表示 xlabel('Population size at t(n)') %設定圖的橫坐標名稱 ylabel('Population size at t(n+1)') %設定圖的縱坐標名稱 title('Graphical analysis:steady states') %設定圖的標題

總結

以上是生活随笔為你收集整理的人口增长模型 源代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。