📄 连分式与无源梯形网络：障板活塞与号筒阻抗的系统构造
评分：6.2/10 | 创新 1.2/2 | 技术严谨 1.2/1.5 | 实验充分 1/1.5 | 清晰度 0.8/1 | 影响力 0.7/1.5 | 开源 0/1.5 | 可复现 0.3/0.5 | 工程/实践 1/1.5
👥 作者与机构
- Ronald M. Aarts：机构信息未能在本次解析中与作者逐人可靠绑定
📌 核心摘要
输入是障板圆形活塞的辐射阻抗与带活塞mouth负载的号筒喉阻抗，输出是由正值电阻、质量与顺性构成的被动梯形网络，难点是原阻抗为超越函数且带分支与振荡波纹，不能用截断幂级数全局逼近。方法先对零点泰勒级数做维斯科瓦托夫算法得到C连分式，再做偶收缩得到J连分式并读出梯形元件，系数为正则直接对应被动可实现网络，截断后嵌套扩展且不改动已有元件。号筒部分把每段均匀波导传输矩阵看作莫比乌斯映射，从mouth向throat回代即为连分式，尾端直接接入活塞网络而保持整体被动，mouth网络误差一一传递到喉阻抗。以最大绝对误差|Z-Zexact|/ρ0cS为指标、在ka≤κ评测设置下，14元件嵌套梯形在κ=2时为1.8e-09、在κ=5时为4.3e-04，后者高于前者。以宽带误差为指标、在0.1≤ka≤50范围内评测时，四元件基网络为0.181、匹配高频1/s项的两点近似为0.153，后者低于前者。在同一四元件拓扑下以电容参数C为比较量，贝拉尼克拟合值1.89与两点近似封闭值1.85相差约2%、后者低于前者。结论适用于需电路仿真与时域仿真的集总换能器模型，超截止波纹与无限指数号筒通带截断不收敛仍需端接处理，原文未披露训练、推理或部署成本。
🔗 开源与复现资源
本次未形成可展示的已核验资源记录，开放状态尚未核实。
可达状态仅表示本次链接检查结果，不代表许可证、本文权重或运行复现已验证。
🧭 深度解读
为什么要用集总网络表示声阻抗
刚进入音频领域的研究生首先要分清两种建模路线。有限元与边界元把几何剖分成大量自由度，精度高但物理图像分散；集总元件网络只用质量、顺性与电阻描述输入端看进去的阻抗，低频时靠近端口的储能元件主导，高频时只剩电阻，两个极限可以直接从电路读出。论文研究的对象正是后者：希望对障板上圆形活塞的辐射阻抗和号筒喉阻抗，给出由正值恒定元件构成、无需优化即可任意加阶、且在电路仿真器与时域中直接可用的网络。
输入是阻抗的幂级数与高频渐近展开，目标是得到无源、可嵌套扩展的梯形电路，输出是元件值表与误差随频率和阶数的变化。必须保留的信息包括归一化方式、误差定义、适用 ka 范围与无源性成立的阶数上限。论文没有训练神经网络，也没有下载权重，全部数字来自随附 Python 包的高精度计算。ka 是波数与活塞半径的乘积，是频率的无量纲度量，误差越小越好。
已有网络综合路线与本篇的位置
经典网络综合已经解决有理正实阻抗的实现问题。Foster 给出无损耗两类元件阻抗的典范实现，Cauer 给出梯形连分式形式，Brune 处理含理想变压器的一般情形，Bott-Duffin 去掉变压器但元件数远超麦克米伦阶。论文明确本篇属于另一类问题：不是实现给定的有理函数，而是对超越的活塞辐射阻抗构造正实有理逼近。相关实践包括电路模拟中的渐近波形估计、经由 Lanczos 的帕德方法、保持无源的 PRIMA，以及用测量数据拟合极零点的向量拟合。
教学例子 arctan 与对数函数用来说明单点展开只在原点附近好、两点匹配才能兼顾无穷远，这些例子不能直接当成活塞阻抗的精度承诺。Beranek 四元件网络是经典选择，低频和高频极限正确但在 ka 约 2.5 附近的实部过冲与后续波纹拟合不好，论文把它作为四元件家族的一个成员重新解释，而不是在同条件下贬低其实用价值。
具体要解决的两个声学对象是什么
第一个对象是无限大障板上圆形活塞的辐射阻抗，按 ρ0cS 归一化后实部含贝塞尔函数 J1、虚部含斯特鲁夫函数 H1，时间约定为 ejωt。它几乎出现在每个障板换能器的集总模型中，电路与状态空间仿真需要把它变成有限网络。第二个对象是号筒喉阻抗，号筒被切成多段均匀波导，每段用传输矩阵级联，末端接活塞口阻抗。论文要回答的是：能否只用正电阻、质量与顺性搭出梯形网络，使每 1 阶都无源、加阶不改旧元件、且误差随阶数系统下降；以及口网络的误差如何传递到喉端。
评价时绝对误差指|Z-Zexact|/ρ0cS，只在汇总时同时给出相对误差以照顾极低频段，指标方向都是越小越好。号筒算例另需区分模型误差与近似误差，前者是平面波无损耗假设与真实号筒的差距，不在评估之列，后者才是离散与网络带来的偏离。
从幂级数到嵌套梯形的总流程
论文为活塞设计的路线可以概括为五步。先把阻抗在 s=0 处展开成处处收敛的幂级数，首项给出辐射质量 8/3π，次项给出低频辐射电阻系数；再用 Viskovatov 算法把级数变成 C 连分式；发现其系数从第三项起符号交替，不能直接对应正元件；于是做偶收缩，把相邻两层合并成 J 连分式，得到全为正的 dk 与 ek。
最后按导纳与阻抗交替解读画出梯形电路，截断即得有限网络，可选地重调末端一两个元件改善宽带误差。下面这张流程图是官方原图，先看箭头主路径再看可选分支，图中方框的文字与箭头走向都是像素可见的，不要猜测框外的内容。
看图路径： 1. 从左上阻抗泰勒级数出发，沿箭头数出四个方框的先后顺序；2. 注意偶收缩后系数为正才进入梯形网络分支；3. 区分主路的截断嵌套与下方的可选终端重调
论文图 1。原论文 Figure 1:：“The route followed for the piston: from the power series of the impedance to a nested ladder of positive elements.”。
图中左上是 s=0 处泰勒级数，中间是符号交替的 C 连分式，右上是系数为正的偶收缩，右下是 M、G、C、R 梯形，左下截断后得到嵌套无源网络，最左是可选的末端重调。正定性对 k≤100 经区间算术认证，这是后文无源性断言的边界，超出该范围只是猜想而非已证事实。整个构造是一次性高精度制表，仿真时只是普通电路求解。
连分式如何从函数算出来
连分式把函数写成嵌套分式，每多一层就多一个有理逼近。计算上有 3 条常用路径。Viskovatov 法从幂级数出发，每步做一次级数除法产生一个部分分母；qd 算法从系数比出发用菱形规则直接产生 Stieltjes 分式的系数；若只有采样点则用 Thiele 反差商插值。
论文用对数函数说明：幂级数在|x|<1 外发散，而连分式在除去分支割线的全平面收敛，x=3 处多项式部分和已发散到很大，连分式十余项即接近 ln4。等价性指分子分母同乘非零常数不改变值，这解释了不同书系数写法不同，也允许逐节缩放网络而不改变阻抗。
\[f=b_{0}+\cfrac{a_{1}}{b_{1}+\cfrac{a_{2}}{b_{2}+\cfrac{a_{3}}{b_{3}+\cdots}}}\]上式是连分式的一般形状，ak、bk 为部分分子分母，截断到 n 层即得 n 阶逼近，其分子分母满足 2 阶递推。实际计算中自顶向下可用修正 Lentz 算法避免溢出，自底向上则需预先定深。论文强调 Viskovatov 与 qd 在双精度下都会丢精度，活塞例子中把阶乘舍入到双精度就会出现虚假的负系数，因此高阶构造必须用精确或上 1,000 位浮点运算，但这只是一次性制表代价。
连分式 × 梯形网络： 连分式是把函数写成逐层分式嵌套的形式，每截断一层就得到一个有理逼近；梯形网络是串联阻抗与并联导纳交替级联的电路，其输入阻抗按从输入向后看的顺序自然展开成同样形状的分式。两者分工是：连分式提供系数的计算与截断规则，梯形网络提供每个系数对应的物理元件，组合后数学上的加阶就等于电路末端追加一节，系数为正则元件为正并保证无源。
下面这张梯形示意图是官方原图，先确认左右贯通的串联支路，再看每个节点下接地的并联支路，输入端向右看进去的顺序就是分式从外到内的顺序。
看图路径： 1. 确认串联 Z 与并联 Y 交替出现的位置；2. 从输入端向右看，体会分式中每一层倒数对应的并联关系
论文图 3。原论文 Figure 3:：“A ladder network. Its input impedance is the continued fraction (11).”。
该图串联 Z1、Z3 与并联 Y2、Y4 交替出现，输入阻抗正是连分式形状。读图时先确认左右贯通的串联支路，再看每个节点下接地的并联支路，输入端向右看进去的顺序就是分式从外到内的顺序，倒数对应并联关系。
活塞阻抗的级数与高频行为
活塞阻抗经归一化后时间尺度取 a/c，s=jka。拉普拉斯域形式把 J1 与 H1 换成修正贝塞尔 I1 与修正斯特鲁夫 L1，显示它是 s 的整函数，没有有限极点，只在无穷远有本性奇点。s=0 处级数系数分奇偶给出，首两项对应辐射质量与低频电阻。s=无穷处有 1/s 整数幂的发散渐近级数，首项趋于 1，对应高频极限为 ρ0cS，这正是两点帕德需要的形式。但频率轴上还有 J1 型振荡项，按(ka)-3/2 衰减，形成实部虚部的波纹，任何极点远离频率轴的有理函数都无法复现无限多个起伏，这构成所有有理方法的误差下限。延迟表示把阻抗写成对活塞弦长的连续延迟积分，有限支撑是直径传播时间，但有理网络只能近似而不能精确复现。
\[\frac{Z(\omega)}{\rho_{0}cS}=1-\frac{2J_{1}(2ka)}{2ka}+\mathrm{j}\,\frac{2\mathbf{H}_{1}(2ka)}{2ka},\]上式是归一化活塞阻抗的定义，S=πa2，J1 为贝塞尔函数，黑体 H1 为斯特鲁夫函数。本篇不涉及分数模型与价值评估模型，只有声阻抗、连分式与梯形网络，后文保持同一称呼，不套用判别器等未涉及角色。
单点与两点帕德如何分配匹配条件
设分子分母同为 n 次，分母常数项归一后有 2n+1 个自由系数，在 0 处匹配 k 个条件、在无穷远匹配剩余条件，两组都是线性方程。逐阶交替提高匹配阶数即得 M 连分式，系数由两点 qd 算法从两端级数算出。arctan 例子中 n=2、k=3 的闭式在 0 处斜率为 1、在无穷远极限为 π/2，两点误差在两端为零、中段隆起，随 n 增大在半轴上均匀下降；而单点在远端趋于 0 或发散，大 x 处相对误差很大。
实践要点是无穷远展开须化成 1/x 整数幂，含对数或根号时先变量代换或分解出渐近因子，k 取中间附近均匀误差最好，线性系统须高精度求解。初学者易误以为阶数越高处处越好，实际上两点误差峰值位于中段波纹区，活塞的高频波纹同样把峰值钉在中高频。
单点帕德逼近 × 两点帕德逼近： 单点帕德逼近只在 s=0 处匹配泰勒系数，适合极低频但不管无穷远极限；两点帕德逼近把自由系数分配到 0 和无穷远两端，同时匹配低频质量与阻尼行为和高频极限与 1/s 修正。两者搭配的理由是活塞阻抗低频趋于 0 而高频趋于特征阻抗，组合后低阶网络也能在两端都正确，代价是中间频段仍受波纹限制，均匀误差下降变慢。
梯形即连分式与活塞梯形的画法
梯形输入阻抗按串并交替展开即得连分式，Zk=sLk、Yk=sCk 时每个部分分母就是一个物理元件。论文用有理函数做 Cauer 综合演示：从高次幂长除得到 Cauer I 低通梯形，从低次幂长除得到 Cauer II 高通梯形，两种实现精确相等，拓扑直接反映展开点选在无穷还是零点。能否实现取决于商是否全为正，这与用偶奇部比值判稳的 Routh-Hurwitz 判据同构。RC 网络的驱动点阻抗是 Stieltjes 函数，极零交替在负实轴，S 连分式系数为正，截断仍是正 RC 梯形。RC 线的结论不能直接搬到活塞，因为活塞极点为复数、实部有过冲，不是 Stieltjes 类，这正是需要偶收缩的原因。
\[Z_{\mathrm{in}}=Z_{1}+\cfrac{1}{Y_{2}+\cfrac{1}{Z_{3}+\cfrac{1}{Y_{4}+\cdots}}}.\]上式是梯形输入阻抗的连分式写法，Z1 为首串联阻抗，Y2 为首并联导纳，依此类推。活塞梯形按导纳与阻抗交替解读，第 0 节为并联质量与电导，之后奇数节串联 C-R、偶数节并联 M-R，每元件多匹配一个泰勒系数。
Cauer 综合 × Stieltjes 函数： Cauer 综合指对有理阻抗在 s=无穷或 s=0 做辗转相除，每一步商对应一个电感电容或质量顺性元件，从而画出梯形电路；Stieltjes 函数指极零点交替在负实轴上的函数类，其 S 连分式系数全为正。分工是 Cauer 给出从展开位置到电路拓扑的对应，Stieltjes 给出系数为正从而截断仍可实现的充分条件，组合后解释了 RC 线为何任意截断都无源，而活塞需要偶收缩后才出现正系数。
下面这张活塞梯形是官方原图，先找到最左侧的并联质量，再向右数串联与并联的交替，最后想象在最右端加一节而不动左侧。
看图路径： 1. 找到最左侧并联的质量与电阻作为第 0 节；2. 向右区分奇数节串联 C-R 与偶数节并联 M-R 的交替规律；3. 想象在最右开路端追加一节而不动左侧
论文图 6。原论文 Figure 6:：“The piston ladder in the mechanical-impedance analogy (inductor=mass, capacitor=compliance, resistor=mechanical resistance).”。
该图左侧 Z 引出端下接并联质量与电导，第 0 节之后交替出现串联 C-R 与并联 M-R。读图时注意每加 1 阶只在最右开路端加一节，左侧值不动，这就是嵌套的电路含义。元件换算到物理单位需乘相应声学量，归一化值本身无量纲。
没有神经网络训练时实际计算了什么
本研究没有训练任何模型，因此不存在冻结参数、梯度路径、监督来源与重置时刻。实际计算分为 3 类。第一类是构造计算：用 1000 位浮点执行 Viskovatov 递推得到 C 连分式系数，再按收缩公式换算成 M、C、R、G 值，每元件对应多匹配一个泰勒系数。第二类是认证计算：在 4000 位区间算术下用整数阶乘与 π 区间包络输入，外向舍入保证每个除法区间不含零，对 k≤100 验证 dk、ek 下界为正，耗时约 1 分钟，从而保证至多 202 个元件的梯形全为正值与无源。
第 3 类是拟合与约化计算：宽带终端只在对数坐标下拟合末节两个元件以最小化宽带最大误差，Foster 型拟合用向量拟合找极点再精修，因非凸而多起点取最优，所报误差是上界而非全局最优。缺失的信息是全 k 正性的解析证明与 M 连分式全正系数的存在性，论文明确标为猜想与未来工作。无训练不等于确定性求解，高精度制表只保证本次系数正确，不保证双精度下重算仍正确。
误差口径频段与号筒算例的条件
活塞实验默认用绝对误差|Z-Zexact|/ρ0cS，只在汇总表同时给出相对误差以照顾极低频段。嵌套梯形考察 ka≤2、5、10、20 的最大误差与误差低于 1% 的最大 ka；宽带比较统一用 0.1≤ka≤50；时域考察单位速度阶跃力在较晚时刻的最大偏离与初值的振荡。号筒算例取喉半径 25.4mm、长 0.8m、口半径 30cm 的指数号筒，截止频率 168 Hz，面积比 140，口 ka 在 50 到 5000 Hz 对应 0.27 到 27，参考解是指数段 Webster 方程闭式传输矩阵，终端口负载为精确活塞。
喉误差以 ρ0c/St 归一化，区分模型误差与近似误差：平面波刚性无损耗指数号筒加障板圆活塞只是模型，其与真实号筒的差距不在评估之列；离散段数与口网络引起的偏离才是方法误差。热黏壁损失经复波数估算在截止附近最大 0.013，以上对比均为无损耗，条件不一致时不能直接比较大小。
嵌套梯形在低频的精度如何随阶数增长
嵌套梯形是 s=0 处对角单点帕德的网络形式，每元件多匹配一个泰勒系数，因此极低频精度极高而截止频率随元件数线性外移。泰勒多项式虽处处收敛但截断后数值完全不可用，说明收敛不等于实用。下表比较同为嵌套截断、同为绝对误差口径时，加阶带来多少低频精度与可用带宽，比较问题是低频段内加阶是否系统改善，指标越小越好。
| 条件 | 元件数 N | 最大绝对误差 | 上限 ka | 误差口径 |
|---|---|---|---|---|
| 嵌套低频 | 14 | 4×10-4 | ka=5 | 绝对误差 |
| 嵌套低频 | 26 | 1.1×10-4 | ka=10 | 绝对误差 |
同样以精确活塞为参考、同一绝对误差口径计算，14 元件在 ka 达 5 时误差已很小，26 元件把同样量级的误差推到 ka 达 10，误差随 N 近似指数下降而带宽线性增长。该对比只反映单点展开的低频优势，不能外推为宽带最优，宽带还需看终端与拟合的对比，高 ka 段仍受波纹限制。
嵌套 × 宽带终端： 嵌套指加阶时只在开路末端追加新节，已有元件值不变，便于按需提高精度而不重做电路；宽带终端指保留嵌套主体，只重算或拟合最后 1 个电阻或最后两个元件以兼顾高频极限与中高频平均误差。两者分工是嵌套保证构造的系统性与无源性，终端用最小的非嵌套代价换取宽带误差明显下降，组合后 14 元件宽带误差从 0.068 降到 0.023。
下面这组阻抗与误差曲线是官方原图，左面板灰色为精确值，彩色虚线为不同元件数，右面板为误差随 ka 的变化，读图时先看低 ka 重合段再看截止点右移。
看图路径： 1. 左图比较 8、14、22 元件曲线与灰色精确曲线的贴合范围；2. 右图沿横轴 ka 增大，观察不同元件数误差曲线的截止点右移；3. 注意高 ka 段误差不再单调下降而是受波纹限制
论文图 7。原论文 Figure 7:：“Left: exact impedance (grey) and nested ladders of 8, 14 and 22 elements (dashed). Right: absolute error for ladders of 6 to 50 elements.”。
图中可见低 ka 段虚线与灰线重合，高 ka 段先出现过冲偏差再进入波纹平均；右图每条曲线在截止前陡降、截止后形成平台，平台高度约等于当地波纹幅度。未胜出项是泰勒截断，尽管理论收敛但数值完全不可用，这提醒不要把收敛域当成精度保证。
宽带终端与四元件家族带来多少改善
低阶时单点梯形高频极限错误，因此需两点思想。2 阶闭式展开为 M0 并联 Ra 与 Rb 串 C，Ra 复现低频电阻，Ra 并 Rb 固定高频极限。该拓扑与把嵌套末电阻按高频极限重算的基网络以及 Beranek 网络同属单参数家族，仅时间常数不同。Beranek 电容与两点值相差约 2%，可视为 2 阶两点帕德的图形精度实现。下表比较同为 14 元件、同为宽带最大误差时，保留主体只改末节带来多少改善，指标越小越好。
| 条件 | 元件数 N | 嵌套误差 | 终端 a 误差 | 终端 b 误差 |
|---|---|---|---|---|
| 宽带 14 元件 | 14 | 0.068 | 0.046 | 0.023 |
同样在宽带条件下，终端 a 用一个非嵌套电阻换来约三成改善，终端 b 用两个非嵌套元件换来约 3 倍精度，但剩余误差主要是衰减波纹，有限阶有理函数只能跟随至某 ka 而后取平均。终端 b 破坏完全嵌套，且中高频下限由波纹幅度钉住，加阶只能缓慢下降。下表比较同为四元件、同一宽带口径时，不同时间常数选择的位置，指标越小越好。
| 条件 | 口模型 | 宽带误差 | 元件数 N | 误差口径 |
|---|---|---|---|---|
| 四元件族 | Beranek | 0.155 | 4 | 绝对误差 |
| 四元件族 | 基网络 | 0.187 | 4 | 绝对误差 |
| 四元件族 | 终端 b | 0.102 | 4 | 绝对误差 |
同样四元件预算下，基网络低频最好，两点与 Beranek 宽带接近，极小极大最优更低但低频稍差。更低宽带误差需要放弃高频极限的精确固定，这是论文明确指出的代价。
拟合方法与时域表现是否改变结论
若不需要电路，直接近似斯特鲁夫函数的三常数公式加精确实部在全范围达 0.14%，远优于任何网络，这是分解出已知结构这一常用补救的胜利。需要电路时，Foster 型拟合在宽带以 13 元件达 0.014，优于同尺寸嵌套终端，但在低频以 14 元件两者相当，而 22 元件时梯形远超拟合，说明拟合低阶宽带占优、嵌套高阶低频占优。下表比较同一时域阶跃力、同一较晚时刻最大偏离时，加阶带来多少时域改善，指标越小越好。
| 条件 | 元件数 N | 时域最大偏离 | 时刻条件 | 误差口径 |
|---|---|---|---|---|
| 阶跃 8 元件 | 8 | 0.048 | 较晚时刻 | 绝对偏离 |
| 阶跃 14 元件 | 14 | 0.018 | 较晚时刻 | 绝对偏离 |
| 阶跃 22 元件 | 22 | 0.008 | 较晚时刻 | 绝对偏离 |
同样以精确阶跃响应为参考，加阶使时域偏离系统下降，有限支撑被密切近似但不能精确实现，初值在 1 附近振荡。未胜出项是 AAA 与 Thiele，前者不保证无源，后者依赖节点选择，都不是嵌套无源梯形。号筒方面均匀管级联与集总 T 节都以 2 阶收敛，但后者常数大 10 倍，需更多节数才达同样喉误差。喉误差几乎等于口误差，完整被动梯形规模很大，虽可进电路仿真但需权衡规模。
哪些断言尚未被证明与尚未评测
已证部分包括连分式与梯形结构、M0 与相关闭式、四元件家族与 2 阶两点闭式，以及经区间算术认证的 k≤100 即 202 元件内正系数从而无源。数值验证部分包括各频段误差、带宽拟合与号筒传递，拟合误差因非凸只是上界，更长搜索可能更低但不会更高。猜想部分是全 k 正性，常用正测度表示因等效 Jacobi 对角为虚数而不直接适用，需对渐近给出显式余项界，或经延迟积分给出倒数阻抗的矩表示。
方法局限有三：展开点在 s=0 故高频极限未强制，终端以一两个非嵌套元件补救，完全嵌套的两点 M 连分式仍待构造；Viskovatov 递推病态，高阶必须高精度一次性制表，双精度会产生虚假负系数；无限指数号筒在通带截断不收敛，需以精确尾阻抗终端。未评测范围包括互阻抗、高阶斯特鲁夫量、模态耦合与有限障板衍射，以及实测喉阻抗对平面波与活塞假设的检验，这些都超出本篇无损耗模型的声称范围，缺证据不等于技术有错。
复现第一步与关键参数如何取
复现应从活塞幂级数出发，先用精确有理数或上 1,000 位浮点实现 Viskovatov 得到 C 连分式，再收缩成正系数并换算元件值，计数含辐射质量与电阻，奇节串联 C-R、偶节并联 M-R。第一步只做至 14 个元件并在 ka≤5 检查绝对误差是否达相应量级，同时检查时域较晚时刻偏离是否约 0.018；通过后再试宽带终端重算末电阻使高频极限为 1，或在对数坐标拟合末节两元件最小化宽带最大误差。关键信息条件是保留归一化换算、误差口径与频段：质量与顺性电阻各有对应乘子。
比较拟合与嵌套时必须同频段同口径。代码方面论文称所有数字表格由随附包产生，但本次未收到可验证的作者资源链接，不得声称代码已公开或可下载。计算预算主要是高精度构造的一次性成本，推理与仿真本身只是普通电路求解，不宜把制表成本当成每次仿真延迟。
什么时候值得尝试与还需补什么验证
当需要把障板活塞辐射直接放进电路或状态空间、且希望加阶不重做电路并在每阶都无源时，值得尝试嵌套梯形；当工作频段明确上限如 ka 达 5 且要极高低频精度时，14 元件已足够；当只有四元件预算且要兼顾高频极限时，应在单参数家族中按时间常数选择，两点值与 Beranek 值几乎重合，极小极大最优更低但低频稍差；当宽带平均误差优先且可接受非嵌套时，可用末节拟合或 Foster 型拟合，但须接受多起点非凸搜索与无源需额外保证。
号筒应用中若口 ka 覆盖较低到二十多且用平面波级联，应先保证口网络精度，因为喉误差几乎等于口误差，离散段数另按 2 阶收敛选取。还需补的验证是全 k 正性的解析证明或显式余项界、完全嵌套两点连分式的正系数构造，以及用实测喉阻抗检验平面波与活塞口假设本身。常见误解是把区间认证的 202 元件内无源当成对任意阶的证明，把拟合上界当成可部署保证，或把对模型的近似误差当成对真实号筒的预测误差。
论文把精确闭式、认证区间、数值验证与猜想分开陈述，复述时也应保持这种分层。
喉阻抗 × 口阻抗： 喉阻抗是号筒小端看进去的声阻抗，决定驱动器负载；口阻抗是号筒大端向外的辐射阻抗，本篇取障板圆活塞阻抗。分工是号筒级联把口阻抗经逐段传输矩阵变换到喉端，相当于以口阻抗为尾的连分式，组合后口网络的误差几乎一对一传递到喉阻抗，因此口模型必须先做好，号筒离散本身另有二阶收敛误差。
📎 论文与评分元数据
排名：前50% | 文档类型：方法研究 | arXiv 原文
⚖️ 评分明细
评分属于系统判断，不是论文实验结果；八维数值与总分见页首，原始审计记录保留在后端。
评分规则：type-aware-v1
评分模型：muse-spark-1.3-contributor
评分请求协议：openai_responses