跳至正文

L-系统与植物分形:代码长出的树

🔵 数值验证 📅 2026年3月 ⏱ 阅读约12分钟

树木的许多枝条,都像是整棵树的微缩版。蕨叶的每一片小叶,都像是整张叶片的翻版。这种”局部复刻整体”的现象,用一个词可以概括:分形自相似性。但植物真的是按照分形规则生长的吗?如果是,那这套规则究竟写在哪里?

20世纪60年代末,匈牙利生物学家阿里斯蒂德·林登迈耶(Aristid Lindenmayer)给出了一个出人意料的回答——他把植物生长写成了一种字符串重写系统,称为 L-系统(L-system)。从此,树木的分枝、蕨类的卷叶、花序的排列,都可以用几行简洁的文法规则来生成。这不仅仅是数学游戏,而是在说:植物的形态信息,本质上是局部规则的迭代放大

本文将带你进入这个规则的内部,看清 L-系统是如何用字符串推导出三维植物形态的,理解它与分形几何的深层联系,并梳理研究者们如何把这套纯数学框架推进到真实植物的建模与应用之中。

📑 本文目录

一、从细胞分裂到字符串:L-系统的文法逻辑

L-系统最初并不是为了画树。林登迈耶创建它,是为了描述藻类细胞的分裂行为。他的核心直觉是:植物器官的生长,可以看作一种并行字符串重写过程——每个字符代表一个细胞或器官状态,每一步所有字符同时按规则替换,生成下一个发育阶段。

🔑 核心概念:L-系统的三要素
  • 字母表 V:符号集合,每个符号代表一种细胞/器官类型
  • 公理 ω:起始字符串,相当于种子状态
  • 产生式规则 P:每个符号在下一步如何展开的替换规则

最经典的例子是藻类细胞的分裂模型。设字母表 V = {A, B},公理 ω = A,产生式规则为:

A → AB
B → A

翻译成人话:A 类细胞每代分裂成一个 A 和一个 B;B 类细胞每代成熟为 A 类。从单个 A 开始迭代,得到字符串序列:A → AB → ABA → ABAAB → ABAABABA……每一代的字符串长度,正好是斐波那契数列。

这不是巧合。植物界的斐波那契模式(叶序角137.5°、螺旋数的相邻斐波那契比)正是类似迭代规则在二维空间展开的结果。L-系统第一次把这种联系写成了可计算的形式。

二、乌龟作图:把字符串变成三维枝条

字符串本身是一串符号,要变成看得见的植物形态,需要一套几何解释机制。研究者引入了计算机图形学中的经典概念:乌龟几何学(Turtle Geometry)。

📐 乌龟状态

想象一只乌龟在三维空间中爬行,它的状态由位置向量和朝向决定:

(x, y, z, H, L, U, δ)
符号含义
x, y, z当前位置坐标
H朝向向量(Heading)
L左方向向量(Left)
U上方向向量(Up)
δ每步转动角度

翻译成人话:乌龟的”爬行日志”就是字符串。遇到 F 就向前走一步画线,遇到 + 就左转δ度,遇到 – 就右转δ度,遇到 [ 就把当前状态压栈(相当于”记住这个岔路口”),遇到 ] 就弹栈回到上一个岔路口。这样,带括号的字符串就自然生成了”先走主干,再回去长侧枝”的树形结构。

一个生成二叉树的简单产生式如下:

公理: F
规则: F → F[+F]F[-F]F

翻译成人话:每一段枝条 F,在下一步展开为:继续延伸,同时向左发出一根侧枝,继续延伸,再向右发出一根侧枝,最后还有主干延伸。迭代几步之后,一棵层级分明的分形小树就出现了——每一层的分枝结构,都是上一层的微缩重复。

药用植物三维虚拟建模的研究验证了这套机制:研究者基于准二叉树(quasi-binary tree)分析植物三维分枝,用 L-系统实现器官和枝条结构的程序化生成,最终得到与真实植物外观高度吻合的三维模型[1]

