シミュレーション論Ⅰ 第１０回 様々なシミュレーション：金利とローン返済. 第９回のレポート ポアソン分布に従う乱数列（乱数表）から乱数を記入 する。 乱数値をその日の客数として、仕入部数が 8 ～ 12 のと きの利益を記入する。 10 日分のシミュレーションをおこない、総売上と最も 利益の高かった仕入部数を調べる。

Presentation on theme: "シミュレーション論Ⅰ 第１０回 様々なシミュレーション：金利とローン返済. 第９回のレポート ポアソン分布に従う乱数列（乱数表）から乱数を記入 する。 乱数値をその日の客数として、仕入部数が 8 ～ 12 のと きの利益を記入する。 10 日分のシミュレーションをおこない、総売上と最も 利益の高かった仕入部数を調べる。"— Presentation transcript:

1 シミュレーション論Ⅰ 第１０回 様々なシミュレーション：金利とローン返済

2 第９回のレポート ポアソン分布に従う乱数列（乱数表）から乱数を記入 する。 乱数値をその日の客数として、仕入部数が 8 ～ 12 のと きの利益を記入する。 10 日分のシミュレーションをおこない、総売上と最も 利益の高かった仕入部数を調べる。 仕入れ価格 c = 80 販売価格 a = 120 1 日の客数 x （乱数表から決定） 仕入量 y （ 8 部、 9 部、 10 部、 11 部、 12 部） 1 日の利益

3 第９回のレポート 回答例（ 9 部・ 11 部も入れた場 合） ポアソン乱数表を用いて日々の客数を記入し、売り上げを計算する 上の例では「９部仕入れ」の場合に総利益がもっとも高くなった 日数乱数 利益（ 8 部仕 入） 利益（ 9 部仕 入） 利益（ 10 部仕 入） 利益（ 11 部仕 入） 利益（ 12 部仕 入） 18320240160800 29320360280200120 36800-80-160-240 413320360400440480 513320360400440480 613320360400440480 716320360400440480 813320360400440480 96800-80-160-240 1011320360400440360 総利益 27202760268026002400

4 復習： EOQ 公式の導出

5 今回の内容 数式によって厳密に定義され、かつ解析的に解ける問 題のシミュレーションについて考える 金利計算とローン返済のモデルを通じて、確定的なシ ミュレーションの意義を知る Excel を用いて前回の新聞売り子問題のシミュレー ションをおこなう

6 雑学：曽呂利新左衛門の褒美 昔、羽柴秀吉の家臣（御伽衆）に曽呂利新左衛門と いう男がいました。ある日、将棋に負けた秀吉が褒 美の希望を聞いたところ、 「今日は米１粒、明日は２粒、翌日はその倍の 4 粒、 その翌日は８粒というように３０日間いただきた い」 と答えたということです。さて、３０日後に秀吉は 何粒の米を与えることになったでしょうか？

7 雑学：曽呂利新左衛門の褒美（２） １日目・１粒が 10 日目には 512 粒、 20 日目には 524,288 粒となり約 15kg 、 22 日目には 2,097,152 粒で 米俵 60kg ・ 1 俵分となる。 30 日目には、何と 536,870,912 粒・米俵 256 俵 (100 石 の殿様 ) にもなる。 解析的に解くにはどうすればいいか考えてみよう。

8 確定的なモデルのシミュレーション 数式によって厳密に定義され、かつ解析的に解ける問 題のシミュレーションは無意味だろうか？ 先ほどの例のように、モデル、数式が分かっていても 「数式だけでは分かりにくい」、「単純な予想を超え る」、「様々な場合を比較したい」場合など、シミュ レーションをおこなうことにより理解を助けることが できる。 身近なところでは、金利の計算（利子・利息）やロー ン返済額、年金額のシミュレーションなどが見受けら れる。

9 単純な数値計算としてのシミュレーション 例：ローン返済のシミュレーション 銀行などからお金を借りるとして、どのような返済方法がよい か？ 様々な場合をあらかじめ試したり、分かりやすく相手に示した りできる。

10 金利と残高 金利の計算方法には大きく分けて「単利」と「複利」 がある。 単利：最初に預けられた（借りた）元金に対してのみ 利息を 計算する方法 複利：一定期間の利息を元金に加え、その合計を新た な元 金として利息を計算する方法

11 単利 単利の元利合計： 元利合計＝元本 × （１＋年利率 × 預入年数） 例） 1 万円を年利率 1 ％の単利で預金したとすると 1 年後： 10,000×(1+0.01×1) = 10,100 2 年後： 10,000×(1+0.01×2) = 10,200 3 年後： 10,000×(1+0.01×3) = 10,300 ・・・ ※利息は元本の 1 万円についてのみ計算される ＝毎年同じ利息がつく ※単利の元利合計は等差数列になる

12 複利 複利の元利合計： 元利合計＝元本 × （１＋利率） 預入期間 複利の利率と預入期間： １年複利 → 利率は年利率、預入期間は１年を１期間とす る。 半年複利 → 利率は（年利率 ÷ ２）、預入期間は半年を１期 間とする。（１年は２期間） １ヶ月複利 → 利率は（年利率 ÷ １２）、預入期間は１ヶ月 を１期間とする。（１年は１２期間）

13 複利（２） 例） 1 万円を年利率 2 ％の半年複利で預金したとすると 半年あたりの利率 ＝ 2÷2 ＝ 1 (%) 半年後： 10,000× （ 1+0.01) 1 ＝ 10,100 1 年後： 10,000× （ 1+0.01) 2 ＝ 10,201 1 年半後： 10,000× （ 1+0.01) 3 ≒ 10,303 ・・・ ※利息は一定期間ごとに（元本＋利息）を新たな元本と して 計算される＝利息が期間ごとに増えていく ※複利の元利合計は等比数列になる

