En este artículo se analiza un método alternativo para calcular los parámetros modales en fibras ópticas, como las constantes de propagación, las distribuciones transversales y la anisotropía, debido a fenómenos lineales y no lineales que actúan como perturbaciones causadas por regiones de sílice dopadas. Este método se basa en una formulación hamiltoniana de la ecuación de Helmholtz y en la teoría de perturbaciones estacionarias, que permite una descripción vectorial completa de los componentes del campo eléctrico cuando se incluyen las inhomogeneidades lineales anisotrópicas y la no linealidad de Kerr. Se pueden hallar los parámetros lineales y no lineales para cada modo de propagación, y su precisión se ha probado con éxito al compararla con los cálculos numéricos del método de elementos finitos vectoriales, y los resultados se publican en la literatura. Este método facilita el cálculo de los efectos de perturbación distribuidos espacialmente en los componentes individuales del campo eléctrico para cada modo de propagación.