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

 
 
> Усрднение значений фазы
amusin
сообщение Aug 20 2010, 06:38
Сообщение #1


Частый гость
**

Группа: Участник
Сообщений: 120
Регистрация: 2-09-05
Из: Екатеринбург
Пользователь №: 8 165



Есть измерения фазы в диапазоне 0...359 градусов.
Допустим, показания болтаются около отметки 0 градусов:
половина значений в первом квадранте, половина - в четвертом.
Среднее арифметическое даст примерно 180 градусов, что неверно.

Пока видится одно решение:
сохраняем массив измерений,
доворачиваем измерения настолько, чтобы сектор "болтанки" не лежал на угле 0 градусов.
Решение не нравится тем, что надо копить данные в ОЗУ.

Есть более красивые варианты?

Сообщение отредактировал amusin - Aug 20 2010, 06:38
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
_4afc_
сообщение Aug 20 2010, 07:25
Сообщение #2


Профессионал
*****

Группа: Свой
Сообщений: 1 262
Регистрация: 13-10-05
Из: Санкт-Петербург
Пользователь №: 9 565



Цитата(amusin @ Aug 20 2010, 10:38) *
Есть измерения фазы в диапазоне 0...359 градусов.
Допустим, показания болтаются около отметки 0 градусов:
половина значений в первом квадранте, половина - в четвертом.
Среднее арифметическое даст примерно 180 градусов, что неверно.

Пока видится одно решение:
сохраняем массив измерений,
доворачиваем измерения настолько, чтобы сектор "болтанки" не лежал на угле 0 градусов.
Решение не нравится тем, что надо копить данные в ОЗУ.

Есть более красивые варианты?


Вариант 1:

Если, по по вашему мнению, доворачивание помогает - доворачивайте сразу перед сохранением массива.
Тогда, как я понимаю массив в вашем алгоритме не понадобится.

Если бы мне нужна была точность в 1 градус - я бы завёл массив из 360 счётчиков и увеличивал их необходимое время.
Далее хоть максимумы берите, хоть распределения и дисперсии смотрите.
-----------------------------------------------------------------------------------------
Вариант 2:

2.0 заводим 3 переменных S,N,F

2.1 изменить представление диапазона с 0...359, на -179...180
2.2 за 0 принять результат предыдущего измерения фазы (F)
2.3 преобразовать измеренное значение в представление 2.1+2.2

2.4 Цикл S=S+измеренное; N++.

2.5 F=F+S/N;

Но лучше ещё добавить скользящее окно.
Тогда к этим трём переменным добавится массив длинной в окно (кольцевой буфер), а N превратится в указатель текущего отсчёта в кольцевом буфере.

Соответственно S=S+измеренное-B[N];B[N]=измеренное;(N++)%длинна окна. F=F+S/длинна окна;

Для борьбы с накоплением ошибки - я бы добавил стекание типа F=F-F/1024;

PS: Никогда не усреднял фазу, но вроде должно сработать smile.gif
Go to the top of the page
 
+Quote Post
_4afc_
сообщение Aug 20 2010, 08:56
Сообщение #3


Профессионал
*****

Группа: Свой
Сообщений: 1 262
Регистрация: 13-10-05
Из: Санкт-Петербург
Пользователь №: 9 565



Цитата(_4afc_ @ Aug 20 2010, 11:25) *
Для борьбы с накоплением ошибки - я бы добавил стекание типа F=F-F/1024;


Вернее S=S-S/1024; конечно же wink.gif

Цитата(Oldring @ Aug 20 2010, 12:17) *
Мне кажется, что если я после двух прочтений так и не понял, что именно вы хотите сделать - это всё слишком сложно, и поэтому, малореализуемо. laughing.gif


Не вижу сложностей в алгоритме, возможно вас пугает детальное описание - имелось ввиду следующее:

Я предлагаю работать не со значениями фазы, а с разницей между измеренным значением фазы (А) и значением фазы вычисленное в предыдущий момент (F) и выводимое например на индикатор.
Т.е. B=f(A,F). Функция f - не просто B=A-F. Она ещё приводит значение B в диапазон -179...+180.
Числа B усредняем скользящим окном. Накопление ошибок в скользящем окне устраняем стеканием.F - выводим на экран.
Go to the top of the page
 
+Quote Post
Oldring
сообщение Aug 20 2010, 11:11
Сообщение #4


Гуру
******

Группа: Свой
Сообщений: 3 041
Регистрация: 10-01-05
Из: Москва
Пользователь №: 1 874



Цитата(_4afc_ @ Aug 20 2010, 12:56) *
Не вижу сложностей в алгоритме, возможно вас пугает детальное описание - имелось ввиду следующее:

Я предлагаю работать не со значениями фазы, а с разницей между измеренным значением фазы (А) и значением фазы вычисленное в предыдущий момент (F) и выводимое например на индикатор.
Т.е. B=f(A,F). Функция f - не просто B=A-F. Она ещё приводит значение B в диапазон -179...+180.
Числа B усредняем скользящим окном. Накопление ошибок в скользящем окне устраняем стеканием.F - выводим на экран.


