Skip to main content

マルコフ行列の中の著者達: どの著者がもっとも人々に影響を与えたのか? (17)


前回は隣接行列の拡張として Markov matrix を導入した.ここで具体的な Markov matrix \(M\) の例を以下のように考える.
\begin{eqnarray*} M &=& \left[ \begin{array}{cc} 0.8 & 0.3 \\ 0.2 & 0.7 \\ \end{array} \right] \end{eqnarray*} \(M\)の要素の意味を具体的に書くと,以下のようになる.
\begin{eqnarray*}  M &=&  \left[   \begin{array}{cc}    \mbox{Stay Berlin}       & \mbox{P $\rightarrow$ B} \\    \mbox{B $\rightarrow$ P} & \mbox{Stay Potsdam}      \\   \end{array}  \right] \end{eqnarray*}
ここで,\(\mbox{B $\rightarrow$ P}\) は Berlin から Potsdam に引っ越す人の割合,\(\mbox{P $\rightarrow$ B}\) は Potsdam から Berlin に引っ越す人の割合を示す.つまり,一年たって,Berlin に 8 割の人がそのまま Berlin に住み,2 割の人は Berlin から Potsdam に移動する.全ての人に関して考えているので,カラムは合計 1となるように(0.8 + 0.2 = 1.0) なっている.Potsdamの場合も同じである.7割の人は一年後も Potsdam に住み,3割の人がBerlin に引っ越す.この場合にも,カラムは合計 1 となっている.

「何割の人が移動する」ということを示す Matrix なので,要素は全て 0 以上であり,1 以下である.たとえば,マイナスの割合の人が移動するということはない.また,カラムの合計が 1 になることは既に見たが,それは住民全員を考えているからである.これが 1 を越えることがないのは,引っ越す人間の合計が住んでいる人の合計よりも多くなることはないからである.この行列では,今年の状態から来年の状態が完全に決定する.このような性質を持つ Matrix を Markov matrix と呼ぶ.

ここでちょっといろいろと予想してみたい.たとえば,最初に Berlin に 900人の人が住んでいて,Potsdam には 100 人の人が住んでいるとする.合計1000人の人がこの2つの街には常にいるとしよう.一年たったら,人口がどのように推移するか octave で計算してみる.

octave:29> M = [0.8 0.3; 0.2 0.7];
octave:30> p = [900 100]';
octave:31> M * p
   750
   250

二年後には次のようになる.

octave:32> M^2 * p
   675
   325

一年,二年,というのを考えてみたが,では,長い年経た後はどうなるであろうか,ちょっと考えてみたい.とはいえ,何故長い年月を考えるのか疑問に思う読者もいるだろう.長い年月の後というのはつまり未来はどうなるかということである.私はこれに興味がある.

多くの科学がそうであるのだが,知識を集める目的は「未来を知りたい」からである.数学が有用であるのは,数学が未来を知る助けになるからである.人間が未来を知りたいのは,それによって生き延びる可能性を広げてきたからではないかと私は思う.未来などどうでもよいという人々もいるだろうが,そういう人々で構成されている社会は滅ぶ可能性が高い.例えば,冬の食料を心配せずに消費してしまう人々はその年の冬に滅んでしまう.未来に備えて子供達を教育しない社会は発展する可能性が低いと思う.我々が長いこと生き延びてきたのは,未来に興味を持ち,それに対して備えることを学んだからであろう.そうでない人々は生きのびてこなかったので,結果として未来に興味を持つ人々が多いのは自然のように思う.

次回はこれまでの仮定から未来を予想してみよう.

Comments

Popular posts from this blog

共有メモリによるプロセス間通信

