"I'm trying to add an 'if' statement to analyze each result, but Matlab is 'ignoring' the condition and not performing the addition."

1 次查看(过去 30 天)
"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 个)

the cyclist
the cyclist 2023-10-19
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 个评论
Walter Roberson
Walter Roberson 2023-10-19
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.
Daniel santiago Umaña
"The values of c are attached."
I need to take these values, and if they are less than 0, d = rad2deg(theta) + 180, if they are greater than zero, d = rad2deg(theta)

请先登录,再进行评论。


Torsten
Torsten 2023-10-19
移动:Torsten 2023-10-19
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 个评论
Torsten
Torsten 2023-10-20
编辑:Torsten 2023-10-20
We don't know the size of the array "theta", but most probably you mean
for i = 1:size(c,1)
if c(i) < 0
d(i) = rad2deg(theta(i)) + 180
else
d(i) = rad2deg(theta(i))
end
end
or shorter
d = rad2deg(theta);
d(c<0) = d(c<0) + 180;

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by