Este trabalho consistiu da realização da análise de regressão periódica dos totais mensais de precipitação de oito localidades do Estado do Paraná. Foram derivados modelos matemáticos que descrevem o transcurso desta variável, considerando um período de doze meses. A análise foi desenvolvida através da álgebra matricial e realizada a expansão em série de Fourier de uma função f(t), periódica e definida para o intervalo. Os dados originais foram ajustados para meses de 30 dias, transformados para normalizar a distribuição. Posteriormente foram determinados todos os componentes harmônicos das funções referentes a cada um dos oito locais. Após a verificação da normalidade dos componentes harmônicos, realizou-se a análise de variância, o que permitiu identificar os componentes harmônicos significativos, que foram escolhidos para fazer parte da equação de regressão correspondente a cada local. Estas equações explicaram de 89,56% a 99,60% da variação devida a meses, considerando o conjunto das localidades estudadas. A probabilidade de ocorrência de um ano conforme o modelo, expressa em porcentagem, variou de 14,19% a 68,42%, enquanto a probabilidade de ocorrência de um ano conforme a média variou de 0,02% a 1,87%.