差商

テンプレート:出典の明記 数学における差商（さしょう、テンプレート:Lang-en-short; 分割差分、差分商）は、差分商をとる操作を再帰的に繰り返すことで与えられる。歴史的には対数表や三角函数表の計算に用いられ、チャールズ・バベッジの階差機関（初期の機械式計算機）はこれを実装するものとして設計された。

差商はニュートン補間における補間多項式の計算に用いることができる。

定義

テンプレート:Math 個の節点 $(x_0,y_0),\dots ,(x_n,y_n)$ に対する、前進差商は $$\begin{array}{cc}[y_{\nu }]& :=y_{\nu }(\nu \in \{0,\dots ,n\}),\\ [y_{\nu },\dots ,y_{\nu +j}]& :=\frac{[y_{\nu +1},\dots ,y_{\nu +j}]-[y_{\nu },\dots ,y_{\nu +j-1}]}{x_{\nu +j}-x_{\nu }}(\nu \in \{0,\dots ,n-j\},j\in \{1,\dots ,n\})\end{array}$$ で定義される。同様に後退差商は $$\begin{array}{cc}[y_{\nu }]& :=y_{\nu }(\nu \in \{0,\dots ,n\}),\\ [y_{\nu },\dots ,y_{\nu -j}]& :=\frac{[y_{\nu },\dots ,y_{\nu -j+1}]-[y_{\nu -1},\dots ,y_{\nu -j}]}{x_{\nu }-x_{\nu -j}}(\nu \in \{j,\dots ,n\},j\in \{1,\dots ,n\})\end{array}$$ と定義される。

以下本項では主に前進差商のみを扱い、それを単に差商と呼ぶ。

等間隔の場合

テンプレート:Main 節点（の テンプレート:Mvar-座標）が等間隔に並んでいるときには、前進差分（あるいは後退・中央の各有限差分）の定める差分商によって差商を記述することができる。この場合は、より一般の差商よりも計算が楽になる。注目すべきは、差商の分母が既知であるならば、差分商から差商を容易に回復することができることである。

すなわち、上と同じ節点のデータが与えられ、適当な数 テンプレート:Math に対して テンプレート:Math テンプレート:Math となっているとき、前進差分は $$\begin{array}{cc}{\mathrm{\Delta }}^{(0)}{y}_i& :=y_i,\\ {\mathrm{\Delta }}^{(k)}{y}_i& :={\mathrm{\Delta }}^{(k-1)}{y}_{i+1}-{\mathrm{\Delta }}^{(k-1)}{y}_i(k\ge 1)\end{array}$$ で定義される。前進差分と差商との関係は $$f[x_0,x_1,\dots ,x_k]=\frac{1}{k!h^k}{\mathrm{\Delta }}^{(k)}f(x_0)$$ と書ける。

ニュートンの三角形

隣り合う二項から次の階数の項が次々に作られる様子は、パスカルの三角形の類似で、差商を三角形状に並べるとわかりよい。

$\begin{array}{ccccc}{x}_0& y_0=[y_0]& & & \\ & & [y_0,y_1]& & \\ x_1& y_1=[y_1]& & [y_0,y_1,y_2]& \\ & & [y_1,y_2]& & [y_0,y_1,y_2,y_3]\\ x_2& y_2=[y_2]& & [y_1,y_2,y_3]& \\ & & [y_2,y_3]& & \\ x_3& y_3=[y_3]& & & \end{array}$ あるいは $\begin{array}{cccc}{y}_0& & & \\ & \mathrm{\Delta }{y}_0& & \\ y_1& & {\mathrm{\Delta }}^2{y}_0& \\ & \mathrm{\Delta }{y}_1& & {\mathrm{\Delta }}^3{y}_0\\ y_2& & {\mathrm{\Delta }}^2{y}_1& \\ & \mathrm{\Delta }{y}_2& & \\ y_3& & & \end{array}$

記法

節点が函数 テンプレート:Mvar を用いて $(x_0,f(x_0)),\dots ,(x_n,f(x_n))$ と与えられているときには、$$\begin{array}{cc}f[x_{\nu }]& :=f(x_{\nu }),\nu \in \{0,\dots ,n\},\\ f[x_{\nu },\dots ,x_{\nu +j}]& :=\frac{f[x_{\nu +1},\dots ,x_{\nu +j}]-f[x_{\nu },\dots ,x_{\nu +j-1}]}{x_{\nu +j}-x_{\nu }}(\nu \in \{0,\dots ,n-j\},j\in \{1,\dots ,n\})\end{array}$$ と書くこともある。同様に $$[x_0,\dots ,x_n]f,[x_0,\dots ,x_n;f],D[x_0,\dots ,x_n]f$$ などとも書く。

