のぺーっと分布した超関数
ディラックのデルタ関数
はい。まあもちろん知らない人が大半です。なのでちょこっと説明します。
デルタ関数の定義
定義はこちら。
これを満たす をデルタ関数といいます。
ちなみに関数ではありません(!)。
関数に似てるけど、関数としては性質がふざけてるので、「超関数」という分類になってます。
これに を代入すると
確率密度関数列の極限としてのデルタ関数
実はデルタ関数は他の確率密度関数から類推してやると理解しやすいです。
2つ紹介します。
正規分布の分散→0にしたもの
分散が となってます。
こいつの を大きくしていくと…
あ、1つめの関数がどの位置にあるかとか、nが大きくなる時の変化の幅とかは適当です。
要するにどういう動きをするのかってことだけわかれば良いので。以下同様。
一様分布の幅を→0にしたもの
一様分布の範囲の幅が になってます。
この を大きくしていくと…
という感じ。
どちらの場合も、面積=1 をキープしたまま、y軸にピッタリと張り付いていく様子がわかりますね。
確率密度関数としてのデルタ関数
例えばですね、これを「x軸上で、ある粒子がどの位置にあるか、という確率を表す」と解釈することが出来るわけですよ。
そしてデルタ関数の場合、「100%絶対に の地点にある」と断言できるわけです。
普通、x軸という連続な範囲において、ある厳密な一点に粒子がある確率は0%です。ほんのちょっとでもズレてたらアウトなので、そういうことになるんです。
「 の範囲にある確率」みたいに範囲で考えて初めて、意味のある確率が出てくるわけです。
なのに、デルタ関数の場合は、厳密に の地点にあると断言できるので、気持ち悪いわけですね。
このあたりからも、デルタ関数が普通の関数ではなくて超関数であることがわかります。
のぺーっとした超関数を考えよう!
ここまでデルタ関数を紹介してきました。
これに対して、のぺーっとした超関数はないのか、というのが今回の疑問。やっと本題ですよ!!
やっとって言っても、texで式書いたりgif動画作るのに時間がかかったのは僕の話であって、読んでると一瞬かもしれないですけどね!!
はい。
「のぺーっとした関数」でググっても、ツイッターで聞いてみても、全然それっぽい話がヒットしなかったので、調べるのは諦めて自分で考えることにしたのが事の発端です。
確率密度関数列の極限としてのイプシロン関数
まずはデルタ関数の時と同様に、確率密度関数の極限を使ってイメージを掴もう!
一様分布の幅を→∞にしたもの
一様分布の範囲の幅が になってます。
この を大きくしていくと…
どちらも面積=1をキープしたまま、x軸にのぺーっと張り付いていく様子が見て取れます。
関数列の収束について
このイプシロン関数を考えている時に不安だったのが、
これ、完全に に一致してしまって、イプシロン関数なんてものを考えることはできないんじゃないか?
ということです。
0.999999.....が1と等しいみたいな。1と別個に0.99999.....という数字を考えることはできないわけです。
実際、この とか は、 に一様収束します。
各点収束もするし、一様収束もするんです。
一様収束は、各点収束よりも割りとガッツリ収束していくので、極限とったら完全に同じものになるんじゃないか。
…と思ってたんですけど。
これ、nをどれだけ大きくしても積分はずっと1なんですよ。一方、 は積分したら0です。
これは、はっきり違いますね。全然一致してないです。
先に極限とって積分するのと、積分してから極限とるのとで、値が違うわけですね。
積分と極限の順序を入れ替えてもいいのは「一様収束」「積分範囲が有限」という2つの条件がありましたが、
今回は、一様収束はしても、積分範囲が無限だから、積分と極限は順序を入れ替えたら値が変わりうるわけですね。
なるほど。
ということで。
どういう関数(超関数)のことを考えているか、大体伝わったと思います。
こいつの性質を調べよう!!!
確率密度関数としてのイプシロン関数
デルタ関数は
「x軸上で、ある粒子が100%絶対に に存在する。」
という確率を表していました。
それに対して。
イプシロン関数は
「x軸上で、ある粒子が存在する確率が、x軸上のありとあらゆる点で同様に確からしい」
と解釈するのが自然でしょう!!!!
…いや、「!」マークまで付けて強調しましたが、それが自然だと感じているのが僕だけだったらごめんなさい。
累積分布関数
イプシロン関数を確率密度関数と考えるなら、その累積分布関数を考えるのは自然な流れです。これは自然な流れです。
で、累積分布関数を考えてみるんですが…
これが、よくわからないんですよ!!
累積分布関数である以上、マイナス∞の極限で0、プラス∞の極限で1になります。これはイプシロン関数の性質「積分すると1」と合致します。
そして、その (-∞,0) と (∞,1) の2点を「水平な直線」で結んだものが、累積分布関数になると思われます。
…うん。いや、飛躍してますかね。どうですかね。
イプシロン関数の値自体はずっと(ほぼ?)0なので、その積分である累積分布関数は傾き(ほぼ?)0の線なわけです。つまり(ほぼ?)直線です。
「(ほぼ?)」といってるのは、厳密に0なのかどうなのか僕がよくわかってないからです。
いやだって厳密に0だったら積分が1にならないじゃない!でも有限の値とってかつx軸全体で同様に確からしかったら、積分が発散しちゃうのでやっぱり厳密に0だよなあ…。
まあいい。
とにかく、水平な直線ということで納得しましょう。
イプシロン関数を定義しよう
とまあ、感覚的な話が続いてしまいましが
数学の記事を書いてるのに感覚的な話に終始していてはいけません。
累積分布関数がどこを通るかわからないのは、きっとイプシロン関数をちゃんと数式で定義していないからなのでしょう。
ということで、定義を試みます。
みたいに、他の関数と掛け算して積分した時にどうなればいいか、っていう観点で定義できそうかな?
どうなればいいんでしょうね。
また感覚的でアレですが、この積分は「 の平均値」を表していると考えられます。
それも普通の平均ではなく、(-∞,∞)という範囲における平均です。
(-∞,∞)という範囲において、イプシロン関数は「同様に確からしいという確率」を表すわけですから、それとを掛け算して積分したら、f(x)の平均値がでるんです。
「値×確率」を全部足すって奴です。説明が超雑ですが、ここで確率論の説明してると長くなるのでこの程度にします。
さて。
まず、 の時は、この積分は1になってなきゃいけません。
というわけです。イプシロン関数に要求している性質です。まあそれはいい。
が定数の時は、その定数がそのまま積分値になりますね。まあそれもいい。
の場合はどうでしょう?
この積分は、x軸上で粒子が発見される位置の「期待値」を表します。「値×確率」を全部足すっていう奴です。
この積分はどうなるでしょうか?
の左半分と右半分で相殺するから、積分結果は0…?
…。
はい、違いますね。
その考え方は、 を「真ん中」と捉えている場合の考え方です。
を「真ん中」としてその左右が相殺すれば、積分値は3になるわけです。
結局、この積分が「わからない」んですよね。宇宙に中心はないわけです。
ということで、
の「???」の部分、ひいては
の「???」の部分が埋められないので、これうまく定義できないですね。
いや、方法はあるのかもしれないけど、今の僕ではうまく扱えません。
そして、実際、うまく扱う方法がなくて何の役にも立たないから、この関数は検索してもヒットしなかったのかもしれません。
そんな気がします。
…
……
………
いや、まだだ!!!!!!!
デルタ関数だって、 を特別な点として扱っていたじゃないか!!
イプシロン関数でも、 を特別な点として扱ってみたらどうだ?
主観的に中心を決めてしまえば良いわけですよ。
だって、あなたがこの宇宙にある何かを感じる時、どうしたってそれはあなた一人が感じているわけで
他の人とまったく同じ感じ方をしていることはありえないわけで
相対性理論だって観測者によって時間の感じ方とか質量とかが変わるって言ってるわけで
まあ要するに
宇宙の中心はある!!
それはお前自身だ!!
いや、むしろ、この文を書いている俺からしたら、宇宙の中心はお前ですらない
宇宙の中心は、常にこの俺だ!!!
というわけでですね。
まったく飛躍のない極めて論理的な推論によって、 を特別扱いすることが根拠付けされました。
今度こそイプシロン関数を定義する
はい。いきなりでアレですが、これでしょうね。
右辺の積分範囲がポイントです。
プラス側とマイナス側に、 を「真ん中」として、同じ速度で積分範囲が大きくなるようになっています。
そして、その積分範囲の幅で積分値を割ることで、 の の範囲での平均値を出します。
最後に、その を にふっ飛ばすわけです。
f(x)にいろいろ入れてみる
f(x)にいろいろ関数をいれてみて、
それらの関数の「x=0を『真ん中』とした平均値」を求めてみましょう。
2次式
これ計算すると、無限大に発散します。
aが正なら正、負なら負の無限大に発散します。
多項式
計算しませんが、奇数次の項しかなければ0だし、2次以上の偶数次の項があれば正か負の∞に発散しますね。
結局、二次以上の項は絶対値がどこまでも大きくなるので、プラスとマイナスで相殺しない限り、その平均値は無限大になっちゃうわけですね。
指数関数
これも発散します。二次式の時と話は同じです。
累積分布関数を確認する
ヘヴィサイドの階段関数 - Wikipediaに書いてあるのと同じようにやってみます。
まず「指示関数」とかいうものを持ってきます。
これです。範囲内にあれば1、なければ0。
これを使って、累積分布関数、つまり「マイナス無限大からxまでの積分」を
と定義します。
xより大きい所では、指示関数が0になってくれるので、結局右辺と左辺は同じ値になるはずですね。
ただし、あくまでも「マイナス無限大から『ある有限の値xまでの積分』」を定義しただけです。この点、覚えておいてください。
ところでこいつは、イプシロン関数の定義より
いきなりg(s)とかいうのが出てきましたが、指示関数 の不定積分(原始関数)です。
x以上の所では定数cで、xより小さい所では傾き1の直線(s-x+定数c)です。
xの関数ではなくてsの関数の話をしているので気をつけてください。この部分ではxは定数です。
ということで、やっと答えが出ましたが、累積分布関数はなんと です!!
まあ、上でも「水平な直線」と言ってたので、わかってたっちゃわかってたんですが、xに依存しないんですね!!
どっか適当なx軸上の点を持ってきて、「その点より左にある確率」を考えると、いつだって1/2だということです。持ってきた点が宇宙の真ん中になるわけですよ多分。おもしろっ。
x→無限大にしても1にならないけど?
累積分布関数はx→無限大とした場合1にならなければいけません。
でもこの 1/2 とかいう累積分布関数は、x→無限大にしてもずっと 1/2 ですね。
ここで思い出してほしいのが、上の定義はあくまで、マイナス無限大から『ある有限の値xまでの積分』」を定義しただけだったということです。
x→無限大とか、x→マイナス無限大の時については何も言ってないわけですよ。
有限の時しか成り立たないですよ、と断って話をしていたので、「無限にしたら矛盾するけど?」とか言われても「知らん。」というわけです。
じゅう無限の時はどうするのって話なんですけど、
x→無限大の時は1、x→マイナス無限大の時は0、って、定義すればいいんじゃないかな!それで完璧に条件を満たしますよ!!
極限とるタイミングをかなり気をつけて扱わなきゃいけない関数の完成ですね。これも超関数なのかな?
粒子の位置を表す関数として解釈する
デルタ関数が「x=0に粒子がある」という確率を表していたので、イプシロン関数についても改めて解釈してみましょう。
…といっても、イプシロン関数から読み取れることは、「x軸全体のどこかにある」という情報だけです。「どこにある確率が高い」とかそういうことは何も言えません。
というのはまあ、そうなんですが
ここで、次の値を計算してみます。
この計算は、「粒子が存在する位置の期待値」を表します。
「値×確率」を全部足してるわけです。何回言ったっけこれ。
ただこれ、xって単なる一次式なので、上で計算したように、答え0なんですよね。
平均値としては、x=0になるんです。
宇宙の中心はあなたというわけですよ。つまりね。そういうこと。
ただ、決して「x=0に存在する確率が高い」とは言ってないことに注意してくださいね。
そうなってくると、「平均はx=0」の意味って何なんでしょうね。何か言えるんですかね。
xが正とxが負の部分で、どちらも確率が1/2なので、まあ、観測を何度も何度もやったら、平均はx=0に収束していくっていうことですかね。
ところがこれすごいのが、x=0の地点にいるあなたがやってもそうなるんですけど、
x=10の地点にいる高橋君が同じことをやってもそうなるんですよ。高橋くんにとってもやっぱり、自分を中心に粒子が発見されるわけです。
x=2億の地点にいるロドリゲスにとっても同じです。ロドリゲス中心に粒子が発見される。
わけわかんないですけど、そういうことなんですよ多分。
あなたにとって世界の中心はあなただけれども、高橋君にとっては世界の中心は高橋君なんです。あなたと高橋君は違う地点に立っているにも関わらず、です。
まあ、そういう前提を置いたんだから、この結果は当然なんですけどね。多分。
原点以外のところへずらす
ここから「EMANの物理学」のデルタ関数の記事を参考にしてます。
ネットで物理とか数学とか調べたことのある人はきっと知っている有名サイトです。
いや、まあ、ここまでも参考にしてたけど。
つーことで。
ε(x-a)について考えると
x'=x-a です。
まあ、そうですね、x=aを「真ん中」とした時の f(x) の平均値がでます。上の高橋君とかロドリゲスとかの話の根拠ですね。
偶関数に似た性質
がんがん省略しますが、
定義式
の右辺について
が言えるのと、左辺について
が言えるのとかを使って、引き算したりすると
が言えます。ほえ〜。
「偶関数『に似た』」と言っている意味合いについては元サイトを御覧ください。
変数のスケールを変える
について考えよう。
つってもまあ、普通に、横方向に 1/|a| 倍されるイメージで良くて
です。|a|が1より大きいなら、面積が横方向に縮んでるわけです。普通の関数のイメージと変わらないですね。
「
これは変だ。横から押し潰された結果としてなぜか高さが縮むことを表しているように思える。普通の などのグラフと比べてしまうと逆のことが起きているイメージなのだが、こう考えないと辻褄が合わないので受け入れるしかない。
」
という記述があります(閲覧日2017-03-18)が、これ、積分記号を省略してるからですね。
積分記号を省略せずに解釈すれば、あくまでも横に縮むイメージのままで大丈夫です。
微分してみる
数式でちゃんと言えないか試してみよう。
…。
はい。
これ0じゃないですね。
こんな感じになりました。元サイトと同じ議論です。
はい。
ε'(x)をf(x)にかけて積分すると、なぜか、f(x)を微分した関数の、全体で平均とった奴の、符号を反転させたものがでてきました。
多分これが、ε'(x)の本質ですね。
δ'(x)と同じで、奇関数的な性質も持っているのでしょう。
……。
を0とした所だけ、ちょっと議論が怪しいか…??
…ちょっとどころじゃないな。
f(x)が発散するとき、イプシロン関数かけても0になる保証はどこにもないじゃないか。
微分に関してはもっと考える必要があるな。
まとめ
ということで、ディラックのデルタ関数の逆バージョンみたいなものを定義することができましたね。
なんでこれ全然言及されてないんでしょうね?????
使い道がないのか…?
いやでも、量子力学で、粒子の運動量が確定してる時の位置の確率分布は、まさにこのイプシロン分布になるでしょう。
運動量と位置のどちらかがデルタ関数の時、もう一方はこのイプシロン関数になるわけですよ。なるよねえ?
あくまで俺の検索が糞なのか…??
おわり
おわりです。
「EMANの物理学」のデルタ関数の記事ではもっといろいろやってますが、いったん、またの機会にします。
f(x)やε(x)に複素数とかデルタ関数とかイプシロン関数とか突っ込んだりするとまだいろいろ調べられますね。多分ね。定義のやり直しもいるかも?
あとは、デルタ関数の応用先を調べて、同じ所でイプシロン関数も応用が効くかどうか調べると面白そう。
ε(δ(x)) と δ(ε(x)) の性質を調べるくらいは、今日帰ったらやってみようかな。
はい。
それにしても。
頭の中でバーっと考えるのにかかった時間を1とすると、
この記事を書くのに40くらいの時間がかかってます。バランス糞ですね。
アウトプットって大変だ。