基于物理的小型程序中的无限循环

Question

这是一个模拟网球从 50 米高的建筑物一侧抛出的程序。 程序应输出每个时间步的 x、y 和速度值。 但是，我似乎陷入了无限循环。

 #include<stdio.h>
 #include<math.h>

 int main() {

     //Intial values 
     float ax = 0; //acceleration in the horizontal direction
     float ay = -9.8; //acceleration in the downward direction
     float x = 0; //top of building at position 0 
     float y = 50; //building is height 50 m
     float vx = 10*cos(30); //velocity in the horizontal direction = 10 m/s * cos(30); 
     float vy = 10*sin(30); //velocity in the vertical direction = 10 m/s * sin(30);     
     int time = 0; //time starts at 0 seconds
     float deltaTime = 0.001; //increment time by .001 each iteration

     //while ball is greater than 0, or above the ground which is at position 0
     while(y > 0) {

     time = time + deltaTime;
     vx = vx + ax*deltaTime;
     vy = vy + ay*deltaTime;
     x = x + vx*deltaTime + (1/2*ax*deltaTime*deltaTime);
     y = y + vy*deltaTime + (1/2*ay*deltaTime*deltaTime);     

     printf("x = %f, y = %f, vx = %f, vy = %f, time = %d, ", x,y,vx,vy,time);

     }
     system ("PAUSE"); 
     return 0;

 }

我的猜测是 y 永远不会变得小于 0，但由于我有限的物理知识，我不知道如何解决它。

Answer 1

1/2 == 0 而不是 0.5

由于 1 和 2 都是整数，因此使用整数除法截断为最接近的整数。使用 0.5f 获取 float 或仅使用 0.5 获取双精度值。

Answer 2

time 是一个 int，而不是 float；那么它不会永远保持为零吗？

Answer 3

将时间声明为 float，而不是 int：因此，它在当前代码中根本不会改变。
中的三角函数接受弧度，而不是度数。

Answer 4

看一下这个词

y = y + vy*deltaTime + (1/2*ay*deltaTime*deltaTime);   

Answer 5

如果您有兴趣了解您的代码如何在 C# 中进行很少的修改而无需杀手循环即可完成：

            float ax = 0; //acceleration in the horizontal direction
            float ay = -9.8f; //acceleration in the downward direction
            float x = 0; //top of building at position 0 
            float y = 50; //building is height 50 m
            float vx = 10f * (float)Math.Cos(30); //velocity in the horizontal direction = 10 m/s * cos(30); 
            float vy = 10 * (float)Math.Sin(30); //velocity in the vertical direction = 10 m/s * sin(30);     
            float time = 0; //time starts at 0 seconds
            float deltaTime = 0.001f; //increment time by .001 each iteration

            //while ball is greater than 0, or above the ground which is at position 0
            while (y > 0)
            {

                time = time + deltaTime;
                vx = vx + ax * deltaTime;
                vy = vy + ay * deltaTime;
                x = x + vx * deltaTime + (1 / 2 * ax * deltaTime * deltaTime);
                y = y + vy * deltaTime + (1 / 2 * ay * deltaTime * deltaTime);

                Console.WriteLine("x = {0}, y = {1}, vx = {2}, vy = {3}, time = {4}, ", x, y, vx, vy, time);

            }
            Console.ReadKey();

我所做的唯一修改是对 Cos 和 Sin 的强制转换以及更改浮动时间。
并在一些初始值（如转换）后添加“f”。

这可能不是 C 的真正答案，但它可能是一条线索？