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