三、自相似的数学结构:L-系统与分形维度

L-系统生成的形态为什么看起来像分形?这不只是视觉错觉,而有数学依据。

📐 分形维度与自相似

对于一个自相似结构,若将其缩小到原来的 1/r,可以拼出 N 个副本,则其相似维度 D 满足:

D = log(N) / log(r)

翻译成人话:整数维度的物体很直觉——线是1维,面是2维,体是3维。但分形介于整数维之间。一棵自然生长的树,其分枝结构既不完全是二维的(平面)也不完全是三维的(填满空间),其分形维度通常介于整数维度之间。L-系统生成的植物形态,正是通过迭代规则精确控制了这个”非整数维度”。

L-系统与植物分形生成融合研究中,研究者将 L-系统与迭代函数系统(IFS)结合,精确生成具有预设分形维度的树木结构,验证了 L-系统产生的层级自相似性可以用分形几何精确刻画[2][3]

💡 类比理解:俄罗斯套娃里藏着什么

L-系统就像一套俄罗斯套娃的制作规则:你只需要知道”大娃娃里放什么样的小娃娃”,剩下的由规则自动嵌套。每迭代一次,就多嵌套一层。植物的分枝层级,正是大自然用自己的”L-系统规则”嵌套出来的。区别在于:套娃的层数有限,而理想分形可以无限迭代。

四、从确定性到随机性:L-系统的变体谱系

经典 L-系统的产生式规则是确定性的——同一个符号每次展开结果完全相同。但真实植物不是这样的:同一棵树上,每根侧枝的角度、长度都有细微差异。为了更贴近现实,研究者发展出了多种 L-系统变体。

L-系统变体谱系
类型特点适用场景
D0L-系统确定性,无上下文理论分析、规律形态
随机L-系统(stochastic)规则按概率选择自然树木的随机变异
参数化L-系统符号携带数值参数枝长、角度精细控制
上下文敏感L-系统规则依赖邻居状态信号传导、极性建立
开放L-系统与外部环境交互光照、重力响应建模

谷类作物分枝建模综述指出,要真实刻画分蘖动态,需要在 L-system 框架中融合概率方法(stochastic rules)、剂量-反应方法和机制方法,三者各自对应确定性规则、参数化规则和环境响应规则的不同组合[4]

参数化 L-系统的关键在于把数值参数嵌入符号,让规则能”感知”器官的当前状态。例如,一个带参数的展开规则可写为:

F(l, w) : l > 0 → F(l·r, w·q) [+(δ)F(l·e, w·q)] [-(δ)F(l·e, w·q)]
参数含义
l当前节间长度
w当前节间直径
r主干延伸比
q直径缩小比
e侧枝相对于主干的长度比
δ分枝角

翻译成人话:这条规则在说——每一段枝条,在下一步会按比例r延伸自己的主干,同时向左和向右各发出一根按比例e收缩的侧枝,所有枝条的直径都按比例q变细。只需调整这几个参数,就能生成从铅笔松到垂柳的各种树形。

L-系统植物结构建模改进研究正是沿着这条路线推进,通过增强参数化规则与随机性,显著提升了从简单规则到复杂真实植物形态的拟合能力[5]

五、根分生组织的细胞谱系:L-系统进入生物实验

L-系统不只是纸面上的数学游戏。在植物发育生物学中,它被用来重建真实细胞的谱系路径——也就是追踪”这个细胞的祖先是谁,后代会是谁”。

对玉米初生根分生组织的研究将细胞分裂序列与”bootstrap” L-系统结合,用形式化规则压缩原本极其庞杂的可能分裂路径空间。结果出人意料:在理论上可能的海量分裂路径中,真实发育只采用了极少数受限路径[6]

🔬 实验发现:细胞分裂的路径约束

玉米根分生组织中,每类细胞包群的分裂模式可以用一套约束规则完整编码。这意味着植物分生组织的发育,并不是随机探索所有可能路径,而是被局部规则严格约束在少数”允许路径”上。L-系统的价值在于,它为这种约束提供了可计算的形式语言[6]

