In a 'for loop' not getting difference of two consecutive iteration

1 次查看（过去 30 天）
feeroz babu
feeroz babu 2020-11-18
评论： feeroz babu 2020-11-19
format long
l=1/2;
u=1/3;
I=[1,0,0;0,1,0;0,0,1];
A=[4,3,3;2,4,1;3,1,2];
B=det(A)*inv(A);
G_1=[1/3,0,0;0,1/2,0;0,0,1];
G_2=[1/4,0,0;0,1/5,0;0,0,1/6];
T=[1,0,0;0,1,0;0,0,1];
F = [1/2,0,0;0,1/2,0;0,0,1/2];
i=1;
x(1) =1;
y(1) =1;
z(1) =1;
X_i=[x(i);y(i);z(i)];
for i = 1:100
a_i = 1/i^(1/2);
t = 1/5;
Z_i = X_i-t*[(I-inv(I+l*G_1))*X_i + B*((I-inv(I+u*G_2)))*A*X_i];
i = i+1;
X_i = a_i*F*X_i+(1-a_i)*T*Z_i
E(i) = norm(X_i);
L(i) = ((x(i+1) - x(i))^3 + (y(i+1) - y(i))^3 + (z(i+1) - z(i))^3)^(1/3); % main problem is here
end
n=[2:1:100];
plot(n,E(n))
  2 个评论
Geoff Hayes
Geoff Hayes 2020-11-18
feeroz - where in the above code are you trying to get the difference between two consecutive iterations? How should this difference be calculated?
feeroz babu
feeroz babu 2020-11-18
I am runing this code in GNU Octave, version 4.2.2. See we get X_i = [x(i); y(i); z(i)] in each iteration. So I want L(i). When I run this code, I got
error: x(3): out of bound 1

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Geoff Hayes
Geoff Hayes 2020-11-18
feeroz - you never update your x, y, and z after their initialization, so perhaps you are intending to use X_i instead. Try the following which may be what you want
X_i = zeros(3,101);
X_i(:,1) = [1; 1; 1];
for i = 1:100
a_i = 1/i^(1/2);
t = 1/5;
Z_i = X_i(:,i)-t*[(I-inv(I+l*G_1))*X_i(:,i) + B*((I-inv(I+u*G_2)))*A*X_i(:,i)];
X_i(:,i+1) = a_i*F*X_i(:,i)+(1-a_i)*T*Z_i;
E(i) = norm(X_i(:,i+1));
L(i) = ((X_i(1,i+1) - X_i(1,i))^3 + (X_i(2,i+1) - X_i(2,i))^3 + (X_i(3,i+1) - X_i(3,i))^3)^(1/3); % main problem is here
end
In the above, we pre-size the X_i array which we updated on each iteration of the loop. We then use that instead of the x, y and z when updating L(i). Note that you could pre-size E and L too.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Whos 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by