第六章地理信息系統(tǒng)應(yīng)用模型

1、一、一、 GISGIS應(yīng)用模型的分類應(yīng)用模型的分類 1.1.根據(jù)所表達(dá)的空間對(duì)象分為：根據(jù)所表達(dá)的空間對(duì)象分為：理論模型，又稱為數(shù)學(xué)模型：是應(yīng)用數(shù)學(xué)分析方法理論模型，又稱為數(shù)學(xué)模型：是應(yīng)用數(shù)學(xué)分析方法建立的數(shù)學(xué)表達(dá)式，反映地理過程本質(zhì)的理化規(guī)律，建立的數(shù)學(xué)表達(dá)式，反映地理過程本質(zhì)的理化規(guī)律，如地表徑流模型、海洋和大氣環(huán)流模型等如地表徑流模型、海洋和大氣環(huán)流模型等; ;經(jīng)驗(yàn)?zāi)Ｐ停夯谧兞恐g的統(tǒng)計(jì)關(guān)系或啟發(fā)式關(guān)系經(jīng)驗(yàn)?zāi)Ｐ停夯谧兞恐g的統(tǒng)計(jì)關(guān)系或啟發(fā)式關(guān)系的模型，是通過數(shù)理統(tǒng)計(jì)方法和大量觀測實(shí)驗(yàn)建立的模型，是通過數(shù)理統(tǒng)計(jì)方法和大量觀測實(shí)驗(yàn)建立的模型，如水土流失模型、適宜性分析模型等的模型，如水土

2、流失模型、適宜性分析模型等; ;混合模型：模型中既有基于理論原理的確定性變量，混合模型：模型中既有基于理論原理的確定性變量，也有應(yīng)用經(jīng)驗(yàn)加以確定的不確定性變量，如資源分也有應(yīng)用經(jīng)驗(yàn)加以確定的不確定性變量，如資源分配模型、位置選擇模型等。配模型、位置選擇模型等。 按照研究對(duì)象的瞬時(shí)狀態(tài)和發(fā)展過程，按照研究對(duì)象的瞬時(shí)狀態(tài)和發(fā)展過程，可將模型分為靜態(tài)、半靜態(tài)和動(dòng)態(tài)三類?？蓪⒛Ｐ头譃殪o態(tài)、半靜態(tài)和動(dòng)態(tài)三類。靜態(tài)模型用于分析地理現(xiàn)象及要素相互靜態(tài)模型用于分析地理現(xiàn)象及要素相互作用的格局作用的格局; ;半靜態(tài)模型用于評(píng)價(jià)應(yīng)用目半靜態(tài)模型用于評(píng)價(jià)應(yīng)用目標(biāo)的變化影響；動(dòng)態(tài)模型用于預(yù)測研究標(biāo)的變化影響；動(dòng)態(tài)模型

3、用于預(yù)測研究目標(biāo)的時(shí)空動(dòng)態(tài)演變及趨勢。目標(biāo)的時(shí)空動(dòng)態(tài)演變及趨勢。 1. GIS1. GIS環(huán)境內(nèi)的模型建造。指應(yīng)用者利用環(huán)境內(nèi)的模型建造。指應(yīng)用者利用GISGIS軟軟件的宏語言件的宏語言( (如如ARC/INFOARC/INFO的的AMLAML和和System9System9的的ATPATP等等) )發(fā)展各自所需的空間分析模型。這種模發(fā)展各自所需的空間分析模型。這種模型法是將由型法是將由GISGIS軟件支持的功能看作模型部件，軟件支持的功能看作模型部件，按照分析目的和標(biāo)準(zhǔn)，對(duì)部件進(jìn)行有機(jī)的組合。按照分析目的和標(biāo)準(zhǔn)，對(duì)部件進(jìn)行有機(jī)的組合。因此，這種構(gòu)模法能充分地利用因此，這種構(gòu)模法能充分地利用G

