喬一娜　楊偉　石孟迪　張錫唐　秦于倩　郝新忠　石長(zhǎng)春

摘要 為研究榆林沙區(qū)樟子松胸徑-樹(shù)高生長(zhǎng)模型，填補(bǔ)榆林沙區(qū)樟子松在胸徑-樹(shù)高生長(zhǎng)模型研究的空白，為榆林沙區(qū)樟子松林分生長(zhǎng)預(yù)測(cè)、生產(chǎn)經(jīng)營(yíng)管理提供參考依據(jù)，以榆林沙區(qū)樟子松人工林為研究對(duì)象，選用15個(gè)常用的胸徑-樹(shù)高生長(zhǎng)模型，通過(guò)5個(gè)評(píng)價(jià)指標(biāo)（R、決定系數(shù)（RMSE）、平均絕對(duì)殘差（MAE）、殘差平均和（SSE）、Akaike信息量準(zhǔn)則（AIC））進(jìn)行對(duì)比分析擬合，從而進(jìn)一步確定，適宜的榆林沙區(qū)樟子松胸徑-樹(shù)高生長(zhǎng)模型。結(jié)果表明：除了Gompertz（1825）模型（14號(hào)模型）無(wú)法輸出參數(shù)，剩余的模型均可以。在剩余的14個(gè)模型中，雙曲線型（2號(hào)模型）、混合型（5號(hào)模型）、二次多項(xiàng)式（6號(hào)模型）、Korf（10號(hào)模型）、修正Veibull（11號(hào)模型）這5個(gè)模型評(píng)價(jià)指標(biāo)較優(yōu)，擬合效果較好。綜合比較分析可知，這5種基礎(chǔ)胸徑-樹(shù)高模型中，修正Veibull模型可以更好地?cái)M合樟子松胸徑-樹(shù)高的關(guān)系，精度也比較準(zhǔn)確，建議選用修正Veibull模型。

關(guān)鍵詞榆林沙區(qū)；樟子松；胸徑-樹(shù)高生長(zhǎng)模型

中圖分類號(hào)S791.253文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)0517-6611（2023）08-0137-05

doi：10.3969/j.issn.0517-6611.2023.08.031開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Study on DBH-Height Growth Model of Pinus sylvestris Plantation in Yulin Sandy Area

QIAO Yi-na，YANG Wei，SHI Meng-di et al（Shaanxi Academy of Forestry，Xi'an，Shaanxi 710082）

AbstractStudying the DBH-tree height growth model of Pinus sylvestris in Yulin sandy area can fill the gap in the research on DBH-tree height growth model of Pinus sylvestris in Yulin sandy area.It can also provide a reference for the growth prediction and production management of Pinus sylvestris in the sandy area of Yulin.This study takes the Pinus sylvestris plantation in Yulin sandy area as the research object，selects 15 commonly used DBH-tree height growth models.Through 5 evaluation indicators （R，coefficient of determination （RMSE），mean absolute residual error （MAE），residual mean sum （SSE），Akaike information content criterion （AIC）） for comparative analysis and fitting.Therefore，the most suitable DBH-tree height growth model of Pinus sylvestris in the sandy area of Yulin is further determined.The result shows：except for the Gompertz （1825） model （model number 14），which cannot output parameters，the rest of the models can.Among the remaining 14 models，these five model evaluation indicators are better，and the fitting effect is better，including the hyperbolic model （Model No.2），the hybrid model （Model No.5），the quadratic polynomial model （Model No.6），the Korf model （Model No.10），the modified Veibull model （Model No.11））.Finally，a comprehensive comparative analysis shows that the modified Veibull model among the five basic DBH-tree height models can better fit the relationship between the DBH and tree height of Pinus sylvestris，and the accuracy is relatively accurate. It is recommended to use the modified Veibull model.

Key wordsYulin Sand District;Pinus sylvetris;DBH-tree height growth model

