clear?all;
close?all;

f0=5;

fi0=pi/2;
fi1=pi/4;

NumCos=16;
TPulse=NumCos*（1/f0）;

NxBinary=8;
xBinary=randint（NxBinary12）;

ts=0.001;
t=ts:ts:TPulse*NxBinary;


for?k=1:NxBinary
????B（（k-1）*length（ts:ts:TPulse）+1:k*length（ts:ts:TPulse））=xBinary（k）;
end

A=1;

y=（1+B.*cos（2*pi*f0*t+fi0））.*cos（2*pi*2*f0*t）;


figure（1）;
plot（ty）;

y_fft=fft（y）;

figure（2）;
plot（abs（y_fft））;


y_fft（1:257）=0;
y_fft（25345:length（y_fft））=0;

figure（3）;
plot（abs（y_fft））;

y_ifft=ifft（y_fft）;
f