4、ISGIS軟件本身軟件本身所具有的資源，模型建造和開發(fā)的效率比較高。所具有的資源，模型建造和開發(fā)的效率比較高。2.GlS2.GlS外部的模型建造。這種方法是基于應(yīng)用外部的模型建造。這種方法是基于應(yīng)用GISGIS的空間數(shù)據(jù)庫和輸出功能，而模型分析功的空間數(shù)據(jù)庫和輸出功能，而模型分析功能則主要是利用其他應(yīng)用領(lǐng)域的軟件。這種構(gòu)能則主要是利用其他應(yīng)用領(lǐng)域的軟件。這種構(gòu)模法雖然運(yùn)行效率受到很大影響，但實(shí)現(xiàn)了軟模法雖然運(yùn)行效率受到很大影響，但實(shí)現(xiàn)了軟件的嫁接，無需在件的嫁接，無需在GISGIS環(huán)境中重編分析軟件，環(huán)境中重編分析軟件，并具有廣泛的適用性。并具有廣泛的適用性。3.3.混合型的模型建造。這是上述

5、兩種建模法的混合型的模型建造。這是上述兩種建模法的結(jié)合，即盡可能利用結(jié)合，即盡可能利用GlSGlS提供的功能，最大限提供的功能，最大限度地減少用戶自行開發(fā)的壓力，又不失具有外度地減少用戶自行開發(fā)的壓力，又不失具有外部建模法的靈活效果，例如自然語言命令或地部建模法的靈活效果，例如自然語言命令或地圖代數(shù)法，就屬于這類方法。圖代數(shù)法，就屬于這類方法。這類方法常用的自然語言命令，例如這類方法常用的自然語言命令，例如: :Avarage Exposure Times 2 Plus Steepness Times 3 For Cost Analysis(益本分析益本分析);Subtract V Minus

6、 U For Change Analysis(動(dòng)態(tài)動(dòng)態(tài)分析分析);Renumber Landuse For Water Assigning 0 To 13, 56, And 1 To 4 (聚合分析聚合分析);Add Overlay l To Overlay 2 For Overlay3(疊加疊加分析分析);Spred Roads To X Meters(緩沖區(qū)分析緩沖區(qū)分析)。適宜性分析是指土地針對(duì)某種特定開發(fā)活適宜性分析是指土地針對(duì)某種特定開發(fā)活動(dòng)的分析，這些開發(fā)活動(dòng)包括農(nóng)業(yè)應(yīng)用、動(dòng)的分析，這些開發(fā)活動(dòng)包括農(nóng)業(yè)應(yīng)用、城市化選址、作物類型布局、道路選線、城市化選址、作物類型布局、道路選線、選

7、擇重新造林的最適宜的土地等。選擇重新造林的最適宜的土地等。 首先確定具體的開發(fā)活動(dòng)，其次選擇其影首先確定具體的開發(fā)活動(dòng)，其次選擇其影響因子，然后評(píng)判某一地域的各個(gè)因子響因子，然后評(píng)判某一地域的各個(gè)因子對(duì)這種開發(fā)活動(dòng)的適宜程度，以作為土對(duì)這種開發(fā)活動(dòng)的適宜程度，以作為土地利用規(guī)劃決策的依據(jù)。地利用規(guī)劃決策的依據(jù)。 設(shè)有某項(xiàng)評(píng)價(jià)目標(biāo)或開發(fā)活動(dòng)設(shè)有某項(xiàng)評(píng)價(jià)目標(biāo)或開發(fā)活動(dòng)T,T,該該T T對(duì)應(yīng)一組影響對(duì)應(yīng)一組影響因素因素X1,X2,XmX1,X2,Xm； 每個(gè)因素對(duì)應(yīng)一組參評(píng)因子每個(gè)因素對(duì)應(yīng)一組參評(píng)因子 每每個(gè)因子有一組屬性個(gè)因子有一組屬性v1,v2, ,vm; v1,v2, ,vm; 因此，每個(gè)因素

