範例：Thiele 內插

函數 > 數據分析 > 內插與預測 > 範例：Thiele 內插

範例：Thiele 內插

在數據集上執行 Thiele 連續分數內插。

Thielecoeff 與 Thiele

使用 Thielecoeff 與 Thiele 函數內插來自有理函數的數據集。

1. 定義內插點數。

2. 定義有理函數，例如 Lorentzian 函數。

3. 取樣函數並將其 x 與 y 值記錄到向量中。

4. 使用 Thielecoeff 函數求解擬合樣本點之連續分數的係數。

5. 呼叫 Thiele 函數執行內插。

<region id="ID0ENHKU" top="556.8" left="38.4" width="218.09" height="25.6" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="88">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id label="None" xml:space="preserve">f_int</id>
        <boundVars>
          <id label="None" xml:space="preserve">x</id>
        </boundVars>
      </function>
      <apply>
        <id label="Function" xml:space="preserve" label-is-contextual="true">Thiele</id>
        <sequence>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">x_data</id>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">c</id>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">x</id>
        </sequence>
      </apply>
    </define>
    <resultFormat>
      <decimal precision="round-trip-decimal-places" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

<region id="ID0E2HKU" top="652.8" left="38.4" width="300.62" height="25.6" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="89">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <function>
        <id label="None" xml:space="preserve">fratint</id>
        <boundVars>
          <id label="None" xml:space="preserve">x</id>
        </boundVars>
      </function>
      <apply>
        <id label="Function" xml:space="preserve" label-is-contextual="true">rationalint</id>
        <sequence>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">x_data</id>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">y_data</id>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">x</id>
        </sequence>
      </apply>
    </define>
    <resultFormat>
      <decimal precision="round-trip-decimal-places" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

Thiele 函數成功，但 rationalint 函數在此內插範圍中的某幾點上失敗。這是因為數據對稱，因而導致在 x 的某些中間值上發生除以零的錯誤：

6. 繪製原始函數、樣本點及內插曲線。

<region id="ID0EFLKU" top="883.2" left="19.2" width="589.606666666667" height="312" 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="92">
          <traceStyle color="#FF0000FF" symbol="rectangle" line-weight="1" line-style="None">線條</traceStyle>
        </trace>
        <trace resultRef="93" num-of-points="500">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="DashDotDot">lines</traceStyle>
        </trace>
        <trace resultRef="94" num-of-points="500">
          <traceStyle color="#FF008000" symbol="none" line-weight="1" line-style="Dot">線條</traceStyle>
        </trace>
      </traces>
      <graph-size width="451.206666666667" height="216" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="1">
          <axisLine position="bottom" positionticmark="-1" legendWidth="139.833333333333" />
          <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="99">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">x_data</id>
              </math>
              <math resultRef="100">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="101">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">x</id>
              </math>
              <math resultRef="102">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue resultRef="96">
              <real xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue resultRef="98">
              <real xmlns="http://schemas.mathsoft.com/math50">1</real>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0.005" end="1.005">
          <axisLine position="origin" positionticmark="0" legendWidth="109.803333333333" />
          <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="103">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">y_data</id>
              </math>
              <math resultRef="104">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="105">
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve">f</id>
                  <id xml:space="preserve">x</id>
                </apply>
              </math>
              <math resultRef="106">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation>
              <math resultRef="107">
                <apply xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve">f_int</id>
                  <id xml:space="preserve">x</id>
                </apply>
              </math>
              <math resultRef="108">
                <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>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

7. 繪製表示 Thiele 內插之殘差的圖形。