Unix の共有メモリを使ったプロセス間通信について調べて実験をしてみた.対象は1つのホスト上での複数のプロセスである.ネット上でいくつか例題はないかと探したが,どうも良い例となるコードが見当たらなかった.結局はある解説記事と,Stack Overflow の議論と,man page を見て作ってみたものになったので,例をここに置くのも有用かと考え,この記事を書く.(もしかしたら探し方が悪くて良いコード例をみつけられなかっただけかもしれない.) mmap を使うかどうかという話がいくつもでていたが,POSIX の方向としては,shmem_open と mmap を使うという方向があるということだったので,それを信じてその形での実装を試してみた. 基本的なコードの流れは次のようになる. 共有メモリ領域を1つのプロセスが shm_open() を使って作成する.その際に,プロセス間で共通の文字列を識別子(``identifier'')とする.(Linux ではこれが /dev/shm/identifier のように見える.) 共有メモリ領域を mmap() でメモリにマップする.共有メモリポインター (shared_ptr)が得られる. shared_ptr を使って複数のプロセスで通信をする. 利用終了後は munmap() をつかってマップを消す. 共有メモリオブジェクトを shm_unlink() によって消す. 以下に示すプログラムは,server と client の2つのプロセスが共有メモリを使って通信をするものである.ここで,server プロセス数と client プロセス数は共に 1 を仮定する.server と client は自分の領域にしか値を書き込まないことで,ロックを避けている.互いに相手の値を読み,それよりも1大きい数を一定の期間ごとに自分の領域に書くという例題である.シンプルではあるが,共有メモリで通信をする基本としては十分なものだと思う.ソースコード(shmem_test.cpp)を以下に付加する.ソースコードのコメントにコンパイル方法とどのように利用するかを書いておく. /*   Shared memory inter process communication minimal exa

複数の線を持つ線グラフを Jenkins の plot plugin で描く方法

私は毎夜のソフトウェアテストを自動化するために Jenkins というツールを使っています.今回は, valgrind  を使ってメモリーリークのテストを自動化することにし ました.その際,エラーの数の結果をグラフとして表そうと思って, Plot plugin  を使うことにしました. Plot plugin の例図からは,複数のデータラインを描くことができるのは明らかなのですが,どうやったらいいのかは参照のページや,例としてあった Perl script,plugin 中の help からは私にはよくわからなかったのです. ここで重要な考えは,それぞれのデータラインにはそれぞれの出力ファイルが必要ということでした.私はこれを誤解していました. 例として,ビルドの時に次の property データファイルを出力します.それぞれのファイルが1つのデータラインを表します. valgrind_trunk_result.definitely.property valgrind_trunk_result.indirectly.property valgrind_trunk_result.possibly.property それぞれのデータの中身は1行のデータ点です.たとえば, valgrind_trunk_result.definitely.property ファイルの中身は次のような1行 です. YVALUE=0 このファイルを ${WORKSPACE} ディレクトリ以下に出力します.ここで," WORKSPACE " は jenkins が提供する環境変数です. 図1が私の plot plugin の設定を示しています.これは jenkins の config 画面です.3つの data series があって,それぞれにデータファイルがあります. Figure 1: Plot plugin configuration in Jenkins 図2が結果です.複数の線が描かれているのがわかります.(実際には 3 本の線がありますが,最初の線と2番目の線が同じデータなので,重ねって見えません.) Fugure 2: Plot data with multiple data lines

ソニーのカメラ (α 5000) の 30 分のビデオ録画時間の制限を外す方法

私は Sony の Alpha 5000 を気にいって使っています。しかし一つだけ問題がありました。それはビデオの録画時間の制限が 30 分というものです。 今日,ちょっと気になって探したらこの制限を解除できることがわかりました。以下のビデオがその紹介です。 https://youtu.be/7cstA_PuRIg このビデオの作者によれば,ほとんどのソニーのカメラのビデオの制限はなくせるそうです。ただし私が試したのは,Alpha 5000 のみです。 手順 カメラ側 スイッチ On Menu -- Setup --- USB connection を MTP にする スイッチ Off and On USB ケーブルでカメラをコンピュータに接続する (以下接続したままにする) コンピュータ側でソフトのダウンロードとインストール (私は Windows 10 で試しました) 次の URL に行く https://sony-pmca.appspot.com/apps ただし,Internet Explorer か Safari のみサポートということでした。Chrome では上手くいきませんでした。私が試したのは Windows 10,Internet Explore 11 です。 注意事項: このサイトは Sony のサイトですが,ここにあるソフトウェアは Sony のものとは限らないので保証はありません。御自分でリスクを判断してご利用下さい。当方も何も責任を負えません。 上記の URL から,OpenMemories のページに移動する。 このページにある PMCADownloader plugin (PMCADownloader.msi) をダウンロードする PMCADownloader をインストールする 私はいちどここでページを閉じてもう一度 https://sony-pmca.appspot.com/apps を開き,OpenMemories のページに移動しました ここで log に Loading plugin Plugin loaded と表示されます。PMCADownloader の Install がされていない時には,``Plugin loaded'&#