8、對(duì)因此，每個(gè)因素對(duì)應(yīng)一個(gè)屬性集應(yīng)一個(gè)屬性集Vi:Vi: 顯然，每個(gè)因素的屬性集都是一個(gè)對(duì)指定的顯然，每個(gè)因素的屬性集都是一個(gè)對(duì)指定的T T從從優(yōu)到劣的全序集，且滿足：優(yōu)到劣的全序集，且滿足： 各個(gè)參評(píng)因子及其屬性值的取得由數(shù)據(jù)庫提取或由各個(gè)參評(píng)因子及其屬性值的取得由數(shù)據(jù)庫提取或由GISGIS空間分析軟件生成，各個(gè)因素按其屬性集的優(yōu)劣，空間分析軟件生成，各個(gè)因素按其屬性集的優(yōu)劣，可用下列矩陣表示可用下列矩陣表示: : 式中式中: P: P為為XiXi對(duì)對(duì)T T的貢獻(xiàn)函數(shù)值的貢獻(xiàn)函數(shù)值; ;W W為為XiXi對(duì)對(duì)T T的權(quán)重值。的權(quán)重值。P P值的確定方法為值的確定方法為: : 將各因子最適宜的指

9、標(biāo)值將各因子最適宜的指標(biāo)值定為貢獻(xiàn)函數(shù)值定為貢獻(xiàn)函數(shù)值l00l00，將各因子最不適宜的指標(biāo)值，將各因子最不適宜的指標(biāo)值定為貢獻(xiàn)函數(shù)值定為貢獻(xiàn)函數(shù)值0 0，在這之間，指標(biāo)值與函數(shù)值按，在這之間，指標(biāo)值與函數(shù)值按線性關(guān)系計(jì)算和確定。線性關(guān)系計(jì)算和確定。有了上述矩陣數(shù)據(jù)和有了上述矩陣數(shù)據(jù)和GlSGlS功能的支持，可以求功能的支持，可以求取基于柵格單元的評(píng)價(jià)分值：取基于柵格單元的評(píng)價(jià)分值： 然后，根據(jù)使然后，根據(jù)使G(Pji)=1-|R(T)-Pji/100|(6-3) G(Pji)=1-|R(T)-Pji/100|(6-3) 的值為最大時(shí)的的值為最大時(shí)的P P所對(duì)應(yīng)的所對(duì)應(yīng)的j,j,即為所求的即為所

10、求的某個(gè)某個(gè)T T的適宜級(jí)的適宜級(jí)SjSj。顯然，當(dāng)適宜級(jí)為。顯然，當(dāng)適宜級(jí)為S2S2或或S3S3時(shí)，必須同時(shí)確定其限制性因子。時(shí)，必須同時(shí)確定其限制性因子。限制性因子的計(jì)算公式為限制性因子的計(jì)算公式為 基于基于GISGIS的土地評(píng)價(jià)方法和過程。的土地評(píng)價(jià)方法和過程。(1)(1)評(píng)價(jià)對(duì)象，玉米種植用地的土地適宜性評(píng)價(jià)，通評(píng)價(jià)對(duì)象，玉米種植用地的土地適宜性評(píng)價(jià)，通過評(píng)價(jià)將研究區(qū)分出不同的適宜性等級(jí)過評(píng)價(jià)將研究區(qū)分出不同的適宜性等級(jí):S1(:S1(最適宜最適宜) )、S2(S2(次適宜次適宜) )、S3(S3(臨界適宜臨界適宜) )和和N(N(不適宜不適宜) )。(2)(2)評(píng)價(jià)方法。采用基于評(píng)價(jià)

