How can I implement Frequency shift Filter or LPTV filter using matalb

Question

Suppose i have signal r(n) = 1+sin(2pin/20), where i have three cyclic frequencies shift a = {a0,a1,a2} = [-1,0,1]2pi/20 such that my signal x(n) = r(n)exp(j