开放科学（资源服务）标志码（OSID）：

0 概述

随着生物电子技术的发展，人工神经网络的研究也不断取得进展。文献[1]根据哺乳动物视觉皮层神经元信号传导特性，构建一种新的神经网络模型。文献[2]将神经网络模型改进为用于图像处理领域的模型，并命名为脉冲耦合神经网络（Pulse Coupled Neural Network，PCNN）。传统的PCNN模型可分为耦合连接和无耦合连接两种模式，被广泛应用于图像分割^[3-4]、图像融合^[5-6]、图像去噪^[7-8]、图像检索^[9-10]以及图像边缘检测^[11-12]等领域。通常将无耦合连接模式下的网络称为无耦合PCNN。无耦合PCNN模型虽然没有耦合连接项，但是保留了PCNN模型的核心子系统部分，具有网络参数少的特点，近年来，逐渐被应用于图像处理领域^[13]。

研究人员通过对图像进行信息提取或变换域处理，融合PCNN模型进行参数自适应设定，并应用于特定的图像处理领域。文献[14]通过SIST（Shift Invariant Shearlet Transform）方法获取图像的高频子带和低频子带频率，使用SOM（Self-Origanizing Maps）网络进行初步分类，之后融合改进的PCNN模型，获得较优的图像分割效果。文献[15]利用非下采样轮廓波变换对图像进行多尺度和多方向分解，获得低通子带和多个带通方向子带，并将其与改进后的PCNN模型相结合用于多聚焦图像融合。文献[16]基于图像NSCT（Non-Subsampled Contourlet Transform）域变换提取图像的边缘特征，提出一种结合自适应PCNN的红外与可见光图像融合算法。文献[17]通过建立待分割图像统计特性与PCNN神经元动态特性之间的直接关系，完成简化PCNN的参数自适应设置。文献[18]通过计算分割图像的最小交叉熵，提出一种新的基于循环迭代的自适应图像分割算法，用于改进PCNN阈值分割机制。此外，研究人员通过对网络模型进行理论分析，以解释网络本质特性并进行参数的自适应设定。文献[19]分析改进后PCNN网络的动态行为，明确PCNN的动态行为特性受网络参数变化的影响，并提出网络参数确定准则，使得所有神经元只点火一次。文献[20]通过分析PCNN中神经元的脉冲发放周期，定义了动态比，以得出近似的脉冲发放周期公式。文献[21]根据PCNN工作机理的时域分析，提出一种在点火神经元没有完全灭火时不允许灭火神经元重新点火的网络工作模式，并用于图像的自适应分割。然而，上述PCNN模型的研究都只是从时域方面或图像本身所含信息方面分析参数与模型之间的关系，在网络本质的特性挖掘方面还存在着一定的局限性和片面性。

本文从频域角度对无耦合PCNN的网络参数与模型特性之间的关系进行分析。通过推导动态门限子系统的频率响应函数，揭示动态门限子系统具有低通特性，同时分析参数对动态门限子系统频率响应的影响，推导无耦合连接时PCNN有效点火时刻方程、动态门限衰减频率和脉冲发放频率的方程，解释无耦合PCNN的网络参数对网络频域特性的影响。

1 基本概念 1.1 数学模型

无耦合PCNN模型由输入子系统、动态门限子系统和输出子系统组成，分别如式（1）~式（3）所示^[22]：

$ {U}_{ij}\left(n\right)={S}_{ij} $

(1)

$ {E}_{ij}\left(n\right)={{\mathrm{e}}}^{-{a}_{E}}{E}_{ij}(n-1)+{v}_{E}{Y}_{ij}(n-1) $

(2)

$ {Y}_{ij}\left(n\right)=\varepsilon \left({U}_{ij}\right(n)-{E}_{ij}(n\left)\right) $

(3)

其中：$ i $、$ j $为本神经元像素坐标点；$ n $为迭代次数；$ {S}_{ij} $为系统的输入，常为归一化之后的像素灰度值；$ {a}_{E} $、$ {v}_{E} $为控制动态门限子系统的2个网络参数。模型的输出子系统由单位阶跃函数构成，决定该神经元的脉冲发放状态。当神经元的输入值大于动态门限值时，神经元发放脉冲又称为神经元点火；否则，神经元灭火。

