特殊函数导数近似值计算：Airy函数Aix导数实现与分析.pdfVIP

s–特殊函数的近似值s17ajc

nagairyaideriv(s17ajc)

1.目的

nagairyaideriv(s17ajc)返回Airy函数Ai(x)的导数值。

2.规范

#includenag.h#

includenags.h

___

doublenagairyaideriv(doublex,NagError*fail)

3.描述

此函数计算Airy函数Ai(x)的导数的近似值。它是基于多个Chebyshev展开式。

对于大的负参数，由于无法以任何精度计算振荡函数的结果，因此函数评估将失败。这发生在x

√

(π/)4/7，其中是机器精度。

0

对于大的正参数，由于Ai以本质上指数方式衰减，存在下溢的风险，因此函数必须失败。

4.参数

x输入：函数的参数x。failNAG错误参数，参见NAGC库的基本介绍。

5.错误指示和警告

NEREALARGGT

输入时，x大于hvaluei：x=hvaluei。x太大且为正数。

函数返回零。

NEREALARGLT

输入时，x小于hvaluei：x=hvaluei。x值太大且

为负数。函数返回零。

6.进一步评论

6.1.精度

对于负参数，该函数是振荡的，因此绝对误差是适当的度量。在正区域内，该函数本质上具

有指数特征，这里需要相对误差。绝对误差，E，和相对误差，，原则上与参数的相对误差，δ，

220

通过E|xAi(x)|δ,|xAi(x)/Ai(x)|δ相关。

实际上，近似等式是所能期望的最佳结果。当δ、或E处于机器精度的量级时，结果中的误差会稍

大一些。

对于小的x，无论是正还是负，这些误差都会被函数强烈衰减，因此将大致被机器精度所限制。

√

对于中等到较大的负数x，误差像函数一样是振荡的；然而，误差的幅度像|x|7/4/π一样增长。

√

因此，如果|x|7/4π/δ，则无法准确计算该函数。

√

3

对于大的正数x，相对误差放大相当严重，/δx。

Howe然而，非常大的参数由于存在下溢的风险而不可能。因此在实际中e

误差放大受到限制。

s–ApproximationsofSpecialFunctionss17ajc

nagairyaideriv(s17ajc)

1.Purpose

nagairyaideriv(s17ajc)returnsavalueofthederivativeoftheAiryfunctionAi(x).

2.Specification

#includenag.h

#i