A method introducing fictious variables for resolving the linear model of reactive power and voltages control in power systems is presented. The method theoritically guarantees the resolvability and optimizition of linear programming model. This method has been proved successful by IEEE 6 -node test system and the several practical power systems such as 49-node, 155 -node systems.