ざるごのブログ

ブラックホール撮影にも使える「スパースモデリング」とは?【機械学習】

投稿日時: 2023/03/08 17:11

最終更新日時: 2024/06/30 04:58

どうもざるご(@zalgo3)です.

世界初のブラックホール撮影の成功例が出たようです.

ブラックホールの撮影に成功 世界初 一般相対性理論を証明 - 毎日新聞

今回のブラックホール撮影は,スパースモデリングという機械学習技術を取り入れたことによる貢献が大きいようです.

天文学に計算機科学の知識が取り入れられて,大きな成果が出たというのは,驚くべきことだと思います.

今回はそんな大成功を巻き起こした「スパースモデリング」について解説していきます.

スパースモデリングとは

スパースモデリングとは,誤解を恐れずにざっくりいうと解けない連立一次方程式を無理やり解くための仕組みです.

次のような連立一次方程式を考えます.

Ax=yAx = y

この方程式は,AAが逆行列を持つときに解くことができて,

x=A1yx = A^{-1} y

となります.では,AAが逆行列を持たないときはどうなるでしょうか?

高校や大学1年生で習った数学では,そういう連立方程式は解けないと習ったと思います.「不能」とか「不定」とかいうやつですね.

ただ,実用上の観点からいうと,せっかく式があるのにハイ解けませんで片付けてしまうのは少々もったいないところがあるんですよね.

そこで登場するのがスパースモデリングです.スパースモデリングでは,以下の2つを仮定することによって解なしの連立一次方程式でも無理やり解を求めてしまいます.

  1. Ax=yAx = yが完全に満たされなくてもよく,AxyAx \approx yと近似的に成り立てばよいものとする.
  2. xxが**スパース(疎)**であるとする.

スパース(疎)というのは人によってはあまり馴染みのない表現かもしれません.

スパースというのは,という字からも連想されるように,データ(0以外の要素)がまばらであることをいいます.

ベクトルxxがスパースであるというのは,xxの要素のほとんどが0で,0以外の要素がちょっとだけある,つまり

x=(0,,0,a,0,,0,b,0,)x = (0, \dots, 0, a, 0 , \dots, 0, b, 0, \dots )^\top

のようなベクトルであるということです.

実世界で観測されるデータには,スパースなものが多いということが知られています.

今回の研究で扱われた宇宙の観測データも,スパースなデータの1つです.

広大な宇宙は、ほとんど真っ暗な世界です。この中の限られた部分だけが天体として観測されるわけで、そうした意味でスパース性を仮定することは妥当だといえるでしょう.

では具体的にこれをどう解くか,というところを順を追って解説していきます.まず1から.

なるべくAxAxyyに近づくようなxxを手に入れたいわけです.そこで,次の関数を最小化することを考えます.

Axy22||Ax - y||^2_2

x2||x||_2という記号は見慣れない人もいるかも知れませんが,高校で習うベクトルの絶対値x|\vec{x}|と同じです.つまり,上の問題の意味は,AxyAx - yの各成分の二乗和(二乗誤差)を最小化するという意味になります.

二乗和は必ず0以上となるので,この関数を最小化すれば0に近い値が得られるはずです.したがって,AxyAx \approx yとなるようなxxが得られるという理屈になります.

この方法を**最小二乗法 (Least Squares Method)**といいます.

しかし,これだけではxxがスパースになってくれるとは限りません.そこでさっきの関数の代わりに次の関数を考えます.

Axy22+cx0||Ax - y||^2_2 + c ||x||_0

c>0c > 0は適当な正の定数を表します.

x0||x||_0

というのは,xxの非ゼロ要素の個数を表す記号です(L0ノルムという).たとえばx=(0,2,0)x = (0, 2, 0)^\topのとき,x0=1||x||_0 = 1となります.

最小二乗法と違って,二乗誤差と一緒に非ゼロ要素の個数を最小化することになるので,非ゼロ要素の数が少ないスパースな解が得られます.

しかし,「非ゼロ要素の数」というのは連続的に変化せず,数学的に結構扱いづらい性質を持っています.そこでこれと似たような問題として,

Axy22+cx1||Ax - y||^2_2 + c ||x||_1

の最小化を考えることがあります.

x1||x||_1は,xxの各成分の絶対値の総和を表します(L1ノルムという).

この問題も,解がスパースになるということがよく知られています.

さっきのL0ノルムと違って,L1ノルムは連続で扱いやすいため,こちらの問題が実用上はよく扱われます.

この問題を,**LASSO (Least Absolute Shrinkage and Selection Operator)**といいます.

さて,ここまでざっくりとしたアイデアを説明してきましたが,次は**このアイデアをブラックホール観測にどう活かしたのか?**ということについて説明していきます.

ブラックホールとスパースモデリング

アインシュタインの相対性理論や,重力波の観測などの実験結果から,ブラックホールが存在することはほぼ確実とされています.

そこで近年では,ブラックホールが作り出す影(ブラックホールシャドウ)を撮影することで,ブラックホールの存在を立証しようという計画が進められてきました.

ブラックホール観測を成功させるための秘訣は,望遠鏡の空間分解能をあげることにあります.

「空間分解能」というのは,近接した2点をきちんと2点として認識できる能力,という意味です.

たとえば道路に家が2軒並んでいるような状況を考えましょう.