1.2 第1次有效点火时刻

在前期研究中^[23]已经得知无耦合PCNN的第1次神经元点火时刻不受网络参数的影响，因此重点分析了网络的第2次点火时刻，又称为第1次有效点火时刻（所有输入值在第1次迭代时发生的点火称为无效点火，其后各次点火为有效点火），通过推导所得第1次有效点火时刻$ {N}_{1} $如式（4）所示：

$ {N}_{1}={\mathrm{c}}{\mathrm{e}}{\mathrm{i}}{\mathrm{l}}\left(2+\frac{1}{{a}_{E}}{\mathrm{l}}{\mathrm{n}}\left(\frac{{v}_{E}}{{S}_{ij}}\right)\right) $

(4)

由于网络数学模型为差分方程，因此当理论点火时刻为小数时，实际点火时刻为下一个整数时刻。在第1次有效点火时刻计算中引入向上取整函数$ {\mathrm{c}}{\mathrm{e}}{\mathrm{i}}{\mathrm{l}} $。由式（4）可知，当灰度值$ {S}_{ij}=0 $时，无法发生第1次有效点火。本文主要针对网络输入的有效点火时刻及输入灰度值$ {S}_{ij} $为1/255~1时，网络所表现出的特性。

2 无耦合PCNN频域分析

在实际图像的处理过程中，输入子系统的输入为图像对应像素点的灰度值，是简单的线性系统。本文对式（1）表示的输入子系统不展开分析。由于无耦合网络的时域特性完全受参数影响，因此本文对于网络的频域分析，根据动态门限子系统的频域特性分析参数对动态门限子系统的影响，将式（1）、式（2）和式（3）联系起来，分别分析网络中动态变量$ {Y}_{ij}\left(n\right) $和$ {E}_{ij}\left(n\right) $的频域特性与参数之间的关系。

2.1 动态门限子系统的频域特性 2.1.1 系统函数

动态门限子系统是无耦合PCNN模型的核心子系统，由差分方程式（2）表示，对式（2）进行$ z $变换可得式（5）：

$ {E}_{ij}\left(z\right)={{\mathrm{e}}}^{-{a}_{E}}{z}^{-1}{E}_{ij}\left(z\right)+{v}_{E}{z}^{-1}{Y}_{ij}\left(z\right) $

(5)

动态门限子系统的系统函数如式（6）所示：

$ {H}_{E}\left(z\right)=\frac{E\left(z\right)}{Y\left(z\right)}=\frac{{v}_{E}{z}^{-1}}{1-{{\mathrm{e}}}^{-{a}_{E}}{z}^{-1}}=\frac{{v}_{E}}{z-{{\mathrm{e}}}^{-{a}_{E}}} $

(6)

根据系统函数$ {H}_{E}\left(z\right) $可知，动态门限子系统没有零点，只有1个极点$ z={{\mathrm{e}}}^{-{a}_{E}} $。此时系统函数的收敛域为$ \left|z\right|\in [{{\mathrm{e}}}^{-{a}_{E}}, {\mathrm{\infty }}] $。当$ {a}_{E} $ > 0时，子系统的收敛域包含单位圆，此时的系统是一个因果稳定系统。当$ {a}_{E} $≤0时，系统不稳定。因此，无耦合PCNN参数$ {a}_{E} $的取值范围为$ (0, +{\mathrm{\infty }}) $。

从零极点对系统影响可知，子系统函数无零点，且极点$ {{\mathrm{e}}}^{-{a}_{E}}\in \left({\mathrm{0, 1}}\right) $在虚轴右边，因此动态门限子系统具有低通特性。由于子系统的极点在单位圆内，因此动态门限子系统是一个衰减系统。$ {{\mathrm{e}}}^{-{a}_{E}} $为衰减因子，且$ {a}_{E} $越大，极点离单位圆越远，系统的衰减速度越快，$ {a}_{E} $越小极点距离单位圆越近，系统的衰减速度越慢。无耦合PCNN模型的网络参数$ {a}_{E} $通过控制极点的大小，影响动态门限$ {E}_{ij}\left(n\right) $的衰减速度。

