用途:给定折断点坐标信息,插值出折线图上对应点的y值
例如:
输入: cut_dots=np.array([[0,0.5],[8,0.65],[18,0.996],[63,0.996]]) #need min/max
输出: 0.5 0.519 0.538 0.556 0.575 0.594 0.612 0.631 0.65 0.685 0.719 0.754 0.788 0.823 0.858 0.892 0.927 0.961 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996 0.996
完整代码
import numpy as np
import matplotlib.pyplot as plt
def write2txt(txt_path,write_array,type):
with open(txt_path,'w') as f:
for i in range(write_array.size):
if type=='float':
print(np.around(write_array[i],3),end=' ',file=f)
print(np.around(write_array[i],3),end=' ')
elif type=='int':
print(write_array[i].astype(int),end=' ',file=f)
print(write_array[i].astype(int),end=' ')
def line_mode(cut_dots):
cut_dots_T = cut_dots.T
x = cut_dots_T[0]
y = cut_dots_T[1]
table_x = np.linspace(x[0],x[-1],x[-1]-x[0]+1)
for i in range(x.size-1):
if i==0:
table_y = np.linspace(y[0],y[1],x[1]-x[0]+1)
else:
table_y = np.concatenate((table_y, np.linspace(y[i],y[i+1],x[i+1]-x[i]+1)[1:]))
plt.plot(table_x, table_y)
plt.show()
return(table_x,table_y)
if __name__ == '__main__':
# 1. 定义折点,须包括第一个和最后一个点(务必确保最小和最大的x值),定义输出值类型
cut_dots=np.array([[0,0.5],[8,0.65],[18,0.996],[63,0.996]]) #need min/max
type = 'float'#'int'
# 2. 绘制多段折线图,并返回拟合的y值
table_x,table_y = line_mode(cut_dots)
# 3. 输出&输出为txt,指定输出文件名
txt = r'table_y.txt'
write2txt(txt,table_y,type)