We herein use a (vectorial) phase-field model description of the evolution of ferroelectric domains to be coupled with the equations of electroelasticity, with the aim of building a simulation framework for electromechanically active materials. The governing equations of the coupled model are discretized in strong form by means of isogeometric collocation and numerically solved by a staggered explicit approach. This is the first time that isogeometric collocation is used for such a complex problem, comprising a vectorial phase-field model for polarization coupled with the equations describing the electrical and the mechanical response of the system. Several numerical experiments are carried out to test the behavior of the adopted simulation framework. The obtained results are excellent and propose isogeometric collocation as an inexpensive but very accurate alternative to standard finite element discretizations also for complex coupled problems.