import numpy as npimport matplotlib.pyplot as pltimport 数据a_time = 数据.a_timea_flow = 数据.a_flowc_time = 数据.c_timec_flow = 数据.c_flowdef lex(x,flow,time):#三项插值(抛物插值) lnx = flow[i] * (x - time[i + 1]) * (x - time[i + 2]) / ( (time[i] - time[i + 1]) * (time[i] - time[i + 2])) + \ flow[i + 1] * (x - time[i]) * (x - time[i + 2]) / ( (time[i + 1] - time[i]) * (time[i + 1] - time[i + 2])) + \ flow[i + 2] * (x - time[i]) * (x - time[i + 1]) / ( (time[i + 2] - time[i]) * (time[i + 2] - time[i + 1])) return lnxx1 = list()y1 = list() #a插值数据x2 = list()y2 = list() #c插值数据 #a流量for i in range(8): for x in np.arange(a_time[i], a_time[i + 2], 0.01): lx = lex(x,a_flow,a_time) y1.append(lx) x1.append(x)for i in range(len(x1)): #a插值平均 x1[i] = float('%.2f'%x1[i])for j in range(len(x1)): for k in range(j+1,len(x1)): if x1[j] == x1[k]: y1[j] = (y1[j] + y1[k]) / 2 y1[k]
不是很难,主要是你贴的太乱了,你点击插入代码段,然后复制,来个有格式的,而且你的是有一个 import 数据???,是有个文件叫 数据.py 吗?为什么要用中文名字?里面又是什么内容?
太复杂了,看不出来
您好,我是有问必答小助手,您的问题已经有小伙伴解答了,您看下是否解决,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632