Отправлено GroundCtrl 16 марта 2003 г. 18:43 В ответ на: Ничо не понял. Показания с акселерометра разложить FFT. Взять 0ую и 1ую гармоники. Сдлелать над ними обратное FFT и проинтегрировать? отправлено goshka 16 марта 2003 г. 18:25

Щас проделал все выкладки. Вот что получилось.

Где N - число отсчётов, Sn[x(n)] - сумма x(n) при n от 0 до N-1,
W(nk) = exp(j*2*pi*n*k/N), a(n) - отсчёты ускорения, Т - период выборки ускорения.

Приращение скорости:

dV = T*Sn[a(n)]

ДПФ от N отчётов ускорения:

A(n) = Sk[a(k)*W(-nk)]
a(n) = (1/N)*Sk[A(k)*W(nk)]

Приращение скорости через ДПФ

dV = (T/N)*Sn[Sk[A(k)*W(nk)]]=(T/N)*Sk[A(k)*Sn[W(nk)]]

т.к. W(nk) = N, если k = m*N и нулю в остальных случаях, то

dV = (T/N)*A(0)*N = T*Sk[a(k)]

Т.е. никакие искусственные изменения в отсчётах ДПФ на интеграл не повлияют. Единственные способ избавиться от периодической вибропомехи, в данном случае, это сделать окно кратным периоду этой помехи, и достаточно длинным. И ДПФ, в чистом виде, оказывается ненужна. Надо просто суммировать длинную последовательность отсчётов.

• Исправление — GroundCtrl   (16.03.2003 18:44, 56 байт)