Se investiga la estabilidad de los sistemas bilineales utilizando técnicas espectrales como el análisis modal selectivo. Se aplican modelos predictivos de sistemas bilineales basados en conocimiento inductivo extraído por técnicas de minería de datos grandes con búsqueda asociativa de patrones estadísticos. Se desarrolla un método y un algoritmo para la solución por elementos de la ecuación de Lyapunov matricial generalizada para sistemas bilineales discretos. El método se basa en calcular la secuencia de valores de un elemento fijo de la matriz solución, que depende del producto de los autovalores de la matriz de dinámica de la parte lineal y los elementos de las matrices de no linealidad. Se obtiene una condición suficiente para la convergencia de todas las secuencias, que también es una condición de estabilidad de sistemas BIBO (entrada acotada salida acotada) para el sistema bilineal.