Iterative solution for 6x6 system
Show older comments
Hello,
I'm trying to solve a system of equations. There are 6 equations 6 variables but one variable is dependant to two variables. Hence, I try to give an initial guess of '0' to a variable. Then I try to iteratively solve the system until the specified tolerance is reached. But when code ran I get only one value for the asked variable instead of recording all iterations.
(Constants and intermediate variables written above of this code):
for i=1:629 if s14dot(i)<0 f14(i)=200; else f14(i)=-100; end
for j=1:1:100
f34prime=f34;
A= [1 0 -mu -mu 0 0;
0 1 1 -1 0 0;
0 0 d-s14(i)-mu*(c/2) -(d+a-s14(i))+mu*(c/2) 1 0;
-a3*sin(theta13(i)) a3*cos(theta13(i)) 0 0 -1 0;
a2*sin(theta12(i)) -a2*cos(theta12(i)) 0 0 0 1;
0 0 0 0 1 0] ;
C=[-f14(i); 0; -(b+c/2)*f14(i); 0; 0; r34*mu*f34*sign(theta13dot(i))];
B=inv(A)*C;
f34x=B(1);
f34y=B(2);
g14prime=B(3);
g14dprime=B(4);
m34=B(5);
T=B(6);
f34=sqrt(f34x^2+f34y^2);
if f34prime==0
elseif (f34-f34prime)/f34prime<0.0001
break
end
end
Answers (1)
Torsten
on 31 Mar 2016
My guess is that you want to insert
F34(i)=f34
before the "end" of the loop over i:
...
if f34prime==0
elseif (f34-f34prime)/f34prime<0.0001
break
end
F34(i)=f34
end
Best wishes
Torsten.
Categories
Find more on Communications Toolbox in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!