基于FFT曲线在MATLAB中辨识传递函数的方法验证与阶次疑问

用户持有如下FFT曲线：

希望在MATLAB中从该曲线辨识传递函数，已执行以下代码：

np = 2; % 极点数量
nz = 1; % 零点数量
iodelay = 0; % 传输延迟
Ts = 0.1; % 采样时间 [s]
input = freq'; % FFT曲线的x轴频率数据
output = xdft'; % FFT曲线的y轴频域响应数据

sysd = tfest(input,output,np,nz,iodelay,'Ts',Ts);
[num, den] = tfdata(sysd);
num = cell2mat(num);
den = cell2mat(den);
transfer_function = tf(num,den);

得到的传递函数如下：

现咨询两个核心问题：

1. 当前操作是否合理正确？
2. 在未知原传递函数的情况下，如何选择零点数量？

一、当前操作的合理性分析 #

1. 整体流程符合规范：使用MATLAB的tfest函数从频域数据辨识传递函数是系统辨识领域的标准方法之一，你的代码逻辑框架没问题——导入频域数据、指定模型阶次与采样时间、拟合得到离散传递函数、提取分子分母，这个流程是正确的。
2. 需要验证的细节点：
   • 频率单位匹配：确认freq是**角频率（rad/s）**而非普通频率（Hz），因为tfest默认输入为角频率，若freq是Hz，需先乘以2π转换，否则拟合结果会出现系统性偏差。
   • 采样时间的准确性：Ts=0.1s必须与实际数据采集的采样间隔一致，或与FFT的频率分辨率匹配（频率分辨率Δf=1/(N*Ts)，N为采样点数），错误的Ts会导致离散传递函数的动态特性失真。
   • 延迟项的假设合理性：从FFT的相频曲线看，低频段无明显线性相位滞后，假设iodelay=0是合理的，但后续可通过残差分析或相位拟合验证是否存在隐藏延迟。
   • 拟合质量验证：建议添加compare(sysd, input, output)命令，对比辨识模型的频域响应与原始FFT数据，查看拟合误差；同时关注tfest输出的拟合优度（Goodness of Fit），通常该值超过80%说明拟合效果可接受。

二、未知原模型时的零点数量选择方法 #

在没有真实传递函数参考的情况下，零点数量可通过以下方式确定：

• 从频域曲线特征判断：
  • 幅频曲线：若幅频曲线出现明显凹陷（ notch），通常对应零点；如果是单调衰减趋势，零点数量可能为0或1。你的FFT幅频曲线在中频段有小波动，选择1个零点是合理的尝试。
  • 相频曲线：零点会带来相位超前，极点带来相位滞后。若相频曲线在某频段有明显的相位超前趋势，说明需要增加零点数量来匹配。
• 逐步试探+指标评估：
  先固定极点数量（可从系统动态特性预估，比如阶跃响应的振荡次数、频域截止频率），从零开始逐步增加零点数量，每次拟合后对比拟合优度和残差平方和。当增加零点后，拟合优度提升幅度小于5%、残差不再显著减小时，停止增加零点——遵循“奥卡姆剃刀”原则，用最简模型达到足够精度。
• 使用自动阶次选择工具：
  打开MATLAB的ident图形界面，导入频域数据后，利用Model Selection工具自动尝试不同阶次的模型，对比各模型的AIC（赤池信息准则）、BIC（贝叶斯信息准则），这两个指标越小，说明模型在拟合精度与复杂度之间的平衡越好。
• 结合系统物理意义：
  若知道系统的物理结构（如机械系统、电气系统），可根据原理预估零点数量。比如RLC电路可能存在1个零点，带减速器的伺服系统通常零点数量为0或1。

内容的提问来源于stack exchange，提问作者Jonathan Frutschy