樟子松（Pinus sylvetris var.mongolica Litv.）是沙地上最穩(wěn)定最持久，生產(chǎn)力最高的樹(shù)種之一，在榆林地區(qū)廣泛栽培，現(xiàn)治沙造林面積已達(dá)到6.67余萬(wàn)hm，部分林區(qū)已郁閉成林，表現(xiàn)出很強(qiáng)的抗逆性和適應(yīng)性，顯現(xiàn)出極強(qiáng)的防風(fēng)固沙效果［1-3］。由此可見(jiàn)，研究榆林地區(qū)樟子松林對(duì)改善榆林沙區(qū)生態(tài)環(huán)境，提高林分質(zhì)量，增強(qiáng)防護(hù)功能等起著重要作用［4-6］。胸徑和樹(shù)高作為林分?jǐn)?shù)量和質(zhì)量調(diào)查的重要指標(biāo)［7］，在林業(yè)調(diào)查各因素之間有著密切的聯(lián)系［8］。胸徑能夠反映出林木生長(zhǎng)情況［9］，樹(shù)高則能夠描述林分生長(zhǎng)變化及當(dāng)前環(huán)境立地質(zhì)量的好壞［10］，因此胸徑和樹(shù)高是計(jì)算植物生物量不可缺少的因素［11］。但在實(shí)際中，樹(shù)高數(shù)據(jù)相對(duì)來(lái)說(shuō)較難獲取，尤其在林分密度較大的情況下，由于視野受限，測(cè)定的數(shù)據(jù)未必準(zhǔn)確，存在較大誤差。因此，研究者通過(guò)建立不同立地類型不同樹(shù)種胸徑-樹(shù)高生長(zhǎng)預(yù)測(cè)模型進(jìn)行了樹(shù)高的估測(cè)［12-14］，以解決實(shí)際測(cè)量的難題，為不同地區(qū)不同樹(shù)種的林木生長(zhǎng)預(yù)測(cè)、森林經(jīng)營(yíng)與管理提供參考依據(jù)。

目前已有一些對(duì)樟子松胸徑-樹(shù)高生長(zhǎng)模型的研究，但總體來(lái)說(shuō)相對(duì)較少。其中，李善堯［15］對(duì)內(nèi)蒙古紅花爾基樟子松天然林胸徑和樹(shù)高進(jìn)行了研究，應(yīng)用了7個(gè)常見(jiàn)胸徑-樹(shù)高模型，得出最適合樟子松天然林胸徑-樹(shù)高模型為冪函數(shù)曲線模型。羅玲等［16］用Logistic方程進(jìn)行了榆林沙區(qū)樟子松樹(shù)高和胸徑生長(zhǎng)過(guò)程的模擬，并推導(dǎo)出不同立地類型條件下樟子松的生長(zhǎng)量方程，擬合精度較好。吳夢(mèng)婉［17］對(duì)章古臺(tái)地區(qū)樟子松人工林的樹(shù)高、胸徑和材積生長(zhǎng)過(guò)程進(jìn)行了擬合，發(fā)現(xiàn)Logistic方程對(duì)樟子松人工林?jǐn)M合效果較好。雷澤勇等［18］研究發(fā)現(xiàn)，樟子松的生長(zhǎng)能夠運(yùn)用Richards模型很好地模擬，符合實(shí)際情況。然而，對(duì)于全國(guó)最早引種樟子松的榆林沙區(qū)樟子松相關(guān)生長(zhǎng)模型研究卻鮮見(jiàn)報(bào)道。鑒于此，筆者研究了榆林沙區(qū)樟子松胸徑-樹(shù)高生長(zhǎng)模型，構(gòu)建了準(zhǔn)確的胸徑-樹(shù)高關(guān)系，以期為榆林沙區(qū)樟子松的生長(zhǎng)推廣、經(jīng)營(yíng)管理提供理論依據(jù)。

1材料與方法

