Отправлено -=ВН=- 19 июля 2005 г. 12:13 В ответ на: Намекните на алгоритм: Как корректнее посчитать RMS от 50Гц-сигнала , при условии что на АЦП сигнал приходит со смещением в середину диапазона? отправлено bp 19 июля 2005 г. 10:41

Ну а так:
SW=SUM(W[n])
SW2=SUM(W[n]^2)
S1=SUM(W[n]*X[n])/SW
S2=SUM((W[n]*X[n])^2)/SW2
RMS=SQRT(S2-S1^2).

SW и SW2 считаются раз и заранее.
Для S1 и S2 налаживается легко рекурсия.

W - окно, обычно из серии "какое вспомнилось, или нашлось". Что- нибудь из классики :-) Ну или Гегенбауэра. Ужасно легко запоминается, 1 (адын) параметр, ужасно легко генерируется. А в зависмости от парметра оч. близкое к классике можно получить.
W[n]=[(n+1)*(N-n)*4/(N+1)]^p. N - длина, p - параметр.

Либо, если считаете, что можно пренебречь погрешностью за счет некратности времени измерения и периода 50 Гц (сеть наверное), то окно прямоугольное. W[n]=1. SW=SW2=N.

• А корректно ли это накладывать окно на сигнал со смещением? — bp   (19.07.2005 12:40 82.140.78.21, пустое)
  • А почему нет? Для Вашей задачи? У Вас вариантов ведь всего 2. — -=ВН=-   (19.07.2005 12:58 194.190.181.231, 494 байт)
    • Посчитал в матлабе для стационарного процесса, довольно хорошо сошлось. Буду в железе реализовывать. Спасибо за совет! — bp   (19.07.2005 13:06 82.140.78.21, пустое)
• А корректно ли это накладывать окно на сигнал со смещением? — bp   (19.07.2005 12:40 82.140.78.21, пустое)