紅楓分形維數中的物理機制

顏揚杰，韋建衛(wèi)，程榮祿，夏雪峰

(重慶理工大學，重慶 400054)

紅楓分形維數中的物理機制

顏揚杰，韋建衛(wèi)，程榮祿，夏雪峰

(重慶理工大學，重慶 400054)

采用Matlab程序計算與實驗測量相結合的方法，以重慶理工大學校園內隨機選取的3棵紅楓為研究對象，分別測量并計算了植物整體和部分的三維空間分形維度，得出該種植物分形維度與陽光等物理因素之間的關系，以及植物生長過程中的維度變化規(guī)律。研究結果表明:植物的維度大小與其生長的物理環(huán)境具有邏輯關系。這些關系有助于深入研究高等植物的生長規(guī)律，進而為人工裁剪和控制的自然原理提供理論支持。

紅楓;分形維度;計盒維數

傳統(tǒng)的幾何學不能夠描述大自然的復雜性和多樣性，從而使我們無法了解這些復雜性和多樣性背后的機制。美國科學怪杰Mandelbrot于1975年首先提出分形幾何的概念。主要用于描寫云彩﹑山嶺﹑海岸線或樹木等自然而不規(guī)則的形狀。而分形幾何學在描述和理解大自然的復雜性和多樣性方面具有重要作用［1－3］。十幾年來，分形迅速發(fā)展，應用范圍涉及自然科學的各個領域。在分形幾何理論中，分形維數是最重要基本概念之一。關于分形維數的定義，主要包括Hausdorff維數、信息維數、關聯維數、計盒維數、相似維數等。高等植物為構件生物，是由構件不斷地生長增加而形成的，如由葉、芽、枝、花、根的反復不斷地生長增加而形成整個植物體。高等植物的形態(tài)結構具有自相似特征，摘下其中任意一支枝干，將其插在地上，可看成縮小的植株。但并非植株的每一部分都生長得同樣茂盛，這與它所在的生長環(huán)境(如陽光、空間分布等物理因素)有密切關系。近年來，運用分形理論對植物的各種具體特征的研究越來越多［4－7］。

分形早期用來描述一些復雜但具有嚴格數學規(guī)則的特殊結構，如柯赫曲線(Koch curve)、柯赫雪花(Koch snowflake)、康托塵埃(Cantor dust)、明可夫基香腸(Minkowski sausage)、謝爾賓斯基方毯(Sierpinski carpet)等等。這些圖形具有嚴格的自相似性，將圖形中任意部分放大都能得到與整體相似的形狀，但是這些結構因其無窮小、無限多的迭代而無法嚴格存在于自然界中。同樣的，自然界中的物體也難以進行嚴格明確的數學定義。英國數學家Falconer［8］參考生物學家的做法，通過列出分形的具體特性來給分形下定義。

本文運用分形幾何的理論對重慶理工大學校園常見的一種植物的生長特性進行了研究，通過計算植物不同部分的計盒分形維度，分析得到的分形維度與植物所在的空間環(huán)境關系，為研究高等植物的生長規(guī)律提供一定的理論支持。

1 分形、分形維數與計盒維數原理

在經典的歐幾里得幾何中，我們已經習慣于用拓撲維度來表示空間:如零維的點、一維的線、二維的平面、三維的立方形，其所表示的維度都是整數。然而，大自然中的山峰、樹木、云朵等物體因其復雜和不規(guī)則性并不能用簡單的線段、矩形、三角形和圓等規(guī)則的圖形來描述，分形因此誕生。它能很好地描述這些復雜的物體，所得到的維度并非整數，而是分數，稱為分形維數。例如計算出柯赫曲線的計盒維數為1.261 8，它大于線的拓撲維度1。原因是柯赫曲線可無限進行下去，將幾乎覆蓋整個平面。

試想有一條單位長度的線段，用長為r的小線段去覆蓋它，需要1/r個小正方體;如果是單位邊長的正方形，則需要1/r2個小正方形;如果是單位邊長的立方體，覆蓋它需要1/r3個小正方體。值得注意的是r的指數和被覆蓋的對象的維度是相同的，這并非偶然。一般我們可以用“盒子”來覆蓋任意形狀的圖形。

