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

 
 
> Определение постоянной составляющей синусоиды...
vladimir_m
сообщение Nov 21 2008, 13:03
Сообщение #1


Участник
*

Группа: Новичок
Сообщений: 17
Регистрация: 9-07-08
Пользователь №: 38 830



Задача такая.
Есть сигнал - синусоида с постоянной составляющей, которая может медленно изменяться.
Синусоида оцифровывается и передаетя в комп буферами. Период синусоиды большой, заранее не известный, намного больше фиксированного размера буфера. Необходимо отслеживать изменение постоянной составляющей по приходу каждого буфера (уточнять с каждым буфером).
Как такое можно сделать?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
mvb
сообщение Nov 27 2008, 00:07
Сообщение #2


Участник
*

Группа: Участник
Сообщений: 72
Регистрация: 7-06-08
Пользователь №: 38 128



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

Сначала сформулирую проблему, которую решал: есть N отсчетов, в которых присутствует синусоида + постоянная составляющая + шум. Период синусоиды заведомо в разы больше N, начальная фаза случайная. Шум -- норм. распределенная случ. величина с нулевым мат. ожиданием.

Я попробовал два озвученных решения: ФНЧ и Станислава. ДПФ и АКФ не смотрел, потому что не понимаю как они здесь помогут.

ФНЧ брал 1й Чебышева. Варьировал все возможные параметры как мог, но не смог добиться того, чтобы на выходе как-то обозначилась искомая величина. Да, выход по времени рассматривал достаточно длинный.

Для реализации метода, предложенного Станиславом взял Curve Fitting toolbox из матлаба. При бесконечном ОСШ он точно выдает искомую величину. В присутствии шума решение начинает гулять. Только надо выставлять область ограничений для решения не бесконечную, иначе минимумов у целевой функции получается несколько.

На все это я потратил 1 час 7 минут.

Вопросы к знающим: 1. как именно вы предлагаете применять фнч?
2. Как здесь можно использовать дпф или акф?

