Neste trabalho, a propagação de sólitons em fibras óticas é analisada através de simulação numérica da equação não linear de Schrödinger, a qual descreve a propagação de pulsos óticos em fibras monomodo do tipo degrau. Uma vez que soluções analíticas para esta equação só podem ser obtidas em alguns casos específicos, implementaram-se dois métodos numéricos, possibilitando a análise da evolução de diferentes formas de pulsos incidentes ao longo de fibras com propriedades diversas de atenuação, dispersão e não linearidades. O método da propagação de Raios, cujo desempenho mostrou-se superior ao do método da série de Fourier, foi o escolhido para a obtenção dos resultados aqui apresentados. Várias características do sóliton fundamental, dos sólitons de ordens superiores, dos sólitons escuros e do fenômeno da interação entre pulsos adjacentes são apresentadas e discutidas, levando-se em consideração as possíveis implicações no desempenho de sistemas óticos. Contrabalançando os efeitos da dispersão da fibra com os efeitos não lineares da automodulações de fase, o que permite sua propagação sem alteração de forma, os sólitons encontram potencial aplicação na transmissão de altas taxas a longas distâncias.