How to take the average of different initial conditions

1 visualización (últimos 30 días)
Samson
Samson el 28 de Ag. de 2019
Comentada: M el 10 de Sept. de 2019
I have a code and I need to run the code for different initial conditions, then save and finally take the average. See code below
a = 0.2; b = 0.2; c = 5.7; sigma = 16; beta = 4; rho = 45.92;
% initial condition
x = rand (9.1);
g = (0: 5: 60);
% g = 0.5;
% computing the trajectory
dt = 0.01;
tspan = 100,000;
xinput = x;
X = zeros (9, tspan);
deltaDR = zeros (length (g));
deltaRA = zeros (length (g));
for j = 1: length (g)
X (:, 1) = x;
for i = 1: tspan
xoutput = rk4angelstepb (@ attractor, dt, xinput, a, b, c, sigma, beta, rho, g (j));
X (:, i) = xoutput;
xinput = xoutput;
end

Respuestas (1)

M
M el 28 de Ag. de 2019
It is not clear what you want to do from the code you posted, especially the first loop :
for j = 1: length (g)
...
end
Could you give a bit more details and tell us where is the difficulty exactly ?
  2 comentarios
Samson
Samson el 28 de Ag. de 2019
The first loop is my coupling values which I want to use for the different initial conditions. I need to run my code for atleast 100 set of initial condition which I have as
x = rand (9.1)
I was thinking if I need 100 initial values, I will have it as
rand([9,1,10])
But how do I construct the loop and then average all. Also, I believe I will have to redefine
deltaDR = zeros (length (g));
deltaRA = zeros (length (g));
This is my issue
M
M el 10 de Sept. de 2019
If you want to run the code for at least 100 iterations, you should not loop over length(g), because dimension of g is 1 x 13 .
Then, rand([9,1,10]) outputs a vector of dimensions 9 x 1 x 10, so you won't have your 100 sets of initial conditions.

Iniciar sesión para comentar.

Categorías

Más información sobre Programming en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by