FFT, Интерпретация результата преобразования

framer
Сообщения: 5
Зарегистрирован: 23 фев 2021, 13:14

Re: FFT, Интерпретация результата преобразования

Сообщение framer »

Alexey-Rus. писал(а):
27 фев 2021, 16:41
Сергей, не могли бы Вы объяснить эту часть выражения: "sum(s0."? sum, как я понимаю, это сумма выражения в скобках при t от 0 до N. А как трактовать "s0."? Еще мне непонятно, что такое "-2i" и "angle(S1/S2)".
Спасибо.
Позволил себе поправить код под SciLab . Можно проверить и поэкспериментировать online https://cloud.scilab.in/ .

Код: Выделить всё

N = 128
#t = [0:0.1:N-1];
t = (0:N-1);

k = 11;
A = 1;
s0 = A * cos(2*%pi*k*t/N);

B = 0.7;
phi = 1.0;

dt = 0.0;
s1 = B * cos(2*%pi*k*(t - dt)/N + phi);

s0c = s0 .* exp(-%i * 2 * %pi * k * t / N);
s1c = s1 .* exp(-%i * 2 * %pi * k * t / N);

S0 = sum( s0c);
S1 = sum( s1c);

result = S1/S0;

gain = abs(result)
phase = atan(imag(result) / real(result))  


.* - Перемножение таблиц
exp(-%i * 2 * %pi * k * t / N); - Комплексная синусоида

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1053
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: FFT, Интерпретация результата преобразования

Сообщение Бахурин Сергей »

2i это 2*i, i - комплексная единица.
s0.* это поэлементное умножение вектора s0 на вектор комплексной экспоненты. По сути это один отсчёт fft соответствующий частоте k.
angle(S1/S0) функция расчёта угла (фазового сдвига) между комплексными спектальными отсчётами исходного и принятого сигналов.

Alexey-Rus.
Сообщения: 4
Зарегистрирован: 26 фев 2021, 08:49

Re: FFT, Интерпретация результата преобразования

Сообщение Alexey-Rus. »

Спасибо! Буду разбираться.

Ответить