4.8.11 · Maths › Numerical Methods
Intuition Badi picture kya hai
Jab hum ek anjaan function f ( x ) ko ek aisi polynomial p n ( x ) se replace karte hain jo n + 1 data points se guzarti hai, toh hum actually points ke beech mein function ko guess kar rahe hote hain. Error humein batata hai guess kitna bura ho sakta hai .
Key insight yeh hai: error bilkul uss agli term jaisi dikhti hai jo tumne Taylor-jaisi expansion mein include nahi ki — yeh control hoti hai (1) x nodes se kitna door hai, aur (2) f kitni "curvy" hai (uske higher derivative).
Definition Interpolation error
Maano f ek aisi function hai jise n + 1 distinct nodes x 0 , x 1 , … , x n par sample kiya gaya hai, aur p n ( x ) woh unique polynomial hai jiska degree ≤ n hai aur p n ( x i ) = f ( x i ) . Kisi point x par error hai
E ( x ) = f ( x ) − p n ( x ) .
Nodes par E ( x i ) = 0 ; unke beech mein yeh generally nonzero hota hai.
Woh central result jo hum derive karenge:
Hume dikhana hai ki E ( x ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ω ( x ) .
Step 1 — Ek point fix karo aur ek helper function banao.
Koi bhi aisa point x ˉ lo jo node nahi hai (agar node hai toh error 0 hai aur kaam khatam). Hume abhi wahan error nahi pata; ek unknown constant K define karo jo is tarah se defined hai:
f ( x ˉ ) − p n ( x ˉ ) = K ω ( x ˉ ) .
Yeh step kyun? Hum error ko (constant)× ω ki form mein force kar rahe hain. Agar hum prove kar saken ki K = ( n + 1 )! f ( n + 1 ) ( ξ ) , toh formula done hai. Yeh Feynman move hai: answer ki shape assume karo, phir usse pin down karo.
Step 2 — Ek nayi variable t mein ek auxiliary function define karo:
g ( t ) = f ( t ) − p n ( t ) − K ω ( t ) .
Yeh step kyun? Hume ek aisi function chahiye jo bahut saare points par vanish kare taaki hum Rolle repeatedly apply kar sakein.
Step 3 — g ke zeros count karo.
Har node par: f ( x i ) = p n ( x i ) aur ω ( x i ) = 0 , isliye g ( x i ) = 0 . Yeh n + 1 zeros hue.
t = x ˉ par: K ki definition se, g ( x ˉ ) = f ( x ˉ ) − p n ( x ˉ ) − K ω ( x ˉ ) = 0 . Ek aur zero.
Toh g ke kam se kam n + 2 distinct zeros hain [ a , b ] mein.
Step 4 — Rolle's theorem baar baar apply karo.
Kyun? Rolle kehta hai g ke do zeros ke beech g ′ ka ek zero hoga.
g ke ≥ n + 2 zeros hain ⇒ g ′ ke ≥ n + 1 zeros hain.
g ′′ ke ≥ n zeros hain.
… g ( n + 1 ) ka ≥ 1 zero hai. Usse ξ bolo.
Step 5 — g ko exactly n + 1 baar differentiate karo.
g ( n + 1 ) ( t ) = f ( n + 1 ) ( t ) − p n ( n + 1 ) ( t ) − K ω ( n + 1 ) ( t ) .
p n ka degree ≤ n hai ⇒ p n ( n + 1 ) ≡ 0 .
ω ( t ) = t n + 1 + ( lower ) , isliye ω ( n + 1 ) ( t ) = ( n + 1 )! (ek constant hai).
Isliye
g ( n + 1 ) ( t ) = f ( n + 1 ) ( t ) − K ( n + 1 )! .
Step 6 — ξ par evaluate karo.
Kyunki g ( n + 1 ) ( ξ ) = 0 :
0 = f ( n + 1 ) ( ξ ) − K ( n + 1 )! ⇒ K = ( n + 1 )! f ( n + 1 ) ( ξ ) .
Step 7 — Wapas substitute karo. Step 1 se,
f ( x ˉ ) − p n ( x ˉ ) = K ω ( x ˉ ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ω ( x ˉ ) . ■
feel kyun sahi lagta hai
Taylor se compare karo: f ( x ) − T n ( x ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ( x − x 0 ) n + 1 . Interpolation sirf ek hi repeated node ( x − x 0 ) n + 1 ko spread-out nodes ∏ ( x − x i ) se replace karta hai. Wahi ( n + 1 )! f ( n + 1 ) ( ξ ) factor rehta hai — interpolation hai "Taylor jisme derivatives ko data points se trade kiya gaya hai."
Hume ξ rarely pata hota hai, isliye hum isse bound karte hain. Maano M n + 1 = max [ a , b ] ∣ f ( n + 1 ) ( x ) ∣ . Tab
∣ E ( x ) ∣ ≤ ( n + 1 )! M n + 1 ∣ ω ( x ) ∣ , ∣ ω ( x ) ∣ ≤ ∏ i = 0 n ∣ x − x i ∣.
Intuition Error ko control karne ke do knobs
f ( n + 1 ) — wiggly functions (bade high derivatives wali) ka interpolation bura hota hai. Nodes choose karke tum yeh fix nahi kar sakte; yeh function ki galti hai.
ω ( x ) — nodes tum choose karte ho. Unhe achhe se cluster karo aur ∣ ω ∣ chhota rahega. Isliye Chebyshev nodes (ends ke paas zyada dense) worst-case max ∣ ω ∣ ko minimise karte hain, jabki equally spaced nodes endpoints ke paas blow up karte hain (Runge phenomenon ).
Worked example Example 1 — Linear interpolation error (
n = 1 )
Do nodes x 0 , x 1 jahan h = x 1 − x 0 hai. Tab
E ( x ) = 2 ! f ′′ ( ξ ) ( x − x 0 ) ( x − x 1 ) .
Yeh step kyun? n + 1 = 2 toh hume f ′′ chahiye aur ω ( x ) = ( x − x 0 ) ( x − x 1 ) .
Factor ( x − x 0 ) ( x − x 1 ) size mein midpoint x = 2 x 0 + x 1 par maximum hota hai, jahan ∣ ( x − x 0 ) ( x − x 1 ) ∣ = 4 h 2 .
Kyun? Yeh ek downward parabola hai roots x 0 , x 1 ke saath; vertex midpoint par hai.
∴ max ∣ E ∣ ≤ 8 h 2 max ∣ f ′′ ∣.
Yeh classic trapezoidal-rule-style O ( h 2 ) error hai.
Worked example Example 2 — Numbers:
f ( x ) = ln x ko [ 1 , 2 ] par linearly interpolate karo
f ′′ ( x ) = − 1/ x 2 , isliye max [ 1 , 2 ] ∣ f ′′ ∣ = ∣ − 1/ 1 2 ∣ = 1 (sabse bada x = 1 par).
h = 1 . Bound: ∣ E ∣ ≤ 8 h 2 ⋅ 1 = 8 1 = 0.125.
x = 1.5 par check karo: sahi ln 1.5 = 0.405465 ; line p ( x ) = ln 1 + ( ln 2 − ln 1 ) ( x − 1 ) = 0.693147 ( x − 1 ) deti hai p ( 1.5 ) = 0.346574 .
Actual error = 0.405465 − 0.346574 = 0.058891 , jo waqai ≤ 0.125 hai. ✓
Bound kyun loose hai: ξ near 1.5 deta hai f ′′ ≈ − 1/1. 5 2 = − 0.444 , worst-case − 1 nahi.
Worked example Example 3 — Quadratic interpolation (
n = 2 )
Teen equally spaced nodes x 0 , x 0 + h , x 0 + 2 h .
E ( x ) = 6 f ′′′ ( ξ ) ( x − x 0 ) ( x − x 0 − h ) ( x − x 0 − 2 h ) .
∣ ω ∣ ko interval par maximise karne se (cubic par calculus) milta hai max ∣ ω ∣ = 3 3 2 h 3 .
Yeh step kyun? ω differentiate karo, zero set karo, interior extremum dhundo, wapas plug karo. Isliye
max ∣ E ∣ ≤ 3 3 ⋅ 6 2 h 3 max ∣ f ′′′ ∣ = 9 3 h 3 max ∣ f ′′′ ∣ = O ( h 3 ) .
Pattern: n + 1 nodes ⇒ error O ( h n + 1 ) smooth f ke liye.
Common mistake "Zyada nodes matlab hamesha chhoti error."
Kyun sahi lagta hai: zyada data ⇒ better fit, jaise zyada pixels add karna.
Kyun galat hai: factor hai ( n + 1 )! f ( n + 1 ) . Halanki ( n + 1 )! badhta hai, max ∣ f ( n + 1 ) ∣ tez badhegi (jaise f ( x ) = 1 + 25 x 2 1 ), aur equally spaced nodes ke saath max ∣ ω ∣ endpoints ke paas explode karta hai. Result: Runge phenomenon — error badhti hai.
Fix: Chebyshev nodes use karo ya piecewise (spline) interpolation; blindly n mat badhao.
ξ ek single known number hai jo main compute kar sakta hoon."
Kyun sahi lagta hai: formula ek ξ likhta hai.
Kyun galat hai: ξ = ξ ( x ) , x par depend karta hai aur sirf exist karne ki guarantee hai (existence theorem), explicitly given nahi.
Fix: practical bounds ke liye, f ( n + 1 ) ( ξ ) ko [ a , b ] par uske maximum se replace karo.
( n + 1 )! bhool jaana, E = f ( n + 1 ) ( ξ ) ω ( x ) likhna.
Kyun sahi lagta hai: leading remainder term jaisa dikhta hai.
Fix: n + 1 Rolle differentiations se exactly ω ( n + 1 ) = ( n + 1 )! denominator mein aata hai (Step 5). Isse kabhi mat choodo.
Recall Feynman: ek 12-saal ke bachche ko explain karo
Socho tum ek graph par dots ko ek smooth curve se connect kar rahe ho. Tum sirf dots ki height jaante ho. Dots ke beech mein tum guess kar rahe ho. Guess kitna galat ho sakta hai? Do cheezein matter karti hain: tum nearest dots se kitna door ho (ω part — jitne zyada dots paas mein, utna safe), aur asli line kitni twisty hai (derivative part — seedhi road guess karna easy hai, roller-coaster nahi). Formula bas "twistiness" ko "dots se distance" se multiply karta hai aur ek bade factorial se divide karta hai jo usse soften karta hai. Surprisingly, dots ko edges ke paas zyada tight pack karna (Chebyshev) worst guess ko bahut better banata hai.
Mnemonic Formula yaad karo
"Next Derivative over Next Factorial, times the product of (x minus nodes)."
Symbol shape: ( n + 1 )! f ( n + 1 ) ( ξ ) ω ( x ) — ( n + 1 ) recurring star hai. Socho: "apni degree se ek zyada, har jagah."
Error formula aur smoothness condition state karo.
Kaunsa factor tum control karte ho, aur kaise?
n badhane se kabhi kabhi cheezein worse kyun hoti hain?
Degree-n interpolation ke liye interpolation error formula kya hai? E ( x ) = f ( x ) − p n ( x ) = ( n + 1 )! f ( n + 1 ) ( ξ ) ∏ i = 0 n ( x − x i ) , jahan ξ ∈ ( a , b ) .
Error formula ke liye smoothness condition kya chahiye? f ∈ C n + 1 [ a , b ] , yaani f ke n + 1 continuous derivatives hain.
Node polynomial ω ( x ) kya hai? ω ( x ) = ∏ i = 0 n ( x − x i ) , jo har interpolation node par vanish karta hai.
Derivation mein p n ( n + 1 ) ≡ 0 kyun hota hai? Kyunki p n ka degree ≤ n hai, n + 1 baar differentiate karne par yeh zero ho jaata hai.
Denominator mein ( n + 1 )! kahan se aata hai? ω ( t ) = t n + 1 + … isliye ω ( n + 1 ) ( t ) = ( n + 1 )! , jo Rolle's theorem n + 1 baar apply karne par milta hai.
Helper g ( t ) ke kitne zeros hote hain, aur kyun? Kam se kam n + 2 : n + 1 nodes mein se har ek par ek, plus chosen point x ˉ par ek.
Spacing h ke saath linear interpolation ka error bound? ∣ E ∣ ≤ 8 h 2 max ∣ f ′′ ∣ .
Do nodes ke liye ∣ ω ∣ ka maximum kahan hota hai? Midpoint par, value h 2 /4 hoti hai.
Runge phenomenon kya hai? Equally spaced nodes ke saath, degree badhane par interpolation error endpoints ke paas blow up ho sakta hai (jaise 1/ ( 1 + 25 x 2 ) ke liye).
Interpolation error ki size ko kaunse do factors control karte hain? High derivative f ( n + 1 ) (function ki curviness, control nahi hoti) aur ω ( x ) (node placement, control hoti hai).
Chebyshev nodes se faida kyun hota hai? Yeh endpoints ke paas cluster hote hain aur max ∣ ω ( x ) ∣ ko minimise karte hain, worst-case error kam karte hain aur Runge oscillations avoid karte hain.
n + 1 equally spaced nodes ke liye, error h mein kis order ka hota hai?Sufficiently smooth f ke liye O ( h n + 1 ) .
Lagrange Interpolation — woh p n banata hai jiska yeh error hai.
Newton Divided Differences — error term equals next divided difference f [ x 0 , … , x n , x ] ω ( x ) .
Taylor's Theorem with Remainder — limiting case jab saare nodes coalesce ho jaate hain.
Rolle's Theorem / Mean Value Theorem — derivation ka engine.
Runge Phenomenon aur Chebyshev Nodes — node-choice ke consequences.
Numerical Integration — Newton–Cotes error is ω factor ko inherit karta hai.
p_n interpolating polynomial
Error E(x) = f(x) - p_n(x)
Helper g(t) = f - p_n - K·omega
p_n^(n+1)=0 and omega^(n+1)=(n+1)!
Error formula: f^(n+1)(xi)/(n+1)! · omega(x)