如何使用MATLAB估算无零点2极点传递函数的m、d参数

你之前调用tfest报错的核心原因是：tfest是用于输入输出配对的动态系统辨识的函数，不能仅传入输出端的位置时间序列，需要同时提供输入推力数据、输出位置数据以及采样时间，构造标准的辨识数据对象。

完全适配你需要的连续域模型需求，步骤如下：

• 首先准备3组基础数据：
  • 输出位置序列x（就是你已经有的数据，长度为N）
  • 对应时间点的输入推力序列u（和x长度一致，每个时刻给系统的推力值）
  • 采样时间间隔Ts（你采集数据的时间步长，单位为秒）
• 构造MATLAB辨识专用的iddata对象：

data_id = iddata(x, u, Ts);

• 调用tfest辨识，同时约束传递函数分母常数项为0（匹配你的模型结构1/(m s² + d s)的特性）：

% 设定传递函数结构：2阶分母，0阶分子，无直通项，强制分母常数项为0
sys = tfest(data_id, 2, 0, 'Feedthrough', false, 'DenominatorConstant', 0);

• 提取参数：辨识得到的分母为[m_est, d_est, 0]，直接取前两项就是m和d的估计值

den_coef = sys.Denominator{1};
m_est = den_coef(1);
d_est = den_coef(2);

因为你已经明确知道系统的物理模型结构，用灰盒辨识可以避免黑盒辨识出现不符合物理意义的结果，直接输出m和d的估计值：

• 第一步定义模型函数，保存为friction_mass_model.m：

function [A,B,C,D] = friction_mass_model(params, Ts)
    m = params(1);
    d = params(2);
    % 连续状态空间模型，对应传递函数1/(m s² + d s)
    A = [0 1; 0 -d/m];
    B = [0; 1/m];
    C = [1 0];
    D = 0;
    if Ts > 0 % 采样时间不为0的话自动离散化
        [A,B,C,D] = c2d(A,B,C,D,Ts);
    end
end

• 第二步构造灰盒模型对象，设定参数初始猜测值：

% 初始猜测m=1，d=0.1，可以根据你的物理场景调整初始值加快收敛
params0 = [1, 0.1];
% 设定参数为非负（物理上质量、阻尼系数不可能为负）
lb = [0, 0];
ub = [inf, inf];
% 创建连续时间灰盒模型
sys_grey = idgrey(@friction_mass_model, params0, 'c', [], lb, ub);

• 第三步用之前构造的data_id对象估计参数：

sys_est = greyest(data_id, sys_grey);
% 直接提取估计得到的m和d
m_est = sys_est.Parameters(1).Value;
d_est = sys_est.Parameters(2).Value;

• 如果你的推力输入是阶跃信号（比如从0时刻开始加恒定推力），辨识效果会更好，尽量保证输入信号有足够的激励性
• 如果没有采集到输入推力数据，需要先做盲辨识，但精度会远低于有输入数据的情况，不推荐

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