Как из fft собрать сигнал в некоторой точке (не по-быстрому а в лоб)?
Период гармоник равен T/1 T/2 T/3, начальная фаза равна арктангенсу, амплитуда - магнитуде комплексного числа? Не получается вычислить даже в нуле. Что я делаю не так (octave)?

T=32; f=(0:T-1)-(T-1)/2; fi=fft(f); k=0; g=0; for j=2:T/2-1, t0=atan2(imag(fi(j)),real(fi(j))); t=T/(j-1); g=g+abs(fi(j))*cos(2*pi*(k*t+t0)); end; g
figure; plot(f); figure; stem(real(fi),'r'); hold on; stem(imag(fi),'b')