# ode45 goes busy endlessly when using the sign of a state variable

조회 수: 2(최근 30일)
Diaa 2021년 4월 23일
답변: Diaa 2021년 4월 23일
This term "- fr*Mass*9.8*sign(y(2))" makes ode45 busy forever specifically when using 'sign(y(2))'
I just need to make -fr*Mass*9.8=0 only when y(2)=0
Any explanations are highly welcomed.
clearvars
clc
close all
[t,y] = myexercise;
plot(t,y(:,2))
function [t,y] = myexercise
tspan = [0 , 10];
y0 = [ 0 ; 5 ];
Mass = 1300;
CD = 0.3;
fr = 0.015;
Af = 1.6;
traction = @(t) (t>3&t<5)*5*t;
[t,y] = ode45(@(t,y) deceleration(t,y) , tspan , y0);
function yp = deceleration(t,y)
yp = [ ...
y(2) ;...
( -0.5*1.23*CD*Af*y(2)^2 + traction(t) - fr*Mass*9.8*sign(y(2)) )/Mass ...
];
end
end
##### 댓글 수: 2표시숨기기 이전 댓글 수: 1
Diaa 2021년 4월 23일
It is a moving vehicle subjected to traction force and resistive aerodynamic drag and rolling resistance. Your comment made me realize that I forgot dividing the whole second equation by Mass. After fixing it, fr*Mass*9.8)*(y(2)>0) was enough to make it work.

댓글을 달려면 로그인하십시오.

### 채택된 답변

Diaa 2021년 4월 23일
Thanks to the explanation of Star Strider, the use of sign() can be replaced by fr*Mass*9.8)*(y(2)>0) to make it work.

댓글을 달려면 로그인하십시오.

### 추가 답변(1개)

Walter Roberson 2021년 4월 23일
You really should read ballode to see how to use event functions. https://www.mathworks.com/help/matlab/math/ode-event-location.html

댓글을 달려면 로그인하십시오.

### Community Treasure Hunt

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

Start Hunting!