李驕　許彥紅　吉靈波

摘要：利用Matlab建模平臺，建立基于Lotka-Volterras生態(tài)學(xué)種群競爭的禿杉數(shù)學(xué)模型，并對其進行Simulink仿真模擬。以云南省騰沖縣禿杉林為研究對象，分析林分內(nèi)禿杉種群與其他樹種種群競爭的關(guān)系，并進行數(shù)值模擬與Logisitic回歸診斷分析，得出3種不同的平衡點競爭模式。分析結(jié)果表明：最理想競爭平衡模式下，林分內(nèi)禿杉種群在樹齡17年之后占據(jù)優(yōu)勢地位，達到了營林的預(yù)期目標；激烈競爭平衡模式下，幼齡林處于第一組平衡的均勢模式；最差競爭平衡模式下，禿杉林在樹齡20年之后處于劣勢地位。本研究改進了傳統(tǒng)競爭指數(shù)的計算方法，從宏觀上對禿杉林種群競爭進行建模仿真分析，并根據(jù)研究結(jié)果對禿杉林后期經(jīng)營提出了應(yīng)對策略及建議，為計算機建模仿真技術(shù)在林學(xué)領(lǐng)域提供了廣闊的應(yīng)用空間。

關(guān)鍵詞：禿杉；Matlab；種群競爭模型；Simulink仿真

中圖分類號： S718.54+1 文獻標志碼： A 文章編號：1002-1302（2014）07-0175-04

收稿日期：2014-03-25

基金項目：國家公益性行業(yè)科研專項（編號：201104053）。

作者簡介：李驕（1982—），男，重慶合川人，碩士研究生，主要從事森林資源管理與評價研究。E-mail：48965058@qq.com。

通信作者：許彥紅，男，云南易門人，碩士，副教授，主要從事森林測計學(xué)、森林資源經(jīng)營管理與評價方面的研究。E-mail：xyhong64@163.com。禿杉（Taiwania flousiana）屬裸子植物杉科臺灣杉屬，是中國特有的世界珍稀植物，現(xiàn)為我國一級保護植物，天然分布于我國湖北省西南部、貴州省東南部及云南省西部等地，垂直分布海拔高度800～2 500 m[1]。目前國內(nèi)對于禿杉林競爭研究只是簡單地通過Hegyi單木競爭方程計算其競爭指數(shù)，并未從宏觀上對禿杉林與不同伴生樹種間的種群競爭進行數(shù)學(xué)建模與分析，無法直觀得出相應(yīng)的競爭干擾程度結(jié)果，從而不能準確地對禿杉林進行經(jīng)營策略指導(dǎo)。本研究引入 Lotka-Volterra 生態(tài)學(xué)數(shù)學(xué)競爭模型對禿杉林種群建模并進行Simulink仿真分析，采用系統(tǒng)、科學(xué)、客觀的分析方法，利用Matlab平臺對實地調(diào)查所得的原始數(shù)據(jù)進行因子指標分類，然后回歸診斷分析得出禿杉與其競爭種群間在不同競爭平衡模式下的內(nèi)在聯(lián)系，為進一步確定合理的禿杉林經(jīng)營密度，制定合理的經(jīng)營管理措施，提高禿杉林經(jīng)營管理水平提供科學(xué)的理論依據(jù)，對禿杉人工林的可持續(xù)經(jīng)營具有十分重要的意義。同時，本研究是計算機建模仿真技術(shù)在林業(yè)領(lǐng)域應(yīng)用的重要發(fā)展方向。

1Matlab建模平臺與Simulink仿真

Matlab是matrix與laboratory 2個詞的組合，意為矩陣工廠（矩陣實驗室），是主要面對科學(xué)計算、可視化以及交互式程序設(shè)計的高科技計算環(huán)境。借助于Matlab強大的數(shù)據(jù)處理、圖形處理能力，可以方便快捷高效地解決數(shù)學(xué)建模中的各種問題[2]。Simulink是Matlab中的一種可視化仿真工具，是一種基于Matlab的框圖設(shè)計環(huán)境，被廣泛應(yīng)用于線性系統(tǒng)、非線性系統(tǒng)、數(shù)字控制及數(shù)字信號處理的建模和仿真中[3]。Simulink可以用連續(xù)采樣時間、離散采樣時間或2種混合的采樣時間進行建模。Simulink與Matlab緊密集成，可以直接訪問Matlab大量的工具來進行算法的研發(fā)、仿真的分析和可視化、批處理腳本的創(chuàng)建、建模環(huán)境的定制以及信號參數(shù)和測試數(shù)據(jù)的定義[4]。

