資源描述：
《阻滯增長模型.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、實驗三一：阻止增長模型（1）設(shè)為x的線性函數(shù)，即。為人口容量，當(dāng)時，代入（1）得：解得（2）1.模型的參數(shù)估計：為了估計（2）中的參數(shù)r和，將其表達(dá)為：算出s=0.0008，r=0.2650且。由公式（2）得到圖一：圖一注：從圖一我們可以看出第8和9個點擬合的不太好，所以我們試著去掉這兩個異常點。若去掉兩個異常點后我們重新算的r=0.2875,s=0.0007.若只去掉一個點（如第8個）得到r=0.2620,s=0.0007.因此，我們只討論去掉一個點。經(jīng)過比較得出第8個點算出的值誤差更小。r=0.2620,s=0.0007且再由公式（2）得到如下圖：圖二附錄：（1）
2、x=[31.4,38.6,50.2,62.9,76.0,92.0,106.5,123.2,131.7,150.7,179.3,204.0,226.5,251.4];y=diff(x);t=x(1:13);y1=y./tP=polyfit(t,y1,1)（2）圖一：x=(0:1:13);y=[31.438.650.262.976.092.0106.5123.2131.7150.7179.3204.0226.5251.4];y1=331.25./(1+9.5494.*exp(-0.2650.*x))plot(x,y,'r*',x,y1)title('阻滯增長模型擬合（18
3、60-1990）')（3）去掉異常點得到的r和sx2=[31.4,38.6,50.2,62.9,76.0,92.0,106.5,123.2,150.7,179.3,204.0,226.5,251.4];y=[7.211.612.713.11614.516.713.751928.624.722.524.9];y2=y./x2;P=polyfit(x2,y2,1)（3）圖二x=(0:1:13);y=[31.438.650.262.976.092.0106.5123.2131.7150.7179.3204.0226.5251.4];y1=331.25./(1+9.5494.
4、*exp(-0.2650.*x))plot(x,y,'r*',x,y1)title('阻滯增長模型擬合（1860-1990）')
