Full-wave electromagnetic (EM) simulation tools have become ubiquitous in the design of microwave components. In some cases, e.g., miniaturized microstrip components, EM analysis is mandatory due to considerable cross-coupling effects that cannot be accounted for otherwise (e.g., by means of equivalent circuits). These effects are particularly pronounced in the structures involving slow-wave compact cells and their numerical optimization is challenging due to expensive simulations and large number of parameters. In this paper, a novel gradient-based procedure with numerical derivatives is proposed for expedited optimization of compact microstrip impedance matching transformers. The method restricts the use of finite differentiation which is replaced for selected parameters by a rank-one Broyden updating formula. The usage of the formula is governed by an acceptance parameter which is made dependent on the parameter space dimensionality. This facilitates handling circuits of various complexities. The proposed approach is validated using three impedance matching transformer circuits with the number of parameters varying from ten to twenty. A significant speedup of up to 50 percent is demonstrated with respect to the reference algorithm. . His research interests include CAD and modeling of microwave and antenna structures, simulation-driven design, surrogate-based optimization, space mapping, circuit theory, analog signal processing, evolutionary computation and numerical analysis. Anna PIETRENKO-DABROWSKA (corresponding author) received the M.Sc. and Ph.D. degrees in Electronic Engineering from Gdansk University of Technology, Poland, in 1998 and 2007, respectively. Currently, she is an Associate Professor with Gdansk University of Technology, Poland. Her research interests include simulationdriven design, design optimization, control theory, modeling of microwave and antenna structures, numerical analysis.