娱乐 每 日 数 学 挑 战（伪）【第-1期】

[可洛]

这不是经典力学吗？我记得以前用mathematica专门画过这个。

（空气阻力的话一般和速度平方成正比，想再复杂一点用cv^2+kv）

[可洛]

2 小时前, Mr.K 018 说道:

是的

空气阻力的话这次先这样吧2333

hhhh这道题目让我想起了当年被经典力学支配的恐惧...

要不下一次出题考虑一下火箭？结合空气阻力和燃料消耗（真正的牛二dp/dt）的问题解起来还挺让人头疼的

补充本次题目的解（伪）：

花了一小会稍微拿matlab弄了一下，凑合着看就行。

具体数学上的解就不算了...实在是太麻烦了，已经不想再体会当年的恐惧了....

剧透

1：如图所示，k = 0.9,初速度50m/s

2：楼中楼放matlab代码，自定义一下v0, theta和k跑一下，看x_h就知道了，懒得写comment了www

剧透

clear all;
dt = 0.001;
v0 = 50;
theta = 30/180*pi();
k = 0.9;
v_h = v0*cos(theta);
v_v = v0*sin(theta);
a = k*v0;
x_h = v_h*dt;
x_v = v_v*dt;
i = 1;
while x_v > 0
    theta = atan(v_v/v_h);
    a_h = a*cos(theta);
    a_v = a*sin(theta);
    v_h = v_h - a_h*dt;
    v_v = v_v - a_v*dt - 9.81*dt;
    x_h = x_h + v_h*dt;
    x_v = x_v + v_v*dt;
    a = k*sqrt(v_h^2+v_v^2);
    a_h_rec(i) = a_h;
    a_v_rec(i) = a_v;
    v_h_rec(i) = v_h;
    v_v_rec(i) = v_v;
    x_h_rec(i) = x_h;
    x_v_rec(i) = x_v;
    i = i + 1;
end
plot(x_h_rec,x_v_rec);
axis equal;

3：这个和k相关，k越大炮台和地面倾角越小...同理，不求解了，数字解的话上面代码改一改挺好找的。

 

五月 13, 2020，由可洛修改

[可洛]

6 小时前, duya 说道:

在粘稠液体中，阻力是和v成正比的

所以只要在蜂蜜里开炮就行了（真·糖衣炮弹

无法反驳www

（其实当年学经典力学的时候给的bv的例子就是蜂蜜...）