例: 関数の線形結合

linfit関数を使用して、項の線形和によってデータセットをモデル化します。各項は、適合する関数と linfit によって計算される係数から成ります。

1. 次の行列を定義します。

2. 適合する関数のベクトルを定義します。

モデル関数は y(x) = c0 * ln(x) + c1 * sqrt x +c3 です。

3. linfit 関数を呼び出して、各関数の係数を計算します。これらの係数をベクトル c に格納します。

4. 適合関数 f を、vf 内の関数に c 内の係数を掛け合わせた項の線形結合として定義します。

5. データ点と適合関数をプロットします。

<region id="ID0EULLT" top="1248" left="19.2" width="600" height="277" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot origin-positioning="true" background-type="white" show-axis-expressions="true">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="30">
          <traceStyle color="#FF00008B" symbol="x" line-weight="1" line-style="None">lines</traceStyle>
        </trace>
        <trace resultRef="31">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="484.8" height="219.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="1" end="31">
          <axisLine position="bottom" positionticmark="-1" legendWidth="64.2033333333333" />
          <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="32">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">X</id>
              </math>
              <math resultRef="33">
                <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="4.05" end="8.05">
          <axisLine position="left" positionticmark="-1" legendWidth="85.5766666666667" />
          <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="36">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">Y</id>
              </math>
              <math resultRef="37">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="38">
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <vectorize />
                  <apply>
                    <id xml:space="preserve">f</id>
                    <id xml:space="preserve">X</id>
                  </apply>
                </apply>
              </math>
              <math resultRef="39">
                <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 resultRef="35">
              <real xmlns="http://schemas.mathsoft.com/math50">8</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

関数の線形結合はデータに良好に適合しています。このことは、次に示すように相関係数の値が 1 に近いことによって確認できます。