Fs = 48000; F = 1000; t = 0 : 1/Fs: 0.1; S0 = 1; phase = 0; Sre = 0.5 * S0 * cos(2*pi*F*t+phase); Sim = 0.5 * S0 * sin(2*pi*F*t+phase); S = complex(Sre, Sim); Sp = S0*exp(1i*2*pi*F*t+phase); Fc = 15e3; Sm = Sre.*cos(2*pi*Fc*t+phase);