A new approach to experimental design for identification of closed-loop models is presented. The method considers the design of an experiment by minimizing experimental cost, subject to probabilistic bounds on the input and output signals, and quality constraints on the identified model. The input and output bounds are common in many industrial processes due to physical limitations of actuators. The aforementioned constraints make the problem non-convex. By assuming that the experiment is a realization of a stationary process with finite memory and finite alphabet, we use results from graph-theory to relax the problem. The key feature of this approach is that the problem becomes convex even for non-linear feedback systems. A numerical example shows that the proposed technique is an attractive alternative for closed-loop system identification.