Построение генератора случайного процесса методом формирующего фильтра

Цель работы –построение и проверка статистических характеристик генератора случайных чисел с заданной корреляционной функцией.

Постановка задачи

При статистическом имитационном моделировании на основе математических, полунатурных и других моделей возникает задача имитации внешних воздействий на систему, имеющих форму случайных процессов с определенными характеристиками. Эта задача решается путем построения генераторов случайных процессов.

Рассмотрим задачу имитации одномерного случайного процесса X(t). Получаемые реализации должны подчиняться закону распределения с заданной ПРВ – Построение генератора случайного процесса методом формирующего фильтра - student2.ru и иметь заданную корреляционную функцию Построение генератора случайного процесса методом формирующего фильтра - student2.ru . Генератор случайного процесса с заданными характеристиками обычно строится на основе генератора белого шума.

Метод формирующего фильтра основан на использовании закономерностей преобразования линейным динамическим звеном спектральной плотности случайного сигнала. Если на вход динамического звена поступает белый шум со спектральной плотностью Построение генератора случайного процесса методом формирующего фильтра - student2.ru , спектральная плотность выходного сигнала Построение генератора случайного процесса методом формирующего фильтра - student2.ru будет определяться через частотную передаточную функцию звена Построение генератора случайного процесса методом формирующего фильтра - student2.ru следующим образом:

Построение генератора случайного процесса методом формирующего фильтра - student2.ru (5.1)

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

Модель случайного процесса с заданной корреляционной функцией Построение генератора случайного процесса методом формирующего фильтра - student2.ru можно построить на основе формирующего фильтра, используя в качестве его входного сигнала белый шум. Необходимая передаточная функция формирующего фильтра определяется из соотношения:

Построение генератора случайного процесса методом формирующего фильтра - student2.ru (5.2)

где Построение генератора случайного процесса методом формирующего фильтра - student2.ru .

Порядок выполнения

1. Построить программный генератор случайных чисел с заданной корреляционной функцией. Использовать метод формирующего фильтра. Значения дисперсии D, параметра a и шага дискретизации h заданы в табл.5.

2. Рассчитать и построить графики заданной и полученной корреляционных функций для интервала Построение генератора случайного процесса методом формирующего фильтра - student2.ru .

3. Рассмотреть различные способы повышения точности построения процесса с заданной корреляционной функцией.

Содержание отчета

1. Графическое и аналитическое отображение заданной корреляционной функции.

2. Расчет параметров формирующего фильтра.

3. Случайный процесс, полученный в результате работы генератора.

4. Корреляционная функция, соответствующая полученному случайному процессу.

5. Текст программы на языке SciLab.

Условия моделирования

Табл.5.

№ варианта D Построение генератора случайного процесса методом формирующего фильтра - student2.ru h № варианта D Построение генератора случайного процесса методом формирующего фильтра - student2.ru h
0.001 0.001
0.1 0.05 0.001
0.5 0.002 0.1 0.2 0.01
0.1 0.001 0.001
0.01 0.1 0.01
0.001 0.001
0.2 0.01 0.002
0.001 0.001
0.01 0.01
0.5 0.02 0.001
0.0011 0.001
0.2 0.01 0.1 0.01
0.001 0.01
0.1 0.01 0.001
0.2 0.001 0.002
0.001 0.5 0.01
0.001 2.5 0.001
0.005 0.001
0.001 0.2 0.01
0.01 0.4 1.5 0.001
0.1 0.02 0.3 0.01
0.5 0.001 0.002
0.001 0.001
0.5 0.01 0.75 0.01

Приложение 1 – программа моделирования по лаб.1

clear all, close all;

global k l m n Kt i1 i2 S V ht x

global delta_max delta_max_st

% определение констант

k=1;

l=10;

m=2;

n=8;

Kt=100;

i1=10;

i2=1;

S=100;

V=800;

delta_max=0.5;

delta_max_st=2.6;

% определение начальных условий

x10=1;

x20=1;

x30=0;

x40=0;

x50=0;

X0=[x10,x20,x30,x40,x50];

ht=0.001; % шаг вычислений

Tk=10; % конечное время моделирования

time=0:ht:Tk;

[t,x]=ode45('system_diff_24',time,X0);

subplot(5,1,1),plot(t,x(:,1)),grid on;

subplot(5,1,2),plot(t,x(:,2)),grid on;

subplot(5,1,3),plot(t,x(:,3)),grid on;

subplot(5,1,4),plot(t,x(:,4)),grid on;

subplot(5,1,5),plot(t,x(:,5)),grid on;

----------------------------------------------------------------

function dx=sysmem_diff_24(t,x)

global k l m n Kt i1 i2 S V

global delta_max

teta=(10000-x(5))/(22000-V*t);

if abs(x(4))<= delta_max

x4_zv=x(4);

else

x4_zv=delta_max*sign(x(4));

end

% система дифференциальных уравнений

dx(1)=k*x(2)-k*x(1);

dx(2)=x(3);

dx(3)=l*x(1)-l*x(2)-m*x(3)+n*x(4);

dx(4)=-Kt*x4_zv-i1*x(2)-i2*x(3)+S*(teta-x(2));

dx(5)=V*sin(x(1));

dx=dx';

Наши рекомендации