For Loop Help Needed

2 次查看(过去 30 天)
Kevin Egan
Kevin Egan 2020-3-19
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 个评论
Ameer Hamza
Ameer Hamza 2020-3-19
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 2020-3-19
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.

请先登录,再进行评论。

回答(1 个)

Harsha Priya Daggubati
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!

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by