1.1研究區(qū)概況該研究區(qū)域位于陜西省榆林市榆陽(yáng)區(qū)（109°42′E，38°21′N），海拔560～1 907 m，平均海拔為1 000 m。其地勢(shì)東北高，中南低，以明長(zhǎng)城為界，長(zhǎng)城以北為風(fēng)沙草灘區(qū)，長(zhǎng)城以南為丘陵溝壑區(qū)。地貌類型主要為固定沙地，沙丘起伏不大，為7～10 m。氣候?qū)贉貛О敫珊荡箨懶约撅L(fēng)氣候，春季干燥，風(fēng)大風(fēng)多，秋、冬、春3季多為西北風(fēng)；夏季炎熱干燥，多為東南風(fēng)；秋季地面迅速冷卻，上熱下冷，空氣日趨穩(wěn)定，天氣比較晴朗涼爽；冬季漫長(zhǎng)、干燥、寒冷、少雪。全年降水量400 mm左右，主要集中在7—9月，年平均蒸發(fā)量為2 388.7 mm，一般為降水量的5倍。年均氣溫7.9 ℃左右，無(wú)霜期150～180 d，全年日照時(shí)數(shù)一般在2 700 h以上。

土壤類型主要為風(fēng)沙土，植被類型主要喬木樹(shù)種為樟子松（Pinus sylvestris var. mongolica Litv.），林下灌木主要有沙柳（Salix cheilophila）、花棒（Corethrodendron scoparium Fisch. et Basiner）、紫穗槐（Amorpha fruticosa Linn.）、檸條錦雞兒（Caragana korshinskii Kom.）、臭柏（Sabina vulgaris Ant.），主要草本植物有羊柴（Hedysarum mongdicum Turcz Var.）、沙蒿（Artemisia desertorum Spreng. Syst. Veg.）、苜蓿（Medicago Sativa Linn）等。

1.2調(diào)查樣地情況該調(diào)查樣地選在陜西省樟子松良種基地，采取典型抽樣的方法進(jìn)行標(biāo)準(zhǔn)樣地調(diào)查，樣地面積為20 m×20 m，共調(diào)查10塊樣地進(jìn)行每木檢尺。使用測(cè)樹(shù)鋼圍尺、塔尺測(cè)量每株樟子松的胸徑（D）和樹(shù)高（H），在進(jìn)行模型構(gòu)建時(shí)，胸徑和樹(shù)高測(cè)量數(shù)據(jù)精確到0.01 m。每塊樣地的樟子松年齡均為14年，株數(shù)密度為500株/hm，共測(cè)量300株樟子松的生長(zhǎng)量情況。樣地測(cè)量數(shù)據(jù)具體詳見(jiàn)表1，樟子松的胸徑-樹(shù)高散點(diǎn)圖見(jiàn)圖1。

1.3建模方法根據(jù)國(guó)內(nèi)外關(guān)于胸徑-樹(shù)高模型的研究［19-21］，從中篩選出15個(gè)常用的胸徑-樹(shù)高模型進(jìn)行擬合，通過(guò)對(duì)比分析比較并從中遴選出適宜榆林沙區(qū)樟子松林的胸徑-樹(shù)高模型，并通過(guò)SPSS 20.0進(jìn)行模型構(gòu)建和參數(shù)檢驗(yàn)（表2）。

選用5個(gè)評(píng)價(jià)指標(biāo)進(jìn)行模型評(píng)價(jià)及檢驗(yàn)（表3），分別為決定系數(shù)（R）、均方根誤差（RMSE）、平均絕對(duì)誤差（MAE）、殘差平均和（SSE）、Akaike信息量準(zhǔn)則（AIC），通過(guò)5個(gè)指標(biāo)綜合評(píng)價(jià)篩選出擬合效果最佳的模型。

1.4數(shù)據(jù)處理及分析運(yùn)用Excel 2010進(jìn)行數(shù)據(jù)整理與統(tǒng)計(jì)，運(yùn)用SPSS 20.0進(jìn)行回歸分析及模型構(gòu)建。

2結(jié)果與分析

