# 'Index exceeds matrix dimensions' error

1 Ansicht (letzte 30 Tage)
Soho am 12 Jul. 2017
Kommentiert: Moe_2015 am 12 Jul. 2017
Hi, I'm trying to write an iteration program that runs a while loop until the variable 'theta' is almost the same as the output x (within a certain accuracy). I keep getting this error, but when I try to change the code to run, it only runs for the first loop. Can anyone help? My code is below (i've been using 0.01 as thetaInitial because i know this is pretty close to the correct answer)
function x=iteration(thetaInitial)
theta(1)=thetaInitial;
x(1)=0.0248;
error=1e-8;
Re=10000;
i=1;
while abs(x(i)-theta(i))>error
x(i)=(((6.25*log(2))+1.5)*log(Re*(theta(i)^0.5))+0.09)^-1;
i=i+1;
end
x=x(1:i+1);
plot(x)
[x(end) i]
##### 0 Kommentare-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

### Akzeptierte Antwort

Moe_2015 am 12 Jul. 2017
theta(i) is not defined after the first index. You set theta(1) but then never do anything with theta again. So when you say theta(i) it is exceeding the matrix dimension (which is 1). Also, you should move i = i+1 as the first line inside the while loop. Also outside the while loop x =x(1:i+1) will fail because you do not reach an x with length (i+1)... your x will be of length i
##### 1 KommentarKeine anzeigenKeine ausblenden
Moe_2015 am 12 Jul. 2017
So something like this (note: set theta(i) to whatever it should be):
function x=iteration(thetaInitial)
theta(1)=thetaInitial;
x(1)=0.0248;
error=1e-8;
Re=10000;
i=1;
while abs(x(i)-theta(i))>error
i=i+1;
theta(i) = 0; %set this to whatever it is supposed to be
x(i)=(((6.25*log(2))+1.5)*log(Re*(theta(i)^0.5))+0.09)^-1;
end
x=x(1:i); % this line isn't really necessary but if you want it this is what it should be
plot(x)
[x(end) i]

Melden Sie sich an, um zu kommentieren.

### Kategorien

Mehr zu Loops and Conditional Statements finden Sie in Help Center und File Exchange

### Community Treasure Hunt

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

Start Hunting!