如何在Python 2.6中将简谐振荡器的时间与位置数据保存到TXT文件

你的代码存在几个关键问题导致无法保存数据：

• 无限循环while 1==1永远不会终止，后续保存逻辑根本执行不到
• 就算执行到，w模式打开文件会覆盖之前内容，无法累积数据
• 运动学迭代逻辑错误，计算出的位置和时间数据不准确

以下是修正后的代码，实现时间t和位置x两列数据的正确保存：

from visual.graph import *

m = 10
k = 10
x = 2.0  # 用浮点数避免整数除法问题
v = 0.0
t = 0.0
dt = 0.1

# 使用with语句管理文件，自动处理关闭
with open("arq02.txt", "w") as f:
    # 写入表头（可选，方便后续查看）
    f.write("时间t,位置x\n")
    # 循环到t达到20时停止
    while t < 20:
        # 计算加速度
        a = -(k / m) * x
        # 欧拉-Cromer法更新速度和位置（比普通欧拉法更稳定）
        v += a * dt
        x += v * dt
        t += dt
        # 将t和x格式化为字符串，按列写入，每行一条数据
        # Python2.6不支持f-string，用传统格式化语法
        f.write("%0.2f,%0.4f\n" % (t, x))

核心改动说明： #

• 把循环条件改为while t < 20，确保循环在时间达标后终止
• 用with语句操作文件，无需手动调用close()，避免资源泄漏
• 移除了无用的x1和v1变量，修正运动学迭代逻辑，采用更适合简谐运动的欧拉-Cromer法
• 写入时添加表头，并用逗号分隔两列数据，生成的文件可直接用Excel、表格工具打开

如果想要生成制表符分隔的TSV文件（兼容性更好），只需修改写入的分隔符：

f.write("%0.2f\t%0.4f\n" % (t, x))

若数据量较大，可先将所有数据存入列表，最后一次性写入文件（效率更高）：

from visual.graph import *

m = 10
k = 10
x = 2.0
v = 0.0
t = 0.0
dt = 0.1

data = []
# 先收集所有数据
while t < 20:
    a = -(k / m) * x
    v += a * dt
    x += v * dt
    t += dt
    data.append((t, x))

# 一次性写入文件
with open("arq02.txt", "w") as f:
    f.write("时间t\t位置x\n")
    for t_val, x_val in data:
        f.write("%0.2f\t%0.4f\n" % (t_val, x_val))

内容的提问来源于stack exchange，提问作者Edinardo Rodrigues