A Screening Rule-Based Iterative Numerical Method for Observability Analysis