%% Impostazioni generali della simulazione n_sim = 20; % numero di simulazioni n_passi = 2000; % numero di passi per ogni simulazione dt = 0.01; % passo temporale T = dt * n_passi; % tempo finale di simulazione %% Definizione del processo: parametri e funzioni a e b %% Commentare o decommentare i diversi valori di sigmas e sigmav per ottenere diversi risultati % parametri S0 = 0; % valore di partenza del processo V0 = 1; k = 1; %sigmas = 0.01; %sigmav = 0.01; %sigmas = 0.1; %sigmav = 0.1; sigmas = 0.6; sigmav = 0.6; S = nan(n_sim, n_passi + 1); S(:,1) = S0*ones(n_sim, 1); V = nan(n_sim, n_passi + 1); V(:,1) = V0*ones(n_sim, 1); dWS = sqrt(dt) * randn(n_sim, n_passi); dWV = sqrt(dt) * randn(n_sim, n_passi); Tvec = linspace(dt, T, n_passi); for i = 1:n_passi S(:, i+1) = S(:, i) + V(:, i) * dt + sigmas * dWS(:, i); V(:, i+1) = V(:, i) - k * S(:, i) * dt + sigmav * dWV(:, i); end plot(S'); axis ([0 n_passi]);