We consider a stochastic sequence $ξ(m)$ with periodically stationary generalized multiple increments of fractional order which combines cyclostationary, multi-seasonal, integrated and fractionally integrated patterns. The filtering problem is solved for this type of sequences based on observations with a periodically stationary noise. When spectral densities are known and allow the canonical factorizations, we derive the mean square error and the spectral characteristics of the optimal estimate of the functional $Aξ=\sum_{k=0}^{\infty}{a}(k) ξ(-k)$. Formulas that determine the least favourable spectral densities and the minimax (robust) spectral characteristics of the optimal linear estimate of the functional are proposed in the case where the spectral densities are not known, but some sets of admissible spectral densities are given.