14 例：単利と複利の比較 年利率 5 ％の単利と１年複利で 10 万円を銀行に預けた場 合、５年後までの毎年の利息と元利合計を計算してみ よう。 単利、複利それぞれの元利合計を X T, X F とし、年数を n とすると

15 例：単利と複利の比較 Excel で計算して１０年後までの結果をグラフにすると 年数単利の利息単利の元利合計複利の利息複利の元利合計 15000105000\5,000\105,000 25000110000\5,250\110,250 35000115000\5,513\115,763 45000120000\5,788\121,551 55000125000\6,078\127,628

16 例：複利の比較 年利率 10 ％の複利で 10 万円を銀行に預ける。１年複利 と半年複利の場合について５年後まで計算してみよう。 １年複利、半年複利それぞれの元利合計を X 1, X 0.5 とし、 年数を n 、半年の期間を m とすると

17 例：複利の比較 年数年数 １年複利の 元利合計 半年の期間 数 半年複利の 元利合計 1 \105,000 11100002 ＝１年 \110,250 3 \115,763 21210004 ＝２年 \121,551 5 \127,628 31331006 ＝３年 \134,010 7 \140,710 41464108 ＝４年 \147,746 9 \155,133 516105110 ＝５年 \162,889

18 ローン返済：元利均等返済 元利均等返済方式： 毎回の返済額（元金，利息の合計）を均等にした返 済方式 。 ローンで最も普及した返済方式で、裁判所の調停で は一般にこの返済方式が用いられている。

19 元利均等返済のシミュレーション １０万円を年利１２％の１ヶ月複利（つまり月 1 ％の複 利）で借り入れ、元利均等返済をする。 ６ヶ月で返す場合と１２ヶ月で返す場合のそれぞれに ついて、毎回の返済金額を計算せよ。 Excel で以下のように入力すると、６回（６ヶ月）で返 済する場合の毎回の返済額が分かる（実際の額は小数 点以下を切り捨てる）。 返済回数を変えて１２回で返済する場合についても試 してみよう。

20 参考：小数点以下の切捨て Excel で小数点以下を切り捨てて整数にするには、 =INT() 関数を用いる 先程の「毎回の返済額」の数式を =INT() のカッコ内に 入れてやるとよい

21 新聞売り子問題のシミュレーション ポアソン分布に従う乱数値をその日の客数、仕入部数 を１～ 20 部として新聞売り子問題のシミュレーション を行う。 10 日分のシミュレーションを繰り返しおこない、最も 平均利益の高かった仕入部数を調べる。 仕入れ価格 c = 80 販売価格 a = 120 1 日の客数 x （乱数表から決定） 仕入量 y （ 8 部、 10 部、 12 部） 1 日の利益

22 ポアソン乱数の近似 Excel でポアソン分布に従う乱数（ポアソン乱数）を生 成するのは結構面倒（ VBA マクロやポアソン分布の表 から作成する方法がある） 平均値 λ が比較的大きい場合、正規分布によってポア ソン分布の近似ができる（平均値 λ 、標準偏差を とする）

23 乱数の生成と客数の決定 以下のような表を作成し、正規乱数でポアソン乱数を 近似する（平均 λ 、標準偏差は λ の平方根）。 まれに負の値が出るので、 MAX 関数と INT 関数を使って 0 以上の整数値に直す。 できたら下へコピー

24 仕入れ部数の設定 仕入れ部数を 1 ～ 20 部としてシミュレーションするため の枠を作成する D1 ～ W1 まで、 1 ～ 20 の数値を入れる 12 行目に利益の合計を計算する欄を作成しておく

25 利益の計算 IF 関数を使って、その日の客数と仕入れ部数から利益を 算出する 客数＜仕入れ部数・・・客数 ×120 －仕入れ部数 ×80 客数≧仕入れ部数・・・仕入れ部数 × （ 120 － 80 ） 入力できたら横・縦へコピーして 10 日分のシミュレー ションを完成させる

26 総利益の計算 SUM 関数を使って、仕入れ部数ごとの総利益を計算す る 入力できたら右へ（ W 列まで）コピーしておく

27 集計部分の作成（１） 14 ～ 15 行に繰り返し回数・総利益の合計・総利益の平 均を記入する欄を作成する

28 集計部分の作成（２） 繰り返し回数、総利益の合計、総利益の平均を計算す る 循環参照のエラーが出るが、キャンセルを押すこと 入力できたら、総利益の合計と平均について右へコ ピーしておく

29 グラフの作成 循環参照を許可する前に、グラフを作成しておく 総利益の平均値を 1 ～ 20 部まで選択し、縦棒グラフを作成 する まだ数値が入っていないので棒が出ないが構わない

30 反復計算を許可する 反復計算を許可する前にファイルに名前をつけて保存 しておく（何か失敗したらそこからやり直せる） 「ツール」 → 「オプション」 → 「計算方法」で「反復計 算」にチェックを入れる。 計算方法は「手動」にしておくと F9 キーだけで計算で きる。 最大反復回数は 1 とする。

31 完成 F9 キーを押して繰り返し計算をしてみよう 仕入れ部数によって総利益（の平均）がどのように異 なるかが分かる

32 第１０回のレポート （１） 先の元利均等返済の問題について、返済回数を 6 回、 12 回としたときの毎月の返済額を記入せよ （２） 作成したシミュレーションを 1000 回繰り返し、仕 入れ部数 が 5 ～ 15 部のときの平均利益をそれぞれ記入し、利 益が 最大となる仕入れ部数を調べよ ノート PC のない方は以下の課題 別課題の課題１、課題２をやって解答を記入