已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
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文件

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

文号 / 917908

名片发私信
学术分 0
总主题 8 帖总回复 113 楼拥有证书:进士 机友
注册于 2020-08-28 12:52最后登录 2024-07-20 11:52
主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:广东

个人简介

火箭爱好者

文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传,正在处理中
空空如也~
处理中...
处理失败
加载中...
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
插入公式
评论控制
加载中...
文号:{{pid}}
加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}
ID: {{user.uid}}