"I'm trying to add an 'if' statement to analyze each result, but Matlab is 'ignoring' the condition and not performing the addition."
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Daniel santiago Umaña
el 19 de Oct. de 2023
Comentada: Daniel santiago Umaña
el 20 de Oct. de 2023
"I'm trying to add an 'if' statement to analyze each result, but Matlab is 'ignoring' the condition and not performing the addition."
I need that when some value of c (not all) is less than 0, d becomes rad2deg(theta) + 180, when any of the values of c is greater than or equal to 0, d becomes rad2deg(theta)
"The values of c are attached."
I think I'm not managing to explain my problem, I'll try again:
theta = acos(UV ./ Uv);
c = Y(:,2); where Y(:,2) is generated by ode45, meaning there are several values (a total of 547).
I need to check the values of "c." When the value of "c" is < 0, the value of "d" is d = rad2deg(theta) + 180; when not, d = rad2deg(theta).
So, for each value analyzed in "c," a value of "d" is calculated (if <0, d = rad2deg(theta) + 180; if not, d = rad2deg(theta)), meaning, in the end, I will have 547 values of "d" based on the rules above.
Thank you very much for your attention.
[T,Y] = ode45(@myprog,tSpan,s,options);
.
.
.
c = Y(:,2);
if c < 0
d = rad2deg(theta) + 180;
else
d = rad2deg(tetha);
end
2 comentarios
Walter Roberson
el 19 de Oct. de 2023
Is theta a scalar or a vector? Is it calculated from c? Is it possibly the same as c ?
Respuestas (2)
the cyclist
el 19 de Oct. de 2023
Is c a vector? If so, then all elements of c would have to be less than zero, to enter the if condition. This is explained in the documentation.
2 comentarios
Walter Roberson
el 19 de Oct. de 2023
The second output of ode45 is guaranteed to have at least three rows (unless there are perhaps some strange combinations of circumstances.) Just maybe if there was an event function that always signalled termination you might be able to get it down to two rows.
Torsten
el 19 de Oct. de 2023
Movida: Torsten
el 19 de Oct. de 2023
c is a vector. Do you want d to be rad2deg(theta) + 180 if all elements of the vector c are < 0 ? Then your code is correct. Otherwise you have to loop over the elements of c:
for i = 1:size(c,1)
if c(i) < 0
...
else
...
end
end
10 comentarios
Ver también
Categorías
Más información sobre Numerical Integration and Differential Equations en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!