Как вывести угол на время в Мупаде

Question

Так что у меня есть неприятная функция с синусами и косинусами, которая представляет положение некоторой точки в определенной системе. Теперь, когда я знаю местоположение точки, зависит от угла бета. Я хочу получить функцию, чтобы найти скорость. Проблема в том, что mupad считает, что бета-версия является константой, когда вы пытаетесь вывести ее на время. Очевидно, что производная бета - это угловая скорость. Но как мне сказать это Мупаду?

Это код, который я до сих пор.

reset();
eq:=(a/cos(Beta))^2=(a/cos(Alpha))^2+d^2-2*a/cos(Alpha)*d*sin(Alpha);

Ex:=-a+Lb*cos(Beta);

a:=simplify(solve(eq,a,Real,IgnoreAnalyticConstraints) assuming PI/2<Alpha<PI and 0<Beta<PI/2)[2][1]

Ex:=((combine(Ex,IgnoreAnalyticConstraints)) assuming Beta in R_ and 0<Beta<PI/2)

plot(Ex|[d=0.122,Alpha=11.5*PI/180,Lb=0.8],Beta=0..PI/2,Scaling=Constrained);

Ey:=a*tan(Alpha)+Lb*sin(Beta)

plot(Ey|[d=0.122,Alpha=11.5*PI/180,Lb=0.8],Beta=0..PI/2,Scaling=Constrained)

d:=0.09:
Alpha:=11.5*PI/180:
Lb:=0.85:

plot1:=plot::Curve2d([Ex,Ey],Beta=0..PI/2):
plot2:=plot::Point2d([Ex,Ey],Beta=0..PI/2,Frames=200):
plot(plot1,plot2,x=0.65,ViewingBoxXRange=0..0.7, ViewingBoxYRange=0..0.9, Scaling=Constrained);
delete(d);delete(Alpha);delete(Lb);


pointE:=matrix([Ex,Ey])

diff(pointE,t) assuming diff(Beta,t)=Omega

Я пробовал что-то в последней строке, но это явно не сработало.

Answer 1

Я нашел решение этой проблемы:

V_E:=diff(puntE|[Beta=Beta(t)],t)