2.1.2 频率响应与参数影响分析

动态门限子系统的系统函数$ {H}_{E}\left(z\right) $收敛域包含单位圆，根据$ z $变换与傅里叶变换之间的关系可得子系统的频率响应，如式（7）所示：

$ {H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)=\begin{array}{cc}\left(\frac{{v}_{E}}{{{\mathrm{e}}}^{{\mathrm{j}}\omega }-{{\mathrm{e}}}^{-{a}_{E}}}\right), {a}_{E} > 0& \begin{array}{cc}& \end{array}\end{array} $

(7)

根据频率响应函数$ {H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}w}\right) $可得其幅频响应，如式（8）所示：

$ \left|{H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)\right|=\left|\frac{{v}_{E}}{{{\mathrm{e}}}^{{\mathrm{j}}\omega }-{{\mathrm{e}}}^{-{a}_{E}}}\right|=\left|\frac{{v}_{E}}{\sqrt{1-2{{\mathrm{e}}}^{-{a}_{E}}{\mathrm{c}}{\mathrm{o}}{\mathrm{s}}\omega +{{\mathrm{e}}}^{-2{a}_{E}}}}\right| $

(8)

由式（8）可知，参数$ {a}_{E} $和$ {v}_{E} $共同控制子系统的幅值特性。改变参数$ {a}_{E} $和$ {v}_{E} $，在$ \omega \in (0, \pi ) $范围内观察到的子系统$ {H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right) $的幅频特性曲线变化如图 1所示。

从图 1（a）可以看出，参数$ {a}_{E} $主要影响幅频响应峰值与谷值之间的差值，$ {a}_{E} $越小，差值越大，此时系统低频分量的幅值较大，高频分量的幅值较小，使得子系统具有较优的低通特性。从图 1（b）可以看出，参数$ {v}_{E} $主要影响幅频特性的幅值大小，$ {v}_{E} $越大，幅频响应的幅值也将增大。因此，无耦合PCNN模型的动态门限子系统具有低通特性，不同的网络参数使得动态门限系统表现出不同的低通特性。子系统的低通特性主要受参数$ {a}_{E} $影响，参数$ {v}_{E} $主要影响子系统的幅值大小。

2.2 网络脉冲发放频域特性 2.2.1 有效点火脉冲方程

无耦合PCNN是一个非线性系统，其神经元的输入与神经元的脉冲发放频率之间存在非线性关系。当参数$ {a}_{E} $=0.5，$ {v}_{E} $=5时，在不同灰度值下神经元的动态门限和点火时刻的变化曲线如图 2所示。从图 2可以看出，单个神经元的脉冲点火过程和动态门限衰减过程都具有周期性，在输入值$ {U}_{ij} $的激励下，输出$ {Y}_{ij}\left(n\right) $和$ {E}_{ij}\left(n\right) $以一定的周期发放脉冲。

因此，输入值与神经元点火周期和动态门限衰减存在非线性映射关系。为分析这种非线性关系，在单个神经元的一个有效点火周期内建立关系，如式（9）所示：

$ {E}_{ij{\mathrm{m}}{\mathrm{i}}{\mathrm{n}}}={U}_{ij}\left(n\right) $

(9)

当单个神经元的输入值与动态门限最小值相等时，神经元发生点火。由式（2）可知，动态门限的最大值是脉冲点火时刻后下一次的门限值，$ {E}_{ij{\mathrm{m}}{\mathrm{a}}{\mathrm{x}}}={{\mathrm{e}}}^{-{a}_{E}}{S}_{ij}+{v}_{E} $，最小值为$ {E}_{ij{\mathrm{m}}{\mathrm{i}}{\mathrm{n}}}={S}_{ij} $，设两次脉冲点火时刻的间隔为$ I $，则由式（9）可得式（10）：

$ {S}_{ij}=({v}_{E}+{S}_{ij}{{\mathrm{e}}}^{-{a}_{E}}){{\mathrm{e}}}^{-{a}_{E}I} $

(10)

间隔$ I $如式（11）所示：

$ I=\frac{1}{{a}_{E}}{\mathrm{l}}{\mathrm{n}}\left(\frac{{v}_{E}}{{S}_{ij}}+{{\mathrm{e}}}^{-{a}_{E}}\right) $