11、方法。采用基于GlSGlS的土地質(zhì)量評(píng)價(jià)法，即的土地質(zhì)量評(píng)價(jià)法，即將玉米作物生長有關(guān)主導(dǎo)生態(tài)條件與土地質(zhì)量將玉米作物生長有關(guān)主導(dǎo)生態(tài)條件與土地質(zhì)量( (供水、供水、供肥等供肥等) )相比照，從而評(píng)定土地的適宜性等級(jí)。相比照，從而評(píng)定土地的適宜性等級(jí)。3 3）評(píng)價(jià)過程。如下。）評(píng)價(jià)過程。如下。評(píng)價(jià)對(duì)象生態(tài)條件的調(diào)查。評(píng)價(jià)對(duì)象玉米屬于評(píng)價(jià)對(duì)象生態(tài)條件的調(diào)查。評(píng)價(jià)對(duì)象玉米屬于禾本科，為一年生草本，其主導(dǎo)生態(tài)條件例如禾本科，為一年生草本，其主導(dǎo)生態(tài)條件例如: :性喜性喜高溫，需水量大，要求土壤肥沃和土層疏松，其根高溫，需水量大，要求土壤肥沃和土層疏松，其根系生長要求防止土壤侵蝕等。系生長要求防止土壤侵

12、蝕等。確定評(píng)價(jià)對(duì)象的影響因素和因子。根據(jù)將玉米作物生確定評(píng)價(jià)對(duì)象的影響因素和因子。根據(jù)將玉米作物生長有關(guān)的生態(tài)條件與土地質(zhì)量相比照，除了溫度可通長有關(guān)的生態(tài)條件與土地質(zhì)量相比照，除了溫度可通過季節(jié)調(diào)節(jié)外，其他影響因素和因子如圖過季節(jié)調(diào)節(jié)外，其他影響因素和因子如圖6-16-1所示所示: : 利用利用GISGIS生成影響因素?cái)?shù)據(jù)。例如供水有效性，其影生成影響因素?cái)?shù)據(jù)。例如供水有效性，其影響因子為水源和土層厚度，它們對(duì)供水或保水有效性響因子為水源和土層厚度，它們對(duì)供水或保水有效性的屬性集如表的屬性集如表6-26-2所示所示: : 根據(jù)表根據(jù)表6-26-2和相應(yīng)的空間數(shù)據(jù)，通過如圖和相應(yīng)的空間數(shù)據(jù)，通

13、過如圖6-26-2所示的所示的GISGIS操作過程，可以有效生成基于柵格操作過程，可以有效生成基于柵格單元的供水條件等級(jí)數(shù)據(jù)文件單元的供水條件等級(jí)數(shù)據(jù)文件R1R1。同理，可。同理，可以生成供肥條件等級(jí)數(shù)據(jù)文件以生成供肥條件等級(jí)數(shù)據(jù)文件R2R2、供氧條件、供氧條件等級(jí)數(shù)據(jù)文件等級(jí)數(shù)據(jù)文件R3R3和土壤侵蝕等級(jí)數(shù)據(jù)文件和土壤侵蝕等級(jí)數(shù)據(jù)文件R4R4。這些數(shù)據(jù)文件是土地適宜性評(píng)價(jià)的重要基礎(chǔ)，這些數(shù)據(jù)文件是土地適宜性評(píng)價(jià)的重要基礎(chǔ)，而它們的可靠性又取決于各個(gè)參評(píng)因子及其而它們的可靠性又取決于各個(gè)參評(píng)因子及其屬性集的確定。屬性集的確定。計(jì)算各因素的權(quán)重和貢獻(xiàn)函數(shù)值。由于每計(jì)算各因素的權(quán)重和貢獻(xiàn)函數(shù)值。由