2.1榆林沙區(qū)樟子松基礎(chǔ)樹(shù)高-胸徑生長(zhǎng)模型的構(gòu)建利用SPSS 20.0對(duì)以上15個(gè)基礎(chǔ)胸徑-樹(shù)高模型進(jìn)行參數(shù)求解及模型構(gòu)建，發(fā)現(xiàn)第14個(gè)模型無(wú)法求解，剩余的模型均能求得參數(shù)。14個(gè)生長(zhǎng)模型的參數(shù)估計(jì)和擬合估計(jì)見(jiàn)表4、5。

通過(guò)對(duì)建模數(shù)據(jù)進(jìn)行分析及擬合估計(jì)分析可以看出，這14個(gè)模型擬合效果R均在0.7以上，其中2、5、6、10、11的模型擬合效果較好，均為0.744，7、8的模型擬合效果較差，僅為0.708?？傮w來(lái)說(shuō)，這14個(gè)模型的擬合R差異不大，并由大到小排序?yàn)?、5、6、10、11、12、13、3、4、14、9、1、7、8。對(duì)決定系數(shù)（RMSE）進(jìn)行擬合發(fā)現(xiàn)，14個(gè)模型的RMSE值在0.323～0.346，其中7、8號(hào)模型的RMSE值最高，2、5、6、10、11號(hào)模型的RMSE值最低。14個(gè)模型的平均絕對(duì)殘差（MAE）值為0.018 8～0.020 1，其大小排序與RMSE值相同。對(duì)Akaike信息量準(zhǔn)則（AIC）也進(jìn)行擬合計(jì)算，發(fā)現(xiàn)擬合值均為負(fù)值，且范圍在-526.86～-507.72（表4）。

不同樹(shù)高-胸徑模型對(duì)樟子松樹(shù)高的參數(shù)檢驗(yàn)分析可以看出，12號(hào)模型的c、13號(hào)模型的b的P值分別為0.005、0.017，說(shuō)明這2個(gè)值有顯著差異。剩余模型的值均為0.000，說(shuō)明有極顯著性差異。

通過(guò)分析比較R、RMSE、MAE、SSE、AIC等值，選擇出雙曲線型（2號(hào)模型）、混合型（5號(hào)模型）、二次多項(xiàng)式（6號(hào)模型）、Korf（10號(hào)模型）、修正Veibull（11號(hào)模型）這5個(gè)模型，這5個(gè)模型評(píng)價(jià)指標(biāo)較優(yōu)，擬合效果較好，優(yōu)于其他模型。因此，選用擬合效果較優(yōu)的以上5個(gè)模型進(jìn)行樟子松胸徑-樹(shù)高關(guān)系擬合，擬合效果如圖2所示。從圖2可以看出，這5種基礎(chǔ)胸徑-樹(shù)高模型擬合效果相似，擬合程度幾乎相同（R均為0.744），說(shuō)明這5種模型均可以很好地描述榆林地區(qū)樟子松胸徑-樹(shù)高的關(guān)系。

2.2榆林沙區(qū)樟子松基礎(chǔ)樹(shù)高-胸徑生長(zhǎng)優(yōu)選模型預(yù)測(cè)效果比較為了從這5種基礎(chǔ)胸徑-樹(shù)高模型中篩選出更適合榆林地區(qū)樟子松胸徑的擬合，將以上5種模型進(jìn)行了進(jìn)一步對(duì)比檢驗(yàn)，將通過(guò)R2、RMSE、MAE、SSE、AIC等指標(biāo)值評(píng)價(jià)，從而進(jìn)一步確定最適合榆林地區(qū)樟子松基礎(chǔ)胸徑-樹(shù)高模型（表6）。

由表6可知，這5種模型的R、RMSE、MAE均相同，難以確定出最適合的基礎(chǔ)胸徑-樹(shù)高模型，只有通過(guò)SSE、AIC進(jìn)行分析確定，比較SSE大小可知，修正Veibull模型＞混合型模型＞二次多項(xiàng)式模型＞Korf模型＞雙曲線型模型，比較AIC絕對(duì)值大小可知，同樣表現(xiàn)為修正Veibull模型＞混合型模型＞二次多項(xiàng)式模型＞Korf模型＞雙曲線型模型。綜合比較分析可知，在這5種基礎(chǔ)胸徑-樹(shù)高模型中修正Veibull模型可以更好地?cái)M合樟子松胸徑-樹(shù)高的關(guān)系，精度也比較準(zhǔn)確。

