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

4楼: Originally posted by lwloveflxgg at 2018-11-17 17:17:27
你这个不添加限制条件是不行的。你要的功能实际上实现函数 f：定义域:实数,值域:路径积分等于输入的曲线簇。不对初始速度和角度作限制的话，它就不是函数，只能解方程f(x)-x=0，而且是无穷多解。
...