This is caused by roundoff error. When time vectors are defined using the above technique, depending on the value of the step size used, there will be times when the time steps will not match Simulink’s time-step. To take a closer look at the difference between the time step that was provided and the time steps that Simulink generates, use the hexadecimal representation through the following command:
format hex
You can then compare the time steps specified using the colon operator technique mentioned above, and the time steps that Simulink generated as follows:
>> t
t =
0000000000000000
3fb999999999999a
3fc999999999999a
3fd3333333333333
3fd999999999999a
3fe0000000000000
>> tout
tout =
0000000000000000
3fb999999999999a
3fc999999999999a
3fd3333333333334
3fd999999999999a
3fe0000000000000
Where “t” represents the user time steps and “tout” represents the time steps generated by Simulink. In particular, notice that the fourth element of these two variables is different by one bit. This causes the Inport block to use the Lagrange interpolation to interpolate the input data and match the Simulink’s time steps. In return, certain operations on the input data results in an unexpected behavior.
The best way to represent evenly spaced time steps such that the time steps will match Simulink’s, is to use the following technique:
time = stepSize * [0:N]'
Where “N” represents the number of time steps, similar to:
N = (stopTime - startTime) / stepSize;
Using the attached model and the following code:
dt = 1/10;
time = dt*(0:5)';
portvalues1 = [0 0 0 0 1 0]';
portvalues2 = [0 0 0 0 1 0]';
t = time;
u = [portvalues1, portvalues2];
sim('test', [t(1) t(end)], [],[t u]);
correct_output = u(:,1)./u(:,2)
simulink_output = yout
The correct result is achieved:
correct_output =
NaN
NaN
NaN
NaN
1
NaN
simulink_output =
NaN
NaN
NaN
NaN
1
NaN
Additional information has been updated in the documentation on this subject. You can find this information under "Importing Data to Test a Discrete Algorithm" at the documentation link here;
<http://www.mathworks.com/help/releases/R2011b/toolbox/simulink/ug/bsu_jie.html>
