解释 scipy 差分进化中顶级参数的直觉

2024-02-26

我正在使用差异进化 https://docs.scipy.org/doc/scipy-0.15.1/reference/generated/scipy.optimize.differential_evolution.htmlscipy 中的优化器,我不明白 tol 参数背后的直觉。具体是文档中说的:

tol:浮动,可选

当总体能量的平均值乘以 tol,除以 总体能量的标准差大于 1 求解过程终止:
convergence = mean(pop) * tol / stdev(pop) > 1

从用户角度来看,设置 tol 代表什么?


也许文档中的公式以下面的形式更容易理解(参见第 508 行和 526 行)the code https://github.com/scipy/scipy/blob/v0.15.1/scipy/optimize/_differentialevolution.py#L508):

std(population_energies) / mean(population_energies) < tol

这意味着当群体中每个个体的能量标准差(以平均值为范数)小于给定的容差值时,就达到了收敛。

优化算法是迭代的。每次迭代都会找到更好的解决方案。容差参数用于定义停止条件。停止条件实际上是所有个体(参数集)具有大致相同的能量,即相同的成本函数值。然后,返回给出最低能量的参数集作为解。

它还意味着所有个体在参数空间中彼此相对接近。因此,对于下一代来说,不能指望有更好的解决方案。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

解释 scipy 差分进化中顶级参数的直觉 的相关文章

  • 三次样条内存错误

    在具有 4GB 内存的计算机上 这种简单的插值会导致内存错误 基于 http docs scipy org doc scipy reference tutorial interpolate html http docs scipy org
  • 在Python中计算矩阵乘以其转置(AA^T)的最快方法

    在Python中将矩阵与其转置 AA T 相乘的最快方法是什么 我认为 NumPy SciPy 没有考虑使用例如时涉及的对称性 np dot or np matmul 得到的矩阵总是对称的 所以我可以想象有一个更快的解决方案 None
  • 如何获得凸包中均匀分布的点?

    给定一组点 points np random randn n 3d points 我想均匀地填充由凸包定义的体积 其中它们位于一个列表 np array of shapenx3 的 3d 点 我可以得到凸包 hull scipy spati
  • 如何通过异常值检测方法在周期性或基于序列的数据上生成脉冲作为异常值以进行实验?

    我想对一些时间序列数据进行一些实验KM https scikit learn org stable auto examples cluster plot cluster iris html sphx glr auto examples cl
  • 约束 3D 表面的 RBF 插值以保持曲率

    我的任务是开发一种算法 给定一组表示现有表面测量值的稀疏点 我们就可以计算表面上任何点的 z 坐标 面临的挑战是找到一种合适的插值方法 该方法可以在仅给定几个点的情况下重新创建 3D 表面 并推断出超出包含初始测量值的范围的值 对于许多插值
  • ValueError:对于optimize.curve_fit中所需的数组来说对象太深

    我正在尝试拟合化学系统中四个变量 A B C D 的人口增长和衰退的动力学模型 我正在尝试求解以下一组方程 我已将其以矩阵形式附加 方程的矩阵形式 https i stack imgur com ysEdZ png 其中 t 是时间步长 k
  • Python 或 C 语言中的 Matlab / Octave bwdist()

    有谁知道 Matlab Octave bwdist 函数的 Python 替代品 此函数返回给定矩阵的每个单元格到最近的非零单元格的欧几里得距离 我看到了一个 Octave C 实现 一个纯 Matlab 实现 我想知道是否有人必须用 AN
  • 使用 Numpy 进行多维批量图像卷积

    在图像处理和分类网络中 一个常见的任务是输入图像与一些固定滤波器的卷积或互相关 例如 在卷积神经网络 CNN 中 这是一种极其常见的操作 我已将通用版本任务减少为 Given 一批 N 个图像 N H W D 和一组 K 个滤镜 K H W
  • 使用 python 生成器高效创建 scipy.lil_matrix

    我有一个生成单一维度的生成器numpy arrays 的长度相同 我想要一个包含该数据的稀疏矩阵 行的生成顺序与我希望它们出现在最终矩阵中的顺序相同 csr矩阵优于lil矩阵 但我认为后者在我描述的场景中更容易构建 假设row gen是一个
  • 在 scipy 中创建新的发行版

    我试图根据我拥有的一些数据创建一个分布 然后从该分布中随机抽取 这是我所拥有的 from scipy import stats import numpy def getDistribution data kernel stats gauss
  • 使用 scipy 的solve_bvp 求解 BVP

    我有一个由 3 个微分方程组成的系统 我相信从代码中可以明显看出 具有 3 个边界条件 我设法在 MATLAB 中用一个循环来解决这个问题 一点一点地改变最初的猜测 而不会在程序即将返回错误时终止程序 然而 关于scipy s solve
  • 将 scipy 稀疏矩阵的几行采样到另一个中

    如何对 scipy 稀疏矩阵的某些行进行采样 并从这些采样的行中形成一个新的 scipy 稀疏矩阵 例如 如果我有一个 10 行的 scipy 稀疏矩阵 A 并且我想创建一个新的 scipy 稀疏矩阵 B 其中 A 的第 1 3 4 行 该
  • 为什么我在 python 中的 Spearman 相关性中得到 nan

    我在用scipy来计算相关性 我计算斯皮尔曼相关性的代码如下 from scipy import stats sequence 1 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 sequence 2 0 0
  • 将分布拟合到直方图

    I want to know the distribution of my data points so first I plotted the histogram of my data My histogram looks like th
  • 如何在Python中拟合阶跃函数

    我有一个关于使用 curve fit 等 scipy 例程拟合阶跃函数的问题 我很难将其矢量化 例如 import numpy as np from scipy optimize import curve fit import matplo
  • 从值数组中计算 sympy 表达式

    我正在尝试 sympy 但遇到了一个无法解决的问题 使用 scipy 我可以编写一个表达式并计算 x 值数组 如下所示 import scipy xvals scipy arange 100 100 0 1 f lambda x x 2 f
  • 如何安装 scipy 杂项包

    我已经安装 实际上是重新安装 scipy 10 x86 64 whl 19 8MB 19 8MB downloaded Installing collected packages scipy Successfully installed s
  • 切片稀疏(scipy)矩阵

    我将不胜感激任何帮助 以理解从 scipy sparse 包中切片 lil matrix A 时的以下行为 实际上 我想根据行和列的任意索引列表提取子矩阵 当我使用这两行代码时 x1 A list 1 x2 x1 list 2 一切都很好
  • 使用 scipy curve_fit 拟合噪声指数的建议?

    我正在尝试拟合通常按以下方式建模的数据 def fit eq x a b c d e return a 1 np exp x b c np exp x d e x np arange 0 100 0 001 y fit eq x 1 1 1
  • 使用 matplotlib 设置或固定二元分布值

    I ve animated a bivariate gaussian distribution using matplotlib 我已经计算过这个distribution通过调整COV matrix来考虑特定的变量 我可以提供有关此过程的更

