MATLAB：正規乱数

研究室の輪読でPRMLを使っている．

原著↓

Pattern Recognition and Machine Learning (Information Science and Statistics)

日本語版↓

パターン認識と機械学習 上 - ベイズ理論による統計的予測

パターン認識と機械学習 下 - ベイズ理論による統計的予測

で，そろそろ数式を追いかけるだけじゃなく，プログラムとして実装しようと

第２章のベイズ推論でガウス分布の平均値のガウス分布(ややこしい．．)をプロットするというものを

MATLABでコーディングしていた．

ベイズ推論自体は数式をforとかifとかで書き下せばいいんだけど，

観測データの平均μ，分散σ^2の正規乱数をどう与えればいいのか分からなかったので調べた．

• • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • -

N = 10;% データ数
Avr = 0.8;% 真の平均
Var = 0.1;% 真の分散
x = sqrt(Var).*randn(1,N)+Avr;% 観測データ行ベクトルxを設定

• • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • -

なんでrandnの前の演算子が.*になっているのか分からない．．

別に*でも結果一緒やん．

あ，今思いついたけど，Varが共分散行列って場合も考えてるのかも．．

明日誰かに聞いてみよう！