Например, тут нужен некоторый механизм запуска, формирующий начальное значение F.
"Постоянно обновляемый индикатор" - это тоже домыслы. Человек вроде хотел простое средне по блоку отсчетов.
Да и описано сложно для понимания словами, а не строгим алгоритмом. Куда идет "В" вы так и не описали. Вы с ним собрались зачем-то "работать", но результат нигде не используется в вашем неформальном описании, нужно догадываться, что именно вы хотели сказать.

Цитата(fontp @ Aug 20 2010, 14:54) *
допустим,
mean( unwrap( phi ) )=0, если фаза болтается в окрестности pi и около

Думаю, автор сказал, про "допустим" чтобы разобраться, а не потому, что это априорно известно.
Если априорно известно, где болтается фаза, то нет проблем, я об этом сказал
А вы Найтман, опять нахамили



Неоднократно вам уже писал, что я не Найтман. Поэтому вы - очевидно паранойик, раз не верите этому моему утверждению.
Считать, что никогда не может писать "чушь" может только самоуверенный идиот, считающий, что он не может ошибаться. Чушь пишут все периодически. Выбирайте уж сами на будущее, кем вам быть милее.

А именно чушь вы написали и еще раз. Достаточно запустить в Матлабе команду

Код
mean(unwrap(pi+1e-2*randn(1000,1)))


чтобы в этом убедиться каждому самостоятельно и независимо от моего субъективного мнения.

И еще раз повторю. Знать, какое среднее фазы, априорно совершенно не нужно. Нужно априорно знать, что разброс не чрезмерно велик.


--------------------
Пишите в личку.
Go to the top of the page
 
+Quote Post
fontp
сообщение Aug 20 2010, 11:32
Сообщение #5


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



QUOTE (Oldring @ Aug 20 2010, 15:11) *
CODE
mean(unwrap(pi+1e-2*randn(1000,1)))


чтобы в этом убедиться каждому самостоятельно и независимо от моего субъективного мнения.

И еще раз повторю. Знать, какое среднее фазы, априорно совершенно не нужно. Нужно априорно знать, что разброс не чрезмерно велик.


Посмотрел реализацию unwrap. Здесь был неправ, этот именно тот последовательный алгоритм, что я привёл выше. И так он и будет работать
Я думал, что это просто приведение к стандартному диапазону.
4afc на самом деле предложил примерно то же самое
Go to the top of the page
 
+Quote Post
Oldring
сообщение Aug 20 2010, 11:47
Сообщение #6


Гуру
******

Группа: Свой
Сообщений: 3 041
Регистрация: 10-01-05
Из: Москва
Пользователь №: 1 874



Цитата(fontp @ Aug 20 2010, 15:32) *
Посмотрел реализацию unwrap. Здесь был неправ, этот именно тот последовательный алгоритм, что я привёл выше. И так он и будет работать



Про "распараллеливание" речь вообще не шла. У человека последовательный микроконтроллер, отсчеты он получает последовательно, и даже не хочет их все хранить для дальнейшей обработки. Человек искал именно последовательный алгоритм. Если же кому-то нужно именно хорошо распараллеливаемый алгоритм - unwrap можно привести к таковому, вычислив первую конечную разность и переставив индексы в двойной сумме, в которую он входит.

Цитата(thermit @ Aug 20 2010, 15:35) *
Усреднение fmean(n) = fmean(n-1) + ( f(n) - fmean(n-1) )*k


Это не усреднение, а фильтрация IIR фильтром первого порядка.


--------------------
Пишите в личку.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- amusin   Усрднение значений фазы   Aug 20 2010, 06:38
- - Oldring   Цитата(amusin @ Aug 20 2010, 10:38) Есть ...   Aug 20 2010, 07:05
|- - amusin   Цитата(Oldring @ Aug 20 2010, 13:05) Спа...   Aug 20 2010, 07:20
- - Fat Robot   усрядняете cos(phi) и sin(phi). в конце - 4-х квад...   Aug 20 2010, 07:09
|- - Oldring   Цитата(_4afc_ @ Aug 20 2010, 11:25) Вариа...   Aug 20 2010, 08:17
|- - fontp   QUOTE (_4afc_ @ Aug 20 2010, 12:56) Верне...   Aug 20 2010, 10:11
||- - Oldring   Цитата(fontp @ Aug 20 2010, 14:11) первый...   Aug 20 2010, 10:43
||- - fontp   QUOTE (Oldring @ Aug 20 2010, 14:43) Чушь...   Aug 20 2010, 10:54
- - thermit   Цитатаamusin: Есть измерения фазы в диапазоне 0......   Aug 20 2010, 11:35
- - thermit   ЦитатаOldring: Это не усреднение, а фильтрация IIR...   Aug 20 2010, 11:55
|- - Oldring   Цитата(thermit @ Aug 20 2010, 15:55) Усре...   Aug 20 2010, 11:58
|- - _4afc_   Ребят. Я в научных терминах не силён. Но предложен...   Aug 20 2010, 12:28
|- - Oldring   Цитата(_4afc_ @ Aug 20 2010, 16:28) Код п...   Aug 20 2010, 12:48
|- - _4afc_   Цитата(Oldring @ Aug 20 2010, 16:48) Во-п...   Aug 20 2010, 13:13
|- - Oldring   Цитата(_4afc_ @ Aug 20 2010, 17:13) Обычн...   Aug 20 2010, 13:18
- - Stanislav   Цитата(amusin @ Aug 20 2010, 10:38) Есть ...   Aug 27 2010, 20:42


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

 


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


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