在MATLAB

Question

我试图使用

找到一个情节抛物线的x截距

x0 = interp1(y,x,0)

然而，因为我的抛物线从原点开始，它返回0。

如何找到远离原点的x截距？目前我正在用眼球估计。

情节代码：

global k g

g = 10;
v0 = 150;
theta = pi/4;
m = 1;
k = 0.001;

tspan = [0 22];

IC = [0; v0*cos(theta); 0; v0*sin(theta)];
[t, oput] = ode45(@dataODE, tspan, IC);

x = oput(:,1);
vx = oput(:,2);
y = oput(:,3);
vy = oput(:,4);

figure(1); clf;
plot(x,y)

，其中

function [p] = dataODE( t, x)

global k g

p = zeros(4,1);
p(1) = x(2);
p(2) = -k*sqrt(x(2)^2 + x(4)^2)* x(2);
p(3) = x(4);
p(4) = -g -k*sqrt(x(2)^2 + x(4)^2)* x(4);

Answer 1

您可以将x和y限制为仅包含一个拦截：

x0 = interp1(y(a:b),x(a:b),0)

但您如何找到a和b？一种方法是只使用y之前和之后的点过零（在向下的路上）：

a = find(diff(y > 0) == -1)
b = a+1