- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我的问题是:是否有最小化算法,最好用 Python 实现,可以对一个慢速(~1-10 秒)的函数进行操作,并从实时系统中获取数据,这不会超过几个几个小时才能完成?
我有一个 FPGA,它对一些传感器数据运行一个过滤器,并使用这个过滤器的输出来提高另一个设备的性能。我想找到最佳过滤器。我对系统建模和使用各种信号处理技术的尝试没有产生足够的结果,所以现在我将尝试在实时系统本身上解决这个问题(如果有的话,只是为了证明这样的最佳过滤器是可能的)。
可以通过串行线路对滤波器进行编程,并且可以通过串行线路测量其他设备的性能。
所以我可以构造一个函数:
这意味着我有一个函数可以用作最小化的目标。但问题是:
对滤波器进行编程大约需要 1.5 秒,获取数据以衡量滤波器的优劣大约需要 6 秒。总而言之,每个函数调用将近 8 秒。换句话说,仅仅调用它 500 次就需要一个多小时。即使加快通信和计算速度也可能不会改变这一点。
(注意下面的x
是我的目标函数的参数空间中的一个向量。)
简单来说,x1 == x2
并不意味着 f(x1) == f(x2)
。由于系统的噪声,在其参数空间中的同一点对目标函数 f(x)
进行采样可能会由于系统的噪声而产生不同的结果。
我想到的第一件事是让目标函数实际平均几个测量值,并增加我正在运行的任何最小化例程的容差值。但在查看实际数字时,在最坏的情况下,我可以让 f(x)
的(平均)值在整个参数范围内变化 2.0,但样本标准偏差为 1.6。这意味着如果我想将标准误差 (s/sqrt(n)
) 降低到 0.1,我需要测量同一个点 250 次,这使得每次测量需要 30 分钟。是的。
我可以使用一些技巧来改进这一点,比如在任何给定点以 0.25 的标准偏差在参数范围内获得 ~20 的摆动。但这些技巧在时间上还有其他权衡。
从好的方面来说,在整个优化空间(我已经这样做以确认确实存在全局最小值)上绘制函数(大大平均)表明事情实际上相当平滑并且最小值不是太尖锐了。另一个好的方面是指标只需要优化到两个或三个有效数字。如果它不是那么慢,优化它会很容易。
我已经开始查看 minimisation routines in SciPy ,但由于许多参数未记录或相互依赖,因此有点像在黑暗中行走(每一步都需要几个小时)。
令我震惊的是,我真正需要的是一种已知可以在最少的函数调用次数中工作的优化算法;尽管也许还有另一种我没有考虑过的方法。
最佳答案
包 scikit-optimize (skopt) 正是为这种设置而设计的:缓慢、嘈杂的目标函数。它使用高斯过程对目标函数建模,并在不确定的评估点(以改进模型)和可能好的点之间切换。他们的示例使用约 100 次评估来恢复最小值。甚至还有一个针对物理实验的界面,它会在其中提出试验值,您运行实验,将结果提供给它,然后它会提出更多试验值。
关于python - 最小化缓慢、嘈杂、未明确定义的目标函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/8469045/
C语言sscanf()函数:从字符串中读取指定格式的数据 头文件: ?
最近,我有一个关于工作预评估的问题,即使查询了每个功能的工作原理,我也不知道如何解决。这是一个伪代码。 下面是一个名为foo()的函数,该函数将被传递一个值并返回一个值。如果将以下值传递给foo函数,
CStr 函数 返回表达式,该表达式已被转换为 String 子类型的 Variant。 CStr(expression) expression 参数是任意有效的表达式。 说明 通常,可以
CSng 函数 返回表达式,该表达式已被转换为 Single 子类型的 Variant。 CSng(expression) expression 参数是任意有效的表达式。 说明 通常,可
CreateObject 函数 创建并返回对 Automation 对象的引用。 CreateObject(servername.typename [, location]) 参数 serv
Cos 函数 返回某个角的余弦值。 Cos(number) number 参数可以是任何将某个角表示为弧度的有效数值表达式。 说明 Cos 函数取某个角并返回直角三角形两边的比值。此比值是
CLng 函数 返回表达式,此表达式已被转换为 Long 子类型的 Variant。 CLng(expression) expression 参数是任意有效的表达式。 说明 通常,您可以使
CInt 函数 返回表达式,此表达式已被转换为 Integer 子类型的 Variant。 CInt(expression) expression 参数是任意有效的表达式。 说明 通常,可
Chr 函数 返回与指定的 ANSI 字符代码相对应的字符。 Chr(charcode) charcode 参数是可以标识字符的数字。 说明 从 0 到 31 的数字表示标准的不可打印的
CDbl 函数 返回表达式,此表达式已被转换为 Double 子类型的 Variant。 CDbl(expression) expression 参数是任意有效的表达式。 说明 通常,您可
CDate 函数 返回表达式,此表达式已被转换为 Date 子类型的 Variant。 CDate(date) date 参数是任意有效的日期表达式。 说明 IsDate 函数用于判断 d
CCur 函数 返回表达式,此表达式已被转换为 Currency 子类型的 Variant。 CCur(expression) expression 参数是任意有效的表达式。 说明 通常,
CByte 函数 返回表达式,此表达式已被转换为 Byte 子类型的 Variant。 CByte(expression) expression 参数是任意有效的表达式。 说明 通常,可以
CBool 函数 返回表达式,此表达式已转换为 Boolean 子类型的 Variant。 CBool(expression) expression 是任意有效的表达式。 说明 如果 ex
Atn 函数 返回数值的反正切值。 Atn(number) number 参数可以是任意有效的数值表达式。 说明 Atn 函数计算直角三角形两个边的比值 (number) 并返回对应角的弧
Asc 函数 返回与字符串的第一个字母对应的 ANSI 字符代码。 Asc(string) string 参数是任意有效的字符串表达式。如果 string 参数未包含字符,则将发生运行时错误。
Array 函数 返回包含数组的 Variant。 Array(arglist) arglist 参数是赋给包含在 Variant 中的数组元素的值的列表(用逗号分隔)。如果没有指定此参数,则
Abs 函数 返回数字的绝对值。 Abs(number) number 参数可以是任意有效的数值表达式。如果 number 包含 Null,则返回 Null;如果是未初始化变量,则返回 0。
FormatPercent 函数 返回表达式,此表达式已被格式化为尾随有 % 符号的百分比(乘以 100 )。 FormatPercent(expression[,NumDigitsAfterD
FormatNumber 函数 返回表达式,此表达式已被格式化为数值。 FormatNumber( expression [,NumDigitsAfterDecimal [,Inc
我是一名优秀的程序员,十分优秀!