Rでスペクトル解析。 正弦波以外の周期をもった波が倍音が出力されることに注意 例１）正弦波 ３０日周期 x*1 例２）３０日周期のイベント（倍音がでる） x plot(x) spec.pgram(x,c(10,10)) https://jp.mathworks.com/help/signal/examples/practical-introd…

Rでガウスの超幾何関数を計算する。 hypergeo 関数は一部パラメータでバグがあるのでgslの超幾何関数でおきかえる。 ただし、gslの超幾何関数は複素数に解析接続されていないので、そこはhypergeoを使う。 library("hypergeo") library(gsl) library("BAS") …

RでHurwitz's zeta関数 フルビッツゼータ関数の近似 よいライブラリがないので、とりあえず、つくってみた。 library("Rmpfr") Hurwit_zeta_approx<-function(beta4,x){ #N<-10^7 #k<-0:N #s<-beta4 ans2<-NULL if(beta4==beta4){ for(j in 1:length(x)){ a3…