|
Не могу понять, откуда лезут гармоники. |
|
|
|
Jan 13 2010, 16:55
|

Участник

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

|
Вопрос в следующем: Берем отрезок синусоиды и делаем его преобразование Фурье. Получаем ожидаемую картинку - "палка" на соответствующей частоте. Далее дискретизируем сигнал и применяем окно Кайзера, делаем преобразование Фурье, получаем набор, кроме основной, набор "палок" малой амплитуды.  При этом сигнал такой  Подскажите, откуда они появляются? Я, как понимаю, шум от квантования сигнала должен быть белвм.
|
|
|
|
|
Jan 14 2010, 06:30
|
Частый гость
 
Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531

|
Цитата(Panukov_ag @ Jan 13 2010, 19:55)  должен быть белым. Это не утверждение, а рекомендация для грубой оценки. Обычно рядом приводят картинку с пилообразной формой ошибки квантования. Ничего белого. Для медленно изменяющихся сигналов ошибка выглядит как пилообразный компонент сигнала, частота которого определяется скоростью нарастания исходного сигнала. Приведенные Вами картинки типичны для любого реального прибора и обязательно приведены в документации каждого АЦП. Пр этом у сигма-дельта АЦП имеются дополнительные шумы около гармоник исходного сигнала.
--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
|
|
|
|
|
Jan 14 2010, 13:52
|
Знающий
   
Группа: Участник
Сообщений: 674
Регистрация: 26-08-05
Пользователь №: 7 997

|
Цитата(Panukov_ag @ Jan 14 2010, 13:00)  По оси y децибелы, Вы правы. Понятно, что уровень гармоник мал, но в данном случае вопрос чисто теоретический, поэтому данные гармоники важны. Вероятно, Вы взяли сигнал, частота которого строго совпадает с сеткой FFT. Это идеализирует сигнал под бесконечный. Надо немного сдвинуть частоту сигнала и спектр повеселеет. Сравнивать с кайзером (или чем то еще) нужно прямоугольное окно. То есть, нужны нулевые хвосты сигнала справа и слева.
|
|
|
|
|
Jan 14 2010, 14:24
|

Участник

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

|
Цитата(SPACUM @ Jan 14 2010, 10:30)  Это не утверждение, а рекомендация для грубой оценки. Обычно рядом приводят картинку с пилообразной формой ошибки квантования. Ничего белого. Для медленно изменяющихся сигналов ошибка выглядит как пилообразный компонент сигнала, частота которого определяется скоростью нарастания исходного сигнала. Был бы благодарен, если бы Вы подсказали сабж. Хотелось бы посмотреть на этот пилообразный компонент. Цитата(sup-sup) Вероятно, Вы взяли сигнал, частота которого строго совпадает с сеткой FFT. Напротив, старался от этого избавиться. Число периодов на длительность сигнала не целое.
|
|
|
|
|
Jan 14 2010, 15:15
|
Частый гость
 
Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531

|
Цитата(Panukov_ag @ Jan 14 2010, 17:24)  Был бы благодарен, если бы Вы подсказали сабж. Хотелось бы посмотреть на этот пилообразный компонент. В Вашем случае все правильно, пила получается при малой разрядности и медленном изменении сигнала. Ваш спектр и есть спектр белого шума. Я специально усиливал шумы резистора и смотрел спектр - получилось очень похоже. И у реального прибора шумовой спектр выглядит так же. Только при значительном усреднении множества реализаций можно получить гладкую кривую.
--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
|
|
|
|
|
Jan 14 2010, 15:26
|

Участник

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

|
Цитата(SPACUM @ Jan 14 2010, 18:15)  В Вашем случае все правильно, пила получается при малой разрядности и медленном изменении сигнала. Ваш спектр и есть спектр белого шума. Я специально усиливал шумы резистора и смотрел спектр - получилось очень похоже. И у реального прибора шумовой спектр выглядит так же. Только при значительном усреднении множества реализаций можно получить гладкую кривую. Да, действительно, понизил частоту и получил пилу на уровне в -300. Осталось только строго понять математически...
|
|
|
|
|
Jan 14 2010, 15:32
|
Частый гость
 
Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531

|
Попробуйте сильно ограничить разрядность 24 бит - 16 бит - 8 бит. Во всех случаях картина будет одинаковая только уровень шума больше. И спектры шума будут иметь приблизительно одинаковый вид. Цитата(Panukov_ag @ Jan 14 2010, 18:26)  Осталось только строго понять математически... Это такая "грязная" математика с вероятностями и функциями рпспределения. Наверное есть публикации, но Я не читал.
--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
|
|
|
|
|
Jan 14 2010, 23:46
|

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

|
Подскажите, откуда они появляются? Я, как понимаю, шум от квантования сигнала должен быть белвм. [/quote]
Найдите ссылку на форуме (сам качал вроде из этого раздела) на книгу "Quantization Noise". Глава 12 по шумам fp и 15 для понимания того, что происходит при fft. Для double - SNR=10*log10((2^(53*2)/2)/(1/12)), т.е. 326дБ только для представления сигнала и без учета выполнения других операций.
|
|
|
|
|
Jan 15 2010, 03:26
|

Участник

Группа: Участник
Сообщений: 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
|
|
|
|
|
Jan 15 2010, 09:08
|

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

|
Эта.
|
|
|
|
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|