Python 做曲线拟合和求积分

 这是一个由加油站油罐传感器测量的油罐高度数据和出油体积,根据体积和高度的倒数,用截面积来描述油罐形状,求出拟合曲线,再用标准数据,求积分来验证拟合曲线效果和误差的一个小项目。 主要的就是首先要安装Anaconda  python库,然后来运用这些数学工具。



###最小二乘法试验###
import numpy as np
import pymysql
from scipy.optimize import leastsq
from scipy import integrate
###绘图,看拟合效果###
import matplotlib.pyplot as plt
from sympy import *


path='E:\PythonProgram\oildata.txt'

lieh0 =[]         #初始第一列,油管高度
liev1 =[]         #初始第二列,油枪记录的体积

h_median =[]      # 高度相邻中位数
h_subtract =[]          #相邻高度作差
v_subtract =[]          #相邻体积作差
select_h_subtr =[]           #筛选的高度作差    ########
select_v_subtr =[]           #筛选的体积作差

VdtH=[]                     #筛选的V 和 t 的 倒数。

def loadData(Spath,lie0,lie1):
    with open(Spath,'r') as f0:
         for i in f0:
             tmp=i.split()
             lie0.append(float(tmp[0]))
             lie1.append(float(tmp[2]))
    print ("source lie0",len(lie0))


def connectMySQL():
    db = pymysql.connect(host='10.**.**.**', user='root', passwd='***', db="zabbix", charset="utf8") #  校罐
    cur = db.cursor()

    try:
        # 查询
        cur.execute("SELECT * FROM station_snapshot limit 10 ")
        for row in cur.fetchall():
            #   print(row)
            id = row[0]
         
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值