2種群競爭數(shù)學(xué)模型與仿真結(jié)構(gòu)圖構(gòu)建

2.1種群競爭模型建立原理

4禿杉種群競爭平衡模式對應(yīng)的經(jīng)營策略

通過研究禿杉種群競爭模型仿真與3種種群競爭模式得出了以下禿杉林人工經(jīng)營策略。

（1）要達到最優(yōu)競爭模式（P1平衡點），需要提升禿杉種群N1所占的份額。首先要做好禿杉幼林郁閉前的撫育工作，幼林撫育的內(nèi)容主要應(yīng)從土壤管理入手，通過松土、除草改善土壤的理化性質(zhì)，排除雜草、灌木對幼林的競爭[19]；其次對林木本身進行必要的抑制調(diào)節(jié)，如除蘗、平茬、間苗等，使幼林成林，到提升N1的效果。

（2）針對最差競爭模式（P2平衡點）的情況，需要降低N2所占份額。在營林上可以考慮采用撫育間伐與合理混交2種措施。幼林郁閉后，林木的個體逐漸增大，主要表現(xiàn)在冠幅、樹高、直徑方面的增加，林木與林木、林木與灌木或草本的競爭隨之逐漸增加，這時需要對密度較大的禿杉林分進行適度

撫育間伐。若間伐強度過大，雖然可以減輕禿杉種群與競爭種群之間的競爭，但林分單位面積的生長量會減少；若間伐強度過小，則達不到撫育效果，解決不了林木激烈競爭的趨勢。此外，還要及時修枝，對不利于優(yōu)良木生長的病腐木、被壓木、生長畸型的林木及時清理[20]，使禿杉有一個良好的生長環(huán)境。造林密度和混交樹種及比例影響林分生長發(fā)育全過程。造林密度不同，混交樹種及比例不同，林分郁閉、林木開始激烈競爭的時間及激烈程度不同。密度過大，郁閉越早，競爭時間越提前，會影響林分生長；密度過小，林木個體在單位面積上株數(shù)較少，不能很好利用營養(yǎng)空間。而禿杉與伴生樹種種群混交密度應(yīng)該遵循參數(shù)優(yōu)化后禿杉種群競爭方程組中的標準化種群密度u1和u2來確定。

（3）林地內(nèi)環(huán)境的保護管理。在研究的數(shù)學(xué)模型中，禿杉種群與其他競爭種群競爭所需的前提環(huán)境資源如陽光、水源、土壤等為固定值，所以在林地內(nèi)應(yīng)保護好地被物，禁止人為干預(yù)枯枝落葉層，以免破壞森林內(nèi)營養(yǎng)元素的循環(huán)。如果林內(nèi)沒有枯枝落葉，營養(yǎng)元素的循環(huán)得不到保證，林地涵養(yǎng)水源能力差，導(dǎo)致地表水大量流失，土壤含水量低；同時由于地表徑流、土壤流失大，林木所需的養(yǎng)分和水分不足，從而加劇林木的競爭。

5結(jié)論

隨著計算機技術(shù)的不斷發(fā)展，計算機在林業(yè)領(lǐng)域中的應(yīng)用已經(jīng)從相對簡單的多媒體計算機技術(shù)時代進入了更高層次的計算機建模仿真技術(shù)時代。本研究利用生態(tài)學(xué)數(shù)學(xué)模型研究禿杉種群與其他競爭種群空間資源爭奪結(jié)果，并建立了3種競爭模式。本研究通過Matlab進行數(shù)值模擬與Simulink仿真分析，得出3個種群競爭平衡點模式，并基于模型參數(shù)數(shù)值模擬圖給出了騰沖縣禿杉人工林經(jīng)營策略建議。endprint

