Помощь - Поиск - Пользователи - Календарь
Полная версия этой страницы: Непонятки с "ручным" обратным преобразованием Фурье?...
Форум разработчиков электроники ELECTRONIX.ru > Cистемный уровень проектирования > Математика и Физика
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?, где в коде ошибка?
Leshii
Вопрс снят, ошибка найдена.
Для просмотра полной версии этой страницы, пожалуйста, пройдите по ссылке.
Invision Power Board © 2001-2025 Invision Power Services, Inc.