2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Wolfram Mathematica
Сообщение03.11.2014, 19:10 
Задача такова.
Рассм. плоскость, наклоненную к горизонту под углом альфа. У основания плоскости стоит пушка, которая стреляет под углом бета к горизонту, выпуская заряд со скоростью V_0. Найти расстояние, на котором снаряд попадет в наклонную плоскость, учитывая силу трения, пропорциональную скорости. Построить траекторию снаряда для трех различных значений коэффициента трения.
К построению траектории:
[math]$Clear[[math][math]$ s_1 = DSolve [{my''[t] + mg - ky'[t] == 0, y'[0] == v_0\cos[b], y[0] == 0}, y[t], t] $[math][math]$ y[t_] = y[t] /. s_1$[math][math]$ s_2 = DSolve [{mx''[t] + mg - kx'[t] == 0, x'[0] == v_0\sin[b], x[0] == 0}, x[t], t] $[math][math]$x[t_] = x[t] /. s_2$[math][math]$g = 9.8;$[math][math]$v_0 = 30;$[math][math]$b = \pi/3;$[math][math]$k = 1;$[math][math]$m = 10;$[math][math]$x[t_] = x[t] /. s_2$[math][math]$y[t_] = y[t] /. s_1$[math][math]$f_1 = ParametricPlot[{x[t], y[t]}, {t, 5, 10}]$[math]
пишу вот так, дифуры считает, а график не строится, при этом ошибки не выдаются. Появляются только оси координат.

 
 
 
 Re: Wolfram Mathematica
Сообщение03.11.2014, 22:12 
newfag,
из вашего кода нельзя понять, почему не получается, код просто не читаем, но сразу видно, что уравнения не правильные: 1) по $x$ действует $ m g$, 2) знак у трения не правильный. Вот решение, разбирайтесь:
Код:
{X[t_],Y[t_]}=DSolveValue[
   {
   m x''[t]+k x'[t]==0,
   m y''[t]+m g+k y'[t]==0,
   x[0]==0,
   y[0]==0,
   x'[0]==v Cos[\[Alpha]],
   y'[0]==v Sin[\[Alpha]]
   },
   {x[t],y[t]},
   t
]  ;

Off[Solve::ifun] ;
T = t  /. Solve[Tan[\[Beta]] X[t] == Y[t],t] // First  // Simplify ;

PARS = N[Thread[{m,g,v,\[Alpha],k,\[Beta]}-> {10,9.8,30,Pi/3,1.,Pi/6}]];
Show[
   Plot[
      Evaluate[{Tan[\[Alpha]]x,Tan[\[Beta]]x} /. PARS],
      {x,0,Evaluate[1.1 X[T]/. PARS ]},
      Frame -> True,
      PlotStyle->{{Black,Dashed}},
      PlotLabel -> PARS
   ],
   ParametricPlot[
   Evaluate[{X[t],Y[t]} /. PARS],
   {t,0,Evaluate[T /.PARS]},
   PlotStyle->{Thick,Red}
   ],
   Graphics[{PointSize[Large],Black,Point[Evaluate[{X[T],Y[T]} /. PARS]]}],
   ImageSize -> 600
]

 
 
 
 Re: Wolfram Mathematica
Сообщение12.11.2014, 18:25 
Спасибо большое, с этой темой разобрался:)

 
 
 [ Сообщений: 3 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group