波动为什么会集中到材料边缘，边界怎样改变内部的响应？

边缘很响，可能是整条链在起作用

想象一排相连的声腔。在中间发出声音，每个声腔把振动传给左右邻居。如果向右传递与向左传递的作用不一样，这种小小的偏差经过许多次传递与反射，可能使整条链的振动模式明显偏向一端。

让人意外的是：不必只在边缘放一个特殊缺陷，也不必只留下一个特殊的边界模式。很多原本属于整条链的模式，都可能一起挤向边缘。把首尾接成环，模式与本征频谱又可能大变。

这就是非厄米皮肤效应的一条典型路线。“皮肤”描述的是大量体模式靠近边界的形状；它与普通导体中交流电的趋肤机制、只存在少数拓扑边缘态的情形，需要分别理解。1

本篇始终分清三件事：**如果系统按某个单独模式振动，它在哪里强；从指定位置持续输入声音，实际哪里强；没有输入时，整个系统会衰减还是增长。**这三件事有关，但不是同一个问题。

2026年的声学研究又把这条路线推进一步：只输入一个频率，通过真实时间调制生成新的频率分量，再让不同分量向不同边界聚集。下面先讲懂静态链，再看这项实验究竟新增了什么。3

从左右传递不一样开始

用复数$q_n$表示第$n$个单元的振幅和相位。一个简化模型为：

$t_R$表示左邻居对本格的作用，即向右传递；$t_L$反之。这里先取两者为正实数，$\gamma>0$是均匀衰减率；时间与参数暂用无量纲单位。

若$t_R\ne t_L$，交换输入端和响应端之后，耦合没有保持原样。把整个系统写成$i\dot{\mathbf q}=H\mathbf q$，相应矩阵不满足$H=H^\dagger$，因此称非厄米。损耗也能使有效矩阵非厄米，但只有损耗这个标签，还不足以推出皮肤效应。1

模型以Hatano–Nelson非互易链为基础。1996年的原研究还处理随机势、局域化以及与超导磁通线的映射；本篇用其中的干净有限链机制，不把原论文的全部物理对象压成一排真实声腔。2

实际装置如何实现这种耦合，要另交代。后文的13声腔实验使用麦克风读取某个腔的声音，经放大和调制后，用扬声器注入另一个腔。方向性和时间变化来自真实的电子反馈与外部供能，不能理解为普通连通管道自动具有这条规则。3

开放链的模式，为什么都带着一个指数偏置

先断开首尾，要求链外没有额外单元。令：

代回方程，左右耦合都会变成$\sqrt{t_Rt_L}$。这意味着，在有限开放链上，可以从一个对称链的正弦模式出发，乘上同一个指数因子，得到原链的右本征模：

当$t_R>t_L$时，$r>1$，这些模式偏向右端。不同模式仍有各自的波节与振荡形状，但指数包络共同把它们推向同一侧。取反方向的不对称，聚集方向也会反转。12

本库取31个单元、$t_R=1.5$、$t_L=1$、$\gamma=0.15$实际计算。振幅的指数局域长度为$1/\log r\approx4.93$个格距；这里说的是振幅，强度是振幅平方，不能混用同一个衰减长度。

31个归一化右本征模的强度重心都在第25.5至29.4格之间。计算出的模式满足原始矩阵方程，残差约为$10^{-15}$。这个结果不依赖某个边缘格子的特殊频率，来自整条链持续存在的方向性差别。

为什么接上最后一条连接，会影响整个谱

闭环要求走完一圈后回到同一个振幅。开放链上的指数变换沿一圈会积累$r^N$，通常无法同时满足周期条件。因此不能把开放链的解原封不动塞进闭环。

在周期链上，允许的模式具有$q_n\propto e^{ikn}$的形式，对应：

当$k$绕过一个周期，$E(k)$在复平面上画出一条椭圆。开放链的本征值则落在另一条线段上。这是边界条件参与决定全链解的一个直接例子；实际建立新响应仍有传播、混合与衰减过程。

上面的示例里，开放链所有模式的增长率都是$-0.15$，闭环却可以达到$+0.35$。把边界一改，不仅空间分布不同，稳定性也变了。这个事实需要明确展示，不能只画两条漂亮的频谱轮廓就省掉它。

也不必一下把末端接得和内部一样强。本篇把末端连接从极弱逐步增加，并比较15、31和61格链。有限矩阵的谱随连接强度连续变化，但更长的链可能对更弱的闭合连接敏感。指数变换能说明原因：端到端连接经过变换后，某个方向会带上随链长快速增长的因子。1

模式偏在右边，不保证从中间输入后右边最响

持续驱动要解的是另一个方程。若以频率参数$\omega$输入空间源$\mathbf b$，稳态响应为：

这个逆矩阵包含了激发位置、频率、所有模式之间的叠加以及损耗。非厄米系统的左右本征向量通常不同；若用模式展开，某个右模式被激发多少，还取决于相应左模式与声源的重叠。知道所有右模式的包络，并没有完成驱动问题。