这一发现意义深远:它说明植物形态的”自相似性”不只是几何现象,在细胞层面也存在规则驱动的重复性。对植物器官轴分枝与细胞后代关系的理论分析进一步表明,确定性的细胞谱系模式约束了器官轴的分枝方式,这为”分枝规则可形式化编码”提供了细胞学基础[7]

把细胞测量数据转化为 L-系统规则,是实验验证的另一条路径。水稻形态发育研究用三维数字化仪连续测量水稻从苗期到成熟期的结构变化,再将实测数据转译为可生成的植物建筑规则,包括分蘖、叶片出现与器官几何参数,构建出”3D virtual rice”[8]

六、功能-结构植物模型:让L-系统生理化

早期 L-系统只关心几何形态——枝条怎么长、叶片怎么排。但真实植物的形态,是由生理过程驱动的:光合产物从叶片运往根,激素信号从顶端传向侧枝,水分沿木质部向上流动。形态与功能必须耦合,才能真正理解植物。

这催生了功能-结构植物模型(Functional-Structural Plant Model,FSPM)——L-系统作为结构骨架,生理方程作为驱动引擎,两者深度整合。

📐 源-汇关系的数学耦合

在 FSPM 框架中,碳分配可以写为一个优化问题。设器官 i 的碳需求为 D_i,整株碳供给为 S,则分配给器官 i 的碳量 C_i 满足:

C_i = S · (D_i / Σ_j D_j)

这个方程嵌入 L-系统每一步迭代,使得新生枝条的生长速率由当前碳收支决定,而非固定规则。

翻译成人话:这不再是”按规则画枝条”,而是”按碳账本决定要不要长枝条”。整棵树有多少光合产物,每个器官需要多少,按比例分配——这才是植物真正的生长逻辑。L-系统提供了枝条的几何框架,生理方程决定了每根枝条能否获得足够的资源去生长。

棉花分枝对源-汇关系影响的研究将不同分枝拓扑嵌入虚拟植物框架,比较结构变化对生理资源分配的作用,说明植物分枝结构既是几何问题,也是功能问题[9]。大豆根系结瘤自调控的建模研究则将 L-system 思路扩展到地下根系,对带侧根与结瘤的三维根系建筑进行建模,协调不同时间尺度上的发育与信号过程[10]

蚕豆冠层 ALAMEDA 模型明确采用 Lindenmayer systems 作为结构模型的形式化语言,同时结合形态学与物候学参数实现动态可视化和数值输出,通过详细实测参数化后完成验证,展示了 L-system 如何支持从单株器官排列到群体冠层结构的建模[11]

七、从猕猴桃藤到芒果树:真实树种的L-系统建模

当 L-系统与 FSPM 框架成熟之后,研究者开始将其应用于具体经济树种,验证模型在真实农业场景中的预测能力。

🌍 案例一:猕猴桃藤本的整合建模

猕猴桃功能-结构模型在 L-system 平台上将枝条建筑、碳动态、环境效应与器官位置效应整合。这意味着 L-系统从”纯几何生成器”升级为”承载生理过程的建模框架”——同一套字符串重写规则,不仅控制枝条生长,还同步计算光合、呼吸和碳分配[12]

🌍 案例二:辐射松的木材质量预测

辐射松结构-功能模型用 L-system 语言定义操作单元及其形态发生、光合作用和呼吸规则,不仅重建树体分枝建筑,还连接木材性质与经营情景。这说明 L-系统框架可以跨越”形态”与”材质”两个尺度,为林业决策提供数量化依据[13]

🌍 案例三:芒果树的果实产量模拟

V-Mango 是面向芒果树的整合型结构-功能模型,精细描述营养与生殖发育过程,可模拟不同栽培措施下的产量响应。这代表了 L-system/FSPM 路线在果树方向的成熟应用——复杂分枝结构通过规则化建筑模型得到再现,并具备预测能力[14]

