Many-body localization in large systems: Matrix-product-state approach