1. 项目概述为什么我们要死磕高频SSVEP脑机接口如果你接触过脑机接口BCI大概率听说过稳态视觉诱发电位SSVEP技术。它的原理很直观让人盯着一个以特定频率闪烁的视觉刺激大脑的视觉皮层会产生一个与刺激频率同步的振荡电信号也就是SSVEP。通过检测脑电EEG中哪个频率的响应最强就能判断用户在看哪个目标从而实现“意念”控制。传统的SSVEP-BCI大多工作在低频比如8-15Hz或中频15-30Hz段因为在这些频段SSVEP的响应幅度大容易检测。但低频刺激有个绕不开的“硬伤”视觉疲劳感强长时间使用容易引发不适甚至对部分光敏性癫痫患者存在风险。此外低频段的脑电信号容易与大脑自身的α节律8-13Hz等自发活动重叠产生干扰。于是把工作频率推到30Hz以上的“高频SSVEP-BCI”就成了一个很有吸引力的研究方向。高频刺激更接近日常荧光灯的闪烁频率视觉感受更舒适也更安全。然而理想很丰满现实很骨感。高频SSVEP的振幅会随着频率升高而急剧衰减信号变得非常微弱。同时这个频段又恰恰是肌电EMG伪影的“重灾区”一个不经意的眨眼、皱眉产生的噪声就可能把微弱的SSVEP信号完全淹没。所以很长一段时间里高频SSVEP-BCI的性能准确率和信息传输速率一直上不去实用化门槛很高。我和团队在几年前也深陷这个泥潭。当时的主流方法是“看功率”对EEG信号做频谱分析看哪个刺激频率对应的功率谱峰值最突出。但在高频段这个方法就像在嘈杂的菜市场里听一根针掉地上的声音非常吃力。我们意识到必须换个思路——SSVEP不仅是频率跟随它还是相位锁定的。也就是说大脑响应不仅频率和刺激一致其波形的“起伏时刻”也与闪烁的“亮-灭”时刻保持固定的时间关系。这个锁定的“波形形状”包含了基波和各次谐波之间独特的比例和相位关系是比单纯的功率更稳定、更具个体特异性的特征。如何把这个“形状”特征高效、鲁棒地提取出来就成了破局的关键。这就引出了我们这次项目的核心相位同步梳状滤波器。它不是简单地在频域找峰值而是在时域直接“雕刻”出与刺激同步的响应波形模板从而在强噪声中牢牢抓住我们想要的信号。2. 系统核心设计硬件与算法的“双剑合璧”一个高性能的BCI系统绝不是算法单方面秀操作就能成的。它需要前端刺激呈现的“绝对稳定”和后端信号处理的“精准狙击”紧密配合。我们的系统设计正是围绕这两个核心展开的。2.1 硬件基石BCI Appliance与精准的刺激控制工欲善其事必先利其器。高频SSVEP对刺激的时序精度要求极高微秒级的抖动都可能导致相位信息错乱让后续的同步处理功亏一篑。用普通电脑屏幕通过软件控制闪烁基本不可行。主流显示器的刷新率通常60Hz或144Hz限制了可呈现频率的精度和范围且软件定时会受操作系统调度影响产生无法预测的延迟和抖动。为此我们专门设计并使用了BCI Appliance这套硬件设备。你可以把它理解为一个为BCI定制的“刺激呈现专用机”。它的核心是一个由微控制器MCU直接驱动的LED阵列作为背光。LCD屏幕覆盖在LED之上用于显示静态的字符或图形内容。这样做的好处是实现了“动静分离”动态闪烁硬件保证每个刺激区域的闪烁频率、占空比都由MCU的硬件定时器精确控制独立可调稳定无抖动。我们实验中使用的是30-39Hz共10个频率点每个都能被精准生成。静态内容软件渲染LCD屏幕负责显示“A”、“B”、“删除”等字符或图标这些内容可以根据任务动态变化而底下的LED背光则按照设定好的频率持续、稳定地闪烁。注意这种硬件分离设计是关键。它确保了诱发SSVEP的物理光源LED其时间特性是绝对可靠、可追溯的为后续算法进行严格的相位对齐提供了可能。如果闪烁是由软件渲染图形比如一个白色方块明暗交替产生的其定时精度会大打折扣。为了给算法提供精确的“时间锚点”我们在每个刺激LED旁都固定了一个光电二极管。它的信号被同步记录到EEG放大器的辅助通道中。这样我们从EEG数据中不仅能读到大脑的电信号还能精确知道每一刻每个刺激方块是“亮”还是“灭”以及状态切换的精确时刻。这个“刺激标记通道”是后续所有相位同步处理的基石。2.2 算法核心相位同步梳状滤波器原理拆解有了硬件提供的精准时间戳我们就可以施展算法的魔法了。传统方法可以看作是在“找声音最大的频率”而我们的方法更像是在“听一首熟悉的歌并从嘈杂环境中把它识别出来”。这首歌的“乐谱”就是个体特异性的SSVEP响应波形模板。1. 空间滤波先给信号“降噪”和“聚焦”原始EEG信号来自头皮多个点位我们用了枕叶和顶叶的8个导联里面混有大量与视觉任务无关的脑电活动和其他伪影。第一步我们使用公共空间模式CSP滤波。简单来说CSP会找到一种对多通道EEG信号的线性组合方式使得在这个新组合出的“虚拟通道”信号上目标刺激比如看35Hz时的信号方差最大而非目标刺激看其他频率时的信号方差最小。这相当于把信号的能量聚焦到与任务最相关的成分上同时抑制了无关的脑电背景噪声为后续处理提供了一个信噪比更高的单通道信号。2. 时域模板构建梳状滤波的实质这是最核心的一步。假设用户正在注视一个以频率F例如35Hz周期T≈28.6ms闪烁的刺激。我们从CSP滤波后的信号中截取一段3秒长的数据缓冲区。关键操作根据光电二极管记录的信号我们精确找到这段3秒数据内目标刺激每一次“亮起”上升沿的时刻点k1, k2, ..., km。对齐与平均以每个“亮起”时刻为起点我们从CSP信号中切出长度为Lp例如1秒的一小段信号。这样我们就得到了m个信号片段每个片段都理论上包含了从刺激开始后1秒内的大脑响应。由于刺激是周期性的如果SSVEP响应是相位锁定的那么这m个片段中的响应波形应该是高度相似的。形成模板将这m个片段进行时间对齐后直接平均。随机噪声在多次平均中会相互抵消而与刺激锁定的SSVEP响应波形则会得到增强最终浮现出来。这个平均后的波形就是我们针对频率F提取到的时域响应模板p^F。这个过程在信号处理上等效于一个梳状滤波器。想象一把梳子它的齿间距正好是刺激周期T。当信号通过这个滤波器时所有与刺激周期T同步的成分即相位锁定的SSVEP及其谐波会在每个齿的位置延迟T、2T、3T...被同相叠加从而增强。而那些不与刺激同步的噪声成分则会在不同延迟下以随机相位叠加相互抵消。这个“梳子”的齿就是由硬件同步信号精确设定的延迟点。实操心得Lp模板长度的选择是个权衡。太短如0.5秒包含的周期数少平均降噪效果有限太长如2秒则可能包含更多非稳态的脑电波动稀释了锁定的响应。我们通过实验发现对于30-39Hz的刺激周期约25-33ms1秒的长度包含约30-40个刺激周期能在降噪和稳定性间取得较好平衡。同时用于平均的缓冲区长度Lb3秒要足够长以确保有足够多的刺激周期m值够大用于平均但也不能太长导致系统响应延迟过高。3. 特征提取与分类从“像不像”到“是不是”在训练校准阶段我们让用户依次注视每个已知频率的刺激为每个频率F都构建一个“黄金标准”模板p^F_T。 在线操作时系统实时地对最新的3秒EEG数据用同样的方法CSP滤波 基于各频率闪烁时刻的对齐平均为当前屏幕上所有8个刺激频率分别生成8个“候选模板”p^F_C。 接下来就是“找对象”环节计算每个候选模板p^F_C与其对应的黄金标准模板p^F_T之间的皮尔逊相关系数r。这个相关系数衡量了实时波形与标准波形在“形状”上的相似程度值越接近1说明形状越像。 最后一个简单的朴素贝叶斯分类器会根据这个相关系数r属于“目标”或“非目标”类别的概率来判断用户当前正在注视哪个频率。如果只有一个频率的相关系数被判定为“目标”系统就输出该频率对应的命令。3. 从理论到实践系统搭建与实验全流程实录纸上谈兵终觉浅绝知此事要躬行。下面我详细拆解我们是如何把上述设计落地并拿到那些漂亮数据的。3.1 实验设置与数据采集的魔鬼细节硬件清单与参数BCI Appliance8个独立LED刺激场每个尺寸4.6cm x 3.2cm观看距离约100cm约2.6度视角。LED由STM32系列MCU驱动频率控制精度达0.01Hz。8个光电二极管监测闪烁。EEG采集TMSI Porti7放大器采样率512Hz。电极帽配备8个盐水电极位置覆盖PO8, O2, O1, PO7, POz, PO3, Pz, Cz国际10-20系统参考电极为PO8地电极在左锁骨。这个导联布局集中覆盖了视觉诱发电位产生最强的枕叶和顶叶区域。实验对象15名健康受试者20-30岁5女10男。所有实验均通过伦理审查。注意事项电极阻抗一定要尽可能低尤其是使用盐水电极时在实验开始前需要花时间确保每个电极与头皮接触良好阻抗最好保持在10kΩ以下。高频SSVEP信号本就微弱高阻抗会引入额外的噪声严重影响信号质量。这是我们早期实验踩过的大坑。实验流程三步走 整个实验分为校准和在线拼写两个阶段全程人机交互界面由BCI Appliance的LCD屏幕呈现。第一步校准约8分钟这是为每个用户“量身定做”模型的关键步骤。屏幕上8个方块同时以不同的频率闪烁从30, 31, ..., 39Hz这10个频率中选8个。屏幕顶部状态栏会提示用户当前需要注视哪个方块用“*”标出。每个试次持续5秒注视 1秒休息。每个频率都会作为目标频率出现8次总共80个试次。在这5秒内系统记录EEG和光电二极管信号用于后续计算每个频率的“黄金标准”模板和训练分类器。第二步在线拼写任务校准完成后立即进行在线测试检验系统实时性能。简单七字母拼写器屏幕上只显示当前要拼写的单词如“tablica”所含的7个字母按顺序排列在第1到第7个刺激块第8个块是“删除”命令。用户只需直接注视想输入的字母即可。任务是用这个界面拼出5个7字母的波兰语单词。全字母表拼写器两级菜单这是一个更接近实际应用的界面。第一级7个刺激块显示字母分组如“ABC”“DEF”…第8块是“删除”。用户先注视目标字母所在的分组。选择分组后进入第二级该分组内的字母及一个“返回”命令显示在8个块上用户再注视目标字母完成输入。任务是用这个两级菜单拼出5个长度从3到23个字母不等的短语。实操心得两级菜单的设计虽然增加了认知负荷但它是将有限刺激目标8个扩展为无限指令集整个字母表的经典方法。实验中发现从简单任务切换到复杂任务时用户需要短暂的适应期。在设计这类实验时务必给用户清晰的界面指引和足够的练习试次否则错误率会显著上升。3.2 信号处理流水线代码级解析下面我结合伪代码和参数选择还原在线运行时的核心处理循环。假设我们已从校准数据中得到了1) CSP空间滤波器权重向量w2) 8个频率对应的黄金标准模板ref_templates[8]3) 分类器参数目标/非目标类别的相关系数概率分布。# 伪代码在线分类核心循环 import numpy as np from scipy import signal, stats buffer_length 3 * fs # 3秒缓冲区fs512Hz - 1536个采样点 template_length 1 * fs # 1秒模板长度 - 512点 update_interval 0.2 * fs # 每0.2秒更新一次 - 102点 # 初始化一个先进先出FIFO缓冲区 eeg_buffer np.zeros((num_channels, buffer_length)) aux_buffer np.zeros((8, buffer_length)) # 8个光电二极管信号 # 主循环 while True: # 1. 数据填充由采集线程实时填入新数据 # ... 假设新数据已填入缓冲区末尾并进行了FIFO移位 ... # 每0.2秒处理一次 if new_data_count update_interval: # 2. 预处理当前缓冲区数据 current_eeg eeg_buffer[:, -buffer_length:] # 取最新3秒数据 current_aux aux_buffer[:, -buffer_length:] # a. 去趋势减去均值消除慢漂移 current_eeg current_eeg - np.mean(current_eeg, axis1, keepdimsTrue) # b. 陷波滤波滤除50Hz工频干扰根据地区可能是60Hz b, a signal.iirnotch(50.0, 30.0, fs) # 中心频率50HzQ值30 current_eeg signal.filtfilt(b, a, current_eeg, axis1) # 零相位滤波 # c. CSP空间滤波投影到最优空间方向 # w 是从校准数据计算出的CSP滤波器1 x num_channels 向量 csp_signal np.dot(w, current_eeg) # 得到1维时间序列 # 3. 为每个刺激频率生成候选模板 candidate_templates [] correlation_features [] for freq_idx in range(8): # 从辅助通道获取该频率刺激的“亮起”时刻点 aux_ch current_aux[freq_idx, :] # 检测上升沿得到索引列表 onset_indices onset_indices detect_rising_edges(aux_ch) # 自定义函数 # 确保有足够多的周期用于平均例如至少20个 if len(onset_indices) 20: candidate_templates.append(None) correlation_features.append(0.0) continue # 对齐并平均 epochs [] for onset in onset_indices: if onset template_length buffer_length: epoch csp_signal[onset : onset template_length] epochs.append(epoch) if epochs: avg_template np.mean(epochs, axis0) candidate_templates.append(avg_template) # 4. 计算特征与黄金标准模板的相关系数 ref_temp ref_templates[freq_idx] corr, _ stats.pearsonr(ref_temp, avg_template) correlation_features.append(corr) else: candidate_templates.append(None) correlation_features.append(0.0) # 5. 朴素贝叶斯分类 # 假设我们已从校准数据估计出目标类相关系数分布 P(r|T) ~ N(μ_T, σ_T) # 和非目标类分布 P(r|NT) ~ N(μ_NT, σ_NT) # 先验概率 P(T) 1/8, P(NT) 7/8 (因为8个目标中只有1个是真的) posterior_probs [] for r in correlation_features: # 计算似然 likelihood_T norm.pdf(r, mu_T, sigma_T) likelihood_NT norm.pdf(r, mu_NT, sigma_NT) # 计算后验概率忽略归一化分母因比较相对大小 post_T likelihood_T * prior_T post_NT likelihood_NT * prior_NT # 判断为该频率是目标的“证据比” evidence_ratio post_T / post_NT posterior_probs.append(evidence_ratio) # 6. 决策 # 找出证据比最高的频率 max_idx np.argmax(posterior_probs) max_ratio posterior_probs[max_idx] # 设置一个阈值避免噪声误触发。例如要求证据比至少大于10且其他频率的证据比都远低于它。 threshold 10.0 if max_ratio threshold and all(r max_ratio/5 for i, r in enumerate(posterior_probs) if i ! max_idx): selected_command command_mapping[max_idx] # 映射到具体命令如字母‘A’ send_command_to_interface(selected_command) # 提供反馈如让选中的方块变绿1秒 provide_visual_feedback(max_idx) else: # 无法做出可靠决策等待下一个0.2秒周期 pass # 清空特征缓存准备下一次循环 correlation_features.clear()关键参数与选择理由缓冲区长度buffer_length 3s这是一个权衡。更长的缓冲区能平均更多刺激周期提升模板信噪比但会增加系统指令输出延迟。3秒是一个在响应速度和分类可靠性之间取得较好平衡的点也是很多SSVEP-BCI研究的常用值。更新间隔update_interval 0.2s即系统每0.2秒尝试分类一次。这决定了BCI的“反应速度”。0.2秒已经接近人脑视觉反应和EEG采集的生理极限再短则每次分析的数据量太少分类不可靠。模板长度template_length 1s如前所述包含约30-40个高频刺激周期足以通过平均有效抑制噪声同时避免引入过多非稳态活动。CSP滤波注意这里的CSP是在校准阶段用带通滤波后的数据训练出来的。例如对于目标试次我们用中心频率为f_T的带通滤波器我们用了1阶切比雪夫II型过滤EEG计算协方差矩阵R_T对于非目标试次则用中心频率为某个f_NT的带通滤波器过滤同一段EEG计算R_NT。这样训练出的CSP滤波器能最大化目标频率响应同时最小化其他频率干扰的响应效果比用原始宽带信号训练更好。4. 性能优化与结果深度分析我们做到了什么程度实验做完数据到手是骡子是马拉出来溜溜。我们不仅关心平均性能更想深挖每个数字背后的意义。4.1 特征有效性验证为什么“形状”比“功率”更靠谱在校准阶段我们就对提出的“时域模板相关系数”特征和传统的“信噪比SNR”特征进行了正面较量。SNR是SSVEP分析中的经典特征计算刺激频率处的功率与其相邻频率带平均功率的比值。我们对每个受试者为其8个最佳频率分别计算了两种特征的接收者操作特征曲线下面积AUC。AUC可以理解为分类器区分“目标”和“非目标”试次的能力1.0表示完美区分0.5表示和随机猜测没区别。结果非常明确对于所有受试者时域模板相关系数特征的AUC值显著高于SNR特征的AUC值Wilcoxon符号秩检验p0.001。这意味着利用相位锁定波形“形状”信息比单纯利用频域“功率”信息能更清晰、更鲁棒地将目标响应从背景噪声和非目标响应中分离出来。这从原理上证实了我们方法的优越性它更充分地利用了SSVEP响应中谐波间相位关系的个体化信息。4.2 BCI Literacy不是所有人都能“玩得转”15名受试者中有9人60%成功完成了全部两个在线拼写任务。另外6人即使在校准阶段表现尚可也无法在在线任务中稳定控制。这个比例与文献中报道的BCI“文盲率”约20%-30%的用户无法有效控制大体一致。深入分析校准数据我们发现成功者与失败者的差异一目了然成功者在注视目标频率时生成的候选模板p^F_C彼此高度一致且与黄金标准模板p^F_T高度相关相关系数分布ρ^F_T集中在0.8以上。同时非目标频率的模板与目标模板差异明显ρ^F_NT分布集中在低值。两者的分布几乎不重叠。失败者目标频率的模板一致性差相关系数分布ρ^F_T很分散从0.2到0.8都有且与非目标分布ρ^F_NT存在大量重叠。这意味着他们的大脑要么没有产生稳定、相位锁定的高频SSVEP响应要么响应太弱被噪声淹没导致算法无法可靠识别。经验总结BCI表现存在显著的个体差异。对于高频SSVEP这种差异可能更明显。在实际应用部署前进行快速的校准和筛查是必要的。对于无法产生稳定响应的用户可能需要尝试其他刺激范式如中低频SSVEP、P300或运动想象。4.3 性能数据解读简单与复杂界面的博弈对于9位成功的使用者我们得到了以下核心性能指标任务类型平均准确率平均信息传输率 (ITR)平均单指令选择时间7字母拼写器96% ± 4%47 ± 5 bits/min3.6 ± 0.2 秒全字母表拼写器89% ± 9%36 ± 9 bits/min3.9 ± 0.3 秒1. 准确率 (Accuracy) 简单拼写器的准确率极高很多情况下达到100%。这说明在界面直观、选项直接的情况下我们的信号处理和分类算法本身非常可靠。全字母表拼写器的准确率有所下降但仍接近90%这主要归因于界面复杂度和用户错误。两级菜单增加了操作步骤和认知负荷用户可能选错分组或在第二级犹豫导致注视目标不集中。2. 信息传输率 (ITR) ITR是衡量BCI通信效率的综合指标它同时考虑了准确率P、指令数量N和选择每个指令所需的时间T。计算公式为ITR (60/T) * [log₂N P log₂P (1-P)log₂((1-P)/(N-1))]bits/min。 我们的7字母拼写器达到了47 bits/min这在高频SSVEP-BCI领域是当时报道的最高水平之一。全字母表拼写器ITR下降主要原因是单指令选择时间增加和准确率轻微下降。3. 选择时间分析 全字母表拼写器平均每次选择比简单拼写器慢310毫秒。这310毫秒差在哪我们回放了分类器日志发现用户在反馈选中项变绿结束后往往需要一小段时间重新定位并聚焦到下一个目标刺激上。在复杂的、非顺序排列的两级界面中这个“视觉搜索和注意力聚焦”时间显著延长。这310毫秒几乎完全消耗在用户的认知和视觉运动过程中而非信号处理延迟。核心洞见BCI系统的最终性能是“机器性能”和“人机交互性能”的乘积。我们的算法将“机器性能”从EEG到指令的解码速度和准确率推到了一个很高的水平。但“人机交互性能”——即界面设计是否符合用户直觉、是否能最小化用户的认知负荷和搜索时间——同样至关重要甚至可能成为瓶颈。一个糟糕的界面设计可以轻易抹杀先进的算法带来的增益。5. 避坑指南与未来展望从实验室到实用化的思考基于这个项目从设计、实现到测试的全过程我总结了几条关键的实操经验和未来可以探索的方向。5.1 常见问题与排查技巧实录问题1校准效果很好但在线任务一塌糊涂。可能原因A电极阻抗漂移。校准后到在线任务前如果用户有较大动作或出汗电极阻抗可能发生变化导致信号质量下降。排查与解决在线任务开始前快速检查一遍各通道的原始信号。如果发现某通道噪声明显增大如出现规律的50Hz干扰或大幅基线漂移需要暂停实验重新注入导电膏或调整电极位置。养成在关键任务节点检查阻抗的习惯。可能原因B用户疲劳或注意力涣散。高频SSVEP需要用户持续集中注意力注视闪烁目标。疲劳会导致响应幅度下降。排查与解决安排合理的休息间隔。实验任务不宜过长中间应有强制休息时间。在线拼写任务可以设计成小块进行。问题2分类器频繁误触发或长时间不触发。可能原因A分类器阈值设置不当。在线使用的阈值如伪代码中的threshold和与其他证据比的比值是基于校准数据静态设置的可能不适应在线数据的微小变化。排查与解决实现一个简单的在线自适应机制。例如可以持续监控最近几十次决策中被分类为“目标”的相关系数值分布动态微调mu_T和sigma_T或者调整决策阈值。这能提升系统在不同时段、不同用户状态下的鲁棒性。可能原因B刺激标记光电二极管信号不同步或丢失。这是致命错误会导致模板对齐完全错误。排查与解决在线运行时增加一个实时可视化模块同时显示原始EEG和8个辅助通道的信号。确保每个刺激频率的方波清晰可见且与屏幕闪烁严格同步。如果发现某个辅助通道信号异常如一直是高电平或低电平应立即检查硬件连接。问题3个别频率始终表现很差。可能原因该频率与用户自身的脑电节律如谐波或环境干扰如设备开关电源噪声耦合。排查与解决在频率选择阶段校准后我们计算了每个频率的AUC。务必舍弃AUC值最低的那几个频率即使这意味着你只能使用6个或7个频率来控制8目标的界面。可以通过设计“频率分组”或“编码协议”来弥补指令数量的减少。强用一个响应差的频率会拉低整体系统性能。5.2 系统优化与扩展方向个性化频率调优我们的研究使用了30-39Hz的固定频率集。实际上每个用户对高频刺激的响应存在“最佳频点”。未来系统可以在校准阶段在一个更宽的频带如28-42Hz内进行扫描为每个用户挑选出8个响应最强的个性化频率有望进一步提升信噪比和BCI Literacy。动态模板更新当前的“黄金标准”模板来自一次性的校准。但用户的生理状态会变化。可以引入增量学习机制在在线使用过程中将高置信度的正确试次数据以一定的权重融合到原有模板中使模板能缓慢适应当天的用户状态。混合特征与深度学习我们的方法主要利用了时域波形形状。可以尝试将其与经过精细处理的频域特征如特定谐波的功率比、空域特征多通道CSP组合进行融合输入到一个更强大的分类器如SVM或浅层神经网络中。近年来端到端的深度学习模型如CNN、LSTM直接从原始多通道EEG中提取特征并进行分类在SSVEP-BCI上也展现出潜力但需要大量的训练数据和对计算资源的更高要求。更智能的界面设计正如实验结果所示界面效率直接影响ITR。对于字符输入这类任务可以引入单词预测和智能布局如将常用字母放在更容易注视的位置减少用户需要做出的选择次数和视觉搜索时间从而在解码速度不变的情况下大幅提升实际通信速率。回过头看这个项目最大的收获不仅仅是那几个性能数字而是让我们更深刻地认识到一个成功的BCI系统必须是精密硬件、鲁棒算法和人性化交互设计三者的深度融合。相位同步梳状滤波器为我们提供了一把在高频噪声中精准捕捉大脑信号的“镊子”而如何让用户更自然、更轻松地使用这把“镊子”与外界交流则是通向真正实用化道路上我们接下来需要持续探索的课题。