I have solved an initial value problem and I have gotten the following equation for that:
theta=(c_0/c_1)- (2/c_1)*atan( exp(-a*c_1*t)*tan((c_0-c_1*theta_0)/2) )
c_0=7*pi/6; c_1=0.3; a=0.055; theta_0=0;
I do expect the MATLAB returns theta=0 for t=0. In other words what I expect to see is:
because for t=0, the first equation can be simplified and as a result we have: theta=theta_0 : independent of c_0,c_1(~=0),and a.
but MATLAB returns something else:
I would be grateful if somebody could explain me how I can get what I expect to get?
thanks a lot, Vahid
No products are associated with this question.
All inverse trigonometry functions return to a specific limited range, because trig functions are periodic. Hence, if x = 9*pi/2, then sin(x) will be 1, so asin(sin(x)) will be pi/2, not 9*pi/2. That's what's happening here -- atan returns values between -pi/2 and pi/2 (see doc atan):
(c_0-c_1*theta_0)/2 % ans = 1.8326 > pi/2 tan((c_0-c_1*theta_0)/2) atan(tan((c_0-c_1*theta_0)/2)) atan(tan((c_0-c_1*theta_0)/2)) + pi