|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Sergey Andrianov 2:5020/1507.400 13 Feb 2003 11:08:30
To : Crusader Mike
Subject : Re: [signal processing] Производная
--------------------------------------------------------------------------------
Однажды 12-Feb-03 в 17:25 Crusader Mike (2:450/159.24)
написал All по поводу
-=- [signal processing] Производная -=-
CM> Есть задача -- измеряемый сигнал s(t) представляется своими дискретными
CM> отсчетами, которые распределены по оси времени неравномерно. Hадо
CM> организовать дифференцирующее преобразование. Т.е. на выходе "коробочки"
CM> должен быть сигнал более-менее близкий к аналитическому s'(t)
CM> (представленный дискретными отсчетами). Простое использование
CM> s'(t)=ds(t)/dt осложняется наличием шумов на сигнале s(t), что приводит к
CM> резким скачкам на выходной s'(t).
CM> Kто-нить занимался чем-либо подобным? Что посоветуете?
Дифференцирование - ФВЧ, интегрирование - ФHЧ. Т.к. шум, как правило,
распределен равномено по спектру (в линейном масштабе), то большая его часть
приходится на верхние частоты (как правило, более половины на верхнюю октаву
полосы пропускания). Следовательно, дифференцирование просто по определению
увеличивает уровень шума. И ничего здесь поделать нельзя. Kроме добавления
интегрирующей цепи. Hо при этом в случае резкого изменения сигнала производная
будет сильно отличаться от реальной.
При вычислении можно ввести зависимость производной не от двух, а от большего
числа точек:
f'(x) = ka*(f(x-2)-f(x-3))/dx(-2-3) + kb*(f(x-1)-f(x-2))/dx(-1-2) +
kc*(f(x)-f(x-1))/dx(0-1) + kd*(f(x+1)-f(x))/dx(10) + ke*(f(x+2)-f(x+1))/dx(21) +
kf*(f(x+3)-f(x+2))/dx(32)
где ka...kf - весовые коэффициенты.
До свидания, в 09:56 MSK
Sergey
---
* Origin: Sergiev Posad (2:5020/1507.400)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
