vh = 10.483314774;
v0 = 19.1184278;
w = 5;
g = 9.8;
m = 0.45;
g = 0.01337;
k = 0.5*0.5*1.29*2*3.14*0.11*0.11;
t1 = vh / (g + (k * vh)/ (2 * m));
h1 = vh * vh / (2 * (g + k * vh/(2 * m)));
t2 = sqrt(2 * (h1 -1.5) / (g - k * vh/(2 * m)))
for t=0:0.01:t1 + t2
    x = -(m * exp(log(v0) - (k * t) /m) * (k * cos((g*t*w)/m) - g * w * sin((g * t * w) /m))) / (g^2*w*w+k*k) + (m * k * exp(log(v0))) / (g*g*w*w+k*k);
    y = -(m * exp(log(v0) - (k * t) /m) * (k * sin((g*t*w)/m) + g * w * cos((g * t * w) /m))) / (g^2*w^2+k*k) + (m * exp(log(v0)) * g* w) / (g*g*w*w+k*k);
    z = -(g/exp(k * t) + k * ( g * t + vh/exp(k * t))) / ( k*k) + (g + k * vh) / ( k *k);
    plot3(y, x, z, \\\'*\\\')
    hold on
end
axis equal

3楼: Originally posted by 1014511134 at 2018-11-17 15:12:16
你这是不是相当于知道落点，速度成未知的？