数值微分和数组积分都有用函数值近似拟合的形式,它们可以表述成
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≈(b−a)∑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)dx≈∑akf(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)
∫1−1f(x)dx≈A0f(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]区间,有高斯-勒让德积分
根据权函数的不同还有对应类似的高斯-切比雪夫积分;高斯-埃尔米特积分