We propose a new general method to study critical points (CP) using the finite-size scaling of Lee-Yang zeros (LYZ). We first study the LYZ in the three-dimensional Ising model on finite lattices. We show that the ratios of multiple LYZ (Lee-Yang-zero ratios: LYZR) have useful scaling properties similar to the Binder cumulants, providing us with a novel method to study CP. In numerical simulations of the Ising model, we confirm that this method works well. We then apply the method to analyze the CP in the three-dimensional three-state Potts model and finite-temperature QCD in heavy-quark region, which are believed to belong to the same universality class as the Ising model. In these models, the partition function at complex parameters can be evaluated by the reweighting method, which allows us to determine the LYZ by varying coupling parameters continuously around the CP. We demonstrate that the LYZR method is powerful in determining the location of the CP in these models.