John Michell answered . 2024-05-15 19:57:47

• If data is transient (that is, the input is a step or impulse, or a pulse) and is not very long:

    data = iddata(outputsignal, inputsignal, sampletime); %sampletime = 1/Fs
    model = tfest(data, np); % np: number of poles
    w = linspace(0,pi/sampletime,100); % an example frequency grid (you can pick something different)
    G = freqresp(model,w); % or, simply do "bode(model)"

• If data has periodic inputs:

    data = iddata(outputsignal, inputsignal, sampletime)
    data2 = fft(data) % data in frequency domain
    model = tfest(data2, np); % np: number of poles
    w = linspace(0,pi/sampletime,100);
    G = freqresp(model,w); % or, simply do "bode(model)"

If you have MIMO data, you could replace TFEST with SSEST (state-space estimator) and follow the steps similar to the above.

Not satisfied with the answer ?? ASK NOW