2023-11-01から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のもくもく会がありました。残念ながら雑談の途中で中座しないといけなくなってしまいましたが、とても楽しかったです。この日のもくもく会の目標は「逆行列を計算できるようにすること」でした。実際は途中までしかできなかった…