(11)

脉冲发放间隔应满足$ I > 0 $，由此可得$ {v}_{E} > (1-{{\mathrm{e}}}^{-{a}_{E}}){S}_{ij{\mathrm{m}}{\mathrm{a}}{\mathrm{x}}} $。其中：$ {S}_{ij{\mathrm{m}}{\mathrm{a}}{\mathrm{x}}} $表示输入最大值。当输入的最大值与$ {v}_{E} $和$ {a}_{E} $满足以上关系时，所有的输入都能表现出良好的脉冲周期发放特性。由脉冲发放间隔得到脉冲发放周期，如式（12）所示：

$ {T}_{{s}_{ij}}=I+1=\frac{1}{{a}_{E}}\left({\mathrm{l}}{\mathrm{n}}\left(\frac{{v}_{E}}{{S}_{ij}}+{{\mathrm{e}}}^{-{a}_{E}}\right)\right)+1 $

(12)

不同输入$ {S}_{ij} $的点火周期$ {T}_{{S}_{ij}} $都为正整数，同样的引入$ {\mathrm{c}}{\mathrm{e}}{\mathrm{i}}{\mathrm{l}} $向上取整函数。因此，无耦合PCNN的实际脉冲发放周期近似为：

$ T={\mathrm{c}}{\mathrm{e}}{\mathrm{i}}{\mathrm{l}}\left({T}_{{s}_{ij}}\right)={\mathrm{c}}{\mathrm{e}}{\mathrm{i}}{\mathrm{l}}\left(\frac{1}{{a}_{E}}\left({\mathrm{l}}{\mathrm{n}}\left(\frac{{v}_{E}}{{S}_{ij}}+{{\mathrm{e}}}^{-{a}_{E}}\right)\right)\right)+1 $

(13)

无耦合PCNN在得到能够反映实际脉冲发放的周期后，脉冲发放时刻如式（14）所示：

$ {Y}_{ij}\left(n\right)=\delta (n-{N}_{1}-rT)\text{，}r\in {\mathrm{1, 2}}, \cdots , n $

(14)

其中：单位冲激函数$ \delta \left(\begin{array}{c}n\end{array}\right) $为脉冲发放；$ {N}_{1} $为第1次有效点火时刻；$ r $为正整数。

2.2.2 脉冲发放频率与参数影响分析

由脉冲发放时刻方程可知，无耦合PCNN脉冲点火具有周期性，而周期为N的序列可展开成离散傅里叶级数（Discrete Fourier Series，DFS）。因此，在一个周期内对脉冲发放时刻式（14）进行傅里叶级数展开，如式（15）所示：

$ {Y}_{ij}\left(k\right)={\mathrm{D}}{\mathrm{F}}{\mathrm{S}}\left({Y}_{ij}\right(n\left)\right)=\sum \limits_{n=0}^{N}{Y}_{ij}\left(n\right){{\mathrm{e}}}^{-{\mathrm{j}}\frac{2{\mathrm{\pi }}}{N}kn}=\sum\limits_{n={N}_{1}}^{T+{N}_{1}-1}\delta (n-{N}_{1}){{\mathrm{e}}}^{-{\mathrm{j}}\frac{2{\mathrm{\pi }}}{T}kn}\\={{\mathrm{e}}}^{-{\mathrm{j}}\frac{2{\mathrm{\pi }}}{T}k{N}_{1}} $

(15)

本文对脉冲发放时刻的傅里叶级数表达式进行傅里叶变换（Fourier Transform，FT），可得其脉冲发放频率表达式，如式（16）所示：

$ {Y}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)={\mathrm{F}}{\mathrm{T}}\left({Y}_{ij}\right(k\left)\right)=\frac{2{\mathrm{\pi }}}{N}\sum\limits _{k=-\infty }^{k=\infty }{Y}_{ij}\left(k\right)\delta \left(\omega -\frac{2{\mathrm{\pi }}}{N}k\right)=\\\frac{2{\mathrm{\pi }}}{T}\sum\limits _{k=-{\mathrm{\infty }}}^{k={\mathrm{\infty }}}{{\mathrm{e}}}^{-{\mathrm{j}}\frac{2{\mathrm{\pi }}}{T}k{N}_{1}}\delta \left(\omega -\frac{2{\mathrm{\pi }}}{T}k\right) $

