範例：亂數種子產生

函數 > 統計 > 機率分佈 > 範例：亂數種子產生

範例：亂數種子產生

使用 Seed 函數控制產生新亂數集的流程。

1. 在簡短的程式中使用函數 Seed 與 rnorm，透過重設種子值，為三個疊代產生相同的常態分佈亂數集。

<region id="ID0EIHTCB" top="268.80000000000007" left="38.400000000000006" height="117.15199999999999" width="372.98" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="84">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id labels="VARIABLE" xml:space="preserve">SameNormal</id>
        <boundVars>
          <id labels="VARIABLE" xml:space="preserve">p</id>
          <id labels="VARIABLE" xml:space="preserve">μ</id>
          <id labels="VARIABLE" xml:space="preserve">σ</id>
        </boundVars>
      </function>
      <program>
        <for>
          <id labels="VARIABLE" xml:space="preserve">i</id>
          <range>
            <real>0</real>
            <real>2</real>
          </range>
          <program>
            <apply>
              <id labels="FUNCTION" xml:space="preserve">種子</id>
              <real>1</real>
            </apply>
            <localDefine>
              <id labels="VARIABLE" xml:space="preserve">nums</id>
              <apply>
                <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">rnorm</id>
                <sequence>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">p</id>
                  <id labels="VARIABLE" xml:space="preserve">μ</id>
                  <id labels="VARIABLE" xml:space="preserve">σ</id>
                </sequence>
              </apply>
            </localDefine>
            <localDefine>
              <apply>
                <matcol />
                <id labels="VARIABLE" xml:space="preserve">M</id>
                <id labels="VARIABLE" xml:space="preserve">i</id>
              </apply>
              <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">nums</id>
            </localDefine>
          </program>
        </for>
        <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
      </program>
    </define>
    <resultFormat>
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

• 因為種子設為 1，所以每次都會產生相同的亂數集。

• SameNormal 會傳回有相同三欄的矩陣，每欄各代表一個疊代。

2. 使用函數 hist 與 max，繪製含有 n_bins 個柱狀之第二個亂數集的柱狀圖。新增水平標記以顯示傳回集合的最大值。

<region id="ID0EVOTCB" top="691.2" left="38.400000000000006" height="45.2" width="228.19666666666669" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="94">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <indexer />
        <id labels="VARIABLE" xml:space="preserve">S1</id>
        <id labels="*" xml:space="preserve">q</id>
      </apply>
      <apply>
        <indexer />
        <parens>
          <apply>
            <matcol />
            <apply>
              <id labels="VARIABLE" xml:space="preserve">SameNormal</id>
              <sequence>
                <id labels="VARIABLE" xml:space="preserve">p</id>
                <id labels="VARIABLE" xml:space="preserve">μ</id>
                <id labels="VARIABLE" xml:space="preserve">σ</id>
              </sequence>
            </apply>
            <real>1</real>
          </apply>
        </parens>
        <id xml:space="preserve">q</id>
      </apply>
    </define>
    <resultFormat>
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

<region id="ID0E6PTCB" top="921.60000000000014" left="38.400000000000006" height="235" width="575.781818181818" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="100">
          <traceStyle color="#FF9ABC66" symbol="none" line-weight="1" line-style="Solid">欄</traceStyle>
        </trace>
      </traces>
      <graph-size width="456.581818181818" height="177.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
          <axisLine position="origin" positionticmark="0" legendWidth="92.1466666666667" />
          <axisGrid>
            <gridFrequency>5</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="107">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">range</id>
              </math>
              <math resultRef="108">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="102">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="104">
              <real xmlns="http://schemas.mathsoft.com/math50">5</real>
            </secondTickValue>
            <endValue resultRef="106">
              <real xmlns="http://schemas.mathsoft.com/math50">20</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="150">
          <axisLine position="origin" positionticmark="0" legendWidth="87.9033333333333" />
          <axisGrid>
            <gridFrequency>7</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers>
            <marker resultRef="109" display="true" color="#FFED1D2F" line-style="Dash" value="129">
              <apply xmlns="http://schemas.mathsoft.com/math50">
                <id xml:space="preserve">max</id>
                <id xml:space="preserve">same</id>
              </apply>
            </marker>
          </markers>
          <numberFormat>
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="117">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">same</id>
              </math>
              <math resultRef="118">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="112">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="114">
              <real xmlns="http://schemas.mathsoft.com/math50">25</real>
            </secondTickValue>
            <endValue resultRef="116">
              <real xmlns="http://schemas.mathsoft.com/math50">150</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

