Наиболее часто используется метод формирующего фильтра, т.к. позволяет получить широкий класс сигналов при манипулировании небольшим числом параметров.
Формирование выходного отсчета АРСС-модели случайного процесса происходит по следующей формуле:
,где yi – i-й отсчет выходного процесса; ai – i-й коэффициент авторегрессионного фильтра; ci – i-й коэффициент фильтра скользящего среднего; zi – i-й отсчет входного гауссовского случайного процесса; p и q – порядки фильтров авторегрессии и скользящего среднего соответственно.
Таким образом, для формирования текущего значения используются как значения входного белого шума (СС-фильтр), так и сформированные на предыдущих шагах значения выходного процесса (АР-фильтр). Следует отметить, что если СС-фильтр еще можно спроектировать для работы только с предыдущими значениями входного процесса, то АР-фильтру для правильного функционирования обязательно нужно использовать текущий отсчет входного случайного процесса.
В зависимости от конкретных задач, измерительный сигнал может быть описан как APCC-, AP- или CC-процесс с конечным числом задающих параметров. Надо отметить, что наиболее часто в технических приложениях используются AP-модели случайных процессов не старше второго порядка [3, 4], т.к. зачастую повышение порядка модели не приводит к существенному повышению точностных параметров модели и слабо влияет на процесс принятия решения при проектировании информационно-измерительных систем (ИИС).
Спектральная плотность мощности (СПМ) АРСС-процесса определяется амплитудно-частотной характеристикой формирующего фильтра. Естественно, что на СПМ выходного процесса значительное влияние окажет вид СПМ исходного гауссовского случайного процесса. Для существующих стандартных генераторов характерно случайное, негладкое поведение СПМ, поэтому при проведении исследований прибегают к усреднению либо по входному воздействию, либо по выходной реакции системы.
В том случае, когда на вход фильтра подается идеальный гауссовский белый шум, СПМ выходного процесса можно вычислить как [3]:
,где w – круговая относительная частота; s2 – дисперсия АРСС–процесса.
Учитывая все вышеизложенное и анализируя приведенную выше формулу, можно определить, что при s2 = 1 для АР- и СС-процессов АЧХ формирующих фильтров описывается следующими выражениями
– для АР-фильтра; – для СС-фильтра.Таким образом, используя метод формирующего фильтра, можно получить случайный сигнал с заданными спектральными свойствами при минимальных аппаратно-программных затратах. Однако следует учитывать и тот факт, что неидеальные параметры используемого генератора белого гауссовского шума могут существенно снизить практическую значимость метода формирующего фильтра.
Примечание: При выполнении лабораторной работы следует использовать блоки библиотек пакета Simulink «Simulink» и «DSP blockset». Так, например, в качестве элемента задержки следует использовать блок «DSP Blockset ® Signal operation ® Integer delay». Вывод гистограмм осуществляется так, как показано на рис. 3. При этом блок «MATHLAB Function» взят из библиотеки «User-Defined Function». В качестве параметра блока следует указать функцию hist(u), а размерность выходной величины установить равной нулю. Блок буфера берется из библиотеки «DSP Blockset ® Signal management ® Buffers». Размер буфера выбирается равным размеру выборки, по которой производится построение гистограммы.
Наиболее важные блоки, использующиеся при выполнении лабораторных работ, приведены в приложении.
Вариант | Параметры генератора равномерно распределенных случайных величин | Параметры гауссовского случайного процесса | Параметры АРСС-фильтра случайного процесса | |||
Разрядность слова | Разрядность регистра | Математическое ожидание | Дисперсия | Тип фильтра | Коэффициенты | |
1 | 5 | 6 | 0,1 | 1,5 | АР | а1= 0,2; а2= -0,7 |
2 | 6 | 9 | 0,5 | 2,2 | АР | а1= 0,6; а2= -0,2 |
3 | 4 | 5 | 0,2 | 1,1 | СС | с1= 0,3; с2= -0,9 |
4 | 5 | 7 | 0,3 | 4,5 | АР | а1= -0,4; а2= 0,7 |
5 | 4 | 9 | 0,0 | 2,7 | СС | с1= -0,2; с2= -0,7 |
6 | 5 | 11 | 0,4 | 1,8 | СС | с1= -0,5; с2= 0,4 |
7 | 4 | 7 | 0,5 | 3,6 | СС | с1= 0,7; с2= -0,8 |
8 | 6 | 6 | 0,3 | 0,6 | АР | а1= -0,9; а2= 0,6 |
9 | 5 | 9 | 0,4 | 2,9 | СС | с1= 0,5; с2= 0,1 |
10 | 4 | 7 | 0,1 | 4,7 | АР | а1= 0,7; а2= -0,1 |
11 | 6 | 11 | 0,0 | 2,5 | АР | а1= 0,4; а2= -0,3 |
12 | 5 | 7 | 0,3 | 1,6 | СС | с1= 0,6; с2= 0,3 |
13 | 6 | 9 | 0,5 | 0,7 | АР | а1= 0,6; а2= 0,2 |
14 | 4 | 7 | 0,8 | 3,0 | СС | с1= -0,5; с2= 0,2 |
15 | 5 | 11 | 0,2 | 2,2 | АР | а1= -0,2; а2= 0,7 |
Контрольные вопросы
1. Чем объясняется широкое применение моделей случайных сигналов как моделей информационных процессов при проектировании и отладке информационно-измерительных систем?
2. По какому принципу работают генераторы псевдослучайных последовательностей с равномерным распределением?