1. 项目概述：当雷达信号处理遇上脑启发计算

在自动驾驶、无人机导航和智能安防这些前沿领域，雷达传感器扮演着“眼睛”的角色，它不惧风雨、穿透力强，能提供目标的速度、距离和方位信息。然而，随之而来的海量雷达频谱数据，却给后端处理系统带来了巨大压力。传统的深度神经网络（DNN）虽然识别精度高，但其动辄数百万的参数、复杂的反向传播训练过程，对计算资源和能耗的“胃口”实在太大，让许多追求实时、低功耗的边缘设备望而却步。这就好比让一台小巧的嵌入式单片机去运行一个大型的图像识别模型，结果往往是算力捉襟见肘，功耗直线飙升，响应速度也跟不上。

正是在这种“高精度”与“低功耗”的矛盾中，一种名为“超维计算”（Hyperdimensional Computing, HDC）的脑启发计算范式进入了我们的视野。它不像DNN那样依赖精密的浮点矩阵运算，而是另辟蹊径，用成千上万个维度的随机向量（称为超向量）来表示信息。这种方法的魅力在于其惊人的计算效率和鲁棒性。想象一下，你不是用一串复杂的数字序列来记住一个概念，而是用一个由数千个随机正负1组成的、独一无二的“高维指纹”来代表它。运算过程也变成了简单的元素级加减乘除和比较。最近，一项名为“DynHD”的研究工作，将这种范式与雷达频谱分类任务深度结合，通过引入动态编码机制，让HDC在处理复杂雷达数据时表现出了媲美DNN的精度，同时在训练速度和推理效率上实现了数量级的提升。这为在资源受限的边缘设备上实现实时、高效的雷达信号理解，打开了一扇新的大门。

2. 超维计算（HDC）核心原理拆解：从“高维指纹”到“脑启发运算”

要理解DynHD的巧妙之处，我们必须先吃透HDC这套独特的世界观和“语法”。它和我们熟悉的神经网络有本质区别，其核心思想源于一个神经科学启发：大脑中的信息可能是以高维、分布式、全息的方式存储和处理的。

2.1 超向量：信息的“高维身份证”

HDC的基本数据单元是“超向量”（Hypervector）。这是一个维度通常在1000到10000之间的随机向量，每个元素的值通常是+1或-1（二值化），或者是+1, 0, -1（三值化）。关键特性在于，随机生成的不同超向量之间几乎是正交的，即它们的点积（相似度）接近于零。这就像在万维空间中，随机指向的两个箭头几乎是垂直的。而一个超向量与自身的相似度则为1（完全一致）。

注意 ：这里的“高维”是核心。维度足够高（通常>1000）才能保证随机向量之间近似正交的概率极高，这是后续所有操作具备鲁棒性的数学基础。维度太低，向量间容易产生非预期的相关性，导致系统混乱。

2.2 HDC的四大基础运算：绑定、捆绑、置换与相似性

HDC通过四种简单的运算来构建复杂的表示和逻辑：

1. 绑定（Binding，符号 ⊙） ：表示“关联”或“组合”。操作是元素级的乘法（XOR的变体）。例如，将代表“红色”的超向量 H_red 和代表“苹果”的超向量 H_apple 绑定，得到 H_red_apple = H_red ⊙ H_apple 。关键性质是，绑定后的新向量与原来的两个向量都不相似（点积≈0），但它唯一地编码了这个组合。这类似于将两个概念“加密”成一个新的、不可逆推的令牌。

2. 捆绑（Bundling，符号 ⊕） ：表示“集合”或“叠加”。操作是元素级的加法。例如，将 H_apple 、 H_banana 、 H_orange 捆绑，得到 H_fruit_set = H_apple ⊕ H_banana ⊕ H_orange 。捆绑后的向量与集合中的每个成员向量都有一定的相似性（点积>0）。这就像把多个声音混音，你仍能从混合音中识别出每个原始音的成分。

