frequenza = 0:0.01:20;
t = 0:0.01:10;
ampiezza = 2 * abs(sinc(frequenza));
fase = frequenza*pi*(-1) - pi*(sinc(frequenza)<0);

rect = ampiezza * cos(frequenza'*t + fase'*ones(1,1001));
plot(t, rect);