龍谷大学>理工学部>数理情報学科>樋口>担当科目>2011年>計算科学☆演習II>01回め
目次 前回 次回 略解
計算科学☆演習 II
樋口さぶろお
*1 配布: 2011-04-15 Fri更新: Time-stamp: ”2012-04-09 Mon 14:08 JST hig”はじめに
講義ののり
必ず計算科学演習
IIと同時に履修してね
.■成績計算
科目の成績
100ピーナッツは
• 10
ピーナッツ
:毎回授業での
quiz.授業の最後にやります
.• 30
ピーナッツ
:プチテスト いまのところ
2011-05-27を予定
• 50
ピーナッツ
:ファイナルトライアル
• 10
ピーナッツ
:授業時間外の予習復習 シラバスから変更するかも 平常点
,予習復習的 な何か
.今日のところは
eラーニングシステムでの予習復習はありません
.■資料
授業で配布
.授業後に欲しい人は
http://hig3.netから各自ダウンロード
. 1-503前のレターボックスに残ってることも
.■欠席届
去年とポリシー違います
.専用用紙に事情を説明する書類を貼って
,授業前後各
5分に提出
(欠席の前でも後でも可
.ただしファイナルトライアルが締切
)演習ののり
かならず計算科学
IIと同時に履修してね
.実習室に行ったら
, http://hig3.net→計算科学☆演習
IIへ
.別紙も参照
.*1 Copyright c⃝2011Saburo HIGUCHI. All rights reserved.
, http://hig3.net(講義のページもここからたどれます), へや:1
号館5階502.
担当者ののり
•
なまえ
:樋口さぶろお
•
へや
: 1-502•
オフィスアワー
:木昼
,金
4.その他
,質問歓迎な時間
:月昼
,火昼
. (水昼
,金昼は授業 や会議の前後
)お弁当持参可
.お湯あげます
.• Web
ページ
:携帯
/PC http://hig3.net今日の演習の指示や
,プチテストのスケ ジュールもここから
.1 ランダムウォークと乱数
今日の目標
•
ランダムウォークとふつうの数列の違いが説明できる
• C
言語で疑似乱数を生成する方法を説明できる
1.1 quiz: 疑似乱数の使いかた
サ ン プ ル プ ロ グ ラ ム を 参 考 に
,関 数
double uniform()と ラ イ ブ ラ リ 関 数
void srand(int)は与えられたとして
, scanfでシードを取得
,設定し
, −1,0,+1をそれぞれ
1/3の確率で出力することを
100回繰り返すプログラムを書こう
.20 40 60 80 100 t
-6 -4 -2 2 4 6 x
2
疑似乱数の正しい / 間違ったプログラム
ソースコード1 乱数
1 /*
2 r a n d 1 . c - - [0 ,1)疑 似 乱 数 の サ ン プ ル プ ロ グ ラ ム 3 Time - s t a m p : " 20 11 -0 4 -1 5 Fri 1 2 : 4 8 JST hig "
4 */
5 # d e f i n e _ C R T _ S E C U R E _ N O _ W A R N I N G S // お ま じ な い 6 # i n c l u d e < s t d i o . h >
7 # i n c l u d e < s t d l i b . h > /* s r a n d () , r a n d () を 使 う の に 必 要 */
8
9 d o u b l e u n i f o r m (); /* 関 数 プ ロ ト タ イ プ 宣 言 */
10
11 int m a i n (){
12 int s e e d ; /* 疑 似 乱 数 の シ ー ド */
13 d o u b l e s ;
14 int t ; /* カ ウ ン タ */
15 int t m a x = 1 0 0 ; /* 疑 似 乱 数 を 得 る 回 数 */
16
17 s c a n f ("% d " ,& s e e d );
18
19 s r a n d ( s e e d ); /* シ ー ド の 設 定 */
20 for ( t =0; t < t m a x ; t + + ) {
21 /* s r a n d ( s e e d ); */ /*こ こ に 置 く と ど う な る の? */
22 s = u n i f o r m ();
23 p r i n t f ("% f \ n " , s );
24 }
25 r e t u r n 0;
26 } 27
28 /** [0 ,1) 疑 似 乱 数 を 返 す */
29 d o u b l e u n i f o r m (){
30 r e t u r n r a n d ( ) / ( R A N D _ M A X + 1 . 0 ) ; 31 }
履修要項の該当ページのコピー
(PDF版では省略
)目次 前回 次回 略解
3