設A是Rn空間的任意非空有界子集，對任意的一個r＞0，Nr(A)表示用來覆蓋A所需的邊長為r的n維立方體(盒子)的最小數目。如果存在一個數d，使得當r→0時，有

那么稱d為A的計盒維數(簡稱盒維數)。盒維數為d，當且僅當存在一個正數k使得

由于方程兩邊都為正，因此可以對方程兩邊取對數，得

這里舍去了logk這一項，因為它是常數項。當r→0時，分母趨于無窮大。

在實際計算中，可以使用不同邊長的立方體來覆蓋圖形，得出不同r值下的Nr(A)，再以logr為橫坐標、logNr(A)為縱坐標的雙對數坐標系中描出點(logri，logNri(A))，最后由這些點擬合出來的直線的斜率的負值便是圖形A的維數。

2 樣本植株

圖1為測量維度的紅楓樣本植株。圖1(a)～(c)分別為樣本1、樣本2、樣本3。

圖1 測量維度的紅楓樣本植株

3 測量方法及所得數據

首先在空間中建立三維直角坐標，選取適當的原點，本次測量都將原點取在離植株較遠的位置，確保植株所有枝干都能在坐標系的第一象限中。定好原點，確定好主枝干的坐標位置，便可開始測量每一級分叉點的空間坐標，即分叉點在xz平面和y－z平面投影坐標，并記下是第幾級分叉點。一人測量，一人記錄。測量工具為米尺，其精度為 mm，本次研究的精度為 cm，且數據取整數。

4 結果及分析

根據測量所得數據，已知2個分叉點的坐標便可在空間確定一條直線，該直線即為樣本枝干的模擬。利用Matlab［9］軟件建立一個6行n列的矩陣，n為所測量的樣本的枝干數目，前3行為一個分叉點的坐標，記為x1，y1，z1;后3行為另一個分叉點的坐標，記為x2，y2，z2。接著在z1到z2之間每隔0.01取一個數，在x－z平面利用直線的點斜式公式便可算得該線段上的其他點的x坐標，同理，在y－z平面得到該線段上其他點的y坐標。這樣一系列的點在空間構成一條線段，重復以上步驟得到三維立體數字樣本模型(見圖2)，利用前面所得的模型可計算其盒維數。首先將空間分割成一個個“盒子”，并為每個“盒子”按其空間位置編號，如(2 3 4)表示“盒子”的x坐標在r到2r之間;y坐標在2r到3r之間;z坐標在3r到4r之間，其中r為“盒子”邊長。每條線段上的點依次判斷在哪個“盒子”里，最后去掉編號相同的“盒子”，所得的便是該樣本模型在空間中占據的“盒子”數，記為Nr。改變“盒子”邊長重復計算得到在不同邊長時占據的“盒子”數。最后在雙對數坐標下描出對應的點，并運用最小二值法擬合直線(見圖3)，所得直線斜率的負值即為盒維度。

圖2 三維立體數字樣本模型

圖3為3個樣本數據在雙對數坐標中擬合的直線。擬合圖縱坐標為盒子數的對數，橫坐標為盒子邊長對數。data1，data2對應樣本1;data3，da－ta4對應樣本2;data5，data6對應樣本3。

圖3 3個樣本數據在雙對數坐標中擬合的直線

表1 3個樣本的維度結果

表2 3個樣本在不同尺寸下的盒子覆蓋數

由于同種植物應有相同的生長規(guī)律，所得的整體維數應該接近。計算結果表明:3棵紅楓的整體維數接近，符合實際情況。

在向陽的方向得到的維數都略大于背陽的方向。從得到的圖像中也可以看出:向陽方向的圖像更為復雜，在實際中表明植物向陽方向將生長得更為茂密，植物有向光方向生長的趨勢。向陽部分的維度結果表明:陽光能夠提供更多的生長動力，使其更多地進行分叉、抽枝，以占領更多空間區(qū)域而獲得更多能量，進而促進整個植株的整體優(yōu)勢地位。

植物的樹冠部分與靠近地面的樹干部分的維數有較大的差異，在實際中表明植物越往上生長，其結構會更為復雜。在底部主枝干分級少，往上生長分級越多，各級枝干越多，生長茂密，橫向寬度也增大。該植物樣本早期時維度較小，隨著生長發(fā)育的進行維度逐漸增加，當成年后基本確定下來。因此，想要使得植物生長較為高大，則需要不斷地修剪多余的分叉，使得植株在維度較低情況下繼續(xù)分叉生長。

