x= [1 1 0 0];
Tb= 0.1;
t=[];
m=[];
fs=1000;
dt= 1/fs;
Ac=2;
fc=10;
ns=length(x);
for index1=1:ns
temp=(index1-1)*Tb:dt:index1-Tb;
if (x(index1)==1)
sbits=ones(1,length(temp));
else
sbits=zeros(1,length(temp));
end
t=[t,temp];
m=[m,sbits];
end
car=Ac*sin(2*3.14*fc*t);
yout=m.*car;
plot (t,ns);