реклама на сайте
подробности

 
 
> БПФ, "грязная" Im составляющая
AndreyVN
сообщение Mar 18 2015, 18:34
Сообщение #1


Знающий
****

Группа: Свой
Сообщений: 754
Регистрация: 29-06-06
Из: Volgograd
Пользователь №: 18 458



Всем привет!

Отлаживаю алгоритм комплексного БПФ, формально работает корректно, функцию туда-сюда гоняет без искажений, для синуса пики приблизительно на частоте N/(2*Pi), N-N/(2*Pi).
Не нравится как выглядит аргумент результата преобразования - с какими-то шумами (первый рисунок). Для сравнения второй рисунок - результат преобразования чистого синуса.

красная линия (пики слева и справа) - модуль комплексных к-ов Фурье образа;
зеленая линия - аргумент комплексных к-ов Фурье образа;
синяя линия - действительная часть обратного преобразования Фурье (мнимая часть =0).

Прикрепленное изображение
Прикрепленное изображение
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
ViKo
сообщение Mar 19 2015, 05:51
Сообщение #2


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Аргумент - это же и есть фаза?
Думаю, проблема вылазит именно из-за кусочно-линейной интерполяции, из-за резких изломов. Амплитуда этих частот мала, незаметна, а скачки фазы - явно видны.
Go to the top of the page
 
+Quote Post
AndreyVN
сообщение Mar 19 2015, 09:45
Сообщение #3


Знающий
****

Группа: Свой
Сообщений: 754
Регистрация: 29-06-06
Из: Volgograd
Пользователь №: 18 458



Цитата(ViKo @ Mar 19 2015, 09:51) *
Аргумент - это же и есть фаза?
Думаю, проблема вылазит именно из-за кусочно-линейной интерполяции, из-за резких изломов. Амплитуда этих частот мала, незаметна, а скачки фазы - явно видны.


Поэксперементировал еще немного, шумы фазы от гладкости действительно зависят, но почему-то сохраняются на гладкой функции тоже. Убрал кусочно-линейную интерполяцию и скормил БПФ Лоренцеву кривую F = 1/(1+X^2). Шумы сильно уменьшились, но совсем не исчезли.

Исходные данные сохранил в виде массива текстовых значений, разделенных запятыми, может кто-то найдет минутку, возьмет преобразование Фурье от этих данных? Все-таки есть опасения, что у меня какой-то, особенный, склонный к шумам алгоритм.

Прикрепленное изображение

Прикрепленный файл  Lorents.zip ( 2.71 килобайт ) Кол-во скачиваний: 77
Go to the top of the page
 
+Quote Post



Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 22nd July 2025 - 16:18
Рейтинг@Mail.ru


Страница сгенерированна за 0.01402 секунд с 7
ELECTRONIX ©2004-2016