"I'm trying to add an 'if' statement to analyze each result, but Matlab is 'ignoring' the condition and not performing the addition."
1 view (last 30 days)
Show older comments
"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 Comments
Walter Roberson
on 19 Oct 2023
Is theta a scalar or a vector? Is it calculated from c? Is it possibly the same as c ?
Answers (2)
the cyclist
on 19 Oct 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 Comments
Walter Roberson
on 19 Oct 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
on 19 Oct 2023
Moved: Torsten
on 19 Oct 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
See Also
Categories
Find more on Numerical Integration and Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!