計算結果如表1和2所示，3棵紅楓植株的整體分形維數分別為1.395 1，1.482 7，1.308 0。本次研究所得的圖像是在三維空間中用描點的形式模擬實際中紅楓植株的形狀，忽略了植株枝干的粗細，所得圖像在三維空間中占少量的體積，得到的維度比實際數據略有偏低，但同時也考慮了樣本1入直徑數據的計算情況。結果發(fā)現:各個部分的維度都有增加，但是它們之間的相對關系并無根本性的變化，因此并不影響結論，故認為此次維數計算是合理的。

5 結束語

本次研究對植株分叉點的坐標采取人工測量，測量過程對樹葉的空間占據有所忽略，但據實物情況來看，葉子數基本正比于枝干的長度，枝干的維度也就基本反映了葉子的各部分情況，所以定性結果是較為準確的。若能采用大型三維投影設備可得到植株的全部信息，增加樣本數量，再對程序進行改進，加上植株的實際直徑，得到的結果將更科學、更有說服力。

［1］ 陳紅，夏青，左婷，等.基于紋理分析的香菇品質分選方法［J］.農業(yè)工程學報，2014，30(3):281－287.

［2］ 單立山，李毅，任偉，等.河西走廊中部兩種荒漠植物根系構型特征［J］.應用生態(tài)學報，2013，24(1):25－31.

［3］ 李志臣，姬長英.基于圖像分析的雜草分形維數計算［J］.農業(yè)工程報，2006，11(11):175－178.

［4］ 朱華，姬翠翠.分形理論及其應用［M］.北京:科學出版社，2011.

［5］ 溫芝元，曹樂平.椪柑果實病蟲害的傅里葉頻譜重分形圖像識別［J］.農業(yè)工程學報，2013，29(23):159－164.

［6］ 丁雪嬌，韓紅，龐彩菊，等.砂砧薹草和假牛鞭草種群空間分布格局的分形特征［J］.浙江農林大學學報，2013，30(2):220－225.

［7］ 袁津生，姚宇飛.基于分形維度的葉片圖像識別方法［J］.計算機工程與設計，2012，33(2):670－673.

［8］ 張濟忠.分形［M］.北京:清華出版社，2011.

［9］ 張波，沈火明.基于Matlab的殼體有限元分析［J］.重慶理工大學學報:自然科學版，2010，24(12):77－81.

(責任編輯 劉 舸)

Physical Mechanism of Fractal Dimensionality of Red Maple

YAN Yang－jie，WEI Jian－wei，CHENG Rong－lu，XIA Xue－feng
(Chongqing University of Technology，Chongqing 400054，China)

Combining experimental measurement and the calculation of Matlab code，the physical factors of plant growth have been analyzed through fractal dimensionality data.Three red maples were selected randomly and measured.The space coordinates of the branches were recorded and loaded into computer.The fractal dimensionalities were calculated and analyzed with Matlab code through boxcounting method.We find the relevance between the sunlight and the fractal dimensionalities.The law of the dimensionality changing of the plant has also been studied.The results indicate that there is an inevitable logical relationship between physical environment and dimensionality of the plant.The relationship will assist the investigation of plant growth rhythm.Furthermore，it will benefit to the natural principle of artificial clip.

red maple;fractal dimensionality;box－counting dimension

O29

A

1674－8425(2014)07－0127－05

10.3969/j.issn.1674－8425(z).2014.07.025

2014－02－20

國家自然科學青年基金資助項目(11204391);重慶市教委項目(KJ130831);2013年校級科研立項項目(GD15)

顏揚杰(1992—)，福建泉州人，主要從事計算機仿真模擬數值計算研究;通訊作者韋建衛(wèi)(1980—)，男，河南焦作人，博士，副教授，主要從事納米材料與器件方面的研究。

顏揚杰，韋建衛(wèi)，程榮祿，等.紅楓分形維數中的物理機制［J］.重慶理工大學學報:自然科學版，2014(7): 127－131.

format:YAN Yang－jie，WEI Jian－wei，CHENG Rong－lu，et al.Physical Mechanism of Fractal Dimensionality of Red Maple［J］.Journal of Chongqing University of Technology:Natural Science，2014(7):127－131.