辛普森近似求解
- 公式证明
- 任意一个对称区间的一元二次函数定积分
- 拆分
- 求和
- :strawberry: 总结 : 如果我们把六分之一乘进去
- 我们只不过在指定的区间采集数据六个求平均，乘以采集数据区间的微元宽度（历史上不少的手稿用h，翻译为微元高度），之后求和
- 除了首尾项，奇数项权重为4，偶数项权重为2
公式证明
任意一个对称区间的一元二次函数定积分
∫ α β a x 2 + b x + c = a 3 x 3 + b 2 x 2 + c x ∣ α β = a 3 β 3 + b 2 β 2 + c β − ( a 3 α 3 + b 2 α 2 + c α ) = a 3 ( β 3 − α 3 ) + b 2 ( β 2 − α 2 ) + c ( β − α ) = β − α 6 [ 2 a ( β 2 + α β + α 2 ) + 3 b ( β + α ) + 6 c ] = β − α 6 [ ( a α 2 + b α + c ) + ( a α 2 + 2 a α β + a β 2 + 2 b α + 2 b β + 4 c ) + ( a β 2 + b β + c ) ] = β − α 6 [ y ( α ) + 4 y ( α + β 2 ) + y ( β ) ] = β − α 6 ( y 0 + 4 y 1 + y 2 ) \LARGE \begin{aligned} \int_{\alpha}^{\beta}ax^2+bx+c&=\frac{a}{3}x^3+\frac{b}{2}x^2+cx\Bigg|_{\alpha}^{\beta} \\ &=\frac{a}{3}\beta^3+\frac{b}{2}\beta^2+c\beta -(\frac{a}{3}\alpha^3+\frac{b}{2}\alpha^2+c\alpha )\\ &=\frac{a}{3}(\beta^3-\alpha^3)+\frac{b}{2}(\beta^2-\alpha^2)+c(\beta-\alpha) \\ &=\frac{\beta-\alpha}{6}\left[{2a}(\beta^2+\alpha\beta+\alpha^2)+{3b}(\beta+\alpha)+6c\right] \\ &=\frac{\beta-\alpha}{6}\left[(a\alpha^2+b\alpha+c)+(a\alpha^2+2a\alpha\beta+a\beta^2+2b\alpha+2b\beta+4c)+(a\beta^2+b\beta+c)\right] \\ &=\frac{\beta-\alpha}{6}\left[y(\alpha)+4y(\frac{\alpha+\beta}{2})+y(\beta)\right] \\ &=\frac{\beta-\alpha}{6}\left(y_0+4y_1+y_2\right) \end{aligned} ∫αβax2+bx+c=3ax3+2bx2+cx αβ=3aβ3+2bβ2+cβ−(3aα3+2bα2+cα)=3a(β3−α3)+2b(β2−α2)+c(β−α)=6β−α[2a(β2+αβ+α2)+3b(β+α)+6c]=6β−α[(aα2+bα+c)+(aα2+2aαβ+aβ2+2bα+2bβ+4c)+(aβ2+bβ+c)]=6β−α y(α)+4y(2α+β)+y(β) =6β−α(y0+4y1+y2)
拆分
🍐 辛普森要求拆分为偶数项，我们拟合一次抛物线近似会占用两个微元面积，如果我们把区间 ( a , b ) (a,b) (a,b)拆分为n份，那么2份就是:
2 ( b − a ) n = β − α \dfrac{2(b-a)}{n}=\beta-\alpha n2(b−a)=β−α
求和
🍇 主要利用定积分的累加性,拆分为n份，自变量从 x 0 ∼ x n x_0\thicksim x_n x0∼xn,每次拟合的函数都不同，只不过我们每次取三个点模拟二次函数
∫ a b ψ ( x ) d x ≈ ∫ x 0 x 2 f 1 ( x ) d x + ∫ x 2 x 4 f 2 ( x ) d x + ⋯ + ∫ x n − 2 x n f n 2 ( x ) d x = b − a 3 n [ ( y 0 + 4 y 1 + y 2 ) + ( y 2 + 4 y 3 + y 4 ) + ⋯ + ( y n − 2 + 4 y n − 1 + y n ) ] \LARGE \begin{aligned} \int_a^b\psi(x)dx &\approx \int_{x_0}^{x_2}f_{1}(x)dx+\int_{x_2}^{x_4}f_{2}(x)dx+\cdots+\int_{x_{n-2}}^{x_n}f_{\frac{n}{2}}(x)dx\\\\ &=\frac{b-a}{3n}\left[(y_0+4y_1+y_2)+(y_2+4y_3+y_4)+\cdots+(y_{n-2}+4y_{n-1}+y_n)\right] \end{aligned} ∫abψ(x)dx≈∫x0x2f1(x)dx+∫x2x4f2(x)dx+⋯+∫xn−2xnf2n(x)dx=3nb−a[(y0+4y1+y2)+(y2+4y3+y4)+⋯+(yn−2+4yn−1+yn)]
🍓 总结 : 如果我们把六分之一乘进去
∫ a b ψ ( x ) d x ≈ ∫ x 0 x 2 f 1 ( x ) d x + ∫ x 2 x 4 f 2 ( x ) d x + ⋯ + ∫ x n − 2 x n f n 2 ( x ) d x = 2 ( b − a ) n [ ( y 0 + 4 y 1 + y 2 ) 6 + ( y 2 + 4 y 3 + y 4 ) 6 + ⋯ + ( y n − 2 + 4 y n − 1 + y n ) 6 ] \LARGE \begin{aligned} \int_a^b\psi(x)dx &\approx \int_{x_0}^{x_2}f_{1}(x)dx+\int_{x_2}^{x_4}f_{2}(x)dx+\cdots+\int_{x_{n-2}}^{x_n}f_{\frac{n}{2}}(x)dx\\\\ &=\frac{2(b-a)}{n}\left[\frac{(y_0+4y_1+y_2)}{6}+\frac{(y_2+4y_3+y_4)}{6}+\cdots+\frac{(y_{n-2}+4y_{n-1}+y_n)}{6}\right] \end{aligned} ∫abψ(x)dx≈∫x0x2f1(x)dx+∫x2x4f2(x)dx+⋯+∫xn−2xnf2n(x)dx=n2(b−a) 6(y0+4y1+y2)+6(y2+4y3+y4)+⋯+6(yn−2+4yn−1+yn)
我们只不过在指定的区间采集数据六个求平均，乘以采集数据区间的微元宽度（历史上不少的手稿用h，翻译为微元高度），之后求和
除了首尾项，奇数项权重为4，偶数项权重为2
- (2*CLOSE+HIGH+LOW+OPEN)/5， 这是股市中收盘价权重为2，最高价，最低价，开盘价 权重都是1，我们采集了五次。如果微元宽度是天数，将会是累积多空双方的一次很好应用。
- 辛普森法的启示在于采集数据的规律是 14242…41。