本研究是國家公益性行業(yè)科研專項“禿杉大徑材人工培育關(guān)鍵技術(shù)研究（201104053）”的內(nèi)容之一，研究成果已應(yīng)用于試驗示范區(qū)建設(shè)。隨著數(shù)字化林業(yè)、林業(yè)信息化的推廣，計算機仿真模擬技術(shù)將會在林業(yè)行業(yè)得到更加廣泛的應(yīng)用。

參考文獻：

[1]陳強，袁明，劉云彩，等. 禿杉的物種確立、天然林種群特征、保護、引種和種源選擇研究[J]. 西部林業(yè)科學(xué)，2012，41（2）：1-16.

[2]席偉，微分方程方向場MATLAB仿真工具箱設(shè)計[J]. 信息安全與技術(shù)，2012（11）：65-67.

[3]Vahidi B，Esmaeeli E. MATLAB-SIMULINK-based simulation for digital differential relay protection of power transformer for educational purpose[J]. Computer Applications in Engineering Education，2013，21（3）：475-483.

[4]劉斐. 基于Simulink的通信系統(tǒng)仿真[J]. 電腦知識與技術(shù)，2011，7（7）：1593-1594.

[5]李海峰，李秋英，陸民燕.基于Logistic測試覆蓋率函數(shù)的軟件可靠性建模研究[J]. 計算機研究與發(fā)展，2011，48（2）：232-240.

[6]黃光球，趙魏娟，陸秋琴.基于3種群Lotka-Volterra模型的種群動力學(xué)函數(shù)優(yōu)化算法[J]. 計算機科學(xué)，2013，40（8）：214-219.

[7]常振海，劉薇. Logistic回歸模型及其應(yīng)用[J]. 延邊大學(xué)學(xué)報：自然科學(xué)版，2012，38（1）：28-32.

[8]王金良，李慧鳳. Logistic種群演化模型的漸近加權(quán)周期性[J]. 應(yīng)用數(shù)學(xué)學(xué)報，2011，34（3）：496-501.

[9]李敏. 真實模型Matlab仿真的過程控制系統(tǒng)實驗研究[J]. 浙江工業(yè)大學(xué)學(xué)報，2011，39（2）：168-173.

[10]張金鋮，李傳東，李超輩. 基于Matlab的憶感器建模仿真及應(yīng)用[J]. 固體電子學(xué)研究與進展，2012（3）：239-245.

[11]Zuliani P，Platzer A，Clarke E M. Bayesian statistical model checking with application to Stateflow/Simulink verification[J]. Formal Methods in System Design，2013，43（2，SI）：338-367.

[12]郝連秀，王浩軒，張亮. 基于Matlab/Simulink的頻率合成器模型設(shè)計[J]. 計算機與現(xiàn)代化，2013（5）：206-210.

[13]劉佩云. 禿杉人工林立地指數(shù)的研究[J]. 林業(yè)科技，2011，36（6）：16-18.

[14]李典，楊永. 點隱式龍格-庫塔方法的應(yīng)用研究[J]. 航空計算技術(shù)，2011，41（3）：66-70.

[15]韓慎友，鐘青. 三階龍格-庫塔時間分裂顯式算法的誤差分析[J]. 氣象學(xué)報，2012，70（1）：119-127.

[16]Medina A，Segundo-Ramirez J，Ribeiro P，et al. Harmonic analysis in frequency and time domain[J]. IEEE Transactions on Power Delivery，2013，28（3）：1813-1821.

[17]趙峰. 貴州雷公山禿杉優(yōu)勢種群的生態(tài)位特征[J]. 中國農(nóng)學(xué)通報，2012，28（1）：17-23.

[18]Zhou M. Using logistic regression models for selection in non-replicated sugarcane breeding populations[J]. Euphytica，2013，191（3）：415-428.

[19]雒志學(xué)，郭金生. 具有擴散和年齡結(jié)構(gòu)競爭種群的最優(yōu)收獲[J]. 數(shù)學(xué)進展，2009，38（2）：209-219.

[20]陳志陽，楊寧，姚先銘，等. 貴州雷公山禿杉種群生活史特征與空間分布格局[J]. 生態(tài)學(xué)報，2012，32（7）：2158-2165.endprint

