文献 [9] 首次将协同滤波算法用于混沌信号去噪并且提出混沌信号协同滤波去噪算法，算法包括三个步骤：相似块分组、协同滤波和聚合重构。

1) 相似块分组

设参考块为R，将任一块S与R的相似度通过它们之间的归一化距离进行量化，其中使用二范数来表示 ${\Vert \cdot \Vert }_2$，并以块的宽度为单位；d代表相似块间的距离。当此距离d较小时，表明块S与参考块R之间的相似度较高，可表示为：

$d\left(S,R\right)={\Vert S-R\Vert }_2^2/\omega$(1)

在设置搜索窗口长度为l的条件下，以参考块R为中心，寻找与之距离最小的m个块进行分组。这些选中的块被存储为 $m\ast \omega$的二维数组形式。随着参考块R从观测信号的起始端逐步以 $\delta$为步长向终端移动，搜索窗口也相应地移动。在每一个位置，都会记录下形成的相似块组 $group\left(R\right)$，并为分组内的每个块标注其具体位置。

2) 协同滤波

协同滤波算法的过程包括三个主要步骤。首先是对时域信号进行变换(比如离散余弦变换)，得到频域上的值；其次，给定阈值对频域值进行过滤，最后将过滤之后的频域值通过逆变换转换为时域值，从而达到滤波的效果。

在对混沌信号进行处理时，在算法的初始阶段，每个分组都会进行二维变换，以获得变换系数矩阵 $G=\left[g_{i,j}\right]\in I^{m\times \omega }$。见公式(2)：

$G\left(R\right)=T\left\{group\left(R\right)\right\}$(2)

其中，T[·]指的是二维离散余弦变换(DCT)， $G\left(R\right)$为分组块对应在变换频域上的值。对频域上的值进行阈值化处理，类似于高通滤波，大于阈值的值被保留，小于阈值的值被舍弃，见公式(3)：

$H_T\left(g\right)=\{\begin{array}{l}g\text{,}\left|g\right|\succ \lambda \left(R\right)\\ o\text{,}\left|g\right|\le \lambda \left(R\right)\end{array}$(3)

上述公式中的阈值是由文献 [11] 提出的VisuShrink方法所求得的，其定义见公式(4)：

$\lambda \left(R\right)=\sigma \sqrt{2\log \left(m\times \omega \right)}$(4)

其中， $\sigma$代表观测信号中噪声的标准偏差，该值是基于系数矩阵G(R)的中位数绝对偏差来估算的。估算表达式见公式(5)：

$\hat{\sigma }\left(R\right)=median\left(G\left(R\right)/0.6745\right)$(5)

最后，将得到的阈值化处理结果进行离散余弦变换的逆变换，将信号从频域还原为到时域，见公式(6)：

$group\left(R\right)=T^{-1}\left\{H_T\left(G\left(R\right)\right)\right\}$(6)

其中， $T^{-1}$表示离散余弦变换的逆变换。

3) 聚合重构

在聚合重构过程中，考虑到相似块之间的重叠性，一个特定的信号点往往会出现在多个相似块中。该信号点的重构值是通过计算所有包含此点的相似块在该位置上的滤波输出的算术平均值来实现的，见公式(7)：

$x\left(n\right)=\frac{{\displaystyle \underset{R}{\sum }{\displaystyle \underset{S\in group\left(R\right)}{\sum }{F}_S\left(n\right)x_{S,R}\left(n\right)}}}{{\displaystyle \underset{R}{\sum }{\displaystyle \underset{S\in group\left(R\right)}{\sum }{F}_S\left(n\right)}}}$(7)

其中， $x_{S,R}\left(n\right)$表示在信号点n处，分组内相似块S的滤波输出结果， $F_S\left(n\right)$表示第n个信号点是否包含在S块中，见公式(8)：

$F_S\left(n\right)=\{\begin{array}{l}1\text{,}n\in S\\ 0\text{,}n\notin S\end{array}$(8)

本文在文献 [9] 的基础上，设计了一种错位搜索的方法来优化相似块的选取，并采用DTW方法对聚合重构方法进行改进。改进后的自适应协同滤波流程图如 图1 所示。

1) 基于排列熵的参数优化

