2017-03-02 1 views
0

그래서 특정 시스템에서 특정 지점의 위치를 ​​나타내는 사인과 코사인을 사용하는 pritty 불쾌한 기능이 있습니다. 이제 각도 베타에 의존하는 지점의 위치를 ​​알았습니다. 나는 속도를 찾기위한 함수를 도출하고자한다. 문제는 mupad가 베타가 일정 시간 동안 파생하려고 시도 할 때 상수라고 생각한다는 것입니다. 우연히 베타의 유도체는 각속도입니다. 하지만 어떻게 이것을 mupad에게 말합니까?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 

마지막 줄에서 뭔가를 시도했지만 분명히 작동하지 않았습니다.

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

답변

0

나는이 문제에 대한 해결책을 발견
관련 문제