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

 
 
13 страниц V  « < 11 12 13  
Reply to this topicStart new topic
> Измерение частоты основной гармоники (50 Гц) с точностью 0.01 Гц
Alex11
сообщение Jun 15 2017, 17:13
Сообщение #181


Гуру
******

Группа: Свой
Сообщений: 2 028
Регистрация: 23-10-04
Из: С-Петербург
Пользователь №: 965



Форма спектра будет совпадать с видом окна. Для хорошей точности нужна аппроксимация соответствующей функцией. Так что квадратичная в этом случае не подойдет. Здесь нужен косинус. Можно попробовать полином (первых три члена от разложения в ряд), но нужно меделировать и смотреть какая точность получится. И я не помню наизусть каких степеней там ряд.
С децимацией фильтр, безусловно лучше. Но еще нужно посмотреть, даст ли он выигрыш в счете по сравнению с Фурье на весь массив сразу. Последний вариант должен быть лучше по точности и усреденения и фильтров.
Go to the top of the page
 
+Quote Post
Ruslan1
сообщение Jun 15 2017, 22:01
Сообщение #182


Гуру
******

Группа: Свой
Сообщений: 2 267
Регистрация: 6-03-06
Из: Кишинев
Пользователь №: 15 025



Цитата(Alex11 @ Jun 15 2017, 19:13) *
Форма спектра будет совпадать с видом окна. Для хорошей точности нужна аппроксимация соответствующей функцией. Так что квадратичная в этом случае не подойдет. Здесь нужен косинус. Можно попробовать полином (первых три члена от разложения в ряд), но нужно моделировать и смотреть какая точность получится. И я не помню наизусть каких степеней там ряд.
С децимацией фильтр, безусловно лучше. Но еще нужно посмотреть, даст ли он выигрыш в счете по сравнению с Фурье на весь массив сразу. Последний вариант должен быть лучше по точности и усреднения и фильтров.

Спасибо! понял, поиграюсь. Благо, реальные данные уже есть в матлабе, есть материал для экспериментов.

Про FFT без децимации сразу на весь массив- сам бы так хотел. Но это у меня просто еще одна дополнительная опция в приборе, желательно бы обойтись без сотрясения основ уже сделанного. А CMSIS-овская уже используемая библиотека только FFT4096 умеет, вот и приходиться децимировать.
Мне бы вообще-то очень бы подошло что-то типа FFT длиной 16К или даже 32К точек float32, но это уже свой код применять вместо встроенных библиотек, про это тоже думал. Снижение быстродействия даже в 5 раз еще приемлемо, по сравнению с CMSIS. Тоже интересно подумать.
Go to the top of the page
 
+Quote Post
Alex11
сообщение Jun 15 2017, 22:44
Сообщение #183


Гуру
******

Группа: Свой
Сообщений: 2 028
Регистрация: 23-10-04
Из: С-Петербург
Пользователь №: 965



Попробуйте, я не смотрел, как это будет оптимизироваться на АРМ'е, но на DSP ложится прилично.
Прикрепленный файл  FFT32k.zip ( 1010 байт ) Кол-во скачиваний: 33
Go to the top of the page
 
+Quote Post
Ruslan1
сообщение Jun 16 2017, 06:40
Сообщение #184


Гуру
******

Группа: Свой
Сообщений: 2 267
Регистрация: 6-03-06
Из: Кишинев
Пользователь №: 15 025



Цитата(Alex11 @ Jun 16 2017, 00:44) *
Попробуйте, я не смотрел, как это будет оптимизироваться на АРМ'е, но на DSP ложится прилично.
Прикрепленный файл  FFT32k.zip ( 1010 байт ) Кол-во скачиваний: 33

Вы будете таки смеяться, но именно этот код я использую уже лет 5, про него и говорил. sm.gif
Работает с небольшими вариациями в разных системах, начиная от мелкого 8-битного Майкрочипа, где данные храню во внешней микросхеме сериальной RAM. Но именно в АРМе я его тоже не пробовал.
Причем, действительно, работает всегда и везде. Нашел в Интернете.
Ну раз и Вы его тоже рекомендуете, значит он действительно неплох, не буду искать что-то еще.
Спасибо!
Go to the top of the page
 
+Quote Post
ivan219
сообщение Jun 16 2017, 14:05
Сообщение #185


Местный
***

Группа: Участник
Сообщений: 345
Регистрация: 16-11-08
Пользователь №: 41 680



Интерполяцию окном Hanna смотри http://calculator2006.narod.ru/articles/haan.htm

Сообщение отредактировал ivan219 - Jun 16 2017, 14:06
Go to the top of the page
 
+Quote Post
Ruslan1
сообщение Jun 16 2017, 14:26
Сообщение #186


Гуру
******

Группа: Свой
Сообщений: 2 267
Регистрация: 6-03-06
Из: Кишинев
Пользователь №: 15 025



Цитата(ivan219 @ Jun 16 2017, 16:05) *
Интерполяцию окном Hanna смотри http://calculator2006.narod.ru/articles/haan.htm

отлично, спасибо!!!!
Скопировал себе на диск и пдф-ку сделал, а то в интернете бывает так- сегодня есть, а завтра -"404" sm.gif

То есть расчет идет по трем точкам.
И странно там в тексте как-то окно задано, формула 1-cos(2pi*n/N) правильная? Она же двойку дает в максимуме, даже с рисунком в том же тексте не совпадает.
Думаю, просто опечатка, должно быть все-таки 0.5*(1-cos(2pi*n/N)), собственно так окно в первоисточнике и задается.
Go to the top of the page
 
+Quote Post
ivan219
сообщение Jun 16 2017, 16:47
Сообщение #187


Местный
***

Группа: Участник
Сообщений: 345
Регистрация: 16-11-08
Пользователь №: 41 680



Будет спектр чуть выше.
И да скорей всего опечатка я у же и не помню почему так вышло.
Go to the top of the page
 
+Quote Post
Alex11
сообщение Jun 16 2017, 17:00
Сообщение #188


Гуру
******

Группа: Свой
Сообщений: 2 028
Регистрация: 23-10-04
Из: С-Петербург
Пользователь №: 965



Боюсь, что для Вашей жуткой точности не хватит интнрполяции по трем точкам. Лучше бы 5-7 точек или даже больше.
Go to the top of the page
 
+Quote Post

13 страниц V  « < 11 12 13
Reply to this topicStart new topic
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 21st November 2017 - 08:22
Рейтинг@Mail.ru


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