シグモイド形状の特殊関数
数学において、 誤差関数 ( ガウス誤差関数とも呼ばれる)は、 erf と表記され 、次のように定義される関数である 。 [1]
e
r
ふ
:
C
→
C
{\displaystyle \mathrm {erf} :\mathbb {C} \to \mathbb {C} }
エルフ
ず
=
2
π
∫
0
ず
e
−
t
2
d
t
。
{\displaystyle \operatorname {erf} z={\frac {2}{\sqrt {\pi }}}\int _{0}^{z}e^{-t^{2}}\,\mathrm {d} t.}
ここでの積分は複素 線積分 であり、複素平面全体で 正則で ある ため経路に依存しません 。多くのアプリケーションでは、関数の引数は実数であり、その場合、関数の値も実数です。
exp
(
−
t
2
)
{\displaystyle \exp(-t^{2})}
C
{\displaystyle \mathbb {C} }
いくつかの古いテキスト [2] では、
誤差関数は の因子なしで定義されています 。この 非基本的な積分は、 確率 、 統計 、 偏微分方程式 でよく見られる シグモイド 関数です 。
2
/
π
{\displaystyle 2/{\sqrt {\pi }}}
統計学では、 x の非負の実数値に対して 、誤差関数は次のように解釈されます。 平均 0 および 標準偏差 で 正規分布する 実数 ランダム変数 Y に対して、 erf x はY が範囲 [− x 、 x ] に含まれる確率です 。
1
/
2
{\displaystyle 1/{\sqrt {2}}}
2つの密接に関連した関数は 相補誤差関数 であり、次のように定義されます。
e
r
f
c
:
C
→
C
{\displaystyle \mathrm {erfc} :\mathbb {C} \to \mathbb {C} }
erfc
z
=
1
−
erf
z
,
{\displaystyle \operatorname {erfc} z=1-\operatorname {erf} z,}
そして 虚数誤差関数は 次のように定義される。
e
r
f
i
:
C
→
C
{\displaystyle \mathrm {erfi} :\mathbb {C} \to \mathbb {C} }
erfi
z
=
−
i
erf
i
z
,
{\displaystyle \operatorname {erfi} z=-i\operatorname {erf} iz,}
ここで、 i は 虚数単位 です 。
名前
「誤差関数」という名前とその略語 erf は、1871 年に JWL Glaisher によって、 「確率論、特に 誤差 理論」との関連から提案されました。 [3] 誤差関数補数についても、同年に Glaisher は別の出版物で論じました。 [4]誤差 密度 が
( 正規分布 )で与えられる
誤差の「容易さの法則」について、Glaisher は p と q の 間にある誤差の確率を 次のように計算します。
f
(
x
)
=
(
c
π
)
1
/
2
e
−
c
x
2
{\displaystyle f(x)=\left({\frac {c}{\pi }}\right)^{1/2}e^{-cx^{2}}}
(
c
π
)
1
2
∫
p
q
e
−
c
x
2
d
x
=
1
2
(
erf
(
q
c
)
−
erf
(
p
c
)
)
.
{\displaystyle \left({\frac {c}{\pi }}\right)^{\frac {1}{2}}\int _{p}^{q}e^{-cx^{2}}\,\mathrm {d} x={\tfrac {1}{2}}\left(\operatorname {erf} \left(q{\sqrt {c}}\right)-\operatorname {erf} \left(p{\sqrt {c}}\right)\right).}
Mathematica 13.1 関数 ComplexPlot3D で作成した、-2-2i から 2+2i までの複素平面における誤差関数 Erf(z) のカラー プロット
アプリケーション
一連の測定結果が 標準偏差 σ と 期待値 0の 正規分布 で記述される場合、 erf( 1つの / σ√2 )は、 a が 正の場合に、 単一の測定の誤差が − a と + a の間にある確率です。これは、たとえば、 デジタル通信システムの
ビット誤り率を決定するのに役立ちます。
誤差関数と相補誤差関数は、例えば、 境界条件が ヘヴィサイドのステップ関数 によって与えられた ときの 熱方程式 の解に発生します。
誤差関数とその近似値は、 高い確率で 成立する結果や低い確率で成立する結果を推定するために使用できます。ランダム変数 X ~ Norm[ μ , σ ] (平均 μ と標準偏差 σ の正規分布 ) と定数 L > μ が与えられている場合、置換積分によって次のように示されます。
Pr
[
X
≤
L
]
=
1
2
+
1
2
erf
L
−
μ
2
σ
≈
A
exp
(
−
B
(
L
−
μ
σ
)
2
)
{\displaystyle {\begin{aligned}\Pr[X\leq L]&={\frac {1}{2}}+{\frac {1}{2}}\operatorname {erf} {\frac {L-\mu }{{\sqrt {2}}\sigma }}\\&\approx A\exp \left(-B\left({\frac {L-\mu }{\sigma }}\right)^{2}\right)\end{aligned}}}
ここで、 A と Bは 特定の数値定数です。L が平均から十分に離れている場合、具体的には μ − L ≥ σ √ ln k の 場合、次のようになります。
Pr
[
X
≤
L
]
≤
A
exp
(
−
B
ln
k
)
=
A
k
B
{\displaystyle \Pr[X\leq L]\leq A\exp(-B\ln {k})={\frac {A}{k^{B}}}}
したがって、 k → ∞ のにつれて確率は 0 になります 。
X が区間 [ L a , L b ] にある確率は 次のように導出できる。
Pr
[
L
a
≤
X
≤
L
b
]
=
∫
L
a
L
b
1
2
π
σ
exp
(
−
(
x
−
μ
)
2
2
σ
2
)
d
x
=
1
2
(
erf
L
b
−
μ
2
σ
−
erf
L
a
−
μ
2
σ
)
.
{\displaystyle {\begin{aligned}\Pr[L_{a}\leq X\leq L_{b}]&=\int _{L_{a}}^{L_{b}}{\frac {1}{{\sqrt {2\pi }}\sigma }}\exp \left(-{\frac {(x-\mu )^{2}}{2\sigma ^{2}}}\right)\,\mathrm {d} x\\&={\frac {1}{2}}\left(\operatorname {erf} {\frac {L_{b}-\mu }{{\sqrt {2}}\sigma }}-\operatorname {erf} {\frac {L_{a}-\mu }{{\sqrt {2}}\sigma }}\right).\end{aligned}}}
プロパティ
erf (− z ) = −erf z という特性は 、誤差関数が 奇関数であることを意味します。これは、被積分関数 e − t 2 が偶関数 である という事実から直接生じます (原点でゼロである偶関数の原始関数は奇関数であり、その逆も同様です)。
誤差関数は、任意の 複素数 z に対して、実数を実数に変換する 完全な関数 であるため、
z
は z の 複素共役 です 。
erf
z
¯
=
erf
z
¯
{\displaystyle \operatorname {erf} {\overline {z}}={\overline {\operatorname {erf} z}}}
積分関数 f = exp(− z 2 ) と f = erf z は、右の図では複素 z 平面上に 領域の色分け とともに示されています。
+∞ での誤差関数は ちょうど 1 です ( ガウス積分を 参照)。実軸では、 erf z は z → +∞ で 1 に近づき、 z → −∞ で −1に近づきます。虚軸では ± i ∞ に近づきます 。
テイラー級数
誤差関数は 整関数 です。特異点は存在せず(無限大を除く)、 テイラー展開は 常に収束します。ただし、 x >> 1 の場合、主要項が相殺されるため、テイラー展開は実用的ではありません。
定義積分は、 基本関数 に関する 閉じた形式 では評価できません( リウヴィルの定理 を参照)が、被 積分関数 e − z 2 をその マクローリン級数 に展開し、項ごとに積分することで、誤差関数のマクローリン級数が得られます。
これは、すべての 複素数 z
に対して成り立ちます。分母項は、 OEIS のシーケンス A007680 です 。
erf
z
=
2
π
∑
n
=
0
∞
(
−
1
)
n
z
2
n
+
1
n
!
(
2
n
+
1
)
=
2
π
(
z
−
z
3
3
+
z
5
10
−
z
7
42
+
z
9
216
−
⋯
)
{\displaystyle {\begin{aligned}\operatorname {erf} z&={\frac {2}{\sqrt {\pi }}}\sum _{n=0}^{\infty }{\frac {(-1)^{n}z^{2n+1}}{n!(2n+1)}}\\[6pt]&={\frac {2}{\sqrt {\pi }}}\left(z-{\frac {z^{3}}{3}}+{\frac {z^{5}}{10}}-{\frac {z^{7}}{42}}+{\frac {z^{9}}{216}}-\cdots \right)\end{aligned}}}
上記の級数の反復計算には、次のような代替定式化が有用である可能性が
ある 。
erf
z
=
2
π
∑
n
=
0
∞
(
z
∏
k
=
1
n
−
(
2
k
−
1
)
z
2
k
(
2
k
+
1
)
)
=
2
π
∑
n
=
0
∞
z
2
n
+
1
∏
k
=
1
n
−
z
2
k
{\displaystyle {\begin{aligned}\operatorname {erf} z&={\frac {2}{\sqrt {\pi }}}\sum _{n=0}^{\infty }\left(z\prod _{k=1}^{n}{\frac {-(2k-1)z^{2}}{k(2k+1)}}\right)\\[6pt]&={\frac {2}{\sqrt {\pi }}}\sum _{n=0}^{\infty }{\frac {z}{2n+1}}\prod _{k=1}^{n}{\frac {-z^{2}}{k}}\end{aligned}}}
−(2 k − 1) z 2 / (2k + 1 ) は、 k番目の項を ( k + 1) 番目の項に 変換するための乗数を表します ( z を 最初の項と見なします)。
虚数誤差関数には非常によく似たマクローリン級数があり、これ
はすべての 複素数 z
に対して成り立ちます 。
erfi
z
=
2
π
∑
n
=
0
∞
z
2
n
+
1
n
!
(
2
n
+
1
)
=
2
π
(
z
+
z
3
3
+
z
5
10
+
z
7
42
+
z
9
216
+
⋯
)
{\displaystyle {\begin{aligned}\operatorname {erfi} z&={\frac {2}{\sqrt {\pi }}}\sum _{n=0}^{\infty }{\frac {z^{2n+1}}{n!(2n+1)}}\\[6pt]&={\frac {2}{\sqrt {\pi }}}\left(z+{\frac {z^{3}}{3}}+{\frac {z^{5}}{10}}+{\frac {z^{7}}{42}}+{\frac {z^{9}}{216}}+\cdots \right)\end{aligned}}}
微分と積分
誤差関数の導関数は定義から直ちに導かれる。
このことから、虚数誤差関数の導関数も直ちに導かれる。
誤差関数の 原始微分
は 部分積分 で得られ 、次のようになる
。虚数誤差関数の原始
微分も部分積分で得られ、次のようになる。
高次の導関数は次のように与えられる
。
ここで H は物理学者の エルミート多項式 である。 [5]
d
d
z
erf
z
=
2
π
e
−
z
2
.
{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}\operatorname {erf} z={\frac {2}{\sqrt {\pi }}}e^{-z^{2}}.}
d
d
z
erfi
z
=
2
π
e
z
2
.
{\displaystyle {\frac {d}{dz}}\operatorname {erfi} z={\frac {2}{\sqrt {\pi }}}e^{z^{2}}.}
z
erf
z
+
e
−
z
2
π
.
{\displaystyle z\operatorname {erf} z+{\frac {e^{-z^{2}}}{\sqrt {\pi }}}.}
z
erfi
z
−
e
z
2
π
.
{\displaystyle z\operatorname {erfi} z-{\frac {e^{z^{2}}}{\sqrt {\pi }}}.}
erf
(
k
)
z
=
2
(
−
1
)
k
−
1
π
H
k
−
1
(
z
)
e
−
z
2
=
2
π
d
k
−
1
d
z
k
−
1
(
e
−
z
2
)
,
k
=
1
,
2
,
…
{\displaystyle \operatorname {erf} ^{(k)}z={\frac {2(-1)^{k-1}}{\sqrt {\pi }}}{\mathit {H}}_{k-1}(z)e^{-z^{2}}={\frac {2}{\sqrt {\pi }}}{\frac {\mathrm {d} ^{k-1}}{\mathrm {d} z^{k-1}}}\left(e^{-z^{2}}\right),\qquad k=1,2,\dots }
ビュルマンシリーズ
テイラー展開よりも x のすべての実数値に対してより速く収束する展開 [6]は、 ハンス・ハインリッヒ・ビュルマン の定理 [7] を使用することによって得られる。
ここで sgnは 符号関数 である。最初の2つの係数のみを保持し、 c 1 = を 選択することによって、
erf
x
=
2
π
sgn
x
⋅
1
−
e
−
x
2
(
1
−
1
12
(
1
−
e
−
x
2
)
−
7
480
(
1
−
e
−
x
2
)
2
−
5
896
(
1
−
e
−
x
2
)
3
−
787
276480
(
1
−
e
−
x
2
)
4
−
⋯
)
=
2
π
sgn
x
⋅
1
−
e
−
x
2
(
π
2
+
∑
k
=
1
∞
c
k
e
−
k
x
2
)
.
{\displaystyle {\begin{aligned}\operatorname {erf} x&={\frac {2}{\sqrt {\pi }}}\operatorname {sgn} x\cdot {\sqrt {1-e^{-x^{2}}}}\left(1-{\frac {1}{12}}\left(1-e^{-x^{2}}\right)-{\frac {7}{480}}\left(1-e^{-x^{2}}\right)^{2}-{\frac {5}{896}}\left(1-e^{-x^{2}}\right)^{3}-{\frac {787}{276480}}\left(1-e^{-x^{2}}\right)^{4}-\cdots \right)\\[10pt]&={\frac {2}{\sqrt {\pi }}}\operatorname {sgn} x\cdot {\sqrt {1-e^{-x^{2}}}}\left({\frac {\sqrt {\pi }}{2}}+\sum _{k=1}^{\infty }c_{k}e^{-kx^{2}}\right).\end{aligned}}}
31 / 200 そして c 2 = − 341 / 8000 、得られた近似値は x = ±1.3796 で相対誤差が最大となり、0.0036127 未満になります。
erf
x
≈
2
π
sgn
x
⋅
1
−
e
−
x
2
(
π
2
+
31
200
e
−
x
2
−
341
8000
e
−
2
x
2
)
.
{\displaystyle \operatorname {erf} x\approx {\frac {2}{\sqrt {\pi }}}\operatorname {sgn} x\cdot {\sqrt {1-e^{-x^{2}}}}\left({\frac {\sqrt {\pi }}{2}}+{\frac {31}{200}}e^{-x^{2}}-{\frac {341}{8000}}e^{-2x^{2}}\right).}
逆関数
逆誤差関数
複素数 zが与えられたとき、 erf w = z を 満たす 唯一の 複素数 w は 存在しない ので、真の逆関数は多値となる。しかし、 −1 < x < 1の場合、 erf −1 x で表される唯一の実数 が 存在し 、
erf
(
erf
−
1
x
)
=
x
.
{\displaystyle \operatorname {erf} \left(\operatorname {erf} ^{-1}x\right)=x.}
逆誤差関数は 通常 、ドメイン (−1,1) で定義され、多くのコンピュータ代数システムではこのドメインに制限されています。しかし、マクローリン級数 [8] を使用して、複素平面のディスク | z | < 1 に拡張することができ
、 c 0 = 1 および
erf
−
1
z
=
∑
k
=
0
∞
c
k
2
k
+
1
(
π
2
z
)
2
k
+
1
,
{\displaystyle \operatorname {erf} ^{-1}z=\sum _{k=0}^{\infty }{\frac {c_{k}}{2k+1}}\left({\frac {\sqrt {\pi }}{2}}z\right)^{2k+1},}
c
k
=
∑
m
=
0
k
−
1
c
m
c
k
−
1
−
m
(
m
+
1
)
(
2
m
+
1
)
=
{
1
,
1
,
7
6
,
127
90
,
4369
2520
,
34807
16200
,
…
}
.
{\displaystyle {\begin{aligned}c_{k}&=\sum _{m=0}^{k-1}{\frac {c_{m}c_{k-1-m}}{(m+1)(2m+1)}}\\[1ex]&=\left\{1,1,{\frac {7}{6}},{\frac {127}{90}},{\frac {4369}{2520}},{\frac {34807}{16200}},\ldots \right\}.\end{aligned}}}
したがって、級数展開は次のようになります(分子と分母から共通因数が消去されています)。
(消去後、分子と分母の値はそれぞれ OEIS : A092676 と OEIS : A092677 になります。消去しない場合、分子項は OEIS : A002067の値になります。) ±∞ における誤差関数の値は ±1 に等しくなります 。
erf
−
1
z
=
π
2
(
z
+
π
12
z
3
+
7
π
2
480
z
5
+
127
π
3
40320
z
7
+
4369
π
4
5806080
z
9
+
34807
π
5
182476800
z
11
+
⋯
)
.
{\displaystyle \operatorname {erf} ^{-1}z={\frac {\sqrt {\pi }}{2}}\left(z+{\frac {\pi }{12}}z^{3}+{\frac {7\pi ^{2}}{480}}z^{5}+{\frac {127\pi ^{3}}{40320}}z^{7}+{\frac {4369\pi ^{4}}{5806080}}z^{9}+{\frac {34807\pi ^{5}}{182476800}}z^{11}+\cdots \right).}
| z | < 1 の場合 、 erf(erf −1 z ) = z となります。
逆相補誤差関数は 次 のように定義されます
。
実数 xに対して、 erfi(erfi −1 x ) = x を 満たす唯一の 実数 erfi −1 x が存在します 。 逆虚数誤差関数は erfi −1 x として定義されます 。 [9]
erfc
−
1
(
1
−
z
)
=
erf
−
1
z
.
{\displaystyle \operatorname {erfc} ^{-1}(1-z)=\operatorname {erf} ^{-1}z.}
任意の実数 x に対して、 ニュートン法を 使用して erfi −1 x を 計算できます。また、 −1 ≤ x ≤ 1 に対して 、次のマクローリン級数が収束します。
ここで、 c k は 上記のように定義されます。
erfi
−
1
z
=
∑
k
=
0
∞
(
−
1
)
k
c
k
2
k
+
1
(
π
2
z
)
2
k
+
1
,
{\displaystyle \operatorname {erfi} ^{-1}z=\sum _{k=0}^{\infty }{\frac {(-1)^{k}c_{k}}{2k+1}}\left({\frac {\sqrt {\pi }}{2}}z\right)^{2k+1},}
漸近展開
大きな実数x に対する相補誤差関数 (したがって誤差関数も) の 便利な 漸近展開は、 次の式です。
ここで、 (
2 n − 1)!!は、 (2 n − 1) の 二重階 乗で、これは、 (2 n − 1) までのすべての奇数の積です 。この級数は、すべての有限 xに対して発散し、漸近展開としての意味は、任意の整数 N ≥ 1 に対して、次の式が得られる
ことです。
ここで、剰余は、
帰納法、書き方
、および部分積分によって簡単に従うことができます。
erfc
x
=
e
−
x
2
x
π
(
1
+
∑
n
=
1
∞
(
−
1
)
n
1
⋅
3
⋅
5
⋯
(
2
n
−
1
)
(
2
x
2
)
n
)
=
e
−
x
2
x
π
∑
n
=
0
∞
(
−
1
)
n
(
2
n
−
1
)
!
!
(
2
x
2
)
n
,
{\displaystyle {\begin{aligned}\operatorname {erfc} x&={\frac {e^{-x^{2}}}{x{\sqrt {\pi }}}}\left(1+\sum _{n=1}^{\infty }(-1)^{n}{\frac {1\cdot 3\cdot 5\cdots (2n-1)}{\left(2x^{2}\right)^{n}}}\right)\\[6pt]&={\frac {e^{-x^{2}}}{x{\sqrt {\pi }}}}\sum _{n=0}^{\infty }(-1)^{n}{\frac {(2n-1)!!}{\left(2x^{2}\right)^{n}}},\end{aligned}}}
erfc
x
=
e
−
x
2
x
π
∑
n
=
0
N
−
1
(
−
1
)
n
(
2
n
−
1
)
!
!
(
2
x
2
)
n
+
R
N
(
x
)
{\displaystyle \operatorname {erfc} x={\frac {e^{-x^{2}}}{x{\sqrt {\pi }}}}\sum _{n=0}^{N-1}(-1)^{n}{\frac {(2n-1)!!}{\left(2x^{2}\right)^{n}}}+R_{N}(x)}
R
N
(
x
)
:=
(
−
1
)
N
π
2
1
−
2
N
(
2
N
)
!
N
!
∫
x
∞
t
−
2
N
e
−
t
2
d
t
,
{\displaystyle R_{N}(x):={\frac {(-1)^{N}}{\sqrt {\pi }}}2^{1-2N}{\frac {(2N)!}{N!}}\int _{x}^{\infty }t^{-2N}e^{-t^{2}}\,\mathrm {d} t,}
e
−
t
2
=
−
(
2
t
)
−
1
(
e
−
t
2
)
′
{\displaystyle e^{-t^{2}}=-(2t)^{-1}\left(e^{-t^{2}}\right)'}
ランダウ記法 における剰余項の漸近的挙動は、
x → ∞ です
。これは次のように求められます。
x
の値が十分に大きい場合、 erfc x の良好な近似値を得るには、この漸近展開の最初の数項のみが必要です(一方、 x の値がそれほど大きくない場合は 、上記の 0 でのテイラー展開により、非常に高速に収束します)。
R
N
(
x
)
=
O
(
x
−
(
1
+
2
N
)
e
−
x
2
)
{\displaystyle R_{N}(x)=O\left(x^{-(1+2N)}e^{-x^{2}}\right)}
R
N
(
x
)
∝
∫
x
∞
t
−
2
N
e
−
t
2
d
t
=
e
−
x
2
∫
0
∞
(
t
+
x
)
−
2
N
e
−
t
2
−
2
t
x
d
t
≤
e
−
x
2
∫
0
∞
x
−
2
N
e
−
2
t
x
d
t
∝
x
−
(
1
+
2
N
)
e
−
x
2
.
{\displaystyle R_{N}(x)\propto \int _{x}^{\infty }t^{-2N}e^{-t^{2}}\,\mathrm {d} t=e^{-x^{2}}\int _{0}^{\infty }(t+x)^{-2N}e^{-t^{2}-2tx}\,\mathrm {d} t\leq e^{-x^{2}}\int _{0}^{\infty }x^{-2N}e^{-2tx}\,\mathrm {d} t\propto x^{-(1+2N)}e^{-x^{2}}.}
連分数展開
相補誤差関数の連分数展開はラプラスによって発見された : [ 10 ] [ 11]
erfc
z
=
z
π
e
−
z
2
1
z
2
+
a
1
1
+
a
2
z
2
+
a
3
1
+
⋯
,
a
m
=
m
2
.
{\displaystyle \operatorname {erfc} z={\frac {z}{\sqrt {\pi }}}e^{-z^{2}}{\cfrac {1}{z^{2}+{\cfrac {a_{1}}{1+{\cfrac {a_{2}}{z^{2}+{\cfrac {a_{3}}{1+\dotsb }}}}}}}},\qquad a_{m}={\frac {m}{2}}.}
階乗級数
逆階乗級数:
Re( z2 )>0
で収束する 。ここで
zn は 上昇階乗 、 s ( n , k )は 第一種 符号付きスターリング数を表す 。 [ 12] [13] 二重階乗 を
含む無限和による表現も存在する :
erfc
z
=
e
−
z
2
π
z
∑
n
=
0
∞
(
−
1
)
n
Q
n
(
z
2
+
1
)
n
¯
=
e
−
z
2
π
z
[
1
−
1
2
1
(
z
2
+
1
)
+
1
4
1
(
z
2
+
1
)
(
z
2
+
2
)
−
⋯
]
{\displaystyle {\begin{aligned}\operatorname {erfc} z&={\frac {e^{-z^{2}}}{{\sqrt {\pi }}\,z}}\sum _{n=0}^{\infty }{\frac {\left(-1\right)^{n}Q_{n}}{{\left(z^{2}+1\right)}^{\bar {n}}}}\\[1ex]&={\frac {e^{-z^{2}}}{{\sqrt {\pi }}\,z}}\left[1-{\frac {1}{2}}{\frac {1}{(z^{2}+1)}}+{\frac {1}{4}}{\frac {1}{\left(z^{2}+1\right)\left(z^{2}+2\right)}}-\cdots \right]\end{aligned}}}
Q
n
=
def
1
Γ
(
1
2
)
∫
0
∞
τ
(
τ
−
1
)
⋯
(
τ
−
n
+
1
)
τ
−
1
2
e
−
τ
d
τ
=
∑
k
=
0
n
(
1
2
)
k
¯
s
(
n
,
k
)
,
{\displaystyle {\begin{aligned}Q_{n}&{\overset {\text{def}}{{}={}}}{\frac {1}{\Gamma {\left({\frac {1}{2}}\right)}}}\int _{0}^{\infty }\tau (\tau -1)\cdots (\tau -n+1)\tau ^{-{\frac {1}{2}}}e^{-\tau }\,d\tau \\[1ex]&=\sum _{k=0}^{n}\left({\frac {1}{2}}\right)^{\bar {k}}s(n,k),\end{aligned}}}
erf
z
=
2
π
∑
n
=
0
∞
(
−
2
)
n
(
2
n
−
1
)
!
!
(
2
n
+
1
)
!
z
2
n
+
1
{\displaystyle \operatorname {erf} z={\frac {2}{\sqrt {\pi }}}\sum _{n=0}^{\infty }{\frac {(-2)^{n}(2n-1)!!}{(2n+1)!}}z^{2n+1}}
数値近似
基本関数による近似
アブラモウィッツとステグンは、 精度の異なるいくつかの近似値(式7.1.25~28)を提示しています。これにより、特定のアプリケーションに適した最速の近似値を選択できます。精度が上がる順に並べると、次のようになります。
(最大誤差:
erf
x
≈
1
−
1
(
1
+
a
1
x
+
a
2
x
2
+
a
3
x
3
+
a
4
x
4
)
4
,
x
≥
0
{\displaystyle \operatorname {erf} x\approx 1-{\frac {1}{\left(1+a_{1}x+a_{2}x^{2}+a_{3}x^{3}+a_{4}x^{4}\right)^{4}}},\qquad x\geq 0}
5 × 10 −4 )
ここで、 a 1 = 0.278393 、 a 2 = 0.230389 、 a 3 = 0.000972 、 a 4 = 0.078108である。
erf
x
≈
1
−
(
a
1
t
+
a
2
t
2
+
a
3
t
3
)
e
−
x
2
,
t
=
1
1
+
p
x
,
x
≥
0
{\displaystyle \operatorname {erf} x\approx 1-\left(a_{1}t+a_{2}t^{2}+a_{3}t^{3}\right)e^{-x^{2}},\quad t={\frac {1}{1+px}},\qquad x\geq 0}
(最大誤差: 2.5 × 10 −5 )
ここで p = 0.47047 、 a 1 = 0.3480242 、 a 2 = −0.0958798 、 a 3 = 0.7478556である。
erf
x
≈
1
−
1
(
1
+
a
1
x
+
a
2
x
2
+
⋯
+
a
6
x
6
)
16
,
x
≥
0
{\displaystyle \operatorname {erf} x\approx 1-{\frac {1}{\left(1+a_{1}x+a_{2}x^{2}+\cdots +a_{6}x^{6}\right)^{16}}},\qquad x\geq 0}
(最大誤差: 3 × 10 −7 )
ここで、 a 1 = 0.0705230784 、 a 2 = 0.0422820123 、 a 3 = 0.0092705272 、 a 4 = 0.0001520143 、 a 5 = 0.0002765672 、 a 6 = 0.0000430638
erf
x
≈
1
−
(
a
1
t
+
a
2
t
2
+
⋯
+
a
5
t
5
)
e
−
x
2
,
t
=
1
1
+
p
x
{\displaystyle \operatorname {erf} x\approx 1-\left(a_{1}t+a_{2}t^{2}+\cdots +a_{5}t^{5}\right)e^{-x^{2}},\quad t={\frac {1}{1+px}}}
(最大誤差: 1.5 × 10 −7 )
ここで、 p = 0.3275911 、 a 1 = 0.254829592 、 a 2 = −0.284496736 、 a 3 = 1.421413741 、 a 4 = −1.453152027 、 a 5 = 1.061405429である。
これらの近似値はすべて x ≥ 0 の場合に有効です。負の xに対してこれらの近似値を使用するには、 erf x が 奇関数である という事実、つまり erf x = −erf(− x ) を使用します。
指数関数的境界と相補誤差関数の純粋指数近似は [14]で与えられる。
erfc
x
≤
1
2
e
−
2
x
2
+
1
2
e
−
x
2
≤
e
−
x
2
,
x
>
0
erfc
x
≈
1
6
e
−
x
2
+
1
2
e
−
4
3
x
2
,
x
>
0.
{\displaystyle {\begin{aligned}\operatorname {erfc} x&\leq {\frac {1}{2}}e^{-2x^{2}}+{\frac {1}{2}}e^{-x^{2}}\leq e^{-x^{2}},&\quad x&>0\\[1.5ex]\operatorname {erfc} x&\approx {\frac {1}{6}}e^{-x^{2}}+{\frac {1}{2}}e^{-{\frac {4}{3}}x^{2}},&\quad x&>0.\end{aligned}}}
上記はN の指数
関数の和に一般化されており [15]、 N に関して精度が増す ため、 erfc xは 2 Q̃ ( √2 x ) で正確に近似または制限される 。
特に、数値係数 {( a n 、 b n )}を解くため の体系的な方法論がある。
Q
~
(
x
)
=
∑
n
=
1
N
a
n
e
−
b
n
x
2
.
{\displaystyle {\tilde {Q}}(x)=\sum _{n=1}^{N}a_{n}e^{-b_{n}x^{2}}.}
1 です 密接に関連する Q関数の ミニマックス 近似または境界値を与える : Q ( x ) ≈ Q̃ ( x ) 、 Q ( x ) ≤ Q̃ ( x ) 、または Q ( x ) ≥ Q̃ ( x ) ( x ≥ 0 )。係数 {( a n 、 b n )} 1 です 指数近似の多くのバリエーションと N = 25 までの境界が包括的なデータセットとしてオープンアクセスに公開されています。 [16]
x ∈ [0,∞)
に対する相補誤差関数の厳密な近似は Karagiannidis & Lioumpas (2007) [17] によって与えられており 、適切なパラメータ { A , B } を
選択すると { A , B } = {1.98,1.135} となり
、すべての x ≥ 0 に対して良好な近似値が得られることが示されています 。特定のアプリケーションに合わせて精度を調整したり、式を厳密な境界に変換したりするために、代替係数も利用できます。 [18]
erfc
x
≈
(
1
−
e
−
A
x
)
e
−
x
2
B
π
x
.
{\displaystyle \operatorname {erfc} x\approx {\frac {\left(1-e^{-Ax}\right)e^{-x^{2}}}{B{\sqrt {\pi }}x}}.}
一項下限値は [19]
であり
、パラメータ βは 近似の望ましい区間における誤差を最小化するように選択することができる。
erfc
x
≥
2
e
π
β
−
1
β
e
−
β
x
2
,
x
≥
0
,
β
>
1
,
{\displaystyle \operatorname {erfc} x\geq {\sqrt {\frac {2e}{\pi }}}{\frac {\sqrt {\beta -1}}{\beta }}e^{-\beta x^{2}},\qquad x\geq 0,\quad \beta >1,}
セルゲイ・ウィニツキは「大域パデ近似」を用いて別の近似値を与えている: [20] [21] : 2–3
ここで、
これは0近傍と無限大近傍で非常に正確になるように設計されており、 相対誤差はすべての実数 x に対して0.00035未満である 。代替値 a ≈ 0.147 を使用すると、最大相対誤差は約0.00013に減少する。 [22]
erf
x
≈
sgn
x
⋅
1
−
exp
(
−
x
2
4
π
+
a
x
2
1
+
a
x
2
)
{\displaystyle \operatorname {erf} x\approx \operatorname {sgn} x\cdot {\sqrt {1-\exp \left(-x^{2}{\frac {{\frac {4}{\pi }}+ax^{2}}{1+ax^{2}}}\right)}}}
a
=
8
(
π
−
3
)
3
π
(
4
−
π
)
≈
0.140012.
{\displaystyle a={\frac {8(\pi -3)}{3\pi (4-\pi )}}\approx 0.140012.}
この近似値を反転すると、逆誤差関数の近似値が得られます。
erf
−
1
x
≈
sgn
x
⋅
(
2
π
a
+
ln
(
1
−
x
2
)
2
)
2
−
ln
(
1
−
x
2
)
a
−
(
2
π
a
+
ln
(
1
−
x
2
)
2
)
.
{\displaystyle \operatorname {erf} ^{-1}x\approx \operatorname {sgn} x\cdot {\sqrt {{\sqrt {\left({\frac {2}{\pi a}}+{\frac {\ln \left(1-x^{2}\right)}{2}}\right)^{2}-{\frac {\ln \left(1-x^{2}\right)}{a}}}}-\left({\frac {2}{\pi a}}+{\frac {\ln \left(1-x^{2}\right)}{2}}\right)}}.}
最大誤差が 1.2 × 10 −7 は 任意の実引数に対して次の式で表される:
[
23]
および
erf
x
=
{
1
−
τ
x
≥
0
τ
−
1
x
<
0
{\displaystyle \operatorname {erf} x={\begin{cases}1-\tau &x\geq 0\\\tau -1&x<0\end{cases}}}
τ
=
t
⋅
exp
(
−
x
2
−
1.26551223
+
1.00002368
t
+
0.37409196
t
2
+
0.09678418
t
3
−
0.18628806
t
4
+
0.27886807
t
5
−
1.13520398
t
6
+
1.48851587
t
7
−
0.82215223
t
8
+
0.17087277
t
9
)
{\displaystyle {\begin{aligned}\tau &=t\cdot \exp \left(-x^{2}-1.26551223+1.00002368t+0.37409196t^{2}+0.09678418t^{3}-0.18628806t^{4}\right.\\&\left.\qquad \qquad \qquad +0.27886807t^{5}-1.13520398t^{6}+1.48851587t^{7}-0.82215223t^{8}+0.17087277t^{9}\right)\end{aligned}}}
t
=
1
1
+
1
2
|
x
|
.
{\displaystyle t={\frac {1}{1+{\frac {1}{2}}|x|}}.}
絶対値で 最大相対誤差が 未満の の近似値は、
については [24] 、
については
erfc
{\displaystyle \operatorname {erfc} }
2
−
53
{\displaystyle 2^{-53}}
(
≈
1.1
×
10
−
16
)
{\displaystyle \left(\approx 1.1\times 10^{-16}\right)}
x
≥
0
{\displaystyle x\geq 0}
erfc
(
x
)
=
(
0.56418958354775629
x
+
2.06955023132914151
)
(
x
2
+
2.71078540045147805
x
+
5.80755613130301624
x
2
+
3.47954057099518960
x
+
12.06166887286239555
)
(
x
2
+
3.47469513777439592
x
+
12.07402036406381411
x
2
+
3.72068443960225092
x
+
8.44319781003968454
)
(
x
2
+
4.00561509202259545
x
+
9.30596659485887898
x
2
+
3.90225704029924078
x
+
6.36161630953880464
)
(
x
2
+
5.16722705817812584
x
+
9.12661617673673262
x
2
+
4.03296893109262491
x
+
5.13578530585681539
)
(
x
2
+
5.95908795446633271
x
+
9.19435612886969243
x
2
+
4.11240942957450885
x
+
4.48640329523408675
)
e
−
x
2
{\displaystyle {\begin{aligned}\operatorname {erfc} \left(x\right)&=\left({\frac {0.56418958354775629}{x+2.06955023132914151}}\right)\left({\frac {x^{2}+2.71078540045147805x+5.80755613130301624}{x^{2}+3.47954057099518960x+12.06166887286239555}}\right)\\&\left({\frac {x^{2}+3.47469513777439592x+12.07402036406381411}{x^{2}+3.72068443960225092x+8.44319781003968454}}\right)\left({\frac {x^{2}+4.00561509202259545x+9.30596659485887898}{x^{2}+3.90225704029924078x+6.36161630953880464}}\right)\\&\left({\frac {x^{2}+5.16722705817812584x+9.12661617673673262}{x^{2}+4.03296893109262491x+5.13578530585681539}}\right)\left({\frac {x^{2}+5.95908795446633271x+9.19435612886969243}{x^{2}+4.11240942957450885x+4.48640329523408675}}\right)e^{-x^{2}}\\\end{aligned}}}
x
<
0
{\displaystyle x<0}
erfc
(
x
)
=
2
−
erfc
(
−
x
)
{\displaystyle \operatorname {erfc} \left(x\right)=2-\operatorname {erfc} \left(-x\right)}
実数値引数の簡単な近似は、
絶対差を維持する 双曲線関数 によって
行うことができます 。
erf
(
x
)
≈
z
(
x
)
=
tanh
(
2
π
(
x
+
11
123
x
3
)
)
{\displaystyle \operatorname {erf} \left(x\right)\approx z(x)=\tanh \left({\frac {2}{\sqrt {\pi }}}\left(x+{\frac {11}{123}}x^{3}\right)\right)}
|
erf
(
x
)
−
z
(
x
)
|
<
0.000358
,
∀
x
{\displaystyle \left|\operatorname {erf} \left(x\right)-z(x)\right|<0.000358,\,\forall x}
値の表
相補誤差関数
相補誤差関数 ( erfc )は 次のように定義される。
Mathematica 13.1 関数 ComplexPlot3D で作成した、-2-2i から 2+2i までの複素平面における相補誤差関数 Erfc(z) のカラー プロット
erfc
x
=
1
−
erf
x
=
2
π
∫
x
∞
e
−
t
2
d
t
=
e
−
x
2
erfcx
x
,
{\displaystyle {\begin{aligned}\operatorname {erfc} x&=1-\operatorname {erf} x\\[5pt]&={\frac {2}{\sqrt {\pi }}}\int _{x}^{\infty }e^{-t^{2}}\,\mathrm {d} t\\[5pt]&=e^{-x^{2}}\operatorname {erfcx} x,\end{aligned}}}
これはまた 、 スケールされた相補誤差関数 erfcxも定義します [25] ( 算術アンダーフローを 避けるために erfc の代わりに使用できます [25] [26] )。 x ≥ 0 に対する erfc x の別の形式は 、発見者にちなんでクレイグの公式として知られています。 [27]この式は x
が正の値の場合にのみ有効ですが、 erfc x = 2 − erfc(− x ) と組み合わせて使用すると、負の値に対する erfc( x ) を取得できます 。 この形式の利点は、積分範囲が固定され有限であることです。 この式を 2 つの非負変数の和の erfcに拡張すると、次のようになります。 [28]
erfc
(
x
∣
x
≥
0
)
=
2
π
∫
0
π
2
exp
(
−
x
2
sin
2
θ
)
d
θ
.
{\displaystyle \operatorname {erfc} (x\mid x\geq 0)={\frac {2}{\pi }}\int _{0}^{\frac {\pi }{2}}\exp \left(-{\frac {x^{2}}{\sin ^{2}\theta }}\right)\,\mathrm {d} \theta .}
erfc
(
x
+
y
∣
x
,
y
≥
0
)
=
2
π
∫
0
π
2
exp
(
−
x
2
sin
2
θ
−
y
2
cos
2
θ
)
d
θ
.
{\displaystyle \operatorname {erfc} (x+y\mid x,y\geq 0)={\frac {2}{\pi }}\int _{0}^{\frac {\pi }{2}}\exp \left(-{\frac {x^{2}}{\sin ^{2}\theta }}-{\frac {y^{2}}{\cos ^{2}\theta }}\right)\,\mathrm {d} \theta .}
虚数誤差関数
虚数誤差関数 ( erfi と表記) は次のように定義される。
Mathematica 13.1 関数 ComplexPlot3D で作成した、-2-2i から 2+2i までの複素平面における虚数誤差関数 Erfi(z) のカラー プロット
erfi
x
=
−
i
erf
i
x
=
2
π
∫
0
x
e
t
2
d
t
=
2
π
e
x
2
D
(
x
)
,
{\displaystyle {\begin{aligned}\operatorname {erfi} x&=-i\operatorname {erf} ix\\[5pt]&={\frac {2}{\sqrt {\pi }}}\int _{0}^{x}e^{t^{2}}\,\mathrm {d} t\\[5pt]&={\frac {2}{\sqrt {\pi }}}e^{x^{2}}D(x),\end{aligned}}}
ここで D ( x ) は ドーソン関数( 算術オーバーフローを 避けるために erfi の代わりに使用できる [25] )。
「虚数誤差関数」という名前にもかかわらず、 xが実数の場合、 erfi x は 実数です 。
誤差関数が任意の 複素 引数 z に対して評価される場合、結果として得られる 複素誤差関数は通常、 Faddeeva関数 としてスケールされた形式で議論されます 。
w
(
z
)
=
e
−
z
2
erfc
(
−
i
z
)
=
erfcx
(
−
i
z
)
.
{\displaystyle w(z)=e^{-z^{2}}\operatorname {erfc} (-iz)=\operatorname {erfcx} (-iz).}
累積分布関数
誤差関数は、 標準 正規累積分布関数 (Φ と表記、 一部のソフトウェア言語では norm ( x )とも呼ば れる ) と本質的に同一であり、スケーリングと変換のみが異なる。実際、
複素平面上にプロットされた正規累積分布関数
Φ
(
x
)
=
1
2
π
∫
−
∞
x
e
−
t
2
2
d
t
=
1
2
(
1
+
erf
x
2
)
=
1
2
erfc
(
−
x
2
)
{\displaystyle {\begin{aligned}\Phi (x)&={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{x}e^{\tfrac {-t^{2}}{2}}\,\mathrm {d} t\\[6pt]&={\frac {1}{2}}\left(1+\operatorname {erf} {\frac {x}{\sqrt {2}}}\right)\\[6pt]&={\frac {1}{2}}\operatorname {erfc} \left(-{\frac {x}{\sqrt {2}}}\right)\end{aligned}}}
または erf と erfc 用に並べ替えます。
erf
(
x
)
=
2
Φ
(
x
2
)
−
1
erfc
(
x
)
=
2
Φ
(
−
x
2
)
=
2
(
1
−
Φ
(
x
2
)
)
.
{\displaystyle {\begin{aligned}\operatorname {erf} (x)&=2\Phi \left(x{\sqrt {2}}\right)-1\\[6pt]\operatorname {erfc} (x)&=2\Phi \left(-x{\sqrt {2}}\right)\\&=2\left(1-\Phi \left(x{\sqrt {2}}\right)\right).\end{aligned}}}
その結果、誤差関数は標準正規分布の裾の確率である Q関数 とも密接に関係している。Q関数は誤差関数を使って次のように表現できる。
Q
(
x
)
=
1
2
−
1
2
erf
x
2
=
1
2
erfc
x
2
.
{\displaystyle {\begin{aligned}Q(x)&={\frac {1}{2}}-{\frac {1}{2}}\operatorname {erf} {\frac {x}{\sqrt {2}}}\\&={\frac {1}{2}}\operatorname {erfc} {\frac {x}{\sqrt {2}}}.\end{aligned}}}
Φ の 逆 関数は 正規分位関数 または プロビット 関数として知られており 、逆誤差関数で次のように表される。
probit
(
p
)
=
Φ
−
1
(
p
)
=
2
erf
−
1
(
2
p
−
1
)
=
−
2
erfc
−
1
(
2
p
)
.
{\displaystyle \operatorname {probit} (p)=\Phi ^{-1}(p)={\sqrt {2}}\operatorname {erf} ^{-1}(2p-1)=-{\sqrt {2}}\operatorname {erfc} ^{-1}(2p).}
標準正規分布の CDF は確率と統計でより頻繁に使用され、誤差関数は数学の他の分野でより頻繁に使用されます。
誤差関数は ミッタク・レフラー関数 の特殊なケースであり、 合流型超幾何関数 (クンマー関数)として表現することもできます。
erf
x
=
2
x
π
M
(
1
2
,
3
2
,
−
x
2
)
.
{\displaystyle \operatorname {erf} x={\frac {2x}{\sqrt {\pi }}}M\left({\tfrac {1}{2}},{\tfrac {3}{2}},-x^{2}\right).}
これはフレネル積分 を使って簡単に表現できます 。 [ さらに説明が必要 ]
正規化されたガンマ関数 P と 不完全ガンマ関数 に関して 、
sgn x は 符号関数 です 。
erf
x
=
sgn
x
⋅
P
(
1
2
,
x
2
)
=
sgn
x
π
γ
(
1
2
,
x
2
)
.
{\displaystyle \operatorname {erf} x=\operatorname {sgn} x\cdot P\left({\tfrac {1}{2}},x^{2}\right)={\frac {\operatorname {sgn} x}{\sqrt {\pi }}}\gamma {\left({\tfrac {1}{2}},x^{2}\right)}.}
相補誤差関数の反復積分
相補誤差関数の反復積分は [29]で定義される。
i
n
erfc
z
=
∫
z
∞
i
n
−
1
erfc
ζ
d
ζ
i
0
erfc
z
=
erfc
z
i
1
erfc
z
=
ierfc
z
=
1
π
e
−
z
2
−
z
erfc
z
i
2
erfc
z
=
1
4
(
erfc
z
−
2
z
ierfc
z
)
{\displaystyle {\begin{aligned}i^{n}\!\operatorname {erfc} z&=\int _{z}^{\infty }i^{n-1}\!\operatorname {erfc} \zeta \,\mathrm {d} \zeta \\[6pt]i^{0}\!\operatorname {erfc} z&=\operatorname {erfc} z\\i^{1}\!\operatorname {erfc} z&=\operatorname {ierfc} z={\frac {1}{\sqrt {\pi }}}e^{-z^{2}}-z\operatorname {erfc} z\\i^{2}\!\operatorname {erfc} z&={\tfrac {1}{4}}\left(\operatorname {erfc} z-2z\operatorname {ierfc} z\right)\\\end{aligned}}}
一般的な漸化式は
2
n
⋅
i
n
erfc
z
=
i
n
−
2
erfc
z
−
2
z
⋅
i
n
−
1
erfc
z
{\displaystyle 2n\cdot i^{n}\!\operatorname {erfc} z=i^{n-2}\!\operatorname {erfc} z-2z\cdot i^{n-1}\!\operatorname {erfc} z}
これらは、対称性の性質に従う
べき級数を持ち
、
i
n
erfc
z
=
∑
j
=
0
∞
(
−
z
)
j
2
n
−
j
j
!
Γ
(
1
+
n
−
j
2
)
,
{\displaystyle i^{n}\!\operatorname {erfc} z=\sum _{j=0}^{\infty }{\frac {(-z)^{j}}{2^{n-j}j!\,\Gamma \left(1+{\frac {n-j}{2}}\right)}},}
i
2
m
erfc
(
−
z
)
=
−
i
2
m
erfc
z
+
∑
q
=
0
m
z
2
q
2
2
(
m
−
q
)
−
1
(
2
q
)
!
(
m
−
q
)
!
{\displaystyle i^{2m}\!\operatorname {erfc} (-z)=-i^{2m}\!\operatorname {erfc} z+\sum _{q=0}^{m}{\frac {z^{2q}}{2^{2(m-q)-1}(2q)!(m-q)!}}}
i
2
m
+
1
erfc
(
−
z
)
=
i
2
m
+
1
erfc
z
+
∑
q
=
0
m
z
2
q
+
1
2
2
(
m
−
q
)
−
1
(
2
q
+
1
)
!
(
m
−
q
)
!
.
{\displaystyle i^{2m+1}\!\operatorname {erfc} (-z)=i^{2m+1}\!\operatorname {erfc} z+\sum _{q=0}^{m}{\frac {z^{2q+1}}{2^{2(m-q)-1}(2q+1)!(m-q)!}}.}
実装
実引数の実関数として
POSIX 準拠のオペレーティングシステムでは 、ヘッダー math.h で を宣言し、数学ライブラリで libm 関数 erfおよび erfc( 倍精度 )と、それらの 単精度 および 拡張精度の 対応関数 erff、 erfl、 erfcfを提供する必要があります erfcl。 [30]
GNU 科学ライブラリは 、、、 erfスケール 誤差関数 を提供している。 [31] erfclog(erf)
複雑な引数の複雑な関数として
libcerfは、複素誤差関数の数値Cライブラリであり、 MIT Faddeevaパッケージに実装されている Faddeeva関数 に基づいて、約13〜14桁の精度で 複素関数 cerf、、 および実関数を提供します 。 cerfccerfcxerfierfcx
参考文献
^ アンドリュース、ラリー C. (1998)。エンジニアのための数学の特殊機能。SPIE プレス。p. 110。ISBN 9780819426161 。
^ ウィテカー、エドマンド・テイラー 、 ワトソン、ジョージ・ネヴィル (2021年)。 モル、ヴィクトル・ユーゴー (編)。 現代分析講座 ( 第5版改訂)。 ケンブリッジ大学出版局 。p.358。ISBN 978-1-316-51893-9 。
^ Glaisher, James Whitbread Lee (1871年7月). 「定積分クラスについて」. ロンドン、エディンバラ、ダブリン Philosophical Magazine and Journal of Science . 4. 42 (277): 294–302. doi :10.1080/14786447108640568 . 2017年 12月6日 閲覧。
^ Glaisher, James Whitbread Lee (1871年9月). 「定積分のクラスについて。パートII」。 ロンドン、エディンバラ、ダブリンの哲学雑誌と科学ジャーナル 。4. 42 (279): 421–436. doi :10.1080/14786447108640600。 2017年 12月6日 閲覧 。
^ ワイスタイン、エリック・W. 「Erf」。 マスワールド 。
^ Schöpf, HM; Supancic, PH (2014). 「Bürmannの定理と線形および非線形 熱 伝達および拡散の問題へのその応用について」. The Mathematica Journal . 16. doi : 10.3888/tmj.16-11 .
^ Weisstein, Eric W. 「Bürmann の定理」 。MathWorld 。
^ Dominici, Diego (2006). 「逆誤差関数の導関数の漸近解析」. arXiv : math/0607230 .
^ Bergsma, Wicher (2006). 「新しい相関係数、その直交分解および関連する独立性検定について」. arXiv : math/0604627 .
^ Pierre-Simon Laplace 、 Traité de mécanique céleste 、第 4 巻 (1805)、livre X、255 ページ。
^ アニー・アム・カイト ;ピーターセン、ヴィグディス B.ヴェルドンク、ブリギッテ。ワーデランド、ホーコン。ジョーンズ、ウィリアム B. (2008)。 特殊関数用の連分数ハンドブック 。スプリンガー・フェルラーク。 ISBN 978-1-4020-6948-2 。
^ シュレーミルヒ、オスカー・ザビエル (1859)。 「ウーバー ファカルテテンレイヘン」。 Zeitschrift für Mathematik und Physik (ドイツ語)。 4 : 390–415。
^ ニールソン、ニールス (1906)。 Handbuch der Theorie der Gammafunktion (ドイツ語)。ライプツィヒ:BG・トイブナー。 p. 283等量3 . 2017 年 12 月 4 日 に取得 。
^ Chiani, M.; Dardari, D.; Simon, MK (2003). 「フェージング チャネルにおけるエラー確率の計算のための新しい指数境界と近似値」 (PDF) . IEEE Transactions on Wireless Communications . 2 (4): 840–845. CiteSeerX 10.1.1.190.6761 . doi :10.1109/TWC.2003.814350.
^ Tanash, IM; Riihonen, T. (2020). 「指数関数の合計によるガウスQ関数のグローバルミニマックス近似と境界」. IEEE Transactions on Communications . 68 (10): 6514–6524. arXiv : 2007.06939 . doi :10.1109/TCOMM.2020.3006902. S2CID 220514754.
^ Tanash, IM; Riihonen, T. (2020). 「指数関数の合計によるガウスQ関数のグローバルミニマックス近似と境界の係数[データセット]」。 Zenodo . doi :10.5281/zenodo.4112978。
^ Karagiannidis, GK; Lioumpas, AS (2007). 「ガウスQ関数の改良近似値」 (PDF) . IEEE Communications Letters . 11 (8): 644–646. doi :10.1109/LCOMM.2007.070470. S2CID 4043576.
^ Tanash, IM; Riihonen, T. (2021). 「Karagiannidis–Lioumpas近似の係数の改良とガウスQ関数の境界」. IEEE Communications Letters . 25 (5): 1468–1471. arXiv : 2101.07631 . doi :10.1109/LCOMM.2021.3052257. S2CID 231639206.
^ Chang , Seok-Ho; Cosman, Pamela C. ; Milstein , Laurence B. (2011 年 11 月)。「ガウス誤差関数の Chernoff 型境界」。IEEE Transactions on Communications。59 ( 11): 2939–2944。doi :10.1109/TCOMM.2011.072011.100049。S2CID 13636638 。
^ Winitzki, Sergei (2003). 「超越関数の均一近似」 . 計算科学とその応用 – ICCSA 2003 . コンピュータサイエンスの講義ノート。第 2667 巻。Springer、ベルリン。pp. 780–789。doi :10.1007/3-540-44839- X_82。ISBN 978-3-540-40155-1 。
^ Zeng, Caibin; Chen, Yang Cuan (2015). 「一般化ミッタク・レフラー関数とその逆関数の全体的パデ近似」. 分数計算と応用解析 . 18 (6): 1492–1506. arXiv : 1310.5592 . doi :10.1515/fca-2015-0086. S2CID 118148950. 実際、Winitzki [32] は、いわゆる全体的パデ近似を提供した。
^ Winitzki, Sergei (2008 年 2 月 6 日)。「誤差関数とその逆関数の便利な近似値」。
^ Press, William H. (1992). Numerical Recipes in Fortran 77: The Art of Scientific Computing . Cambridge University Press. p. 214. ISBN 0-521-43064-X 。
^ Dia, Yaya D. (2023). 「近似不完全積分、相補誤差関数への応用」. SSRN電子ジャーナル . doi :10.2139/ssrn.4487559. ISSN 1556-5068.
^ abc Cody, WJ (1993 年 3 月)、「アルゴリズム 715: SPECFUN - 特殊関数ルーチンとテスト ドライバのポータブル FORTRAN パッケージ」 (PDF) 、 ACM Trans. Math. Softw. 、 19 (1): 22–32、 CiteSeerX 10.1.1.643.4394 、 doi :10.1145/151271.151273、 S2CID 5621105
^ Zaghloul, MR (2007 年 3 月 1 日)、「Voigt 線プロファイルの計算について: 減衰正弦関数による単一の固有積分」、 Monthly Notices of the Royal Astronomical Society 、 375 (3): 1043–1048、 Bibcode :2007MNRAS.375.1043Z、 doi : 10.1111/j.1365-2966.2006.11377.x
^ John W. Craig、「2次元信号コンステレーションのエラー確率を計算するための新しい、シンプルで正確な結果」 Wayback Machine で2012年4月3日にアーカイブ、1991 IEEE軍事通信会議の議事録、第2巻、571~575ページ。
^ Behnad, Aydin (2020). 「CraigのQ関数公式の新しい拡張とデュアルブランチEGCパフォーマンス分析へのその応用」. IEEE Transactions on Communications . 68 (7): 4117–4125. doi :10.1109/TCOMM.2020.2986209. S2CID 216500014.
^ Carslaw, HS ; Jaeger, JC (1959). 固体の熱伝導 (第2版). Oxford University Press. p. 484. ISBN 978-0-19-853368-9 。
^ 「math.h - 数学的宣言」. opengroup.org . 2018年. 2023年 4月21日 閲覧 。
^ 「特殊関数 – GSL 2.7 ドキュメント」。
さらに読む
外部リンク