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

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

Какая- то симметрия лишняя - если "дальше симметрично"
ФУРЬЕ.JPG

Андрей33
Сообщения: 12
Зарегистрирован: 01 окт 2015, 16:29

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

Сообщение Андрей33 »

Прогнал на STM-е, все нормально считает, один пик на 16 отсчете, но там библиотека другая.
Либо я функции неправильно входные данные даю, либо ошибка в библиотеке.
Пробовал функцию rfft - такой же результат.

Есть еще предположение:
fft в этой библиотеке использует таблицу синусов, область памяти которой в процессоре нужно вручную прописать.
Если таблица уместилась не вся, такой эффект может наблюдаться?

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

Так, надо немного остановиться и подумать...
У нас ведь выходной вектор тоже комплексный - и я надеюсь, что имеет такую структуру, как и входной, т.е.
1-2 действительная часть частоты 0;мнимая часть частоты 0;
3-4 действительная часть частоты 125 Гц;мнимая часть частоты 125 Гц;
5-6 действительная часть частоты 250 Гц;мнимая часть частоты 250 Гц;
7-8 действительная часть частоты 375 Гц;мнимая часть частоты 375 Гц;
.....
17-18 действительная часть частоты 1000 Гц;мнимая часть частоты 1000 Гц;
.....
64-65 действительная часть частоты 3850 Гц;мнимая часть частоты 3850 Гц;
65-66 действительная часть частоты 4000 Гц;мнимая часть частоты 4000 Гц; Частота Найквиста.
67-68 действительная часть частоты -125 Гц;мнимая часть частоты -125 Гц;
69-70 действительная часть частоты -250 Гц;мнимая часть частоты -250 Гц;
71-72 действительная часть частоты -375 Гц;мнимая часть частоты -375 Гц;
....
81-82 действительная часть частоты -1000 Гц;мнимая часть частоты -1000 Гц;
....
127-128 действительная часть частоты -3850 Гц;мнимая часть частоты -3850 Гц;

То есть 2 "пика" должно быть на 1000 Гц и на -1000 Гц

Андрей33
Сообщения: 12
Зарегистрирован: 01 окт 2015, 16:29

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

Сообщение Андрей33 »

Разобрался.
Таблицу синусов я размещал во флеш памяти прцессора, а с ней он походу не успевал просчитывать результат.
Переместил эту секцию вручную в ОЗУ и все заработало.
Спасибо за помощь.

Николай72
Сообщения: 5
Зарегистрирован: 01 дек 2015, 11:06

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

Сообщение Николай72 »

Добрый день!
ЦОС и Matlab пока только изучаю.
Попробовал запустить программу.
Раскомментированы строки кроме % fftshift(512);
Изменено значение 256 на 1024 (по замечанию Santik)
Результат работы скрипта прилагается

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

Николай72 писал(а): Результат работы скрипта прилагается
Где? :)

Николай72
Сообщения: 5
Зарегистрирован: 01 дек 2015, 11:06

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

Сообщение Николай72 »

Прошу прощения. Не разобрался. PDF-файл не прикрепляется.
Вложения
nuvsignal1.jpg

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

А в чем вопрос? Если частота Найквиста посередине - то всё нормально!

Николай72
Сообщения: 5
Зарегистрирован: 01 дек 2015, 11:06

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

Сообщение Николай72 »

То есть картина такая, какой ей положено быть в примере автора темы?

Аватара пользователя
Santik
Сообщения: 609
Зарегистрирован: 28 дек 2010, 08:04
Откуда: Мирный (Якутия)
Контактная информация:

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

Сообщение Santik »

Я не понял вопроса. У автора темы были неправильные данные. Или эти данные именно из этого примера?

Ответить