(16)

幅值谱函数如式（17）所示：

$ \left|{Y}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)\right|=\frac{2{\mathrm{\pi }}}{T}\delta \left(\omega -\frac{2{\mathrm{\pi }}}{T}k\right) $

(17)

由式（17）可知，无耦合PCNN的脉冲发放幅值谱是由一系列幅值为$ 2{\mathrm{\pi }}/T $和间隔为$ 2{\mathrm{\pi }}/T $的冲击序列构成。当选取输入$ {S}_{ij}=20/255 $时，$ \omega \in (0, {\mathrm{\pi }}) $范围幅频特性随参数的$ {a}_{E} $和$ {v}_{E} $变化趋势如图 3和图 4所示。根据离散周期信号频谱特性可知，频域谱线越密，时域脉冲周期越大，频域谱线越疏，时域脉冲周期越小。从图 3和图 4可以看出，参数变化不会引起脉冲发放频谱函数幅值的变化。但是参数$ {a}_{E} $的增大使得脉冲发放频谱函数的谱线间隔增大，从而增大脉冲点火时刻周期。参数$ {v}_{E} $的增大反而使得脉冲发放频谱函数的谱线间隔变小，脉冲点火时刻周期减小。

因此，在无耦合PCNN模型中，参数$ {a}_{E} $的增大使得脉冲发放时刻周期变小，参数$ {v}_{E} $的增大会使得脉冲发放时刻周期变大，但不影响脉冲发放的强度。

2.2.3 网络动态门限衰减频率与参数影响分析

动态门限子系统是一个衰减系统，因此将该子系统的输出动态变量$ {E}_{ij}\left(n\right) $的频率称为动态门限衰减频率。动态门限衰减频率与脉冲发放频率存在以下关系，如式（18）所示：

$ {E}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)={Y}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)\times {H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right) $

(18)

由于式中$ {H}_{E}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right) $和$ {Y}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right) $已由式（7）和式（15）给出，因此动态衰减频率的频域方程如式（19）所示：

$ {E}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)=\frac{2{\mathrm{\pi }}}{T}{{\mathrm{e}}}^{-{\mathrm{j}}\frac{2{\mathrm{\pi }}}{T}k{N}_{1}}\delta \left(\omega -\frac{2{\mathrm{\pi }}}{T}k\right)\left(\frac{{v}_{E}}{{{\mathrm{e}}}^{{\mathrm{j}}\omega }-{{\mathrm{e}}}^{-{a}_{E}}}\right) $

(19)

频谱函数如式（20）所示：

$ \begin{array}{l}\left|{E}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)\right|=\frac{2{\mathrm{\pi }}}{T}\delta \left(\omega -\frac{2{\mathrm{\pi }}}{T}k\right)\\ \frac{{v}_{E}}{\sqrt{1-2{{\mathrm{e}}}^{-{a}_{E}}{\mathrm{c}}{\mathrm{o}}{\mathrm{s}}\omega +{{\mathrm{e}}}^{-2{a}_{E}}}}\end{array} $

(20)

当选取输入$ {S}_{ij}=20/255 $时，参数$ {a}_{E} $和$ {v}_{E} $增大在$ \omega \in (0, {\mathrm{\pi }}) $范围幅频特性随参数变化趋势如图 5和图 6所示。

从式（20）可知，动态门限幅频函数$ \left|{E}_{ij}\left({{\mathrm{e}}}^{{\mathrm{j}}\omega }\right)\right| $的频谱谱线间隔与脉冲发放频谱函数的谱线间隔相同，参数对动态门限周期与脉冲发放周期的影响程度相同。从图 5和图 6可以看出，与脉冲发放频谱相比，随着参数$ {a}_{E} $的增大，动态门限幅频函数的频谱主要由高频分量组成，使得动态门限的衰减频率变大，衰减速度变快。参数$ {v}_{E} $的增大反而会使得动态门限幅频函数的频谱主要由低频分量组成，动态门限$ {E}_{ij}\left(n\right) $的衰减频率变低。

