clear; clf;
n=[0:1:1000]; n2=[0:1:1000];
x=cos(0.48*pi*n)+cos(0.52*pi*n);
psd(x,1024,2,kaiser(128,pi),0,0.95);
title('using function psd'); xlabel('frequency in pi units');