雑感等

音楽,数学,語学,その他に関するメモを記す.

2019-02-01から1ヶ月間の記事一覧

TeX Live2018のセットアップ(Windows10)とソースの例

TeXWorksの設定(タイプセット) 名前:MAIN_upLaTeX (ptex2pdf) プログラム:C:/texlive/2018/bin/win32/ptex2pdf.exe 引数 -l -u -ot -kanji=utf8 $synctexoption -no-guess-input-enc -kanji=utf8 $fullname 名前:CITE_upBibTeX プログラム:C:/texlive/…

pythonでレスラー方程式の解軌道を描画(ルンゲクッタ法)

pythonでルンゲクッタ法を実装した. レスラー方程式を与えて描画する.プログラム内関数のdxdt, dydt, dzdtを書き換えれば別の微分方程式も解ける. 下記プログラムはローレンツ方程式を解くようにも書き換えられる. import numpy as np def dxdt(a, b, c,…

pythonで自己MIC(相互情報量)関数を計算

自己相関関数では,ある信号sig(t)とその信号を時間kだけ遅延させた信号sig(t+k)との相関係数を計算するが, ここでは相関係数の代わりに相互情報量(MIC)を求める関数automicを示す. MIC自体はminepyで計算する.下記プログラムを実行すると,グラフが2つ表…

matlabで振幅スペクトルの表示(fft)

matlabで高速フーリエ変換し,振幅スペクトルを表示する. function showfft(sig,fs) %sig:信号.1次元配列 %fs:信号のサンプリング周波数 N=length(sig); % 信号のサンプル数 plot(fs*(0:N-1)/N,abs(fft(sig)/N)) %振幅スペクトルの描画:縦軸が振幅,横軸が…