果树作物模型的更新综述回顾了这一领域的整体进展,明确总结了 L-system、pipe model 和 Markov 模型在分枝、芽命运与节间伸长模拟中的分工:L-system 负责生成结构骨架,pipe model 约束导管横截面积与叶面积比,Markov 模型处理芽的命运概率[15]。被子植物树体结构的分子遗传与激素基础综述则从生物学角度解释了为何植物会表现出重复、自相似、可规则化的分枝模式,为算法生成提供了机制背景[16]

八、与IFS融合:分形生成技术的协奏

L-系统并非唯一的分形生成工具。在计算机图形学领域,另一类强大的方法是迭代函数系统(Iterated Function System,IFS)。两者的逻辑不同,但可以互补。

📐 IFS的核心方程

IFS 通过一组仿射变换的迭代来生成分形。设有 n 个收缩映射 {w₁, w₂, …, wₙ},每个 wᵢ 形如:

wᵢ(x) = Aᵢx + bᵢ
符号含义
Aᵢ线性变换矩阵(旋转、缩放、反射)
bᵢ平移向量
x当前点坐标

反复对任意初始点应用这组映射,最终吸引到的极限集合,就是 IFS 的吸引子——即分形图案本身。

翻译成人话:IFS 就像是一套”折纸规则”——把一张纸按特定方式缩小、旋转、移动,重叠放置,反复操作下去,最终”折”出一个分形图案。著名的巴恩斯利蕨就是用4个仿射变换的IFS生成的,视觉上几乎与真实蕨叶无法区分。

L-系统负责基于文法的生长规则,IFS 负责几何自相似变换——两者结合可以生成单独使用任何一种方法都难以达到的复杂层级结构。三维分形模拟研究直接展示了这种融合:L-system 定义全局分枝拓扑,IFS 在局部补充自相似细节,两套机制在不同尺度上分工协作[2]。植物仿真融合研究进一步将这种方法推广到树木、枝条和叶片的层级细节生成,获得了更逼真的分枝结构与视觉自相似性[3]

植物生长可视化研究则专注于如何把形式文法规则转化为连续的三维生长动画与几何结构,解决了”规则正确但视觉不自然”的工程难题[17]

九、基因组遇见L-系统:植物建筑的未来

L-系统的参数(分枝角、节间长、生长速率)从哪里来?在传统模型中,这些参数由研究者手动测量或拟合。但植物的建筑形态是由基因决定的——能否把基因组信息直接”注入”L-系统的参数?

这个方向已经有了实质性进展。苹果树结构-功能模型 MAppleT 与全基因组信息相结合,模拟一年生苹果树的建筑发育,将器官发生、形态发生和遗传决定因子联结起来——植物分枝结构不仅可由规则生成,还能与可测遗传参数耦合[18]

🚀 前沿方向:从基因组到虚拟植株

如果每个 L-系统参数都对应可测的遗传变异,理论上可以从个体基因型直接预测其三维植株建筑。这将彻底改变育种流程:不再需要种出来才知道长什么样,而是在计算机中先”生长”出来。数值植物模型与表型相关空间的结合研究探索了这一路线,把植物性状修改投影到表型相关空间中,评估理想性状组合是否在发育上可达[19]。遗传算法与功能-结构植物模型结合的水稻研究则以光合作用输出为目标函数,动态优化虚拟水稻的叶角、株高、分蘖数等表型参数,探索植物建筑与性能之间的设计空间[20]

这条路线的核心逻辑是:L-系统是连接基因型空间表型空间的桥梁。基因决定规则参数,规则迭代产生形态,形态影响生理表现,生理表现反过来约束基因的选择压力。把这个闭环计算化,就是下一代植物科学的核心工具之一。

📜 历史回顾:从藻类到数字森林
  • 1968年:林登迈耶提出 L-系统,用于描述藻类细胞分裂
  • 1970-80年代:引入乌龟几何学,L-系统生成第一批植物分形图像
  • 1990年代:参数化、随机化、上下文敏感变体相继成熟
  • 2000年代:与功能生理模型整合,FSPM框架形成
  • 2010年代:进入真实树种建模,与基因组信息耦合
  • 2020年代:迈向”从基因到三维植株”的直接预测

