Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

区分求積法(Quadratures By Parts)

区分求積法と定積分

関数 f(x) は区間 [a, b] で連続であるとし、f(x) と直線 x=a, x=b および x 軸によって囲まれる面積を S とします。

区分求積法Quadratures

図のように区間 [a, b]n+1 個の点
(1)a=x0, x1, x2, , xk1, xk, , xn=b
で分割して、それぞれの点から
(2)f(x0), f(x1), f(x2), , f(xk1), f(xk), , f(xn)
の各点に引いた垂線によって S を帯状に分けます。さらに xk1xk の間に点 pk をとり、各々の帯の面積を底辺 Δxk=xkxk1, 高さ f(pk) の長方形の面積で近似します。すべての長方形を足し合わせると
(3)Sn=k=1nf(pk)Δxk
と表されます。ここで分割数 n を限りなく小さくすると、Δxk も限りなく小さくなります。このとき SnS に限りなく近づくので、その極限値
(4)limnk=1nf(pk)Δxk
f(x)定積分 (definite integral) と定義して、
(5)abf(x)dx=limnk=1nf(pk)Δxk
のように表します。a, b をそれぞれ積分下限、積分上限といいます。(5) は 不定積分 の記号 dx とは全く別個に定義したものなので、現段階では不定積分との間の関連性は明らかではありません。

定積分の基本性質

定義 (5) より、定積分に関して次のような公式が成り立つことがわかります。
(6)ab{f(x)+g(x)}dx=abf(x)dx+abg(x)dxabkf(x)dx=kabf(x)dx
両式を1つにまとめると、
(7)ab{kf(x)+lg(x)}dx=kabf(x)dx+labg(x)dx
と表せます。また、被積分関数が f(x)=c (c は定数) であるときには、
(8)abcdx=c(ba)
となります。区間 [a, b]f(x)0 ならば
(9)abf(x)dx0
が成り立つことも明らかです。区間 [a, b]f(x)g(x) が成り立つなら、f(x)g(x)0 なので、(9) より
(10)ab{f(x)g(x)}dx0
が成り立つので、
(11)abf(x)dxabg(x)dx
が成立します。また、積分の下限と上限が等しいとき、すなわち b=a のとき、
(12)aaf(x)dx=0
が成り立ちます。

平均値の定理

微分で表された平均値の定理においては、区間 [a, b] において傾きが f(c) (a<c<b) となるような c が必ず存在するとされていました。この記事では長方形 f(c)(ba) が定積分 abf(x)dx と等しくなるような c が存在すると主張する、平均値の定理の積分形を導きます。

積分で表された平均値の定理

いま、f(x) は区間 [a, b] において連続で、区間内における最小値と最大値がそれぞれ m, M であるとします。先ほど証明した定積分の性質
[a, b]
において
f(x)g(x)abf(x)dxabg(x)dx
を用いると、

(13)abmdxabf(x)dxabMdx
という不等式が成り立ちます。すなわち
(14)m(ba)abf(x)dxM(ba)
各辺を ba で割ると
(15)m1baabf(x)dxM
よって、mAM を満たすある定数を A として
(16)1baabf(x)dx=A
となります。中間値の定理により、
(17)f(c)=A(a<c<)
となるような c が必ず存在するので、
(18)1baabf(x)dx=f(c)
を満たす c が存在します。c=a+s(ba) とおくと
(19)abf(x)dx=(ba)f(a+s(ba))(0<s<1)
となり、これが積分で表された平均値の定理です。

微分積分学の基本定理

平均値の定理を使うと、不定積分と定積分の関係を明らかにすることができます。積分上限を変数 x とする関数
(20)F(x)=axf(t)dt
を定義して、変数を x から Δx だけ動かしたときの定積分を計算すると
(21)F(x+Δx)=ax+Δxf(x)dx=axf(t)dt+ax+Δxf(t)dt=F(x)+ax+Δxf(t)dt
ここで先ほどの平均値の定理
abf(x)dx=(ba)f(a+s(ba))(0<s<1)
を使うと、
(22)F(x+Δx)=F(x)+(x+Δxx)f(x+sΔx)=F(x)+Δxf(x+sΔx)
となるので、
(23)F(x+Δx)F(x)Δx=f(x+sΔx)
ここで Δx0 の極限をとると、
(24)F(x)=f(x)
が得られます。これは 微分積分学の基本定理 とよばれています。この定理を少し言い換えると
(25)F(x)=axf(t)dt
f(x) の原始関数のひとつです。いま、F(x)=f(x) を満たすような原始関数のひとつ F0(x) が見つかったとします。F0(x) に任意定数 C を加えてもやはり原始関数なので、
(26)axf(t)dt=F0(x)+C
x=a とおくと左辺は 0 になるので
(27)0=F0(a)+C
したがって、C=F0(a) と定まり、
(28)axf(t)dt=F0(x)F0(a)
x=b と置き直すと
(29)abf(t)dt=F0(b)F0(a)
F0(x) は任意の原始関数なので F(x) と置き直します。また、左辺の積分変数 tx に置き換えると
(30)abf(x)dx=F(b)F(a)
となります。すなわち、下限 a から上限 b の定積分は F(b) から F(a) を引いて得られることになります。積分下限と積分上限を入れ替えると、
(31)abf(x)dx=F(a)F(b)
となるので、
(32)baf(x)dx=abf(x)dx
という公式も導かれます。

≫ Python高水準パッケージ SciPy による数値計算法

エクセルや数学に関するコメントをお寄せください