14、于每個(gè)影響因素或土地質(zhì)量個(gè)影響因素或土地質(zhì)量Qi(i=1,2,3,4)Qi(i=1,2,3,4)對(duì)評(píng)價(jià)對(duì)評(píng)價(jià)對(duì)象的貢獻(xiàn)大小和適宜程度不同，因此不同對(duì)象的貢獻(xiàn)大小和適宜程度不同，因此不同因素有不同的權(quán)重因素有不同的權(quán)重WiWi，不同級(jí)別的因素對(duì)評(píng)，不同級(jí)別的因素對(duì)評(píng)價(jià)對(duì)象的貢獻(xiàn)函數(shù)值價(jià)對(duì)象的貢獻(xiàn)函數(shù)值PjiPji也不相同。也不相同。 首先，將土地質(zhì)量按其對(duì)評(píng)價(jià)對(duì)象重要性的首先，將土地質(zhì)量按其對(duì)評(píng)價(jià)對(duì)象重要性的不同，進(jìn)行由小到大的排隊(duì)；其次，確定不同，進(jìn)行由小到大的排隊(duì)；其次，確定Qi+1Qi+1對(duì)對(duì)QiQi的重要性程度的重要性程度Vi+1,iVi+1,i，用倍數(shù)表，用倍數(shù)表示，并令示，并令V1,

15、0=1V1,0=1，其結(jié)果如表，其結(jié)果如表6-36-3所示所示; ;最最后，按照下列公式計(jì)算各土地質(zhì)量的權(quán)重后，按照下列公式計(jì)算各土地質(zhì)量的權(quán)重: : 式中，式中，U1=V1,0, U2=V2,1U1=V1,0, U2=V2,1V1,0, Ui=Vi,i-1V1,0, Ui=Vi,i-1Vi-1,i-Vi-1,i-2 2V1,0V1,0，所以，和土地質(zhì)量，所以，和土地質(zhì)量Qi(i=1,2,3,4)Qi(i=1,2,3,4)對(duì)應(yīng)的權(quán)重分別對(duì)應(yīng)的權(quán)重分別為為W1=0.67W1=0.67，W2=0.17W2=0.17，W3=0.08W3=0.08，W4=0.08W4=0.08。貢獻(xiàn)函數(shù)值可以按照影響

16、因素的級(jí)別來確定。例如，如果將貢獻(xiàn)函數(shù)值可以按照影響因素的級(jí)別來確定。例如，如果將各因子最適宜的指標(biāo)值定為貢獻(xiàn)函數(shù)值各因子最適宜的指標(biāo)值定為貢獻(xiàn)函數(shù)值100100，將各因子最不適宜，將各因子最不適宜的指標(biāo)值定為貢獻(xiàn)函數(shù)值的指標(biāo)值定為貢獻(xiàn)函數(shù)值0 0，當(dāng)影響因素分為，當(dāng)影響因素分為4 4級(jí)時(shí)，則各因素的級(jí)時(shí)，則各因素的P P值分別為值分別為P1i=100, P2i=67, P3i=33, P4i=0P1i=100, P2i=67, P3i=33, P4i=0根據(jù)各因素的權(quán)重與貢獻(xiàn)函數(shù)值的關(guān)系式，可以建立各因素根據(jù)各因素的權(quán)重與貢獻(xiàn)函數(shù)值的關(guān)系式，可以建立各因素的評(píng)價(jià)指標(biāo)表的評(píng)價(jià)指標(biāo)表( (表表6

17、 64)4)，該指標(biāo)表與不同土地質(zhì)量的數(shù)據(jù)文件，該指標(biāo)表與不同土地質(zhì)量的數(shù)據(jù)文件RiRi相結(jié)合，為土地適宜性評(píng)價(jià)提供依據(jù)。相結(jié)合，為土地適宜性評(píng)價(jià)提供依據(jù)。計(jì)算機(jī)適宜性評(píng)級(jí)。通過將不同土地質(zhì)量計(jì)算機(jī)適宜性評(píng)級(jí)。通過將不同土地質(zhì)量數(shù)據(jù)文件數(shù)據(jù)文件RiRi的等級(jí)，切換為與表的等級(jí)，切換為與表6-46-4相同等相同等級(jí)對(duì)應(yīng)的指標(biāo)值，便可計(jì)算出基于柵格單級(jí)對(duì)應(yīng)的指標(biāo)值，便可計(jì)算出基于柵格單元的評(píng)價(jià)分值元的評(píng)價(jià)分值R(T):R(T): 然后，根據(jù)使然后，根據(jù)使 G(Pji)=1-|R(T)-Pji/100| G(Pji)=1-|R(T)-Pji/100| 的值為最大時(shí)的的值為最大時(shí)的P P所對(duì)應(yīng)的所對(duì)應(yīng)

