シミュレーション物理7 乱数.

Slides:



Advertisements
Similar presentations
画像処理・実習 第十四回:パターン認識 東海大学 情報理工学部 情報メディア学科 濱本和彦. 今回の内容 5. パターン認識 5.1 マッチングの原理 5.2 テンプレートマッチング 実習 相互相関とテンプレートマッチング.
Advertisements

シミュレーション論Ⅰ 第 7 回 待ち行列のシミュレーション(2). 第 6 回のレポート(解答例) 乱数表より乱数を記入し、到着間隔・サービス時間にした がってグラフを作成する 例) 最大待ち人数:2人 最大待ち時間:5分 平均待ち時間:3分.
エクセルと SPSS による データ分析の方法 社会調査法・実習 資料. 仮説の分析に使う代表的なモデ ル 1 クロス表 2 t検定(平均値の差の検定) 3 相関係数.
ロボットビジョン(ロボットの視覚能力)のための デジタル画像処理
シミュレーション論Ⅰ 第6回 待ち行列のシミュレーション.
1.1 C/C++言語 Hello.ccを作りコンパイルしてa.outを作り出し実行する
運動方程式の方法: 惑星の軌道 出席のメール(件名に学生番号と氏名)に,中点法をサブルーチンを使って書いたプログラムを添付
第4回 (10/16) 授業の学習目標 先輩の卒論の調査に協力する。 2つの定量的変数間の関係を調べる最も簡単な方法は?
心理統計学 II 第7回 (11/13) 授業の学習目標 相関係数のまとめと具体的な計算例の復習 相関係数の実習.
大数の法則 平均 m の母集団から n 個のデータ xi をサンプリングする n 個のデータの平均 <x>
シミュレーション物理5 運動方程式の方法: サブルーチンの使い方.
心理学情報処理法Ⅰ やってみよう:Excelを使ってみよう.
放射線の計算や測定における統計誤差 「平均の誤差」とその応用(1H) 2項分布、ポアソン分布、ガウス分布(1H) 最小二乗法(1H)
回帰分析.
数理統計学  第8回 第2章のエクササイズ 西山.
情報教育論 第9回 仮定文の仕組み 政策・メディア研究科 岡田 健.
応用統計学の内容 推測統計学(inferential statistics)   連続型の確率分布   標本分布   統計推定   統計的検定.
問題 1 キーボードから入力した数の合計を計算するプログラムを 作成せよ。最初に、何個の数を入力するかその数を入力 するようにする。
寺尾 敦 青山学院大学社会情報学部 エクセルでの正規分布の グラフの描き方 寺尾 敦 青山学院大学社会情報学部
湘南工科大学 2013年12月10日 プログラミング基礎1 湘南工科大学情報工学科 准教授 小林 学.
ネットワーク性能評価.
寺尾 敦 青山学院大学社会情報学部 エクセルでの正規分布の グラフの描き方 寺尾 敦 青山学院大学社会情報学部
ー 第1日目 ー 確率過程について 抵抗の熱雑音の測定実験
概要 Boxed Economy Simulation Platform(BESP)とその基本構造 BESPの設計・実装におけるポイント!
プロジェクト演習Ⅱ インタラクティブゲーム制作 イントロダクション2
情報工学Ⅱ (第9回) 月曜4限 担当:北川 晃.
プログラミング演習(2組) 第8回
応用統計学の内容 推測統計学(inferential statistics)   連続型の確率分布   標本分布   統計推定   統計的検定.
混合ガウスモデルによる回帰分析および 逆解析 Gaussian Mixture Regression GMR
モデルの逆解析 明治大学 理工学部 応用化学科 データ化学工学研究室 金子 弘昌.
寺尾 敦 青山学院大学社会情報学部 エクセルでの正規分布の グラフの描き方 寺尾 敦 青山学院大学社会情報学部
東北大学 大学院情報科学研究科 応用情報科学専攻 田中 和之(Kazuyuki Tanaka)
湘南工科大学 2013年10月22日 プログラミング基礎1 湘南工科大学情報工学科 准教授 小林 学.
デジタル画像とC言語.
情報処理 タイマの基礎 R8C タイマの基礎.
プロジェクト演習Ⅱ インタラクティブゲーム制作
岩村雅一 知能情報工学演習I 第9回(後半第3回) 岩村雅一
情報理論2 第3回 小林 学 湘南工科大学 2011年10月25日 〒 神奈川県藤沢市辻堂西海岸1-1-25
情報処理Ⅱ 第2回:2003年10月14日(火).
or-4. モンテカルロシミュレーション (オペレーションズリサーチを Excel で実習するシリーズ)
東京工科大学 コンピュータサイエンス学部 亀田弘之
生物統計学・第3回 全体を眺める(2) クラスタリング、ヒートマップ
疑似乱数, モンテカルロ法によるシミュレーション
ex-8. 平均と標準偏差 (Excel 実習シリーズ)
「データ学習アルゴリズム」 第3章 複雑な学習モデル 報告者 佐々木 稔 2003年6月25日 3.1 関数近似モデル
ベイズ最適化 Bayesian Optimization BO
地域情報学 C言語プログラミング 第2回 変数・配列、型変換、入力 2017年10月20日
シミュレーション論Ⅰ 第7回 シミュレーションの構築と実施.
東京工科大学 コンピュータサイエンス学部 亀田弘之
最尤推定・最尤法 明治大学 理工学部 応用化学科 データ化学工学研究室 金子 弘昌.
シミュレーション物理 スピングラスなどについて.
プログラミング入門2 第13回、14回 総合演習 情報工学科 篠埜 功.
情報工学Ⅱ (第9回) 月曜4限 担当:北川 晃.
数理統計学 西 山.
データ解析 静岡大学工学部 安藤和敏
小標本に関する平均の推定と検定 標本が小さい場合,標本分散から母分散を推定するときの不確実さを加味したt分布を用いて,推定や検定を行う
or-4. モンテカルロシミュレーション (オペレーションズリサーチを Excel で実習するシリーズ)
or-6. 待ち行列シミュレーション (オペレーションズリサーチを Excel で実習するシリーズ)
プログラミング論 相関
ex-8. 平均と標準偏差 (Excel を演習で学ぶシリーズ)
or-4. モンテカルロシミュレーション (オペレーションズリサーチを Excel で実習するシリーズ)
情報工学Ⅱ (第8回) 月曜4限 担当:北川 晃.
cp-15. 疑似乱数とシミュレーション (C プログラミング演習,Visual Studio 2019 対応)
統計解析 第11回.
情報処理Ⅱ 第2回 2004年10月12日(火).
「速さ」で学んだことを用いて 課題を解決しよう
確率的フィルタリングを用いた アンサンブル学習の統計力学 三好 誠司 岡田 真人 神 戸 高 専 東 大, 理 研
プログラミング入門2 第5回 配列 変数宣言、初期化について
第1章 文字の表示と計算 printfと演算子をやります.
混合ガウスモデル Gaussian Mixture Model GMM
Presentation transcript:

