842
个编辑
更改
跳到导航
跳到搜索
第1行：
第1行： − + − − + − − |- − | colspan="1" rowspan="1" |'''中文名称''' − | colspan="1" rowspan="1" |菲茨休-南云模型 − |- − | colspan="1" rowspan="1" |'''英文名称''' − | colspan="1" rowspan="1" |FitzHugh–Nagumo Model（FHN Model） − + − − + − − + − − + − − + − − + − − + − − HH 模型虽精确，却因变量多、参数繁而难以做几何与解析分析。菲茨休受范德波尔（van der Pol）弛豫振荡器启发，于 1961 年在《生物物理杂志》（''Biophysical Journal''）提出仅含两个变量的模型，证明可兴奋性的核心特征并不依赖具体离子通道。1962 年，南云仁一等用隧道二极管搭建出具有相同动力学的等效电路，使这一抽象模型获得物理实现，模型由此得名。它也常被称为邦赫费尔-菲茨休模型，以纪念更早提出类似可兴奋系统思想的德国物理化学家卡尔-弗里德里希·邦赫费尔。+ + + − FHN 模型把霍奇金-赫胥黎模型的四维动力学压缩到两个变量：快变量 v 代表膜的兴奋程度（对应[[膜电位]]与钠激活的综合效应），慢变量 w 代表恢复过程（对应[[钠失活]]与[[钾激活]]）。快变量沿一条三次型零倾线运动，具有静息、兴奋平台两个稳定支和中间一个不稳定支；慢变量则沿一条线性零倾线缓慢调节。两者在相平面上的相互追赶，使系统在阈下刺激时回到静息态、在阈上刺激时绕相平面做一次大幅"远足"（对应一个锋电位），或在参数合适时进入周期性弛豫振荡。+ + + + + + + − 一种经过无量纲化、并相对原始表述把恢复变量符号反号的写法为：+ − + − $$\dfrac{dw}{dt}=\varepsilon\,(v-a-bw$$+ + − 其中 v、w 分别为无量纲的兴奋与恢复坐标，I(t) 为重新标度后的刺激；a、b 决定恢复零倾线的位置，$$0<\varepsilon\ll$$ 表示两个变量之间强烈的时间尺度分离。令右端为零可得两条零倾线：v 零倾线为三次曲线 $$w=v-v^{3}/3+$$，w 零倾线为直线 $$w=(v-a)/$$。需要说明，菲茨休 1961 年的经典原始约定中恢复项写作 $$\varepsilon(v+a-bw$$，与上式仅相差恢复变量的符号约定，二者相平面拓扑完全等价。+ − 在相平面上，三次零倾线的左支对应静息（恢复）状态、右支对应兴奋平台、中间支不稳定。当两条零倾线的交点落在中间支上时，静息态失稳，系统被推向围绕不稳定点的稳定[[极限环]]，产生周期锋电位；交点落在外侧稳定支时，系统在小扰动后回到唯一稳定[[不动点]]。下图给出阈下与阈上两种刺激的对比，以及相平面上的零倾线、向量场与极限环。+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + − + − FHN 相平面（自绘，采用经典原始符号约定）：三次 v 零倾线与线性 w 零倾线、向量场及稳定极限环；轨迹沿左支恢复、在膝点跳转至右支兴奋平台。+ + − [[文件:时间响应.png|无框|772x772像素]]+ + + + + − FHN 时间响应（自绘数值仿真）：左为阈下刺激后回到稳定静息态，右为阈上刺激下的弛豫振荡，快变量 v 形成周期性锋电位+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + − * '''可兴奋性的几何解释'''：阈值、全或无响应和兴奋后的缓慢恢复，都可由零倾线的形状与快慢时间尺度直接读出。+ − * '''弛豫振荡'''：快慢两变量交替主导，形成波形不对称、周期主要由慢过程决定的振荡。 − * '''分岔分析'''：随刺激 I 改变，不动点可经[[霍普夫分岔]]失稳并产生极限环，也可出现双稳态与[[迟滞]]。 − * '''波形传播'''：把单胞模型空间耦合，可得到[[可兴奋介质]]中的[[行波]]与[[螺旋波]]，被广泛用于神经组织、心肌乃至化学反应（[[别洛乌索夫-扎博京斯基反应]]）。 − + − '''优势'''：仅两个变量即可完整呈现兴奋性与振荡转换，相平面与分岔分析直观可行，计算量极小，非常适合定性研究、教学和大规模耦合网络的理论分析。+ − '''局限'''：变量是无量纲的抽象坐标，不对应可直接测量的膜电位或[[门控概率]]；模型不保留钠、钾通道的特异性解释，锋电位波形、离子电流和[[不应期]]的定量细节无法从中直接获得，与实验数据对接时需要重新标定。+ − == 与其他模型的关系 ==+ − + − # FitzHugh, R. (1961). Impulses and Physiological States in Theoretical Models of Nerve Membrane. ''Biophysical Journal'', 1(6), 445–466. DOI: 10.1016/S0006-3495(61)86902-6+ − # Nagumo, J., Arimoto, S., & Yoshizawa, S. (1962). An Active Pulse Transmission Line Simulating Nerve Axon. ''Proceedings of the IRE'', 50(10), 2061–2070. − ''说明：本文相平面配图采用菲茨休 1961 年经典原始符号约定（恢复项为 v+a−bw），与正文给出的反号改写形式拓扑等价；配图均为自绘示意性数值仿真。''+
无编辑摘要
'''菲茨休–南云模型'''（FitzHugh–Nagumo model，简称 '''FHN 模型'''）是描述可兴奋系统的二维非线性模型。它用一个快速的兴奋变量和一个较慢的恢复变量，刻画神经膜的兴奋、恢复及周期发放。作为[[霍奇金-赫胥黎模型]]的定性简化，FHN 模型保留了部分关键动力学机制，但不逐一描述钠、钾等离子通道的电流。<ref name="FitzHugh_1961">{{cite journal|author1=Richard FitzHugh|title=Impulses and Physiological States in Theoretical Models of Nerve Membrane|journal=Biophysical Journal|year=1961|volume=1|issue=6|pages=445–466|doi=10.1016/S0006-3495(61)86902-6|url=https://doi.org/10.1016/S0006-3495(61)86902-6}}</ref>
'''[[菲茨休-南云模型]]（FitzHugh–Nagumo Model，简称 FHN 模型）'''是[[霍奇金-赫胥黎模型]]的二维简化模型，用一个快速的"[[兴奋变量]]"和一个缓慢的"[[恢复变量]]"刻画神经膜的阈值、[[全或无式响应]]与[[弛豫振荡]]。它由美国生物物理学家理查德·菲茨休（Richard FitzHugh）于 1961 年提出，后由南云仁一（Jin-ichi Nagumo）等人于 1962 年用电子电路（[[隧道二极管]]）实现，因而合称菲茨休-南云模型。该模型不追踪具体[[离子通道]]，却以[[相平面]]几何清晰地揭示了[[可兴奋性]]的本质，是研究兴奋、振荡与[[分岔]]的经典低维范式。
== 基本信息 ==
== 基本信息 ==
{| class="wikitable"
{| class="wikitable"
! colspan="1" rowspan="1" |项目
! 项目 !! 内容
! colspan="1" rowspan="1" |内容
|-
|-
| colspan="1" rowspan="1" |'''又称'''
| 中文名称 || 菲茨休–南云模型
| colspan="1" rowspan="1" |[[邦赫费尔-菲茨休模型]]（Bonhoeffer–FitzHugh）
|-
|-
| colspan="1" rowspan="1" |'''建模尺度'''
| 英文名称 || FitzHugh–Nagumo model（FHN model）
| colspan="1" rowspan="1" |单细胞（低维抽象）
|-
|-
| colspan="1" rowspan="1" |'''数学形式'''
| 数学形式 || 二维常[[微分方程]]组，可推广为反应–扩散方程或耦合网络
| colspan="1" rowspan="1" |二维常[[微分方程]]组（[[快慢系统]]）
|-
|-
| colspan="1" rowspan="1" |'''状态变量'''
| 状态变量 || 兴奋变量 <math>v</math>、恢复变量 <math>w</math>
| colspan="1" rowspan="1" |兴奋变量 v、恢复变量 w（均无量纲）
|-
|-
| colspan="1" rowspan="1" |'''关键参数'''
| 主要参数 || 外部输入 <math>I(t)</math>、恢复参数 <math>a,b</math>、时间尺度参数 <math>\varepsilon</math>
| colspan="1" rowspan="1" |刺激 I；恢复[[零倾线]]参数 a、b；[[时间尺度分离]]参数 ε≪1
|-
|-
| colspan="1" rowspan="1" |'''提出者 / 年份'''
| 历史来源 || FitzHugh（1961）的 Bonhoeffer–van der Pol 模型，Nagumo 等（1962）的等效电路
| colspan="1" rowspan="1" |R. FitzHugh（1961）；J. Nagumo 等（1962，电路实现）
|-
|-
| colspan="1" rowspan="1" |'''代表现象'''
| 典型行为 || 静息、可兴奋响应、弛豫振荡，以及特定参数下的多稳态与空间传播
| colspan="1" rowspan="1" |阈值兴奋、全或无[[锋电位]]、弛豫振荡、[[双稳态]]
|}
|}
== 历史溯源 ==
== 历史溯源 ==
1961 年，理查德·菲茨休（Richard FitzHugh）在范德波尔振荡器的基础上加入恢复项和刺激项，构造了能够表现稳定静息态、兴奋响应与重复发放的二维模型。他将其称为 Bonhoeffer–van der Pol 模型，简称 BVP 模型。这一工作的目的，是用较简单的数学结构分析神经膜的兴奋与恢复，而非证明具体离子通道对真实神经活动不重要。<ref name="FitzHugh_1961" />
1962 年，南云仁一（J. Nagumo）、S. Arimoto 和 S. Yoshizawa 用隧道二极管等元件构建了模拟神经轴突的有源脉冲传输线。该电路提供了相关兴奋动力学及脉冲传播的物理实现，后来这一模型通常被称为 FitzHugh–Nagumo 模型。<ref name="Nagumo_1962">{{cite journal|author1=J. Nagumo|author2=S. Arimoto|author3=S. Yoshizawa|title=An Active Pulse Transmission Line Simulating Nerve Axon|journal=Proceedings of the IRE|year=1962|volume=50|issue=10|pages=2061–2070|doi=10.1109/JRPROC.1962.288235|url=https://doi.org/10.1109/JRPROC.1962.288235}}</ref>
== 模型概述 ==
== 模型概述 ==
=== 方程与变量 ===
本文采用一种常见的无量纲写法：<ref name="Izhikevich_FitzHugh_2006">{{cite journal|author1=Eugene M. Izhikevich|author2=Richard FitzHugh|title=FitzHugh-Nagumo model|journal=Scholarpedia|year=2006|volume=1|issue=9|pages=1349|doi=10.4249/scholarpedia.1349|url=https://www.scholarpedia.org/article/FitzHugh-Nagumo}}</ref>
{{NumBlk|:|<math>\begin{aligned}
\frac{\mathrm{d}v}{\mathrm{d}t} &= v-\frac{v^3}{3}-w+I(t),\\
\frac{\mathrm{d}w}{\mathrm{d}t} &= \varepsilon\left(v+a-bw\right).
\end{aligned}</math>|{{EquationRef|1}}}}
其中，<math>t</math> 是无量纲时间。<math>v</math> 是与膜电位相关的兴奋坐标，<math>w</math> 是汇总恢复作用的坐标，二者均为标量。<math>I(t)</math> 表示经过标度变换的外部输入。<math>a\in\mathbb{R}</math> 决定恢复方程的偏置，<math>b>0</math> 调节恢复项，<math>\varepsilon>0</math> 控制两变量的相对变化速度。当其他量处于相近尺度且 <math>0<\varepsilon\ll1</math> 时，<math>v</math> 通常比 <math>w</math> 变化更快。无量纲变量不能直接按毫伏或毫秒解读，若要与实验比较，需要明确标度及参数标定。
<nowiki>$$\dfrac{dv}{dt}=v-\dfrac{v^{3}}{3}-w+I(t$$</nowiki>
三次非线性使兴奋变量能够产生再生性变化，恢复变量则提供较慢的负[[反馈]]。这种二维描述便于观察完整相轨迹，但不是将 HH 模型的四个变量直接合并后得到的精确等价模型。<ref name="Izhikevich_FitzHugh_2006" />
=== 不同文献的符号约定 ===
文献中存在恢复项写成 <math>\varepsilon(v-\tilde a-bw)</math> 的形式。若其他项保持不变，它与式（1）的对应关系是 <math>\tilde a=-a</math>，即参数的符号约定不同。只把 <math>w</math> 反号会同时改变两条方程中的相关项，不能用“恢复变量反号”解释只改动 <math>a</math> 的写法。同一个正数分别代入这两种约定，一般会得到不同的系统。
FitzHugh 的原始 BVP 方程使用 <math>x,y</math> 及不同时间尺度，不能将式（1）称为未经变换的原始公式。比较文献或程序时，应同时核对变量、时间和参数的变换。<ref name="FitzHugh_1961" />另一些文献把三次项写成 <math>-v^3</math>，相应的数值参数与稳定性公式也需调整，不能直接套用本文含 <math>-v^3/3</math> 的结果。<ref name="CebrianLacasa_2024">{{cite journal|author1=Daniel Cebrián-Lacasa|author2=Pedro Parra-Rivas|author3=Daniel Ruiz-Reynés|author4=Lendert Gelens|title=Six decades of the FitzHugh–Nagumo model: A guide through its spatio-temporal dynamics and influence across disciplines|journal=Physics Reports|year=2024|volume=1096|pages=1–39|doi=10.1016/j.physrep.2024.09.014|url=https://arxiv.org/abs/2404.11403}}</ref>
== 动力学机制 ==
== 动力学机制 ==
=== 零倾线与平衡点 ===
以下平衡点和分岔分析假定输入恒定，即 <math>I(t)=I</math>。两条零倾线分别由某一变量的时间导数为零定义。由式（1）得到
{{NumBlk|:|<math>\begin{aligned}
\dot v=0:&\quad w=v-\frac{v^3}{3}+I,\\
\dot w=0:&\quad w=\frac{v+a}{b}.
\end{aligned}</math>|{{EquationRef|2}}}}
它们的交点是平衡点，也称[[不动点]]。记其坐标为 <math>(v_*,w_*)</math>，代入式（2）可得
{{NumBlk|:|<math>\frac{b}{3}v_*^3+(1-b)v_*+a-bI=0,\qquad w_*=\frac{v_*+a}{b}.</math>|{{EquationRef|3}}}}
当 <math>0<b<1</math> 时，式（3）左端对 <math>v_*</math> 的导数为 <math>bv_*^2+1-b>0</math>，因此每个恒定输入对应唯一平衡点。当 <math>b>1</math> 时，某些 <math>a,I</math> 可对应三个不同平衡点。平衡点的个数本身不决定稳定性。
零倾线并不是轨迹必须沿着运动的曲线。在零倾线上的某一瞬间，只有对应的速度分量为零。对于小 <math>\varepsilon</math>，部分轨迹会在三次零倾线的吸引支附近缓慢运动，这是一种快慢系统近似。
=== 平衡点的局部稳定性 ===
由式（1）在平衡点处线性化，得到雅可比矩阵
{{NumBlk|:|<math>\mathbf{J}_*=\begin{pmatrix}
1-v_*^2 & -1\\
\varepsilon & -\varepsilon b
\end{pmatrix}.</math>|{{EquationRef|4}}}}
其中 <math>\mathbf{J}_*</math> 为二阶矩阵，其迹和行列式为
{{NumBlk|:|<math>\begin{aligned}
\operatorname{tr}\mathbf{J}_* &= 1-v_*^2-\varepsilon b,\\
\det\mathbf{J}_* &= \varepsilon\left(1-b+bv_*^2\right).
\end{aligned}</math>|{{EquationRef|5}}}}
在非临界情形下，<math>\det\mathbf{J}_*>0</math> 且 <math>\operatorname{tr}\mathbf{J}_*<0</math> 表示平衡点局部渐近稳定。行列式为负时，它是鞍点。行列式为正而迹为正时，它不稳定。迹或行列式为零的临界情形需要进一步分析。
[[文件:相平面.png|无框|522x522像素]]
固定 <math>w</math> 考察快子系统时，三次零倾线的外支 <math>|v|>1</math> 为吸引支，中支 <math>|v|<1</math> 为排斥支。这一结论依据快方程的导数 <math>1-v^2</math>，不能代替完整二维系统的稳定性判据。尤其是，交点位于中支并不必然意味着平衡点不稳定，也不单独保证存在稳定[[极限环]]。
=== Hopf 分岔与周期发放 ===
当控制参数变化时，平衡点可能发生 Hopf 分岔。由式（5）的零迹和正行列式条件，得到候选位置
{{NumBlk|:|<math>\begin{aligned}
v_{\mathrm H} &= \pm\sqrt{1-\varepsilon b},\\
I_{\mathrm H} &= \frac{v_{\mathrm H}+a}{b}-v_{\mathrm H}+\frac{v_{\mathrm H}^3}{3},\\
\det\mathbf{J}_{\mathrm H} &= \varepsilon\left(1-\varepsilon b^2\right)>0.
\end{aligned}</math>|{{EquationRef|6}}}}
式（6）中的两个非零候选值要求 <math>1-\varepsilon b>0</math>，并满足所列的正行列式条件。要确认为通常的非退化 Hopf 分岔，还需检查特征值实部随控制参数穿过零的横截条件及相应非线性条件。新生周期轨道是否稳定，也需要非线性分析。平衡点的局部失稳条件与稳定大幅周期发放的参数范围不能简单画等号。有关局部和全局[[分岔]]的讨论可参见相关综述。<ref name="CebrianLacasa_2024" />
=== 快慢运动与配图 ===
当参数支持稳定周期轨道且时间尺度分离较强时，轨迹通常在三次零倾线的两条外支附近缓慢运动，并在接近折点时快速跳转，形成弛豫振荡。对式（1），快子系统的折点位于 <math>v=\pm1</math>。有限 <math>\varepsilon</math> 下的实际轨迹及跳转位置不必与零倾线或折点完全重合。
[[文件:相平面.png|thumb|left|520px|图1：原词条相平面示意图，图中标注恒定输入 <math>I=0.8</math>。横轴为兴奋变量 <math>v</math>，纵轴为恢复变量 <math>w</math>。红色三次曲线和蓝色直线分别为式（2）的两条零倾线，灰色箭头表示向量场，绿色闭合轨迹表示图示参数下的周期运动。图中的恢复零倾线使用 <math>w=(v+a)/b</math>，与本文约定一致。“中支不稳定”指固定恢复变量时快子系统的性质。图源：站内原文件，由 Matthew 上传。]]
<div style="clear:both;"></div>
[[文件:时间响应.png|thumb|left|760px|图2：原词条时间响应示意图。横轴为无量纲时间 <math>t</math>，蓝线为 <math>v(t)</math>，红线为 <math>w(t)</math>。左图标注恒定输入 <math>I=0.1</math>，显示轨迹趋于静息态。右图标注恒定输入 <math>I=0.8</math>，显示持续振荡。这是两个恒定输入情形的比较，不能直接解释为短时阈下、阈上脉冲的比较。图源：站内原文件，由 Matthew 上传。]]
<div style="clear:both;"></div>
两幅原图的文件页未给出完整参数、初始条件和数值求解设置，因此这里保留其定性示意用途，不据此推断精确周期、阈值或分岔位置。
== 神经科学意义与典型现象 ==
== 神经科学意义与典型现象 ==
=== 可兴奋响应与恢复 ===
在可兴奋参数区域，系统具有稳定静息态。小扰动只产生较小响应，适当的较大扰动可以触发一次明显的相平面远足，随后回到静息态。恢复变量的慢变化会影响短时间内再次激发的难易程度。有限时间尺度分离下，响应振幅可随刺激连续变化，出现部分振幅的脉冲，因此不宜把 FHN 模型概括为严格的“全或无”模型。<ref name="Izhikevich_FitzHugh_2006" />
=== 持续发放与兴奋阻滞 ===
短时刺激触发的单次兴奋响应，应与恒定输入支持的持续周期发放区分。某些参数下，增大恒定输入可使系统从静息进入振荡。进一步增大输入又可使系统停留在较高兴奋水平的稳定平衡态，持续发放终止，这称为兴奋阻滞。因此，输入越强并不意味着始终产生更多脉冲。<ref name="Izhikevich_FitzHugh_2006" />
=== 参数示例 ===
为说明式（3）和式（5）的用法，取 <math>a=0.7</math>、<math>b=0.8</math>、<math>\varepsilon=0.08</math>。这是一组常见的示例参数。<ref name="Izhikevich_FitzHugh_2006" />按本文方程计算得到以下结果，数值保留四位小数：
{| class="wikitable"
! 恒定输入 <math>I</math> !! <math>v_*</math> !! <math>w_*</math> !! <math>\operatorname{tr}\mathbf{J}_*</math> !! <math>\det\mathbf{J}_*</math> !! 平衡点的局部性质
|-
| 0.1 || −1.1375 || −0.5469 || −0.3579 || 0.0988 || 渐近稳定
|-
| 0.8 || −0.2729 || 0.5339 || 0.8615 || 0.0208 || 不稳定
|-
| 1.7 || 1.1692 || 2.3364 || −0.4309 || 0.1035 || 渐近稳定
|}
式（6）给出的两个 Hopf 候选输入约为 <math>I_{\mathrm H}^{-}=0.3313</math> 和 <math>I_{\mathrm H}^{+}=1.4187</math>。表格仅说明平衡点的局部性质，不用于证明周期轨道的存在、稳定性或吸引范围。这里完整给出的参数也不等同于对原配图参数的确认。
== 空间推广与应用 ==
单个 FHN 单元只描述随时间变化的局部状态。若要研究兴奋信号在介质中的传播，可以加入扩散项，得到反应–扩散形式：<ref name="CebrianLacasa_2024" />
{{NumBlk|:|<math>\begin{aligned}
\frac{\partial v}{\partial t} &= v-\frac{v^3}{3}-w+I+D_v\nabla^2v,\\
\frac{\partial w}{\partial t} &= \varepsilon\left(v+a-bw\right)+D_w\nabla^2w.
\end{aligned}</math>|{{EquationRef|7}}}}
式（7）中的 <math>v=v(\mathbf{r},t)</math>、<math>w=w(\mathbf{r},t)</math> 为标量场，<math>\mathbf{r}</math> 为位置向量，<math>\nabla^2</math> 为拉普拉斯算子。<math>D_v>0</math>、<math>D_w\geq0</math> 为扩散系数。神经电生理的常见简化取 <math>D_w=0</math>，即只有兴奋变量发生空间扩散。具体解还取决于参数、空间维度、初始条件和边界条件。<ref name="CebrianLacasa_2024" />
在适当条件下，空间模型能够产生传播脉冲、波列和二维螺旋波。离散耦合的 FHN 单元则可用于研究神经网络中的同步及其他集体动力学。模型及其变体也用于心肌兴奋、化学可兴奋介质等问题，主要提供动力学层面的定性描述。空间脉冲的传播不要求每个未耦合单元本身处于持续振荡状态。<ref name="Nagumo_1962" /><ref name="CebrianLacasa_2024" />
== 优势与局限 ==
== 优势、局限与其他模型的关系 ==
FHN 模型变量少，便于用相平面分析兴奋和恢复过程，也便于构建空间系统或耦合网络。它适合研究静息、兴奋、振荡之间的转换及其几何机制。<ref name="FitzHugh_1961" />
模型中的恢复变量不等同于某个单独的离子通道门控概率。标准无量纲形式也不直接给出真实神经元的离子电流、动作电位持续时间或不应期长度。它可以表现某些恢复与不应性的定性特征，但定量预测仍需标定，并检验所选变体是否适合目标细胞或组织。与 HH 模型相比，FHN 保留的生物物理细节较少。<ref name="Izhikevich_FitzHugh_2006" />
从数学历史看，FHN 与范德波尔振荡器的联系体现在三次非线性、恢复动力学和弛豫振荡。从神经建模看，它与 HH 模型的联系是对兴奋行为的简化描述。仅凭其他模型同样具有兴奋或抑制变量，不能将其认定为 FHN 的直接离散化或严格等价形式。
FHN 是 HH 模型的低维几何化简化，上承 HH 的兴奋机制、下启各类低维可兴奋系统。[[Rulkov 映射]]可视为其离散时间对应物；[[威尔逊-考恩模型]]在群体层面沿用了"快兴奋—慢恢复/抑制[[反馈]]"的思想。在只需判断"是否兴奋、是否振荡、在何种参数下转换"的研究问题中，FHN 往往比 HH 更高效。
== 参考文献 ==
== 参考文献 ==
<references />
== 署名 ==
杨明哲