Size: 2301
Comment:
|
Size: 2603
Comment:
|
Deletions are marked like this. | Additions are marked like this. |
Line 15: | Line 15: |
* $\nu$が既知の信号に対してモデリングを行う | * $\nu$が既知の信号に対してモデリングを行う -> 高いquantileの方が$\nu$の再現がよい(定常乱数の場合) |
Line 20: | Line 20: |
* quantile以外の基準の検討<<BR>><<BR>> | * quantile以外の基準の検討 -> 分布全体でフィッティング(最小二乗法)<<BR>><<BR>> |
Line 24: | Line 24: |
* /frames以下のデータでテスト -> 16kHzサンプリング128秒のデータを20秒程度で処理できている<<BR>> | |
Line 35: | Line 36: |
2. 各周波数ビンに対してStudent-tでフィッティング:(求根アルゴリズム) | 2. 各周波数ビンごとにヒストグラム化:Cern ROOT 3. ヒストグラムのフィッティング:最少二乗法(GSL) |
Student-t noise modeling
Project description
検出器雑音の非ガウスさを定量化し、干渉計診断の1つの指標とする
非ガウスさ指標はStudent-t分布のフィッティングパラメータ$\nu$を用いる
オンラインDetChar(GUIシステムからのオフラインDetChar)として動かす
Project schedule and Task List
1. 非ガウスさ$\nu$推定ルーチンの実装
参考文献の再現 (〜6/10)
quantileを基準に最適な$\nu$を推定
$\nu$の探索範囲は[4:64] ($\nu$が大きいと解析パフォーマンス改善の恩恵は無い)
シミュレーション信号を用いたフィッティングテスト (〜6/17)
$\nu$が既知の信号に対してモデリングを行う -> 高いquantileの方が$\nu$の再現がよい(定常乱数の場合)
参考文献の改善 or 新たな方法の考案 (〜7/7)
$\nu$の決定精度はどの程度必要か
dT, dFをどのように決めるか
どのquantileを基準とするか
quantile以外の基準の検討 -> 分布全体でフィッティング(最小二乗法)
2. リアルタイムノイズモデリング
実データを用いた準リアルタイムモデリングのテスト (〜7/14)
/frames以下のデータでテスト -> 16kHzサンプリング128秒のデータを20秒程度で処理できている
CLIO, TAMAデータを用いる
シェルスクリプト等でT秒おきにファイルを食わせて走らせる
計算コストの見積もり、削減によりリアルタイム化 (〜7/21)
計算自体はほとんどがSFTとSort (要Vector化)
プログラム全体としてはIO, ROOT起動が遅い?
3. DetCharシステムに組み込み
HasKAL GUIから動かす (〜7/22)
Analysis method
- 時系列雑音をTF平面に変換:SFT
- 各周波数ビンごとにヒストグラム化:Cern ROOT
- ヒストグラムのフィッティング:最少二乗法(GSL)
Results
Presentation
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