youtu.be いま、スーパーカブ110プロ（JA07）とCRF1100L Adventure Sports ES Dual Clutch Transmission（2BL-SD10）に乗っています。アフリカツインに乗り始めて、もうすぐ2年になろうとしています。 スーパーカブ110プロ（JA07） CRF1100L Adventure Sport…

昨日「相互相関の計算 - 丸井綜研」で「circshiftのあたりが面倒な気がしたので」と書きましたが、そんなに面倒ではなかったので書き残しておきます。 (defun circshift (arry k) "Circular shift the array ARRY by amount K. Shift towards right if K is …

相互相関 (cross-correlation) を計算するCommon Lispプログラムを作りました。とは言っても、畳み込み演算ができれば簡単に実装できます。Matlabのxcorr関数とそれなりの互換性があるものが欲しかったので、相互相関値だけでなくラグも返すようにしました。…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-06-22 13278.4 1.99 165.83 48.39 3.43 2022-07-27 13405.0 2.21 157.92 57.29 2.76 2022-08-31 13546.7 2.92 156.16 48.53 3.22 2022-11-05 13699.4 3.21 157.01 47.57…

「二項検定と二種類の過誤 - 丸井綜研」で、Common Lispで二項検定ができるようになりました。今日は1群のt検定を実装したいと思います。 t分布 t検定をするには、まずt分布の確率密度関数と累積分布関数を作る必要があります（検定だけなら累積分布関数だけ…

6月25日(日)にShibuya.lispのもくもく会がありました。14:30にオンラインで集合して、自己紹介と質問タイム（勉強になった！）があり、15時頃～17時頃にもくもくと作業をしていきました。（もくもく会の後の雑談は予定を1時間近くオーバーして盛り上がりまし…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 2860 12.74 160 18.21 8.78 2022-07-31 3128 13.63 157 19.66 7.99 2022-08-27 3401 16.63 157 16.42 9.56 2022-11-03 372…

リスト内の数値の合計を計算しようとして (defun sum (lst) (apply #'+ lst)) (sum '(1 2 3 4 5)) ;=> 15 というような関数を書きました。普通に使うぶんには全く問題なかったのですが、長いリスト（具体的には480,000個のdouble-float）を引数にしたら「ス…

先日、二項分布の累積分布関数を計算できるようになりました。 marui.hatenablog.com これを使って、簡単な二項検定をやってみたいと思います（あまりLISPとは関係ない内容ですが）。題材としては、ジャンケンにやたら強い人がいて「俺には天から与えられた…

統計に使える関数の実装を続けています。2月のLISPもくもく会では二項分布の確率質量関数（dbinom）を実装しました。今度は累積分布関数（pbinom）が欲しいところです。 marui.hatenablog.com 原理どおりに動くのであれば確率質量関数を累積加算すると累積分…

千葉でのキャンプのときに大荷物を乗せてのバイク走行が大変だったので、積載の方法を考え直さないといけないなーと思っているうちに、福島でのキャンプ予定が決まってしまいました。せっかくの大型連休ですしね。 走行の当日は雨になる予報だったので、東京…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-04-10 2283 14.23 163 19.40 8.40 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 2860 12.74 160 18.21 8.78 2022-07-31 3128 13.63 157 19.66 7.99 2022-08-27 340…

日本のAmazonで「Julia」と検索すると（「18歳以上ですか？」に「いいえ」と答えても）セクシーな画像が多く並んでしまうので、これから購入したい本のメモも兼ねて一覧を作りました。とは言っても出版された件数はまだまだ少ないので作成は楽ちん。リンクは…

youtu.be 給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-04-10 2283 14.23 163 19.40 8.40 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 2860 12.74 160 18.21 8.78 2022-07-31 3128 13.63 157 19.66 7.99 2022-…

以下のエントリに、二項分布の確率質量関数をCommon Lispで計算する関数を紹介したのですが、どうも意図した通りに動いてくれないときがありました。少し時間を作って調べてみたら、コードに間違いがあるのが明らかになりました。 marui.hatenablog.com Load…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39 3.43 2022-07-27 13405.0 2.21 157.92 57.29…

ごくたまにオーディオシステムの準備をすることがありますが、パワーアンプからスピーカへの結線は赤黒ケーブルが逆になってしまったり、左右スピーカへの接続を逆にしてしまったりということがあります。そこで、ステレオスピーカへの結線がうまくいってい…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-23 12910.1 2.55 158.04 37.57 4.21 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39…

【2023-04-17】このページの最下部に載っているコードでは、特定の引数の時に間違った答えが出力されたり、無限ループに入ってしまうことが分かりました。たとえば(dbinom 4 10 0.5) ;=>0.557となってしまいます（正しくは0.205）。元のアルゴリズムには問題…

車輪の再発明を続けています。今日は画像の書き出しです。 スペクトログラムを計算できるようになったものの、それを画面表示するためにはなんらかの画像にしないといけないんですね。画面にウィンドウを作って表示するよりもファイルにするほうが簡単なので…

Common Lispで振幅スペクトルの計算ができるようになり、Hann窓もかけられるようになったので、次はスペクトログラムの実装にいきたいと思います。 スペクトログラムは、振幅あるいはパワーの強さを時間と周波数ごとに計算したものです。以下の画像ではAudac…

macOSにはafで始まる名前の、オーディオファイルを操作するコマンドラインアプリが入っています（afclip、afconvert、afhash、afida、afinfo、afplayの6つ）。僕がよく使うのはafinfoとafconvertです。 afinfoは以下のようにオーディオファイルの情報を表示…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-23 12910.1 2.55 158.04 37.57 4.21 2022-03-06 13058.0 2.83 160.07 52.26 3.06 2022-05-26 13182.1 2.46 158.10 50.45 3.09 2022-06-22 13278.4 1.99 165.83 48.39…

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-02-05 1740 18.22 158 18.77 8.42 2022-03-07 2007 16.46 162 16.22 9.99 2022-04-10 2283 14.23 163 19.40 8.40 2022-06-26 2628 18.97 166 18.19 9.13 2022-07-18 286…

非負整数nの階乗n!はn×(n-1)×(n-2)×(n-3)×...×2×1と定義されています。また、0!=1とします。例えば6!であれば、6×5×4×3×2×1=720です。これを、整数以外にも拡張したのがガンマ関数で、nが整数の時にはΓ(n+1) = n!となります。例えばΓ(6)=120、Γ(7)=720となり…

目次 正規分布の確率密度関数 累積分布関数の近似式（Abramowitz and Stegun） 累積分布関数の逆関数の近似式（Abramowitz and Stegun） 累積分布関数と逆関数を誤差関数から計算する（Numerical Recipes） Common Lispを音だけでなく、統計にも使いたいです…

youtu.be キャンプ道具を揃えたシリーズの番外編です。本編は以下の3つ。 キャンプ道具を揃えた2022（バッグと小物類） キャンプ道具を揃えた2022（調理器具） キャンプ道具を揃えた2022（テーブルとチェア） 最初のキャンプ道具を揃えた2022（バッグと小物…

信号にフィルタをかけられるようになり、そして、レシピ通りのフィルタを作れるようになりました。 IIRフィルタをかける関数（Lisp Advent Calendar 2022） - 丸井綜研 双二次フィルタの係数（Lisp Advent Calendar 2022） - 丸井綜研 次に、フィルタの周波…

この記事はLisp Advent Calendar 2022の25日目です。Lispで音プログラミングに向けた、さらなる一歩。 目次 はじめに RBJ Cookbookによる低域フィルタ フィルタを使ってみる 他の形状のフィルタ はじめに 先日フィルタをかける関数を作りました（IIRフィルタ…

Julia Advent Calendar 2022の18日目の記事です。毎年Juliaで音に関係したことを書いていて、これまでにこのブログで書いたアドカレ記事は以下の通りです。 楽器音の分析とスペクトログラム (Julia Advent Calendar 2018) - 丸井綜研 パワースペクトル計算の…