Este estudio presenta modelos analíticos para simular las propiedades térmicas de sistemas lineales triatómicos, utilizando el oscilador de Rosen-Morse modificado y el potencial del oscilador armónico para representar los modos vibracionales. Los modelos emplean funciones de partición existentes para derivar las funciones termodinámicas para los modos de flexión simétricos, asimétricos y degenerados de 2 veces. Estas funciones termodinámicas se aplican a moléculas triatómicas gaseosas como BO, HCN, N y SiN. Los resultados demuestran una alta precisión, con desviaciones absolutas porcentuales medias (MPAD) de menos del 0.17% para la entropía molar y la energía libre de Gibbs. Para la entalpía y la capacidad calorífica, los valores de MPAD están por debajo del 2% en comparación con los datos del Instituto Nacional de Estándares y Tecnología (NIST). Los hallazgos están en fuerte acuerdo con la literatura existente sobre moléculas triatómicas gaseosas, confirmando la fiabilidad de los modelos propuestos.