|
Усрднение значений фазы |
|
|
|
 |
Ответов
|
Aug 20 2010, 07:25
|

Профессионал
    
Группа: Свой
Сообщений: 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: Никогда не усреднял фазу, но вроде должно сработать
|
|
|
|
|
Aug 20 2010, 08:56
|

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

|
Цитата(_4afc_ @ Aug 20 2010, 11:25)  Для борьбы с накоплением ошибки - я бы добавил стекание типа F=F-F/1024; Вернее S=S-S/1024; конечно же  Цитата(Oldring @ Aug 20 2010, 12:17)  Мне кажется, что если я после двух прочтений так и не понял, что именно вы хотите сделать - это всё слишком сложно, и поэтому, малореализуемо.  Не вижу сложностей в алгоритме, возможно вас пугает детальное описание - имелось ввиду следующее: Я предлагаю работать не со значениями фазы, а с разницей между измеренным значением фазы (А) и значением фазы вычисленное в предыдущий момент (F) и выводимое например на индикатор. Т.е. B=f(A,F). Функция f - не просто B=A-F. Она ещё приводит значение B в диапазон -179...+180. Числа B усредняем скользящим окном. Накопление ошибок в скользящем окне устраняем стеканием.F - выводим на экран.
|
|
|
|
|
Aug 20 2010, 11:11
|

Гуру
     
Группа: Свой
Сообщений: 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))) чтобы в этом убедиться каждому самостоятельно и независимо от моего субъективного мнения. И еще раз повторю. Знать, какое среднее фазы, априорно совершенно не нужно. Нужно априорно знать, что разброс не чрезмерно велик.
--------------------
Пишите в личку.
|
|
|
|
|
Aug 20 2010, 11:32
|

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

|
QUOTE (Oldring @ Aug 20 2010, 15:11)  CODE mean(unwrap(pi+1e-2*randn(1000,1))) чтобы в этом убедиться каждому самостоятельно и независимо от моего субъективного мнения. И еще раз повторю. Знать, какое среднее фазы, априорно совершенно не нужно. Нужно априорно знать, что разброс не чрезмерно велик. Посмотрел реализацию unwrap. Здесь был неправ, этот именно тот последовательный алгоритм, что я привёл выше. И так он и будет работать Я думал, что это просто приведение к стандартному диапазону. 4afc на самом деле предложил примерно то же самое
|
|
|
|
|
Aug 20 2010, 11:47
|

Гуру
     
Группа: Свой
Сообщений: 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 фильтром первого порядка.
--------------------
Пишите в личку.
|
|
|
|
Сообщений в этой теме
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
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|