% The modulation function:
function y=mod(data)
fc = 2000;
Ac = 1;
t = 0:1/4:(length(data)-1)/4;

%Generate carrier wave
carrier = Ac*sin(2*pi.*t);
%Generate BPSK signal using the generated binary data
y=data.*carrier + .01*randn(1,length(data));