Моделирование фазового шума во временной области, Как сгенерировать синус с фазовым шумом в MatLab/Octave/SciLab??? Опции

[:-)]

Стало интересно как можно сгенерировать неущее колебание с наложенным фазовым шумом для моделирования система ФАПЧ. Google/yandex неохотно подсказывают как это сделать.

Вот один из найденных вариантов: http://www.dsprelated.com/showcode/246.php

Если я правильно понял идею, то сначала генерируется белый шум, от него берется БПФ, потом на полученный спектр накладывается требуемая форма спектра фазовых шумов и в конце берется ОБПФ.

Неясно как это использовать, если хочется дополнительно наложить меняющийся во времени Доплер на несущую. Есть ли ещё какие способы моделирования фазового шума во временной области?