на главнуюВсе эхи RU.ALGORITHMS
войти ?

выделение основных частот сигнала

От Nickita A Startcev (2:5022/49.61) к Alex Mizrahi

В ответ на Заголовок предыдущего сообщения в треде (Имя Автора)


Привет, Alex !


14 Jun 10 , 13:32 Alex Mizrahi писал к Nickita A Startcev:

NAS>> Можно сделать так:
NAS>> сделать БПФ, посчитать модули амплитуд (сумму квадратов
NAS>> проекций) выделить самый высоки пик. Сворачивать исходный сигнал
NAS>> с синусом и косинусом этой частоты, а потом чуток шевелить ее
NAS>> (частоту),
NAS>> ориентируясь на максимизацию амплитуды.

AM> Амплитуду?

Интенсивность суммы квадратов свёрток (то что у тебя ниже обозначено как A^2+B^2).

AM> Я считал интеграл свёрток

AM> A := int(f(x)*sin(f*x), x=-Pi..Pi);
AM> B := int(f(x)*cos(f*x), x=-Pi..Pi);
AM> А потом максимизировал A^2 + B^2.

не понял. f(x) и f*x - это разные f?

AM> Максимум оказался совсем не на нужной частоте.

AM> Более того, если у вектора съехавшая фаза, то максимум вообще
AM> находится в нуле!

AM> Я подозреваю проблема вот в чём -- т.к. у меня период функции больше
AM> чем исследуемый интервал,

Ага. Проблема именно в этом. Hа исследуемом интервале должно помещаться несколько периодов.


NAS>> ps: сам я плотник, если где сильно прогнал - жду комментариев от
NAS>>
AM> спецов.

AM> Что-то мне кажется тут практически никого уже не осталось :(
AM> А жаль...

Посмотри в ЖЖ. Там Санитар Женя часто бывает.

. С уважением, Hикита.
icq:240059686, lj-user:nicka_startcev
... Добросовестно и усердно веселиться?!

--- GoldED+/LNX 1.1.5-
* Origin: Люди Билли не любили... (c) (2:5022/49.61)

Ответы на это письмо:

From: Username
Заголовок следующего сообщения в треде может быть длинным и его придется перенести на новую строку

From: Username
Или коротким

FGHI-url этого письма: area://RU.ALGORITHMS?msgid=2:5022/49.61+4c1674bc