本篇保持31格链和中间同一单位源不变，只增加均匀损耗。这个改动给所有本征值加上相同的负虚部，不改变任何右本征模的形状。结果却是：

损耗较小时，方向性传递可以把响应推到边界；损耗增大后，声音在传递过程中被更快消耗，最强响应回到输入附近。于是，是否在实验图上看到边缘最亮，还取决于实际驱动与衰减，而不只取决于本征模。

长期稳定，也可能先放大再衰减

为什么第一组开放链的所有模式都在衰减，短时间内的总振幅却仍可增加？因为这些模式通常不是彼此正交的。初始状态中，不同模式可以相互抵消；演化后抵消关系发生变化，总振幅就可能暂时变大，之后才衰减。

在本篇选定的振幅范数下：

决定“有没有某种状态能立刻增长”的，是右边这个厄米矩阵是否有正方向；决定长期指数稳定性的，则是$H$本征值的虚部。二者可以不同。

本库实际计算了传播矩阵$e^{-iHt}$的最大奇异值。$\gamma=0.15$时，虽然所有本征增长率均为负，在取样时间$t=13$附近仍得到约2345的最大范数平方增益；这是对不同时间分别寻找最有利初态的包络，不代表任意初态或某个固定声源都能达到。

当$\gamma=0.7$时，上式的厄米矩阵已经负定，相应范数从一开始就单调衰减。这里计算的是明确的振幅范数，不是测得的电声效率；若真实装置出现放大，供能与饱和必须由其具体实现来承担。

时间调制怎样从一个输入频率生成多个分量

现在回到2026年的声腔实验。输入频率为$f_0$，把耦合以频率$\Omega$周期改变，就会产生$f_0+m\Omega$的分量，其中$m$为整数。3

论文把这些分量称为谐波。它们在这里是相对于调制产生的边带，不必是输入频率$f_0$的整数倍。例如输入4000 Hz、调制85 Hz，可以出现3915和4085 Hz。

这不需要用材料本身的非线性去倍频。即使振幅方程保持线性，系数随时间变化，也会把原来单一频率的振动混合到别的频率。能量的交换要包括控制耦合的调制装置，不能只盯着最初的声源。

把时间耦合写成傅里叶级数，并把各边带振幅分别记为$\mathbf q_m$，它们满足一组相互连接的线性方程：

只有$m=0$这一层直接接到外部单频源；其他层通过调制耦合获得振幅。各频率又位于静态参考谱的不同位置，因此它们可以有不同的空间传递倾向。34

对论文使用的单带周期链，每个$k$处的哈密顿量只是一个标量，时间演化可直接积分，因此其准频率谱与时间平均哈密顿量存在简单关系。换成有限开放链，或换成多带系统，就不能据此无条件地用时间平均代替全部动力学。本篇的实际声腔响应仍解完整谐波方程，并另用分段时间演化核验。

复频谱的绕行，提供了什么信息

在这些一维模型中，可以把输入的实频率看作复平面上的一个参照点，再问周期谱绕它转了几圈、朝哪个方向转。这就是谱绕数的一种直观读法；要求参照点没有落在谱线上。

在所讨论的单带链中，谱的绕行方向与开边界下的方向性响应有关。它是有条件的诊断，不是单独保证某个有限样品在任何损耗下都能把声场放大的万能规则。

一条椭圆主要提供一个方向；加入下一近邻耦合后，谱可以形成方向相反的两个瓣。于是，不同频率可以落进不同的瓣，对应不同的聚集方向。3

13个声腔的实验，怎样实现左右分流

第一种装置用相邻声腔之间的被动管道提供约75 Hz的双向耦合，再叠加朝左的方波主动耦合。声腔中心频率约4000 Hz，损耗参数约38 Hz；方波高值约120 Hz，占空比为一半，平均主动耦合约60 Hz。3

注意，75 Hz和120 Hz在这里是耦合强度的频率单位，不是额外声源的音高。麦克风、放大器、移相和扬声器共同实现从一个腔到另一个腔的反馈；参数来自独立的二腔或三腔传输标定，并非完全无拟合的几何预测。作者在评审回复里还说明了声源归一化和组件筛选方法。45

第二种装置重新连接反馈通路，加入跨过一个声腔的耦合：一个方向的静态参数约为$-28$ Hz，另一个方向的方波高值约56 Hz。负号表示相位关系，不表示负声能。两种装置都包含13个声腔，从第7腔输入。

本篇采用这些四舍五入后的参数，分别解谐波矩阵和方波每段的时间演化。后者先求出一个周期后的状态映射，找到周期稳态，再对时间信号作傅里叶分解。主要结果如下：

这些方向和结构特征与原文及附录的实验比较相符。本库没有重新拟合原始声压数据；公开资料没有提供可逐点重算全部实验拟合的原始表格，因此不把这种相符写成逐点数值复现实验。34

