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

 
 
> минимизация погрешности, при плавающей частоте
TigerSHARC
сообщение Sep 4 2009, 10:48
Сообщение #1


Знающий
****

Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195



Есть сигнал. Частота основной гармоники 50 гц.
Используется алгоритм ДПФ. При смещении частоты на 1-5 Гц. Появляется дополнительная погрешность по модулю ДПФ. Оконное сглаживание не годится.
частота дискретизации 1200гц задана жёстко и не меняется, выборок строго 24 за период.
Слышал про алгоритм интерполяции, но так ничего конкретного не нашёл.
Может кто сталкивался или известен источник.????
Необходимо минимизировать погрешность БПФ при смещении частоты сигнала при постоянной частоте дискретизации.

Сообщение отредактировал TigerSHARC - Sep 4 2009, 10:49
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Oldring
сообщение Oct 21 2009, 06:27
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 3 041
Регистрация: 10-01-05
Из: Москва
Пользователь №: 1 874



Цитата(TigerSHARC @ Sep 4 2009, 14:48) *
Есть сигнал. Частота основной гармоники 50 гц.
Используется алгоритм ДПФ. При смещении частоты на 1-5 Гц. Появляется дополнительная погрешность по модулю ДПФ. Оконное сглаживание не годится.
частота дискретизации 1200гц задана жёстко и не меняется, выборок строго 24 за период.
Слышал про алгоритм интерполяции, но так ничего конкретного не нашёл.
Может кто сталкивался или известен источник.????
Необходимо минимизировать погрешность БПФ при смещении частоты сигнала при постоянной частоте дискретизации.


Вы забыли написать про время измерения. 24 точки на период, а сколько всего точек?


--------------------
Пишите в личку.
Go to the top of the page
 
+Quote Post
TigerSHARC
сообщение Oct 21 2009, 19:38
Сообщение #3


Знающий
****

Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195



Цитата(Oldring @ Oct 21 2009, 10:27) *
Вы забыли написать про время измерения. 24 точки на период, а сколько всего точек?


... алгоритм выполняется постоянно, оценка БПФ производиться каждый период сигнала(20 мс), оценка может производится сколь угодно долго, главное - 24 точечное БПФ, которое должно выполниться за 20 мс, получается так... каждые 20 мс вычисляется БПФ... так вот, если частота начнёт гулять, нужно что бы данные полученные от БПФ (модуль) не были искажены...
Go to the top of the page
 
+Quote Post
alexkok
сообщение Oct 22 2009, 01:16
Сообщение #4


Знающий
****

Группа: Участник
Сообщений: 609
Регистрация: 3-03-07
Из: San Jose
Пользователь №: 25 837



Цитата(TigerSHARC @ Oct 21 2009, 23:38) *
... алгоритм выполняется постоянно, оценка БПФ производиться каждый период сигнала(20 мс), оценка может производится сколь угодно долго, главное - 24 точечное БПФ, которое должно выполниться за 20 мс, получается так... каждые 20 мс вычисляется БПФ... так вот, если частота начнёт гулять, нужно что бы данные полученные от БПФ (модуль) не были искажены...

А зачем в таком случае БПФ?
Тут наверное что-то следящее больше подходит, типа цифровой ФАПЧ.


--------------------
Go to the top of the page
 
+Quote Post
TigerSHARC
сообщение Oct 27 2009, 15:10
Сообщение #5


Знающий
****

Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195



Цитата(alexkok @ Oct 22 2009, 05:16) *
А зачем в таком случае БПФ?
Тут наверное что-то следящее больше подходит, типа цифровой ФАПЧ.


...цифровой ФАПЧ? нужно программное решение...

