C 练习实例20 – 小球自由下落：从物理模型到代码实现
在学习 C 语言的过程中，我们常常会遇到一些看似简单却内涵丰富的编程练习。C 练习实例20 – 小球自由下落，正是这样一个经典题目。它不仅帮助你掌握循环、浮点数运算和函数封装等核心语法，更重要的是，它让你第一次用代码“模拟真实世界”——让一个小球从空中落下，计算它的位置、速度和总路程。
这就像你在物理课上用公式推导自由落体，但这次，你不是在纸上写，而是在键盘上敲出一个“数字世界里的小球”。
问题描述与物理背景
想象一下，你从一栋高楼的顶端，轻轻松开手，一个金属小球开始自由下落。忽略空气阻力，它的运动遵循经典力学中的匀加速直线运动规律。
根据物理公式：
- 速度公式：v = g × t
- 位移公式：s = ½ × g × t²
- 其中 g 是重力加速度，通常取 9.8 m/s²
我们的目标是：编写一个 C 程序，模拟这个小球从静止开始下落，每秒记录一次它的位置和速度，直到落地为止。
这个练习看似简单，但背后涉及多个编程概念：变量定义、循环控制、浮点数精度、函数抽象和输出格式控制。我们一步步来。
设计程序结构：模块化思维
在写代码之前，先思考“做什么”和“怎么做”。这是专业程序员的基本素养。
我们可以把整个程序拆分为几个功能模块：
- 定义重力加速度常量
- 初始化小球状态（初始高度、初始速度）
- 循环计算每一秒的运动状态
- 输出结果
- 结束程序
这种“分而治之”的思想，就像搭积木：先准备好每一块，再拼成完整的模型。
创建变量与常量
#include <stdio.h>
int main() {
// 定义重力加速度，单位：m/s²
// 使用 const 保证数值不会被意外修改
const double gravity = 9.8;
// 初始高度，单位：米
double initial_height = 100.0;
// 初始速度，小球从静止释放，所以为 0
double velocity = 0.0;
// 当前高度，随时间变化
double current_height = initial_height;
// 时间变量，单位：秒
int time = 0;
// 输出表头
printf("时间(s)\t高度(m)\t速度(m/s)\n");
printf("----------------------------------------\n");
// 主循环开始
while (current_height > 0) {
// 计算当前时刻的速度：v = g * t
velocity = gravity * time;
// 计算当前高度：h = h0 - 0.5 * g * t²
current_height = initial_height - 0.5 * gravity * time * time;
// 输出当前状态
printf("%d\t\t%.2f\t\t%.2f\n", time, current_height, velocity);
// 时间递增
time++;
}
// 小球落地后，补充输出落地瞬间的状态
printf("小球在第 %d 秒落地，速度为 %.2f m/s\n", time - 1, velocity);
return 0;
}
代码注释说明：
const double gravity = 9.8;：使用const修饰重力加速度，确保它在程序运行中不会被修改，提升代码安全性和可读性。double类型用于存储浮点数，因为高度和速度可能不是整数。current_height = initial_height - 0.5 * gravity * time * time;：这是自由落体的核心公式，表示经过 t 秒后，小球下落的距离。while (current_height > 0)：循环条件，只要小球还没落地（高度大于 0），就继续计算下一秒的状态。printf("%d\t\t%.2f\t\t%.2f\n", ...);：使用格式化输出，\t表示制表符，使表格对齐；%.2f表示保留两位小数。
优化输出：让结果更清晰
原始输出虽然正确，但对齐略显混乱。我们可以通过调整格式，让表格更美观。
优化后的输出代码：
// 在主函数中替换原来的 printf 输出部分
printf("%-6s %-10s %-12s\n", "时间(s)", "高度(m)", "速度(m/s)");
printf("%-6s %-10s %-12s\n", "--------", "----------", "------------");
while (current_height > 0) {
velocity = gravity * time;
current_height = initial_height - 0.5 * gravity * time * time;
// 使用 %-6s 等左对齐，确保表格整齐
printf("%-6d %-10.2f %-12.2f\n", time, current_height, velocity);
time++;
}
printf("小球在第 %d 秒落地，速度为 %.2f m/s\n", time - 1, velocity);
输出效果示例：
时间(s) 高度(m) 速度(m/s)
-------- ---------- ------------
0 100.00 0.00
1 95.10 9.80
2 80.40 19.60
3 55.90 29.40
4 21.60 39.20
5 -22.90 49.00
小球在第 4 秒落地，速度为 39.20 m/s
提示：
%-6d表示左对齐，宽度为 6 个字符，避免右侧空格导致表格错位。
深入理解：为什么用 double 而不是 float？
在 C 语言中，float 和 double 都是浮点数类型，但精度不同：
float：单精度，约 6~7 位有效数字double：双精度，约 15~16 位有效数字
在本例中，我们使用的是 double，因为：
- 重力加速度 9.8 本身是浮点数
- 计算过程中涉及平方运算（time * time），数值可能变大
- 若使用
float，在多次迭代后可能出现精度误差，导致“小球在负高度时仍在运行”
例如，若 current_height 变成 -0.00001，程序可能还会继续循环一次，造成不准确。
所以，在涉及物理模拟、金融计算等场景时，优先使用 double。
增强功能：添加“落地判断”逻辑
当前程序中，我们通过 current_height > 0 判断是否落地，但实际落地时高度可能略小于 0（由于浮点误差）。我们可以加入一个“安全阈值”判断：
// 改为：
if (current_height <= 0.001) { // 高度接近 0，视为落地
printf("小球在第 %d 秒落地，速度为 %.2f m/s\n", time, velocity);
break; // 立即退出循环
}
这样可以避免小球“穿透地面”后继续计算。
扩展思考：从静态模拟到动态动画
虽然 C 语言不支持图形界面（除非引入第三方库如 graphics.h），但我们可以通过输出“ASCII 动画”来模拟下落过程。
例如，每秒输出一个小球的位置，用 * 表示：
// 在循环中加入：
int ground_level = 0;
int ball_position = (int)(initial_height - current_height); // 球离地面的高度（像素级）
printf("第 %d 秒：", time);
for (int i = 0; i < initial_height; i++) {
if (i == ball_position) {
printf("*"); // 小球位置
} else {
printf(" ");
}
}
printf("\n");
虽然不能播放动画，但可以让你看到小球“从高处落下”的视觉效果，非常直观。
常见问题与调试技巧
| 问题 | 原因 | 解决方法 |
|---|---|---|
| 小球下落速度一直为 0 | 未正确更新 velocity | 检查 velocity = gravity * time; 是否在循环内 |
| 程序无限循环 | 判断条件错误或浮点误差 | 增加 <= 0.001 判断，或使用 break 机制 |
| 输出高度为负数 | 计算公式错误或时间过长 | 检查 0.5 * gravity * time * time 是否正确 |
| 编译报错“undefined reference” | 缺少头文件 | 确保包含 #include <stdio.h> |
总结：从练习到思维升级
C 练习实例20 – 小球自由下落，表面上是一个简单的数学模拟，实则是对编程思维的一次全面训练。
你学会了：
- 如何用 C 语言描述物理现象
- 如何合理使用
double类型避免精度问题 - 如何通过
while循环模拟连续过程 - 如何用格式化输出构建清晰的表格
- 如何通过调试发现并修复逻辑错误
更重要的是，你开始用“程序”去理解世界——不是被动接受公式，而是主动构建模型。
当你看到一个小球在屏幕上“落地”的那一刻，你不仅在运行代码，更在验证物理定律。
这，就是编程的魅力。
下一次，不妨试试让小球反弹一次，加入弹性系数，看看它如何“弹跳”？那将是一个更有趣的挑战。