FourierSeries[(1/2 (SquareWave[t/(4pi)] + 1)Sin[t])]