Questions tagged «error-estimation»

4
数值误差的科学标准
在我的研究领域中,实验误差的规范被普遍接受,未能提供实验误差的出版物受到了强烈批评。同时,我经常发现提供的数值计算结果没有考虑任何数值误差,尽管(或者可能是因为)经常有可疑的数值方法正在起作用。我说的是由离散化和数值计算的有限精度等导致的误差。当然,这些误差估计值并不总是很容易获得,例如在流体动力学方程的情况下,但是我似乎常常是由于懒惰而导致的,我相信数值误差估计的规范应该与实验结果一样是标准的。因此,我的问题是:

5
为什么等距点的表现不好?
实验说明: 在拉格朗日插值中,精确的方程在NNN个点(多项式N−1N−1N - 1)处采样,并在101个点处插值。这里从2变化到64。每个时间,和制备误差图。可以看出,当在等距点对函数进行采样时,误差开始下降(直到小于约15时才会发生),然后误差随着进一步增加而上升。NNNL1L1L_1L2L2L_2L∞L∞L_\inftyNNNNNN 而如果初始采样是在Legendre-Gauss(LG)点(Legendre多项式的根)或Legendre-Gauss-Lobatto(LGL)点(Lobatto多项式的根)进行的,则误差降至机器级别,并且不会当进一步增加时增加。NNN 我的问题是 等距点到底会发生什么? 为什么多项式阶数增加会导致误差在某个点之后上升? 这是否还意味着如果我将等距点用于WENO / ENO重建(使用拉格朗日多项式),那么在平滑区域中,我会得到错误?(嗯,这些只是假设的问题(以我的理解),对于WENO方案,重建15或更高阶数的多项式确实是不合理的) 额外细节: 功能近似: X∈[-1,1]f(x)=cos(π2 x)f(x)=cos⁡(π2 x)f(x) = \cos(\frac{\pi}{2}~x),x∈[−1,1]x∈[−1,1]x \in [-1, 1] ñxxx分为等距点(以后称为LG)点。每次在101点插值该函数。NNN 结果: a)等距点(插值): N=65N=65N = 65 b)等距点(误差图,对数刻度): a)LG点(插值): N=65N=65N = 65 b)LG点(误差图,对数刻度):

4
带导数的数字正交
大多数求积分的数值方法都将被积数视为黑盒函数。如果我们有更多信息怎么办?特别是,如果了解被积物的前几个导数,我们可以从中得到什么好处?还有哪些其他信息可能有价值? 特别是对于导数:基本正交(矩形/梯形/辛普森规则)的误差估计密切相关。也许有一种方法可以预先选择采样分辨率,而不是依靠动态适应性? 我对单变量和多维情况都感兴趣。


1
与Adams-Bashforth算法相比,使用Adams-Moulton有哪些相对优势?
我正在计算两个空间维度和时间上两个耦合的PDE的系统。由于函数评估很昂贵,因此我想使用多步方法(使用Runge-Kutta 4-5初始化)。 使用五个先前的函数求值的Adams-Bashforth方法的全局误差为(在s = 5的情况下O(h5)O(h5)O(h^5)s=5s=5s=5在下面引用的Wikipedia文章下),并且每步需要一个函数求值(每个PDE)。 另一方面,Adams-Moulton方法每个步骤需要进行两次功能评估:一项用于预测步骤,另一项用于校正步骤。再一次,如果使用五个函数求值,则全局误差为。(s = 4O(h5)O(h5)O(h^5)s=4s=4s=4 Wikipedia文章中的) 那么,在Adams-Bashforth上使用Adams-Moulton的背后原因是什么?它具有相同数量级的错误,功能评估的次数是两倍。从直觉上讲,预测校正方法应该是有利的,但是有人可以对此进行定量解释吗? 参考:http : //en.wikipedia.org/wiki/Linear_multistep_method#Adams.E2.80.93Bashforth_methods

4
估计硬件错误概率
假设我在http://www.nersc.gov/users/computational-systems/edison/configuration上的100k内核上运行了4个小时的超级计算机计算,通过网络交换了大约4 PB的数据,并执行了大约4 TB的I /哦 计算都是整数,所以结果是对还是错(没有中间数值错误)。 假设代码正确,我想估计由于硬件故障而导致计算错误的可能性。有什么好的方法可以解决这个问题?是否有足够的资源可以估算出所需的数字?

2
有限元
w ^1 ,∞w ^1个,∞W^{1,\infty}∥ ü′H- 你′∥∞‖üH′-ü′‖∞\|u'_h - u'\|_\infty (来自MathOverflow的交叉发布,在这里几乎没有兴趣,但是也许在这里我可以找到更多具有FEM背景的人员。)

2
线性PDE的这种简单误差估计如何?
令为的凸多边形有界Lipschitz域,令。- [R 2 ˚F ∈ 大号2(Ω )ΩΩ\Omega[R2R2\mathbb R^2F∈ 大号2(Ω )f∈L2(Ω)f \in L^2(\Omega) 然后的狄利克雷问题的解决方案在,上具有独特的解决方案并且被适定的,即对于某一常数我们有。Δ û = ˚FΔu=f\Delta u = fΩΩ\Omega跟踪u = 0trace⁡u=0\operatorname{trace} u = 0∂Ω∂Ω\partial\OmegaH2H2H^2CCC∥ ü ∥H2≤ ç∥ ˚F∥大号2‖u‖H2≤C‖f‖L2\|u\|_{H^2} \leq C \|f\|_{L^2} 对于某些有限元逼近üHuhu_h,例如在均匀网格上具有节点元素的情况下,我们得到误差估计 ∥ ü - üH∥H1个≤ ç^ h ∥ ü ∥H2‖u−uh‖H1≤Ch‖u‖H2\| u - u_h \|_{H^1} \leq C h \| u …
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.