範例：產生亂數

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

範例：產生亂數

建立在區間上各自均等分佈、常態分佈或指數分佈的亂數向量。

均勻分佈

1. 呼叫函數 runif 產生在高低之間均等分佈的 n_set 亂數。

2. 使用函數 hist 產生含有 n_bins 個柱狀的均等集合柱狀圖，然後計算此集合的均數。

3. 繪製柱狀圖並顯示數據集的預期均數：

<region id="ID0E56BCB" top="480.00000000000006" left="38.400000000000006" height="216" width="576" 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="202">
          <traceStyle color="#FF31ADC2" symbol="none" line-weight="1" line-style="Solid">欄</traceStyle>
        </trace>
        <trace resultRef="203">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Dash">線條</traceStyle>
        </trace>
      </traces>
      <graph-size width="437.6" height="158.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
          <axisLine position="origin" positionticmark="0" legendWidth="92.1466666666667" />
          <axisGrid>
            <gridFrequency>11</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="208">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">range</id>
              </math>
              <math resultRef="209">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="205">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue resultRef="207">
              <real xmlns="http://schemas.mathsoft.com/math50">20</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="75">
          <axisLine position="origin" positionticmark="0" legendWidth="111.033333333333" />
          <axisGrid>
            <gridFrequency>6</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="216">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">uniform</id>
              </math>
              <math resultRef="217">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="218">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">mean</id>
              </math>
              <math resultRef="219">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="211">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="213">
              <real xmlns="http://schemas.mathsoft.com/math50">15</real>
            </secondTickValue>
            <endValue resultRef="215">
              <real xmlns="http://schemas.mathsoft.com/math50">67</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

若要取得新的亂數集，請在「計算」標籤上，按一下「控制項」群組的「計算」。每次重新計算時，都使用新的種子值產生亂數。

常態分佈

1. 呼叫函數 rnorm 產生含有均數為 μ 且標準差為 σ 之常態分佈的 n_set 亂數。

2. 使用函數 floor、ceil、min 及 max 計算數據集值域的最小值與最大值，然後計算每個柱狀的寬度：

3. 設定繪圖的水平縮放：

4. 呼叫函數 dnorm 計算數據集的預期常態分佈 F，然後繪製含有 n_bins 個柱狀的常態集合柱狀圖。

<region id="ID0EBNCCB" top="1344.0000000000002" left="38.400000000000006" height="36.800000000000004" width="227.07666666666668" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="247">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <indexer />
        <id labels="*" xml:space="preserve">F</id>
        <id labels="*" xml:space="preserve">n</id>
      </apply>
      <apply>
        <mult />
        <apply>
          <mult />
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">n_set</id>
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">w</id>
        </apply>
        <apply>
          <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">dnorm</id>
          <sequence>
            <apply>
              <indexer />
              <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">int</id>
              <id xml:space="preserve">n</id>
            </apply>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">μ</id>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">σ</id>
          </sequence>
        </apply>
      </apply>
    </define>
    <resultFormat>
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

<region id="ID0E4NCCB" top="1401.6000000000001" left="38.400000000000006" height="222.00000000000003" width="575" 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="251">
          <traceStyle color="#FF662D91" symbol="none" line-weight="1" line-style="Solid">欄</traceStyle>
        </trace>
        <trace resultRef="252">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Dash">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="436.6" height="164.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="-7" end="7">
          <axisLine position="origin" positionticmark="0" legendWidth="72.65" />
          <axisGrid>
            <gridFrequency>21</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="259">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">int</id>
              </math>
              <math resultRef="260">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="254">
              <real xmlns="http://schemas.mathsoft.com/math50">-7.0</real>
            </startValue>
            <secondTickValue resultRef="256">
              <real xmlns="http://schemas.mathsoft.com/math50">-6.3</real>
            </secondTickValue>
            <endValue resultRef="258">
              <real xmlns="http://schemas.mathsoft.com/math50">7.0</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="150">
          <axisLine position="origin" positionticmark="10" legendWidth="102.85" />
          <axisGrid>
            <gridFrequency>6</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="267">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">normal</id>
              </math>
              <math resultRef="268">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="269">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">F</id>
              </math>
              <math resultRef="270">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="262">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="264">
              <real xmlns="http://schemas.mathsoft.com/math50">30</real>
            </secondTickValue>
            <endValue resultRef="266">
              <real xmlns="http://schemas.mathsoft.com/math50">150</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

指數分佈

1. 呼叫函數 rexp 產生一組按機率 r 以指數方式分佈的亂數。

2. 呼叫函數 dexp 計算數據集的預期指數分佈 F，然後繪製含有 n_bins 個柱狀的指數集合柱狀圖。

<region id="ID0E5VCCB" top="2131.2000000000003" left="38.400000000000006" height="235" width="576" 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="301">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">欄</traceStyle>
        </trace>
        <trace resultRef="302">
          <traceStyle color="#FF2E3192" symbol="none" line-weight="1" line-style="Dash">線條</traceStyle>
        </trace>
      </traces>
      <graph-size width="418.4" height="177.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
          <axisLine position="origin" positionticmark="0" legendWidth="72.65" />
          <axisGrid>
            <gridFrequency>11</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="307">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">int</id>
              </math>
              <math resultRef="308">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="304">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue resultRef="306">
              <real xmlns="http://schemas.mathsoft.com/math50">20</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="400">
          <axisLine position="origin" positionticmark="0" legendWidth="133.973333333333" />
          <axisGrid>
            <gridFrequency>9</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="315">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">exponential</id>
              </math>
              <math resultRef="316">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="317">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">F</id>
              </math>
              <math resultRef="318">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="310">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="312">
              <real xmlns="http://schemas.mathsoft.com/math50">50</real>
            </secondTickValue>
            <endValue resultRef="314">
              <real xmlns="http://schemas.mathsoft.com/math50">400</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

3. 呼叫函數 pexp 與 qexp 分別計算及繪製值 x 的累積機率分佈，以及值 p 的累積機率逆分佈。

<region id="ID0EOZCCB" top="2611.2000000000003" left="38.400000000000006" height="254.20000000000002" width="576" 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="334">
          <traceStyle color="#FF2E3192" symbol="none" line-weight="1" line-style="Dash">線條</traceStyle>
        </trace>
        <trace resultRef="335">
          <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Dash">lines</traceStyle>
        </trace>
        <trace resultRef="336">
          <traceStyle color="#FF068149" symbol="none" line-weight="1" line-style="Dash">線條</traceStyle>
        </trace>
      </traces>
      <graph-size width="476" height="177.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
          <axisLine position="origin" positionticmark="0" legendWidth="72.65" />
          <axisGrid>
            <gridFrequency>11</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="341">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">int</id>
              </math>
              <math resultRef="342">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="338">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue resultRef="340">
              <real xmlns="http://schemas.mathsoft.com/math50">20</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="1000">
          <axisLine position="origin" positionticmark="0" legendWidth="71.24" />
          <axisGrid>
            <gridFrequency>11</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="349">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">F</id>
              </math>
              <math resultRef="350">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="351">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">Fp</id>
              </math>
              <math resultRef="352">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="353">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">Fq</id>
              </math>
              <math resultRef="354">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="344">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <secondTickValue resultRef="346">
              <real xmlns="http://schemas.mathsoft.com/math50">100</real>
            </secondTickValue>
            <endValue resultRef="348">
              <real xmlns="http://schemas.mathsoft.com/math50">1000</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

函數 qexp 的第一個引數是機率且必須小於 1，因此會將 int 除以 100。

相關主題

均等分佈

常態分佈

指數分佈

這是否有幫助?