% repliment de spectre
%---------------------
FS = 11025;
h = 1/FS;
t = 4;
th = 0:h:t;

x = .98*cos(2*pi*(1000*th+2000*th.^2));
sound(x,FS);
specgram(x,[],FS);
