State estimation (SE) is an algorithmic function of an energy management system (EMS). SE provides an actual-time monitoring and control of modern electrical power grids. State Estimation can be worked with sufficiency using Phasor Measurement Units optimally placed within a power grid. This paper concerns the implementation of proper algorithms embedded in optimization solvers to the optimal PMU localization problem solving globally. The optimization model is formulated as a 0 - 1 nonlinear minimization problem. The problem is transformed to a polyhedron using linearization methods and B&B tree. In this model, we use a linear cost function under polynomial constraints and binary restrictions on the design variables in a symbolic format. This mathematical model is programmed in the YALMIP environment which is fully compatible with MATLAB. The 0 - 1 Nonlinear Programming (NLP) model is suitable for getting concisely global optimal solutions. The optimal solution is given by a wrapped optimization engine including a local optimizer routine performing together with a mixed-Integer-Linear Programming routine. The solution is achieved within a zero-gap precisely encountered during the iterative process. This tolerance criterion is a necessity for a successful implementation of the B&B tree because it ensures global optimality with an acceptance relative gap. The minimization model is implemented in a YALMIP code fully compatible with MATLAB in two stages. Initially, an objective function with one term is minimized to discover a number of sensors for wide-area monitoring, control and state estimator applications. Then, an extra product is considered in the objective to suffice maximum reliability for observing the network buses. The numerical minimization models are applied to standard power networks in the direction to be solved globally.