4.10.6 · D3 · Maths › Advanced Topics (Elite Level) › Residue theorem — computing real integrals
Intuition Yeh page kya hai
Parent note ne machine banai. Yahan hum usse har tarah ke terrain mein drive karte hain . Har example se pehle tum ek forecast karo (answer ki shape guess karo), phir hum steps grind karte hain, phir verify karte hain. End tak tumne har sign, har quadrant mein pole placement, har degenerate case, aur do "gotcha" cases dekh liye honge jo naive recipe ko tod dete hain.
Definition Is poori page par use hone wala notation
Yahan ke zyaadatar integrands rational functions hain: ek ratio Q ( x ) P ( x ) jahan ==P numerator polynomial hai aur Q denominator polynomial== hai. Symbol deg P ka matlab hai P ki degree — usmein x ki sabse badi power (jaise deg ( x 2 + 1 ) = 2 ). Gap ek single number hai:
gap = deg Q − deg P .
Yeh measure karta hai ki denominator, numerator se door jaake kitni tezi se badhta hai — aur yahi decide karta hai ki hamare contour ka bada arc khatam hoga ya nahi. Gap 2 ya usse zyada hone ka matlab hai Q P ki shrinking ∣ z ∣ 2 1 ya usse tez hoti hai.
Jab ek oscillating factor e ia z present ho, hum integrand ko R ( z ) e ia z mein split karte hain, jahan ==R wo rational part hai jo exponential ko peel karne ke baad bachti hai==. R wahi hai jo Jordan's Lemma maangta hai ki infinity par vanish ho.
Is chapter ke har real integral ko inhi cells mein se ek mein rakha ja sakta hai. Neeche ka har worked example us cell ke saath labelled hai jo woh fill karta hai.
Cell
Kya special hai
Chhupi hui danger
Example
A Basic rational
gap = deg Q − deg P ≥ 2 , simple poles
galat half-plane choose karna
Ex 1
B Higher-order pole
denominator mein repeated factor hai
residue ke liye differentiate karna padega
Ex 2
C Oscillatory (cos/sin)
integrand mein cos a x ya sin a x hai
e ia z use karo, cos z nahi
Ex 3
D Odd integrand → sine
"imaginary part" branch
real part 0 hai, sine answer hai
Ex 4
E Trig over [ 0 , 2 π ]
already ek loop hai → unit circle
mandatory i z 1 factor
Ex 5
F Degenerate: pole on the real axis
contour ek spike se guzarti hai
indent karo — principal value lo
Ex 6
G Gap = 1 with no oscillation
arc nahi marta
recipe fail — ise pehchano
Ex 7
H Word problem
physics/engineering dressing
ise ek jaane-pehchane cell mein translate karo
Ex 8
Ab hum sabhi aath par jaate hain. Labels dekhte rehna.
I = ∫ − ∞ ∞ x 4 + 1 x 2 d x
Forecast: Integrand even aur positive hai, isliye answer ek positive real number hona chahiye, aur (aane wale algebra se) usme 2 hoga. Aage padhne se pehle guess karo.
Step 1 — check karo ki arc marta hai ya nahi. Yeh step kyun? Boxed formula tabhi legal hai jab bada semicircle kuch contribute na kare. Yahan P = x 2 hai toh deg P = 2 , aur Q = x 4 + 1 hai toh deg Q = 4 ; gap = 4 − 2 = 2 ≥ 2 . ✓ Arc marta hai, machine legal hai.
Step 2 — poles dhundho. Kyun? Hume har wo spike chahiye jo contour trap kar sake. Solve karo z 4 + 1 = 0 ⇒ z 4 = − 1 = e iπ . Chaar fourth-roots hain:
z k = e i ( π + 2 π k ) /4 , k = 0 , 1 , 2 , 3 ,
yaani angles 4 5 ∘ , 13 5 ∘ , 22 5 ∘ , 31 5 ∘ .
Figure mein kya dekhna hai: chaar poles ek circle par evenly spaced baithe hain. Yellow path haara contour hai — real axis plus upper arc Γ ρ . Do pink dots (4 5 ∘ aur 13 5 ∘ ) loop ke andar hain (positive imaginary part), toh ye count hote hain; do blue crosses (22 5 ∘ , 31 5 ∘ ) axis ke neeche hain, loop ke bahar , toh inhe ignore karte hain. Jo do bachte hain:
z 0 = e iπ /4 = 2 1 + i , z 1 = e i 3 π /4 = 2 − 1 + i .
Step 3 — residues (simple poles, g / h ′ use karo). g / h ′ kyun? Quartic ko factor karna ugly hai; formula Res = g ( z 0 ) / h ′ ( z 0 ) use karne se bach jaate hain. Yahan g = z 2 , h = z 4 + 1 , h ′ = 4 z 3 , toh
Res z k = 4 z k 3 z k 2 = 4 z k 1 .
Step 4 — har UHP pole par evaluate karo. Kyunki z k unit circle par hai, z k 1 = z k = e − i a r g z k . Toh
Res z 0 = 4 1 e − iπ /4 , Res z 1 = 4 1 e − i 3 π /4 .
Step 5 — sum karo aur 2 π i se multiply karo. Kyun? Yahi residue theorem ka verdict hai. Har exponential ko rectangular form mein likho taaki add kar sakein:
e − iπ /4 = 2 1 − i , e − i 3 π /4 = 2 − 1 − i .
Yeh kyun help karta hai: real parts + 2 1 aur − 2 1 cancel ho jaate hain, aur imaginary parts − 2 i add ho jaate hain, ek purely imaginary sum bachta hai:
e − iπ /4 + e − i 3 π /4 = 2 − 2 i = − i 2 .
Isliye residue sum hai 4 1 ( − i 2 ) = − 4 i 2 , aur
I = 2 π i ⋅ ( − 4 i 2 ) = 4 2 π 2 = 2 π = 2 π 2 .
(Note karo 2 π i ⋅ ( − i ) = 2 π , kyunki − i 2 = 1 — yahi imaginary sum ko real answer mein badalta hai.)
Verify: Answer positive hai aur 2 carry karta hai — forecast se match karta hai. Numerically π / 2 ≈ 2.221 , aur height ≤ 2 1 wale bump ka ∫ kuch units ki width ke saath plausibly ≈ 2.2 hai. ✓
I = ∫ − ∞ ∞ ( x 2 + 4 ) 2 d x
Forecast: Even, positive, aur poles ± i ki jagah ± 2 i par hain — answer ∫ ( x 2 + 1 ) 2 d x = π /2 case se chota hona chahiye kyunki bump flatter hai. π ka ek fraction guess karo.
Step 1 — arc check. P = 1 toh deg P = 0 ; Q = ( x 2 + 4 ) 2 toh deg Q = 4 ; gap = 4 − 0 = 4 ≥ 2 . Arc marta hai. ✓
Step 2 — poles. ( z 2 + 4 ) 2 = ( z − 2 i ) 2 ( z + 2 i ) 2 . UHP pole hai z = 2 i , order m = 2 . Order 2 kyun? Factor ( z − 2 i ) squared appear karta hai.
Step 3 — order-2 residue formula. Ek baar differentiate kyun? Order-m pole ke liye hum ( z − 2 i ) 2 f ka m − 1 = 1 derivative lete hain:
Res 2 i = 1 ! 1 lim z → 2 i d z d [ ( z − 2 i ) 2 ⋅ ( z − 2 i ) 2 ( z + 2 i ) 2 1 ] = d z d ( z + 2 i ) − 2 2 i .
Step 4 — derivative karo.
d z d ( z + 2 i ) − 2 = ( z + 2 i ) 3 − 2 .
z = 2 i par: ( 2 i + 2 i ) 3 = ( 4 i ) 3 = 64 i 3 = − 64 i , toh Res 2 i = − 64 i − 2 = 32 i 1 .
Step 5 — assemble karo.
I = 2 π i ⋅ 32 i 1 = 32 2 π = 16 π .
Verify: π /16 ≈ 0.196 , jo π /2 ≈ 1.571 se kaafi chota hai — "flatter bump" forecast se match karta hai. Sanity scaling: ∫ ( x 2 + 4 ) 2 d x mein x = 2 u substitute karne se 2 3 1 ∫ ( u 2 + 1 ) 2 d u = 8 1 ⋅ 2 π = 16 π milta hai. ✓
I = ∫ − ∞ ∞ x 2 + 4 cos ( 3 x ) d x
Forecast: Positive (tail of cos , 1/ x 2 se damp hai aur 0 ke paas peak dominate karta hai jahan cos ≈ 1 ). Decay factor e − a ⋅ 2 = e − 6 ise bahut chota banata hai. Kuch 2 π e − 6 jaisa guess karo.
Step 1 — cos 3 x ko e i 3 z se replace karo. cos z kyun nahi? Kyunki cos ( 3 z ) = 2 1 ( e i 3 z + e − i 3 z ) aur e − i 3 z UHP mein explode karta hai (uska modulus e + 3 y hai). Sirf e i 3 z wahan decay karta hai (Jordan's Lemma ). Toh set karo f ( z ) = z 2 + 4 e i 3 z , a = 3 .
Step 2 — Jordan's Lemma legality. Yahan exponential ko peel karne ke baad jo rational leftover bachta hai woh hai R ( z ) = z 2 + 4 1 (notation box se yaad karo: R woh non-oscillating part hai). Yeh ∣ z ∣ → ∞ par → 0 hai, aur a = 3 > 0 . Jordan kehta hai arc integral → 0 — chahe gap sirf 2 ho, hume bonus milta hai ki yeh gap 1 par bhi kaam karta. ✓
Step 3 — UHP pole & residue. z = 2 i (simple). g = e i 3 z , h ′ = 2 z :
Res 2 i = 2 ( 2 i ) e i 3 ( 2 i ) = 4 i e − 6 .
e − 6 kyun? i ⋅ 3 ⋅ 2 i = 6 i 2 = − 6 .
Step 4 — full complex integral, phir real part lo.
∫ − ∞ ∞ x 2 + 4 e i 3 x d x = 2 π i ⋅ 4 i e − 6 = 2 π e − 6 .
Yeh already real hai, toh Real part = 2 π e − 6 aur (imag part 0 ) ∫ x 2 + 4 sin 3 x d x = 0 — consistent hai kyunki woh integrand odd hai.
Verify: 2 π e − 6 ≈ 1.5708 × 0.00248 ≈ 0.00389 — positive aur tiny, exactly as forecast. ✓
I = ∫ − ∞ ∞ x 2 + 1 x sin x d x
Forecast: Note karo ki x sin x even hai (odd × odd), toh integral nonzero aur positive hai. Yeh woh case hai jahan sine (imaginary part) answer carry karta hai. Kuch π / e -jaisa guess karo.
Step 1 — complex integrand banao. Kyun? sin x = Im ( e i x ) , toh use karo f ( z ) = z 2 + 1 z e i z , a = 1 .
Step 2 — arc / Jordan check. Exponential peel karne ke baad rational leftover hai R ( z ) = z 2 + 1 z → 0 jab ∣ z ∣ → ∞ (gap = 2 − 1 = 1 ), aur a = 1 > 0 . Jordan's Lemma apply hota hai — arc exactly isliye marta hai kyunki exponential usse kill karta hai chahe gap 1 hi ho. Issi liye hume Jordan chahiye tha, naive gap-2 rule nahi.
Step 3 — UHP pole & residue. z 2 + 1 = 0 ⇒ z = i (UHP). Simple pole, g = z e i z , h ′ = 2 z :
Res i = 2 i i e i ⋅ i = 2 e − 1 .
Kyun: e i ⋅ i = e i 2 = e − 1 , aur g mein i upar-neeche 2 i ke saath cancel hokar 2 1 banta hai.
Step 4 — complex integral assemble karo.
∫ − ∞ ∞ x 2 + 1 x e i x d x = 2 π i ⋅ 2 e − 1 = e π i .
Step 5 — imaginary part lo. Imaginary kyun? Kyunki hume sin chahiye, aur sin x = Im ( e i x ) , toh
I = Im ( e π i ) = e π .
Real part 0 hai — yeh match karta hai ki x 2 + 1 x cos x odd hai, toh uska integral zero hota hai.
Verify: π / e ≈ 1.156 , forecast ke anusaar positive. Cross-check: pairing "real part → cos → 0 (odd), imag part → sin → answer" exactly wahi parity bookkeeping hai. ✓
I = ∫ 0 2 π 5 + 4 cos θ d θ
Forecast: Denominator 5 ± 4 = [ 1 , 9 ] ke beech ghoomta hai, hamesha positive, toh I > 0 . Kyunki 5 + 4 c o s θ 1 ka average 2 π length par 1 se kam hai, expect karo I < 2 π ≈ 6.28 . Kuch 2 π /3 ke aas-paas guess karo.
Step 1 — z = e i θ substitute karo. Kyun? Integral already ek loop hai; unit circle par yeh ek contour integral ban jaata hai jise hum residues se attack kar sakte hain. Use karo (parent ):
cos θ = 2 z + z − 1 , d θ = i z d z .
i z 1 kyun? Kyunki d z = i e i θ d θ = i z d θ . Ise bhoolna Cell-E ka classic trap hai.
Step 2 — rewrite karo.
I = ∮ ∣ z ∣ = 1 5 + 2 ( z + z − 1 ) 1 ⋅ i z d z = ∮ i z ( 5 + 2 z + 2 z − 1 ) d z = ∮ i ( 2 z 2 + 5 z + 2 ) d z .
z − 1 kyun vanish ho gaya? Numerator aur denominator ko z se multiply karne se woh clear ho gaya; note karo isse z = 0 par ek pole aa sakta hai — check karo! Yahan z = 0 pe 2 z 2 + 5 z + 2 = 2 = 0 , toh origin par koi extra pole nahi .
Step 3 — ∣ z ∣ = 1 ke andar poles. Solve karo 2 z 2 + 5 z + 2 = 0 : z = 4 − 5 ± 25 − 16 = 4 − 5 ± 3 , jisse milte hain z = − 2 1 aur z = − 2 .
Figure mein kya dekhna hai: yellow unit circle haara naya contour hai. z = − 2 1 par pink dot andar hai (distance 0.5 < 1 ), toh yeh count hota hai; z = − 2 par blue cross bahar hai (distance 2 > 1 ), toh isse discard karte hain. Sirf ek residue matter karta hai.
Step 4 — z 0 = − 2 1 par residue (simple). g = i 1 , h = 2 z 2 + 5 z + 2 , h ′ = 4 z + 5 ke saath:
Res − 1/2 = 4 ( − 2 1 ) + 5 1/ i = 3 1/ i = 3 i 1 .
Step 5 — assemble karo.
I = 2 π i ⋅ 3 i 1 = 3 2 π .
Verify: 3 2 π ≈ 2.094 , positive aur 2 π se kam — forecast se match. General formula check: ∫ 0 2 π a + b c o s θ d θ = a 2 − b 2 2 π se milta hai 25 − 16 2 π = 3 2 π . ✓
I = P.V. ∫ − ∞ ∞ x sin x d x
Forecast: Yeh famous "sinc" integral hai; integrand even hai, average-par-positive hai, aur answer π ka ek clean multiple hai. π guess karo.
Step 1 — degeneracy. Yeh special kyun hai? f ( z ) = z e i z ka pole z = 0 par hai, jo haare contour par andar nahi, par hi hai. Seedha recipe illegal hai; hum contour ko ek choti semicircle γ ε of radius ε ke saath indent karte hain jo origin ke upar se detour karti hai (dekho Principal Value Integrals ).
Figure mein kya dekhna hai: yellow path real axis ke saath chalti hai lekin origin ke aas-paas ek chota pink half-loop γ ε banati hai, phir bade upper arc Γ ρ se close hoti hai. Kyunki detour z = 0 ke around step karta hai, koi pole enclosed nahi — poora loop zero integrate ho jaata hai. Answer tiny arc se aayega, enclosed residue se nahi.
Step 2 — closed contour koi pole enclose nahi karta. Kyun? Sirf pole z = 0 tiny detour se exclude hai, aur koi strictly inside nahi hai. Toh Cauchy's Integral Theorem & Formula se:
∫ real, indented z e i z d z + ∫ γ ε z e i z d z + ∫ Γ ρ z e i z d z = 0.
Step 3 — bada arc marta hai (Jordan). R ( z ) = 1/ z → 0 , a = 1 > 0 ⟹ ∫ Γ ρ → 0 .
Step 4 — chota arc ek half residue contribute karta hai. Half kyun? Tiny semicircle simple pole z = 0 ke around π angle (half a full loop) clockwise sweep karta hai. Simple pole ke liye, ϕ angle ka chota arc i ϕ ⋅ Res contribute karta hai; yahan ϕ = − π (clockwise). Res 0 z e i z = e 0 = 1 , toh
∫ γ ε z e i z d z → − iπ ⋅ 1 = − iπ .
Step 5 — collect karo. Indented real integral principal value P.V. ∫ x e i x d x ki taraf tend karta hai:
P.V. ∫ − ∞ ∞ x e i x d x − iπ + 0 = 0 ⇒ P.V. ∫ x e i x d x = iπ .
Step 6 — imaginary part lo. sin x = Im ( e i x ) :
∫ − ∞ ∞ x s i n x d x = Im ( iπ ) = π .
(Real part deta hai P.V. ∫ x c o s x d x = 0 — sahi hai, woh integrand odd hai.)
Verify: π ≈ 3.1416 , forecast se match karta hai. Standard result: ∫ 0 ∞ x s i n x d x = 2 π , evenness se double karke = π . ✓
I = ∫ − ∞ ∞ x 2 + 1 x d x — kya machine legal hai?
Forecast: Integrand odd hai, toh principal value 0 hai, lekin ordinary integral diverge karta hai (large ∣ x ∣ par ∫ x d x jaisa). Naive residue recipe tumhe ek galat finite number degi — fraud pakdo.
Step 1 — pehle degree/decay check. Pehle kyun? Kyunki yeh guardrail hai. P = x toh deg P = 1 ; Q = x 2 + 1 toh deg Q = 2 ; gap = 2 − 1 = 1 . Koi oscillating e ia z nahi hai, toh Jordan apply nahi karta. Toh arc integral vanish nahi hota — boxed formula yahan illegal hai.
Step 2 — dikhao ki arc sach mein survive karta hai. Γ ρ par, z = ρ e i θ :
z 2 + 1 z ∼ ρ 1 , arc length = π ρ ⇒ bound ∼ ρ 1 ⋅ π ρ = π → 0.
Estimate ek constant par ruk jaata hai — arc contribute karta hai aur drop nahi ho sakta.
Step 3 — jo IS sach hai. Principal value oddness se exist karta hai:
P.V. ∫ − ∞ ∞ x 2 + 1 x d x = lim R → ∞ ∫ − R R x 2 + 1 x d x = lim R → ∞ [ 2 1 ln ( x 2 + 1 ) ] − R R = 0 ,
kyunki integrand odd hai aur limits symmetric hain. Lekin ordinary two-sided limit diverge karta hai.
Verify: Antiderivative 2 1 ln ( x 2 + 1 ) even hai, toh symmetrically evaluate karne par 0 deta hai; upper aur lower limits ko independently infinity par bhejo toh blow up ho jaata hai. Lesson: residues sum karne se pehle hamesha arc/decay check karo. ✓ (P.V. = 0 .)
Worked example Ek damped resonance ki signal energy
Ek electrical filter ka response hai H ( ω ) = ω 2 + ω 0 2 1 resonance parameter ω 0 > 0 ke saath. Total spectral weight hai I = ∫ − ∞ ∞ H ( ω ) d ω = ∫ − ∞ ∞ ω 2 + ω 0 2 d ω . Ise dhundho, aur ω 0 = 2 ke liye evaluate karo.
Forecast: Zyaada damping (bada ω 0 ) curve ko spread aur lower karna chahiye, toh I ko decrease karna chahiye jab ω 0 badhta hai — expect karo I ∝ 1/ ω 0 .
Step 1 — cell pehchano. Kyun? Rename karo ω → x : yeh Cell A hai, basic rational integral. Yahan P = 1 (deg P = 0 ) aur Q = ω 2 + ω 0 2 (deg Q = 2 ), toh gap = 2 ≥ 2 — machine legal hai. Poles ± i ω 0 par hain.
Step 2 — arc check. Gap = 2 . ✓
Step 3 — UHP pole & residue. z = i ω 0 (simple), h ′ = 2 z :
Res i ω 0 = 2 i ω 0 1 .
Step 4 — assemble karo.
I = 2 π i ⋅ 2 i ω 0 1 = ω 0 π .
Step 5 — ω 0 = 2 plug karo. I = 2 π .
Verify: I = π / ω 0 sach mein shrink karta hai jab ω 0 badhta hai — "more damping → lower weight" forecast se match. Units: π / ω 0 mein 1/ ω 0 scaling hai, ek broader, lower Lorentzian ke saath consistent. Numerically π /2 ≈ 1.571 , aur elementary check ∫ ω 2 + 4 d ω = 2 1 arctan 2 ω − ∞ ∞ = 2 1 ⋅ π = 2 π . ✓
Recall Humne matrix ki har cell cover ki
A basic rational ::: Ex 1 (π / 2 )
B higher-order pole ::: Ex 2 (π /16 )
C cosine via e ia z ::: Ex 3 (2 π e − 6 )
D odd → sine, imaginary part ::: Ex 4 (π / e )
E trig over [ 0 , 2 π ] ::: Ex 5 (2 π /3 )
F pole on the axis, principal value ::: Ex 6 (π )
G arc survives — recipe illegal ::: Ex 7 (P.V. = 0 )
H word problem ::: Ex 8 (π / ω 0 )
Mnemonic Chaar sawaalon mein decision tree
Oscillating hai? haan → e ia z + Jordan (Cells C, D, F). nahi → gap check karo.
Gap ≥ 2 ? haan → semicircle machine (A, B, H). nahi & oscillation bhi nahi → recipe illegal, P.V. socho (G).
Pole axis par hai? haan → indent karo, half-residue i ϕ Res (F).
Loop over [ 0 , 2 π ] ? → unit circle, yaad rakho i z d z (E).