- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
借助 SciPy 中的 curve_fit
函数,我能够确定代表下图中所示曲线的系数。
def func2(t, tau):
return np.exp(-t / tau)
t2 = np.linspace(0, 4, 50)
y2 = func2(t2, 1.2)
y2_noise = 0.2 * np.random.normal(size=t2.size)
y2_curve_noise = y2 + y2_noise
popt2, pcov2 = curve_fit(func2, t2, y2_curve_noise)
tau2, = popt2
y2_fit = func2(t2, tau2)
我想用类似的函数来表示一些数据点。但是,我无法使用这种方法来拟合数据点,如下所示。
def func4(t, a, tau, c):
return a * np.exp(-t / tau) + c
t4 = np.array([15445.1, 15445.6, 15446.1, 15446.6, 15447.1, 15447.6, 15448.1,
15448.6, 15449.1, 15449.6, 15450.1, 15450.6, 15451.1, 15451.6,
15452.1, 15452.6, 15453.1, 15453.6, 15454.1, 15454.6, 15455.1,
15455.6, 15456.1, 15456.6, 15457.1, 15457.6, 15458.1, 15458.6,
15459.1, 15459.6, 15460.1, 15460.6, 15461.1, 15461.6, 15462.1,
15462.6, 15463.1, 15463.6, 15464.1, 15464.6, 15465.1, 15465.6,
15466.1, 15466.6, 15467.1, 15467.6, 15468.1, 15468.6, 15469.1,
15469.6, 15470.1, 15470.6, 15471.1, 15471.6, 15472.1, 15472.6,
15473.1, 15473.6, 15474.1])
y4 = np.array([4.129, 4.125, 4.123, 4.121, 4.119, 4.118, 4.116, 4.115, 4.114,
4.113, 4.112, 4.11, 4.11, 4.109, 4.108, 4.108, 4.106, 4.105,
4.105, 4.104, 4.103, 4.102, 4.102, 4.101, 4.1, 4.1, 4.099,
4.098, 4.098, 4.097, 4.097, 4.096, 4.095, 4.095, 4.094, 4.094,
4.093, 4.092, 4.092, 4.091, 4.091, 4.09, 4.09, 4.089, 4.089,
4.088, 4.088, 4.087, 4.087, 4.086, 4.086, 4.085, 4.085, 4.084,
4.084, 4.084, 4.083, 4.083, 4.082])
popt4, pcov4 = curve_fit(func4, t4, y4, p0=(4.129, 1.2, 4.082))
a4, tau4, c4 = popt4
y4_fit = func4(t4, a4, tau4, c4)
如何在 SciPy 中应用 curve_fit
来拟合数据点?或者我应该使用不同的曲线拟合方法吗?我也不确定初始猜测使用什么值 p0
。我只是根据数据选择了一些数字,但显然这对拟合没有帮助。
最佳答案
问题是 exp(-15000)
必须用大得离谱的 a
值来平衡,问题变得非常严重,所以优化例程失败。
规范化 t
使它们从 0 变为 1 有助于解决缩放问题。合理的初始猜测可以是:tau 为 1,c 为最小 y 值,a 为最大和最小 y 值之差。
t4_norm = (t4 - t4[0])/(t4[-1] - t4[0]) # normalized
c_0 = y4[-1]
tau_0 = 1
a_0 = (y4[0] - y4[-1])
popt4, pcov4 = curve_fit(func4, t4_norm, y4, p0=(a_0, tau_0, c_0))
a4, tau4, c4 = popt4
y4_fit = func4(t4_norm, a4, tau4, c4)
plt.plot(t4, y4, 'r.')
plt.plot(t4, y4_fit, 'b')
plt.show()
找到参数后,可以根据原来的t重新计算。的确,目前得到的曲线是
y = a*exp(- (t4 - t4[0])/(t4[-1] - t4[0]) / tau) + c
可以重写为
y = a*exp(t4[0]/(t4[-1] - t4[0]) / tau) * exp(-t4/(t4[-1] - t4[0]) / tau) + c
这意味着原始变量的参数是
a_orig = a*exp(t4[0]/(t4[-1] - t4[0]) / tau)
tau_orig = (t4[-1] - t4[0]) * tau
c_orig = c
关于python - 使用给定数据点曲线拟合 Python 中的指数衰减函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/49565152/
在下面的代码中,我得到一个 uninitialized value警告,但仅限于第二个 given/when例子。为什么是这样? #!/usr/bin/env perl use warnings; u
整个“开关”功能是否已成为实验性的?在没有 Perl 的 future 版本破坏我的代码的情况下,我可以依赖其中的某些部分吗?一般来说,将稳定功能更改为实验性的政策是什么? 背景use feature
有没有办法在一个条件语句中写出如下语句? a和b不能同时等于5。 (a可以是5,b可以是5,但是a AND b不能是5) 最佳答案 正如克里斯指出的那样,您要查找的是逻辑异或,相当于逻辑不等于 !=:
我正在寻找一种算法来找到给定 n 条线段的所有交点。以下是来自 http://jeffe.cs.illinois.edu/teaching/373/notes/x06-sweepline.pdf 的伪
数组中有 N 个元素。我可以选择第一项最多 N 次,第二项最多选择 N-1 次,依此类推。 我有 K 个 token 要使用并且需要使用它们以便我可以拥有最大数量的项目。 arr = [3, 4, 8
我正在尝试修复法语文本中的语法性别,想知道是否有办法从某个词条中获取所有单词的列表,以及是否可以在此类列表中进行查找? 最佳答案 尝试: import spacy lemma_lookup = spa
我正在为 Win32 编写一个简单的自动化测试应用程序。它作为一个单独的进程运行,并通过 Windows API 访问目标应用程序。我可以阅读窗口层次结构,查找标签和文本框,并通过发送/发布消息等来单
在 nodeJs 中使用 Sequelize 时,我从 Sequelize 收到此错误,如下所示: { [SequelizeUniqueConstraintError: Validation erro
本文https://arxiv.org/pdf/1703.10757.pdf使用回归激活映射 (RAM) - 而不是类激活映射 (CAM) 来解决问题。有几篇文章描述了如何实现 CAM。但是我找不到
我正在研究 Mach 动态链接器 dyld。这个问题适用于所有 Apple 平台,但很高兴得到特定于平台的答案;我正在使用 ObjC,但如果对你有用的话,我也很乐意翻译 Swift。 The rele
我有一个包含数千个 Instagram 用户 ID 的列表。我如何获得他们的 Instagram 用户名/句柄? 最佳答案 你必须使用这个 Instagram API: https://api.ins
我在下面的代码: def main(args: Array[String]) { val sparkConf = new SparkConf().setAppName("Spark-Hbase").s
我有一个表格,其中包含从 1 到 10 的数字。(从 D2 到 M2) 假设A1中有03/09/2019 并且在B1中有06/09/2019 并且在C1中有Hello 在A 列中,我有多个系列的单词,
我想在给定服务对应的 URI 的情况下检索服务的注释(特别是 @RolesAllowed )。这是一个例子: 服务: @GET @Path("/example") @RolesAllowed({ "B
我看到 OraclePreparedStatementexecuteQuery() 表现出序列化。也就是说,我想使用相同的连接对 Oracle 数据库同时运行两个查询。然而,OraclePrepare
import java.util.Scanner; public class GeometricSumFromK { public static int geometricSum(int k,
我创建了一个抽象基类Page,它说明了如何构建动态网页。我正在尝试想出一种基于作为 HttpServletRequest 传入的 GET 请求生成 Page 的好方法。例如... public cla
我的字符串是一条短信,采用以下两种格式之一: 潜在客户短信: 您已收到 1 条线索 标题:我的领导 潜在客户 ID:12345-2365 警报设置 ID:890 短信回复: 您已收到 1 条回复 标题
我在 python 中有以下代码: class CreateMap: def changeme(listOne, lisrTwo, listThree, listFour, listfive):
这是在 Hibernate 上运行的 JPA2。 我想检索相同实体类型的多个实例,给定它们的 ID。其中许多已经在持久性上下文和/或二级缓存中。 我尝试了几种方法,但似乎都有其缺点: 当我使用 ent
我是一名优秀的程序员,十分优秀!