🗊Презентация Програмна реалізація фрактальної моделі броунівського руху

Модуль inceigen.m function [M,Rxx,Sxx,G,rho,iFault]=inceigen(H,N,G,maxG); %Если значение G не было задано, то оно задается так, чтобы % 2AG>=2(N-1) if G==0 rho=log(2*(N-l))/log(2); GTemp=round(rho); if Gtemp<rho GTemp=GTemp+1; end; else GTemp=G; end; %Проверка GTemp перед началом цикла if GTemp>maxG іFault=2; return; end; iFault=0; rho=l; eig=0; sum=0; М=2^GTemp; MHalf=M/2; Rxx=zeros(1,M); for i=l:MHalf Rxx(i)=inccov(і—1,H); end; Rxx(MHalf+1:end)=Rxx(MHalf:-1:1); %Вычисление ДПФ Sxx=fft(Rxx);

Робота програми Запуск програми здійснюється шляхом набору у командному рядку MATLAB імені програми. Файл з текстом програми повинен знаходитись у поточному каталозі системи MATLAB. Вхідними параметрами програми є: показник Херста (Н [0,1]) та розмір реалізації процесу, що моделюється. Розглянуті 3 випадки моделювання для різних значень показника Херсту: Н<0.5 (кореляція негативна), Н=0.5 (кореляція відсутня) та Н>0.5 (кореляція позитивна).  Рисунок 1 - вигляд коваріаційної функції та спектральної щільності процесу Рисунок 2 містить реалізацію фрактальної броунівської функції та фрактального шуму. Рисунки 3 (додаток Б) ілюструє властивості масштабної інваріантності узагальненого броунівського руху:а) відображує процес, де елементарний крок часу дорівнює 4t, б) містить реалізацію вихідного процесу. Рисунок 4 (додаток Б) - приведений процес приростів фрактальної броунівської функції для реєстрації на кожному кроці t і 4t відповідно

Випадок 1: Н<0.5 Моделювання фрактального броунівського руху Введіть показник Херста Н: 0.1 Введіть необхідну кількість відліків процесу N:2400 Встановлене значення G: 13. Рисунок 1 а) Коваріаційна функція б) Спектральна щільність

Випадок 2: Н=0.5 » main Моделювання фрактального броунівського руху Введіть показник Херста Н: 0.5 Введіть необхідну кількість відліків процесу N:10000 Встановлене значення G: 15. Рисунок 1 - Характеристики процесу: а) Коваріаційна функція б) Спектральна щільність

Випадок 3: Н>0.5 » main Моделювання фрактального броунівського руху Введіть показник Херста Н: 0.9 Введіть необхідну кількість відліків процесу N:4000 Встановлене значення G: 13. Рисунок 1 - Характеристики процесу: а) Коваріаційна функція б) Спектральна щільність