例: 正弦回帰

sinfit関数を使用して正弦回帰を実行します。この関数は一定周期で変動するデータに特に適しています。

1. データセットを定義します。

2. 推定値のベクトルを入力します。

場合によっては、最良適合を得るために推定値を調整する必要があります。1 つ目の値にはデータの振幅とほぼ等しい値、2 つ目の値には位相シフトと等しい値、3 つ目の値にはオフセットと等しい値を指定する必要があります。

3. sinfit 関数を呼び出して、正弦適合のパラメータを求めます。

これらのパラメータは次の正弦方程式に対応しています。

4. データと正弦適合をプロットします。

<region id="ID0E2EHT" top="883.20000000000016" left="0" height="312" width="622" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot origin-positioning="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="66">
          <traceStyle color="#FFFF0000" symbol="x" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
        <trace resultRef="67">
          <traceStyle color="#FF000000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="483.6" height="216" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="50">
          <axisLine position="origin" positionticmark="5" legendWidth="92.1466666666667" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="68">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">vx</id>
              </math>
              <math resultRef="69">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="70">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">range</id>
              </math>
              <math resultRef="71">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-10" end="10">
          <axisLine position="origin" positionticmark="0" legendWidth="113.52" />
          <axisGrid>
            <gridFrequency>11</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="76">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">vy</id>
              </math>
              <math resultRef="77">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="78">
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve">f</id>
                  <id xml:space="preserve">range</id>
                </apply>
              </math>
              <math resultRef="79">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="73">
              <real xmlns="http://schemas.mathsoft.com/math50">-10</real>
            </startValue>
            <endValue resultRef="75">
              <real xmlns="http://schemas.mathsoft.com/math50">10</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

適合がほぼ完璧であることが次の相関係数によって確認できます。