MATLAB Answers

For Loop Help Needed

1 view (last 30 days)
Kevin Egan
Kevin Egan on 19 Mar 2020
clear all, close all, clc
figpath = '../figures/';
addpath('./utils');
%% generate Data
polyorder = 5;
usesine = 0;
sigma = 10; % Lorenz's parameters (chaotic)
beta = 8/3;
rho = 28;
n = 3;
x0=[-8; 7; 27];
range2 = [((10^3.6)*.001), ((10^3.7)*.001) , ((10^3.8)*.001) , ((10^3.9)*.001) , ((10^4)*.001)]
for i = range2
tspan=[.001:.001:i];
N = length(tspan);
options = odeset('RelTol',1e-12,'AbsTol',1e-12*ones(1,n));
[t,x]=ode45(@(t,x) lorenz(t,x,sigma,beta,rho),tspan,x0,options);
for j=1:length(x)
dx(j,:) = lorenz(0,x(j,:),sigma,beta,rho);
end
dx_n0 = (dx);
Theta_n0 = poolData(x,n,polyorder,usesine);
m = size(Theta_n0,2);
total_n0(i,:) = [dx_n0 Theta_n0];
end
I currently have the above code but am receiving the error "Index in position 1 is invalid. Array indices must be positive integers or logical values.". Is it possible to run this code through the range of values and save each resulting dataframe separately?

  4 Comments

Show 1 older comment
Kevin Egan
Kevin Egan on 19 Mar 2020
function dy = lorenz(t,y,sigma,beta,rho)
dy = [
sigma*(y(2)-y(1));
y(1)*(rho-y(3))-y(2);
y(1)*y(2)-beta*y(3);
];
From running the code without total_n0(i,:) = [dx_n0 Theta_n0]; there doesn't seem to be an error. It appears that the code runs through the entire range and saves the final (10^4)*0.001 from tspan
Ameer Hamza
Ameer Hamza on 19 Mar 2020
I also don't have function poolData, but from your description in comment, the error seems to be related to dimensions of dx_n0 and Theta_n0.
Kevin Egan
Kevin Egan on 19 Mar 2020
Yes, if I could possibly align the results side by side, that would be ideal. Although the length of the matrix would be unequal.
The resulting total_n0 would initially have 59 columns, therefore I'd like to add the next matrix to the 60th column.

Sign in to comment.

Answers (1)

Harsha Priya Daggubati
Harsha Priya Daggubati on 23 Mar 2020
It would be more helpful, if you could copy the entire error message you get while posting the question in the forum. I guess the issue might be with the value of 'i'. You can make use of breakpoints to debug the issue much easily.
If that's no the case other workaround would be to possibly define 'total_n0' to be a cell array, to store variable length arrays.
total_n0{i} = [dx_n0 Theta_n0];
Hope this helps!

  0 Comments

Sign in to comment.

Tags