前面两期，我们回溯了声表面波滤波器的理论起源和声波滤波器的演进路径；本次我们将解读声表面波领域经典综述《Coupling-of-Modes Analysis of SAW Devices》（Victor Plessky、Julius Koskela，2000）。Victor Plessky与Julius Koskela在文章中系统梳理了COM理论的基础原理、模型体系、改进扩展、参数提取与工程应用路径，解析了COM模型的核心思想、数学结构、参数物理意义与工程落地方法，是声表面波领域中一篇必读之作。本文对原文中公式详细推导进行了省略，仅梳理其逻辑，鉴于笔者水平有限，不可避免出现部分疏漏，欢迎指正。
01 引言：周期结构波传播研究百年回顾
如果拆开一部现代智能手机的主板，在射频前端电路中，你会看到几颗米粒大小的黑色芯片。它们就是支撑整个移动通信系统的核心元器件——声表面波（Surface Acoustic Wave, SAW）滤波器与体声波（Bulk Acoustic Wave, BAW）滤波器。在纷繁复杂的电磁信号里，它们如同精细的滤网，只放行指定频段的通信信号，把各类干扰尽数阻隔。正是这些小小的元器件，默默守护着手机通话、上网的稳定，避免信号串扰、网络卡顿。
过去二十年间，移动通信产业的爆发式增长，推动SAW滤波器进入了高频、低损、小型化的全新阶段。与传统横向耦合滤波器不同，现代SAW器件的电极厚度与声波波长已处于同一量级，电极的质量加载会显著降低声速，同时产生强烈的声波反射——正是这些反射，将声能束缚在器件内部，成为高性能SAW器件的工作基础。
设计这类高性能器件，离不开精准高效的仿真工具。目前行业内存在两类主流设计方案：一类是基于有限元法（FEM）、边界元法（BEM）的纯数值仿真，它直接从材料常数和器件几何推导器件特性，理论精度最高，但计算量极大，无法直接用于器件迭代优化；另一类是以COM模型为代表的维象模型，它用简洁的数学框架抓住SAW传播、反射、激发的核心物理规律，计算速度极快，可直接嵌入优化算法，是工业界SAW设计的主流标准。
正如Plessky和Koskela在论文开篇所言：“The coupling-of-modes formalism is a particular branch of the highly developed theory of wave propagation in periodic media, which has an exciting history of more than 100 years.”
Victor Plessky和Julius Koskela
1.1 理论溯源：从基础物理到COM框架的诞生
COM 理论并非声表面波领域的独有理论，它根植于百年以来学界对周期介质中波传播特性的基础物理研究，凝聚了数代学者的研究成果。
1. 理论源头：布拉格反射的发现 1913年，英国物理学家布拉格父子在研究X射线晶体衍射时，提出了著名的布拉格定律：当波入射到周期排列的结构上时，若周期p等于入射波半波长的整数倍（p=nλ/2），不同周期单元的反射波会相位一致、相干增强，形成极强的总反射。这个规律是所有波的共性：无论是X射线、电磁波、光波还是声波，只要遇到周期结构，就会出现布拉格反射。即使单个周期单元的反射强度极弱，只要周期单元数量足够多，总反射率也能逼近100%，最终形成波无法自由传播的“禁带（stopband）”，这是COM理论最核心的物理源头。值得一提的是布拉格父子二人也因此获得1915年诺贝尔物理学奖；
William Henry Bragg与William Lawrence Bragg
2.数学基石：从Mathieu方程到Floquet-Bloch理论
1868年，法国数学家Émile Mathieu提出了Mathieu方程，首次用数学形式描述了周期介质中的波动行为，成为周期系统波动分析的基础方程。
Duhem, P. (1892). "Émile Mathieu, his life and works". Bull. Amer. Math. Soc. 1 (7): 156–168. doi:10.1090/s0002-9904-1892-00067-5.
1883年，法国数学家Gaston Floquet提出了Floquet定理：线性周期微分方程的解，可以表示为一个周期函数与指数函数的乘积。1928年，瑞士物理学家Felix Bloch将这一定理推广到固体物理领域，提出了描述晶体中电子运动的Bloch定理。
Achille Marie Gaston Floquet与Felix Bloch
Floquet-Bloch理论指出：任何周期结构中的波，都可以分解为一系列离散的空间谐波（Floquet谐波），这些谐波之间存在相互耦合。在绝大多数情况下，只有少数几个谐波的振幅足够大，其余高次谐波的贡献可以忽略不计——这个结论，正是C这套理论的核心思想非常简洁：在布拉格反射占主导的频率附近，整个波场可以用两个最主要的谐波来近似描述——一个正向传播的波，一个反向传播的波，两者通过周期结构相互耦合，其余高次谐波的贡献可以完全忽略，这也是 “耦合模” 命名的由来。随后，麻省理工学院的Hermann A. Haus教授成为耦合模理论的集大成者，将这套理论系统化、严谨化，为其后续进入声学领域铺平了道路。OM理论的核心数学依据。
3.COM的诞生：电磁学与光学的理论积累
“耦合模”这个名称，最早出现在1950年代的光学与电磁学研究中。1960年代，行波管领域的先驱J.R. Pierce等人，将周期结构中波的耦合行为系统化，发展出了一套完整的耦合模理论框架。
John R. Pierce Holding An Early Traveling Wave Tube.
这套理论的核心思想非常简洁：在布拉格反射占主导的频率附近，整个波场可以用两个最主要的谐波来近似描述——一个正向传播的波，一个反向传播的波，两者通过周期结构相互耦合，其余高次谐波的贡献可以完全忽略，这也是 “耦合模” 命名的由来。随后，麻省理工学院的Hermann A. Haus教授成为耦合模理论的集大成者，将这套理论系统化、严谨化，为其后续进入声学领域铺平了道路。
Hermann A. Haus
1.2 周期结构中的波传播基础
要理解COM理论，首先要回到周期介质中波传播的通用物理模型。对于周期加载的波，我们可以用经典的波动方程（质量负载下的）描述：
其中k0=ω/ν0是自由衬底的波数，ζ(x)是x方向周期为p的加载密度，描述金属电极带来的电负载与机械负载。
根据Floquet定理，周期结构中的本征模满足
因此波场可以展开为无穷多Floquet谐波的叠加：
这些谐波并非独立存在，而是通过周期负载相互耦合。在布拉格频率附近，只有正向入射波和反向反射波两个谐波的振幅足够大，其余高次谐波可以完全忽略——这就是耦合模近似的核心逻辑。
1.3 耦合模理论的基础形式
耦合模理论有两种等价的表达形式，分别适用于不同的分析场景：
代数形式：通过忽略弱耦合高阶谐波，直接从Floquet谐波的耦合方程截断得到，保留了完整的物理精度，是严格数值分析的基础。
微分形式：假设波的振幅在单个周期内变化极慢，将代数方程线性化得到一阶微分方程组，是工程化COM模型的基础。
两种形式都能推导出COM理论最核心的成果：色散关系
其中q是波数，δ是失谐因子，κ是反射系数。这个公式清晰地描述了禁带的形成：当∣δ∣<∣κ∣时，q为纯虚数，波呈指数衰减，无法自由传播，对应布拉格禁带；禁带的相对宽度完全由反射系数决定：Δf/f0=∣κp∣/π。
耦合模色散关系
需要特别说明的是：SAW领域的“耦合模”，大多数情况下指的是周期栅格中反向传播的两个相同模式之间的耦合，而非不同模式之间的耦合，“分布式反射模型”其实是更准确的描述，但“耦合模”已成为行业通用术语。当然，SAW器件领域也存在不同类型表面声波模式间的耦合现象，比如在周期性电极阵列中，纵声表面波转换为瑞利波，或是瑞利波与兰姆波模式发生耦合，这类问题同样可以利用耦合模理论进行有效分析。
02 SAW器件的COM模型
耦合模（COM）理论与声表面波器件的结合，最早始于 20 世纪 70 年代中期。1976年，由Y. Suzuki及其合作者在IEEE超声研讨会上发表的论文《Some Studies on SAW Resonators and Multiple-Mode Filters》中率先完成。紧接着，Hermann A. Haus于1977年在Journal of Applied Physics上发表《Modes in SAW grating resonators》，从波动方程的本征模分析出发，系统推导了周期性栅格中SAW模式的特性，为COM方法的数学形式化奠定了坚实基础。二者殊途同归，正式拉开了 COM 理论应用于声表面波领域的序幕。历经数十年深耕完善，如今 COM 理论已发展为一套体系完备、技术成熟的声表面波器件建模方法。
Y. Suzuki, H. Shimizu, M. Takeuchi, K. Nakamura and A. Yamada, "Some Studies on SAW Resonators and Multiple-Mode Filters," 1976 Ultrasonics Symposium, Annapolis, MD, USA, 1976, pp. 297-302, doi: 10.1109/ULTSYM.1976.196685.
2.1 COM核心方程与P矩阵形式
SAW器件的COM模型，最终凝练成了三个一阶微分方程组，用最简洁的形式描述了声波传播、反射、激发与电流生成的全部物理过程：
同步均匀声表面波谐振器
物理含义：R(x)：正向传播声波振幅；S(x)：反向传播声波振幅；V：外加电压；I：电流；δ=2π(f−f0)/ν − iγ ：失谐因子。
1979年，学者G. Tobolka提出了具有里程碑意义的P矩阵（P-matrix） 概念，彻底解决了复杂SAW结构的级联问题。P矩阵是一个3×3的线性变换矩阵，它将SAW器件的三个端口——两个声学端口（左右方向传播的声波）和一个电学端口（电压与电流）的输入输出关系完整描述：
P矩阵示意图
其中：
P11、P22 是左右端口的声反射系数，P12=P21 是声透射系数；
P13、P23 是电压激发声波的效率，P31=−2P13、P32=−2P23 是声波产生电流的效率；
P33 是器件的总导纳。
G. Tobolka, "Mixed Matrix Representation of SAW Transducers," in IEEE Transactions on Sonics and Ultrasonics, vol. 26, no. 6, pp. 426-427, Nov. 1979, doi: 10.1109/T-SU.1979.31128.
对于叉指换能器、反射栅等各类 SAW 基本单元，P 矩阵都能实现精准描述，相当于为每一个功能模块建立了标准化的数学模型。依托成熟的矩阵级联规则，工程人员可像拼接模块一样组合不同单元的 P 矩阵，直接求解整体器件特性，省去了反复求解复杂波动方程的步骤，也为 COM 模型走向工程实用化筑牢了根基。
2.2 均匀周期结构的解析解
众多学者的持续探索，推动 COM 模型从理论研究全面落地工程设计。1982年，C. S. Hartmann、P. V. Wright、R. J. Kansy与E. M. Garber在Ultrasonics Symposium上发表《An analysis of SAW interdigital transducers with internal reflections and the application to the design of single-phase unidirectional transducers》，首次在COM框架内包含换能器内部多次反射效应，并给出了均匀结构的解析闭式解。
在均匀结构中，COM参数均为常数。对于这一极其重要的特例，COM方程组构成一个非齐次的一阶线性微分方程组。此时，通解可表示为对应齐次方程组解与特定解的线性组合——前者对应短路栅本征模，后者对应激发场。
短路栅与开路栅
1.短路栅本征模：对应短路栅（电极接地）中自由传播的正向波与反向波，满足色散关系
其禁带的上下边缘为：
2.开路栅本征模：对应电极悬浮的插指结构，此时电流密度为零，电压由声波振幅决定，等效为修正后的失谐因子与反射系数，形成独立的开路禁带。
对比开路栅与短路栅：开路栅的等效反射系数包含了一个额外的"电气反射"项；这解释了为什么：短路栅的禁带位置主要由机械反射决定，开路栅的禁带位置受机电耦合强度影响。
3.激发场解：对应外加电压驱动下的强制声场，描述换能器的电声转换过程。
基于这些解，均匀换能器、反射栅的所有P矩阵分量都可以通过代数公式直接计算，无需复杂的数值求解，极大降低了COM理论的使用门槛。
2.3 五个核心参数：决定器件声学特性的关键要素
整个COM模型，完全由五个独立的核心参数支配，它们决定了SAW器件的全部性能：
1.声速ν（单位：m/s）：声表面波在周期电极栅格中的传播速度。金属电极的质量加载与电加载通常会使声波速度低于自由衬底的声速。
2.反射系数κ（单位：1/m）：单位长度的反射强度，衡量周期栅格将正向波耦合为反向波的能力，直接决定了禁带的宽度。
3.换能系数α（单位：Ω−1/2）：衡量压电效应的效率——外加电压激发声波，以及声波反向产生电流的能力。
4.单位长度电容C（单位：F/m）：叉指换能器单位长度的静电电容，描述器件的电学特性。
5.衰减系数γ（单位：Neper/m）：声波传播过程中的能量损耗，包括热声子散射、缺陷散射、体波辐射等所有损耗机制。
1985年，D.P. Chen与Haus发表于IEEE Transactions on Sonics and Ultrasonics的论文《Analysis of metal-strip SAW gratings and transducers》中，证明了这五个参数并非完全独立：它们全部可以由开路栅格和短路栅格的色散关系唯一确定，彻底解决了“COM参数从哪里来”的核心问题。
为了消除器件尺寸的影响，行业通用以中心频率波长λ0为基准的归一化参数体系：归一化反射系数κp、归一化换能系数αp=αλ0、归一化电容Cp=Cλ0、归一化衰减γp=γλ0。这些参数是尺度不变的，不同频率、不同尺寸的SAW设计，可以直接通过归一化参数进行性能对比与优化。
Analysis of Metal-Strip SAW Gratings and Transducers
2.4 模型的固有局限性
COM模型的简洁性，来自于一系列合理的近似假设，这也带来了它的固有边界：
假设声波的空间结构不随位置、波数、频率变化，这一假设对水平剪切型声波严重不成立；
假设周期扰动足够弱，忽略高次Floquet谐波的贡献；
仅适用于均匀周期结构，周期突变处（如换能器与反射栅的间隙）的相位变化、体波散射等效应无法直接描述；
本质是一维模型，无法描述二维波导效应、总线声辐射、三维寄生效应。
2.5 与其他SAW模型的对比
目前SAW领域的主流维象学模型还包括混合矩阵模型、等效电路模型等。如果参数准确，所有模型的计算结果几乎一致，但COM模型具有不可替代的优势：
物理图像最清晰，直接对应声波的振幅、相位等物理量，而非等效电路中的虚拟量；
扩展性最强，寄生效应、二次效应可以很方便地纳入模型框架；
计算速度快，可直接嵌入优化算法，支持数千个频率点的快速仿真。
03 COM模型的改进与扩展
经典COM模型对瑞利型SAW的建模非常出色，但面对实际工程中的各种非理想效应、特殊波型，需要进行针对性的改进与扩展。这些改进分为两类：一类是不改变核心方程的“补丁式改进”，另一类是针对水平剪切波的根本性扩展。
3.1 补丁式改进：兼容实际工程非理想效应
3.1.1 电极电阻效应
电极电阻率会带来欧姆损耗、电压沿器件分布不均、额外色散等问题。1989年，P.V. Wright提出了广义COM形式，通过简单的变量替换V →V – dI(x)/dx，将电极电阻率完整纳入COM方程，无需改变模型的核心结构。
包括电极电阻IDT等效模型
3.1.2 速度色散与体波辐射修正
经典COM假设声速为常数，但实际中声速会随相对电极厚度h/λ变化，导致宽频范围内的频率偏移。V. P. Plessky、D. P. Chen 和 C. S. Hartmann在《‘Patch’ improvements to COM model for leaky waves》文章中通过引入频变声速ν(f) ，可以很好地修正这一”琴状偏差”。
C. S. Hartmann
对于高频下的体波辐射损耗，可以通过在衰减系数中加入频变修正项，以及在导纳中加入体波激发的经验项，实现宽频范围内的精准建模。
3.2 水平剪切波：针对特殊波型的关键技术突破
经典COM理论的最大挑战，来自水平剪切（SH）型声波——包括表面横波（STW）和漏声表面波（LSAW），其中旋转Y切割LiTaO₃和LiNbO₃衬底上的LSAW，是目前手机射频SAW滤波器的绝对主力波型。
3.2.1 SH波的核心难题
SH波有两个经典COM无法描述的核心特性：
速度强色散：SH波的传播速度会随频率发生显著变化；
能量分布的频率相关性：SH波的能量在表面的深度会随频率和栅格加载强度显著变化。当频率升高到一定阈值，SH波会从束缚在表面的声波，转变为向衬底内部辐射的体波，能量快速衰减——经典COM假设波的空间结构固定，完全无法描述这个过渡过程。
3.2.2 Plessky的二参量模型
1993年，Plessky在论文中《A two parameter coupling-of-modes model for shear horizontal type SAW propagation in periodic gratings》提出了专门针对SH型SAW的二参量COM模型，首次解决了这个难题。Plessky将无限周期栅格中的色散关系近似表示为：
式中引入ε和η，其中，ε表示周期栅格中传播时与反方向波的耦合；η表示漏波与快切变体波的相互作用。该方法与经典耦合模（COM）模型的主要区别在于：它能够精准描述波的空间分布形态，且本征模可由表面波与体波叠加构成。模型可自然体现高频下体波散射的产生过程，以及随之而来的波衰减效应。
3.2.3 STW-COM模型：框架的融合
1995年，Abbott和Hashimoto在论文《A coupling-of-modes formalism for surface transverse wave devices》将Plessky的色散模型与标准COM框架结合，提出了STW-COM模型。他们引入了5个新参量模型（ΔV，ΔB，κ，κB，ηB），在保留COM原有P矩阵级联、解析解等工程优势的同时，正确描述了STW和LSAW的色散、衰减与体波辐射特性。下式是STW-COM色散公式：
Benjamin P. Abbott与Prof. Ken-ya Hashimoto
3.2.4 未竟的挑战
即使是STW-COM模型，也并非完美：模型在远离禁带的频率区域，渐近行为与严格数值结果存在偏差；SH波的换能系数α会随频率显著变化，与经典COM的常数假设不符；换能器端面、间隙等不连续处产生的额外体波辐射，也尚未完全纳入模型。直到今天，SH波的高精度COM建模，仍是SAW领域的活跃研究方向。
04 COM参数的提取
COM模型是维象学模型，它的精度完全取决于核心参数的准确性。经过几十年的发展，学界形成了实验测量与数值仿真两大类成熟的参数提取方法。
4.1 实验提取：从测试结构到参数映射
最可靠的参数提取方法，来自对专门设计的测试结构的电学测量，其中最经典的是同步单端口谐振器法。
4.1.1 同步单端口谐振器
1993年，Plessky和Hartmann首次提出：同步单端口谐振器（由中间的叉指换能器和两侧的反射栅组成，所有单元周期相同）是最理想的COM参数提取测试结构。
同步单端口谐振器
通过该结构的导纳曲线的关键特征，与COM参数有一一对应的关系实现COM参数的提取。如果COM理论足够精确，那么通过严格计算得到的电容以及短路栅和开路栅的禁带边缘，可以轻松提取出速度、耦合系数和换能系数。在对称结构中，开路栅和短路栅的禁带会有一个重合的边缘；而在单向结构中，禁带的边缘会彼此分离。
(a)/(b)/(c)双向结构的短路栅（sc）和开路栅（oc）本征模的色散曲线，（d）单向结构的短路栅（sc）和开路栅（oc）本征模的色散曲线
远离禁带的低频/高频区域的导纳水平，直接确定单位长度电容C；谐振峰的品质因子（Q值），反映衰减系数γ的大小。在所有参数中，衰减γ的提取是最棘手的。它受多种因素共同影响，还很难与电极电阻带来的损耗区分，因此衰减参数的提取通常需要多特征的综合拟合权衡。
4.1.2 谐波导纳自动化提取
1999年，Koskela、Plessky和Salomaa提出了基于谐波导纳的自动化参数提取方法。该方法的核心是：对无限长周期栅格施加不同相位的周期性电压激励，计算其谐波导纳函数，其中开路栅格的色散对应谐波导纳的零点，短路栅格的色散对应谐波导纳的极点。通过拟合这些零点和极点，就可以逆向提取出所有COM参数，整个过程可以通过算法自动完成，无需人工干预。
4.2 数值仿真提取：灵活的参数扫描
与实验提取相比，数值仿真提取的优势在于灵活性：研究者可以独立改变电极厚度、占空比、材料常数等参数，系统研究每个设计变量对COM参数的影响，这是实验方法难以实现的。
目前主流的数值提取方法，是基于边界元法（BEM）或有限元法（FEM）仿真周期栅格的谐波导纳，再与COM理论的谐波导纳表达式进行拟合。通过调整COM参数，让理论导纳与仿真导纳最优拟合，即可得到完整的参数集。
数值提取也存在局限性：电极形状的精确建模、材料参数的不确定性、短器件与无限周期结构的场分布差异，都会影响提取精度。因此，数值提取得到的参数，通常还需要实验测试数据进行最终校准。
05 COM理论的器件设计与工程落地
COM理论的最大价值，在于它将复杂的SAW器件设计，转化为标准化的模块化工程流程，支撑了现代射频SAW滤波器的大规模量产。
5.1 SAW设计的基础构建模块
复杂的SAW器件，都可以拆解为几个基础功能模块，每个模块用对应的P矩阵描述，通过级联实现完整器件的仿真：
叉指换能器（IDT）：实现电声转换的核心单元，强反射电极下的IDT响应接近LC谐振器，是谐振型滤波器的基础。
反射栅：实现声波反射的布拉格光栅，长反射栅的禁带内反射率可接近100%；通过加权反射系数，可以有效抑制反射系数的旁瓣，优化滤波器的带外特性。
间隙：换能器与反射栅之间的自由衬底区域，用于调整声波相位，其带来的体波散射损耗需要作为维象参数纳入模型。
5.2 单端口SAW谐振器
单端口谐振器是所有阻抗型滤波器的基础单元，COM模型可以精准仿真其全部特性。
1.同步均匀谐振器：最经典的谐振器结构，也是COM参数提取的标准测试结构。其谐振发生在禁带的低频边缘，具有最大的压电耦合系数，是梯型滤波器的核心单元。
2.同步“hiccup”谐振器：在换能器中心引入1/4波长间隙，将谐振频率移到禁带中心，大幅降低器件对工艺偏差的敏感度，主要用于振荡器应用。
同步“hiccup”谐振器
3.SAW标签：由换能器和多个反射栅组成，通过COM模型可以精准仿真其频域与时域响应，实现无线传感与识别功能。
5.3 双端口谐振器与耦合谐振滤波器
5.3.1 耦合谐振滤波器（CRF）
耦合谐振滤波器是射频SAW滤波器的主流结构之一，典型的三换能器CRF由中心输入换能器、两侧输出换能器与两端反射栅组成。COM模型可以分别构建每个单元的P矩阵，通过级联得到完整的Y参数矩阵，仿真滤波器的插入损耗、带外抑制、矩形度等关键指标，优化换能器指条比、间隙尺寸等核心设计参数。
耦合谐振滤波器（CRF）
5.3.2 横向耦合滤波器
由两个平行的窄孔径谐振器通过声场重叠实现耦合，具有极低的通带损耗和极高的带外抑制，主要用于窄带中频滤波器。COM模型可以很好地描述其波导模式的耦合特性，实现精准的设计优化。
横向耦合滤波器
5.4 阻抗元滤波器（IEF）
梯型阻抗元滤波器是目前手机射频前端应用最广泛的SAW滤波器结构，它由串联臂和并联臂的单端口谐振器交替连接组成。
Ladder Filter
COM模型可以分别优化串联臂和并联臂谐振器的谐振、反谐振频率：通带内，串联臂阻抗极低，信号无损耗通过；阻带内，并联臂将信号短路到地。COM模型还可以分析谐振器的串并联组合，实现50Ω系统阻抗匹配，优化滤波器的插入损耗与功率耐受能力。
除此之外，基于阻抗元还可以衍生出平衡桥滤波器（BBF）、陷波滤波器等结构，满足不同的应用需求。与其他滤波器结构相比，阻抗元滤波器具有更高的功率承受能力、更宽的通带和更低的损耗，是5G时代射频滤波器的主流技术路线。
5.5 工艺敏感性分析：量产良率的保障
COM理论的另一个关键工程价值，是评估制造工艺偏差对器件性能的影响：
· 电极厚度：主要改变波速与反射率，对换能系数、电容影响微弱；横波机电耦合强度会随厚度增大而提升。厚度增加会使多数基底材料的波速下降、反射率上升，进而拉低谐振频率。
· 电极宽度 / 金属化率：会同时影响波速、反射率、换能系数与电容，改变谐振、反谐振频率。
· 电极形貌：影响复杂，其中反射率对形貌变化最敏感；波速由电极总质量决定，换能系数和电容主要取决于电极表面积。
通过COM模型分析工艺参数的敏感性，工程师可以提前优化设计参数，降低器件对工艺偏差的敏感度，大幅提升大规模生产的良率。
06 边界与未来
在论文的最后，Plessky和Koskela坦诚讨论了COM模型仍未解决的核心挑战，也对其未来发展做出了精准的预言。
6.1 现存的核心挑战
1.水平剪切波的精确建模：SH波在不连续处的体波辐射、换能参数的频变特性，仍未完美纳入COM框架；
2.波导效应的描述：强压电衬底上SAW的强各向异性，导致二维波导行为十分复杂，目前的一维COM模型无法完全精准描述；
3.三维寄生效应：封装寄生参数、电极端面散射、总线声辐射等三维效应，仍需结合等效电路或数值方法补充建模。
6.2 未来的发展方向
Plessky和Koskela在论文中预言：“作为对不同物理本质的波在周期介质中传播的通用简化描述，COM模型将会被一直使用下去。但对于精确设计，趋势将是使用数值FEM/BEM方法——至少是用于COM参数的确定。”
他们同时指出，未来滤波器技术的核心是“融合”：滤波器与射频系统的融合、不同声学技术路线的融合、理论模型与材料工艺的深度融合。进入5G/6G时代，TF-SAW、IHP SAW、XBAR等新型声学器件不断涌现，但无论器件形态如何演变，COM理论作为周期结构声波传播的核心基础，仍将持续发挥不可替代的作用。
07 结束语
COM理论的百年发展史，是一段从纯粹物理探索走向成熟工程落地的经典历程。该理论最早起源于19世纪物理学家对周期介质中波传播规律的基础性研究，历经电磁学、光学领域的层层迭代与理论完善，在二十世纪七十年代正式落地应用于声表面波器件研究领域。伴随着移动通信技术的飞速迭代革新，COM理论逐步成为核心器件研发的关键理论根基，支撑起数十亿人日常移动通信的正常运转。
自该文献发表以来，COM模型始终处在持续优化与拓展的发展进程中。研究人员在传统一维基础模型之上，逐步衍生出二维COM模型、多模耦合COM模型等优化版本，极大弥补了经典模型的适用局限，适配了更为复杂的器件结构与工作场景。与此同时，依托计算机算力与数值仿真技术的飞速发展，FEM、HCT-FEM等数值计算方法逐步成熟，在声表面波器件仿真、结构优化等开发环节的应用频次持续提升。即便各类新型仿真算法与计算模型不断涌现，COM理论依旧是声表面波器件设计、机理分析与性能预判的核心基础理论与主流设计工具。其兼具物理机理清晰、计算效率高、工程适配性强的核心优势，能够精准刻画声表面波的激发、传播与耦合损耗特性，至今仍是行业内器件正向设计、参数优化与产品迭代不可或缺的核心理论支撑，在声表面波器件研发体系中占据着不可替代的核心地位。
参考文献：
[1] V. Plessky and J. Koskela, "Coupling-of-Modes Analysis of SAW Devices," Int. J. High Speed Electron. Syst., vol. 10, no. 4, pp. 867-947, 2000.
[2] W. L. Bragg, "The Diffraction of Short Electromagnetic Waves by a Crystal," Proc. Camb. Philos. Soc., vol. 17, pp. 43-57, 1913.
[3] É. Mathieu, "Mémoire sur le mouvement vibratoire d'une membrane de forme elliptique," J. Math. Pures Appl., vol. 13, pp. 137-203, 1868.
[4] P. Duhem, "Émile Mathieu, his life and works," Bull. Amer. Math. Soc., vol. 1, no. 7, pp. 156-168, 1892. doi:10.1090/s0002-9904-1892-00067-5.
[5] G. Floquet, "Sur les équations différentielles linéaires à coefficients périodiques," Ann. Sci. Éc. Norm. Supér., vol. 12, pp. 47-88, 1883.
[6] F. Bloch, "Über die Quantenmechanik der Elektronen in Kristallgittern," Z. Phys., vol. 52, no. 7-8, pp. 555-600, 1928.
[7] J. R. Pierce, Traveling Wave Tubes. New York: D. Van Nostrand Company, 1950.
[8] H. A. Haus, "Coupled-mode theory," Proc. IEEE, vol. 79, no. 10, pp. 1505-1518, 1991.
[9] Y. Suzuki, H. Shimizu, M. Takeuchi, et al., "Some studies on SAW resonators and multiple-mode filters," in Proc. 1976 IEEE Ultrason. Symp., 1976, pp. 297-302.
[10] H. A. Haus, "Modes in SAW grating resonators," J. Appl. Phys., vol. 48, no. 12, pp. 4955-4961, 1977.
[11] H. A. Haus, "Bulk scattering loss of SAW grating cascades," IEEE Trans. Sonics Ultrason., vol. SU-24, no. 4, pp. 259-267, 1977.
[12] G. Tobolka, "Mixed matrix representation of SAW transducers," IEEE Trans. Sonics Ultrason., vol. SU-26, no. 6, pp. 426-428, 1979.
[13] C. S. Hartmann, P. V. Wright, R. J. Kansy, and E. M. Garber, "An analysis of SAW interdigital transducers with internal reflections and the application to the design of single-phase unidirectional transducers," in Proc. 1982 IEEE Ultrason. Symp., 1982, pp. 40-45.
[14] D. P. Chen and H. A. Haus, "Analysis of metal-strip SAW gratings and transducers," IEEE Trans. Sonics Ultrason., vol. SU-32, no. 3, pp. 395-408, 1985.
[15] P. V. Wright, "A new generalized modeling of SAW transducers and gratings," in Proc. 43rd Annu. Symp. Freq. Control, 1989, pp. 596-605.
[16] B. P. Abbott, "A derivation of the coupling-of-modes parameters based on the scattering analysis of SAW transducers and gratings," in Proc. 1991 IEEE Ultrason. Symp., 1991, pp. 5-10.
[17] V. P. Plessky, "A two parameter coupling-of-modes model for shear horizontal type SAW propagation in periodic gratings," in Proc. 1993 IEEE Ultrason. Symp., 1993, pp. 195-200.
[18] B. P. Abbott and K. Hashimoto, "A coupling-of-modes formalism for surface transverse wave devices," in Proc. 1995 IEEE Ultrason. Symp., 1995, pp. 239-245.
[19] J. Koskela, V. P. Plessky, and M. M. Salomaa, "SAW/LSAW COM parameter extraction from computer experiments with harmonic admittance of a periodic array of electrodes," IEEE Trans. Ultrason., Ferroelect., Freq. Control, vol. 46, no. 4, pp. 806-817, 1999.
[20] N. Kimura and M. Takeuchi, "Direct numerical analysis of SAW coupling-of-modes equations and its applications to weighted single-phase unidirectional transducers and weighted reflectors," in Proc. 1999 IEEE Ultrason. Symp., 1999, pp. 203-208.
[21] V. P. Plessky and C. S. Hartmann, "Characteristics of leaky SAWs on 36-LiTaO3 in periodic structures of heavy electrodes," in Proc. 1993 IEEE Ultrason. Symp., 1993, pp. 1239-1242.
[22] S. Datta and B. J. Hunsinger, "First-order reflection coefficient of surface acoustic waves from thin-strip overlays," J. Appl. Phys., vol. 50, no. 9, pp. 5661-5665, 1979.
[23] K. Hashimoto and M. Yamaguchi, "Precise simulation of surface transverse wave devices by discrete Green function theory," in Proc. 1994 IEEE Ultrason. Symp., 1994, pp. 253-258.
[24] P. Ventura, J. Desbois, and L. Boyer, "A mixed FEM/analytical model of the electrode mechanical perturbation for SAW and PSAW propagation," in Proc. 1993 IEEE Ultrason. Symp., 1993, pp. 205-208.
[25] K. Hashimoto, G. Q. Zheng, and M. Yamaguchi, "Fast analysis of SAW propagation under multi-electrode-type gratings with finite thickness," in Proc. 1997 IEEE Ultrason. Symp., 1997, pp. 279-284.
[26] G. Kovacs, M. Anhorn, H. E. Engan, G. Visintini, and C. C. W. Ruppel, "Improved material constants for LiNbO₃ and LiTaO₃," in Proc. 1990 IEEE Ultrason. Symp., 1990, pp. 435-438.
[27] C. S. Hartmann, D. P. Chen, and J. Heighway, "Modeling of SAW transversely coupled resonator filters using coupling-of-modes modeling technique," in Proc. 1992 IEEE Ultrason. Symp., 1992, pp. 39-43.
[28] S. V. Biryukov, G. Martin, V. G. Polevoi, and M. Weihnacht, "Consistent generalization of COM equations to three-dimensional structures and the theory of the SAW transversely coupled waveguide resonator filter," IEEE Trans. Ultrason., Ferroelect., Freq. Control, vol. 42, no. 4, pp. 612-618, 1995.
[29] O. Ikata, T. Miyashita, T. Matsuda, T. Nishihara, and Y. Satoh, "Development of low-loss band-pass filters using SAW resonators for portable telephones," in Proc. 1992 IEEE Ultrason. Symp., 1992, pp. 111-115.
[30]Ken-Ya Hashimoto. (2000). Surface acoustic wave devices in telecommunications: modelling and simulation.Springer Berlin Heidelberg.
