Phasor Measurement Units (PMUs) are the backbone of smart grids that are able to measure power system observability in real-time. The deployment of synchronized sensors in power networks opens up the advantage of real-time monitoring of the network state. An optimal number of PMUs must be installed to ensure system observability. For that reason, an objective function is minimized, reflecting the cost of PMU installation around the power grid. As a result, a minimization model is declared where the objective function is defined over an adequate number of constraints on a binary decision variable domain. To achieve maximum network observability, there is a need to find the best number of PMUs and put them in appropriate locations around the power grid. Hence, maximization models are declared in a decision-making way to obtain optimality satisfying a guaranteed stopping and optimality criteria. The best performance metrics are achieved using binary integer, semi-definite, and binary polynomial models to encounter the optimal number of PMUs with suitable PMU positioning sites. All optimization models are implemented with powerful optimization solvers in MATLAB to obtain the global solution point.