9つのものを袋に分ける分け方は何通りあるか？という問題について、「数え上げ理論」本に書いてあったアルゴリズムで、友人がJAVAで書いてくれました。答えは手で数えたものと合っています。とりあえずこのままJ言語に翻訳してみようと思う(予定)。 int p(in…

数学の本でベクトルというと、 ● 向きを持った量(矢印)のこと ● 1行もしくは1列に並べて書かれた順序をもった実数の組 という2種類の定義があって混乱していました。本を読むとたいていどちらかしか書いてない。ぜんぜん関係ないのかと思うとそうでもないよ…

「急がばまわれ式・堅実で一番効率的な英語の勉強法」という記事が話題になっています。 内容はとてもまっとうなことが書いてあるのですが、冒頭の自己紹介で、TOEIC975点と書いてあって、私と同じ点数なので自分が書いたのかと思いました。 先日も朝コーヒ…

まず最初に言葉の定義ですが、頭の中に絵を描いてください。 正方形の上に三角が載った家の形です。正方形の対角線を結んでバッテン(×)を描いてください。 この絵で点が5個あります。ちゃんと数えてくださいね。 辺の数はまんなかのバッテンを入れて8本です…

「数え上げ理論」というブルーバックスの真ん中あたりを開いたら、9つのものを袋に分ける分け方は何通りあるか？という問題が載っていました。 本を真ん中から読むというのはよくやります。 推理小説の最初と最後だけ読むというのもたまにやります。 9つのも…

表題のようなフレーズが頭に浮かびました。 今の人生は、犬だったころにまじめに生きたことのごほうび 輪廻転生とか信じているわけではありませんが、かと言って神とか魂の存在を否定しているわけでもない。 世界には神様を信じている人の数の方が信じていな…

アマゾンでシャープのPC-Z1を購入しました。 今まで富士通のLOOX Uにつないでいたイーモバイル(D02HW)をつないで使っています。 はじめちょっとトラブリましたがその後快適。 一日中どこにいても2秒でインターネットにつながります。 画面も広いし、読みやす…

江戸川土手のサイクリングロードがあります。葛西臨海公園駅スタートでどこまでさかのぼってもいいのですが、関宿を目的にする例がネットにありました。56キロを往復するのはわたしのようなヘタレサイクリストにはつらいので、輪行を考えます。ちょうど真ん…

とある本で行列の積の定義があって、例としてつぎのような行列A、Bの積を行っていました。 B 0 3 1 4 2 5 A 6 8 7 9 B mp A 21 27 34 44 47 61 さて、通常J言語で3x2(3行2列)の行列を作ろうと思うと、次のようにします。 i.3 2 0 1 2 3 4 5 これが例のBと微…

ベクトルa = (a1, a2)とベクトルb = (b1, b2)の成す角度θを求めよ。 ベクトルの内積(dot product)はa1*b1 + a2*b2である。 これをa dp b = a1*b1+a2*b2と書くものとする。 ベクトル自身の内積はベクトルの長さの二乗になる。 a dp a = a1*a1 + a2*a2 b dp b …

固有値を解くのが楽しいので、ネットで「固有値 例題」でヒットしたサイトの例題を解いてみます。 http://next1.msi.sk.shibaura-it.ac.jp/MULTIMEDIA/linearalg02/node14.htm ]a=:3 3$3 0 0 0 2 _5 0 1 _2 3 0 0 0 2 _5 0 1 _2 3 Naigen a _3 1 _3 1 p. 3 N…

中野先生の固有値の論文を順々に読んでいったら、ル・ヴェリエ、ファデーエフ法というのが出てきました。 これがとても優秀のようです。 とりあえず転載。 NB. Leverrier-Faddeev Algorithm char=: 3 : 0 X=.I=.=@i.n=.#y[p=.1 for_k. >:i.n do. X=.y +/ . *…

固有値問題の実例をやってみます。 ]a=:4 4$1 3 5 7 1r3 1 1 5 1r5 1 1 3 1r7 1r5 1r3 1 1 3 5 7 1r3 1 1 5 1r5 1 1 3 1r7 1r5 1r3 1 この4x4の行列はExcelで学ぶAHP入門という本の最初の例で、スポーツクラブを選択する際に、費用、施設・環境、交通の便、…

BD-1という自転車があります。日記を調べたら、女房が2002年に購入していた。 乗ってもよい、というのでたまに乗るのですがとても乗り心地が良い。 見かけは風変わりな自転車で、折りたたみだし、車輪がこれは18インチかな。 いわゆるミニベロというのでしょ…

自分がいまどんな風に考えているかの魚拓です。 経済問題はあとになってみると以前どう考えていたかがわからなくなるので、定期的に魚拓を取っておく必要があります。 さて、10000円前後、89円前後の相場が続いて株高に振れるでもない、円高に振れるでもない…

dcde ged2 egag8a8 <d>bag ega2 <dcd2> egag e.g8d2 a<cd2> <cd>ag ag8e8d2</cd></cd2></dcd2></d>

つまらないことでひっかかっていましたが、ようやく動くようになりました。 すごい！ あとで実例はお見せしますが、下記のプログラムで多次元の固有値問題を解いてしまいます。 ひっかかっていたのはプログラム中にタブがあるとエラーになってしまうことに気…

数学や統計の本を読んでいると固有値というのが出てきますが、APL/J言語でどうやって解くのだろうと思っていました。英語ではなぜかドイツ語みたいな単語でeigenvalueと呼ぶらしいのですが、英語でも日本語でも調べていてしっくりこないでいたところ、ネット…

ビクターが音楽レッスン用のボイスレコーダーのようなものを発売しました。 ポータブルデジタルレコーダーXA-LM1希望小売価格26250円 高音質レコーダー+クロマチックチューナー+メトロノーム で、購入予定なのですが、おもしろいところを突いてきたなと思っ…

J言語では有理数をあつかえてとても気持ちがいいのですが、このときx:(エックスコロン、extended precision)という動詞を使います。 2 % 3 0.666667 x: 2 % 3 2r3 ここで2r3とかが分数(fraction)です。分数のままで普通の計算ができます。自動的に約分されて…

(=/; +./; *./; ~:/; -.@=/; =&.-./)~ 0 1 +---+---+---+---+---+---+ |1 0|0 1|0 0|0 1|0 1|0 1| |0 1|1 1|0 1|1 0|1 0|1 0| +---+---+---+---+---+---+ =(イコール)は同じかどうかのテストですのでブール代数としてもそのまま使える。 +.(プラスドット)はO…