n = (0:199).';
    s = 4.7*cos(0.12*pi*n-1.3) + ...
        3.8*cos(0.19*pi*n+0.8) + ...
        5.1*cos(0.23*pi*n+2.4) + ...
        2.0*randn(size(n));
    plot(n,s)