本研究是國家公益性行業(yè)科研專項“禿杉大徑材人工培育關(guān)鍵技術(shù)研究（201104053）”的內(nèi)容之一，研究成果已應(yīng)用于試驗示范區(qū)建設(shè)。隨著數(shù)字化林業(yè)、林業(yè)信息化的推廣，計算機仿真模擬技術(shù)將會在林業(yè)行業(yè)得到更加廣泛的應(yīng)用。

參考文獻：

[1]陳強，袁明，劉云彩，等. 禿杉的物種確立、天然林種群特征、保護、引種和種源選擇研究[J]. 西部林業(yè)科學(xué)，2012，41（2）：1-16.

[2]席偉，微分方程方向場MATLAB仿真工具箱設(shè)計[J]. 信息安全與技術(shù)，2012（11）：65-67.

[3]Vahidi B，Esmaeeli E. MATLAB-SIMULINK-based simulation for digital differential relay protection of power transformer for educational purpose[J]. Computer Applications in Engineering Education，2013，21（3）：475-483.

[4]劉斐. 基于Simulink的通信系統(tǒng)仿真[J]. 電腦知識與技術(shù)，2011，7（7）：1593-1594.

[5]李海峰，李秋英，陸民燕.基于Logistic測試覆蓋率函數(shù)的軟件可靠性建模研究[J]. 計算機研究與發(fā)展，2011，48（2）：232-240.

[6]黃光球，趙魏娟，陸秋琴.基于3種群Lotka-Volterra模型的種群動力學(xué)函數(shù)優(yōu)化算法[J]. 計算機科學(xué)，2013，40（8）：214-219.

[7]常振海，劉薇. Logistic回歸模型及其應(yīng)用[J]. 延邊大學(xué)學(xué)報：自然科學(xué)版，2012，38（1）：28-32.

[8]王金良，李慧鳳. Logistic種群演化模型的漸近加權(quán)周期性[J]. 應(yīng)用數(shù)學(xué)學(xué)報，2011，34（3）：496-501.

[9]李敏. 真實模型Matlab仿真的過程控制系統(tǒng)實驗研究[J]. 浙江工業(yè)大學(xué)學(xué)報，2011，39（2）：168-173.

[10]張金鋮，李傳東，李超輩. 基于Matlab的憶感器建模仿真及應(yīng)用[J]. 固體電子學(xué)研究與進展，2012（3）：239-245.

[11]Zuliani P，Platzer A，Clarke E M. Bayesian statistical model checking with application to Stateflow/Simulink verification[J]. Formal Methods in System Design，2013，43（2，SI）：338-367.

[12]郝連秀，王浩軒，張亮. 基于Matlab/Simulink的頻率合成器模型設(shè)計[J]. 計算機與現(xiàn)代化，2013（5）：206-210.

[13]劉佩云. 禿杉人工林立地指數(shù)的研究[J]. 林業(yè)科技，2011，36（6）：16-18.

[14]李典，楊永. 點隱式龍格-庫塔方法的應(yīng)用研究[J]. 航空計算技術(shù)，2011，41（3）：66-70.

[15]韓慎友，鐘青. 三階龍格-庫塔時間分裂顯式算法的誤差分析[J]. 氣象學(xué)報，2012，70（1）：119-127.

[16]Medina A，Segundo-Ramirez J，Ribeiro P，et al. Harmonic analysis in frequency and time domain[J]. IEEE Transactions on Power Delivery，2013，28（3）：1813-1821.

[17]趙峰. 貴州雷公山禿杉優(yōu)勢種群的生態(tài)位特征[J]. 中國農(nóng)學(xué)通報，2012，28（1）：17-23.

[18]Zhou M. Using logistic regression models for selection in non-replicated sugarcane breeding populations[J]. Euphytica，2013，191（3）：415-428.

[19]雒志學(xué)，郭金生. 具有擴散和年齡結(jié)構(gòu)競爭種群的最優(yōu)收獲[J]. 數(shù)學(xué)進展，2009，38（2）：209-219.

[20]陳志陽，楊寧，姚先銘，等. 貴州雷公山禿杉種群生活史特征與空間分布格局[J]. 生態(tài)學(xué)報，2012，32（7）：2158-2165.endprint