把谐波截断从$|m|\le5$扩大到25，三个主要分量相对于独立时间解的复振幅整体误差，从最高约0.26%降到约0.002%；时间网格加倍的变化约为千万分之一到百万分之一量级。所算六组参数的周期映射特征乘子模长均小于1，因此这些有限模型确实有稳定周期响应。

这个检验很重要：求逆矩阵总能在不少参数下给出一组形式上的受迫解，但如果系统本身在增长，那组解未必是实际长期会达到的响应。本篇先检查了这一点，再谈频率分流。

同样的平均耦合，为什么占空比还能改变结果

方波开得更短，但每次开得更强，完全可以保持相同的平均值。论文正是固定$R\kappa_0=\bar\kappa$，再改变占空比$R$，以保持时间平均周期谱不变。3

但平均值相同，并不意味着非零傅里叶分量相同。对以时间原点为中心的实方波：

正负一次分量相同，取决于这里明确的波形与相位约定。$R$变大，第一次边带耦合的幅度下降；边带之间相互转移的权重也会随之改变。只看平均哈密顿量，会漏掉这部分频率转换机制。

本库在3930 Hz输入、170 Hz调制下重算了$R=0.4$到0.8。同侧链的上边带权重约从11.3%降至1.27%；两侧链约从33.6%降至4.78%。这里按论文口径，只在$m=-1,0,+1$三个分量内归一化。3

再把$R$设为1，耦合始终打开、平均值仍然相同，时间调制消失。本库对照中，两侧边带回到数值舍入量级。参考的平均谱还在，频率转换却不再发生；这直接显示了两部分机制各做什么。

这个分母不能被称为“全系统能量效率”。在本篇六组代表性计算中，$|m|\le25$范围内还有约0.2%至4.6%的振幅平方落在三个主要分量之外；此外，实验并未因此测定所有输入电功率、输出声功率与耗散。声压平方的相对分配和器件效率，是不同层面的量。

占空比也不是越小越好。固定平均值时，同侧链的高值会从$R=0.8$时75 Hz升到$R=0.4$时150 Hz。作者说明，更低占空比需要更大的高值，装置会遇到饱和与失效，因此实验只展示0.4至0.8。把某项权重的单调趋势无限外推，会遗漏真实控制器的范围。5

所有皮肤效应都需要主动放大吗

不需要把“非厄米”与“主动放大器”画等号。2026年的另一项研究，用双层声学结构、交错连接与额外吸声损耗，报告了被动的$Z_2$皮肤效应。6

这项实验有21个元胞、84个声腔。额外损耗放在特定连接管中；研究者测量各腔的声压振幅和相位，再把上下层的信号组合成两个“伪自旋”通道。一个通道偏左，另一个偏右。去掉额外损耗的对照中，两个通道的空间响应恢复得更对称。

这里的伪自旋是两层振幅与相位的数学组合，不是声波突然有了电子的真实自旋。实验中的左右分流也需要按这些组合分别读取；只把所有声压强度混在一起，会掩盖所要分辨的方向性。

本文按其简化矩阵做了一个代数核对：两个通道的复传播因子模长分别约为0.775和1.291，一小于1、一大于1，对应相反的指数包络。它们各自的普通行列式谱绕数相反，相乘后的总绕数却为零。这说明前面单带链的一种绕数诊断，不能不加区分地拿来否定多通道皮肤效应。

这只是对原文块对角模型的复核，不是重跑84腔有限元或实验，也没有独立计算其完整Pfaffian拓扑不变量。实际被动装置还含共同衰减；简化相对矩阵里出现正负虚部，并不等于实物提供了增益。该实验的附录和原始数据本篇未取得，限制保留在来源记录。6

怎样判断一张“边缘更亮”的图靠不靠谱

关键不在于标签多新，而在于它有没有把不同问题分开。

边界怎样成为机制的一部分

这类系统提供了一条清楚的因果链：局部方向性耦合在反复传递中积累成空间偏置；开放或闭合的末端条件决定哪些全链模式能够成立；指定的声源、频率与损耗再决定实际响应；若耦合还随时间变化，就进一步决定不同频率之间怎样转移振幅。

这也是本篇较稳的结论。2026年的两类声学实验进一步展示：既可以用实时主动调制把频率分到不同边界，也可以用结构与损耗把不同伪自旋通道分开。它们没有因此证明一种普遍的增益优势，也不能只凭边缘聚集就判断器件效率。

本篇自己的复算覆盖有限矩阵、六组13腔模型的时间与频率两种解，以及被动平台简化矩阵的部分代数关系。原始实验数据、全波声场、电子反馈延迟、噪声和非线性饱和没有被重新验证。把这些边界留下，仍足以理解为什么一条链的最后一处连接、一个损耗参数或同平均值的另一种时间波形，可以改变看得见的整体响应。