import matplotlib.pyplot as plt import math x = [] #创建数组 y = [] x1 =[] #plt.ion() # 开启interactive mode 成功的关键函数 l=1931 r=565/2 for i in range(int(-l/2-r),int(l/2+r+1),10): if i >= l/2: m = r**2-(i-l/2)**2 j = m**0.5 else: if i <= -l/2: m = r**2 - (i+l/2)**2 j = m**0.5 #pow(x,y) x的y次幂 else: j = r y.append(i+915.9) #将i的值逐个写入x中 x.append(round(j,2)+1279.85) #侧板的质心 已经莫名其妙的定位 x1.append(-round(j,2)+1279.85) plt.ylim() plt.xlim() #设置坐标轴的显示范围 plt.plot(x,y) fp1= open(r'F:\Users\Desktop\SCAU\JinPeng\Modeling\数据\大链轮x数据.txt','w') fp2= open(r'F:\Users\Desktop\SCAU\JinPeng\Modeling\数据\大链轮y数据.txt','w') fp3= open(r'F:\Users\Desktop\SCAU\JinPeng\Modeling\数据\大链轮-x数据.txt','w') fp1.truncate() fp2.truncate() fp3.truncate() for i in range(len(x)): fp1.write(str(x[i])+'\n') #str()数字转字符串 fp2.write(str(y[i])+'\n') fp3.write(str(x1[i])+'\n') print(len(x))
06-14
366
03-30
2589
06-01
320