Отдельная просьба: ваше мнение о моих способностях, оставьте при себе.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- vladimir_m   Определение постоянной составляющей синусоиды...   Nov 21 2008, 13:03
- - D.I.M.A   Оценивать мат. ожидание процесса.   Nov 21 2008, 13:16
- - rezident   А постоянная составляющая по какому-то закону изме...   Nov 21 2008, 13:24
- - vladimir_m   ЦитатаА постоянная составляющая по какому-то закон...   Nov 21 2008, 13:34
|- - sysel   По идее, для получения постоянной составляющей над...   Nov 21 2008, 13:49
- - Самурай   Цитата(vladimir_m @ Nov 21 2008, 16:03) З...   Nov 21 2008, 14:12
- - vladimir_m   Большое окно - это хорошо. Расчетом занимается и т...   Nov 21 2008, 14:17
|- - fontp   Цитата(vladimir_m @ Nov 21 2008, 17:17) Б...   Nov 21 2008, 15:07
|- - Serg76   Цитата(fontp @ Nov 21 2008, 18:07) Фильтр...   Nov 26 2008, 16:34
- - ivan219   Зделать FFT и 0 отщёт будет = пост состовляющей.   Nov 21 2008, 19:57
|- - GetSmart   Цитата(ivan219 @ Nov 22 2008, 01:57) Здел...   Nov 24 2008, 22:23
- - ivan219   1   Nov 21 2008, 19:57
- - Stanislav   Цитата(ivan219 @ Nov 21 2008, 22:57) Здел...   Nov 21 2008, 20:36
- - Михаил_K   А чем не подходит это http://electronix.ru/forum/i...   Nov 22 2008, 10:01
|- - Stanislav   Цитата(Михаил_K @ Nov 22 2008, 13:01) А ч...   Nov 22 2008, 12:53
|- - Михаил_K   Цитата(Stanislav @ Nov 22 2008, 15:53) А ...   Nov 26 2008, 04:58
|- - Stanislav   Цитата(Михаил_K @ Nov 26 2008, 07:58) Что...   Nov 26 2008, 09:02
|- - Михаил_K   Цитата(Stanislav @ Nov 26 2008, 12:02) Мд...   Nov 26 2008, 09:04
|- - Stanislav   Цитата(Михаил_K @ Nov 26 2008, 12:04) Буд...   Nov 26 2008, 09:08
- - vladimir_m   Всем спасибо! Лично я остановился на варианте ...   Nov 22 2008, 12:31
|- - Tanya   Цитата(vladimir_m @ Nov 22 2008, 15:31) В...   Nov 22 2008, 12:58
|- - fontp   Цитата(Tanya @ Nov 22 2008, 15:58) Нескол...   Nov 24 2008, 08:17
|- - Stanislav   Цитата(fontp @ Nov 24 2008, 11:17) В прин...   Nov 24 2008, 08:30
|- - fontp   Цитата(Stanislav @ Nov 24 2008, 11:30) Ст...   Nov 24 2008, 08:41
|- - Stanislav   Цитата(fontp @ Nov 24 2008, 11:41) Просто...   Nov 24 2008, 09:01
|- - fontp   Цитата(Stanislav @ Nov 24 2008, 12:01) Кр...   Nov 24 2008, 09:51
|- - alexkok   Цитата(fontp @ Nov 24 2008, 12:51) Про ...   Nov 24 2008, 21:03
|- - Stanislav   Цитата(fontp @ Nov 24 2008, 12:51) Да лад...   Nov 25 2008, 00:31
- - GetSmart   Цитата(Stanislav)Для вычисления параметров синусои...   Nov 22 2008, 17:12
|- - Stanislav   Цитата(GetSmart @ Nov 22 2008, 20:12) На ...   Nov 22 2008, 18:09
- - GetSmart   Цитата(Stanislav)Послушайте, уважаемый Получатель ...   Nov 22 2008, 23:38
|- - Stanislav   Цитата(GetSmart @ Nov 23 2008, 02:38) Тол...   Nov 23 2008, 00:14
- - GetSmart   Цитата(Stanislav)Возможность практической реализац...   Nov 23 2008, 22:24
|- - Stanislav   Цитата(GetSmart @ Nov 24 2008, 01:24) Я у...   Nov 23 2008, 23:20
|- - GetSmart   Цитата(Stanislav @ Nov 24 2008, 05:20) А ...   Nov 23 2008, 23:45
|- - Stanislav   Цитата(GetSmart @ Nov 24 2008, 02:45) Ну ...   Nov 23 2008, 23:49
- - GetSmart   По приходу нового буфера берётся несколько последн...   Nov 24 2008, 23:54
- - Stanislav   ЗЫ. fontp, я готов прекратить спор, если будет при...   Nov 25 2008, 07:30
|- - fontp   Цитата(alexkok @ Nov 25 2008, 00:03) Не о...   Nov 25 2008, 07:54
|- - Stanislav   Цитата(fontp @ Nov 25 2008, 10:54) Fittin...   Nov 25 2008, 08:23
|- - alexkok   Цитата(fontp @ Nov 25 2008, 10:54) Fittin...   Nov 25 2008, 17:52
|- - blackfin   Цитата(alexkok @ Nov 25 2008, 20:52) Fitt...   Nov 26 2008, 05:23
- - GetSmart   Цитата(alexkok)Fitting, как Вам уже ответили, при ...   Nov 25 2008, 23:17
|- - alexkok   Цитата(GetSmart @ Nov 26 2008, 02:17) А э...   Nov 26 2008, 08:10
- - GetSmart   Цитата(blackfin)Для выделения полезного сигнала на...   Nov 26 2008, 07:03
|- - alex_os   Цитата(GetSmart @ Nov 26 2008, 10:03) При...   Nov 26 2008, 08:25
- - RadioJunior   пост не в тему... Stanislav, не обижайтесь, но эт...   Nov 26 2008, 09:11
|- - Stanislav   Цитата(RadioJunior @ Nov 26 2008, 12:11) ...   Nov 26 2008, 21:29
|- - Михаил_K   Цитата(RadioJunior @ Nov 26 2008, 12:11) ...   Nov 28 2008, 05:00
|- - Stanislav   Цитата(mvb @ Nov 28 2008, 00:10) Вот мой ...   Nov 28 2008, 06:55
- - GetSmart   Цитата(RadioJunior)Со мной кто-нибудь согласен по ...   Nov 26 2008, 10:43
|- - Stanislav   Цитата(mvb @ Nov 27 2008, 03:07) Для реал...   Nov 27 2008, 00:23
- - GetSmart   Цитата(mvb)2. Как здесь можно использовать дпф или...   Nov 27 2008, 04:12
- - mvb   Вот мой код, без претензий на робастность и скорос...   Nov 27 2008, 21:10
- - mvb   Если бы у меня был бы фрейм с достаточным количест...   Nov 28 2008, 10:03
|- - Stanislav   Цитата(mvb @ Nov 28 2008, 13:03) Если бы ...   Nov 28 2008, 11:45
- - GetSmart   Цитата(mvb)2. Далее скомпенсировал бы ошибку в ДПФ...   Nov 28 2008, 13:08
|- - Stanislav   Цитата(GetSmart @ Nov 28 2008, 16:08) Да,...   Nov 28 2008, 20:59
- - mvb   Попробовал на досуге определить период синуса в ма...   Nov 30 2008, 09:21


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

 


RSS Текстовая версия Сейчас: 24th July 2025 - 11:44
Рейтинг@Mail.ru


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