对于嘈杂或结构精细的数据,有没有比中点法则更好的正交?
这个长问题的仅前两个部分是必不可少的。其他仅用于说明。 背景 较高阶的正交函数,例如更高阶的牛顿-科茨复合材料,高斯-莱根特雷和罗姆贝格复合材料,似乎主要用于那些可以对函数进行精细采样而不能进行分析集成的情况。但是,对于具有比采样间隔(例如,请参阅附录A)或测量噪声更精细的结构的函数,它们无法与简单的方法(例如中点或梯形规则)(请参见附录B进行演示)竞争。 这有点直观,例如,复合Simpson规则通过分配较低的权重实质上“丢弃”了四分之一的信息。这样的正交对于足够无聊的功能更好的唯一原因是,正确处理边界效应要大于丢弃信息的效应。从另一角度看,从我的直觉上很清楚,对于具有精细结构或噪声的函数,远离积分域边界的样本必须几乎等距且权重几乎相同(对于大量样本而言) )。另一方面,更好地处理边界效果(与中点方法相比)可能会受益于此类函数的正交。 题 假设我希望在数值上整合嘈杂或精细结构的一维数据。 采样点的数量是固定的(由于功能评估成本很高),但是我可以自由放置它们。但是,我(或方法)不能交互地放置采样点,即基于其他采样点的结果。我也不事先知道潜在的问题区域。因此,像Gauß–Legendre(非等距采样点)之类的东西就可以了。自适应正交不是必需的,因为它需要交互式放置采样点。 在这种情况下,是否建议使用任何超出中点方法的方法? 或者:是否有任何证据证明在这种情况下中点方法是最好的? 更笼统地说:关于此问题是否已有任何工作? 附录A:精细结构函数的特定示例 我想为估算: 其中和。一个典型的函数如下所示:∫10f(t)dt∫01f(t)dt\int_0^1f(t)\, \mathrm{d}tf(t)=∑i=1ksin(ωit−φi)ωi,f(t)=∑i=1ksin(ωit−φi)ωi, f(t) = \sum_{i=1}^{k} \frac{\sin(ω_i t-φ_i)}{ω_i},φi∈[0,2π]φi∈[0,2π]φ_i∈ [0,2π]logωi∈[1,1000]logωi∈[1,1000]\log{ω_i} ∈ [1,1000] 我为以下属性选择了此功能: 可以对控制结果进行分析集成。 它在某种程度上具有良好的结构,使得无法使用我正在使用的样本数量()捕获全部。<102<102<10^2 它并不以其精细的结构为主导。 附录B:基准 为了完整起见,这是Python中的基准测试: import numpy as np from numpy.random import uniform from scipy.integrate import simps, trapz, romb, fixed_quad begin = 0 end = 1 …