flickr2006/09/21行列の固有値を求めたい。なにをいまさら感たっぷりなんだけど、理系の出版社でテクニカルレビューっぽいことをしていると、特に統計の本なんかで行列の固有値を検算したい場合がなきにしもあらず。特性多項式を作ってそれを解いているらしく、うまく解が求まりやがらないことがある。実際、maxima/share/matrix/eigen.macを見ると、charpolyという関数で特定多項式を作り、それをsolveという関数で解いている。いう多項式を解く関数がうまくないといったことがMaximaのマニュアルにも書いてある。polynomial幸い、固有値を求めたい行列は対称行列ばっかりなので、自分でJacobi法をナイーブに実装してみた。対称行列だとうれしい理由やJacobi法については『プログラミングのための線形代数』の第5章で。2.5472427345288424)行列のデータをどういうふうに表現すべきなのか、ひどく悩んだ。本のようにベクトルを使うのがセオリーなんだろう。ベクトルなら行列のij要素を取り出したり変更したりするのも簡単だ。行と列を順繰りになめて各要素をちこちこ掛けたり足したりする処理ばっかり書くことになりそうで、ちょとブルー。そんなわけで今回はあえてリストのリストとして定義した(行ベクトルが1つのリスト。それを要素に持つ縦ベクトルのリストが行列)。リストなので、行列の積や和なんかも高階関数だけで定義できる。ms)))ちょっと満足。単位行列や回転行列を定義するのに要素を頭からなめるはめになってるんだけどね。(11)GaucheにもRubyみたいなARGFがほしい。正確にいうと、ほしがっているのは僕じゃないんだけ...荒れまくっていたウィキペディアの実数の項目が素敵な解説に直されて先週復活していた。編集履歴を見ると、...行列の固有値を求めたい。なにをいまさら感たっぷりなんだけど、理系の出版社でテクニカルレビューっぽい...ときどきむしょうにハンダ付けがしたくなる。サーバがアップデートされていた。ありがとうご...会社を休んで日本ソフトウエア科学会のPPLサマースクールを受講してきた。へとへとになって日暮里まで帰...いよいよネタ切れか。内部でプロシージャを定義するのを止めて、ついでに末尾再帰にだけしておこ...急がないと歯医者に間に合わないけど、ひきつづきクイックソート。partitionかあ。いちおう毎日ちょっとでも続けたいんだけど正直ネタ切れ。多値でクイックソート。やってみようと思った。出版される前の本の内容は、通常は著者や編集者に代表される「制作サイド」の人間にしか読まれない。専門性の高い本であれば査読とか監修といったプロセスを有識者にお願いすることはあるけど、そうしたお願いをするときには有償だったりカバーや袖に名前を出したりすることが多いので、これも「...職務経歴書いろいろあって職務経歴書を作成してみることにしました。機能を使うのが一番お手軽です。編集をお手伝いしました。Blogger.自己紹介Unknown詳細プロフィールを。