本日の給油(スーパーカブ110プロ/JA07)

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-08-31 13546.7 2.92 156.16 48.53 3.22 2022-11-05 13699.4 3.21 157.01 47.57 3.30 2023-02-06 13848.1 3.26 157.06 45.61 3.44 2023-03-02 13973.9 2.60 156.92 48.38…

信号の包絡線を計算する (Common Lisp)

信号の包絡線 (signal envelope) は「信号の頂点にシーツをかぶせたときのような概形」を持つ波形です。 MatlabやRで書いた信号包絡線の計算プログラムをCommon Lispに移植しました。信号包絡線の計算方法は振幅の絶対値に対して、(1)低域フィルタをかける、…

Aseprite買った

以前から欲しいと思っていたAsepriteというドット絵作成ツールが半額になっていたので、思わず購入してしまった。半額セールはたぶん2月13日まで。 www.aseprite.org で、初めて描いてみたのが以下の絵。16×16ピクセルのBrüel & Kjær Type 4128-Cさんを想像…

本日の給油(スーパーカブ110プロ/JA07)

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 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 3.30 2023-02-06 13848.1 3.26 157.06 45.61…

本日の給油(CRF1100Lアフリカツイン/2BL-SD10)

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-08-27 3401 16.63 157 16.42 9.56 2022-11-03 3727 17.20 158 18.95 8.34 2023-02-04 3987 16.75 157 15.52 10.12 2023-04-20 4165 10.37 157 17.16 9.15 2023-04-30 45…

キャンプ道具を揃えた2023(更新版)

キャンプ道具を揃えたシリーズの2023年更新版です。これまで購入したものは以下。 キャンプ道具を揃えた2022(バッグと小物類) キャンプ道具を揃えた2022(調理器具) キャンプ道具を揃えた2022(テーブルとチェア) キャンプ道具を揃えた2022(テント) そ…

ピンクノイズを生成するLISP関数

以前PythonやJSFX用のものも作ったのですが、ピンクノイズはいろいろなところで活躍するので、Common Lisp用にも作っておきます。 ピンクノイズを作るPython関数 - 丸井綜研 REAPERのJSFXでピンクノイズを生成してみた - 丸井綜研 Voss-McCartney法 DSP gene…

同じスペクトルを持つのにクレストファクタが違う波形(Julia Advent Calendar 2023)

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

アフリカツインにパニアケース購入

少しずつ快適なキャンプになるようにと考えていると、なぜかキャンプ道具が増えていきます。容量45リットルのシートバッグの上に、容量60リットルのソフトな防水シートバッグを載せて使っていましたが、車体重心が高くなってしまうことがネックでした。下の…

重回帰分析 (Lisp Advent Calendar 2023)

Lisp Advent Calendar 2023の第1日目の記事です。Lispについては全くのヘナチョコなので1日目は恐縮なのですが、カレンダーの先頭が空になっているのも幸先が悪いので……。 さて、車輪の再発明を続けて、ようやく逆行列の計算までたどりつきました。今回は逆…

逆行列の計算

Common Lispで逆行列の計算をするための準備をしてきました。 行列式の計算(Common Lisp) - 丸井綜研 LU分解(Common Lisp) - 丸井綜研 LUP分解(Common Lisp) - 丸井綜研 連立方程式を解く(Common Lisp) - 丸井綜研 今日はついに逆行列の計算です。と…

本日の給油(スーパーカブ110プロ/JA07)