因此，在无耦合PCNN模型中，参数$ {a}_{E} $和参数$ {v}_{E} $的变化会影响动态门限衰减周期及动态门限衰减频率的高低。

3 实验仿真与验证

本文通过式（7）、式（16）和式（19）描述了无耦合PCNN模型的频域特性，通过分析参数对无耦合PCNN模型频域特性的影响，得出参数影响网络运行的结论。在主频2.8 GHz（双核），内存12 GB的计算机上利用Python3.8平台完成数值模拟实验，验证了无耦合PCNN频域相关公式推导及理论分析结论的正确性。

3.1 实际网络脉冲发放时刻

为验证对脉冲发放时刻推导的正确性，本文构建一幅灰度值已知的人工测试灰度图与灰度值，如图 7所示。

当网络参数$ {a}_{E} $=0.9，$ {v}_{E} $=15时，无耦合PCNN实际迭代过程如图 8所示。从图 8可以看出，第1次迭代为网络的无效点火时刻，之后的各次点火均为有效点火。灰度值在110~250之间的像素第1次有效点火时刻为第6次迭代，点火周期$ T=5 $。灰度值在50~100之间的像素第1次有效点火时刻为第7次迭代，点火周期$ T=6 $。灰度值在20~40之间的像素第1次有效点火时刻为第8次迭代，点火周期$ T=7 $。灰度值为$ 10 $的像素第1次有效点火时刻为第9次迭代，点火周期$ T=8 $。根据第1次点火公式（4）可计算出各像素第1次有效点火时刻，根据点火周期公式（13）可计算出各像素点火周期。不同像素的第1次有效点火时刻和点火周期对比如图 9所示。从图 9可以看出，本文所推导公式的脉冲发放过程与人工测试图实际迭代过程相一致。

3.2 实际网络参数对脉冲发放与动态门限的影响

本文分别选取灰度值$ {S}_{ij}=10/255 $，$ {S}_{ij}=110/255 $作为对网络中动态变量$ {Y}_{ij}\left(n\right) $受参数影响的观测对象。当网络迭代第60次时，其神经元脉冲发放过程如图 10所示。从图 10可以看出：在灰度值$ {S}_{ij} $相同的情况下，保持参数$ {v}_{E}=15 $不变，$ {a}_{E} $从0.5增大到0.9会使得网络脉冲发放时刻周期变小；保持参数$ {a}_{E}=0.5 $不变，$ {v}_{E} $从15增大到30使得网络脉冲发放时刻周期变大。因此，在实际的脉冲发放过程中，参数对脉冲发放频率的影响与理论推导分析相一致。

本文分别选取灰度值$ {S}_{ij}=10/255 $，$ {S}_{ij}=110/255 $作为网络中动态变量$ {E}_{ij}\left(n\right) $受参数影响分析的观测对象。当网络迭代第60次时，不同灰度值下网络实际动态门限变化过程如图 11所示。从图 11可以看出：在输入$ {S}_{ij} $相同的情况下，保持参数$ {v}_{E}=15 $不变，$ {a}_{E} $从0.5增大到0.9，使得动态门限$ {E}_{ij}\left(n\right) $衰减变快，衰减周期变小；保持参数$ {a}_{E}=0.5 $不变，$ {v}_{E} $从15增大到30，使得动态门限$ {E}_{ij}\left(n\right) $的幅值增益变大，衰减变慢，衰减周期变大。因此，动态变量$ {E}_{ij}\left(n\right) $的迭代过程与之前分析结果一致。

4 结束语

本文分析了无耦合PCNN的频域特性。通过对无耦合PCNN动态门限子系统进行分析，得出动态门限子系统稳定的条件，确定参数的取值范围，并指出动态门限子系统具有低通特性。同时分析单个神经元的脉冲发放频率特性和动态阈值频域特性，明确参数对网络动态频域的影响。仿真实验结果验证理论分析结论的正确性。后续将从图像特性角度分析图像统计特性对网络参数的影响，设计无耦合PCNN模型用于图像分割时的参数自适应设置方法，使无耦合PCNN模型广泛应用于图像处理领域。