Есть гармоника частотой 62 Гц есть выборка 48000 сэмпла ДПФ на 48000 и FFT на 512 и частота дискретизации 48 кГц.

Суть вопроса вот вчём мне нужно из сигнала вычесть гармонику с частотой 62 Гц и зделать FFT на 512 в которой этой гармоники недолжно быть. Я вот что подумал если мы зделаем так: зделаем ДПФ 48000 на гармонику 62 Гц и определим её фазу и амплитуду после чего сгенерируем эту гармонику и из имеющихся 48000 сэмплов вычтем эту гармонику после чего зделаем FFT на 512 будет ли эта операция коректной? Или есть другое решения? Мне нужно максимально точное.

Затрата ресурсов не важна время тоже.

Сообщение отредактировал ivan219 - Sep 17 2010, 19:28

..может я чего на понял ,
но по моему - надо сделать fft 512 точек , потом обнулить ненужную гармонику в результате , сделать обратное fft 512 точек , получим сигнал без вашей гармоники .

з.ы.

всегда думал что ДПФ и fft - одно и тоже, в смысле конечного результата.

Нет если у меня была бы кратная гармоника т.е. 48000 / 512 = 93,75 Гц то да метот подашол. Но гармоника не кратная а значит её спектр растекётся по всем 256 бинам.
И ДПФ это почти БПФ (FFT) гораздо медленне но позволяет вычислить любую гармонику кратную целому числу а не только двойке как в БПФ.

Сообщение отредактировал ivan219 - Sep 17 2010, 20:25

Дааааааааа.... О какой специфике вопроса мы можем здесь говорить, если человек даже не знает, что такое знаки препинания. Про грамматику речи даже не идет. Ужас какой-то!!!!!!!! Не сочтите за грубость.

Сообщение отредактировал Serg76 - Sep 17 2010, 20:33

КаДДа нечего сказать по-теме - лучше промолчать...

Наверное, да

Вычитать сгенерированную гармонику можно, если у вас 62Гц достаточно точно и стабильно в записанном окне. Но при этом нужно:

1. нужно предварительно, перед определением параметров 62 Гц вычесть из сигнала среднее;
2. воспользоваться подходящей оконной функцией;
3. Определить вклад 62 Гц точно алгоритмом Герцеля, а не при помощи FFT, так как 48 кГц не кратно 62 Гц.

Но в целом успех этого предприятия зависит от неописанных вами и поэтому неизвестных мне свойств сигнала и помехи.

А нельзя ли не обращать внимания на эту гармонику?
Если нельзя, я бы просто сделал FIR максимально возможной длины.

62Гц стабильно во всём окне. Среднее значение не подумал новерное правильно. А вот оконная функция зачем онаже исказить сигнал? А заместо алгоритма Герцеля ДПФ на 48000 выборок не подайдёт ведь ресурсы не важны??? И FFT я буду делать уже после того как 62 Гц будут убраны из сигнала.


Думал но он вырежит определённую полосу. А надо только одну гармонику 62 Гц.

Сообщение отредактировал ivan219 - Sep 18 2010, 07:44

А она есть? Гармоника? Это абстракция, условность... Под гармоникой FFT всегда подразумевают ту или иную полосу... Бин DFT/FFT ничто иное как полосовой фильтр с sinc-образной АЧХ... ну если без окон...

Кстати 2-а соседних к 62 Hz бина в принципе можно так исказить, что они будут давить именно 62 Hz. Делаем режекторный фильтр на 62 Hz, имеем его передаточную характеристику... Делаем ее FFT... Смотрим получилось ли - из-за влияния окон может быть хреновой реализацией...

P.S.: а уж если совсем абстрагироваться - то под гармоникой понимают ту или иную корреляцию с той или иной функцией ...

Уидроу Б.,Стринз С.
Адаптивная обработка сигналов.
http://lord-n.narod.ru/walla.html#uidrouAOC

В этой книге рассказано как просто такую решить.

Потому что 48000 на 62 не делится. Оконная функция нужна только для определения параметров гармоники.

Ну как это её може не быть я её же сам и синтезирую. Вот блин как всё сложно голова кругом



Ага нашол адаптивный режекторный фильтр надо почитать, спасибо.


Ну да зато на промежутке в 48000 семпла и частотой дискретизации 48 кГц укладывается целое количество колебаний с частотой 62 Гц по этому я и взял ДПФ 48000 окно длинной в секунду.

Сообщение отредактировал ivan219 - Sep 18 2010, 09:55

А вот так ...

При выборке ограниченной длины точность измерения (вырезания) частоты не будет лучше величины, обратной времени выборки. Максимальная точность (в центре выборки) будет при размере ядра FIR равной времени выборки.