具体的に書き下せば: $$\begin{array}{cc}f[x_0]& =f(x_0)\\ f[x_0,x_1]& =\frac{f(x_0)}{(x_0-x_1)}+\frac{f(x_1)}{(x_1-x_0)}\\ f[x_0,x_1,x_2]& =\frac{f(x_0)}{(x_0-x_1)\cdot (x_0-x_2)}+\frac{f(x_1)}{(x_1-x_0)\cdot (x_1-x_2)}+\frac{f(x_2)}{(x_2-x_0)\cdot (x_2-x_1)}\\ f[x_0,x_1,x_2,x_3]& =\frac{f(x_0)}{(x_0-x_1)\cdot (x_0-x_2)\cdot (x_0-x_3)}+\frac{f(x_1)}{(x_1-x_0)\cdot (x_1-x_2)\cdot (x_1-x_3)}+\frac{f(x_2)}{(x_2-x_0)\cdot (x_2-x_1)\cdot (x_2-x_3)}+\frac{f(x_3)}{(x_3-x_0)\cdot (x_3-x_1)\cdot (x_3-x_2)}\end{array}$$ のようになる。一般に $$f[x_0,\dots ,x_n]=\sum _{j=0}^n\frac{f(x_j)}{\prod _{k=0,\dots n;k\ne j}(x_j-x_k)}$$ となり、これを多項式函数 $q(\xi )=(\xi -x_0)\cdots (\xi -x_n)$ を用いて$$f[x_0,\dots ,x_n]=\sum _{j=0}^n\frac{f(x_j)}{q^{\prime }(x_j)}$$ と書くことができる。また極限を用いて $$f[x_0,\dots ,x_n]=\sum _{j=0}^n\underset{x\to x_j}{\lim }\left[\frac{f(x_j)(x-x_j)}{\prod _{k=0}^n(x-x_k)}\right]$$ とも書ける。

性質

• 線型性: $$\begin{array}{cc}& (f+g)[x_0,\dots ,x_n]=f[x_0,\dots ,x_n]+g[x_0,\dots ,x_n]\\ & (\lambda \cdot f)[x_0,\dots ,x_n]=\lambda \cdot f[x_0,\dots ,x_n]\end{array}$$
• ライプニッツの法則: $$(f\cdot g)[x_0,\dots ,x_n]=f[x_0]\cdot g[x_0,\dots ,x_n]+f[x_0,x_1]\cdot g[x_1,\dots ,x_n]+\dots +f[x_0,\dots ,x_n]\cdot g[x_n]$$
• 対称性: $$f[x_0,\dots ,x_n]=f[x_{\sigma (0)},\dots ,x_{\sigma (n)}](\sigma \in S_n)$$
• 差商に対する平均値の定理より: $$f[x_0,\dots ,x_n]=\frac{f^{(n)}(\xi )}{n!}(\mathrm{\exists }\xi \in (\min x_k,\max x_k))$$

別定義

ペアノ形

次数 テンプレート:Math のBスプライン テンプレート:Mvar を用いて、差商を $$f[x_0,\dots ,x_n]=\frac{1}{n!}\underset{x_0}{\overset{x_n}{\int }}{f}^{(n)}(t)B_{n-1}(t)dt$$と書くことができる。テンプレート:Math は テンプレート:Mvar の テンプレート:Mvar-階導函数である。

これを差商のペアノ形 (Peano form) と言い、テンプレート:Mvar をこの差商のテンプレート:Ill2と呼ぶ。ともに名称はジゼッペ・ペアノに因む。

テイラー形

考えている節点が集積しているならば、ほとんどゼロに近い値での割り算が生じ、相対誤差が大きくなるから、数値計算はおぼつかない。しかしその場合も、差分商を微分商で（あるいはその逆に）近似することはできる:$$\frac{f(y)-f(x)}{y-x}\approx f^{\prime }(x)(x\approx y).$$

この近似はテイラーの定理 $$f(y)=f(x)+f^{\prime }(x)\cdot (y-x)+f^{″}(x)\cdot \frac{(y-x{)}^2}{2!}+f^{‴}(x)\cdot \frac{(y-x{)}^3}{3!}+\cdots$$ が適用できる函数に関しては等式 $$\frac{f(y)-f(x)}{y-x}=f^{\prime }(x)+f^{″}(x)\cdot \frac{y-x}{2!}+f^{‴}(x)\cdot \frac{(y-x{)}^2}{3!}+\cdots$$に変えることができる。また、中央差分を利用すれば テンプレート:Mvar の奇数冪の項を消すことができる: $$\frac{f(y)-f(x)}{y-x}=\frac{f(m+h)-f(m-h)}{2\cdot h}=f^{\prime }(m)+f^{‴}(m)\cdot \frac{h^2}{3!}+\cdots .$$