シミュレーション物理7 乱数

今回の授業の目的 乱数の発生のさせ方を覚える なぜ乱数が重要か? 乱数は数値実験の基本 自然界の複雑な現象を乱数を使ってまねする 温度による雑音,サンプル平均など

乱数の出し方 基本的には大きな整数に大きな整数をかける->次の値が予測がつかなくなることを利用 4byte整数の範囲:232/2= 2147483648なので大きな整数に大きな整数をかけると  -2147483648~2147483647 の間に一様に分布することを利用する。これをプログラムで書くとこうなる。 iran=65539*iran+125654 !大きな整数iranに適当な大きな整数をかける,足す if(iran<0)iran=iran+2147483647+1 !2147483647を越えると負になるので補正 x=iran/2147483648.0 !結果が0~1の間になるように規格化

この乱数には問題がある 一見ランダムに見えて,相関がある。 周期も短い よってよくできたサブルーチンを利用する。

乱数の分布 乱数の分布としてよく使われる二つ 一様分布 0~1の間に同じ確率で分布 正規分布 釣り鐘状のガウス分布,平均は0,分散は1に規格化されているものを用意し,用途に応じて使う 実は上手いやり方があり,一様分布乱数2個から正規分布乱数2個を作ることが出来る(Box-Muller法)

確率変数の変数変換 変数2個では?

Box-Muller法

本日使うプログラム KindNumbers.f90 !double precisionを定義 random.f90 !一様乱数のサブルーチン random_uniform.f90 ![0,1]の一様乱数をサブルーチンを使って発生させるメインプログラム random_gauss.f90 ! 上の正規分布(平均0,分散1)バージョン 講義ノートのページからrandom.tarをPCにダウンロードして、dahlmanにアップロード tar xvf random.tar で解凍する

実行のさせ方 まずモジュールをコンパイル 実行ファイルを作る 次に実行 f90 –c KindNumbers.f90 f90 –c random.f90 実行ファイルを作る f90 –o randomuniform random_uniform.f90 KindNumbers.o random.o f90 –o randomgauss random_gauss.f90 KindNumbers.o random.o 次に実行 randomuniform randomgauss

課題  乱数を1000個発生させ,このヒストグラムを作る。一様乱数,正規分布乱数両方とも。(乱数を画面でなくファイルに書きだし,それをエクセルで解析すればよい。分析ツールを使う。0.2間隔の列を作りそれを範囲として選択する。(列ごとは選ばない。出力は他のシートにする。)) 乱数を使って円周率を求めるプログラムを作る。やり方は黒板で説明。これは提出課題(12月末までに提出)。