Variational study of isotropic spin‐1 chains using matrix‐product states