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

 
 
2 страниц V   1 2 >  
Reply to this topicStart new topic
> Не могу понять, откуда лезут гармоники.
Panukov_ag
сообщение Jan 13 2010, 16:55
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 31
Регистрация: 1-02-08
Из: Омск
Пользователь №: 34 669



Вопрос в следующем:
Берем отрезок синусоиды и делаем его преобразование Фурье. Получаем ожидаемую картинку - "палка" на соответствующей частоте.
Далее дискретизируем сигнал и применяем окно Кайзера, делаем преобразование Фурье, получаем набор, кроме основной, набор "палок" малой амплитуды.


При этом сигнал такой



Подскажите, откуда они появляются? Я, как понимаю, шум от квантования сигнала должен быть белвм.
Go to the top of the page
 
+Quote Post
eugen_pcad_ru
сообщение Jan 14 2010, 06:12
Сообщение #2


Знающий
****

Группа: Свой
Сообщений: 642
Регистрация: 15-11-07
Пользователь №: 32 353



Вы перепутали шумы по оси абсцисс и по оси ординатsmile.gif При дискретизации конечного буфера шумы дискретизации (не квантования, при квантовании все как раз нормально) имеют явно периодический характер - не так? Собственно Вы их и наблюдаете.


--------------------
Правильно сформулированый вопрос содержит в себе половину ответа.
P.S.: Некоторые модераторы в качестве ответа так навязчиво предлагают посетить свой сайт, что иначе как саморекламу такие действия интерпретировать сложно.
Go to the top of the page
 
+Quote Post
SPACUM
сообщение Jan 14 2010, 06:30
Сообщение #3


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

Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531



Цитата(Panukov_ag @ Jan 13 2010, 19:55) *
должен быть белым.

Это не утверждение, а рекомендация для грубой оценки. Обычно рядом приводят картинку с
пилообразной формой ошибки квантования. Ничего белого.
Для медленно изменяющихся сигналов ошибка выглядит как пилообразный компонент сигнала,
частота которого определяется скоростью нарастания исходного сигнала.
Приведенные Вами картинки типичны для любого реального прибора и
обязательно приведены в документации каждого АЦП.
Пр этом у сигма-дельта АЦП имеются дополнительные шумы около гармоник исходного сигнала.


--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
Go to the top of the page
 
+Quote Post
bahurin
сообщение Jan 14 2010, 07:18
Сообщение #4


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



если по оси y децибелы на спектре то я бы даже не парился. А если это линейный масштаб, то тогда вопрос сколько уровней квантования?
Go to the top of the page
 
+Quote Post
Panukov_ag
сообщение Jan 14 2010, 09:00
Сообщение #5


Участник
*

Группа: Участник
Сообщений: 31
Регистрация: 1-02-08
Из: Омск
Пользователь №: 34 669



По оси y децибелы, Вы правы. Понятно, что уровень гармоник мал, но в данном случае вопрос чисто теоретический, поэтому данные гармоники важны.
Go to the top of the page
 
+Quote Post
sup-sup
сообщение Jan 14 2010, 13:52
Сообщение #6


Знающий
****

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



Цитата(Panukov_ag @ Jan 14 2010, 13:00) *
По оси y децибелы, Вы правы. Понятно, что уровень гармоник мал, но в данном случае вопрос чисто теоретический, поэтому данные гармоники важны.

Вероятно, Вы взяли сигнал, частота которого строго совпадает с сеткой FFT.
Это идеализирует сигнал под бесконечный.
Надо немного сдвинуть частоту сигнала и спектр повеселеет.
Сравнивать с кайзером (или чем то еще) нужно прямоугольное окно.
То есть, нужны нулевые хвосты сигнала справа и слева.
Go to the top of the page
 
+Quote Post
Panukov_ag
сообщение Jan 14 2010, 14:24
Сообщение #7


Участник
*

Группа: Участник
Сообщений: 31
Регистрация: 1-02-08
Из: Омск
Пользователь №: 34 669



Цитата(SPACUM @ Jan 14 2010, 10:30) *
Это не утверждение, а рекомендация для грубой оценки. Обычно рядом приводят картинку с
пилообразной формой ошибки квантования. Ничего белого.
Для медленно изменяющихся сигналов ошибка выглядит как пилообразный компонент сигнала,
частота которого определяется скоростью нарастания исходного сигнала.


Был бы благодарен, если бы Вы подсказали сабж. Хотелось бы посмотреть на этот пилообразный компонент.

Цитата(sup-sup)
Вероятно, Вы взяли сигнал, частота которого строго совпадает с сеткой FFT.


Напротив, старался от этого избавиться. Число периодов на длительность сигнала не целое.
Go to the top of the page
 
+Quote Post
bahurin
сообщение Jan 14 2010, 14:54
Сообщение #8


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



