#author("2021-04-09T20:37:38+09:00","default:LabMember","LabMember")
#author("2021-04-10T22:46:37+09:00","default:LabMember","LabMember")
[[PythonTutorial]]

* Pythonを使ったデータ解析その2 [#xc25415b]

#contents

** 準備 [#d44aa4e0]
[[ex2.zip>https://granite.phys.s.u-tokyo.ac.jp/wiki/Lab/?plugin=attach&pcmd=open&file=ex1.zip&refer=PythonTutorialExample1]]
[[ex2.zip>https://granite.phys.s.u-tokyo.ac.jp/wiki/Lab/?plugin=attach&pcmd=open&file=ex2.zip&refer=PythonTutorialExample2]]
をダウンロードして展開し,ex2.ipynbを開く

*** 中身について [#y6c67512]
今回はミッドフリンジに制御したマイケルソン干渉計の制御信号から,元の外乱のスペクトルを推定する.
- ex2.ipynb: メインのデータ解析をするnotebook
- OLTF.csv: 制御のオープンループ伝達関数
- feedback.csv: フィードバック信号のパワースペクトル密度

** アクチュエーター伝達関数の推定 [#se11de68]
外乱を推定するために,測定したオープンループ伝達関数からアクチュエーターの伝達関数を推定する.
そのためにはセンシング効率,フィルターの伝達関数の情報が必要である.

*** フィルターの伝達関数 [#sd51c88c]
ここでは[[前回>PythonTutorialExample1]]の解析で用いたものと同じフィルターを使用している.
つまり,ゲインが100倍,カットオフが1Hzのローパスフィルターである.

*** センシング効率 [#y854e4d2]
マイケルソン干渉計をミッドフリンジに制御した時のセンシング効率は,フリンジの最大値Vmax,最小値Vmin,波長λで求まる(細かい式はnotebookなり誰かの修論・D論なりを参照).
今回フリンジの最大値/最小値は
 Vmax = 0.232 [V]
 Vmin = 0.020 [V]
と与えている(事前の測定で求めた).

***アクチュエーターの伝達関数 [#sa47bdc8]
アクチュエーターの伝達関数は,今回はピエゾで鏡を駆動しているので測定帯域ではほぼ周波数に依存しないとしてよい.
なので,値 A をフィッティングパラメーターとしてオープンループ伝達関数のゲインからフィッティングにより推定する.

** 外乱の推定 [#nb02d06b]
推定したアクチュエーターの伝達関数を用いて外乱のパワースペクトルの推定をする.
細かい式はnotebookなり誰かの修論・D論なりを参照.

トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS