企业官网建设流程全解析

做控制系统的朋友，尤其是搞过电源小信号建模、电机驱动或是在MATLAB里做系统辨识的，一定有过这种瞬间：手里有一个传递函数，想看看它的增益特性，有人说“看中频增益就行”，又有人说“化成尾一标准型，直接读系数K就是增益”。结果你按后一种方法做了，拿到一个数值，和Bode图上明显对不上，当场就懵了。最近就有个做DAB电源的工程师拿着同样的问题来找我，问“中频增益”和“尾一标准型里的K”到底是不是一个东西。这问题问得特别好，因为它是控制理论里一个非常典型、又很容易被讲混的概念。

先给结论：严格说，它们不是同一个东西。“中频增益”描述的是系统幅频特性在某个频段内的代表值，是一个随频率定义域变化的响应值；而“尾一标准型”里的系数K，是传递函数做代数恒等变形后从分子分母常数项里提取出来的一个放大倍数。二者只在某些特定系统里数值碰巧相等，但定义、物理含义、适用场景全都不一样。这篇文章我会把这两个概念彻底掰开揉碎，再用MATLAB的tfest辨识流程给你演示一遍，看完你就能明白它们什么时候相等、什么时候八竿子打不着。

1. 先弄明白：中频增益到底指哪个数值

1.1 频率分三段，“中频”其实是个工程概念

要想搞清中频增益，先得接受一个前提：控制理论里的“低频”“中频”“高频”并没有像数学公式一样精确的分界线。它更像是在说一件事——把频率轴从低到高拉出来，系统的幅频特性通常可以分成三段来看：

• 低频段：一般指远低于第一个转折频率的区域。这个区域里，增益渐近线往往是水平的，数值上等于系统的直流增益（如果没有积分环节），或者以 -20dB/dec、-40dB/dec 的斜率上升/下降（如果有积分/微分项）。
• 中频段：指的是转折频率集中、相位变化最剧烈的区域。对大多数控制系统来说，这就是穿越频率（0dB交点）附近以及相位裕度起作用的频带。在这个区间，增益曲线可能平坦，也可能快速爬升或下降。
• 高频段：远高于最后一个转折频率的区域，受系统阶次和零极点分布影响，增益通常按 -20dB/dec×（极点个数-零点个数）的斜率衰减。

所以“中频”不是一个点，而是一个频带。中频增益，就是系统在这个频带内的“代表增益值”。但问题来了，既然是一个带而不是一个点，取哪个频率上的幅值才算数？这没有绝对标准，实际操作中一般看你想用这个增益去衡量什么：如果是想评估稳定性，就看穿越频率ωc处的增益；如果是设计补偿器，就看平台段的增益；如果是分析某个带通滤波器，就看通带内的平均增益。

1.2 中频增益是“一段频带的代表值”，不是某个代数常数

“代表值”这个词很关键。它意味着中频增益本质上是一个频率响应幅值，记作 |G(jω)| 在某段频率区间内的取值。它带有频率属性，是真实物理系统在正弦激励下的输出/输入幅值比，单位可以是倍数，也可以换算成 dB。

举个例子。一个带通滤波器，中心频率1kHz，通带增益20dB，那么它的中频增益就可以说是20dB（约10倍）。这个20dB并不是传递函数分子上随便一个常数，而是当ω落在通带内时，|G(jω)| 真实达到的高度。换一个频率范围，比如远离通带的100Hz或100kHz，增益可能就只有-20dB甚至更低。

而对于低通系统，比如一个典型的一阶惯性环节 G(s)=K/(Ts+1)，它在低频段是平台K，过了转折频率1/T后开始以-20dB/dec下降。你能说它的“中频增益”是K吗？严格说不能，因为一阶系统根本没有一个像样的、相位变化集中的“中频平台”，它的幅频曲线从低频平台直接过渡到高频衰减，中间只是一个转折点。习惯上我们叫K为低频增益，不叫中频增益。如果非要在Bode图上找一个和“中频”沾边的值，那大概就是转折频率附近的增益，约等于 0.707K，也就是-3dB点。这个值显然不等于K。

1.3 典型例子：带通型网络的平台增益

为了让你对“中频增益”这个概念有画面感，我直接写一个典型的三阶带通型传递函数：

G(s) = K * (Tzs + 1) / [(Tp1s + 1)(Tp2*s + 1)]

