xx=linspace(0,2*pi,33);
x = zeros(1,32);
for i=1:32
    x(i)=xx(i);
end
y=(cos(14.*x)+2*sin(15.*x)+3.*cos(16.*x))/16.;
z=fft(y);
zz = fftshift(z);
z
zz