Objetivo: La máxima extracción de potencia de los sistemas de cosecha de energía basados en el efecto piezoeléctrico es crucial para el desarrollo de diversas tecnologías actuales. Sin embargo, el uso directo del teorema de transferencia de máxima potencia en este tipo de captadores de energía produce modelos de carga inestables y no causales, que impiden su síntesis a través de un circuito electrónico de potencia y reducen la potencia disponible para el hardware específico de la aplicación. En este artículo, se presenta el cálculo de una función de transferencia causal y estable, que aproxima la máxima extracción de potencia de un cosechador piezoeléctrico lineal. Metodología: Dicho modelo se obtiene mediante el algoritmo de Optimización por Enjambre de Partículas (PSO), seleccionado por su sencillez de implementación. En este estudio, se utiliza una variación del algoritmo original para incluir las restricciones de estabilidad. Resultados: El algoritmo se empleó para obtener modelos de 7º, 8º, 16º y 20º grado para evaluar el impacto de la complejidad del modelo en los valores de la función objetivo obtenidos. El mínimo error cuadrático acumulado alcanzó un valor de la función objetivo de aproximadamente 6 · 10−7Ω −2 alcanzado por el modelo de 7º grado. Conclusiones: Se obtuvo un modelo estable y causal mediante una implementación de PSO. Los resultados sugieren una baja dependencia del valor mínimo de la función objetivo con el aumento de complejidad del modelo. No obstante, el mejor modelo coincidió con el grado de la función de transferencia de la admitancia. Finalmente, la implementación del algoritmo PSO sin restricciones produjo mejores resultados, sugiriendo que un algoritmo de optimización diferente puede alcanzar mejores resultados