The intersample behavior of the input signals influences the estimation, simulation and prediction of continuous-time models. A sampled signal is characterized only by its values at the sampling instants. However, when you apply a continuous-time input to a continuous-time system, the output values at the sampling instants depend on the inputs at the sampling instants and on the inputs between these points.
idfrd objects have an
which stores how the input behaves between the sampling instants.
You can specify the behavior between the samples to be piecewise constant
(zero-order hold), linearly interpolated between the samples (first-order
hold) or band-limited. A band-limited intersample behavior of the
input signal means:
A filtered input signal (an input of finite bandwidth) was used to excite the system dynamics
The input was measured using a sampling device (A/D converter with antialiasing) that reported it to be band-limited even though the true input entering the system was piecewise constant or linear. In this case, the sampling devices can be assumed to be a part of the system being modeled.
When input signal is band-limited, the estimation is performed as follows:
Time-domain data is converted into frequency domain
fft and the sample
time of the data is set to zero.
Discrete-time frequency domain data (
= 'frequency' or
idfrd with sample time Ts≠0)
is treated as continuous-time data by setting the sample time Ts to
The resulting continuous-time frequency domain data is used for model estimation. For more information, see Pintelon, R. and J. Schoukens, System Identification. A Frequency Domain Approach, section 10.2, pp-352-356,Wiley-IEEE Press, New York, 2001.
Similarly, the intersample behavior of the input data affects
the results of simulation and prediction of continuous-time models.
InterSample property to choose the right
algorithm for computing model response.
The following example simulates a system using first-order hold (
foh ) intersample behavior for input signal.
sys = idtf([-1 -2],[1 2 1 0.5]); rng('default') u = idinput([100 1 5],'sine',,,[5 10 1]); Ts = 2; y = lsim(sys, u, (0:Ts:999)', 'foh');
iddata object for the simulated input-output data.
data = iddata(y,u,Ts);
The default intersample behavior is zero-order hold (
ans = zoh
Estimate a transfer function using this data.
np = 3; % number of poles nz = 1; % number of zeros opt = tfestOptions('InitMethod','all','Display','on'); opt.SearchOption.MaxIter = 100; modelZOH = tfest(data,np,nz,opt)
modelZOH = From input "u1" to output "y1": -217.2 s - 391.6 --------------------------------- s^3 + 354.4 s^2 + 140.2 s + 112.4 Continuous-time identified transfer function. Parameterization: Number of poles: 3 Number of zeros: 1 Number of free coefficients: 5 Use "tfdata", "getpvec", "getcov" for parameters and their uncertainties. Status: Estimated using TFEST on time domain data "data". Fit to estimation data: 81.38% (simulation focus) FPE: 0.1261, MSE: 0.111
The model gives about 80% fit to data. The sample time of the data is large enough that intersample inaccuracy (using
zoh rather than
foh ) leads to significant modeling errors.
Re-estimate the model using
foh intersample behavior.
data.InterSample = 'foh'; modelFOH = tfest(data, np, nz,opt)
modelFOH = From input "u1" to output "y1": -1.197 s - 0.06843 ------------------------------------- s^3 + 0.4824 s^2 + 0.3258 s + 0.01723 Continuous-time identified transfer function. Parameterization: Number of poles: 3 Number of zeros: 1 Number of free coefficients: 5 Use "tfdata", "getpvec", "getcov" for parameters and their uncertainties. Status: Estimated using TFEST on time domain data "data". Fit to estimation data: 97.7% (simulation focus) FPE: 0.001747, MSE: 0.001693
modelFOH is able to retrieve the original system correctly.
Compare the model outputs with data.
compare(data, modelZOH, modelFOH)
modelZOH is compared to data whose intersample behavior is
foh. Therefore, its fit decreases to around 70%.