文献 [10] 通过引入排列熵，对协同滤波中的参数进行优化，效果显著。本文也采用排列熵的方法，对协同滤波中的参数进行优化。排列熵(PE)是一种评估混沌信号复杂度的有效工具。混沌信号的复杂度越高，其PE值亦随之增大；相反，具有较低复杂度的混沌信号，其PE值也相对较低。当混沌信号遭受噪声干扰时，其复杂度会提升，并且随着噪声水平的提高，复杂度也会相应增加。 表1 展示了原始信号与在不同噪声水平下的预测误差PE情况。

表1 中分别给出了Lorenz系统和Chen系统生成的模拟混沌信号在不同信噪比条件下的PE值。随着噪声水平的增高，含噪混沌信号的复杂度也会相应上升。

2) 基于错位搜索方法的相似块优化

在信号处理的众多应用场景中，如滤波和小波变换等传统技术已被广泛采用以提高信号品质。这些建立在信号频率属性基础上的方法，有效地增强了信号的可解读性和信息精度。本文提出一种基于错位搜索方法的相似数据块匹配及调整的信号处理策略，通过识别并精确调整信号中的相似数据块，提高去噪的效果。

通过错位搜索对相似块进行微小的位置调整，我们可以进一步提高相似度，从而更有效地去除噪声、进行信号压缩或修复图像，其流程图如 图2 所示。

首先将接收到的信号分割成若干个预定尺寸的块，其中块的尺寸基于最佳块宽度进行调整。针对信号内的每一块，本研究计算其与其他所有块的相似程度。基于所得到的相似度，本研究为每个目标块挑选出最为相近的若干块。此环节可能需设置特定阈值或按相似度高低选择若干块作为最匹配块，依据文献 [10] 的方法，本文选取了前30个最为相似的块。

对每一选中的相似块，确定其位置调整的可能范围，该范围既可以是预设固定值，亦可依据块的具体特性进行动态设定。位置优化调整旨在通过细微位移寻找更佳的匹配点，进而提升块间的相似度水平。在规定的调整范围内，对每个相似块的位置进行精细调整。此过程通常包括在目标块周边的有限区域内对相似块进行位移，并在每次位移后重新评估其与目标块的相似度。对于每一次的位置调整，均需计算调整后的块与目标块间的相似度，并通过比较不同调整位置的相似度，选出相似度最高的位置作为相似块的最终位置。如 图3 所示，某些相似块选择上虽然形状相似但是首位差别大互相错开，因此经过优化调整后相似程度得到较大提升。

3) 基于DTW的聚合重构优化

在信号重建的应用场景中，DTW主要被用来优化信号片段间的过渡效果。它通过寻找两个信号片段之间的最优对齐路径，能够生成平滑且自然的过渡效果，进而提升重建信号的整体质量。如 图4 所示，序列A中的节点i可能并不与序列B中的任何节点直接对应，这种情况下的误差可能会显著影响到时间

序列相似度的准确评估。

与传统的欧几里得距离方法不同，DTW通过在时间轴上对时间序列进行伸缩调整，以寻求两个序列之间的最佳匹配点，这一特性使得DTW尤其适合于处理那些在长度、速度或时间节点上有所差异的时间序列数据。动态时间规整简介DTW是一种算法( 图5 )，用于测量两个时间序列之间的相似度，与传统的欧几里得距离不同，DTW能够考虑到时间序列的形状和模式，即使两个序列在时间轴上有所偏移 [12] 。这使得DTW成为处理时间序列数据，特别是在信号重建领域的一个强大工具。

滤波参数的选择受到信号特性、采样频率和噪声水平的影响。为了验证所提出方法的有效性，本研究应用该方法在不同信号特性、采样频率和噪声水平的条件下进行滤波参数的选取。为了便于分析，本文还对排列熵与窗口大小之间的关系以及去噪效能与窗口大小之间的关系进行了图形化展示。此外，采用均方误差(MSE)作为评估去噪效能的标准，即通过对比观测信号与经过滤波处理后的信号之间的差异进行评价。

1) 不同信号特征

