matlab - 如何找到以不同速率采样的输入和输出之间的传递函数？

Question

我知道通常我会ftest()在使用iddata(). 但是，为了iddata()正常工作，我需要以相同的速率对输入和输出数据进行采样。是否有与速率无关的变体iddata()或任何其他方式可以让我完成我需要完成的事情？

Answer 1

我正在解决同样的问题（https://dsp.stackexchange.com/questions/19458/how-to-compute-transfer-function-from-experimental-data），但我不确定我是否找到了解决方法，但我与您分享我的想法，因此我们可能会找到适合我们双方的解决方案（如果您已经找到了方法，请分享）。

方法一

如果您在时域中有信号，则可以同步它们，然后使用该tfestimate功能。

% Define timeseries 
ts_output = timeseries(x,time1,'Name','output');
ts_input = timeseries(y,time2,'Name','input');

% Synchronization
[ts_output,ts_input] = synchronize(ts_output,ts_input,'uniform',...
'interval',delta_t);

% Compute transfer function
Fs = 1/delta_t;
[Txy,W] = tfestimate(ts_input.data,ts_output.data,[],[],[],Fs);

方法二

您可以重新采样以较低频率采样的信号，而不是同步：假设Fs1 > Fs2

[P,Q] = rat(Fs1/Fs2);
y2 = resample(y2,P,Q);