对了，其中图七中发动机名称和生产商为中文，是因为一开始我是这么做的，不过放openrocket里边这俩成乱码了，大概是编码问题，看不下去又懒得搞，所以就换英文了，图六是第二次换了英文之后的。另外由于Python包装后很大，一个简单的功能就33m了，所以还是推荐用源码，我以后看看能不能用c写，这样好很多。最后，如果用源码，而且之前没有安装matplotlib库的，用pip安一下吧，不过不知道为什么直接安不上，打开cmd后要输入

pip install matplotlib -i XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX/simple

可能是外国的源下不了？

顶一个，能上传.py的文件吗

可以的，不过依旧需要下载matplotlib库，以下是.py文件

解算.py 2.28KB PY 4次下载

很好的程序，但这好像有点bug,还有时间应该是横轴，力是竖轴你这是不是搞反了

应该不是的，我自己用没问题，不过这个是针对一种格式的，要不私发下文件，我看看。我最近很忙，有空我会尽量弄出图形化界面，并兼容不同格式文件。

另外这是最近的改进版本，运算速度了不少，并增加了删去试车前后的数据的功能。不过还是针对这种文件，不一样的话可以找我，我有时间改，或者闲下来后做个有图形化界面，可通用的程序。

以下是原本的F-t图

以下是第二版的F-t图

可以看出它能删去无用的数据。

多段推力的火箭也没问题，因为是只删头尾。

其余操作与1.0一样

一下是代码

以下是.py文件

解算2.0.py 2.34KB PY 5次下载

以下是.exe文件（要解压缩）

解算2.0.zip 33.58MB ZIP 2次下载

同样的，我测试过可以用，但我不知道是否会有没被发现的bug，所以有问题麻烦联系我，谢谢

提示下各位，经过试车，本程序可以用，但是我除错质量了，导致比冲计算错误，不过总冲是准的，以下是修改后的代码

import matplotlib.pyplot as plt

<br>

<br>

f = open("TEST_SD.TXT")

a = f.read()

b = a.split("\n")

c = []

time = []

force = []

I = 0

for i in b:

   if (i != "\n")and(i != ""):

       c.append(i)

for i in range(len(c)):

   d = c[i].split(",")

   time.append(float(d[0]))

   force.append(float(d[1]))

while force[0] < 3:

   force.pop(0)

   time.pop(0)#删除试车前的数据

while force[-1] < 3:

   force.pop(-1)

   time.pop(-1)#删除试车后的数据，不直接设定阈值是考虑到多推力峰

t0 = time[0]

for i in range(len(time)):

   force[i]*=0.1#我的实验得到的数据单位是0.1N,我不知道是不是都这样，如果不是的话可以改下代码

   time[i] -= t0

for i in range(len(time)-1):

   dt = time[i+1]-time[i]

   di = force[i]*dt

   I += di#积分算总冲

data = input('''请输入相关数据，格式如下：

发动机名称/直径（mm）/长度（mm）/开伞延时（s）/燃料质量（kg）/总质量（kg）/制造商

若要跳过请输入0，不计算比冲，不形成eng文件\n''')

if data != "0":

   d = data.split("/")

   name = d[0]

   diameter = int(d[1])

   length = int(d[2])

   dertatime = int(d[3])

   fuelmass = float(d[4])

   mass = float(d[5])

   maker = d[6]

else:

   fuelmass = 1

f.close()

f1 = open("时间-推力表.txt", "w")

f1.write("总冲（N.s）：%f\n"%I)

if data != "0":

   isp = (I/fuelmass)/9.8#计算以秒为单位的比冲

   f1.write("比冲（s）：%f\n"%isp)

fmax = force[0]

for i in range(len(force)):

   if fmax < force[i]:

       fmax = force[i]#计算最大推力

f1.write("最大推力（N）：%f\n"%fmax)

f1.write("时间（s）    推力（N）\n")

for i in range(len(time)):

   f1.write("%f    %f\n"%(time[i],force[i]))

f1.close()

fig = plt.figure(figsize=(20, 10), dpi=100)

plt.plot(time, force)

plt.xlabel("F(N)")

plt.ylabel("t(s)")

plt.title("F - t")

plt.show()

fig.savefig("F - t 图")#绘制并保存推力时间曲线

if data != "0":

   f2 = open("导入openrocket的数据.eng","w")

   f2.write("%s %d %d %d %f %f %s\n"%(name, diameter, length, dertatime, fuelmass, mass, maker))

   for i in range(len(time)):

       f2.write("%.3f %f\n"%(time[i], force[i]))

   f2.write(";")

   f2.close()#生成可以直接导入openrocket的.eng文件

import matplotlib.pyplot as plt


f = open("TEST_SD.TXT")
a = f.read()
b = a.split("\n")
c = []
time = []
force = []
I = 0
for i in b:
   if (i != "\n")and(i != ""):
       c.append(i)
for i in range(len(c)):
   d = c[i].split(",")
   time.append(float(d[0]))
   force.append(float(d[1]))
while force[0] < 3:
   force.pop(0)
   time.pop(0)#删除试车前的数据
while force[-1] < 3:
   force.pop(-1)
   time.pop(-1)#删除试车后的数据，不直接设定阈值是考虑到多推力峰
t0 = time[0]
for i in range(len(time)):
   force[i]*=0.1#我的实验得到的数据单位是0.1N,我不知道是不是都这样，如果不是的话可以改下代码
   time[i] -= t0
