如何估计具有 2 个极点且无零点的传递函数的参数？

Question

我是 matlab 新手，我有一个传递函数 1/(ms^2+ds) 用于表示经历推力、阻力和重力的物体的位置。我想估计 m 和 d。我有时域的输出数据，但我不确定如何使用 matlab 估计 m 和 d。

我尝试过以下方法：

x = data            %position in time domain
np = 2;
nz = 0;
sys = tfest(x,np,nz);

sys.Denominator

我收到以下错误

使用 tfest 时出错
“tfest”命令不能用于时间序列数据。考虑使用“ar”、“arx”或“armax”命令。

我知道我不能在 s 域中使用时域数据，但我不知道如何解决这个问题。我了解到armax不支持连续时间模型估计。我该如何解决我的问题？

Answer 1

在调用 tfest 之前，将用于测量输出的输入添加到数据对象中。该功能需要至少一个输入和一个输出数据序列才能运行。 如果输入不可用，我建议要么使用输入进行新的模拟，要么将时域数据转换为频域数据（假设脉冲作为输入）。