python拟合二元一次函数
背景:
参考:
代码:
import numpy as np
from scipy.optimize import leastsq
import pylab as pl
def func(x,p):
"""
数组拟合函数
"""
A,k,theta = p
return A*(x-k)**2+theta
def residuals(p,y,x):
"""
残差
"""
return y-func(x,p)
x = np.linspace(0,2,100)
A,k,theta = 10.,1,2. #真实数据参数
y0 = func(x,[A,k,theta]) #真实数据
y1 = y0 + 2 * np.random.randn(len(x)) #加入噪声序列
p0 = [7.,0.2,1.]
plsq = leastsq(residuals,p0,args = (y1,x))
print("真实参数:",[A,k,theta])
print("拟合参数:",plsq[0]) #试验数据拟合后的参数
pl.plot(x,y0,label = "real")
pl.plot(x,y1,label = "real+noise")
pl.plot(x,func(x,plsq[0]),label = "fitting")
pl.legend()
pl.show()
结果:
(貌似这里的求解方法用了智能算法,因为每次的结果都有细小差异。具体资料没见到,以后有精力再找)
真实参数: [10.0, 1, 2.0]
拟合参数: [ 11.06873474 0.99709893 1.49107462]
python拟合二元一次函数相关推荐
- python二元多次函数拟合_Python实现的拟合二元一次函数功能示例【基于scipy模块】...
本文实例讲述了Python实现的拟合二元一次函数功能.分享给大家供大家参考,具体如下: 背景: 使用scipy拟合一元二次函数. 参考: HYRY Studio-<用Python做科学计算> ...
- python自定义二元一次函数_Python实现的拟合二元一次函数功能示例【基于scipy模块】...
本文实例讲述了Python实现的拟合二元一次函数功能.分享给大家供大家参考,具体如下: 背景: 使用scipy拟合一元二次函数. 参考: HYRY Studio-<用Python做科学计算> ...
- python二元函数图像在线绘制_Python实现的拟合二元一次函数功能示例【基于scipy模块】...
Python实现的拟合二元一次函数功能示例[基于scipy模块] 发布时间:2020-08-31 00:46:24 来源:脚本之家 阅读:121 本文实例讲述了Python实现的拟合二元一次函数功能. ...
- Python解二元一次方程
使用Python解二元一次方程组 用到的库:sympy 下面举一个例子: 解方程组:y = x + 1,y = -x + 1 #导入sympy from sympy import *#定义变量 x = ...
- python拟合统计分布_利用 Python Fitter 判断数据样本符合那种分布函数
# Python拟合数据样本的分布 # 安装fitter # pip install fitter # 生成一段模拟数据 from scipy import stats import numpy as ...
- 数值计算——最小二乘拟合二元一次多项式
数值计算--最小二乘拟合二元一次多项式 最小二乘拟合: 就是根据一系列给定的数据点,求一条曲线使得数据点到曲线的某些(水平.竖直.垂直)距离最短. 推导过程: 1. 设拟合多项式为: 2. ...
- 麻雀搜索算法(SSA)求解二元一次函数(附完整代码)
目标函数: function z =fun( x,y ) z=-20*exp(-0.2*sqrt(0.5*(x.^2+y.^2)))-exp(0.5*(cos(2*pi*x)+cos(2*pi*y)) ...
- 数学建模入门-python拟合曲面
前言 找了好久python拟合曲面的方法,没找到,就借鉴 https://blog.csdn.net/Haipai1998/article/details/85345823 这篇博客,对方法进行封装更 ...
- python定义类方法用来计算整数的减法运算_详解python的二元算术运算,为什么说减法只是语法糖?...
原标题:详解python的二元算术运算,为什么说减法只是语法糖? 原题 | Unravelling binary arithmetic operations in Python 作者 | Brett ...
- 用最小二乘法拟合二元多次曲线
引用 http://blog.sina.com.cn/s/blog_6e51df7f0100thie.html 对代码稍作修改和注释,防止链接失效. 1 ///<summary> 2 // ...
最新文章
- UVa572 Oil Deposits DFS求连通块
- react native redux 梳理
- ByteArrayOutputStream用法
- 《Effective Java》读书笔记 - 5.泛型
- 8. Action过滤
- 宏锦软件 Android 的 ListView 使用详解
- bd3.1 Python 高级
- 动态规划—最长上升子序列(POJ 1458)
- 计算机无法备份,电脑无法备份怎么办
- 取出被正则表达式匹配的值
- MySQL数据库索引优化实战分析
- 第三季-第14课-有名管道通讯编程
- 安装debian文件管理器
- 大型网站架构的演化[转]
- 解决cef加载flash时弹出黑框的问题
- 计算机发明于1946年用英语怎么说,电子计算机发明于哪一年,电子计算机发明与1946年...
- ctfshow密码easyrsa5
- 操作系统有哪些.智能手机的操作系统有哪些?
- C/C++ 判断当前操作系统语言、获取当前操作系统使用语言 GetSystemDefaultLangID
- web--拉灯泡切换黑天与白夜的精美动画
热门文章
- 【智能路由器】轻量级web服务器lighttpd架设——打造家庭影院
- 视频教程-SAP快速入门-其他
- linux内核剖析之traps.c
- 神策分析 Android SDK 之用户路径采集
- mugen linux主程序,MUGEN主程序
- java httpserver stop_java web tomcat服务停止Stopping Coyote HTTP/1.1 on http-8089
- 蜂鸟处理器+OV5640摄像头模块开发
- 高通9008端口刷linux,高通黑砖强制9008_高通救砖9008刷底层
- 【愚公系列】2022年10月 微信小程序-电商项目-收货地址功能实现
- Python(十三)企业微信Api自动化框架搭建