Google Groups no longer supports new Usenet posts or subscriptions. Historical content remains viewable.
Dismiss

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

6 views
Skip to first unread message

Michael Yuzhakov

unread,
Jan 27, 2012, 11:29:47 AM1/27/12
to
Есть кто живой?

Интересует книжка по сабжу. Желательно без особого мозгоклюйства,
применительно к практике.
Буду признателен за наводку.

С большим приветом, Южаков Михаил.
_________________________________________________________
е-мыло: yuzhakov<собака>etel.ru

Valentin Kuznetsov

unread,
Jan 28, 2012, 6:19:22 AM1/28/12
to
Мои пpибоpы сказали `Spam detected`, я пpовеpил, а оказалось, это Michael Yuzhakov пишет All:

MY> Интересует книжка по сабжу. Желательно без особого
MY> мозгоклюйства, применительно к практике.
Л.А.Залманзон. Преобразование Фурье, Уорша, Харра и их применение
Легко яндексится
Проще и понятнее не встречал

Valentin

Michael Yuzhakov

unread,
Mar 18, 2012, 1:27:01 AM3/18/12
to
Джентльмены!

On Saturday January 28 2012 you wrote to Michael Yuzhakov:

VK> Л.А.Залманзон. Преобразование Фурье, Уорша, Харра и их применение
VK> Легко яндексится
VK> Проще и понятнее не встречал

Дошли руки, почитал, спасибо. Впрочем, ничего особо нового в этот вумной
книжке не встретил, просто вспомнил институтский курс.

А теперь практический вопрос. Может, кто подскажет.

Имеется квазипериодический сигнал (с затуханием). Для оценки спектра
делается БПФ N последних отсчетов (N естественно равно 2 в целой степени). Hа
выходе имеем оценку спектра; при этом расстояние между частотами "гармоник"
обратно пропорционально числу отсчетов.
Хотелось бы оценить спектр более точно. Тем более что интересует собственно
не весь спектр, а какая-то его область. N увеличивать так просто не хочется,
ибо увеличение инерционности.
Пока мысли вслух. Или попробовать оценить период, под него подстроить
частоты "гармоник". Hо возможны погрешности из-за шумов. Кроме того, возможна
сумма нескольких квазипериодических сигналов. Или каким-то образом сместить
спектр так, чтобы анализировать интересующую его часть.

Valentin Kuznetsov

unread,
Mar 18, 2012, 4:33:04 PM3/18/12
to
��� �p���p� ������� `Spam detected`, � �p���p��, � ���������, ��� Michael Yuzhakov ����� Valentin Kuznetsov:

MY> ������� ������������������ ������ (� ����������). ���
MY> ������ ������� �������� ��� N ��������� �������� (N
MY> ����������� ����� 2 � ����� �������). H� ������ ����� ������
MY> �������; ��� ���� ���������� ����� ��������� "��������"
MY> ������� ��������������� ����� ��������.
MY> �������� �� ������� ������ ����� �����. ��� ����� ���
MY> ���������� ���������� �� ���� ������, � �����-�� ���
MY> �������. N ����������� ��� ������ �� �������, ��� ����������
MY> �������������.
MY> ���� ����� �����. ��� ����������� ������� ������, ���
MY> ���� ���������� ������� "��������". H� �������� �����������
MY> ��-�� �����. ����� ����, �������� ����� ����������
MY> ������������������ ��������. ��� �����-�� ������� ��������
MY> ������ ���, ����� ������������� ������������ ��� �����.
��� ��� ��� ��� ����� � ����� ���� ��������� ����
���� ���������� ������� �����, �� ������ ������� ���� ��� ������ ������� � �ݣ
����� ��������, � ����� �������� �� ����� �������, ����� ����� �������, �����
����� ��������������� �������� ����� ����� ������� ��� ������ ��������
���������� ������ ������������
������ ������� ������� ���������� ���ޣ��� ���ģ���, � ��� ������������ ��
���� ����� �� "�� ������ ���������"

Valentin

Michael Yuzhakov

unread,
Mar 19, 2012, 1:58:51 AM3/19/12
to
�����������!

On Monday March 19 2012 you wrote to Michael Yuzhakov:

