A new model-reference data-driven approach is presented which uses the frequency response data of a system in order to avoid the problem of unmodeled dynamics associated with low-order parametric models. It is shown that a convex optimization problem can be formulated (in either the H-infinity , H-2 or H-1 sense) to shape the closed-loop sensitivity functions while guaranteeing the closed-loop stability. The effectiveness of the method is illustrated by considering several case studies where the proposed design scheme is applied in both simulation and to a power converter control system for a specific accelerator requirement at CERN.