We investigate the asymptotic behavior of a family of multiple orthogonal polynomials that is naturally linked with the normal matrix model with a monomial potential of arbitrary degree d + 1. The polynomials that we investigate are multiple orthogonal with respect to a system of d analytic weights defined on a symmetric (d + 1)-star centered at the origin. In the first part we analyze in detail a vector equilibrium problem involving a system of d interacting measures (µ1, . . . , µ d ) supported on star-like sets in the plane. We show that in the subcritical regime, the first component µ * 1 of the solution to this problem is the asymptotic zero distribution of the multiple orthogonal polynomials. It also characterizes the domain where the eigenvalues in the normal matrix model accumulate, in the sense that the Schwarz function associated with the boundary of this domain can be expressed explicitly in terms of µ * 1 . The second part of the paper is devoted to the asymptotic analysis of the multiple orthogonal polynomials. The asymptotic results are obtained again in the subcritical regime, and they follow from the Deift/Zhou steepest descent analysis of a Riemann-Hilbert problem of size (d + 1) × (d + 1). The vector equilibrium problem and the Riemann-Hilbert problem that we investigate are generalizations of those studied recently by Bleher-Kuijlaars in the case d = 2.