VK> ��� ��� ��� ��� ����� � ����� ���� ��������� ����
VK> ���� ���������� ������� �����, �� ������ ������� ���� ��� ������ �������
VK> �
VK> �ݣ ����� ��������, � ����� �������� �� ����� �������, ����� �����
VK> �������,
VK> ����� ����� ��������������� �������� ����� ����� ������� ��� ������
VK> �������� ���������� ������ ������������ ������ ������� ������� ����������
VK> ���ޣ��� ���ģ���, � ��� ������������ �� ���� ����� �� "�� ������
VK> ���������"

����� ������ �������, ����� �������� �����/������ ������?
��� ����� �����, � ��, ��� �������� �����, ���������� �� ���������, � ��
������������ ���� ����� ���� �� ���.

� ������� ��������, ������ ������.
_________________________________________________________
�-����: yuzhakov<������>etel.ru

Valentin Kuznetsov

unread,
Mar 20, 2012, 4:47:58 AM3/20/12
to
Мои пpибоpы сказали `Spam detected`, я пpовеpил, а оказалось, это Michael Yuzhakov пишет Valentin Kuznetsov:

VK> Вот как раз для этого в книге дано несколько схем
VK> Если нарисовать сложную схему, на выходе которой есть все нужные сигналы
VK> и
VK> ещё масса ненужных, а потом ненужные из схемы стереть, легко будет
VK> увидеть,
VK> какие части вычислительного процесса можно будет удалить без потери
VK> точности дешифровки нужных составляющих Тоесть хранить большОе количество
VK> отсчётов придётся, а вот обрабатывать их надо будет не "по полной
VK> программе"

MY> Может ткнешь пальцем, какое ключевое слово/раздел
MY> искать?
Следует освоить разницу между "обычными" и "быстрыми" преобразованиями Фурье
По сути математики - это одно и то же, но "быстрое" требует меньше
вычислительной работы из-за рационализации схемы вычислений
Для понимания разницы очень полезна эта самая схема вычислений
Если разница станет понятна, то станет понятна и суть этих самых схем
Далее следует сварганить схему, наиболее подходящую для поставленной задачи

Valentin

Michael Yuzhakov

unread,
Apr 20, 2012, 8:54:34 AM4/20/12
to
On Tuesday March 20 2012 you wrote to Michael Yuzhakov:

VK> Следует освоить разницу между "обычными" и "быстрыми" преобразованиями
VK> Фурье По сути математики - это одно и то же, но "быстрое" требует меньше
VK> вычислительной работы из-за рационализации схемы вычислений Для понимания
VK> разницы очень полезна эта самая схема вычислений Если разница станет
VK> понятна, то станет понятна и суть этих самых схем Далее следует сварганить
VK> схему, наиболее подходящую для поставленной задачи

Hу вот, руки дошли, написал свою процедуру FFT. Берем последовательность,
разделяем на четную / нечетную, для каждой подпоследовательности рекуррентно
вызываем ту же процедуру FFT; получив результат, из преобразования четной и
нечетной последовательности компонуем преобразование текущей
последовательности. Проверил - работает.
Однако просветленным пока себя не чувствую. Или может мы друг друга не
совсем поняли.
Сейчас-то пока идет речь не об оптимизации количества вычислений, пока
занимаюсь моделированием, на скорость плевать.
Речь идет о том, чтобы из преобразования N последних отсчетов сигнала
вытянуть наиболее точную информацию, оценив частоты/фазы гармоник.
В процессе работы пойдет палка о двух концах: с одной стороны, точность
оценки частот зависит от количества отсчетов сигнала, на которых делаем
преобразование; с другой стороны, от этого же зависит инерционность (система
предполагается реалтаймовая).
Конкретные задачи:
1. Предельно точно для данного N оценить частоты гармоник сигнала (которые в
общем случае не совпадают с частотами членов ряда).
2. Предельно точно оценить амплитуду и особенно фазу для найденных в п.1
гармоник.
3. Избавиться от влияния лепестков, которые возникают в первое время начала
сигнала, когда буфер еще частично заполнен нулями.
Вот это и буду прорабатывать. Какие-то сведения нашел в
http://psi-logic.narod.ru; если ткнете еще куда, буду признателен.
0 new messages