Transformación estado-espacio 

Obtenga las ecuaciones gobernantes de un péndulo esférico en coordenadas cartesianas, colóquelas en la forma afín estado-espacio, y conviértalas en coordinadas esféricas.

In[1]:=

X $Line = 0;

Las ecuaciones gobernantes a partir de la función de Lagrange.

In[2]:=

X \[ScriptCapitalL] = 1/2 m (x'[t]^2 + y'[t]^2 + (x[t] x'[t] + y[t] y'[t])^2/( r^2 - x[t]^2 - y[t]^2)) + m g Sqrt[r^2 - x[t]^2 - y[t]^2];

In[3]:=

X Short[eqns = Table[D[D[\[ScriptCapitalL], q'[t]], t] - D[\[ScriptCapitalL], q[t]] == Subscript[F, q][t], {q, {x, y}}]]

Out[3]//Short=

La representación afín con estados y entradas .

In[4]:=

X Short[eqns = Table[D[D[\[ScriptCapitalL], q'[t]], t] - D[\[ScriptCapitalL], q[t]] == Subscript[F, q][t], {q, {x, y}}]]; AffineStateSpaceModel[ eqns, {x[t], y[t], x'[t], y'[t]}, {Subscript[F, x][t], Subscript[F, y][t]}, {Sqrt[r^2 - x[t]^2 - y[t]^2]}, t]; assm = AffineStateSpaceModel[%, {x, y, Subscript[x, d], Subscript[y, d]}] // Simplify

Out[4]=

La transformación entre las coordenadas esféricas y cartesianas.

In[5]:=

X ctrans = Most@ CoordinateTransform[ "Spherical" -> "Cartesian", {r, \[Theta][t], \[Phi][t]}]

Out[5]=

La transformación de los derivados de coordenadas.

In[6]:=

X ctransd = D[ctrans, t] /. {\[Theta]' -> Subscript[\[Theta], d], \[Phi]' -> Subscript[\[Phi], d]}

Out[6]=

La transformación de estado completa.

In[7]:=

X trans = Thread[{x[t], y[t], Subscript[x, d][t], Subscript[y, d][t]} -> Join[ctrans, ctransd]]

Out[7]=

Las expresiones en coordenadas esféricas son mucho más simples.

In[8]:=

X StateSpaceTransform[ assm, {trans, {\[Theta][t], Subscript[\[Theta], d][t], \[Phi][t], Subscript[\[Phi], d][t]}}] // PowerExpand

Out[8]=