Python基于最小二乘法实现曲线拟合示例
作者:Inside_Zhang 时间:2021-08-06 15:47:07
本文实例讲述了Python基于最小二乘法实现曲线拟合。分享给大家供大家参考,具体如下:
这里不手动实现最小二乘,调用scipy库中实现好的相关优化函数。
考虑如下的含有4个参数的函数式:
构造数据
import numpy as np
from scipy import optimize
import matplotlib.pyplot as plt
def logistic4(x, A, B, C, D):
return (A-D)/(1+(x/C)**B)+D
def residuals(p, y, x):
A, B, C, D = p
return y - logisctic4(x, A, B, C, D)
def peval(x, p):
A, B, C, D = p
return logistic4(x, A, B, C, D)
A, B, C, D = .5, 2.5, 8, 7.3
x = np.linspace(0, 20, 20)
y_true = logistic4(x, A, B, C, D)
y_meas = y_true + 0.2 * np.random.randn(len(y_true))
调用工具箱函数,进行优化
p0 = [1/2]*4
plesq = optimize.leastsq(residuals, p0, args=(y_meas, x))
# leastsq函数的功能其实是根据误差(y_meas-y_true)
# 估计模型(也即函数)的参数
绘图
plt.figure(figsize=(6, 4.5))
plt.plot(x, peval(x, plesq[0]), x, y_meas, 'o', x, y_true)
plt.legend(['Fit', 'Noisy', 'True'], loc='upper left')
plt.title('least square for the noisy data (measurements)')
for i, (param, true, est) in enumerate(zip('ABCD', [A, B, C, D], plesq[0])):
plt.text(11, 2-i*.5, '{} = {:.2f}, est({:.2f}) = {:.2f}'.format(param, true, param, est))
plt.savefig('./logisitic.png')
plt.show()
PS:这里再为大家推荐两款相似的在线工具供大家参考:
在线多项式曲线及曲线函数拟合工具:
http://tools.jb51.net/jisuanqi/create_fun
在线绘制多项式/函数曲线图形工具:
http://tools.jb51.net/jisuanqi/fun_draw
希望本文所述对大家Python程序设计有所帮助。
来源:https://blog.csdn.net/lanchunhui/article/details/50358943
标签:Python,最小二乘法,曲线拟合
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
PHP析构函数destruct与垃圾回收机制的讲解
2023-07-11 21:05:36
谷歌浏览器Chrome的javascript引擎
2008-09-04 12:24:00
![](https://img.aspxhome.com/file/UploadPic/20089/4/200894122755641s.gif)
开发保留标准的浏览器功能的AJAX应用
2008-02-03 14:03:00
解读ASP.NET 5 & MVC6系列教程(6):Middleware详解
2023-07-23 22:27:34
基于google图表API的jquery组件 I
2010-03-01 10:20:00
![](https://img.aspxhome.com/file/UploadPic/20103/1/01-93s.png)
论坛首页效果图设计
2009-03-19 13:46:00
![](https://img.aspxhome.com/file/UploadPic/20093/19/newbbs-600x1523-48s.jpg)
PHP自动生成月历代码
2023-11-18 09:22:52
SQL Server连接失败错误及解决
2008-01-28 21:09:00
![](https://img.aspxhome.com/file/UploadPic/20081/28/2008128211959685s.jpg)
ASP.NET 2.0中的数据操作之九:跨页面的主/从报表
2023-07-19 20:28:42
![](https://img.aspxhome.com/file/2023/5/75525_0s.png)
关于document.cookie的使用
2008-03-25 12:07:00
:hover在IE6下的问题
2009-06-18 21:09:00
非常全面的实用JavaScript开发工具列表
2010-01-05 16:44:00
![](https://img.aspxhome.com/file/UploadPic/20101/11/200912311344278032-17s.jpg)
Cookies 欺骗漏洞的防范方法(vbs+js 实现)
2011-03-09 11:09:00
thinkPHP框架实现类似java过滤器的简单方法示例
2023-11-22 12:24:47
如何提高Request集合的使用效率?
2010-06-07 20:53:00
php开启mysqli扩展之后如何连接数据库
2023-11-14 17:36:08
php中ob_flush函数和flush函数用法分析
2023-11-15 06:12:59
thinkphp学习笔记之多表查询
2023-11-15 02:57:15
asp如何选择访问速度最快的站点?
2010-06-10 18:34:00
SQL 判断给定日期值(或时间段)所在星期的星期一和星期天的日期
2011-10-24 20:14:52