【正文】 10 11 12 合計 運(yùn)用MATLAB軟件，分別得到2005—2012年8年海桑的預(yù)測值如下表9和圖6，并預(yù)測2013年和2014年海桑數(shù)量（即年度總凋落物量）。表9：8年海桑2005—2012年實際值和預(yù)測值年份2005年2006年2007年2008年2009年2010年2011年2012年2013年2014年實際值 　　預(yù)測值 圖6：8年海桑2005—2012年實際值和預(yù)測值對比由表9和圖6可知，實際值和預(yù)測值非常接近，誤差極小，因此能夠放心預(yù)測。并且海桑數(shù)量并不是持續(xù)增長，而是漲落有序，由此可以看出，海桑并不會給福田紅樹林濕地生態(tài)系統(tǒng)造成巨大的影響，但是也不應(yīng)該持續(xù)引進(jìn)，并且要嚴(yán)格監(jiān)測其數(shù)量變化。、針對海桑入侵的保護(hù)管理建議作為優(yōu)良速生的引進(jìn)樹種，無瓣海桑在華南沿海防護(hù)林體系建設(shè)中發(fā)揮了積極的作用。目前紅樹林人工林中約90％為無瓣海桑林，由于它可以在短期內(nèi)郁閉成林，與其他人工林種類相比能盡快發(fā)揮紅樹林的各種生態(tài)效應(yīng)。我們在綜合考慮預(yù)警結(jié)果的前提下給出如下保護(hù)與管理建議：（1）鑒于無瓣海桑具有較強(qiáng)的種間競爭能力(且仍處于競爭增強(qiáng)階段)和一定的(局部)天然擴(kuò)散更新能力。建議今后在利用無瓣海桑構(gòu)建紅樹林人工林的同時需加強(qiáng)其長期定點的生態(tài)監(jiān)測；（2）深入評價其引種效果，進(jìn)一步開展區(qū)域性栽培試驗、生產(chǎn)力測定及抗逆性評價等；開發(fā)利用無瓣海桑樹種的其他特性，如木材(紙漿)、果膠提取等，以發(fā)揮其多重經(jīng)濟(jì)效益；（3）研究如何改良林分結(jié)構(gòu)以增強(qiáng)其生態(tài)效應(yīng)，尤其是防風(fēng)消浪方面作用；（4）由于純林存在許多弊病，包括易受病蟲害侵襲、景觀單調(diào)、不穩(wěn)定等，建議逐步對現(xiàn)有無瓣海桑純林進(jìn)行改造，修枝后在林下種植適宜的鄉(xiāng)土樹種以形成更穩(wěn)定的混交結(jié)構(gòu)。六、模型檢驗和評價問題一中主要運(yùn)用DPSIR模型，通過構(gòu)建生態(tài)健康評價指標(biāo)體系，進(jìn)而確定指標(biāo)層次和指標(biāo)因子，接著采用層次分析法計算出各指標(biāo)間的權(quán)重，最后通過綜合評價模型得出：福田紅樹林濕地生態(tài)系統(tǒng)的結(jié)構(gòu)層次簡單且穩(wěn)定性相對較差，易受其他因素影響，其穩(wěn)定性在下降，這與事實相符。問題二中通過建立生態(tài)場模型得出海桑的引進(jìn)對周圍生態(tài)因子有所影響，海桑的引進(jìn)使福田紅樹林原生態(tài)系統(tǒng)得生物多樣性增強(qiáng)，增加了食物鏈能量等級，有利于生態(tài)系統(tǒng)穩(wěn)定性的提高。接著基于2005年2012年的月凋落物均值數(shù)據(jù)，通過BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測出20132014年海桑的數(shù)量，通過對比可得，海桑數(shù)量與以往8年相比基本穩(wěn)定，說明海桑的引入未對原生態(tài)系統(tǒng)構(gòu)成生物入侵，這也與事實中濕地生態(tài)系統(tǒng)常通過引進(jìn)海桑提高系統(tǒng)穩(wěn)定性相符，說明我們的模型建立與求解是正確的和符合實際的。七、模型的評價與推廣、模型的優(yōu)點、DPSIR模型的優(yōu)點作為PSR模型的延伸，DPSIR模型的“驅(qū)動力壓力狀態(tài)影響響應(yīng)”部分以分析性因果關(guān)系的方式來體現(xiàn)所有的信息，很好地區(qū)分原因和影響以及人類為了控制影響所采取的響應(yīng)措施,作為分析研究自然生態(tài)系統(tǒng)一一社會相互作用的基礎(chǔ)，允許使用人在不同的實際情況下自定義主要的指標(biāo)。該模型把紅樹林濕地生態(tài)系統(tǒng)與人類的發(fā)展看做是一個人與自然相互作的綜合系統(tǒng)，把影響生態(tài)系統(tǒng)的各個指標(biāo)劃分到相應(yīng)的各個環(huán)節(jié)，層次清晰且考慮全面。、層次分析法的優(yōu)點本題選用層次分析法計算權(quán)重優(yōu)點在于其與DPSIR模型可以較好的結(jié)合，可以在DPSIR模型的基礎(chǔ)之上進(jìn)行定性和定量分析，該方法具有系統(tǒng)、靈活、簡潔的優(yōu)點。、其他方面優(yōu)點（1）首先對數(shù)據(jù)及元素間的相關(guān)性進(jìn)行分析處理，簡化了計算（2）所建立的模型與實際緊密聯(lián)系，由一些利用簡單的模型就能達(dá)到很好的效果，有很好的通用性和推廣性。（3）對數(shù)據(jù)進(jìn)行了無量綱化處理，使得避免了運(yùn)用數(shù)據(jù)時單位不統(tǒng)一的麻煩。（4）運(yùn)用Matlab軟件進(jìn)行計算，可信度高。（5）論文中圖形與數(shù)據(jù)相結(jié)合更具有說服力。、模型的缺點（1）DPSIR模型模型線性的因果關(guān)系過度簡化了實際情況,僅代表了傳統(tǒng)的響應(yīng)式環(huán)境保護(hù)觀念。而且所選指標(biāo)易受數(shù)據(jù)可獲得性的限制，對綜合評價結(jié)果產(chǎn)生不確定性影響。（2）對于層次分析法，當(dāng)指標(biāo)過多時數(shù)據(jù)統(tǒng)計量大，且權(quán)重難以確定，特征值和特征向量的精確求法比較復(fù)雜。（3）在處理數(shù)據(jù)和求解過程中不可避免的出現(xiàn)各種誤差，在一定也影響到模型求解的精確度。、模型的推廣本模型為福田紅樹林自然保護(hù)區(qū)所在的濕地生態(tài)系統(tǒng)構(gòu)建了新型的生態(tài)系統(tǒng)動態(tài)監(jiān)測和健康評估及預(yù)警科學(xué)管理支撐體系，對于生態(tài)系統(tǒng)的發(fā)展具有極其重要的意義。尤其在當(dāng)今以環(huán)境和經(jīng)濟(jì)并重的社會中，該模型的建立對于解決和描述環(huán)境問題及其與社會發(fā)展關(guān)系的問題具有重要的參考價值，同時利用該模型也可以應(yīng)用于水、土壤、海洋等資源的管理保護(hù)以及環(huán)境管理科學(xué)的決策與實施。八、參考文獻(xiàn)[1]伯蓀，廖寶文，王勇軍，等．深圳灣紅樹林生態(tài)系統(tǒng)及其持續(xù)發(fā)展．北京：科學(xué)出版社，2002．1—362．[2] 李玫，. 第3期(總84期):.[3] 曹紅軍. 淺評DPSIR 模型. (增刊):110111+.[4]郭菊蘭，朱耀軍，武高潔，:.[5] 昝啟杰,王勇軍,王伯蓀. 深圳福田紅樹林無瓣海桑+海桑群落N、P、. 22 ( 4): 331 – .[6] 劉莉娜，陳里娥，韋萍萍，李喻春，廖文波，昝啟杰. ：.[7] 駱永明，鄭茂坤，趙其國，滕應(yīng). DPSIR 體系及其在土壤圈環(huán)境管理中的意義. 土 壤. 第38 卷:.附件主要MATLAB程序A=[1,1/2,2。2,1,3。1/2,1/3,1]disp(39。特征值為：39。)。[V,D]=eig(A)。D1=diag(D)。D1=D1.39。disp(39。對應(yīng)的特征向量為：39。)。Vdisp(39。最大特征值為：39。)。m=max(D(:))[n,z]=size(A)。CI=(mn)/(n1)。RI=[0 0 ]。ri=RI(n)。CR=CI/ridisp(39。是否通過一致性檢驗？39。)。if(CR) disp(39。不通過39。)else disp(39。通過39。) [r,c]=find(m==D)。 v=V(:,c)。 disp(39。權(quán)向量為39。) B=v./sum(v)enddisp (’輸入判斷矩陣’）A=input (’ A=’）；[m, n]=size (A)。RI=[0 0 0. 52 0. 89 1. 12 1. 26 1. 36 1. 41 1. 46 1. 49 1. 52 1. 54 1. 56 1. 58 1. 59]。 R=rank(A)。[V, D]=eig (A)。 tz=max (D)。B=max (tz)。[row, col]=find(D==B)。C=V(:,col)。CI=(Bn)/(n1)。CR=CI/RI(1,n)。if CR0. 10disp (’ CI=’）； disp (CI)。disp ( CR=，)。 disp (CR)。disp (’對比矩陣A通過一致性檢驗，各向量權(quán)重向量Q為：’）Q=zeros(n, 1)。 for i=1 :nQ(i, 1)=C(i, 1)/sum(C(:, 1))。endQelsedisp (’對比矩陣未通過一致性檢驗，需要對對比矩陣A重新構(gòu)造’）；end三維繪圖程序 mx=min(x)。 Mx=max(x)。 my=min(y)。 My=max(y)。 Nx=1000。 Ny=1000。 cx=linspace(mx,Mx,Nx)。 cy=linspace(my,My,Ny)。 cz=griddata(x,y,z,cx,cy39。,39。cubic39。)。 mesh(cx,cy,cz) xlabel(39。X坐標(biāo)39。)。ylabel(39。Y坐標(biāo)39。)。zlabel(39。PH值39。)。 title(39。PH值效應(yīng)場39。) title(39。光效應(yīng)場39。)21