【Python】piecewise与curve_fit分段线性函数拟合
思路:确定所有的分割点(x*,y*)以及首尾的斜率(k*)
参考:两段的分段函数;三段的分段函数;scipy.optimize.curve_fit
代码:
from scipy import optimize
import matplotlib.pyplot as plt
import numpy as npx = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9, 10 ,11, 12, 13, 14, 15], dtype=float)
y = np.array([5, 7, 9, 11, 13, 15, 28.92, 42.81, 56.7, 70.59, 84.47, 98.36, 102.25, 106.14, 110.03])# 一个输入序列,4个未知参数,2个分段函数
def piecewise_linear(x, x0, y0, k1, k2):# x<x0 ⇒ lambda x: k1*x + y0 - k1*x0# x>=x0 ⇒ lambda x: k2*x + y0 - k2*x0return np.piecewise(x, [x < x0, x >= x0], [lambda x:k1*x + y0-k1*x0, lambda x:k2*x + y0-k2*x0])def piecewise_linear3(x,x0,x1,y0,y1,k0,k1):return np.piecewise(x , [x <= x0, np.logical_and(x0<x, x<= x1),x>x1] ,[lambda x:k0*(x-x0) + y0,#根据点斜式构建函数lambda x:(x-x0)*(y1-y0)/(x1-x0)+y0,#根据两点式构建函数lambda x:k1*(x-x1) + y1])# 用已有的 (x, y) 去拟合 piecewise_linear 分段函数
p , e = optimize.curve_fit(piecewise_linear3, x, y,bounds=(0, [16,16,120,120,10,10]))xd = np.linspace(0, 15, 100)
plt.plot(x, y, "o")
plt.plot(xd, piecewise_linear3(xd, *p))
用3段的函数拟合实际2段的数据,可能出现报错
OptimizeWarning: Covariance of the parameters could not be estimated category=OptimizeWarning)
参考:限制参数取值
可以帮助解决报错,但是实际拟合效果和bound设置有关系,bound设置的比较好,拟合就会比较准
拟合比较准的用的bounds=(0, [6,13,20,120,4,6])
没有具体研究curve_fit所以原因还不明确。
【Python】piecewise与curve_fit分段线性函数拟合相关推荐
- Ilog cplex, java 表示分段线性函数 piecewise function
文章目录 1. 什么是分段线性函数 2. Ilog CPLEX 表示 Piecewise Function 分段线性函数 3. java Cplex 表示分段线性函数 4. matlab/gurobi ...
- relu函数为分段线性函数,为什么会增加非线性元素
relu函数为分段线性函数,为什么会增加非线性元素 我们知道激活函数的作用就是为了为神经网络增加非线性因素,使其可以拟合任意的函数.那么relu在大于的时候就是线性函数,如果我们的输出值一直是在大于0 ...
- Python之建模数值逼近篇–最小二乘拟合
Python之建模数值逼近篇–最小二乘拟合 介绍 系数ak的确定 函数rk(x)r_k(x)rk(x)的选取 理解和区别 样例 介绍 曲线拟合问题的提法是,已知一组(二维)数据,即平面上的n个点(x ...
- 分段线性拟合经典案例:计算多年气温最低值和最高值的分段线性变化趋势(附分段线性拟合工具下载)
分段线性回归:是用虚拟变量估计不同数量水平的解释变量对被解释变量的影响.在经济关系中,当解释变量X的值达到某一水平X′之前,与被解释变量之间存在某种线性关系:当解释变量X的值达到或者超过X′以后,与被 ...
- 高级建模方法(Gurobi):线性化、分段线性函数、逻辑变量/约束
文章目录 乘积式 绝对值 max∣x∣\max |x|max∣x∣, big-M max∣x∣\max |x|max∣x∣, SOS-1约束 分段线性函数 SOS2表示法 0-1变量表示法 min ...
- 灰度变换-分段线性函数
目录 1. 分段线性函数介绍 2. 代码实现 3. other 1. 分段线性函数介绍 分段线性函数同样是点运算,基于像素的图像增强,也就是对比度拉伸. 大概的原理就是:将不同灰度区间的灰度值经过不同 ...
- python二元多次函数拟合_对python实现二维函数高次拟合的示例详解
在参加"数据挖掘"比赛中遇到了关于函数高次拟合的问题,然后就整理了一下源码,以便后期的学习与改进. 在本次"数据挖掘"比赛中感觉收获最大的还是对于神经网络的认识 ...
- 【使用华为云MindSpore框架实现线性函数拟合___实验报告(一)】
*************************************************** 码字不易,收藏之余,别忘了给我点个赞吧! *************************** ...
- Python:实现segmented sieve分段筛算法(附完整源码)
Python:实现segmented sieve分段筛算法 import math def sieve(n):"""Segmented Sieve."" ...
最新文章
- 在react next 中使用rem_CSS中的rem:了解和使用rem单位
- 记录我学github的路程(二)
- linux 日志服务器简单用
- axios发送post数据后台收不到_使用axios post 提交数据,后台获取不到
- django批量修改table_django-formset实现数据表的批量操作
- 【nodejs】安装browser-sync 遇到错误提示
- 初探EntityFramework——空Code First模型
- 【教程】VsCode搭建Java开发环境
- 设计模式学习笔记——责任链(Chain of Responsibility)模式
- Linux软件安装及基本概念
- excel导入到mysql
- CRF++总结 (1)(2)【转】
- linux输入法小结
- 【原创】告别恼人的水平滚动条——滚动条宽度到底是多少?
- 2017 计蒜之道 初赛 第五场 A. UCloud 机房的网络搭建
- 【建议收藏】10个适合程序员逛的在线社区
- 【围棋游戏——使用Python实现(纯tkinter gui)】
- 无损连接,函数依赖性判定
- 学习fortran77基础语法
- LRU算法,走迷宫,数根,星际战争