其中令 Tp1 = 1/ωp1，Tz = 1/ωz，Tp2 = 1/ωp2，并且 ωp1 < ωz < ωp2。也就是说，第一个极点先让增益下降，零点把下降的趋势挡停，第二个极点再让增益继续下降。这样幅频曲线的形状就是：低频平台→下降→中频平台→下降，形成一个“凸台”。这个凸台的高度就是中频增益的典型体现。

用渐近线近似算一下，中频平台高度约等于 K * ωp1 / ωz。如果K=10，ωp1=1000rad/s，ωz=10000rad/s，那么中频平台增益约等于1，也就是0dB。你看，K明明是10，但中频增益只有1。这直接说明：K并不等于中频增益，中频增益取决于K与零极点频率的相对位置。

这时候再回头看“中频增益”的定义，你就明白了：它是在一段频率范围内幅频特性的实际高度，和系统零极点布局强相关，不是一个纯粹的代数系数。

2. 尾一标准型里的系数K是怎么来的

2.1 从“首一标准型”到“尾一标准型”

控制理论入门时，大家接触最多的其实是“首一标准型”，也就是把分母最高次项的系数化为1。比如二阶系统标准形式：

G(s) = ωn² / (s² + 2ζωn s + ωn²)

这就是首一标准型，因为分母最高次项s²的系数是1。这种形式的好处是高阶项系数干净，方便分析极点位置和稳定性。

而“尾一标准型”恰好相反，它把分母多项式常数项化为1，使得所有与时间常数相关的项都变成“某某s + 1”的形式。同一个二阶系统，分子分母同除以ωn²，得到：

G(s) = 1 / (1 + 2ζ/ωn * s + s²/ωn²)

再写成更一般的形式，就是：

G(s) = K * (τz1s + 1)(τz2s + 1) / [(τp1s + 1)(τp2s + 1)]

因为每个因子尾巴上的常数项都是1，所以叫“尾一”。这种形式在工程上特别受欢迎——每个τ都直接对应一个转折频率（τ=1/ω），你能一眼看出系统有哪些低频、高频拐点，做Bode图几乎不用动脑子。它有时也被称为“时间常数标准型”或“Bode标准型”，但本质是同一件事。

2.2 K的本质：零频增益、积分增益，还是都不是？

尾一标准型里的系数K是在做恒等变形过程中提取出来的常数因子。它到底代表什么，取决于分母里是否含有积分因子s。

第一种情况，系统分母没有s因子（也就是没有积分环节）。此时如果令s→0，你会发现所有(τ*s+1)都趋近于1，于是：

G(0) = K

所以K就是系统的直流增益，也就是低频渐近线的高度。例如一阶惯性环节 G(s)=K/(Ts+1)，K既是尾一标准型系数，也是直流增益。

第二种情况，系统分母里有积分环节，也就是s的因子。比如一个典型的Type-III补偿器：

Gc(s) = K * (1 + s/ωz1)(1 + s/ωz2) / [s * (1 + s/ωp1)(1 + s/ωp2)]

注意这里分母有一个单独的s，所以不能把分母常数项直接化为1。通常我们把式子写成K/s乘以剩余尾一形式。这时候K有什么物理含义？它不是直流增益（直流增益是无穷大，因为积分环节在s=0处有极点），也不是中频增益，而是决定积分增益大小的一个系数。你可以把它理解为“在ω=1rad/s附近，系统增益大致等于K”的基准参数，但这个基准与中频平台的真实高度之间隔着零点和极点频率的拉扯，不能直接画等号。

第三种情况，分子有零点在原点，比如包含微分环节。此时分子常数项为0，尾一标准型写不出来，需要把s因子提取到K里面变成K*s，K的单位就会变成与s相关的量纲。这种情况下的K既不是直流增益，也不是中频增益，只是一个比例系数。

2.3 尾一标准型并不是“唯一标准型”，K会随写法而变

很多人忽略了一个细节：尾一标准型不是唯一的。它只是把“常数项为1”作为约定，但你可以选择把哪个因子吸收进K里，写法不同，K的数值就不同。

举例说明。系统 G(s) = (0.1s+1)/(0.01s+1)(0.001s+1)，尾一化之后K=1。这是标准做法。但如果我偏要把分子写成 10*(0.01s+0.1) 这种不伦不类的形式，K就会变成10。虽然传递函数没变，K却变了。所以K作为代数系数，本质上依赖于你采用的“标准化约定”。