給油日 オドメーター (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…

連立方程式を解く

先日以下の記事にしたように、LUP分解をすると正則行列𝐴に対して𝑃𝐴 = 𝐿𝑈という分解ができるのでした。𝐿は下三角行列、𝑈は上三角行列、𝑃は並べ替えの情報が入っています。 marui.hatenablog.com これを使って連立方程式𝐴𝑥 = 𝑏を解くプログラムを作ります。LU…

LUP分解

【2023-11-09:入力行列を書き換えながら処理をするようになっていたので、呼び出し元でも行列の内容が書き変わってしまっていました。複製を作って作業するように修正しました。】 昨日LU分解のコードを掲載しましたが、単純なアルゴリズムなので、𝑎₁,₁要素…

LU分解

10月29日(日)にShibuya.lispのもくもく会がありました。残念ながら雑談の途中で中座しないといけなくなってしまいましたが、とても楽しかったです。この日のもくもく会の目標は「逆行列を計算できるようにすること」でした。実際は途中までしかできなかった…

Pure Dataのexternalを作ってみた(Apple Silicon対応版バイナリを作る)

何年か前の夏休みにPure Dataのexternalオブジェクトを作るチュートリアル記事を書きました。以下の2篇です。 marui.hatenablog.com marui.hatenablog.com その当時はIntelチップ搭載のMacに向けて書きましたが、もう最近のmacOSは32 bitをサポートしていま…

行列式の計算

Common Lispで車輪の再発明をし続けています。いろいろと統計に使う関数を作っていましたが、そろそろ逆行列を計算したくなってきました。逆行列と固有値は線形代数の基本ですものね。*1 逆行列の計算をする前には、行列式(determinant)という、逆行列が計…

本日の給油(CRF1100Lアフリカツイン/2BL-SD10)

youtu.be 給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-08-27 3401 16.63 157 16.42 9.56 2022-11-03 3727 17.20 158 18.95 8.34 2023-02-04 3987 16.75 157 15.52 10.12 2023-04-20 4165 10.37 157 17.16 9.15 2023…

本日の給油(スーパーカブ110プロ/JA07)

給油日 オドメーター (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…

上方倍音列と下方倍音列

とある音楽理論の一派では、自然倍音列にヒントを得て作曲に活かしたりするそうです。そんな倍音列の周波数と音名との対応を見られるJuliaプログラムを書いてみました。そもそも「倍音ってなに?」については以下の記事を参照してください。 marui.hatenablo…

コラッツ問題

コラッツ問題という数学の未解決問題があります。自然数nが偶数のときn÷2を計算、nが奇数のとき3n+1を計算し、その結果に同じルールを適用し続けると、最終的に1に至る(あるいは最終的に1→4→2→1の無限ループになる)、というものです。すべての自然数につい…

本日の給油(CRF1100Lアフリカツイン/2BL-SD10)

youtu.be 給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2022-08-27 3401 16.63 157 16.42 9.56 2022-11-03 3727 17.20 158 18.95 8.34 2023-02-04 3987 16.75 157 15.52 10.12 2023-04-20 4165 10.37 157 17.16 9.15 2023…

カレーのグレイビーを省力化チートで作った

印度カリー子さんのスパイスカレーを知ってから、カレー作りが楽しいです。カレー・グレイビーの基本レシピがやたら覚えやすいですし(ニンニク1片、ショウガはニンニクと同じくらい、玉ねぎ1個、トマト1個、ターメリック・クミン・コリアンダー・塩が各小さ…

cl-portaudioを使った音の入出力

前回の記事でCommon Lisp上でPortAudioが使えるようになりました。 marui.hatenablog.com 今回は音の入出力をやってみます。基本的にはcl-portaudioのドキュメントに載っている「入力と出力を15秒間だけつなげる」というデモプログラムを実行するだけです。…

LispからPortAudioでデバイス情報を得る

Common Lispからcl-portaudioを使ってPortAudioを呼び出して音を再生したいのですが、まずはデバイス情報を得るところまでやりました。 cl-portaudioを利用するには例によってQuicklispを使うと便利みたいです。 (ql:quickload :cl-portaudio) とはいえ、bod…

金属の熱伝導率

キャンプ道具には鉄、ステンレス、アルミニウム、チタンなどの金属が使われています。「チタンは熱伝導率が低いから料理には向かない」などと言われることもあったりするので、どういうことなのか調べてみることにしました。 手元に少し前の理科年表(2020年…

国内4大メーカーの二輪車販売台数

国内で二輪車を製造している主要4企業の2022年度の販売台数を調べてみました。二輪車の販売台数順に並べましたが、ホンダとスズキは四輪車事業も付記してあります。ホンダはインドネシアとかインドとかで強いので小型二輪あたりが多いのでしょうが、国内外あ…

アフリカツインを入手するまでの話

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関数とそれなりの互換性があるものが欲しかったので、相互相関値だけでなくラグも返すようにしました。…