🎯 关键要点
  • L-系统用字符串并行重写规则模拟植物生长,本质是局部规则的迭代放大,天然产生分形自相似结构
  • 乌龟几何学将抽象字符串转化为三维枝条,”[“和”]”实现分枝回溯,是L-系统生成树形结构的关键机制
  • L-系统生成的结构具有可计算的分形维度,与IFS融合可在不同尺度上协同生成复杂植物形态
  • 玉米根分生组织实验显示,真实发育路径被严格约束在少数允许路径上,验证了L-系统的生物学基础
  • 功能-结构植物模型将L-系统与生理方程耦合,使形态生成由碳收支、光照和激素信号驱动,而非固定规则
  • 猕猴桃、辐射松、芒果等真实树种的成功建模,证明L-系统框架已超越理论,具备农业与林业应用价值
  • 基因组信息与L-系统参数的耦合,正在开启”从基因型预测三维植株建筑”的新时代

📚 参考文献

  1. Ding D et al. Virtual medical plant modeling based on L-system. African Health Sciences. 2014. PMID: 25834517
  2. Deng Fang et al. An application of L-system and IFS in 3D fractal simulation. WSEAS Transactions on Systems. 2008.
  3. Jinshu Han et al. Plant Simulation Based on Fusion of L-System and IFS. 2007.
  4. Evers J et al. Modeling branching in cereals. Frontiers in Plant Science. 2013. PMID: 24133499
  5. Zhang Shu et al. Improvement of Plant Structure Modeling Based on L-System. Journal of Image and Graphics. 2002.
  6. Lück J et al. Cell genealogies in a plant meristem deduced with the aid of a ‘bootstrap’ L-system. Cell Proliferation. 1994. PMID: 10465023
  7. Barlow P et al. Deterministic cellular descendance and its relationship to the branching of plant organ axes. Protoplasma. 2004. PMID: 15614474
  8. Watanabe T et al. Rice morphogenesis and plant architecture: measurement, specification and the reconstruction of structural development by 3D architectural modelling. Annals of Botany. 2005. PMID: 15820987
  9. Dong Li et al. Modeling branching effects on source-sink relationships of the cotton plant. arXiv: 1012.3276. 2010.
  10. Han L et al. A functional-structural modelling approach to autoregulation of nodulation. Annals of Botany. 2011. PMID: 20826439
  11. Ruiz-Ramos M et al. ALAMEDA, a structural-functional model for faba bean crops: morphological parameterization and verification. Annals of Botany. 2006. PMID: 16390842
  12. Cieslak M et al. A functional-structural kiwifruit vine model integrating architecture, carbon dynamics and effects of the environment. Annals of Botany. 2011. PMID: 20855486
  13. Fernández M et al. A functional-structural model for radiata pine (Pinus radiata) focusing on tree architecture and wood quality. Annals of Botany. 2011. PMID: 21987452
  14. Boudon F et al. V-Mango: a functional-structural model of mango tree growth, development and fruit production. Annals of Botany. 2020. PMID: 32391865
  15. Grisafi F et al. Fruit tree crop models: an update. Tree Physiology. 2022. PMID: 34542149
  16. Hollender C et al. Molecular basis of angiosperm tree architecture. New Phytologist. 2015. PMID: 25483362
  17. Yongxiang Sun et al. Visualization of plant growth based on l-system. 2010 World Automation Congress. 2010.
  18. Migault V et al. Combining Genome-Wide Information with a Functional Structural Plant Model to Simulate 1-Year-Old Apple Tree Architecture. Frontiers in Plant Science. 2016. PMID: 28127302
  19. Victor Picheny et al. Using numerical plant models and phenotypic correlation space to design achievable ideotypes. arXiv: 1603.03238. 2016.
  20. Ding W et al. Genetic algorithm based approach to optimize phenotypical traits of virtual rice. Journal of Theoretical Biology. 2016. PMID: 27179460