Auxiliary model-based interval-varying multi-innovation least squares identification for multivariable OE-like systems with scarce measurements