这一点很重要，因为在实际辨识工作中，tfest返回的是分子分母多项式的数值系数，这个系数可能是首一形式的，也可能是其他归一化形式的。如果你不先做标准化处理就直接把某个系数当“增益”，很容易出笑话。正确做法是：先提取零极点，或者按照“分母常数项化为1”的规则重新整理，再读K。

3. 中频增益和系数K，到底是不是一个东西

3.1 结论先行：只有在特定条件下数值才可能相等

现在可以正面回答最开始那个问题了。

严格说，中频增益是一个频域响应值，它等于|G(jω)|在某个中频频带内的取值；尾一标准型的K是一个代数系数，它是传递函数在标准化约定下的比例因子。两者定义域完全不同，一个是函数在特定自变量下的值，一个是函数表达式的参数。

它们什么时候相等？只有一种常见情形：系统在“低频段到中频段”没有明显转折，且低频平台一直延伸到你要考察的中频频带内。也就是说，如果幅频特性从ω=0到你关心的频段都是平坦的，那么这个平台高度既等于K，也可以勉强称为该频段内的增益。比如一个截止频率远高于工作频带的惯性环节，你在工作频带内测到的增益就是K。

但控制系统的“中频段”通常伴随着转折和相位变化。一旦有零极点落在频带内，中频增益和K就分道扬镳了。下面我用几个例子让你彻底记住。

3.2 用二阶低通系统说明：谐振峰增益≠K

二阶低通系统的首一标准形式是：

G(s) = ωn² / (s² + 2ζωn s + ωn²)

化成尾一标准型：

G(s) = 1 / (1 + 2ζ/ωn * s + s²/ωn²)

此时K=1。直流增益是1（0dB），这个没问题。但如果你把阻尼比设小一点，比如ζ=0.2，Bode图上在ωn附近会出现一个明显的谐振峰。谐振峰的幅值用公式算：

Mr = 1 / (2ζ√(1-ζ²))

代入ζ=0.2，Mr≈2.55，换算成dB大概是8.1dB。这个谐振峰就落在“中频段”附近，你可以把它看成系统在转折频率附近的增益表现。你看，K=1（0dB），中频附近的实际增益却是8.1dB，差了8dB。这时候你还能说系数K就是中频增益吗？显然不能。

这个例子特别经典，因为它说明了一个普遍规律：K给出的是渐近线的基准，而真实幅频响应会被零极点的Q值（阻尼比）抬高或压低。中频增益是“真实响应”，K只是“代数骨架”。

3.3 DAB电源补偿器里的实战辨析

再把视角拉回到电力电子领域。DAB（双有源桥）变换器的电压环通常用Type-III补偿器，传递函数我前面写过：

Gc(s) = K * (1 + s/ωz1)(1 + s/ωz2) / [s * (1 + s/ωp1)(1 + s/ωp2)]

型如这种补偿器，K虽然被很多人称为“增益系数”，但它和环路的中频增益完全不是一回事。环路设计时，我们真正关心的是穿越频率ωc处的开环增益——通常要求它等于1，也就是0dB。为了实现这个目标，K该定多少？需要代入ωc，计算|Gc(jωc)|，再乘上被控对象在ωc处的增益，反解出K。也就是说，K是根据“目标中频增益”反推出来的结果，而不是中频增益本身。

我在调试一个DAB数字电源时遇到过类似的事。仿真扫频得到被控对象在5kHz处的增益是-8dB，Type-III补偿器零点设在1.5kHz和800Hz，极点设在20kHz和80kHz，我一开始随便设K=10000，结果Bode图上穿越频率跑到12kHz，相位裕度只剩15度。后来老老实实按目标穿越频率反算K，才把带宽压回设计值。那一次我彻底记住：K只是数学上的比例因子，中频增益是Bode图上的真实值，两者隔着零极点频率的“搬运”过程。

3.4 区别对照表

为了让你以后查阅方便，我把这两个概念放在一张表里对比：

4. MATLAB实操：从tfest辨识到提取中频增益与K

4.1 先用delayest估计纯滞后时间τ

聊完概念，来点实际能用的。很多帖子都在问“用tfest辨识完传递函数后，怎么确定纯滞后时间τ”。这里我给出我的标准流程：不要等辨识完再猜τ，而是用数据先估计延迟，再把延迟作为已知或待估参数交给tfest。

第一步，拿到输入输出数据后，构造iddata对象，然后用delayest做一个快速估计：

