677
Comment:
|
1659
|
Deletions are marked like this. | Additions are marked like this. |
Line 1: | Line 1: |
<<TableOfContents(2)>> |
|
Line 3: | Line 5: |
= Time Seriesなどの計算小技 = | 工事中 ----- = TimeSeriesなどの計算小技・注意事項 = |
Line 6: | Line 12: |
gwpyのTime SeriesやFrequency Seriesは、四則演算や指数の演算がそのまま使える。横軸の配列(timesやfrequencies)が同一なら、Time Series同士の演算も可能。 | gwpyの[[https://gwpy.github.io/docs/stable/api/gwpy.timeseries.TimeSeries.html#gwpy.timeseries.TimeSeries.psd|TimeSeries]]や [[https://gwpy.github.io/docs/stable/api/gwpy.frequencyseries.FrequencySeries.html#gwpy.frequencyseries.FrequencySeries|FrequencySeries]]は、 四則演算や指数の演算がそのまま使える。 横軸の配列(timesやfrequencies)が同一なら、TimeSeries同士の演算も可能。 |
Line 15: | Line 24: |
== FFT length の注意 == TimeSeriesからスペクトルを計算する際、fftlengthの値が小さいと下図左のように謎の冪スペクトルに埋もれてしまうことがある。(ドリフト+窓関数の影響?)<<BR>> 十分なfftlengthをとればいいのだが、平均回数を稼ぎたいときや短時間のみを見たい場合にはそうもいかない。<<BR>> そんなときは先に適当な[[https://gwpy.github.io/docs/stable/api/gwpy.timeseries.TimeSeries.html#gwpy.timeseries.TimeSeries.highpass|highpass]]フィルターをかけることで、下図右のように回避できる。<<BR>> {{attachment:wo_hp.png}} {{attachment:w_hp.png}} |
numpy, scipyなどのメモ
工事中
TimeSeriesなどの計算小技・注意事項
四則演算、指数
gwpyのTimeSeriesや FrequencySeriesは、 四則演算や指数の演算がそのまま使える。 横軸の配列(timesやfrequencies)が同一なら、TimeSeries同士の演算も可能。
FFT length の注意
TimeSeriesからスペクトルを計算する際、fftlengthの値が小さいと下図左のように謎の冪スペクトルに埋もれてしまうことがある。(ドリフト+窓関数の影響?)
十分なfftlengthをとればいいのだが、平均回数を稼ぎたいときや短時間のみを見たい場合にはそうもいかない。
そんなときは先に適当なhighpassフィルターをかけることで、下図右のように回避できる。