Моделирование эффекта Доплера в КВ-канале без использования встроенных функций MATLAB. Опции

[СФТКСР]

Добрый вечер.
Собственно, в теме вопроса сам вопрос и заключается.

Реализация Doppler Shift вопросов не вызывает, т.к. вроде все понятно и есть конкретная формула, описывающая зависимость сдвига от скорости и угла:
http://picturepush.com/public/12446817

Вопрос: как, не используя встроенные функции "rayleighchan" и "ricianchan", реализовать Doppler Spread?

[andyp]

Вам нужно умножить Ваш сигнал на случайный процесс, моделирующий замирания. В качестве модели замираний используются:

1. Пропущенный через некий фильтр белый гауссовский шум. ИХ фильтра соответствует желаемой АКФ замираний. Здесь все сводится к синтезу фильтра с требуемой АЧХ.

2. Некий детерминированный процесс (обычно это сумма взвешенных синусоид) со случайными параметрами.

Про подход 2 можно почитать например здесь: Matthias Pätzold, Мobile Fading Channels

В обоих случаях, процесс, моделирующий замирания, узкополосный, поэтому, как правило, он синтезируется на низкой частоте дискретизации, а затем интерполируется.

PS В КВ обычно используются гауссовские спектры замираний (ну или бигауссовские (модель Ватерсона)). Есть также стандарт Recommendation ITU-R F.1487, описывающий каналы для тестирования модемов.
В матлабе есть какие-то модули для бигауссовской модели (здесь). Может, там и код посмотреть можно.