for i in range(len(time)-1):
   dt = time[i+1]-time[i]
   di = force[i]*dt
   I += di#积分算总冲
data = input('''请输入相关数据，格式如下：
发动机名称/直径（mm）/长度（mm）/开伞延时（s）/燃料质量（kg）/总质量（kg）/制造商
若要跳过请输入0，不计算比冲，不形成eng文件\n''')
if data != "0":
   d = data.split("/")
   name = d[0]
   diameter = int(d[1])
   length = int(d[2])
   dertatime = int(d[3])
   fuelmass = float(d[4])
   mass = float(d[5])
   maker = d[6]
else:
   mass = 1
f.close()
f1 = open("时间-推力表.txt", "w")
f1.write("总冲（N.s）：%f\n"%I)
if data != "0":
   isp = (I/mass)/9.8#计算以秒为单位的比冲
   f1.write("比冲（s）：%f\n"%isp)
fmax = force[0]
for i in range(len(force)):
   if fmax < force[i]:
       fmax = force[i]#计算最大推力
f1.write("最大推力（N）：%f\n"%fmax)
f1.write("时间（s）    推力（N）\n")
for i in range(len(time)):
   f1.write("%f    %f\n"%(time[i],force[i]))
f1.close()
fig = plt.figure(figsize=(20, 10), dpi=100)
plt.plot(time, force)
plt.xlabel("F(N)")
plt.ylabel("t(s)")
plt.title("F - t")
plt.show()
fig.savefig("F - t 图")#绘制并保存推力时间曲线
if data != "0":
   f2 = open("导入openrocket的数据.eng","w")
   f2.write("%s %d %d %d %f %f %s\n"%(name, diameter, length, dertatime, fuelmass, mass, maker))
   for i in range(len(time)):
       f2.write("%.3f %f\n"%(time[i], force[i]))
   f2.write(";")
   f2.close()#生成可以直接导入openrocket的.eng文件

唉，刚刚太着急，代码的模式搞错了

我的天呐，我昨天还在为这个东西着急,想着要不要一个一个去复制粘贴，没想到啊，碰到了救世主。

最近想起来了，现在已经改了一遍，不过没加图形窗口，我觉得没必要以下是.py、.exe、源码

改动内容：增加了删去错误数据的部分，改动部分代码，增加可读性，改正时间和力标反的错误，改正数据单位的错误，减小exe文件大小

解算2.0.py 2.31KB PY 0次下载

解算2.0.exe 28.77MB EXE 10次下载

(第一行注释是改名了)

#by某倒吊的亚雷斯塔（辰星引力航天工作室）
import matplotlib.pyplot as plt

f = open("TEST_SD.TXT")
c = [i for i in f.read().split('\n') if i != '']
time = [float(i.split(',')[0]) for i in c]
force = [float(i.split(',')[1]) for i in c]
I = 0
for i in force:
    if i < 0:
        i = 0#删除负数的错误数据
while force[0] < 20:
    force.pop(0)
    time.pop(0)#删除试车前的数据
while force[-1] < 20:
    force.pop(-1)
    time.pop(-1)#删除试车后的数据，不直接设定阈值是考虑到多推力峰
t0 = time[0]
for i in range(len(time)):
    force[i]*=0.098#我的实验得到的数据单位是0.098N,我不知道是不是都这样，如果不是的话可以改下代码
    time[i] -= t0
for i in range(len(time)-1):
    dt = time[i+1]-time[i]
    di = force[i]*dt
    I += di#积分算总冲
data = input('''请输入相关数据，格式如下：
发动机名称/直径（mm）/长度（mm）/开伞延时（s）/燃料质量（kg）/总质量（kg）/制造商
若要跳过请输入0，不计算比冲，不形成eng文件\n''')
if data != "0":
    d = data.split("/")
    name = d[0]
    diameter = int(d[1])
    length = int(d[2])
    dertatime = int(d[3])
    fuelmass = float(d[4])
    mass = float(d[5])
    maker = d[6]
else:
    mass = 1
f.close()
f1 = open("时间-推力表.txt", "w")
f1.write("总冲（N.s）：%f\n"%I)
if data != "0":
    isp = (I/mass)/9.8#计算以秒为单位的比冲
    f1.write("比冲（s）：%f\n"%isp)
fmax = force[0]
for i in range(len(force)):
    if fmax < force[i]:
        fmax = force[i]#计算最大推力
f1.write("最大推力（N）：%f\n"%fmax)
f1.write("时间（s）    推力（N）\n")
for i in range(len(time)):
    f1.write("%f    %f\n"%(time[i],force[i]))
f1.close()
fig = plt.figure(figsize=(20, 10), dpi=100)
plt.plot(time, force)
plt.xlabel("t(s)")
plt.ylabel("F(N)")
plt.title("F - t")
plt.show()
fig.savefig("F - t 图")#绘制并保存推力时间曲线
if data != "0":
    f2 = open("导入openrocket的数据.eng","w")
    f2.write("%s %d %d %d %f %f %s\n"%(name, diameter, length, dertatime, fuelmass, mass, maker))
    for i in range(len(time)):
        f2.write("%.3f %f\n"%(time[i], force[i]))
    f2.write(";")
    f2.close()#生成可以直接导入openrocket的.eng文件