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

 
 
 
Reply to this topicStart new topic
> Моделирование шума, по заданной спектр. плотности в полосе
sergunas
сообщение Feb 19 2008, 12:56
Сообщение #1


Местный
***

Группа: Свой
Сообщений: 441
Регистрация: 7-12-04
Пользователь №: 1 373



дана спектр. плотность w в полосе F
Как смоделировать подачу такого шума во времени на вход модели?
Go to the top of the page
 
+Quote Post
fontp
сообщение Feb 19 2008, 12:59
Сообщение #2


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



Цитата(sergunas @ Feb 19 2008, 15:56) *
дана спектр. плотность w в полосе F
Как смоделировать подачу такого шума во времени на вход модели?


Берётся белый шум и подаётся на фильтр h(F)
w(F) = h(F)*h'(F)
Go to the top of the page
 
+Quote Post
sergunas
сообщение Feb 19 2008, 13:27
Сообщение #3


Местный
***

Группа: Свой
Сообщений: 441
Регистрация: 7-12-04
Пользователь №: 1 373



Цитата(fontp @ Feb 19 2008, 15:59) *
Берётся белый шум и подаётся на фильтр h(F)
w(F) = h(F)*h'(F)

1) А белый шум какого размаха по амплитуде берётся?
2) Я так понимаю h(F) - идеальный полосовой фильтр, а
поясните что такое h'(F) и формулу я не понял w(F) = h(F)*h'(F).
Go to the top of the page
 
+Quote Post
fontp
сообщение Feb 19 2008, 13:31
Сообщение #4


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



Цитата(sergunas @ Feb 19 2008, 16:27) *
1) А белый шум какого размаха по амплитуде берётся?
2) Я так понимаю h(F) - идеальный полосовой фильтр, а
поясните что такое h'(F) и формулу я не понял w(F) = h(F)*h'(F).


1) обычно генераторы генерируют с единичной дисперсией, множитель нужно посчитать,
что нибудь типа 1/fd
2) h' комплексно-сопряжение от h
Go to the top of the page
 
+Quote Post
sergunas
сообщение Feb 19 2008, 13:59
Сообщение #5


Местный
***

Группа: Свой
Сообщений: 441
Регистрация: 7-12-04
Пользователь №: 1 373



Цитата(fontp @ Feb 19 2008, 16:31) *
1) обычно генераторы генерируют с единичной дисперсией, множитель нужно посчитать,
что нибудь типа 1/fd
2) h' комплексно-сопряжение от h


Простите, вероятно туплю, поясните поподробнее:
1) вот мне задана спектральная плотность шума w (Вт/Гц), и мне нужно сгенерировать белый шум такой амплитуды, чтобы его спектр был для всех частот на уровне w, так?

2) h' комплексно-сопряжение - понял, а как понять формулу w(F) = h(F)*h'(F)? То есть нужно подобрать h(F) такой, чтобы h(F)*h'(F) равнялось заданной w?
Go to the top of the page
 
+Quote Post
fontp
сообщение Feb 22 2008, 07:41
Сообщение #6


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



Цитата(sergunas @ Feb 19 2008, 16:59) *
Простите, вероятно туплю, поясните поподробнее:
1) вот мне задана спектральная плотность шума w (Вт/Гц), и мне нужно сгенерировать белый шум такой амплитуды, чтобы его спектр был для всех частот на уровне w, так?


ну да. Коэффициент пропорциональности нужно считать c учётом полосы и усиления

Цитата(sergunas @ Feb 19 2008, 16:59) *
2) h' комплексно-сопряжение - понял, а как понять формулу w(F) = h(F)*h'(F)? То есть нужно подобрать h(F) такой, чтобы h(F)*h'(F) равнялось заданной w?


да w = |h|*|h| в части зависимости от F. w - спектральная плотность
Go to the top of the page
 
+Quote Post
SpellDrive
сообщение Feb 29 2008, 11:18
Сообщение #7





Группа: Новичок
Сообщений: 9
Регистрация: 2-07-07
Пользователь №: 28 829



w(F) в данном примере - коэффициент передачи по мощности фильтра ( |h|*|h| = h*h', где h - коэффициент передачи по амплитуде).

А теперь по порядку.
1) Необходимо сгенерировать белый гауссовский шум. Если это Matlab, то это просто randn. Если это физическая модель - можно сделать на регистрах сначала равномерное распределение, а затем проссумировать несколько отсчетов равномерного.
2) Необходимо полученный белый гауссовский шум пропустить через полосовой фильтр - он оставит из "бесконечной" полосы ту часть, что Вам нужна.
3) Промасштабировать до нужной дисперсии. Для математики напоминаю, что sigma^2 = N0 * Fd / 2.
Go to the top of the page
 
+Quote Post
repairDV
сообщение Feb 29 2008, 13:42
Сообщение #8


Знающий
****

Группа: Свой
Сообщений: 578
Регистрация: 7-11-06
Из: Хабаровск
Пользователь №: 22 044



Это функция SinX/X. Я делал по этой теме диплом. Институт делал устройство для анализа канала связи, от меня требовалось сделать схему узла, формирующего сигнал в определённом диапазоне частот с одинаковой амплитудой. Где-то валяется, если не потерялось. Там - все формулы.


--------------------
Маленький нанайца. А-а. А-а. Оморочком плыл. Маленький проточка. Осетра ловил.
Go to the top of the page
 
+Quote Post
SpellDrive
сообщение Mar 1 2008, 10:30
Сообщение #9





Группа: Новичок
Сообщений: 9
Регистрация: 2-07-07
Пользователь №: 28 829



Цитата
Это функция SinX/X


Если Вы о импульсной характеристике идеального полосового фильтра, то почти - sinc(x) * гармонику на центральной частоте... Тут проблема в том, что когда обрезаешь отсчеты этого синка (дабы сделать ЦФ), то появляются большие выбросы... В железе лучше такой результат, полученый преобразованием Фурье П-образной АЧХ, не применять. Есть хорошая теория синтеза цифровых полосовых фильтров. ...для Матлаба же подойдет, если нет желания заморачиваться с синтезом окон или аппроксимирующих функций...
Go to the top of the page
 
+Quote Post

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

 


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


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