3. 置换（Permutation，符号 ρ） ：表示“顺序”或“位置”。操作通常是一个循环移位。例如， ρ(H_apple) 表示“第一个位置的苹果”， ρ(ρ(H_apple)) 表示“第二个位置的苹果”。置换后的向量与原向量几乎正交。这为处理时序数据或序列提供了可能。

4. 相似性（Similarity，符号 δ） ：用于比较和检索。最常用的是余弦相似度或归一化点积。通过计算一个查询超向量与所有已存储的类别超向量的相似度，并取最大值，即可完成分类决策。

2.3 传统HDC编码器的局限：静态的“词汇表”

如何将实际的雷达数据（一连串的浮点数）转换成超向量？传统方法主要有两种：

• 基于记录的编码器（Record-based） ：为数据的每个特征维度（如频谱的每个频率点）分配一个固定的“地址”超向量，为每个可能的特征值（量化后的）分配一个“值”超向量。编码时，将每个特征对应的地址和值向量绑定，再把所有特征的结果捆绑起来。这就像用一本固定的密码本（codebook）来翻译数据。
• N-gram编码器 ：常用于序列数据，为每个特征值分配一个值向量，然后根据该特征在序列中的位置进行相应次数的置换，最后捆绑所有结果。

这两种方法的共同问题是： 编码矩阵（即那些“地址”向量或随机矩阵）是随机生成且固定不变的 。它就像一个预先定义好的、与具体任务无关的“通用词汇表”。对于结构简单、特征明显的任务（如识别简单的手写字母），这个词汇表可能够用。但面对雷达频谱这种高维、复杂、信噪比多变的数据，这个静态词汇表提取特征的能力就非常有限了，导致模型的学习能力和最终分类精度遇到瓶颈。

3. DynHD动态编码机制深度解析：让编码器“学会”雷达的语言

DynHD的核心创新，正是针对上述“静态词汇表”的痛点，提出了一种 动态、可学习的编码机制 。它的目标不是替换HDC的整个框架，而是优化其中最前端、也是最关键的一步——如何将原始数据映射到高维空间。

3.1 基石：非线性HDC编码与随机傅里叶特征

在介绍动态性之前，DynHD采用了一个比传统绑定/捆绑更强大的编码基础：基于随机傅里叶特征（Random Fourier Features, RFF）的非线性编码。其灵感来源于核方法（Kernel Methods），目的是隐式地实现一个径向基函数（RBF）核的近似。

具体操作简洁而巧妙：对于一个d维的输入数据点 x ，其对应的D维超向量 H_x 通过以下方式生成： H_x = cos(ω * x + b) 其中， ω 是一个 D x d 的随机矩阵，每个元素从标准正态分布 N(0, 1) 中采样； b 是一个D维随机向量，每个元素从 [0, 2π] 均匀分布中采样。 cos 是元素级的余弦函数。

实操心得 ：这里的 cos 非线性变换至关重要。它使得编码后的超向量之间的相似度 δ(H_x, H_y) 近似等于 exp(-||x - y||^2 / 2) ，即一个高斯RBF核函数。这意味着，在原始数据空间中相似的两个点，映射到超维空间后依然保持高相似度，反之亦然。这为后续的线性分类器（HDC的捆绑操作本质上是线性的）处理非线性可分的雷达数据提供了理论保障。

3.2 DynHD的动态学习引擎：可训练的编码矩阵变换网络

静态编码的 ω 矩阵是随机的、任务无关的。DynHD的想法是： 我们能否学习一个变换，将这个随机矩阵“扭曲”成一个对当前雷达分类任务更有利的矩阵？

DynHD引入了一个轻量级的编码-解码器结构的神经网络（文中称为变换网络），来实现这一目标。其架构和工作流程如下：

