Непонятки с "ручным" обратным преобразованием Фурье?... Опции

[Leshii]

clear all
t=0:(1024-1);
f=0.4;
A=2.5;
fi=pi/4;
s0=A*exp(sqrt(-1)*2*pi*f*t+fi);
s=[zeros(1, 100) s0 zeros(1,250)];
s=s+(0.5*randn(1,length(s))+sqrt(-1)*0.5*randn(1,length(s)));
t=0:length(s)-1;

mfT=abs(ifft(fft(s)));

sp_mfT=fft(s)/length(s);

ug=unwrap(angle(sp_mfT));
ab=abs(sp_mfT);

sall=zeros(1,length(t));

for N=1:length(sp_mfT)
A_1=ab(N);
fi_1=ug(N);
f_1=(1/length(sp_mfT))*(N-1);
s_1=A_1*exp(sqrt(-1)*2*pi*f_1*t-fi_1);
sall=sall+s_1;
end;

figure, plot(real(mfT),'r'), hold on, plot(abs(sall),'*'), vaxis



Почему не совпадают mfT и sall?, где в коде ошибка?