随机推荐

  • R/LaTeX 表创建包的建议

    我一直在使用xtable很长一段时间以来 并期待着用 R 编写我的第一个包 所以我认为 如果我有一些值得实施的 酷 想法 那么很有可能有人在我之前到达那里 我对专门用于 LaTeX 表创建的函数 包感兴趣 当然是通过 R 我撞上了quant
  • 如何复制 Google 工作表并保留其项目触发器和脚本?

    我目前有一个 Google 表格 用作主模板 也就是说 我会为每个请求制作该模板的副本 我想向我的主模板添加一个 Google 应用程序脚本 当工作表完成后 它会在编辑时发布到我的服务器 该脚本将被复制并为此模板的每个副本运行 我已经尝试从
  • 是否应该确保 IDisposable.Dispose() 可以安全地多次调用?

    应实施IDisposable make Dispose 多次拨打安全吗 或者相反 大多数 NET Framework 类采用什么方法 具体来说 打电话安全吗 System Data Linq DataContext Dispose 多次 我
  • 建立 iTunes 商店的联属链接而不进行重定向?

    苹果公司在 从 iPhone 应用程序启动 App Store http developer apple com library ios qa qa1629 index html 如何建立一个到应用商店的联属链接并在后台处理重定向 这样就不
  • LLDB GUI 的文档

    我最近一直在研究 LLDB 的 GUI 功能 有关此功能的 stackoverflow 链接有描述here https stackoverflow com questions 19820247 is there a tui mode for
  • MongoDB 3.2 身份验证失败

    我使用以下命令集创建一个用户 这应该在两者中创建用户admin数据库以及我的目标数据库 c2d mongo 127 0 0 1 27017 MongoDB shell version 3 2 6 29 g5c19788 connecting
  • 使用epplus创建.xlsm文件

    我正在尝试让网站导出 xlsm 文件 但似乎找不到任何有帮助的内容 这是我用来测试构建 xlsx 文件 有效 的一个简单示例 using OfficeOpenXml div Change file extension to xlsm to
  • 如何在 PHP 中解析日期字符串?

    日期字符串为Apr 30 2010 我怎样才能将字符串解析为2010 04 30使用 PHP 使用 DateTime API 需要 PHP 5 3 dateTime DateTime createFromFormat F d Y Apr 3
  • 表达请求 .query.xyz 作为字符串 |细绳[]

    我正在尝试为我的变量分配一个类型 该变量从express 请求查询中获取其值 req query 的类型为 QueryString ParsedQs 例如 req query accountId 的类型为 string QueryStrin
  • 是否有一个玩笑配置会导致 console.warn 测试失败?

    如何配置笑话测试以因警告而失败 console warn stuff fail test 您可以使用这个简单的覆盖 let error console error console error function message error a
  • R 的 xtsum 命令?

    我们正在处理面板数据 Stata中有一个命令 xtsum 这给出了数据集中变量的内部方差和方差之间 R 是否有类似的命令可以产生干净的输出 我用了一个小函数来做到这一点 函数 XTSUM 接受三个输入 数据 数据集varname xtsum
  • C# 嵌套 Try Catch 语句或方法?

    简单的最佳实践问题 应该嵌套 try catch 语句还是只使用方法 例如 如果您有一个打开文件的方法确实有效并关闭该文件 那么您将在 try catch 之外进行打开和关闭 或者更确切地说 在 finally 块中进行关闭 现在 如果您的
  • 如何为 'a:before' 和 'a:after' 编写 ':hover' 条件?

    我怎样才能写 hover and visited的条件a before 我想a before hover 但它不起作用 这取决于您实际想要做什么 如果您只是想将样式应用到 before伪元素 当a元素匹配一个伪类 你需要写a hover b
  • 将结构体转换为数组是否合法?

    考虑以下 Just a sequence of adjacent fields of same the type repr C derive Debug struct S
  • 使用 MongoDB 与 MySQL 与大量 JSON 字段?

    有一种微博类型的应用程序 两个主要的基本数据库存储是 MySQL 或 MongoDB 我计划对大量数据进行非规范化 即对帖子进行的投票存储在投票表中 并且主帖子表中的计数也会增加 该帖子还涉及其他操作 例如点赞 投票否决 如果我使用 MyS
  • WindowManager:MainActivity 泄露了原始窗口 com.android.internal.policy.impl.PhoneWindow$DecorView@40731aa0

    这个警告是什么意思 我开发了一个 Android 应用程序 我在 logcat 中看到了这个警告 此警告不会导致关闭应用程序 一切都运行没有任何问题或中断 但它是否会导致一些在我的应用程序中看不到的问题 Continue logcat at
  • 为 RSpec 套件中的所有示例设置一次变量(不使用全局变量)

    设置变量一次以供 RSpec 套件中的所有示例使用的传统方法是什么 我目前在中设置了一个全局变量spec helper检查规范是否正在 调试模式 下运行 debug ENV key DEBUG ENV DEBUG casecmp false
  • 将“onclick”事件附加到 D3 图表背景

    我有一个 D3 直方图 我在其中的条形图上附加了一个 onclick 事件 var bar svg selectAll bar data data enter append g attr class bar attr transform f
  • 我可以引用外部 DTD 来了解 DTD 中 ELEMENT 的结构吗?

    我可以定义一个 DTD 使其中的元素的结构在外部 DTD 中定义吗 我的意思是类似以下内容 这是一个我知道在某些方面无效的示例 但希望它能让您了解我正在寻找的内容
  • 解释 scipy 差分进化中顶级参数的直觉

    我正在使用差异进化 https docs scipy org doc scipy 0 15 1 reference generated scipy optimize differential evolution htmlscipy 中的优化