samer@32: % fq2poles - resonant frequencies to poles
samer@32: % fq2poles :: nonneg, nonneg -> [[2]->complex].
samer@32: % fq2poles :: [[1,N]->nonneg], [[1,N]->nonneg] -> [[2,N]->complex].
samer@32: function p=fq2poles(f,q)
samer@32: 
samer@32: 	p1=exp(2*pi*(i-q).*f);
samer@32: 	p=[p1;conj(p1)];