Цитата(sup-sup @ Oct 22 2009, 08:12) *
Дык ведь, если доступен поток данных непрерывный. то и делаем интерполяционный фильтр с какой надо точностью (то есть, с необходимым нам увеличением частоты отсчетов, а потом анализируем. Проблем то нет. Лучше, все равно, не получится. Самое прямое и понятное решение. После интерполяции размер БПФ увеличенный берем, опять, такой какой нам нужен для точности.


можно подробнее про этот подход...

Сообщение отредактировал TigerSHARC - Oct 27 2009, 15:07
Go to the top of the page
 
+Quote Post
sup-sup
сообщение Oct 27 2009, 16:28
Сообщение #6


Знающий
****

Группа: Участник
Сообщений: 674
Регистрация: 26-08-05
Пользователь №: 7 997



Цитата(TigerSHARC @ Oct 27 2009, 18:10) *
...цифровой ФАПЧ? нужно программное решение...



можно подробнее про этот подход...

На пальцах так.
Предположим (для упрощения) что достаточно увеличить число точек в 10 раз.
Делаем интерполяцию любым доступным инструментом (Matlab, Excel, ..).
Для этого добавляем по 9 нулей между каждым отсчетом и "получаем частоту дискретизации" 12000 Hz.
Пропускаем этот "поток" через цифровой фильтр нижних частот (интерполяционный фильтр), задача которого - обязательно пропускать первую гармонику (50 Hz), но достаточно хорошо () подавлять все частоты выше 600 Гц. Можно взять для пробы затухание 40-60 дБ, а там видно будет. Нужно брать FIR, так как он не трогает фазы.
Выходной сигнал будет иметь в 10 раз больше точек, что позволит повысить "точность".
Если есть шумы, а они есть, то применяемый фильтр их подавит (вне первой гармоники), поэтому чем длиннее фильтр, тем до некоторых пор это лучше. Нужно бы знать как быстро "плывет" сеть, чтобы выбрать фильтр, способный это отслеживать, если это надо.
Go to the top of the page
 
+Quote Post
TigerSHARC
сообщение Oct 28 2009, 17:24
Сообщение #7


Знающий
****

Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195



Цитата(sup-sup @ Oct 27 2009, 20:28) *
На пальцах так.
Предположим (для упрощения) что достаточно увеличить число точек в 10 раз.
Делаем интерполяцию любым доступным инструментом (Matlab, Excel, ..).
Для этого добавляем по 9 нулей между каждым отсчетом и "получаем частоту дискретизации" 12000 Hz.
Пропускаем этот "поток" через цифровой фильтр нижних частот (интерполяционный фильтр), задача которого - обязательно пропускать первую гармонику (50 Hz), но достаточно хорошо () подавлять все частоты выше 600 Гц. Можно взять для пробы затухание 40-60 дБ, а там видно будет. Нужно брать FIR, так как он не трогает фазы.
Выходной сигнал будет иметь в 10 раз больше точек, что позволит повысить "точность".
Если есть шумы, а они есть, то применяемый фильтр их подавит (вне первой гармоники), поэтому чем длиннее фильтр, тем до некоторых пор это лучше. Нужно бы знать как быстро "плывет" сеть, чтобы выбрать фильтр, способный это отслеживать, если это надо.

... не совсем понимаю, что даёт интерполяция сама по себе. Пусть даже интерполировал не нулями, а вычисленными значениями (пусть интерполяция сплайнами).
Получил 24 выборки (анализ в течении 0.02с) -> Интерполировал, получил 240 точек на период -> провёл fft...
получил спектр, стало видно, что на погрешность это не влияет(((...
как стало видно увеличение числа точек даёт увеличение полосы спектра, но спектральные отчёты продолжают итди с тем же интервалом...

другое дело - это когда увеличиваем интервал анализа( скажем берём не один "период" (0,02с), а пусть 20 (0.4)), то тут спектральные составляющие идут уже с частотой 2.5 Гц (1/T), а не 50Гц.
Но проблема заключается в том, что нужно получать данные о БПФ каждые 20мс(стало быть по данным одного периода) и получается что априори имеем дело с 24 выборками((...

Сообщение отредактировал TigerSHARC - Oct 28 2009, 17:41
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- TigerSHARC   минимизация погрешности   Sep 4 2009, 10:48
- - fontp   Цитата(TigerSHARC @ Sep 4 2009, 14:48) Ес...   Sep 4 2009, 10:55
|- - Den   Цитата(fontp @ Sep 4 2009, 14:55) Посмотр...   Nov 20 2012, 12:13
|- - fontp   QUOTE (Den @ Nov 20 2012, 15:13) Добрый д...   Nov 20 2012, 13:17
- - thermit   ЦитатаЕсть сигнал. Частота основной гармоники 50 г...   Sep 4 2009, 11:12
|- - fontp   Цитата(thermit @ Sep 4 2009, 15:12) Похож...   Sep 4 2009, 11:16
- - thermit   ЦитатаБудет размазано. Но если априорно известно, ...   Sep 4 2009, 11:27
|- - sup-sup   Интереснео, но не очень понятно. Похоже на передач...   Sep 4 2009, 18:20
- - TigerSHARC   Это нужно для оценки самой сети (в трёхфазной сети...   Sep 5 2009, 12:16
|- - sup-sup   Так что, частота дискретизации задана от кварца, а...   Sep 5 2009, 17:09
|- - alex_os   Цитата(TigerSHARC @ Sep 5 2009, 16:16) Эт...   Sep 5 2009, 17:33
- - bahurin   Цитата(TigerSHARC @ Sep 4 2009, 14:48) Ес...   Sep 6 2009, 13:01
- - thermit   ЦитатаЭто нужно для оценки самой сети (в трёхфазно...   Sep 7 2009, 06:56
- - TigerSHARC   Цитата(fontp @ Sep 4 2009, 14:55) Посмотр...   Oct 15 2009, 19:31
|- - fontp   Цитата(TigerSHARC @ Oct 15 2009, 23:31) ....   Oct 16 2009, 07:20
- - bahurin   Цитата(TigerSHARC @ Sep 4 2009, 14:48) Ес...   Oct 16 2009, 11:17
- - Alex11   Мы тут уже сделали прибор для измерения парметров ...   Oct 16 2009, 18:04
- - AndeyP   Если нужно что то померить и точности FFT не хвата...   Oct 16 2009, 19:04
|- - bahurin   Цитата(AndeyP @ Oct 16 2009, 23:04) Если ...   Oct 17 2009, 05:54
|- - AndeyP   Цитата(bahurin @ Oct 17 2009, 09:54) Доба...   Oct 17 2009, 09:48
|- - fontp   Цитата(AndeyP @ Oct 17 2009, 13:48) Гейзе...   Oct 17 2009, 10:08
|- - AndeyP   Цитата(fontp @ Oct 17 2009, 14:08) Здесь ...   Oct 17 2009, 19:16
|- - GetSmart   Цитата(AndeyP @ Oct 18 2009, 01:16) А есл...   Oct 18 2009, 17:20
||- - TigerSHARC   Цитата(GetSmart @ Oct 18 2009, 21:20) Тип...   Oct 18 2009, 18:49
|- - bahurin   Цитата(AndeyP @ Oct 17 2009, 23:16) А спр...   Oct 18 2009, 19:28
|- - GetSmart   Цитата(bahurin @ Oct 19 2009, 01:28) прин...   Oct 18 2009, 22:08
||- - fontp   Цитата(GetSmart @ Oct 19 2009, 02:08) При...   Oct 19 2009, 07:34
||- - анатолий   Цитата(fontp @ Oct 19 2009, 10:34) В реал...   Oct 19 2009, 10:21
||- - bahurin   Цитата(анатолий @ Oct 19 2009, 14:21) Зде...   Oct 19 2009, 10:37
||- - fontp   Цитата(анатолий @ Oct 19 2009, 14:21) Есл...   Oct 19 2009, 10:43
|- - AndeyP   Цитата(bahurin @ Oct 18 2009, 23:28) прин...   Oct 19 2009, 20:29
|- - fontp   Цитата(AndeyP @ Oct 20 2009, 00:29) Мне к...   Oct 20 2009, 09:07
|- - AndeyP   Выкладываю свою тестовую среду для алгоритмов оцен...   Oct 20 2009, 20:49
|- - GetSmart   Цитата(AndeyP @ Oct 21 2009, 02:49) Может...   Oct 20 2009, 20:58
- - анатолий   Лет 15 назад мы такое делали на 8051 микроконтролл...   Oct 17 2009, 15:36
- - TigerSHARC   Амплитуда мерялась по максимуму интерполяции графи...   Oct 18 2009, 15:15
- - sup-sup   Цитата(TigerSHARC @ Sep 4 2009, 13:48) Ес...   Oct 19 2009, 05:10
|- - bahurin   Цитата(sup-sup @ Oct 19 2009, 09:10)...   Oct 19 2009, 05:56
|- - sup-sup   Цитата(sup-sup @ Oct 19 2009, 08:10)...   Oct 19 2009, 11:04
- - thermit   ЦитатаМогучий принцип. Наверное им хорошо обосновы...   Oct 20 2009, 07:18
- - анатолий   Странные какие-то споры. Сигнал 50 Гц - это, в общ...   Oct 20 2009, 10:50
|- - Oldring   Цитата(TigerSHARC @ Oct 21 2009, 23:38) ....   Oct 21 2009, 20:39
|||- - bahurin   Цитата(TigerSHARC @ Oct 28 2009, 20:24) ....   Oct 28 2009, 17:48
|||- - sup-sup   ... не совсем понимаю, что даёт интерполяция сама ...   Oct 28 2009, 19:04
|||- - TigerSHARC   Цитата(sup-sup @ Oct 28 2009, 23:04)...   Oct 28 2009, 20:36
|||- - sup-sup   Цитата(TigerSHARC @ Oct 28 2009, 23:36) ....   Oct 28 2009, 21:25
|||- - bahurin   Цитата(sup-sup @ Oct 29 2009, 00:25)...   Oct 29 2009, 05:12
|||- - sup-sup   Цитата(bahurin @ Oct 29 2009, 09:12) Вы т...   Oct 29 2009, 06:14
|||- - TigerSHARC   Цитата(sup-sup @ Oct 29 2009, 10:14)...   Oct 29 2009, 18:23
||- - alexkok   Цитата(TigerSHARC @ Oct 27 2009, 18:10) ....   Oct 27 2009, 16:30
|- - sup-sup   Цитата(TigerSHARC @ Oct 21 2009, 22:38) ....   Oct 22 2009, 04:12
|- - bahurin   Цитата(TigerSHARC @ Oct 21 2009, 23:38) ....   Oct 22 2009, 05:11
- - TigerSHARC   ... как ещё можно минимизировать погрешность???   Oct 29 2009, 20:21
|- - fontp   Цитата(TigerSHARC @ Oct 29 2009, 23:21) ....   Oct 30 2009, 06:35
||- - sup-sup   Цитата(fontp @ Oct 30 2009, 10:35) Вообще...   Oct 30 2009, 11:51
|- - AlexU   Цитата(TigerSHARC @ Oct 29 2009, 23:21) ....   Nov 1 2009, 17:52
|- - GetSmart   Цитата(AlexU @ Nov 1 2009, 23:52) На вход...   Nov 1 2009, 18:33
||- - AlexU   Цитата(GetSmart @ Nov 1 2009, 21:33) Это ...   Nov 2 2009, 06:29
||- - fontp   Цитата(GetSmart @ Nov 1 2009, 21:33) Это ...   Nov 2 2009, 09:08
||- - AlexU   Цитата(fontp @ Nov 2 2009, 12:08) "У...   Nov 2 2009, 09:27
|- - GetSmart   Цитата(AlexU @ Nov 1 2009, 23:52) На вход...   Nov 2 2009, 18:49
|- - AlexU   Цитата(GetSmart @ Nov 2 2009, 21:49) Отку...   Nov 3 2009, 08:20
- - Евгений Николаев   TigerSHARC, Вы определитесь - Вам вычисление часто...   Nov 1 2009, 20:12
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 1 2009, 23...   Nov 3 2009, 16:40
- - Евгений Николаев   Если аппаратно, то речь, как я понимаю идёт о ПЛИС...   Nov 5 2009, 08:04
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 5 2009, 11...   Nov 5 2009, 12:27
- - Евгений Николаев   Цитатакоэффицент несинусоидальности посчитать... У...   Nov 5 2009, 15:09
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 5 2009, 18...   Nov 5 2009, 15:30
- - Евгений Николаев   Цитатадля предотвращения наложения используется RC...   Nov 5 2009, 16:16
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 5 2009, 19...   Nov 5 2009, 19:38
- - Евгений Николаев   Все эти приборы работают с окном наблюдения намног...   Nov 6 2009, 13:17
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 6 2009, 16...   Nov 6 2009, 14:00
|- - анатолий   Цитата(Евгений Николаев @ Nov 6 2009, 16...   Nov 7 2009, 15:23
|- - TigerSHARC   Цитата(анатолий @ Nov 7 2009, 18:23) Прав...   Dec 1 2009, 13:56
- - Евгений Николаев   А чем Вас так держит аппаратная частота? Можно мен...   Nov 6 2009, 14:54
|- - TigerSHARC   Цитата(Евгений Николаев @ Nov 6 2009, 17...   Nov 6 2009, 16:07
- - Den   Спасибо. Уважаемый fontp, прочитал много Ваших по...   Nov 21 2012, 05:13
|- - TigerSHARC   Цитата(Den @ Nov 21 2012, 09:13) Спасибо....   Nov 21 2012, 07:40
|- - fontp   QUOTE (Den @ Nov 21 2012, 08:13) Спасибо....   Nov 21 2012, 09:12
|- - Den   Цитата(fontp @ Nov 21 2012, 13:12) Амплит...   Nov 21 2012, 09:33
||- - fontp   QUOTE (Den @ Nov 21 2012, 12:33) А разве ...   Nov 21 2012, 09:36
|- - Den   Цитата(fontp @ Nov 21 2012, 13:12) Я прив...   Nov 21 2012, 09:47
|- - fontp   QUOTE (Den @ Nov 21 2012, 12:47) Большое ...   Nov 21 2012, 09:59
- - Den   128 точечное БПФ у меня высчитывается за 1.1 мкс. ...   Nov 21 2012, 08:32
- - Den   Спасибо. Буду пробывать. С уважением, Den.   Nov 21 2012, 10:09
- - Serg76   Уважаемые знатоки! Для обсуждаемых здесь алгор...   Nov 21 2012, 18:46
|- - fontp   QUOTE (Serg76 @ Nov 21 2012, 21:46) Уважа...   Nov 22 2012, 08:21
- - Serg76   2 fontp Спасибо большое за такой развернутый отве...   Nov 22 2012, 11:14
- - Serg76   2 fontp пока тема не "остыла", хотел ещ...   Nov 22 2012, 17:31
|- - fontp   QUOTE (Serg76 @ Nov 22 2012, 20:31) 2 fon...   Nov 23 2012, 08:31
|- - Serg76   Цитата(fontp @ Nov 23 2012, 11:31) Смысл ...   Nov 23 2012, 11:02
- - Den   Добрый день, уважаемый fontp! Что-то я совсем...   Dec 1 2012, 07:27
|- - fontp   QUOTE (Den @ Dec 1 2012, 10:27) Добрый де...   Dec 2 2012, 11:35
- - Den   Добрый день. Пытаюсь всё разобраться. Вот ещё код...   Dec 20 2012, 10:03


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

 


RSS Текстовая версия Сейчас: 21st June 2025 - 06:27
Рейтинг@Mail.ru


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