為了進(jìn)一步證明修正Veibull模型在樟子松胸徑-樹(shù)高關(guān)系表達(dá)的準(zhǔn)確性，對(duì)其預(yù)測(cè)值與殘差值進(jìn)行作圖。由圖3可以看出，殘差值分布均勻，無(wú)明顯異質(zhì)性，表明用修正Veibull模型擬合效果較好，適合榆林地區(qū)樟子松基礎(chǔ)胸徑-樹(shù)高模型計(jì)算，可以在實(shí)踐中應(yīng)用。

3討論

林分生長(zhǎng)模型是研究林業(yè)和林木生長(zhǎng)的基礎(chǔ)，可為研究林木生長(zhǎng)提供精準(zhǔn)的林分生長(zhǎng)預(yù)測(cè)［22-25］。胸徑和樹(shù)高都是林木調(diào)查中基本的調(diào)查因子，同時(shí)也是衡量林分結(jié)構(gòu)和林木質(zhì)量的重要指標(biāo)［7］，在林業(yè)調(diào)查中具有不可替代的實(shí)際意義。我國(guó)對(duì)于林木基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型的研究較早，相對(duì)比較成熟，且在不同樹(shù)種不同地區(qū)不同樹(shù)高的基礎(chǔ)生長(zhǎng)模型方面有很深的研究，李明華等［26］將上海市主要造林樹(shù)種（池杉、柳樹(shù)、木蘭類、女貞、水杉、楊樹(shù)、榆樹(shù)類）進(jìn)行了模型構(gòu)建，得到各樹(shù)種對(duì)應(yīng)的最優(yōu)模型結(jié)果；婁明華等［27］以寧波地區(qū)常見(jiàn)的石櫟-木荷天然常綠落葉混交林為研究對(duì)象，通過(guò)5個(gè)模型評(píng)價(jià)指標(biāo)分析擬合比較，進(jìn)而確定最適宜的石櫟-木荷天然常綠落葉混交林胸徑-樹(shù)高生長(zhǎng)模型；白天恒［28］通過(guò)對(duì)章古臺(tái)地區(qū)樟子松人工固沙林進(jìn)行模型構(gòu)建和擬合，得出Richards模型的擬合精度最高，并適宜在章古臺(tái)地區(qū)進(jìn)行推廣應(yīng)用。筆者在前人研究的基礎(chǔ)上，在榆林地區(qū)的樟子松良種基地進(jìn)行胸徑、樹(shù)高等試驗(yàn)數(shù)據(jù)的調(diào)查，并選用了15個(gè)基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型進(jìn)行擬合和檢驗(yàn)，以此來(lái)研究榆林地區(qū)樟子松的生長(zhǎng)過(guò)程。通過(guò)對(duì)比分析15個(gè)胸徑-樹(shù)高生長(zhǎng)模型的擬合精度，分析得出修正Veibull模型擬合效果較好。這一結(jié)論與前人的研究結(jié)果不一致［17-18］，原因可能是研究對(duì)象不同所致，但也有可能是研究地區(qū)不一致所致。該研究的研究對(duì)象為榆林地區(qū)樟子松人工林，其胸徑為1.62～6.37 cm，進(jìn)一步說(shuō)明模型的適應(yīng)性可能與環(huán)境、樹(shù)種、經(jīng)緯度等外界因素有關(guān)，這就需要后續(xù)進(jìn)行進(jìn)一步研究。