3. 修改程式以透過增量每個疊代的種子值，產生三個不同的常態分佈亂數集。

<region id="ID0EPRTCB" top="1267.2000000000003" left="38.400000000000006" height="117.15199999999999" width="365.72" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="130">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id labels="VARIABLE" xml:space="preserve">DiffNormal</id>
        <boundVars>
          <id labels="VARIABLE" xml:space="preserve">p</id>
          <id labels="VARIABLE" xml:space="preserve">μ</id>
          <id labels="VARIABLE" xml:space="preserve">σ</id>
        </boundVars>
      </function>
      <program>
        <for>
          <id labels="VARIABLE" xml:space="preserve">i</id>
          <range>
            <real>0</real>
            <real>2</real>
          </range>
          <program>
            <apply>
              <id labels="FUNCTION" xml:space="preserve">種子</id>
              <apply>
                <plus />
                <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">i</id>
                <real>1</real>
              </apply>
            </apply>
            <localDefine>
              <id labels="VARIABLE" xml:space="preserve">nums</id>
              <apply>
                <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">rnorm</id>
                <sequence>
                  <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">p</id>
                  <id labels="VARIABLE" xml:space="preserve">μ</id>
                  <id labels="VARIABLE" xml:space="preserve">σ</id>
                </sequence>
              </apply>
            </localDefine>
            <localDefine>
              <apply>
                <matcol />
                <id labels="VARIABLE" xml:space="preserve">M</id>
                <id labels="VARIABLE" xml:space="preserve">i</id>
              </apply>
              <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">nums</id>
            </localDefine>
          </program>
        </for>
        <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
      </program>
    </define>
    <resultFormat>
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

• 既然種子值會隨著每次疊代而變更，因此每個亂數集都是唯一的。

• 由於兩個程式的疊代 0 都是以種子 1 起始，所以這兩個傳回矩陣的欄 0 會包含相同的集合。

4. 從傳回的矩陣萃取第二個數據集：

<region id="ID0E4VTCB" top="1593.6000000000001" left="38.400000000000006" height="45.2" width="220.93666666666667" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="132">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <indexer />
        <id labels="VARIABLE" xml:space="preserve">S2</id>
        <id labels="*" xml:space="preserve">q</id>
      </apply>
      <apply>
        <indexer />
        <parens>
          <apply>
            <matcol />
            <apply>
              <id labels="VARIABLE" xml:space="preserve">DiffNormal</id>
              <sequence>
                <id labels="VARIABLE" xml:space="preserve">p</id>
                <id labels="VARIABLE" xml:space="preserve">μ</id>
                <id labels="VARIABLE" xml:space="preserve">σ</id>
              </sequence>
            </apply>
            <real>1</real>
          </apply>
        </parens>
        <id xml:space="preserve">q</id>
      </apply>
    </define>
    <resultFormat>
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

5. 使用函數 hist 產生含有 n_bins 個柱狀之亂數集的柱狀圖，然後計算其最大值。

6. 繪製柱狀圖，並使用水平標記顯示其最大值：

<region id="ID0EPZTCB" top="1785.6000000000004" left="38.400000000000006" height="235" width="575.781818181818" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="138">
          <traceStyle color="#FFF47D44" symbol="none" line-weight="1" line-style="Dash">欄</traceStyle>
        </trace>
      </traces>
      <graph-size width="456.581818181818" height="177.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
          <axisLine position="origin" positionticmark="0" legendWidth="92.1466666666667" />
          <axisGrid>
            <gridFrequency>5</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="145">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">range</id>
              </math>
              <math resultRef="146">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="140">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="142">
              <real xmlns="http://schemas.mathsoft.com/math50">5</real>
            </secondTickValue>
            <endValue resultRef="144">
              <real xmlns="http://schemas.mathsoft.com/math50">20</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="175">
          <axisLine position="origin" positionticmark="0" legendWidth="80.32" />
          <axisGrid>
            <gridFrequency>8</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers>
            <marker resultRef="147" display="true" color="#FF662D91" line-style="Dash" value="155">
              <apply xmlns="http://schemas.mathsoft.com/math50">
                <id xml:space="preserve">max</id>
                <id xml:space="preserve">diff</id>
              </apply>
            </marker>
          </markers>
          <numberFormat>
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="155">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">diff</id>
              </math>
              <math resultRef="156">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="150">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="152">
              <real xmlns="http://schemas.mathsoft.com/math50">25</real>
            </secondTickValue>
            <endValue resultRef="154">
              <real xmlns="http://schemas.mathsoft.com/math50">175</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

相關主題

均等分佈

常態分佈

指數分佈

這是否有幫助?