yading@10: function Y = genspecsines_sin(ploc, pmag, pphase, N)
yading@10: y=zeros(N,1);
yading@10: w=blackmanharris(N);
yading@10: w=w./sum(w);
yading@10: w(2:N)=w(2:N)*2;
yading@10: for k=1:length(ploc)
yading@10: f=ploc(k)/N;
yading@10: A=10ˆ(pmag(k)/20);
yading@10: phi=pphase(k);
yading@10: n=[0:1:N-1]';
yading@10: y=A*cos(2*pi*f.*n+phi)+y;
yading@10: end
yading@10: y=y.*w;
yading@10: Y=abs(fft(y));