物理学者(ポスドク)による日々の研究生活のメモ書きです ( python/emacs/html/Japascript/シェルスクリプト/TeX/Mac/C言語/Linux/git/tmux/R/ポケモンGO)
sin波の位相と振幅を調べる必要があった
振幅はsin波のmaxを調べればいい
→np.argmaxを使えばOK
位相はzero crossする時刻から割り出せばよいのかな?と思った
ひとまず、zero crossするindexを割り出す方法を調べてみた
■ 参考 : Efficiently detect sign-changes in python
ここに色々とまとまってた
とりあえずメモだけしておく
ただ、この方法だと、+から-に0を交差する点と-から+に0を交差する点を区別できない
そういうときは具体的に求まったindexでのdataの値を見てみれば良い
それをwhereを使って+から-だけ残したのが以下のコード
mp.whereの使い方は以下
■ 参考 : NumPyで条件に応じた処理を行うnp.whereの使い方
whereの使用例としては、位相を-180~180度の間に制限したいときとかは以下のような感じ
# convert the range [-180:180]
#b= np.where(b<-180.0, b+360.0, b)
#b = np.where(b>180.0, b-360.0, b)
np.diffは微分をするためによく用いられるらしい
逆に積分をするcumsumという関数もある
■ 参考 : 要素の差分、足し合わせを計算するNumPyのdiff関数とcumsum関数の使い方
ツイート
振幅はsin波のmaxを調べればいい
→np.argmaxを使えばOK
位相はzero crossする時刻から割り出せばよいのかな?と思った
ひとまず、zero crossするindexを割り出す方法を調べてみた
■ 参考 : Efficiently detect sign-changes in python
ここに色々とまとまってた
とりあえずメモだけしておく
import numpy as np
a = np.array([0, -2, -1, 0, 1, 2, 0])
zero_crossings = np.where(np.diff(np.signbit(a)))[0]
print(zero_crossings)
# output: [0 2]
signで符号を取り出して、diffで前の符号との差分を調べているa = np.array([0, -2, -1, 0, 1, 2, 0])
zero_crossings = np.where(np.diff(np.signbit(a)))[0]
print(zero_crossings)
# output: [0 2]
ただ、この方法だと、+から-に0を交差する点と-から+に0を交差する点を区別できない
そういうときは具体的に求まったindexでのdataの値を見てみれば良い
それをwhereを使って+から-だけ残したのが以下のコード
index = np.where(np.diff(np.sign(data) > 0))[0]
index2 = index_zc[np.where(data[index] > 0)]
index2 = index_zc[np.where(data[index] > 0)]
mp.whereの使い方は以下
■ 参考 : NumPyで条件に応じた処理を行うnp.whereの使い方
whereの使用例としては、位相を-180~180度の間に制限したいときとかは以下のような感じ
# convert the range [-180:180]
#b= np.where(b<-180.0, b+360.0, b)
#b = np.where(b>180.0, b-360.0, b)
np.diffは微分をするためによく用いられるらしい
逆に積分をするcumsumという関数もある
■ 参考 : 要素の差分、足し合わせを計算するNumPyのdiff関数とcumsum関数の使い方
PR
この記事にコメントする
プロフィール
HN:coffee
職業:物理屋(自称)
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
職業:物理屋
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
カテゴリー
ブログ内検索
リンク
相互リンク募集中です
(Twitterにてお知らせください)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
最新記事
(11/20)
(05/17)
(05/17)
(05/09)
(03/05)
(02/29)
(02/21)
(02/21)
(02/21)
(02/21)