|
Size: 160
Comment:
|
Size: 4988
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| = Student-t noise = | = Student-t noise modeling = ---- |
| Line 3: | Line 4: |
| 重力波望遠鏡における雑音は扱いやすさからガウス性を仮定する事があるが、<<BR>> 現実的にはガウス雑音に比べ、裾を持った非ガウスな分布となっている事が多い。<<BR>> この非ガウスさがどの程度であるかを把握する事は望遠鏡の状態を知る上で非常に重要な要素であり、<<BR>> 重力波探索に置いても非ガウスさ仮定する事でパフォーマンスが向上するという報告もなされている<<BR>> 本研究の最も重要な目的は干渉計雑音の非ガウスさを定量化し、干渉計診断の１つの指標とすることである<<BR>> この際、非ガウスさの指標にはStudent-t分布のパラメータ$\nu$を用い、<<BR>> オンライン／オフライン問わずに非ガウスさを追えるモニタを実装し、<<BR>> 望遠鏡の状態把握、重力波探索などに利用できる指標を提供する。<<BR>> ---- |
|
| Line 4: | Line 16: |
| '''1. 非ガウスさ$\nu$推定ルーチンの実装''' * '''参考文献の再現''' (〜6/10)<<BR>> * ''quantileを基準に最適な$\nu$を推定''<<BR>> * ''$\nu$の探索範囲は[4:100] ($\nu$が大きいと解析パフォーマンス改善の恩恵は無い)''<<BR>> * '''シミュレーション信号を用いたフィッティングテスト''' (〜6/17)<<BR>> * $\nu$が既知の信号に対してモデリングを行う -> 高いquantileの方が$\nu$の再現がよい(定常乱数の場合) * '''参考文献の改善 or 新たな方法の考案''' (〜7/7)<<BR>> * $\nu$の決定精度はどの程度必要か<<BR>> * dT, dFをどのように決めるか<<BR>> * どのquantileを基準とするか<<BR>> * quantile以外の基準の検討 -> 分布全体でフィッティング(最小二乗法)<<BR>><<BR>> ''' 2. リアルタイムノイズモデリング''' * '''実データを用いた準リアルタイムモデリングのテスト''' (〜7/14)<<BR>> * /frames以下のデータでテスト -> 16kHzサンプリング128秒のデータを10秒程度で処理できている<<BR>> * CLIO, TAMAデータを用いる<<BR>> * シェルスクリプト等でT秒おきにファイルを食わせて走らせる<<BR>> * '''計算コストの見積もり、削減によりリアルタイム化''' (〜7/21)<<BR>> * 計算の大半がSFTとソートorヒストグラム化(Vectorで高速化？)<<BR>> * プログラム全体としてはIO, 画面描画が遅い？<<BR>><<BR>> ''' 3. DetCharシステムに組み込み''' * '''HasKAL GUIから動かす''' (〜7/22)<<BR>> ---- |
|
| Line 5: | Line 41: |
| 1. 時系列雑音をTF平面に変換：SFT 2. 各周波数ビンごとにヒストグラム化：Cern ROOT 3. ヒストグラムのフィッティング：最小二乗法(GSL), CDFのquantile, 最尤法 ---- |
|
| Line 6: | Line 47: |
| * シミュレーション信号で与えた$nu$付近が求まる事を確認 [[attachment:sim.png]] * CLIOデータでファイルIOからプロットまでの動作を確認 [[attachment:C_1034554496.png]] * LIGOデータでいくつかのチャンネルデータをプロット ---- == Ploblem == 1. 計算高速化の目処が立たない * 16kHzの主干渉信号に対して、リアルタイム処理が実現できない(Overlapにもよる) * 低サンプリングのチャンネルだと間に合う(〜1kHz) -> Unboxな型を用いる事で解決。4kHzサンプリング、4096秒のデータが40秒弱で処理できる 2. メモリの圧迫 * 1000秒データの解析に十数GBのメモリが必要(遅延評価が原因) * 正格評価に出来れば数MBで済むはず -> Unbox化後、GBを超えない所までは落ちた。 ---- |
|
| Line 7: | Line 64: |
| (2015/3/22) 山本尚弘 他, 「重力波検出器雑音の非ガウス・非定常性雑音成分のキャラクタリゼーション」, 日本物理学会 第70回年次大会, 早稲田大学, 2015年3月<<BR>> (2015/1/10) T.Yamamoto,「Noise characterization of non-gaussianity and non-stationarity for the gravitational wave detector」A04/05 Joint Camp, Niigata, Jan. 2015<<BR>> (2014/12/20) T.Yamamoto, et al.,「Development of Non-Gaussian Noise Characterization System」, 7th Japan-Korea Workshop on KAGRA, Toyama Univ., Dec. 2014<<BR>> (2014/9/20) 山本尚弘 他,「KAGRA Detector Characterization: リアルタイム非ガウス性雑音モデリングシステムの開発」, 日本物理学会 2014年秋期大会, 佐賀大学, 2014年9月<<BR>> ---- |
|
| Line 8: | Line 71: |
| * Reference * [1] C. Rover, LIGO-T1100497 (2011) * [2] C. Rover, ''Phys. Rev.'' D '''84''', (2011) 122004 * [3] C. Rover, ''et. al, Class. Quantum Grav.'' '''28''', (2011) 015010 ---- |
Student-t noise modeling
Project description
重力波望遠鏡における雑音は扱いやすさからガウス性を仮定する事があるが、
現実的にはガウス雑音に比べ、裾を持った非ガウスな分布となっている事が多い。
この非ガウスさがどの程度であるかを把握する事は望遠鏡の状態を知る上で非常に重要な要素であり、
重力波探索に置いても非ガウスさ仮定する事でパフォーマンスが向上するという報告もなされている
本研究の最も重要な目的は干渉計雑音の非ガウスさを定量化し、干渉計診断の１つの指標とすることである
この際、非ガウスさの指標にはStudent-t分布のパラメータ$\nu$を用い、
オンライン／オフライン問わずに非ガウスさを追えるモニタを実装し、
望遠鏡の状態把握、重力波探索などに利用できる指標を提供する。
Project schedule and Task List
1. 非ガウスさ$\nu$推定ルーチンの実装
参考文献の再現 (〜6/10)
quantileを基準に最適な$\nu$を推定
$\nu$の探索範囲は[4:100] ($\nu$が大きいと解析パフォーマンス改善の恩恵は無い)
シミュレーション信号を用いたフィッティングテスト (〜6/17)
$\nu$が既知の信号に対してモデリングを行う -> 高いquantileの方が$\nu$の再現がよい(定常乱数の場合)
参考文献の改善 or 新たな方法の考案 (〜7/7)
$\nu$の決定精度はどの程度必要か
dT, dFをどのように決めるか
どのquantileを基準とするか
quantile以外の基準の検討 -> 分布全体でフィッティング(最小二乗法)
2. リアルタイムノイズモデリング
実データを用いた準リアルタイムモデリングのテスト (〜7/14)
/frames以下のデータでテスト -> 16kHzサンプリング128秒のデータを10秒程度で処理できている
CLIO, TAMAデータを用いる
シェルスクリプト等でT秒おきにファイルを食わせて走らせる
計算コストの見積もり、削減によりリアルタイム化 (〜7/21)
計算の大半がSFTとソートorヒストグラム化(Vectorで高速化？)
プログラム全体としてはIO, 画面描画が遅い？
3. DetCharシステムに組み込み
HasKAL GUIから動かす (〜7/22)
Analysis method
- 時系列雑音をTF平面に変換：SFT
- 各周波数ビンごとにヒストグラム化：Cern ROOT
- ヒストグラムのフィッティング：最小二乗法(GSL), CDFのquantile, 最尤法
Results
シミュレーション信号で与えた$nu$付近が求まる事を確認 sim.png
CLIOデータでファイルIOからプロットまでの動作を確認 C_1034554496.png
- LIGOデータでいくつかのチャンネルデータをプロット
Ploblem
- 計算高速化の目処が立たない
- 16kHzの主干渉信号に対して、リアルタイム処理が実現できない(Overlapにもよる)
- 低サンプリングのチャンネルだと間に合う(〜1kHz)
-> Unboxな型を用いる事で解決。4kHzサンプリング、4096秒のデータが40秒弱で処理できる
- メモリの圧迫
- 1000秒データの解析に十数GBのメモリが必要(遅延評価が原因)
- 正格評価に出来れば数MBで済むはず
-> Unbox化後、GBを超えない所までは落ちた。
Presentation
(2015/3/22) 山本尚弘 他, 「重力波検出器雑音の非ガウス・非定常性雑音成分のキャラクタリゼーション」, 日本物理学会 第70回年次大会, 早稲田大学, 2015年3月
(2015/1/10) T.Yamamoto,「Noise characterization of non-gaussianity and non-stationarity for the gravitational wave detector」A04/05 Joint Camp, Niigata, Jan. 2015
(2014/12/20) T.Yamamoto, et al.,「Development of Non-Gaussian Noise Characterization System」, 7th Japan-Korea Workshop on KAGRA, Toyama Univ., Dec. 2014
(2014/9/20) 山本尚弘 他,「KAGRA Detector Characterization: リアルタイム非ガウス性雑音モデリングシステムの開発」, 日本物理学会 2014年秋期大会, 佐賀大学, 2014年9月
Document
- Reference
- [1] C. Rover, LIGO-T1100497 (2011)
[2] C. Rover, Phys. Rev. D 84, (2011) 122004
[3] C. Rover, et. al, Class. Quantum Grav. 28, (2011) 015010