ваши гармоники на уровне -270 дБ это очень мало. Причин может быть несколько:
1. шумы квантования сигнала
2. шумы как результат округлений при вычислении БПФ
3. шумы при вычислении комплексных экспоент при рассчете БПФ (синусы и косинусы считаются с ошибкой возникают фазовые рассогласования)
4. микроскопическая нелинейность при рассчете исходного сиганала (микроскопическое отличие полученной синусоиды от истинной вываливает кучу гармоник которые между собой комбинируются и привет).
5. возможно ошибки при рассчете окна какие нить вылезли. Окно кайзера не выражается через элементарные функции а это опять ошибки рассчета
6. любая комбинация или все перечисленные варианты.

А вообще на уровне -270 дБ уже квантовая механика рулит там может все что угодно быть, даже то чего в принципе быть не может. Так что думаю лучше успокоится.
Go to the top of the page
 
+Quote Post
SPACUM
сообщение Jan 14 2010, 15:15
Сообщение #9


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

Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531



Цитата(Panukov_ag @ Jan 14 2010, 17:24) *
Был бы благодарен, если бы Вы подсказали сабж. Хотелось бы посмотреть на этот пилообразный компонент.

В Вашем случае все правильно, пила получается при малой разрядности и медленном изменении сигнала.
Ваш спектр и есть спектр белого шума. Я специально усиливал шумы резистора и смотрел спектр - получилось очень похоже.
И у реального прибора шумовой спектр выглядит так же. Только при значительном усреднении множества реализаций можно получить гладкую кривую.


--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
Go to the top of the page
 
+Quote Post
Panukov_ag
сообщение Jan 14 2010, 15:26
Сообщение #10


Участник
*

Группа: Участник
Сообщений: 31
Регистрация: 1-02-08
Из: Омск
Пользователь №: 34 669



Цитата(SPACUM @ Jan 14 2010, 18:15) *
В Вашем случае все правильно, пила получается при малой разрядности и медленном изменении сигнала.
Ваш спектр и есть спектр белого шума. Я специально усиливал шумы резистора и смотрел спектр - получилось очень похоже.
И у реального прибора шумовой спектр выглядит так же. Только при значительном усреднении множества реализаций можно получить гладкую кривую.


Да, действительно, понизил частоту и получил пилу на уровне в -300. Осталось только строго понять математически...
Go to the top of the page
 
+Quote Post
SPACUM
сообщение Jan 14 2010, 15:32
Сообщение #11


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

Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531



Попробуйте сильно ограничить разрядность 24 бит - 16 бит - 8 бит. Во всех случаях картина будет одинаковая только уровень шума больше. И спектры шума будут иметь приблизительно одинаковый вид.

Цитата(Panukov_ag @ Jan 14 2010, 18:26) *
Осталось только строго понять математически...

Это такая "грязная" математика с вероятностями и функциями рпспределения. Наверное есть публикации, но Я не читал.


--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
Go to the top of the page
 
+Quote Post
anton
сообщение Jan 14 2010, 21:37
Сообщение #12


Знающий
****

Группа: Свой
Сообщений: 622
Регистрация: 31-07-04
Пользователь №: 422



матлаб использует 64бит плавучку но и она конечна 270дб это 45 чесных бит.
Go to the top of the page
 
+Quote Post
sdy
сообщение Jan 14 2010, 23:46
Сообщение #13





Группа: Участник
Сообщений: 10
Регистрация: 17-12-04
Из: Санкт-Петербург
Пользователь №: 1 541



Подскажите, откуда они появляются? Я, как понимаю, шум от квантования сигнала должен быть белвм.
[/quote]

Найдите ссылку на форуме (сам качал вроде из этого раздела) на книгу "Quantization Noise". Глава 12 по шумам fp и 15 для понимания того, что происходит при fft. Для double -
SNR=10*log10((2^(53*2)/2)/(1/12)), т.е. 326дБ только для представления сигнала и без учета выполнения других операций.
Go to the top of the page
 
+Quote Post
Panukov_ag
сообщение Jan 15 2010, 03:26
Сообщение #14


Участник
*

Группа: Участник
Сообщений: 31
Регистрация: 1-02-08
Из: Омск
Пользователь №: 34 669



Цитата
Найдите ссылку на форуме (сам качал вроде из этого раздела) на книгу "Quantization Noise". Глава 12 по шумам fp и 15 для понимания того, что происходит при fft.


Вы имеете в виду это?
Название: Quantization noise
Автор: Widrow B. Kollar I.
Издательство: Cambridge university press

Сабж карйне интересен, спасибо за наводку.

Сообщение отредактировал Panukov_ag - Jan 15 2010, 03:48
Go to the top of the page
 
+Quote Post
sdy
сообщение Jan 15 2010, 09:08
Сообщение #15





Группа: Участник
Сообщений: 10
Регистрация: 17-12-04
Из: Санкт-Петербург
Пользователь №: 1 541



Эта.
Go to the top of the page
 
+Quote Post

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

 


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


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