MATLAB Answers

Constraints other than nonnegativity on ode solver?

13 views (last 30 days)
Iris
Iris on 13 Dec 2014
Commented: Iris on 15 Dec 2014
I'm using ode45 to solve the equations of motion of some robot arm. I know I can use the odeset option 'NonNegative' to tell the solver that certain states must be nonnegative, but what if I want certain states to be bounded by different numbers? So say, I have a state q: q = [q1 q2 q3] and I want q1 to be in between -pi/2 and pi/2 and I want q3 to be in between 0 and 1... Is there any way to do this?
Thanks in advance

  0 Comments

Sign in to comment.

Answers (1)

Matt J
Matt J on 13 Dec 2014
Edited: Matt J on 15 Dec 2014
You could try transforming the state variable like
q3=sin(y)/2+1/2
where y is the transformed state variable. The derivatives in the state equation also have to rewritten in terms of y
dq3/dt= (cos(y)/2)*dy/dt
but it looks like the cos(y)/2 part can be accounted for using a mass matrix.

  5 Comments

Show 2 older comments
Iris
Iris on 15 Dec 2014
Unfortunately, this doesn't seem to work. The point is, when z becomes zero, it should stay zero because the arm cannot drop through the floor. It should only go up again because there is a force making it go up. If I use a sine function, I do not get this result which kind of makes sense, as a sine wave is an oscillation. I really hope there is another solution to this problem.. Would be most appreciated!
Matt J
Matt J on 15 Dec 2014
Doesn't the floor exert a force on the arm? And don't your equations of motion account for that force together with any others? If so, why do you need to add separate constraints in the first place? Why don't the applied forces alone (and hence the equations of motion) enforce the constraints?
In any case, did you test the final state trajectory q3(y(t)) that my approach gives you? Does it not satisfy both the equations of motion and the constraints? If so, and yet you think there is a 2nd solution out there that you like better, it means that the solution is non-unique. In turn, that means you haven't provided enough equations to fully specify the solution.
Iris
Iris on 15 Dec 2014
Normally I would add constraints to my equations, but this is an assignment for which I have to use the port-Hamiltonian formalism and I have not included the constraints. But I was just going to try to build in constraints now myself. I'll let you know if it worked. Thanx for your input so far though!

Sign in to comment.

Sign in to answer this question.