The paper presents a method to construct gn improved Lyapunov function for power system stability analysis, based on the generalized multivariable Popov criterion developed by Moore and Anderson. The method has been evaluated in two different power systems. Calculation results show that compared with the traditional Lure Lyapunov functions, the proposed Lyapunov function can detect a stability region much closer to the true power system stability boundary.