通過(guò)以上對(duì)榆林地區(qū)樟子松人工林的基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型進(jìn)行研究，從中篩選出擬合效果最優(yōu)的基礎(chǔ)胸徑-樹(shù)高模型（修正Veibull模型），填補(bǔ)了榆林地區(qū)關(guān)于樟子松人工林基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型的空白，對(duì)把握榆林地區(qū)樟子松林分動(dòng)態(tài)及預(yù)估榆林地區(qū)樟子松林分生長(zhǎng)狀態(tài)具有極為重要的意義。因此，可以利用修正Veibull模型在榆林地區(qū)進(jìn)行研究應(yīng)用，以監(jiān)測(cè)樟子松的正常生長(zhǎng)。

4結(jié)論

該研究以榆林地區(qū)樟子松人工林為研究對(duì)象，利用陜西省樟子松良種基地14年樟子松10塊標(biāo)準(zhǔn)樣地的胸徑、樹(shù)高數(shù)據(jù)進(jìn)行模型構(gòu)建及參數(shù)檢驗(yàn)，共建立了14個(gè)基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型，經(jīng)過(guò)綜合比較分析，從中篩選出擬合效果最優(yōu)的基礎(chǔ)胸徑-樹(shù)高模型（修正Veibull模型）。具體結(jié)論如下：

（1）在15個(gè)基礎(chǔ)胸徑-樹(shù)高模型中，可以求解參數(shù)的模型有14個(gè)，且這14個(gè)模型的P＜0.01，說(shuō)明這14個(gè)模型的參數(shù)值均有顯著統(tǒng)計(jì)學(xué)差異。

（2）通過(guò)對(duì)以上14個(gè)基礎(chǔ)胸徑-樹(shù)高生長(zhǎng)模型進(jìn)行擬合和參數(shù)檢驗(yàn)，分析得出修正Veibull模型擬合效果較好，其模型表達(dá)式為：H=1.3＋5.322（1-exp（-0.149D）），參數(shù)估計(jì)為：a=5.322，b=0.149，c=1.147，適合榆林地區(qū)樟子松基礎(chǔ)胸徑-樹(shù)高模型計(jì)算。

參考文獻(xiàn)

［1］ 趙曉彬，朱建軍，郜超，等. 陜北毛烏素沙地樟子松引種造林綜合分析［J］.防護(hù)林科技，2018（7）：52-54.

［2］ 趙曉彬，劉光哲. 沙地樟子松引種栽培及造林技術(shù)研究綜述［J］.西北林學(xué)院學(xué)報(bào)，2007，22（5）：86-89.

［3］ 劉軒，趙珮杉，高廣磊，等. 沙地樟子松（Pinus sylvestris var.mongolica）物候特征及其對(duì)氣候的響應(yīng)［J］.中國(guó)沙漠，2022，42（2）：25-35.

［4］ 白曉霞，艾海艦. 榆林沙地樟子松人工林土壤養(yǎng)分變化特征［J］.西部林業(yè)科學(xué)，2020，49（3）：80-85.

［5］ 馮偉，李衛(wèi)，楊文斌，等. 毛烏素沙地樟子松固沙林土壤水分對(duì)降雨的動(dòng)態(tài)響應(yīng)［J］.水土保持研究，2019，26（3）：101-107.

［6］ 翟明普，沈國(guó)舫. 森林培育學(xué)［M］.3版.北京：中國(guó)林業(yè)出版社，2016：417.

［7］ 劉春云，方文靜，蔡瓊，等. 中國(guó)落葉松林胸徑-樹(shù)高相關(guān)關(guān)系的探討［J］.北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，53（6）：1081-1088．

［8］ 雷娜慶，鐵牛，劉洋，等. 興安落葉松天然林林分直徑分布和樹(shù)高分布［J］.東北林業(yè)大學(xué)學(xué)報(bào)，2017，45（1）：90-93.

［9］ 張雄清，張建國(guó)，段愛(ài)國(guó)．基于貝葉斯法估計(jì)杉木人工林樹(shù)高生長(zhǎng)模型［J］.林業(yè)科學(xué)，2014，50（3）：69-75.

