数值微分和数值积分

数值微分和数组积分都有用函数值近似拟合的形式,它们可以表述成
x i x_i xi的数组微分(积分) = ∑ a k f ( x k ) =\sum a_kf(x_k) =akf(xk)

数值微分

先插值,再对插值函数求导

三个点的一二阶导是:
f ′ ( x 0 ) = 1 2 h [ − 3 f ( x 0 ) + 4 f ( x 1 ) − f ( x 2 ) ] f'(x_0)=\frac{1}{2h}[-3f(x_0)+4f(x_1)-f(x_2)] f(x0)=2h1[3f(x0)+4f(x1)f(x2)]
f ′ ( x 1 ) = 1 2 h [ − f ( x 0 ) + f ( x 2 ) ] f'(x_1)=\frac{1}{2h}[-f(x_0)+f(x_2)] f(x1)=2h1[f(x0)+f(x2)]
f ′ ( x 2 ) = 1 2 h [ f ( x 0 ) − 4 f ( x 1 ) + 3 f ( x 2 ) ] f'(x_2)=\frac{1}{2h}[f(x_0)-4f(x_1)+3f(x_2)] f(x2)=2h1[f(x0)4f(x1)+3f(x2)]

f ′ ′ ( x i ) = 1 h 2 [ f ( x 0 ) − 2 f ( x 1 ) + f ( x 2 ) ] f''(x_i)=\frac{1}{h^2}[f(x_0)-2f(x_1)+f(x_2)] f′′(xi)=h21[f(x0)2f(x1)+f(x2)]

数值积分

牛顿-柯特斯公式
是拉格朗日插值的积分,它最终形式有
∫ a b f ( x ) d x ≈ ( b − a ) ∑ C k ( n ) f ( x k ) \int^b_af(x)dx≈(b-a)\sum C_k^{(n)}f(x_k) abf(x)dx(ba)Ck(n)f(xk)
它要求 x k x_k xk之间等距
C k ( n ) C_k^{(n)} Ck(n)是一个与 a , b , f ( x ) a,b,f(x) a,b,f(x)无关的柯特斯系数,可以查表得到

n C k ( n ) C_k^{(n)} Ck(n)
1 1 2 \frac{1}{2} 21 1 2 \frac{1}{2} 21
2 1 6 \frac{1}{6} 61 4 6 \frac{4}{6} 64 1 6 \frac{1}{6} 61
3 1 8 \frac{1}{8} 81 3 8 \frac{3}{8} 83 3 8 \frac{3}{8} 83 1 8 \frac{1}{8} 81

n=0 是矩形公式
n=1 是梯形公式
n=2 是抛物线公式,也叫辛普森公式

复合积分
分段求插值近似积分

高斯型积分

高斯型积分也是
∫ a b f ( x ) d x ≈ ∑ a k f ( x k ) \int^b_af(x)dx≈\sum a_kf(x_k) abf(x)dxakf(xk)的形式
不过它的节点 x k x_k xk不是等距的,并且和积分区间有关联

∫ 1 − 1 f ( x ) d x ≈ A 0 f ( x 0 ) + A 1 f ( x 1 ) \int^{-1}_1f(x)dx≈A_0f(x_0)+A_1f(x_1) 11f(x)dxA0f(x0)+A1f(x1)
构造方程,待定系数法求解
f ( x ) = 1 − > A 0 + A 1 = 2 f(x)=1 ->A_0+A_1=2 f(x)=1>A0+A1=2
f ( x ) = x − > A 0 x 0 + A 1 x 1 = 0 f(x)=x ->A_0x_0+A_1x_1=0 f(x)=x>A0x0+A1x1=0
f ( x ) = x 2 − > A 0 x 0 2 + A 1 x 1 2 = 2 / 3 f(x)=x^2 ->A_0x_0^2+A_1x_1^2=2/3 f(x)=x2>A0x02+A1x12=2/3
f ( x ) = x 3 − > A 0 x 0 3 + A 1 x 1 3 = 0 f(x)=x^3 ->A_0x_0^3+A_1x_1^3=0 f(x)=x3>A0x03+A1x13=0

高斯节点是任意n次正交多项式的零点
它的 a k a_k ak x k x_k xk可以查表

它和正交多项式有对应关系
在[-1,1]区间,有高斯-勒让德积分
根据权函数的不同还有对应类似的高斯-切比雪夫积分;高斯-埃尔米特积分

你可能感兴趣的:(数值分析学习笔记,数值分析)