18、的j j，即為所求的對(duì)，即為所求的對(duì)T T的的適宜性等級(jí)適宜性等級(jí)SjSj。例如，當(dāng)某柵格單元。例如，當(dāng)某柵格單元Q1Q1的指標(biāo)值的指標(biāo)值為為67.0,Q267.0,Q2為為11.411.4，Q3Q3為為2.62.6，Q4Q4為為0 0時(shí)，該柵格時(shí)，該柵格單元的單元的R(T)=81.0R(T)=81.0，然后將該值和，然后將該值和PjiPji的值依次代的值依次代入式（入式（6 63 3）進(jìn)行計(jì)算，只有當(dāng)?。┻M(jìn)行計(jì)算，只有當(dāng)取P2i=67P2i=67時(shí)，時(shí)，G(P2i)G(P2i)的值為最大的值為最大(0.86)(0.86)，所以該柵格單元的適，所以該柵格單元的適宜性等級(jí)為宜性等級(jí)為S2(S2(

19、次適宜）。直至研究地區(qū)全部柵次適宜）。直至研究地區(qū)全部柵格單元都獲得相應(yīng)的適宜性等級(jí)，得到玉米作物格單元都獲得相應(yīng)的適宜性等級(jí)，得到玉米作物種植用地的適宜性分級(jí)圖（圖種植用地的適宜性分級(jí)圖（圖6 63 3）。）。 一、一、 一般介紹一般介紹 發(fā)展預(yù)測是運(yùn)用已有的存儲(chǔ)數(shù)據(jù)和系統(tǒng)提供的手發(fā)展預(yù)測是運(yùn)用已有的存儲(chǔ)數(shù)據(jù)和系統(tǒng)提供的手段，對(duì)事物進(jìn)行科學(xué)的數(shù)量分析，探索某一事物在今段，對(duì)事物進(jìn)行科學(xué)的數(shù)量分析，探索某一事物在今后的可能發(fā)展趨勢，并作出評(píng)價(jià)和估計(jì)，以調(diào)節(jié)、控后的可能發(fā)展趨勢，并作出評(píng)價(jià)和估計(jì)，以調(diào)節(jié)、控制計(jì)劃或行動(dòng)，在地理信息研究中，如人口預(yù)測、資制計(jì)劃或行動(dòng)，在地理信息研究中，如人口預(yù)測、

20、資源預(yù)測、糧食產(chǎn)量預(yù)測以及社會(huì)經(jīng)濟(jì)發(fā)展預(yù)測等，都源預(yù)測、糧食產(chǎn)量預(yù)測以及社會(huì)經(jīng)濟(jì)發(fā)展預(yù)測等，都是經(jīng)常要解決的問題。是經(jīng)常要解決的問題。預(yù)測方法通常分為定性、定量、定時(shí)和概率預(yù)測。在預(yù)測方法通常分為定性、定量、定時(shí)和概率預(yù)測。在信息系統(tǒng)中，一般采用定量預(yù)測方法，它利用系統(tǒng)存信息系統(tǒng)中，一般采用定量預(yù)測方法，它利用系統(tǒng)存儲(chǔ)的多目標(biāo)統(tǒng)計(jì)數(shù)據(jù)，由一個(gè)或幾個(gè)變量的值，來預(yù)儲(chǔ)的多目標(biāo)統(tǒng)計(jì)數(shù)據(jù)，由一個(gè)或幾個(gè)變量的值，來預(yù)測或控制另一個(gè)研究變量的取值。這種數(shù)量預(yù)測常用測或控制另一個(gè)研究變量的取值。這種數(shù)量預(yù)測常用的數(shù)學(xué)方法有移動(dòng)平均數(shù)法、指數(shù)平滑法、趨勢分析的數(shù)學(xué)方法有移動(dòng)平均數(shù)法、指數(shù)平滑法、趨勢分析法、時(shí)

21、間序列分析法、回歸分析法，以及灰色系統(tǒng)理法、時(shí)間序列分析法、回歸分析法，以及灰色系統(tǒng)理論模型的應(yīng)用。論模型的應(yīng)用。以下以人口和勞動(dòng)力的預(yù)測為例，說明人口統(tǒng)計(jì)數(shù)據(jù)如何以下以人口和勞動(dòng)力的預(yù)測為例，說明人口統(tǒng)計(jì)數(shù)據(jù)如何在定量預(yù)測模型中的應(yīng)用。在定量預(yù)測模型中的應(yīng)用。 根據(jù)人口預(yù)測模型：根據(jù)人口預(yù)測模型：式中：式中：Pt Pt 第第t t年人口數(shù)年人口數(shù); ;P0 P0 基年人口數(shù)基年人口數(shù); ; 人口出生率人口出生率; ; 人口死亡率人口死亡率; ;t t 時(shí)間時(shí)間( (年份年份) )。設(shè)根據(jù)研究地區(qū)一組人口統(tǒng)計(jì)數(shù)據(jù)的分析，得設(shè)根據(jù)研究地區(qū)一組人口統(tǒng)計(jì)數(shù)據(jù)的分析，得12.5%,12.5%,6.56

22、.5，將基年定為，將基年定為19851985年，并且年，并且P0=612.7P0=612.7萬人。萬人。設(shè)每年凈遷入該研究地區(qū)的人口數(shù)為設(shè)每年凈遷入該研究地區(qū)的人口數(shù)為W=5W=5萬人。則萬人。則 于是可得到規(guī)劃期的人口預(yù)測數(shù)（表于是可得到規(guī)劃期的人口預(yù)測數(shù)（表6 65 5）。）。同理，根據(jù)勞動(dòng)力預(yù)測方程同理，根據(jù)勞動(dòng)力預(yù)測方程: :L(t)=LR(t)L(t)=LR(t)* *L(t-1)+LW(t)L(t-1)+LW(t)式中式中: L(t): L(t)為第為第t t年勞動(dòng)力狀態(tài)向量，即年勞動(dòng)力狀態(tài)向量，即 ; ;LWLW（t t）為第）為第t t年勞動(dòng)力遷移向量，即年勞動(dòng)力遷移向量，即

23、LR(t)LR(t)為勞動(dòng)力存留系數(shù)矩陣，即為勞動(dòng)力存留系數(shù)矩陣，即 式中下標(biāo)式中下標(biāo)18-6018-60表示勞動(dòng)力的年齡表示勞動(dòng)力的年齡; ;r r代表分年齡層的勞動(dòng)力存留比率。代表分年齡層的勞動(dòng)力存留比率。 于是，得到研究地區(qū)規(guī)劃期勞動(dòng)力的預(yù)測數(shù)（表于是，得到研究地區(qū)規(guī)劃期勞動(dòng)力的預(yù)測數(shù)（表6-66-6）。）。有了這些預(yù)測的結(jié)果，將這些結(jié)果與表示每個(gè)鎮(zhèn)、市有了這些預(yù)測的結(jié)果，將這些結(jié)果與表示每個(gè)鎮(zhèn)、市中心點(diǎn)的中心點(diǎn)的x x、y y坐標(biāo)聯(lián)系起來，便得到一組點(diǎn)的數(shù)據(jù)，坐標(biāo)聯(lián)系起來，便得到一組點(diǎn)的數(shù)據(jù)，這組數(shù)據(jù)加上研究地區(qū)的邊界數(shù)據(jù)，輸入這組數(shù)據(jù)加上研究地區(qū)的邊界數(shù)據(jù)，輸入SYMAPSYMAP程

24、程序，通過使用繪制等值線這個(gè)權(quán)，便輸出一幅人口序，通過使用繪制等值線這個(gè)權(quán)，便輸出一幅人口發(fā)展預(yù)測圖發(fā)展預(yù)測圖( (圖圖6-4)6-4) 地學(xué)模擬模型是應(yīng)用計(jì)算機(jī)、數(shù)字模地學(xué)模擬模型是應(yīng)用計(jì)算機(jī)、數(shù)字模擬技術(shù)及綜合分析的方法來模擬許多地理擬技術(shù)及綜合分析的方法來模擬許多地理過程或現(xiàn)象，例如氣候變遷、沙漠化過程、過程或現(xiàn)象，例如氣候變遷、沙漠化過程、土地退化過程、湖泊沼澤化、河道沖淤、土地退化過程、湖泊沼澤化、河道沖淤、沙嘴發(fā)育，以及土壤侵蝕等，使得受幾個(gè)沙嘴發(fā)育，以及土壤侵蝕等，使得受幾個(gè)因素的共同影響，要經(jīng)過若干年才能完成因素的共同影響，要經(jīng)過若干年才能完成的地理過程，采用計(jì)算機(jī)模擬模型，只

25、需的地理過程，采用計(jì)算機(jī)模擬模型，只需幾分鐘就能得出類似的結(jié)果，為資源開發(fā)、幾分鐘就能得出類似的結(jié)果，為資源開發(fā)、國土整治、水土保持、工程論證等提供依國土整治、水土保持、工程論證等提供依據(jù)。據(jù)。 步驟：首先確定土壤侵蝕的數(shù)值分析模型，步驟：首先確定土壤侵蝕的數(shù)值分析模型，根據(jù)模型確定影響土壤侵蝕的因子，這些根據(jù)模型確定影響土壤侵蝕的因子，這些因子必須能夠反映不同的土壤性質(zhì)，不同因子必須能夠反映不同的土壤性質(zhì)，不同的坡面形態(tài)，以及不同的植被條件等，然的坡面形態(tài)，以及不同的植被條件等，然后選擇格網(wǎng)尺寸，建立各個(gè)影響因子的柵后選擇格網(wǎng)尺寸，建立各個(gè)影響因子的柵格數(shù)據(jù)文件，最后將多種信息加以復(fù)合，格數(shù)

26、據(jù)文件，最后將多種信息加以復(fù)合，確定研究地區(qū)土壤侵蝕量的各種不同等級(jí)，確定研究地區(qū)土壤侵蝕量的各種不同等級(jí)，為制定區(qū)域的水土保持規(guī)劃提供依據(jù)。為制定區(qū)域的水土保持規(guī)劃提供依據(jù)。 土壤侵蝕的數(shù)值分析模型隨具體區(qū)域而不同，美國普土壤侵蝕的數(shù)值分析模型隨具體區(qū)域而不同，美國普渡大學(xué)曾根據(jù)渡大學(xué)曾根據(jù)3030余個(gè)觀測站的數(shù)以萬計(jì)的資料，用余個(gè)觀測站的數(shù)以萬計(jì)的資料，用電子計(jì)算機(jī)加以分析，得出下列通用的土壤流失方電子計(jì)算機(jī)加以分析，得出下列通用的土壤流失方程。程。式中：式中：A A為土壤流失量，為土壤流失量， ; ;R R為降雨侵蝕力為降雨侵蝕力; ;K K為土壤可蝕性為土壤可蝕性; ;L L為坡長為坡長; ;S S為坡度為坡度; ;C C為植被覆蓋度為植被覆蓋度: : P P為土壤侵蝕控制措施為土壤侵蝕控制措施。 根據(jù)模型