data = iddata(y, u, Ts); % y和u是列向量，Ts是采样周期 nDelay = delayest(data); % 估计延迟，结果以采样周期为单位 tau0 = nDelay * Ts; % 换算成秒 fprintf('估计纯滞后: %d 个采样周期, 即 %.6f 秒\n', nDelay, tau0);

delayest的内部原理并不复杂，它本质上是在输入输出数据之间做相关性分析，寻找让输入解释输出效果最好的时间偏移量。这个方法对含噪数据有一定鲁棒性，可以给我们一个不错的τ初值。但如果系统本身动态较慢、噪声较大，这个初值可能会有几个采样周期的偏差，没关系，下一步我们可以把τ交给tfest继续精修。

还有一种思路：如果数据是频域响应数据，可以通过相位曲线的线性下降斜率来估计τ。因为纯延迟e^{-τs}的相位是-ωτ，在相位曲线上表现为一条随频率线性下降的线，斜率就是-τ。你可以把无延迟模型（比如初步辨识的模型）的相位与原数据的相位相减，再用最小二乘拟合剩余相位的斜率，换算成τ。这个方法在扫频数据质量好的时候很准，但要注意相位缠绕（unwrap）问题。

4.2 用tfest辨识传递函数（含延迟）

有了τ的初值以后，就可以做正式辨识了。tfest可以直接把纯延迟作为待估参数，做法是给InputDelay传递一个搜索区间。我比较常用的是在tau0附近留出正负几个采样周期的余量：

np = 3; % 期望的极点个数，按你的系统阶次来 nz = 2; % 期望的零点个数 % 方式1: 让tfest在τ的区间内自动搜索 tauLow = max(0, tau0 - 5*Ts); tauHigh = tau0 + 5*Ts; sys1 = tfest(data, np, nz, 'InputDelay', [tauLow, tauHigh]); % 方式2: 固定τ，只辨识零极点，更稳但依赖τ初值准度 sys2 = tfest(data, np, nz, 'InputDelay', tau0);

方式1的好处是τ也会被优化，拟合度通常更好；缺点是如果数据噪声大，延迟和零极点参数之间可能存在耦合，偶尔会收敛到奇怪的局部最优。所以我通常先用方式1快速看结果，再把得到的延迟四舍五入成采样周期的整数倍，用方式2固定延迟重新辨识一遍，对比两者的拟合度（可以用compare(data, sys)看），取更优的。

如果想直接估计连续时间模型，可以在tfest里加一个选项：

sysc = tfest(data, np, nz, 'InputDelay', [tauLow, tauHigh], 'Ts', 0);

设定Ts=0表示让MATLAB在连续时间域里拟合传递函数。做完之后，sysc就包含零极点、增益和纯延迟，可以用bodeplot、step、margin等工具继续分析。

4.3 把辨识结果化成尾一标准型，提取K

辨识得到的模型默认是多项式传递函数，也就是分子分母按s降幂排列。想化成尾一标准型，我的做法是用tfdata提取分子分母系数，然后统一除以分母常数项，再把分子常数项部分提取为K。

[num, den] = tfdata(sysc, 'v'); % 将分母常数项化为1 denTail = den / den(end); numTail = num / den(end); % 分子常数项就是K（前提是没有积分环节） K = numTail(end); % 提取K后，剩余分子常数项为1 numTail = numTail / K; sysTail = K * tf(numTail, denTail); fprintf('尾一标准型系数K = %.6g\n', K);

注意，这里的sysc如果是用方式1辨识的，还带InputDelay属性。尾一标准型通常不包含延迟，所以你要单独记录一下纯延迟：

tauFinal = sysc.InputDelay;

如果你的系统有积分环节，比如分母常数项本来就是0，那直接用上面这段代码会报错。这时需要先把系统转为zpk形式，手动把s因子提取出来，再对剩余部分做尾一化。简单来说就是：

syszpk = zpk(sysc); % 观察zero/pole里是否有s=0或s=∞的因子

辨识模型里有积分环节的情况并不罕见，尤其是处理含内环的开关电源控制对象时。遇到这种情况不要慌，提取积分因子后，K的意义就变成积分增益的系数了，不再是直流增益。

4.4 在Bode图上算中频增益

最后一步，计算中频增益。这取决于你关心的频段是什么。如果你已经确定了穿越频率ωc，那就直接算这个频率点的增益：

% 计算某个频率点的增益 wc = 2*pi*5000; % 以5kHz为例 Gjw = evalfr(sysc, 1i*wc); midGain = abs(Gjw); midGain_dB = 20*log10(midGain); fprintf('在 %.2f kHz 处的增益为 %.2f dB\n', wc/2/pi, midGain_dB);

