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?, где в коде ошибка?