TAIZO(博多)のラーメン

黒い麻油がかかったこってり黒ラーメン
博多に来て一番美味しかった。

pythonを使って固有値・固有ベクトルを求める

Pythonのnumpyとscipyを使って固有値・固有ベクトルを指定した数分、作成するコードを作ってみたのでメモ書き

　
 #encoding:UTF-8

import numpy
import scipy.linalg

hi = 2
lo = 0
#行列のサイズ
w = 3
list = []
for i in range( 10 ):
    list.append(i)

def EigenValue():
    global list
    a = []
    #行列
    a = []
    for i in range(w):
        if i == w:
            break
        else:
            a.append(list[i:i+w])
    A = numpy.matrix(a)
    eigen_value,eigen_vector = scipy.linalg.eigh(A,eigvals=(lo,hi))
    eigen_id = numpy.argsort(eigen_value)[::-1]
    eigen_vector = eigen_vector[:,eigen_id]
    eigen_value = eigen_value[:,eigen_id]
    print eigen_value
    print eigen_vector

if __name__=="__main__":
    EigenValue()

Rのおぼえがき

Rの勉強してて、手順を覚えがき

• 列ごとに平均値を出す

apply(dt[,c(-1:-3)] , 2 , mean)

• 行列操作の整理

コマンド	説明
x[2, ]	2 行目を取り出す
x[ ,2]	2 列目を取り出す
x[1, 2]	1 行 2 列目の成分を取り出す
x[c(1,2), 2]	1,2 行 2 列目の成分を取り出す
x[c(1,2), c(1,3)]	1, 2 行目と 1, 3 列を取り出す
x[ ,-c(1, 3)]	1, 3列を除外した行列を取り出す
x[ , c(T,F,T)]	1, 3列を取り出す
rowSums()	行の総和
colSums()	列の総和
rowMeans()	行の平均
colMeans()	列の平均

Rでお勉強（行列のsum(A)って？）

とある本を読んでいてRの操作で記述されているsum(A)って全要素の和でいいんだっけ？と思い、実際にRで試してみた。

勉強もかねて、テキストファイルからデータを読み込むところもおまけで

まずは、データを読む

> A <- read.table("/home/toku/r_test.txt")

表示

> A
  V1 V2 V3 V4 V5
1  1  2  3  4  5
2  6  7  8  9 10

行と列の最初に、でてるのは、ただの表示ですね。

和の実行
> sum(A)
[1] 55

ふむ、全要素の和ですね。

プランク温度

突然思い出したように物事が不思議になる。ここ一ヶ月くらいは、温度の上限がテーマです。「なにが光速で動いているか」がわからない。そんな世の中を暮らすのに１００％役に立たないことを考えるのが好きです。

wikipediaより

プランク温度は、プランク単位系における温度の単位。

プランク温度T_Pは、次のように定義される。

ここで、m_Pはプランク質量、cは真空中の光速度、はディラック定数、kはボルツマン定数、Gは万有引力定数である。

SI単位では1.41679 × 10³² ケルビンとなる。

1プランク温度は、ビッグバンの瞬間から1プランク時間経過したときの宇宙の温度である。1プランク温度以上の温度で物理的に意味のあるものは知られていない。

ふーん。よく分からない。。。

Rをシェルスクリプトで実行してみた。

スゴく簡単なんだけど、いろいろ書いてシェルスクリプトから実行することを想定して練習してみる。以下を記述して実行。

#!/bin/sh
R -e "sqrt(100);\
round(100);\
log10(100)"

Rのインストール（Ubuntu）

• /etc/apt/sources.listに以下を追記

deb http://cran.md.tsukuba.ac.jp/bin/linux/ubuntu lucid/

• Ubuntu レポジトリの公開鍵を入手

 $ gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9

 $ gpg -a --export E084DAB9 | sudo apt-key add -

• R 本体のインストール

 $ sudo apt-get update

 $ sudo apt-get install r-base

あとは、Synapticで必要そうなものを『r-cranで検索』してインストールする。