1. 初始化 ：首先，我们仍然随机生成一个初始的编码矩阵 ω 和偏置向量 b 。
2. 变换网络 ：这是一个小型全连接网络，输入和输出维度都是超维度D，中间有一个256维的隐藏层。它的参数是可训练的。我们将随机矩阵 ω 的每一行（对应一个超维度）作为一个输入，送入这个网络。
3. 生成动态编码矩阵 ：变换网络输出一个变换后的矩阵 ω' 。这个 ω' 将替代原始的 ω ，用于上述 H_x = cos(ω' * x + b) 的编码过程。此时，编码过程从 H_x = cos(ω * x + b) 变成了 H_x = cos(f_θ(ω) * x + b) ，其中 f_θ 就是我们的变换网络。
4. 交替训练 ：训练过程分为两个交错进行的阶段：
   • 阶段一（固定编码器，训练分类器） ：使用当前的 ω' 编码一批训练数据，得到超向量，然后使用OnlineHD等自适应算法更新类别超向量（分类器）。计算损失（如分类错误）。
   • 阶段二（固定分类器，训练编码器） ：将阶段一的损失，通过分类器反向传播到编码后的超向量，再进一步反向传播到编码矩阵 ω' ， 最终传播到变换网络 f_θ 的参数 。更新变换网络的参数。
   • 更新几次变换网络后，用更新后的网络重新处理原始的随机 ω ，生成新的、更优的 ω' ，然后回到阶段一。

这个过程的关键在于，变换网络学习的是 如何生成一个好的编码矩阵的映射规则 ，而不是直接记忆一个固定的矩阵。这使得模型能够适应数据分布，并且由于变换网络本身很小，引入的开销非常有限。

3.3 与NeuralHD的对比：从“进化”到“学习”

在DynHD之前，NeuralHD也尝试了动态编码。它的思路是“进化”：定期评估每个超维度对分类的重要性（例如，通过计算该类超向量在不同维度上的方差），然后将重要性低的维度标记为“无效”，并用新的随机值替换编码矩阵中对应的行，同时将分类器中该类超向量的对应维度置零。

注意事项 ：NeuralHD的“动态”是一种启发式的、离散的替换。它像是一个维度的“自然选择”过程，淘汰弱的，引入新的随机突变。而DynHD的“动态”是连续的、基于梯度的学习。它通过一个可微分的网络，系统性地调整整个编码空间的“形状”，使其更贴合雷达数据的流形结构。从结果看，DynHD这种“学习”的方式，通常能比“进化”的方式找到更优的编码策略，从而获得更高的精度。

4. 从理论到实践：DynHD在雷达频谱分类上的实现要点

理解了DynHD的原理，我们来看看如何将其应用于实际的雷达频谱分类任务。这里以公开的MSTAR（合成孔径雷达目标识别）和CARRADA（汽车雷达数据集）为例，拆解关键步骤。

4.1 数据预处理与特征工程

雷达原始数据（如ADC采样数据）不能直接送入HDC模型。通常需要经过以下预处理：

1. 时频分析 ：对于调频连续波（FMCW）雷达如CARRADA数据，需要通过2D-FFT（距离-多普勒处理）将原始中频信号转换为距离-多普勒谱（RD谱）。对于合成孔径雷达（SAR）图像如MSTAR数据，本身已是2D图像。
2. 标准化与向量化 ：将得到的频谱或图像数据（可能是2D矩阵）展平为一维向量。进行标准化处理（如减均值、除方差），使输入数据分布在一个相对稳定的范围内。
3. 降维（可选） ：如果原始特征维度（d）非常高（例如数万），而超维度（D）只有几千，直接编码可能信息压缩过度。可以考虑使用PCA等线性方法进行适度降维，保留主要能量成分。

实操心得 ：预处理步骤对HDC模型性能的影响比DNN更敏感。因为HDC的编码相对线性（虽然有cos非线性，但核心是线性投影），特征的质量和区分度至关重要。在CARRADA数据上，直接使用RD谱的幅度值作为特征效果就不错；而在MSTAR图像上，可能需要先提取一些纹理特征（如HOG）再输入，效果会比直接用像素更好。

4.2 模型构建与训练流程

以下是基于PyTorch框架的一个简化实现逻辑：

import torch
import torch.nn as nn
import torch.nn.functional as F
import numpy as np

class DynamicEncoder(nn.Module):
    """DynHD的可学习编码器变换网络"""
    def __init__(self, hyperdim_D):
        super().__init__()
        self.fc1 = nn.Linear(hyperdim_D, 256)
        self.fc2 = nn.Linear(256, hyperdim_D)
        self.relu = nn.ReLU()

    def forward(self, base_omega):
        # base_omega: [D, d]
        # 我们将每一行（一个D维向量）独立地通过网络变换
        # 这里为了效率，可以重塑张量，但概念上是逐行处理
        transformed = self.relu(self.fc1(base_omega))
        transformed = self.fc2(transformed)
        # 输出变换后的omega矩阵，维度不变 [D, d]
        return transformed

class DynHDClassifier:
    """DynHD分类器，包含编码器和类别超向量"""
    def __init__(self, input_dim, hyperdim_D, num_classes):
        self.D = hyperdim_D
        self.d = input_dim
        self.num_classes = num_classes

        # 初始化随机基础矩阵和偏置
        self.base_omega = torch.randn(self.D, self.d)  # 固定，不训练
        self.base_bias = torch.rand(self.D) * 2 * torch.pi  # 固定，不训练

        # 动态编码器网络
        self.encoder_net = DynamicEncoder(self.D)

        # 初始化可训练的编码矩阵（由网络生成）
        with torch.no_grad():
            self.current_omega = self.encoder_net(self.base_omega)

        # 初始化类别超向量为零向量
        self.class_hypervectors = torch.zeros(num_classes, self.D)

        # 优化器只优化编码器网络参数
        self.encoder_optimizer = torch.optim.Adam(self.encoder_net.parameters(), lr=1e-3)

    def encode(self, x):
        """使用当前动态编码矩阵将输入数据编码为超向量"""
        # x: [batch_size, d]
        # 计算 H = cos(omega * x + b)
        projection = F.linear(x, self.current_omega.T, self.base_bias) # [batch_size, D]
        hypervectors = torch.cos(projection) # [batch_size, D]
        return hypervectors

    def train_step(self, data_batch, label_batch, classifier_lr=0.01):
        """
        一个训练步骤，包含分类器更新和编码器更新（每N步一次）
        """
        # 1. 编码数据
        H = self.encode(data_batch) # [batch_size, D]

        # 2. OnlineHD风格更新分类器（简单示例，非完整OnlineHD）
        similarities = F.cosine_similarity(H.unsqueeze(1), self.class_hypervectors.unsqueeze(0), dim=2) # [batch_size, num_classes]
        preds = similarities.argmax(dim=1)

        for i in range(len(data_batch)):
            h = H[i]
            true_label = label_batch[i]
            pred_label = preds[i]

            # 更新正确类别的超向量：加上（1 - 相似度）* 样本向量
            sim_correct = similarities[i, true_label]
            self.class_hypervectors[true_label] += classifier_lr * (1 - sim_correct) * h

            # 如果预测错误，惩罚预测错误的类别：减去（1 - 相似度）* 样本向量
            if pred_label != true_label:
                sim_wrong = similarities[i, pred_label]
                self.class_hypervectors[pred_label] -= classifier_lr * (1 - sim_wrong) * h

        # 3. 计算分类损失（用于更新编码器）
        loss = F.cross_entropy(similarities * 10, label_batch) # 缩放相似度以匹配交叉熵输入范围

        # 4. 每隔一定批次（例如10个批次），更新一次编码器网络
        if self.update_encoder_counter % 10 == 0:
            self.encoder_optimizer.zero_grad()
            loss.backward() # 梯度会通过encode函数流回current_omega，进而流回encoder_net
            self.encoder_optimizer.step()

            # 用更新后的网络重新生成编码矩阵
            with torch.no_grad():
                self.current_omega = self.encoder_net(self.base_omega)

        self.update_encoder_counter += 1
        return loss.item()

4.3 超参数选择与调优经验

• 超维度D ：这是HDC最重要的参数。论文中测试了1k和10k。经验是： 精度要求高、数据复杂，选大D（如4k, 10k）；追求极致速度和低内存，选小D（如1k） 。对于MSTAR/CARRADA这类任务，1k维度已能获得不错效果，4k或10k能逼近DNN精度。
• 编码器网络结构 ：DynHD原文使用了一个两层的MLP（D -> 256 -> D）。这是一个很好的起点。如果输入特征维度d很大，可以考虑在第一层前加入一个从d到某个中间维度的映射，但会增加参数。 保持这个网络尽可能小是保证效率的关键 。
• 学习率与更新频率 ：分类器更新（ classifier_lr ）通常较大（如0.01-0.1），因为是对超向量直接做加法。编码器网络的学习率（ encoder_optimizer 的lr）要小得多（如1e-4到1e-3），因为它是在学习一个缓慢变化的映射。编码器网络的更新频率（如每10个批次）也需要平衡，太频繁则训练不稳定，太慢则动态性不足。
• 量化与部署 ：HDC模型天生对量化友好。因为核心运算（编码时的乘加、相似度计算的点积）都是整数或定点数操作。在部署到嵌入式设备（如ARM Cortex-M系列）时，可以将编码矩阵 ω' 、偏置 b 和类别超向量全部量化为8位甚至4位整数，在推理时使用整数运算单元，能极大降低功耗和延迟。DynHD论文也验证了其在低精度下的鲁棒性。

5. 性能评估与对比分析：为何DynHD是边缘设备的优选？

根据论文中的实验结果，我们可以从几个维度来审视DynHD的优势。

5.1 精度：不逊色于DNN，远超传统机器学习

模型	MSTAR数据集准确率	CARRADA数据集准确率	备注
DNN (3层MLP)	~95%	~92%	传统深度学习基线，参数量大
SVM (RBF核)	~88%	~85%	经典机器学习方法，推理慢
决策树	~82%	~80%	简单，但精度有限
传统HDC (RecordHD)	~85%	~83%	静态编码，性能瓶颈明显
OnlineHD (10k维)	~93%	~90%	自适应训练，但编码静态
NeuralHD (1k维)	~94%	~91%	动态维度进化
DynHD (1k维)	~95%	~92%	动态可学习编码
DynHD (4k维)	~97%	~93%	动态可学习编码

从上表可以看出， DynHD在仅使用1k维度时，就能达到与DNN相当的精度；当维度提升到4k时，其精度甚至能超越DNN 。这充分证明了动态编码机制有效提升了HDC的特征提取能力，使其能够捕捉雷达数据中更细微、更复杂的模式。

5.2 效率：训练与推理的“降维打击”

这是HDC范式，尤其是DynHD，最闪耀的优势所在。

• 训练速度 ：在笔记本电脑CPU上，DynHD-1k的训练时间相比DNN可缩短 8倍以上 。这是因为：

  1. 无反向传播 ：HDC分类器的更新（OnlineHD）是简单的向量加减法，复杂度为O(D)，远低于DNN的梯度反向传播（O(参数数量)）。
  2. 快速收敛 ：HDC模型通常只需遍历训练集少数几次（甚至一次）就能达到不错的精度，而DNN需要数十上百个epoch。
  3. 轻量编码器更新 ：DynHD的编码器网络虽然需要反向传播，但其网络非常小（如D->256->D），且更新频率低（每10个批次），开销远小于训练整个DNN。

• 推理速度与功耗 ：在树莓派这类嵌入式CPU上，优势更为明显。

  • 推理延迟 ：DynHD-1k的推理速度可比DNN快 100倍 。推理过程就是一次编码（矩阵乘加+cos）和一次与所有类别向量的点积比较（O(D * 类别数)），计算极其简单。
  • 功耗与能量 ：论文中测量了训练过程中的总能耗。由于训练时间大幅缩短且单次操作计算量小，DynHD的总能耗可比DNN低 1到2个数量级 。这对于电池供电的边缘设备至关重要。

5.3 鲁棒性：对量化与噪声的天然抵抗力

边缘设备常使用低精度计算（如INT8, INT4）来节省内存和能耗。DNN模型在低精度量化下精度往往骤降。

模型精度	DNN (MSTAR)	DynHD-4k (MSTAR)
FP32 (全精度)	95.1%	96.8%
INT8 (8位整型)	92.5%	96.2%
INT4 (4位整型)	75.3%	94.1%

DynHD（及HDC模型）对量化表现出惊人的鲁棒性。即使在4位整数量化下，精度损失也微乎其微。这是因为HDC的核心运算（绑定是XOR，捆绑是加法）本身对数值精度不敏感，且高维空间的表示本身具有冗余性和容错性。同样，对于传感器数据中常见的噪声，HDC模型也比DNN更具韧性。

6. 常见问题与实战避坑指南

在实际尝试复现或应用DynHD时，你可能会遇到以下问题：

问题1：超维度D到底设多大？是不是越大越好？

• 现象 ：盲目设置D=10000，导致模型内存占用大，推理速度提升不明显。
• 分析与解决 ：D的增加会线性增加内存（存储类别向量）和计算量（编码和相似度计算）。存在一个“收益递减”点。 建议从1k开始，逐步增加（2k, 4k, 10k），观察验证集精度变化 。对于许多任务，4k可能是一个甜点。同时，考虑目标硬件平台的缓存大小，过大的D可能导致缓存频繁失效，反而降低速度。

问题2：动态编码器训练不稳定，精度震荡甚至下降。

• 现象 ：启用DynHD的动态编码训练后，损失函数波动剧烈。
• 分析与解决 ：
  1. 检查学习率 ：编码器网络的学习率可能过高。尝试从更小的值开始（如1e-5），并配合学习率调度器（如ReduceLROnPlateau）。
  2. 调整更新频率 ：不要每个批次都更新编码器网络。论文中建议每10个分类器训练批次更新一次。可以尝试更长的间隔（如20或50批次），让分类器在当前编码下“充分学习”一段时间。
  3. 冻结初期训练 ：在训练初期（如前几个epoch），先冻结编码器网络，只用静态编码训练分类器，待分类器初步稳定后，再解冻编码器进行联合微调。

问题3：在嵌入式设备上部署，编码时的cos非线性函数计算较慢。

• 现象 ：虽然HDC运算简单，但 cos 函数在缺乏硬件浮点单元的MCU上计算开销较大。
• 分析与解决 ：
  1. 查表法 ：预先计算一个 cos 函数查找表。由于输入是 ω*x + b ，其值域有一定范围，可以量化到固定点数并索引查找表。
  2. 多项式近似 ：使用低阶多项式（如泰勒展开前几项）在特定区间内近似 cos 函数。
  3. 简化编码 ：对于极端资源受限场景，可以回归研究更简单的编码方式（如基于量化的绑定/捆绑），虽然可能损失一些精度，但能换来极致的速度。DynHD的动态思想也可以尝试迁移到这些简化编码上。

问题4：如何处理多模态或时序雷达数据？

• 现象 ：雷达数据可能包含距离、速度、角度等多维信息，或者是连续的时间帧序列。
• 分析与解决 ：
  • 多模态融合 ：为每种模态数据（如距离谱、多普勒谱）分别构建一个DynHD编码器，生成各自的超向量 H_range 和 H_doppler 。然后通过 绑定 操作将它们融合： H_fused = H_range ⊙ H_doppler 。绑定操作能很好地表示不同模态特征的联合出现。
  • 时序处理 ：对于连续帧，可以将每一帧编码为一个超向量 H_t 。然后利用 置换 操作来表示时序顺序，例如， H_sequence = ρ^1(H_t) ⊕ ρ^2(H_{t-1}) ⊕ ρ^3(H_{t-2}) ，将最近几帧的信息捆绑在一起，置换赋予了它们时间位置信息。这为处理雷达目标跟踪等时序任务提供了可能。

DynHD为我们展示了一条清晰的路径：通过将脑启发的超维计算与可学习的动态编码相结合，我们能够在资源严格的边缘环境下，实现复杂雷达信号的高效、精准理解。它不仅仅是一个算法，更是一种设计范式——在追求性能的同时，始终将能效和可行性放在核心位置。在实际项目中，当你面临“既要精度又要速度还要低功耗”的三角难题时，不妨将HDC及其变体如DynHD纳入你的候选方案清单，它可能会带来意想不到的突破。