例: カスタム分布関数を使用したモンテカルロシミュレーション

関数 > 実験計画 > モンテカルロシミュレーション > 例: カスタム分布関数を使用したモンテカルロシミュレーション

この例では、カスタム分布関数を montecarlo 関数とともに使用する方法を示します。カスタム分布関数を使用すると、一様分布関数を使用する場合と比較して計算速度が向上します。

1. 一様分布に基づくカスタム分布関数を作成します。

<region id="ID0EYKEX" actualWidth="349.11666666666673" actualHeight="113.60000000000001" top="96.000000000000014" left="19.200000000000003" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="4">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id labels="VARIABLE" xml:space="preserve">CustomDist</id>
        <boundVars>
          <id labels="VARIABLE" xml:space="preserve">n</id>
        </boundVars>
      </function>
      <program>
        <if>
          <test>
            <apply>
              <equal />
              <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">n</id>
              <real>1</real>
            </apply>
          </test>
          <then>
            <program>
              <apply>
                <indexer />
                <apply>
                  <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">runif</id>
                  <sequence>
                    <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">n</id>
                    <real>-0.01</real>
                    <real>0.01</real>
                  </sequence>
                </apply>
                <id labels="SYSTEM" xml:space="preserve" label-is-contextual="true">ORIGIN</id>
              </apply>
            </program>
          </then>
          <else>
            <program>
              <apply>
                <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">runif</id>
                <sequence>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">n</id>
                  <real>-0.01</real>
                  <real>0.01</real>
                </sequence>
              </apply>
            </program>
          </else>
        </if>
      </program>
    </define>
  </math>
</region>

2. シミュレートする関数を定義します。

3. 生成するサンプルの数を定義します。

4. montecarlo を呼び出してサンプルを生成します。

<region id="ID0ENPEX" actualWidth="704.6733333333334" actualHeight="261.62" top="556.80000000000007" left="19.200000000000003" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="7">
    <eval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">montecarlo</id>
        <sequence>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">F</id>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">N</id>
          <matrix rows="1" cols="3">
            <str xml:space="preserve">入力</str>
            <real>0</real>
            <real>1</real>
          </matrix>
          <matrix rows="1" cols="1">
            <id labels="VARIABLE" xml:space="preserve">CustomDist</id>
          </matrix>
          <matrix rows="1" cols="2">
            <real>-0.01</real>
            <real>0.01</real>
          </matrix>
        </sequence>
      </apply>
      <unitOverride>
        <placeholder />
      </unitOverride>
    </eval>
  </math>
</region>