- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我有一段较大的代码,其核心是将函数与数据相匹配。要拟合的数据和函数是动态的。最近我向整个系统附加了一个额外的数据点,现在 curve_fit
总是返回初始猜测(或太接近它的东西),无论我如何选择它。对于非常不同的 y 值和 x 值(十组前者,两组后者)会发生这种情况。
我知道选择起始值很重要,但我以前从未遇到过使用默认值的问题(我的函数通常很简单)并且可以通过取消注释添加附加数据的新代码来恢复到正常工作的状态观点。现在人们会认为新代码显然是问题所在,但在新添加和实际将数据馈送到 curve_fit
之间有相当多的步骤。我已经检查过 curve_fit
的输入类型是否相同:np.ndarray
,在有问题的情况下只有一个元素长。
但是,在创建 MWE 时,我注意到它只是导致问题的确切 x 数组。当我在我的 MWE 中复制主程序的打印 x 向量而不是内部表示时,它完全消失了。因此我只能用外部文件来说明问题:local_z.npy [150kB]
MWE:
import numpy as np
from scipy.optimize import curve_fit
values = np.array([[1.37712972, 1.58475346, 1.78578759, 1.9843099, 1.73393093],
[-0.0155715, -0.01534987, -0.00910744, -0.00189728, -1.73393093],
[1.23613934, 0.76894505, 0.18876817, 0.06376843, 1.1637315 ],
[0.8535248, 0.53093829, 0.13033993, 0.04403058, 0.80352895],
[0.51505805, 0.32039379, 0.0786534, 0.02657018, 0.48488813]])
heights = np.array([ 22.110203, 65.49054, 110.321526, 156.54034, 166.59094])
local_z = np.load('local_z.npy')
print('difference in heights', local_z - heights)
def func(z, a0, a1):
return a0 + a1*z
for v in values:
popt_non_working = curve_fit(func, local_z, v)[0]
print('not working: ', popt_non_working)
popt_working = curve_fit(func, heights, v)[0]
print('working: ', popt_working)
我使用 Python 2.7.6、numpy 1.14.1 和 scipy 1.0.0 的输出:
$ python auxiliary/SUBLIME_fit_test.py
('difference in heights', array([-2.10693358e-07, -4.49218746e-07, -4.26269537e-07, 4.23828126e-06, 2.38281251e-06]))
/home/csag5117/venv/local/lib/python2.7/site-packages/scipy/optimize/minpack.py:785: OptimizeWarning: Covariance of the parameters could not be estimated category=OptimizeWarning)
('not working: ', array([1., 1.]))
('working: ', array([1.35420488, 0.00325281]))
('not working: ', array([1., 1.]))
('working: ', array([ 0.38896878, -0.00714073]))
('not working: ', array([1., 1.]))
('working: ', array([ 1.06301278, -0.00363439]))
('not working: ', array([1., 1.]))
('working: ', array([ 0.73398503, -0.00250946]))
('not working: ', array([1., 1.]))
('working: ', array([ 0.442922 , -0.00151433]))
如您所见,我使用 heights
作为 x 值的版本按预期工作(返回拟合参数),而我使用存储的 `local_z' 的版本没有,甚至尽管两个阵列之间的差异非常小。我只显示多个 y 值以表明这不是百万分之一的故障,可以通过适当的起始值来修复。这也只是一个示例,我还有一个具有相同行为的更多数据点(24 个而不是 5 个)。
为了完整起见,我添加了代码块(当我关闭它时一切正常)。有趣的是,通过在 MWE 中使用 local_z[:-1]
遗漏 local_z
中的最后一个值(这是代码块添加的值)并没有修复问题。
zi_minus_dd -= 1
zf_long = np.append(out.zf, np.squeeze(data.zf[t])[z_mask_full[-1] + 1])
u_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.u, np.squeeze(data.u[t])[z_mask_full[-1] + 1]))
v_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.v, np.squeeze(data.v[t])[z_mask_full[-1] + 1]))
th_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.th, np.squeeze(data.th[t])[z_mask_full[-1] + 1]))
zh_long = np.append(out.zh, np.squeeze(data.zh[t])[z_mask_full[-1] + 1])
uw_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.uw_raw, np.squeeze(data.uw[t])[z_mask_full[-1] + 1]))
vw_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.vw_raw, np.squeeze(data.vw[t])[z_mask_full[-1] + 1]))
tke_zi = np.interp(zi_minus_dd, zf_long,
np.append(out.tke, np.squeeze(data.TKE[t])[z_mask_full[-1] + 1]))
out.zf = np.append(out.zf, zi_minus_dd)
out.u = np.append(out.u, u_zi)
out.v = np.append(out.v, u_zi)
out.th = np.append(out.th, u_zi)
out.zh = np.append(out.zh, zi_minus_dd)
out.uw_raw = np.append(out.uw_raw, u_zi)
out.vw_raw = np.append(out.vw_raw, u_zi)
out.tke = np.append(out.tke, u_zi)
out.zf
和out.zh
是后面做成local_z
的向量。整个代码相当大,并且还依赖于 netCDF 文件(上面代码段中的 data
)。我已经问过了here ,但那是工作代码。
我很困惑,不知道如何解决这个问题,甚至不知道如何继续调试。 copy vs deepcopy 或类似的问题是否可能存在问题?虽然我想知道如何通过存储数组将其传输到 MWE...
最佳答案
追踪这个很有趣。 :-)
这不是值,而是它们的类型。这是一个精度问题:有效的 heights
是 float64,无效的 local_z 只是 float32。
我们有
In [70]: heights
Out[70]: array([ 22.110203, 65.49054 , 110.321526, 156.54034 , 166.59094 ])
In [71]: heights.dtype
Out[71]: dtype('float64')
In [72]: curve_fit(func, heights, v)[0]
Out[72]: array([1.35420488, 0.00325281])
和
In [73]: local_z
Out[73]:
array([ 22.110205, 65.49054 , 110.321526, 156.54034 , 166.59094 ],
dtype=float32)
In [74]: curve_fit(func, local_z, v)[0]
C:\Python\lib\site-packages\scipy\optimize\minpack.py:794: OptimizeWarning: Covariance of the parameters could not be estimated
category=OptimizeWarning)
Out[74]: array([1., 1.])
但如果我们愿意,我们可以让 local_z 工作:
In [75]: curve_fit(func, local_z.astype(np.float64), v)[0]
Out[75]: array([1.35420488, 0.00325281])
或高度失败:
In [76]: curve_fit(func, heights.astype(np.float32), v)[0]
C:\Python\lib\site-packages\scipy\optimize\minpack.py:794: OptimizeWarning: Covariance of the parameters could not be estimated
category=OptimizeWarning)
Out[76]: array([1., 1.])
关于python - Scipy curve_fit 仅针对非常特定的 x 值静默失败,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/52390002/
我仅在 WIN7 PC 上收到此通知,仅使用 IE。 Firefox 总是很好,旧版 Windows 上的 IE 似乎也不错。这让我大吃一惊,我不知道为什么 IE 认为 SSL 证书有问题。有没有人以
概述 对于我产品的新版本 v1.9.0,我创建了一个新的 MSI 安装程序。该应用程序的先前版本是 v1.7.0。 卸载旧版本然后安装新版本工作正常。 但是当我尝试使用 v1.9.0 安装程序更新旧版
该网站有一个全高图像启动。更多内容位于首屏下方,图像底部有一个“滚动”元素,以提示用户发现其余内容。单击后,我成功地使网站向下滚动 300 像素。然而,我想顺利地做到这一点。这是我当前的代码: w
var i = 0; function Myfunc() { var newdiv = document.createElement('div'); var el = document
这纯粹是为了学习目的;我知道 CSS 将是这种情况下的首选方法。 我知道在 JavaScript 中,您可以使用内联事件处理将鼠标悬停在图像上,如下所示: 我知道您可以在您的站点中安装 jQuery
我只想从curl请求中获取 header curl -I www.google.com 一切都很棒。现在我想这样做,但也传递发布数据: curl -I -d'test=test' www.google
以下代码旨在更改一个字段的颜色: Untitled Document var bkColor =
我正在使用 grep 递归搜索目录,并使用以下参数希望只返回第一个匹配项。不幸的是,它返回了不止一个——事实上,我上次查看时返回了两个。似乎我有太多的争论,尤其是没有得到想要的结果。 :-/ # gr
我只想搜索当前目录中的所有文件。我试过这个 grep foo * 但我收到此错误 grep: bar: Is a directory 我也尝试过这个 grep -r foo 但这也在搜索子目录。 最佳
我正在构建一个销售点应用程序,我想打印一张收据。问题是我使用的打印机无法打印纯文本的任何图形,我在 javafx 中只能找到使用 Print API 打印节点或使用像 jasper 这样都包含图形的报
是否有任何操作系统在完全加载时仅提供用于控制台应用程序执行的 java 环境?理想情况下,它会在加载时自动启动程序 最佳答案 这是一个名称为:JavaOS 的东西 从我的角度来看,更好的方法是安装一个
在工作中,我们有一个每晚执行 mysql 数据转储的脚本。对于开发,我们通常需要使用来自最近转储的数据。一段时间以来,我们一直每天都进行数据库还原,但现在我们已经到了每天还原花费近一个小时的地步。有没
我的移动模式菜单有问题。 onClick 它淡出。我想保留此设置,但我不希望它在单击下拉部分时淡出。这是链接:http://jsfiddle.net/zLLzrs6b/3/感谢您的帮助! html:
经过大量研究和反复试验,我谦虚地向各位 CSS 专家寻求帮助。这就是我需要的: 我有两张图片:titlelogo 和 newlogo。 在全屏模式下,newlogo 需要在左边,titlelogo 在
这个问题在这里已经有了答案: Exclusive CSS selector (3 个答案) 关闭 3 年前。 我的文档结构如下: ... ... something something someth
我有一个具有以下要求的表: 所有列的宽度必须可变 所有列的宽度不得超过必要的宽度 所有单元格必须保留空白(white-space:pre/pre-wrap) 当(且仅当)超过最大定义宽度 (1000p
我正在寻找一个正则表达式来仅匹配具有特殊 字符且大小为4+ 的数字 字符串。我对此处发布的问题做了一些评论: 测试网站: http://regexlib.com/RETester.aspx 1- re
我正在为我的元素开发一个纯 CSS 灯箱解决方案。我用谷歌搜索了它,但到目前为止只找到了部分解决方案。 我正在寻找这些功能: 显示任意宽任意高的内容(无固定高/宽) 垂直居中和水平居中 如果内容宽度和
出于各种原因,我目前正在尝试使用 HTML/CSS 创建网格布局(我知道 Bootstrap 等,但在这种情况下没有选择,而且我无法添加标记元素)。 我有以下代码(容器 div,每次都有一个带有 ul
有没有办法使用String.format()格式化 double 以仅获取小数? System.out.println(String.format("%.2f", 1.23456d)); 正如预期的那
我是一名优秀的程序员,十分优秀!