收起
这样改写一下:
import matplotlib.pyplot as plt
import numpy as np
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
m=1575;Cd=0.3;A=2.5;f=0.012;at=0.9;g=9.8;dt=1.1
t=[]
for i in range(1,51):
Pm=49+i
u=np.arange(1e-20,100)
Ft=3600*Pm*at/u
F1=m*g*f+Cd*A*u**2/21.15
a=(Ft-F1)/(dt*m)
b=1/a
avg_b=b[1:100]+np.diff(b)/2
deta_t=np.diff(u)*avg_b
t.append(sum(deta_t)/3.6)
Pm=np.arange(50,100)
plt.plot(Pm,t)
plt.xlabel('功率/kW')
plt.ylabel('加速时间/s')
plt.show()
报告相同问题?