［10］ ARABATZIS A A，BURKHART H E.An evaluation of sampling methods and model forms for estimating height-diameter relationships in loblolly pine plantations［J］.Forest science，1992，38（1）：192-198.

［11］ ZEIDE B，ZAKRZEWSKI W T.Selection of site trees：The combined method and its application［J］.Canadian journal of forest research，1993，23（6）：1019-1025.

［12］ HUANG S，TITUS S J，WIENS D P.Comparison of nonlinear height-diameter functions for major Alberta tree species［J］.Canadian journal of forest research 1992，22（9）：1297-1304.

［13］ FANG Z，BAILEY R L.Height-diameter models for tropical forest on Hainan Island in southern China［J］.Forest ecology and management，1998，110（1/2/3）：315-327.

［14］ HUANG S.Ecoregion-based individual tree height-diameter models for lodgepole pine in Alberta［J］.Western journal of applied forestry，1999，14（4）：186-193.

［15］ 李善堯. 樟子松天然林樹(shù)高-胸徑模型及胸徑分布規(guī)律分析［J］.林業(yè)科技，2018，43（2）：10-13.

［16］ 羅玲，廖超英. 榆林沙區(qū)不同立地條件下樟子松林木個(gè)體生長(zhǎng)的模型擬合［J］.林業(yè)資源管理，2008（1）：49-52，25.

［17］ 吳夢(mèng)婉. 遼寧章古臺(tái)沙地樟子松人工林樹(shù)木生長(zhǎng)及其對(duì)氣候的響應(yīng)［D］.北京：北京林業(yè)大學(xué)，2019.

［18］ 雷澤勇，韓艷剛，趙國(guó)軍，等. 遼寧章古臺(tái)樟子松生長(zhǎng)過(guò)程分析［J］.浙江農(nóng)林大學(xué)學(xué)報(bào)，2018，35（2）：324-330.

［19］ 盧軍，張會(huì)儒，雷相東，等. 長(zhǎng)白山云冷杉針闊混交林幼樹(shù)樹(shù)高-胸徑模型［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2015，37（11）：10-25.

［20］ 鐘譚衛(wèi).神經(jīng)網(wǎng)絡(luò)中的多維 MP 模型［J］.華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2004，25（2）：121-122.

［21］ 馬天曉，趙曉峰，黃家榮，等. 基于人工神經(jīng)網(wǎng)絡(luò)的樹(shù)高曲線模型研究［J］.河南林業(yè)科技，2006，26（1）：4-5，11.

［22］ 徐步強(qiáng)，張秋良，彌宏卓，等. 基于 BP 神經(jīng)網(wǎng)絡(luò)的油松人工林生長(zhǎng)模型［J］.東北林業(yè)大學(xué)學(xué)報(bào)，2011，39（12）：33-35.

［23］ 金燕. 基于神經(jīng)網(wǎng)絡(luò)和機(jī)器視覺(jué)的南方葡萄專家系統(tǒng)研究［D］.長(zhǎng)沙：湖南農(nóng)業(yè)大學(xué)，2009.

［24］ 楊興，朱大奇，桑慶兵. 專家系統(tǒng)研究現(xiàn)狀與展望［J］.計(jì)算機(jī)應(yīng)用研究，2007，24（5）：4-9.

［25］ 鄧紅兵，郝占慶，王慶禮，等. 紅松單木高生長(zhǎng)模型的研究［J］.生態(tài)學(xué)雜志，1999，18（3）：19-22，31.

［26］ 李明華，肖舜禎. 上海市主要造林樹(shù)種的胸徑-樹(shù)高模型研究［J］.江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2019，41（3）：501-511.

［27］ 婁明華，白超，楊同輝. 寧波石櫟-木荷天然常綠闊葉混交林的樹(shù)高-胸徑模型［J］.林業(yè)與環(huán)境科學(xué)，2021，37（4）：46-54.

［28］ 白天恒. 沙地樟子松人工林胸徑生長(zhǎng)的影響因子研究［D］.阜新：遼寧工程技術(shù)大學(xué)，2019.