非均勻梯形值與帶符號誤差

本頁導覽
  1. 目標與範圍
  2. 已知條件與符號
  3. 規則與來源
  4. 逐步主鏈
  5. 局部子步與引理
  6. 逐步反向審查
  7. 錯誤控制與反例
  8. 結論與限制
  9. 驗證及來源綁定

← 導讀目錄 · 原始 MKO

derivation-mko-trapezoidal-quadrature-01 · ocme-new142-reader-adoption-20261021

目標與範圍

對 f(x)=x^2、nodes (0,1/2,2) 計算 T、I-T 與絕對誤差。

允許非均勻分割,完整二次函數給獨立積分;不證一般誤差界或收斂階。

已知條件與符號

  • algebra
  • polynomial
  • trap_def
  • integral_poly
  • G1 — f(x)=x^2,[a,b]=[0,2],nodes (0,1/2,2),values (0,1/4,4)。
  • T — 梯形近似值
  • E=I-T — 固定符號誤差
  • h_i — 每段自己的寬度

規則與來源

  • mko_definition · cited — 來源語義轉述:合法非均勻分割的梯形值逐段以各自寬度乘端點值平均後加總;有限表不決定一般積分。 · mko-trapezoidal-quadrature
  • algebra · prerequisite — 實數域的等式運算;除法前分母必非零,不等式乘正數保持方向。
  • polynomial · prerequisite — 多項式在實數上連續且可逐項微分;d(x^k)/dx=k x^(k-1),常數導數為零。
  • trap_def · prerequisite — 複合梯形值為各段 (x_i-x_(i-1))*(y_(i-1)+y_i)/2 之和;nodes 嚴格遞增、端點值對齊,n 段有 n+1 nodes。
  • integral_poly · prerequisite — 實多項式有限區間積分逐項以 x^(k+1)/(k+1) 原始函數計算。

逐步主鏈

S01 · implies

操作前

0<1/2<2,\quad\#\text{nodes}=3

操作後

n=2,\quad\#\text{values}=3

條件:無

規則:trap_def · 依賴:G1

葉節點依據:trap_def

合法分割有兩段三節點,值陣列長度相符。

S02 · equivalent

操作前

(f(0),f(1/2),f(2))

操作後

(0,1/4,4)

條件:無

規則:algebra · 依賴:G1, S01

葉節點依據:algebra

由完整函數獨立回算全部端點值。

S03 · equivalent

操作前

h_1=x_1-x_0,\quad h_2=x_2-x_1

操作後

h_1=1/2,\quad h_2=3/2

條件:無

規則:algebra · 依賴:G1, S01

葉節點依據:algebra

兩段寬度不同,不能偷換成單一 h。

S04 · implies

操作前

T_1=\tfrac12\tfrac{0+1/4}{2},\quad T_2=\tfrac32\tfrac{1/4+4}{2}

操作後

T_1=1/16,\quad T_2=51/16

條件:無

規則:trap_def · 依賴:S02, S03

葉節點依據:trap_def

每段使用自己的寬度乘自己的端點平均。

S05 · equivalent

操作前

T=T_1+T_2

操作後

T=52/16=13/4

條件:無

規則:algebra · 依賴:S04

葉節點依據:algebra

精確加總兩段。

S06 · implies

操作前

I=\int_0^2x^2\,dx

操作後

I=[x^3/3]_0^2=8/3

條件:無

規則:integral_poly · 依賴:G1

葉節點依據:integral_poly

完整二次函數允許不依樣本值的獨立 reference integral。

S07 · equivalent

操作前

E=I-T=8/3-13/4

操作後

E=-7/12,\quad|E|=7/12

條件:無

規則:algebra · 依賴:S05, S06

葉節點依據:algebra

signed error 固定 exact integral 減梯形值;另外報絕對誤差。

S08 · implies

操作前

\ell(x)=mx+b,\quad u<v

操作後

\int_u^v\ell(x)\,dx=(v-u)\left(m\frac{u+v}{2}+b\right)=(v-u)\frac{\ell(u)+\ell(v)}2

條件:u<v

規則:integral_poly · 依賴:S03

葉節點依據:integral_poly

逐段原始函數恆等式證任意合法分割的 affine exactness,再加總各段;非 affine 本例已給非零誤差。

局部子步與引理

本稿無局部子步。

逐步反向審查

  • S08 · premise_audit · S03 — affine exactness 不是一般二階誤差界或收斂率。
  • S07 · inverse_equivalence · S05, S06 — 負號表示本例梯形值偏大,不把負誤差當負誤差大小。
  • S06 · premise_audit · G1 — 只給節點表時沒有這個完整 integrand 證書。
  • S05 · inverse_equivalence · S04 — T 是 quadrature 值,不先稱 exact integral。
  • S04 · premise_audit · S02, S03 — 中節點在兩段各參與一次,是公式本意。
  • S03 · inverse_equivalence · G1, S01 — 兩寬相加等 2=b-a 作回查。
  • S02 · inverse_equivalence · G1, S01 — 不是未綁函數的樣本表。
  • S01 · premise_audit · G1 — 三個 nodes 不是三個 subintervals。

錯誤控制與反例

每段寬度

先獨立算每個 h_i 再乘端點平均。

誤差符號

固定 I-T,另列絕對值。

  • 三個 nodes 有三段。 → reject — 只有兩對相鄰端點。
  • 本非均勻分割可用同一 h。 → reject — 兩段分別 1/2、3/2。
  • 給三個值就唯一決定積分。 → reject — 可加正倍數 x^2(x-1/2)^2(x-2)^2;節點值不變而積分增加。

結論與限制

T=13/4,I=8/3,I-T=-7/12、絕對誤差 7/12;affine 函數逐段 exact。

    驗證及來源綁定

    3 個本稿宣告算術檢查;只支持所列算式。

    • source commit: b1b2f1dcf8c50d67a93fe52d95f8eef7b3abba9f
    • adopted source SHA-256: 0523238679e2cdf17d8dbeaaa405cec37b031589a6c3c5e271550f3786ebdf06
    • original source SHA-256: 0523238679e2cdf17d8dbeaaa405cec37b031589a6c3c5e271550f3786ebdf06
    • Lean: NOT_RUN · independent AI: SCOPED_LEAD_AND_BILINGUAL_CONTEXT_REVIEW

    同源完整 JSON(含公式、AST、反控制與修訂記錄)