calculate matrix to predict control variables depending on past and prospective changes of measurable disturbance variables using step a response disturbance model
function DisturbanceMatrix input Real[:, :, :] Sd "step response matrix of measurable disturbances"; input Integer Nl "lower prediction horizon"; input Integer Np "upper prediction horizon"; input Integer Nmd "model horizon of measurable disturbances"; input Integer p "number of control systems outputs"; input Integer q "number of measurable disturbance variables"; output Real[p*(Np - Nl + 1), q*(Nmd - 1 + Np)] E; end DisturbanceMatrix;