Constraints other than nonnegativity on ode solver?
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
Antworten (1)
Matt J
am 13 Dez. 2014
Bearbeitet: Matt J
am 15 Dez. 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 Kommentare
Matt J
am 15 Dez. 2014
Bearbeitet: Matt J
am 15 Dez. 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.
Siehe auch
Kategorien
Mehr zu Operating on Diagonal Matrices finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!