読者です 読者をやめる 読者になる 読者になる

Python再開してみています

Python

15年くらい前には中途半端にPythonを使っていたのですが、最近のデータサイエンス業界でのPythonの盛り上がりが気になったので、Pythonに再入門してみました。(15年前は主としてCとJavaをやっていたのでCythonではなくJythonを使っていましたが、いろいろと不満もあってMatlabに移行したのでした)

ここ1〜2週間ほどPythonSoloLearn.comのウェブサイトとiPhoneアプリで勉強しつつ、これまでMatlabやRで書いたコードをPythonに移植しています。

全てのものが一カ所にまとまっているMatlabと違って、PythonではNumPy、SciPy、Matplotlibなど様々なライブラリに散らばっていて面倒だと感じましたが、とりあえずpyenvでanacondaを入れておけば最小公倍数的に必要そうなものが全部そろいますし普段のMac OS X環境を破壊してしまう心配がないので良いと思います。今後、本腰を入れてPythonを使うかどうかは分かりませんが、そうであれば必要なものだけに絞り込んでインストールしたい気がします。

さて、今日はオーディオファイルを読み込んでスペクトルなどを表示するところまで行きました。以下、撥弦楽器のサンプルを分析したものです。

import soundfile as sf
x, fs = sf.read("sound.wav")

# 自前のツール群を読み込み
exec(open("./audioutils.py").read())

# スペクトルを描画(両線形軸)
plot_spectrum(x, fs)

f:id:amarui:20170323131642p:plain

# スペクトルを描画(両線形軸・範囲指定)
plot_spectrum(x, fs, xlim=(0,5000))

f:id:amarui:20170323132810p:plain

# スペクトルを描画(両対数軸)
plot_magnitude(x, fs)

f:id:amarui:20170323132118p:plain

# スペクトル中心(Hz)の時間変化を描画
sc, t = spectral_centroid_moving(x, fs, window_size=4096, hop_size=1024)
import matplotlib.pyplot as plt
plt.plot(t, sc)
plt.xlabel("Time (sec)")
plt.ylabel("Spectral Centroid (Hz)")
plt.grid()
plt.show()

f:id:amarui:20170323155424p:plain

# 実効値(dB)の時間変化を描画
r, t = rms_moving(x, fs, window_size=4096, hop_size=1024)
plt.plot(t, r)
plt.xlabel("Time (sec)")
plt.ylabel("RMS (dB)")
plt.grid()
plt.show()

f:id:amarui:20170323161423p:plain

plt.show()の代わりにplt.savefig("hoge.pdf")などとすると表示せずファイルに保存。

本日の給油

バイク
給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km)
2017-03-18 8559.6 2.76 125.72 45.94 2.74

今日は埼玉県南埼玉郡まで。前回給油が去年10月29日に茨城県つくば市に行ったときだったので、それからは近所への買い物くらいしかしてなかったことになる。真夏用のメッシュジャケットしか持っていないというのが原因なのかも? 来冬には冬用のジャケットやグローブ買おうかな。

本日の給油

バイク
給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km)
2016-02-13 6714.7 2.96 108.78 49.26 2.21
2016-04-03 6847.3 2.57 108.95 51.60 2.11
2016-04-10 7025.7 3.42 111.99 52.16 2.15
2016-05-18 7113.8 1.68 108.93 52.44 2.08
2016-06-01 7286.3 3.28 114.94 52.59 2.19
2016-06-15 7464.2 2.65 118.11 67.13 1.76
2016-07-24 7579.7 2.54 107.87 45.47 2.37
2016-09-09 7740.9 2.87 127.87 56.17 2.28
2016-10-08 7878.6 2.88 105.90 47.81 2.21
2016-10-09 8075.8 3.13 112.14 63.00 1.78
2016-10-13 8238.7 2.89 127.68 56.37 2.27
2016-10-29 8432.8 3.52 123.86 55.14 2.25

4回分まとめて。

48ヶ月点検

バイク

スーパーカブ110プロ(JA07)の4年目の定期点検。春頃から、以前と音が変わったかなーと思ったまま乗り続けていたところ、それがチェーンにけっこうなダメージを与えていたらしいとのこと。今年は調整で何とかしてもらい、来年の点検時にはスプロケット込みで交換の可能性が出てきました。ぐむむむ。

これまでの走行距離はだいたい8000 km。来年の今頃は10,000 kmくらいか。1万キロでチェーン交換なんて聞いたことないよー(自業自得)

本日の給油

バイク
給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km)
2016-02-13 6714.7 2.96 108.78 49.26 2.21
2016-04-03 6847.3 2.57 108.95 51.60 2.11
2016-04-10 7025.7 3.42 111.99 52.16 2.15
2016-05-18 7113.8 1.68 108.93 52.44 2.08
2016-06-01 7286.3 3.28 114.94 52.59 2.19
2016-06-15 7464.2 2.65 118.11 67.13 1.76
2016-07-24 7579.7 2.54 107.87 45.47 2.37
2016-09-09 7740.9 2.87 127.87 56.17 2.28

急にガソリン価格が上がったなー。