2

誤差関数の公式

461
0
全般的に正当化不能な計算を行っているので注意

以前紹介したHankel変換など偶関数に対してしか定義されなかったものの定義域を拡張する試みの一環で、厳密化しようとする過程にある計算なので怪しいものが多数あるが、きれいなので紹介する

誤差関数erf(x)=2π0xet2dt
相補誤差関数erfc(x)=1erf(x)
スケーリング相補誤差関数(scaled complementary error function)
erfcx(x)=ex2erfc(x)
Dawson関数D(x)=π2iex2erf(ix)

1f(x)=0xf(t)dt
のようにの(右)逆元を定義することができる。これを使うと
D(x)=ex20xet2dt=ex21ex2=ex21ex21=(ex2ex2)11=(+2x)11=(e2/4xe2/4)11=e2/4(2x)1e2/41=e14212x
となる。実際、途中にある式が示唆するように、Dは微分方程式
dD(x)dx+2xD(x)=1
を満たす。両辺を微分すると簡単に示せる次の式
et2sin(xt)tdt=π erf(x/2)
を使うと Weierstrass変換 のCauchy主値として
D(x)=12πe(xt)2dtt=e14212x
となっていることがわかる。書き換えると次のようになる。
erf(x)=ex212x2π=1πex2e1421x
次の公式について考える。

12πσ2erf(ax+b)exp((xμ)22σ2)dx=erf(aμ+b1+2a2σ2)

Weierstrass変換を用いると次のように書ける(ただしA=1+2a2σ2)
eσ222erf(ax)=erf(ax/A)
eσ222ea2x2e14a221x=Aea2Ax2eA4a221x
LCT を用いると
T+(10σ2i1)T+(12ia201)T+(10i/(2a2)1)1x=AT+(12ia2/A01)T+(10iA/(2a2)1)1x
という等式に言い換えられる。
(A00A1)=((10σ2i1)(12ia201)(10i/(2a2)1))1(12ia2/A01)(10iA/(2a2)1)
というSL(2,C)の等式を使うことで準同型性から

1x=AT+(A001/A)1x
という式に帰着する。この議論を逆行することで公式をしめすことができるように思えるが、LCTのパラメータがSL(2,R)ではなくSL(2,C)なので、普遍被覆群の基本群に由来する係数の部分を定めることができない。
もとに戻ると、
erfcx(x)=ex2+1πe1421x
となる。Weierstrass変換は積分変換として定義されていたが、指数関数の級数展開を行うとn!x2nのオーダーで発散してしまうため、微分演算子としては通常の計算ができない。しかし、漸近展開としては計算できる。
2nx1=(2n)!x12n=(1/2)n(1)n4nx12n
なので無理やり計算すると
erfcx(x)=ex2+1xπ 2F0(1,12;;1x2)

これは漸近展開の公式
erfcx(x)1xπ 2F0(nterm)(1,12;;1x2)
を示唆している。2F0(nterm)n項までの部分和である。
他の関数としてQ-functionがある。
Q(x)=12erfc(x/2)=12+12πe12x2e1221x
となる。

投稿日:20231120
更新日:20231123
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

赤げふ
赤げふ
92
15293
東工大情報B4 数学,理論物理,Minecraft計算機/微分演算子の記事を書きます/主に表現論,量子群,物理の数理に興味があります

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中