原理的には、テイラー級数および任意のテンプレート:Ill2に差商近似が適用できる。テイラー級数は冪函数に関する無限和で、函数にその差商を対応させる写像 テンプレート:Math は線型汎函数であるから、基底函数にこの汎函数を適用すればよい。

通常の冪函数 テンプレート:Math によって、通常のテイラー級数を $$f=f(0)\cdot p_0+f^{\prime }(0)\cdot p_1+\frac{f^{″}(0)}{2!}\cdot p_2+\frac{f^{‴}(0)}{3!}\cdot p_3+\dots$$ と書けば、差商のテイラー級数は $$f[x_0,\dots ,x_n]=f(0)\cdot p_0[x_0,\dots ,x_n]+f^{\prime }(0)\cdot p_1[x_0,\dots ,x_n]+\frac{f^{″}(0)}{2!}\cdot p_2[x_0,\dots ,x_n]+\frac{f^{‴}(0)}{3!}\cdot p_3[x_0,\dots ,x_n]+\cdots$$ と書けることになる。この最初の テンプレート:Mvar 項は多項式の次数よりも高階の差分だから消えており、それ以降の項も以下のように知ることができる: $$\{\begin{array}{cc}{p}_j[x_0,\dots ,x_n]=0& (j<n)\\ p_n[x_0,\dots ,x_n]=1& \\ p_{n+1}[x_0,\dots ,x_n]& =x_0+\dots +x_n\\ p_{n+m}[x_0,\dots ,x_n]& =\sum _{0\le a_1\le a_2\le \dots \le a_m\le n}\prod _{j\in \{a_1,a_2,\dots ,a_m\}}{x}_j.\end{array}$$

行列表示

多項式の差商はライプニッツの法則の恩恵を受けられる点で興味を持たれる。

行列 テンプレート:Mvar を $$J=\left(\begin{array}{cccccc}{x}_0& 1& 0& 0& \cdots & 0\\ 0& x_1& 1& 0& \cdots & 0\\ 0& 0& x_2& 1& & 0\\ ⋮& ⋮& & \ddots & \ddots & \\ ⋮& ⋮& & & \ddots & 1\\ 0& 0& 0& 0& \cdots & x_n\end{array}\right)$$ とすれば、これは節点 テンプレート:Math に関する恒等写像 テンプレート:Math の差商をすべて包摂したものと考えることができる。特に テンプレート:Mvar は冪函数 テンプレート:Mvar に対する差商が全て現れる。差商の線型性により、多項式 テンプレート:Mvar に対する多項式函数を テンプレート:Mvar に適用したものは、$$\begin{array}{cc}p(J)& =a_0\cdot I+a_1\cdot J+\dots +a_n\cdot J^n\\ & =\left(\begin{array}{ccccc}p[x_0]& p[x_0,x_1]& p[x_0,x_1,x_2]& \dots & p[x_0,\dots ,x_n]\\ 0& p[x_1]& p[x_1,x_2]& \dots & p[x_1,\dots ,x_n]\\ ⋮& \ddots & \ddots & \ddots & ⋮\\ 0& \dots & 0& 0& p[x_n]\end{array}\right)\end{array}$$ と書くことができ、オピッツの公式と呼ばれる。

多項式 テンプレート:Mvar の次数を無限大に飛ばして形式冪級数とすれば、テイラー展開可能な函数 テンプレート:Mvar に対しても、その差商を テンプレート:Mvar を用いた行列表示 $$T_f(x_0,\dots ,x_n)=f(J)=\left(\begin{array}{ccccc}f[x_0]& f[x_0,x_1]& f[x_0,x_1,x_2]& \dots & f[x_0,\dots ,x_n]\\ 0& f[x_1]& f[x_1,x_2]& \dots & f[x_1,\dots ,x_n]\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ 0& 0& 0& \dots & f[x_n]\end{array}\right)$$ によって与えることができる。

節点 テンプレート:Math が全て等しいとき、テンプレート:Mvar はジョルダン細胞であり、ジョルダン標準形を考えることでスカラー函数を行列函数に一般化することができる。

関連項目

• 差分商
• テンプレート:Ill2
• 多項式補間
• 差商に関する平均値定理
• ネールント–ライス積分
• パスカルの三角形

参考文献

• テンプレート:Cite book
• テンプレート:Cite book
• テンプレート:Cite book

外部リンク

• テンプレート:MathWorld
• テンプレート:PlanetMath
• An implementation in Haskell.