Neste trabalho, um modelo de otimização é formulado para se estimar estados harmônicos a partir de um algoritmo robusto à introdução de erros grosseiros. Na abordagem proposta, magnitudes e ângulos de tensões harmônicas são as variáveis de estados a serem determinadas pela solução do problema de otimização proposto. Medições fasoriais sincronizadas são consideradas para a abordagem sendo o sistema IEEE de 14 barras utilizado para validar o método proposto e atestar a eficiência do mesmo considerando cargas não lineares conectadas ao sistema. Comparações são realizadas com um fluxo de potência e o método tradicional dos Mínimos Quadrados Ponderados (MQP). A fim de se melhorar o desempenho computacional do processo da estimativa de estados para cada ordem harmônica, o uso de processamento paralelo é adotado. Palavras-chave: estimação de estados; qualidade de energia; distorção harmônica; fluxo de potência.