本研究是國家公益性行業(yè)科研專項“禿杉大徑材人工培育關(guān)鍵技術(shù)研究（201104053）”的內(nèi)容之一，研究成果已應(yīng)用于試驗示范區(qū)建設(shè)。隨著數(shù)字化林業(yè)、林業(yè)信息化的推廣，計算機仿真模擬技術(shù)將會在林業(yè)行業(yè)得到更加廣泛的應(yīng)用。

參考文獻：

[1]陳強，袁明，劉云彩，等. 禿杉的物種確立、天然林種群特征、保護、引種和種源選擇研究[J]. 西部林業(yè)科學(xué)，2012，41（2）：1-16.

[2]席偉，微分方程方向場MATLAB仿真工具箱設(shè)計[J]. 信息安全與技術(shù)，2012（11）：65-67.

[3]Vahidi B，Esmaeeli E. MATLAB-SIMULINK-based simulation for digital differential relay protection of power transformer for educational purpose[J]. Computer Applications in Engineering Education，2013，21（3）：475-483.

[4]劉斐. 基于Simulink的通信系統(tǒng)仿真[J]. 電腦知識與技術(shù)，2011，7（7）：1593-1594.

[5]李海峰，李秋英，陸民燕.基于Logistic測試覆蓋率函數(shù)的軟件可靠性建模研究[J]. 計算機研究與發(fā)展，2011，48（2）：232-240.

[6]黃光球，趙魏娟，陸秋琴.基于3種群Lotka-Volterra模型的種群動力學(xué)函數(shù)優(yōu)化算法[J]. 計算機科學(xué)，2013，40（8）：214-219.

[7]常振海，劉薇. Logistic回歸模型及其應(yīng)用[J]. 延邊大學(xué)學(xué)報：自然科學(xué)版，2012，38（1）：28-32.

[8]王金良，李慧鳳. Logistic種群演化模型的漸近加權(quán)周期性[J]. 應(yīng)用數(shù)學(xué)學(xué)報，2011，34（3）：496-501.

[9]李敏. 真實模型Matlab仿真的過程控制系統(tǒng)實驗研究[J]. 浙江工業(yè)大學(xué)學(xué)報，2011，39（2）：168-173.

[10]張金鋮，李傳東，李超輩. 基于Matlab的憶感器建模仿真及應(yīng)用[J]. 固體電子學(xué)研究與進展，2012（3）：239-245.

[11]Zuliani P，Platzer A，Clarke E M. Bayesian statistical model checking with application to Stateflow/Simulink verification[J]. Formal Methods in System Design，2013，43（2，SI）：338-367.

[12]郝連秀，王浩軒，張亮. 基于Matlab/Simulink的頻率合成器模型設(shè)計[J]. 計算機與現(xiàn)代化，2013（5）：206-210.

[13]劉佩云. 禿杉人工林立地指數(shù)的研究[J]. 林業(yè)科技，2011，36（6）：16-18.

[14]李典，楊永. 點隱式龍格-庫塔方法的應(yīng)用研究[J]. 航空計算技術(shù)，2011，41（3）：66-70.

[15]韓慎友，鐘青. 三階龍格-庫塔時間分裂顯式算法的誤差分析[J]. 氣象學(xué)報，2012，70（1）：119-127.

[16]Medina A，Segundo-Ramirez J，Ribeiro P，et al. Harmonic analysis in frequency and time domain[J]. IEEE Transactions on Power Delivery，2013，28（3）：1813-1821.

[17]趙峰. 貴州雷公山禿杉優(yōu)勢種群的生態(tài)位特征[J]. 中國農(nóng)學(xué)通報，2012，28（1）：17-23.

[18]Zhou M. Using logistic regression models for selection in non-replicated sugarcane breeding populations[J]. Euphytica，2013，191（3）：415-428.

[19]雒志學(xué)，郭金生. 具有擴散和年齡結(jié)構(gòu)競爭種群的最優(yōu)收獲[J]. 數(shù)學(xué)進展，2009，38（2）：209-219.

[20]陳志陽，楊寧，姚先銘，等. 貴州雷公山禿杉種群生活史特征與空間分布格局[J]. 生態(tài)學(xué)報，2012，32（7）：2158-2165.endprint