为了澄清一点。有两个相关的问题:

最小化功能

拟合模型数据

要拟合模型对观测数据是要找到一个模型的这种参数,减少某种模型数据与观测数据之间的错误。 (x)可以是一个向量)。

F(x) = 0.5 * sum(rho(f_i(x)**2), i = 0, ..., m - 1)

(rho是亏损的功能,默认是rho(x) = x所以不介意它)

least_squares(func, x0)预计调用func(x)将返回一个矢量[a1, a2, a3, ...]为其平方和,将计算:S = 0.5 * (a1^2 + a2^2 + a3^2 + ...)。

least_squares将调整x0以尽量减少S。

因此,为了用它来拟合模型数据,必须构建误差的函数的模型和实际数据之间 - 残差然后最小化残差功能。

import numpy as np

from scipy.optimize import least_squares

x = np.asarray(range(20,401,20))

y = np.random.rand(20)

params = np.array([1,1])

def funcinv(x, a, b):

return b + a/x

def residuals(params, x, data):

# evaluates function given vector of params [a, b]

# and return residuals: (observed_data - model_data)

a, b = params

func_eval = funcinv(x, a, b)

return (data - func_eval)

res = least_squares(residuals, params, args=(x, y))

这给出了一个结果:在你的情况,你可以如下写 print(res)

...

message: '`gtol` termination condition is satisfied.'

nfev: 4

njev: 4 optimality: 5.6774618339971994e-10

status: 1

success: True

x: array([ 6.89518618, 0.37118815])

但是,作为一个残差功能几乎相同的所有的时间(res = observed_data - model_data),有scipy.optimize中的快捷方式叫做curve_fit:curve_fit(func, xdata, ydata, x0)。 curve_fit自动构建残差函数,你可以简单地写:

import numpy as np

from scipy.optimize import curve_fit

x = np.asarray(range(20,401,20))

y = np.random.rand(20)

params = np.array([1,1])

def funcinv(x, a, b):

return b + a/x

res = curve_fit(funcinv, x, y, params)

print(res) # ... array([ 6.89518618, 0.37118815]), ...

python画反比例函数_拟合一个反比例函数相关推荐

  1. python画交互式地图_使用Python构建交互式地图-入门指南

    python画交互式地图 Welcome to The Beginner's Guide to Building Interactive Maps in Python 欢迎使用Python构建交互式地 ...

  2. python画圆形螺旋线_宝宝爱看小猪佩奇,很简单,让我们用python搞定它

    现在很多宝宝喜欢看小猪佩奇,今天就教大家用python的海龟画图画一个乖巧萌萌的小猪佩奇,引导对编程产生浓浓 的兴趣. 画图前引导 1.让我们打开百度,输入python进入python官网 2.在官网 ...

  3. python如何拟合函数_我们如何在Python中拟合一个sigmoid函数?(How do we fit a sigmoid function in Python?)...

    由于可重现性的原因,我正在共享我在这里工作的简单数据集. 为了弄清楚我在做什么--从第2列开始,我读取当前行并将其与前一行的值进行比较.如果它更大,我会不断比较.如果当前值小于前一行的值,我想用当前值 ...

  4. 用python画满天星花朵_如何使用python做一个可以画一朵花的小程序?

    这篇文章可以交给你怎么在windows或macOS与python3环境下做出一个可以根据你的自定义设置画出一朵花的程序. 在开始写正式的逻辑代码之前,我们要先导入一些库: from tkinter i ...

  5. python和报表生成工具_推荐一个利用 python 生成 pptx 分析报告的工具包:reportgen...

    reportgen v0.1.8 更新介绍 这段时间,我对 reportgen 进行了大工程量的修改和更新.将之前在各个文章中出现的函数进行了封装,同时也对现有工具包的一些逻辑进行了调整. 1.rep ...

  6. 用python画机器猫代码_如何用Python画一只机器猫?| 原力计划

    原标题:如何用Python画一只机器猫?| 原力计划 作者 | 人邮异步社区 责编 | 胡巍巍 出品 | CSDN博客 自信心是成功的源泉,对刚入门编程行业的初级程序员来说,多敲代码多做项目就是构建自 ...

  7. python解魔方程序_写一个解二阶魔方的程序

    本文需要读者有一定的魔方基础, 最起码也要达到十秒内还原二阶魔方的水平, 并且手上最好有一个二阶魔方, 否则文中的很多东西理解不了. 另外, 这里使用的算法是我自己写着玩的, 如果你需要更成熟和专业的 ...

  8. python画e指数函数_[宜配屋]听图阁

    这里用Python逼近函数y = exp(x);同样使用泰勒函数去逼近: exp(x) = 1 + x + (x)^2/(2!) + .. + (x)^n/(n!) + ... #!/usr/bin/ ...

  9. python画大学教室_开学季,用Python画大学教室座位神分区图,你中枪了吗?

    这就用Python给大家画这幅有趣的神分区图,看看是不是跟你的情况一样呢? 基本思路:我们将画布分成若干子画布,这些子画布构成绘图区域,在这些绘图区域上分别绘制图形,这就让我们想到函数subplot2 ...

  10. python画简单圣诞树_教你用Python画了一棵圣诞树!圣诞节给女朋友制作个小惊喜...

    如何用Python画一个圣诞树呢? 最简单: height = 5 stars = 1 for i in range(height): print((' ' * (height - i)) + ('* ...

最新文章

  1. Android4.0 Design之UI设计易犯的错误2
  2. Logistic Regression and Classification
  3. Apache Doris : 一个开源 MPP 数据库的架构与实践
  4. P3250 [HNOI2016]网络(利用堆建线段树 + 树剖)
  5. TP-LINK WR841N V8刷OpenWRT
  6. 创维酷开电视能换成android系统吗,创维酷开电视怎么系统升级【图文教程】
  7. LoadRunner11完美破解
  8. ZYNQ研究----(2)基于开发板制作串口测试程序
  9. 奇迹服务端架设教程技术分享探究——学习
  10. SYNOPSYS AXI VIP using note
  11. MT6573驱动开发日志之touchpanel
  12. 存储基础知识--存储网络
  13. CT重建学习笔记(一)
  14. r语言t检验输出检验统计量_R语言:常用统计检验
  15. 「Python海龟画图」利用海龟画笔绘制奥运五环
  16. 图像拼接matlab
  17. STM32自学笔记-7-CAN总线通信
  18. 水井坊VS舍得酒业,谁是A股川酒“老三“?
  19. Open mv识别三角形的办法
  20. html直线箭头,HTML中利用div+CSS实现简单的箭头图标的代码

热门文章

  1. 鸿蒙判后而成翻译,我在洪荒一百四十六石猴出世(1/3)胜神州。东胜神洲,海外有一国土,名曰傲来-金羚文学...
  2. 无论用手工处理还是用计算机进行处理,会计电算化试卷
  3. java微信公众号上传永久素材_微信开放平台永久素材视频文件上传
  4. 转:以Delphi Package架构多人开发应用程序环境
  5. matlab 安装时报错: dsp_doc_en_common 时检测到以下错误
  6. jQuery02(达)
  7. D. Berserk And Fireball(Educational Codeforces Round 91 (Rated for Div. 2))
  8. 关于修复msvcp110.dll丢失的问题
  9. 专访:6位飞桨开发者的「新生代」AI创新故事!
  10. mstar 平台GPIO在线调试