> restart; with(DЕtools):
> DEplot({diff(x(t),t)=y, diff(y(t),t)=x-x^3},
[x(t),y(t)], t=0..20, [[0,1,0.2], [0,0,1],
[0,1,0.4], [0,1,0.75], [0,0,1.5], [0,-0.1,0.7]],
stepsize=0.1, arrows=none, linecolor=black);
2. Мынадай дифференциалдық теңдеулер жүйесінің фазалық портретін төменде берілген бірнеше бастапқы шарттар үшін құру керек болсын:
,
> restart; with(DЕtools):
> sys:=diff(x(t),t)=y, diff(y(t),t)=sin(x):
> DEplot({sys},[x(t),y(t)], t=0..4*Pi, [[0,1,0],
[0,-1,0], [0,Pi,1], [0,-Pi,1], [0,3*Pi,0.2],
[0,3*Pi,1], [0,3*Pi,1.8], [0,-2*Pi,1]],
stepsize=0.1, linecolor=black);
3. Мынадай дифференциалдық теңдеулер жүйесінің фазалық портретін салу керек:
Бастапқы шарттар, айнымалының өзгеру диапазоны және координаталық осьтердің өлшемдері фазалық портреттің көркемдігіне байланысты алынатын болсын.
> restart; with(DЕtools):
> sys:=diff(x(t),t)=3*x+y, diff(y(t),t)=-x+y:
> phaseportrait([sys],[x(t),y(t)],t=-10..10,
[[0,1,-2], [0,-3,-3], [0,-2,4], [0,5,5], [0,5,-3],
[0,-5,2], [0,5,2], [0,-1,2]], x=-30..30,y=-20..20,
stepsize=.1, colour=blue,linecolor=black);
Достарыңызбен бөлісу: |