In R2012a, use the ARMAX function with "IntegrateNoise" attribute specified using name-value pair. See ARMAX documentation for more info.
For example, model = armax(data, [na nc], 'IntegrateNoise', true) estimates an ARIMA model of time series "data".
