Este trabalho trata de funções de Lyapunov lineares por partes para análise de estabilidade de sistemas lineares discretos no tempo com entradas saturáveis. São apresentadas novas condições necessárias e suficientes para funções lineares por partes positivas definidas serem funções de Lyapunov. Um procedimento computacional é proposto para determinação dessas funções de Lyapunov e suas correspondentes regiões de estabilidade assintótica. Comparadas às funções de Minkowski, as funções lineares por partes apresentam desempenho estritamente superior e são naturalmente mais flexíveis e adaptadas ao comportamento dinâmico radialmente variável dos sistemas saturados.