近所から見ると,「家が2軒あるな」ということはすぐにわかります.しかし,これを遠く離れた宇宙から見るとどうでしょうか?

おそらく人間の目では家が2軒あることを認識するのは不可能でしょう.では,望遠鏡で見るとどうでしょうか?

宇宙から撮影された衛星写真を見ると,家が2軒あることくらいは余裕で識別することができます.

この二つの違いを生じさせるのが空間分解能です.望遠鏡は,人間の目に比べて非常に高い空間分解能を持っています.

しかしそんな天体望遠鏡の空間分解能でも,ブラックホールを観測するには全く歯が立ちません.

ブラックホールとブラックホールの距離は,約1光年と言われています.一方,ブラックホールが存在する可能性のある銀河系で,地球から最も近いものでも,239万光年離れていると言われています.

スケールがでかすぎてよくわかりませんが,239km先にある1cm間隔のパチンコ玉を識別できるかと聞かれれば,その無理ゲーさが理解できると思います.

しかし,この無理ゲーを解決してこそ科学者です.そこで使われた解決策が,電波干渉計,そしてスパースモデリングです.

まず電波干渉計について.これは2つ以上の電波望遠鏡を連動させて,波の干渉を用いて空間分解能を向上させる技術です.

「波の干渉」といったら,高校物理を履修していた人は一度は聞いたことがあると思います.ヤングの干渉実験とかが有名なやつですね.

光は波の性質を持っているので,干渉が起こります.電波干渉計のアイデアは,2つの望遠鏡で観測された波が干渉したときの時間差を測定することで,波がどちらの方向から来たかを割り出すというところにあります.

この技術によって,空間分解能は飛躍的に上がります.

しかし,ブラックホールの観測にはこれだけじゃまだ足りないのです.ブラックホール恐るべし.

そこで近年になって注目されてきたのが,スパースモデリングです.

ブラックホール観測にスパースモデリングをどう用いるか,ぼくも完全に理解しているわけではないですが,ざっくりと書きます.

まず,電波干渉計で観測されたデータVVというのは,実際の天体画像IIに,観測行列(フーリエ変換)FFをかけてやったものと考えられます.つまり,

V=FIV = FI

です.

目標は,天体画像IIを,なるべく高解像度で復元してやることとなります.(空間分解能をあげるため)

基本的に望遠鏡の数はそんなにたくさんは使えないので,FFの逆行列を求めて方程式をそのまま解くという方法は使えません.

そこで生きるのがさっきのスパースモデリングです.

次の関数を最小化して(LASSO),天体画像を復元してやります.

VFI22+cI1||V - FI||^2_2 + c ||I||_1

この方法により,実際に研究チームが復元した画像がこちらです1

この研究からさらに様々な改良が重ねられ,今回の報道が合ったような大発見につながったようです.

まとめ

5000文字を超える長文記事になってしまいましたが,今回の大発見のすばらしさ,そしてスパースモデリングの魅力について少しでも伝えられていたら幸いです.

まだぼくもスパースモデリングについては勉強中の身なので,間違っている点等あったら申し訳ないです.ぜひコメントとかで教えてくださると嬉しいです.

また,この記事でスパースモデリングに興味を持った人向けに,入門の教科書のリンクをいくつか貼っておきます.一緒に勉強しよう!

スパース推定法による統計モデリング (統計学One Point)

スパースモデリングで用いられるモデルがざっくりと書いてある本. 細かい証明等はほとんどないが,イメージを掴むのにおすすめ.

スパース推定法による統計モデリング (統計学One Point) ")

スパース推定法による統計モデリング (統計学One Point)

スパースモデリング- 基礎から動的システムへの応用

基礎から応用までしっかり書かれている本. ソースコード例も載っているので,実際にプログラムを実装しながら学びたい人にはおすすめ.

スパースモデリング- 基礎から動的システムへの応用 -

スパースモデリング- 基礎から動的システムへの応用 -

スパース性に基づく機械学習 (機械学習プロフェッショナルシリーズ)

数学的な背景までスパースモデリングをガッツリ知りたい人向け. ぼくのように証明がないと発狂してしまう人類はこちらを.

スパース性に基づく機械学習 (機械学習プロフェッショナルシリーズ) ")

スパース性に基づく機械学習 (機械学習プロフェッショナルシリーズ)

宣伝

学術系の知識とかをおもしろおかしく解説するYouTubeチャンネルをやってます. ↓の動画は「安定マッチング」というグラフ理論の問題について解説した動画です. まだ動画数は少ないですが,今後機械学習やスパースモデリングについても動画を出していこうと思っているのでよければチャンネル登録よろしくお願いします!

【京大式】数学的に最も正しい彼氏彼女の作り方を解説【安定結婚問題】

*1:本間ら「スパースモデリング天文学 — ブラックホール撮像から時間変動減少まで」,科学研究費補助金新学術領域研究「スパースモデリングの深化と高次元データ駆動科学の創成」最終成果報告会 (2017/12/18-20)

Footnotes

  1. 本間ら「スパースモデリング天文学 — ブラックホール撮像から時間変動減少まで」,科学研究費補助金新学術領域研究「スパースモデリングの深化と高次元データ駆動科学の創成」最終成果報告会 (2017/12/18-20)