考虑到不同混沌信号频谱宽度的异质性，选择适当的块宽度需依据各类混沌信号的特征进行调整。本文以Lorenz混沌信号和Chen混沌信号为研究对象，进行了详细分析。相较于Lorenz信号，Chen信号展现了更广泛的频谱范围。在Lorenz和Chen信号上分别施加零均值高斯白噪声，其中信噪比设定为10 db，采样周期设定为0.01 s。通过应用不同宽度的滤波块对这两种带噪声的混沌信号进行降噪处理，并对处理后的信号进行排列熵计算。如 图6 所示，降噪后信号的排列熵反映了滤波块宽度对降噪性能的影响。在特定的采样频率和噪声水平下，理想的滤波块宽度对于Lorenz信号和Chen信号分别约为110和60，这与先前的分析结果一致。

2) 不同采样频率

在较高的采样频率下，混沌信号的协同滤波去噪方法展现出更优的去噪效能，并且对于不同采样频率，所需的最佳滤波块尺寸也存在差异。通过将采样时间(t)设定为0.02、0.01、0.0075、0.005秒，生成Lorenz信号，并向其添加零均值的高斯白噪声，其中信噪比固定为10分贝。然后，采用各种宽度的滤波块对在不同采样频率下的含噪信号执行去噪操作。正如 图7 与 表2 展现的，随着采样频率的升高，降噪性能逐渐提升，与此同时，最优的滤波块宽度也相应增大。

以Lorenz混沌信号为研究对象，通过向Lorenz信号中添加零均值的高斯白噪声，并设置信噪比为

图7. 采样频率对块宽的影响：(a) 不同采样频率对应的PE；(b) 不同采样频率的去噪效果

10分贝，采样时间为0.01秒。基于排列熵算法确定的最佳块宽度作为参数，对相似块进行匹配。本研究分别采用了错位搜索相似块算法和未经优化的相似块匹配算法进行比较分析，其结果展示如 图8 和 表3 所示。

通过对比上述结果的图表和数据，可以明显观察到，采用错位搜索相似块算法优化后的自适应协同滤波算法在去噪性能上有了显著提高。进一步地，在不同频率条件下的对比分析表明，经过错位搜索相似块算法处理后的协同滤波去噪效能显著增强，具体对比结果如 表4 所示。

以Lorenz混沌信号为例，给Lorenz信号加上零均值高斯白噪声，信噪比取10 dB，采样时间取0.01 s，运用不同的聚合重构方法进行仿真实验，试验结果如 表5 所示。

根据上表所示，通过DTW聚合重构方法相较于传统的直接平均聚合重构方法，在精度上实现了显著提高，并且与其他方法相比展现出了明显的优越性。此次仿真实验验证了该方法的实用性和优势所在。

1) 不同信号特性下与滤波效果的对比分析

为全面评价该方法的性能，我们将其与传统的协同滤波去噪技术进行了比较。选取洛伦兹混沌系统作为信号源，对系统产生的信号分别加入了三种不同水平的高斯白噪声，分别为低、中、高三个噪声等级，通过信噪比(SNR)分别设置为30 dB、20 dB和10 dB来模拟不同的噪声条件。

首先，通过数值方法生成混沌系统的时间序列数据，对每个生成的混沌信号，通过加入不同强度的高斯白噪声来模拟不同的噪声环境，噪声强度通过调整噪声的方差来控制。实施了两种去噪策略，一种是采用新提出的优化策略，首先利用和谐相似块法识别信号中的相似块，随后通过DTW重构技术对这些块进行优化重构；另一种则是作为对照的传统协同滤波去噪方法。结果见 表6 。

2) 与其他去噪方法的比较分析

除了与传统的协同滤波进行比较以外，还与几种常用的去噪技术，包括小波去噪、高斯滤波和经验模态分解(EMD)去噪方法进行了比较。为全面评估各去噪技术的表现，我们采用了三个主要评价指标：决定系数(R²)、均方误差(MSE)和均方根误差(RMSE)。实验结果如 表7 所示。

从对比的结果来看，随着信号点个数的增加，本文提出的改进自适应协同滤波方法的去噪效果有明显的提升，并且比经验模态分解方法更好，比高斯滤波和小波滤波方法稍差。但是应该注意到，协同滤波算法的原理比较简单，计算量相对比较小，在滤波精度要求不太高的条件下，具有广泛的应用场景。总体而言，混沌信号的长度越长，自适应协同滤波算法的去噪效果越显著，并且随着信号长度的增加，去噪效果呈现出明显的提升趋势。这些发现不仅为混沌信号去噪领域提供了新的研究视角，也为未来的研究方向和应用实践指明了方向。