如果你想看整个中频段内的增益变化曲线，然后用鼠标在Bode图上直接读，可以这样：

figure; bodeplot(sysc, {2*pi*100, 2*pi*20000}); grid on;

我更推荐在程序里先把频带内的幅值拉出来，然后取平均值或者找平台段，这样便于批量处理数据：

wMid = logspace(log10(2*pi*1000), log10(2*pi*10000), 500); mag = squeeze(abs(freqresp(sysc, wMid))); mag_dB = 20*log10(mag); % 中频段平均增益 midAvg_dB = mean(mag_dB); % 中频段最大增益（如果想看谐振峰） midPeak_dB = max(mag_dB);

你会发现，这样计算出来的中频增益，和前面提取的K在大多数情况下不是同一个数值。尤其是当Bode图上有明显拐点或谐振峰时，K只是低频基准或代数系数，而中频增益才是你真正该关心的环路设计指标。

5. 常见误区与避坑清单

5.1 误区一：以为“K=直流增益”在所有情况下都成立

这是最普遍的误解。尾一标准型系数K只有在系统没有积分环节的前提下才等于直流增益。分母一旦含有s因子，直流增益要么是无穷大（分母含积分），要么是0（分子含积分），这时候K代表的只是某个频段的比例基准。你在MATLAB里输出一个很大的K，别急着说“系统增益很高”，先看有没有积分项再说。

5.2 误区二：拿谐振峰当“中频平台”

有些系统的幅频曲线在中频段不是平坦的，而是有一个尖锐的谐振峰。比如你前面看到的ζ=0.2的二阶低通系统，在ωn附近有8dB左右的峰值。这时候如果你把谐振峰高度当成“中频增益”，设计环路时很容易高估增益，导致实际穿越频率比预期高很多。正确的做法是：定义清楚你评估的频带范围，一般取穿越频率周围0.5~2倍频程，看平均增益而不是峰值增益。如果曲线起伏太大，说明系统Q值太高，优先考虑加阻尼或改补偿器。

5.3 误区三：忽略纯滞后对相位的影响，只盯着幅值增益

纯延迟e^{-τs}的幅值恒为1，所以在Bode幅频曲线上它完全不改变增益形状。于是很多人辨识完模型，看到幅频曲线拟合得很好，就忽略了τ的准确性。问题是，纯延迟对相位的影响是线性的，τ一个采样周期的偏差，在高频段就可能带来几十度的相位差。相位错了，环路稳定性分析全错。所以辨识时一定要像第4节那样认真对待τ，而不是只看幅频拟合。

5.4 误区四：以为tail-1形式只有一种

前面说过，尾一标准型依赖于标准化约定。你从论文里看到某个“K”和从MATLAB辨识出的“分子第一项”不见得是同一个值。为了避免交流歧义，建议每次给出尾一标准型时，都把分子写成∏(τzs+1)的形式，分母写成s^N * ∏(τps+1)的形式，明确写出K的位置。这样别人复制你的模型时，不会因为归一化方式不同而算错。

5.5 实操速查表

6. 最后说点个人体会

我在刚接触控制理论那两年，也被K和“增益”这两个词坑过。当时做一个Buck变换器的建模，把传递函数化成尾一标准型，看到K=12000，脑子里第一反应是“这系统增益真高”，后来用Bode图一看，低频确实是41.6dB，但5kHz附近早就掉到个位数了，哪里还有什么12000。那时候才明白，传递函数是一个函数，不是一个数，它的“增益”天然就是随频率变化的。

现在我做环路设计时给自己立了一个规矩：凡是涉及“增益”的讨论，一律先在Bode图上标注一个具体频率点或频段，然后再说那个值。比如“穿越频率30kHz处的增益是-1.2dB”，而不是笼统地说“增益是K”。K只作为代数系数和中间计算量存在，它帮助我快速搭建渐近线、估算零极点位置、反推补偿器参数，但它永远代替不了幅频曲线本身的含义。你手上有MATLAB也好，用公式推导也好，最可靠的做法始终是：化尾一标准型得到K，然后回到幅频曲线上重新确认中频增益。两者对上了，说明你的系统足够简单；对不上，说明正是那些转折点上的零极点，塑造了系统的真实动态。

这个小习惯，帮我少走了很多弯路。希望这篇文章也能帮你把这两个概念彻底理顺。