Mercurial > hg > pmhd

annotate extra/genspecsines_sin.m @ 13:844d341cf643 tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 function Y = genspecsines_sin(ploc, pmag, pphase, N)
yading@10 2 y=zeros(N,1);
yading@10 3 w=blackmanharris(N);
yading@10 4 w=w./sum(w);
yading@10 5 w(2:N)=w(2:N)*2;
yading@10 6 for k=1:length(ploc)
yading@10 7 f=ploc(k)/N;
yading@10 8 A=10ˆ(pmag(k)/20);
yading@10 9 phi=pphase(k);
yading@10 10 n=[0:1:N-1]';
yading@10 11 y=A*cos(2*pi*f.*n+phi)+y;
yading@10 12 end
yading@10 13 y=y.*w;
yading@10 14 Y=abs(fft(y));