f=0.0000001:0.001:1;
Tb=1;
x=2*pi*f*Tb;
 
plot(f,sp);
%pause;
hold on
sinc=sin(x/2)./(x/2);
sinc_patrat=sinc.*sinc;
spec=sinc_patrat.*sp;
plot(f,spec);