插值计算

# -*-coding:utf-8 -*-
# -*-soc-OCV曲线插值计算 -*-
import numpy as np
from scipy import interpolate
import pylab as pl
import pandas as pd
import xlrd
import xlwings as xw
Filename = input('输入文件名,且在程序同一文件夹下: >> ')
df=pd.read_excel(Filename+'.xlsx',sheet_name=0)
x=df['soc']
y=df['ocv']
print(df)
xnew=np.linspace(0,1,100000)
pl.plot(x,y,"ro")
for kind in ["cubic"]:#插值方式
#"nearest","zero"为阶梯插值
#slinear 线性插值
#"quadratic","cubic" 为2阶、3阶B样条曲线插值
f=interpolate.interp1d(x,y,kind=kind)
# ‘slinear’, ‘quadratic’ and ‘cubic’ refer to a spline interpolation of first, second or third order)
ynew=f(xnew)
pl.plot(xnew,ynew,label=str(kind))
pl.legend(loc="lower right")
pl.savefig('SOC-OCV')
df = pd.DataFrame(ynew,xnew)
print(df)
df.to_excel('SOC-OCV.xlsx')
pl.show()


上一篇:洛谷P2801 分块


下一篇:【OSATE学习笔记】AADL Property 特性