Как правильно использовать interp1 в matlab, чтобы найти несколько решений, как в parabola? - PullRequest
0 голосов
/ 03 марта 2019

Я решил ряд обыкновенных дифференциальных уравнений для достижения решения следующим образом:

f = @(t,x) [(-B/m)*x(1);-g-(B/m)*x(2);x(1);x(2)];
[t,x_sol] = ode45(f,[0 5],[V_0*cos(alpha) V_0*sin(alpha) 0 0]);

и составил одно из решений следующим образом:

plot(t,x_sol(:,4));

enter image description here

Все, что мне нужно сделать, это найти время, когда Y (t) станет равным нулю.Но когда я использую t_f = interp1(x_sol(:,4),t,0), я получаю первое решение из 0. Это может иметь простое решение, но я не могу найти другое решение, где Y становится равным нулю (что, по-видимому, составляет около 4 секунд из графика).На самом деле это проблема траектории, и все, что мне нужно, это просто найти время траектории.

...