Перейти к содержанию
    

Найти экстремумы функции с гармонической составляющей

После чего находим частоту:

 

f = [1/[pi*(t3 - t1)]]*arccos[(S3 + S1)/(2*S2)].

 

Как-то так..

Решение красивое, но что будет, если а) добавить шум, б) добавить интерференцию

Можно привести пример. Из пространственного спектроанализа, но сути это не меняет. В АФАР можно определить направление на цель всего по 2-м сенсорам, оценив задержку сигнала, полученного с одного сенсора относительно другого. Однако эта схема работает с относительно сильными и узкополосными сигналами, не содержащими интерференцию. Если же это условие не выполняется (что происходит довольно часто), то придётся а) наращивать количество сенсоров (т.е. количество сэмплов для анализа) и б) усложнять алгоритм (в отличие от временного спектроанализа, где мы в основном можем увеличить длину выборки, в пространственном размер решётки ограничен, т.е. Фурье с его малым разрешением неприемлим).

 

Марпл , Хайкин... А всё равно кошерных методов "и/или разделения близко расположенных по частоте сигналов" - не существует!

Предлагаю рассмотреть элементарную задачу - протонный магнитометр. За 0.5 - 1 сек предлагается определить частоту прецессии (2-3 кГц) с точностью 0.001 Гц

Давайте рассмотрим аналогичную задачу (но сделаем её ещё сложнее):

post-81866-1438326399_thumb.png

символьная частота = 100 кГц, ОСШ = 60 дБ, частота сигнала = 80,348 Гц, количество символов для анализа = 1024 + 128 (как видно из картинки - менее одного периода), погрешность ~ 0,001 Гц.

Изменено пользователем serjj

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

...Если такой сигнал можно рассматривать как комплексную экспоненту с шумом, то все ваши фантазии будут строго ограничены пределом Крамера-Рао...

Ну мои сексуальные фантазии так далеко не заходят... :-)

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Ну мои сексуальные фантазии так далеко не заходят... :-)

печалька :(, не богатая значит у Вас фантазия, трудно подобрать более простой сигнал, чем exp(jwt) :)

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Ну тогда порадуй старика новым оригиальным решением :-)

все старо как мир, читай выше, кроме того, никто поиск по форуму не отменял

PS. blackfin привел одно из простых решений линейной интерполяцией, существует масса других, более/менее сложных с соответствующей точностью оценки.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

blackfin привел одно из простых решений линейной интерполяцией, существует масса других, более/менее сложных с соответствующей точностью оценки.

Страшно далеки они от народа...(с) :-)

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Страшно далеки они от народа...(с) :-)

«Паситесь, мирные народы!

Вас не разбудит чести клич.

К чему стадам дары свободы?

Их должно резать или стричь.

Наследство их из рода в роды

Ярмо с гремушками да бич.» ©

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Решение красивое, но что будет, если а) добавить шум, б) добавить интерференцию

Считать по этой формуле частоту с относительной точностью ∆f/f можно только если СКО шума меньше чем:

 

σ < (S2/6)*(∆f/f),

 

и при этом:

 

abs(S1 + S3) < abs(S2).

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Присоединяйтесь к обсуждению

Вы можете написать сейчас и зарегистрироваться позже. Если у вас есть аккаунт, авторизуйтесь, чтобы опубликовать от имени своего аккаунта.

Гость
Ответить в этой теме...

×   Вставлено с форматированием.   Вставить как обычный текст

  Разрешено использовать не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отображать как обычную ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставлять изображения напрямую. Загружайте или вставляйте изображения по ссылке.

×
×
  • Создать...