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

 
 
> Как убрать постоянную составляющую из сигнала., Сингал 50Гц
Pat
сообщение Oct 31 2005, 11:12
Сообщение #1


Местный
***

Группа: Свой
Сообщений: 421
Регистрация: 27-05-05
Из: Энергодар
Пользователь №: 5 480



Сигнал 50 Гц сдвинут относительно нуля.
Измерение производим АЦП с частотой выборки 5000Гц.
Надо получить его же только относительно нуля.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
bmf
сообщение Oct 31 2005, 15:23
Сообщение #2


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

Группа: Свой
Сообщений: 146
Регистрация: 6-01-05
Из: Украина
Пользователь №: 1 831



Поcмотри:
простейший HP фильтр
dc_bias - 32 bit integer
clean_sample - 16 bit integer

dc_bias += ((((int32_t) sample << 15) - dc_bias) >> 14);
clean_sample = sample - (dc_bias >> 15);

думаю поможет
Go to the top of the page
 
+Quote Post
bmf
сообщение Oct 31 2005, 17:32
Сообщение #3


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

Группа: Свой
Сообщений: 146
Регистрация: 6-01-05
Из: Украина
Пользователь №: 1 831



Цитата(bmf @ Oct 31 2005, 18:23)
    dc_bias += ((((int32_t) sample << 15) - dc_bias) >> 14);
    clean_sample = sample - (dc_bias >> 15);

примерно так

HP - high pass filter

формула для 1-го порядка:
y(n) = x(n) - x(n-1) + R * y(n-1)

"R" 0.9 .. 1
для -3dB: R = 1 - (pi*2 * frequency /samplerate)

ну а то что в выше - просто реализация что бы небыло дробного умножения

немного теории:
http://ccrma.stanford.edu/~jos/filters/
http://ccrma.stanford.edu/~jos/filters/DC_Blocker.html
Go to the top of the page
 
+Quote Post
bmf
сообщение Oct 31 2005, 17:48
Сообщение #4


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

Группа: Свой
Сообщений: 146
Регистрация: 6-01-05
Из: Украина
Пользователь №: 1 831



Цитата(bmf @ Oct 31 2005, 20:32)
y(n) = x(n) - x(n-1) + R * y(n-1)

может более привычный вид
1-st Order Filter of the Form:
y(n)=b0*x(n)+b1*x(n-1)+a1*y(n-1)
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Pat   Как убрать постоянную составляющую из сигнала.   Oct 31 2005, 11:12
- - moonrock   Вам нужно до АЦП сдвинуть или уже из отсчетов убра...   Oct 31 2005, 11:51
|- - Pat   Цитата(moonrock @ Oct 31 2005, 14:51)Вам нужн...   Oct 31 2005, 12:30
|- - Vic   Можно не накапливать а вычитать на каждый отсчет, ...   Oct 31 2005, 14:23
|- - Pat   Цитата(Vic @ Oct 31 2005, 17:23)Можно не нака...   Oct 31 2005, 14:38
|- - ms1   Цитата(Pat @ Oct 31 2005, 17:38)Цитата(Vic ...   Oct 31 2005, 18:24
|- - Pat   Цитата(ms1 @ Oct 31 2005, 21:24)Лучше всего и...   Oct 31 2005, 18:57
|- - vm1   А какой у Вас АЦП? Есть АЦП с внутренним сдвигом, ...   Oct 31 2005, 21:20
- - BVU   Цитата(Pat @ Oct 31 2005, 15:12)Сигнал 50 Гц ...   Oct 31 2005, 11:57
- - VAI   ЦитатаУровень то известен, но он тоже не совершене...   Oct 31 2005, 15:01
- - Vic1   В алгоритмах ничего лучше вычисления среднего и не...   Oct 31 2005, 15:23
|- - Vic1   Цитата(bmf @ Oct 31 2005, 22:32)Цитата(bmf ...   Nov 1 2005, 05:39
|- - Stanislav   Цитата(bmf @ Oct 31 2005, 20:32)формула для 1...   Nov 1 2005, 19:04
- - Vic1   To bmf А поподробнее для всех и автора темы можно?...   Oct 31 2005, 16:12
- - Vic1   Так это же совсем другой алгоритм (или я ничего не...   Nov 2 2005, 10:27
|- - Stanislav   Цитата(Vic1 @ Nov 2 2005, 13:27)Так это же со...   Nov 2 2005, 13:05
- - Vic1   To Stanislav Спасибо! С пунктами просто замеча...   Nov 2 2005, 14:41
|- - Stanislav   Цитата(Vic1 @ Nov 2 2005, 17:41)To Stanislav ...   Nov 2 2005, 15:26
- - Vic1   Я имела в виду, что помимо погрешности дискретизац...   Nov 2 2005, 15:42
|- - Stanislav   Цитата(Vic1 @ Nov 2 2005, 18:42)Я имела в вид...   Nov 2 2005, 17:08
- - mbc   Цитата(Pat @ Oct 31 2005, 16:12) Сигнал 5...   Nov 30 2005, 19:00


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

 


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


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