4.7.3 · D3 · Maths › Partial Differential Equations › Fourier series — motivation from periodic functions
Yeh page Fourier series ki practice arena hai. Parent note ne recipe derive ki thi; yahan hum isse har tarah ke function par run karte hain jo tum ko diye ja sakte hain — even, odd, neither, constant, discontinuous, aisi period jo friendly 2 π nahi hai, aur ek real heat-flow word problem. Is page ke end tak koi bhi case class tumhe surprise nahi kar sakti.
Intuition "Har scenario" ka matlab kya hai
Ek Fourier problem char sawaalon se decide hoti hai. (1) Kya f even hai, odd hai, ya neither? (decide karta hai ki kaunse coefficients vanish hote hain). (2) Period 2 L kya hai? (decide karta hai ki n x likhte ho ya L nπ x ). (3) Kya f continuous hai ya koi jump hai? (decide karta hai ki series f ko hit karegi ya midpoint ko). (4) Kya yeh pure maths hai ya physics ka starting shape? (decide karta hai ki answer ka matlab kya hai). Hum saare combinations tabulate karte hain, phir har ek ko tackle karte hain.
Neeche ke har example ko same teen formulas mein daalte hain jo parent note se hain. Unhe ek baar yahan state karte hain taaki yeh page apne aap mein complete ho.
Intuition Symmetry shortcut jis par hum baar baar lean karte hain
Even and Odd Functions use karke: agar f even hai (f ( − x ) = f ( x ) ) toh saare b n = 0 (sirf cosines); agar f odd hai (f ( − x ) = − f ( x ) ) toh saare a n = 0 (sirf sines). Reason: symmetric interval [ − L , L ] par odd integrand ka integral 0 hota hai.
Decision-tree mein teen symmetry classes hain (even / odd / neither) jo do continuity classes (continuous / jump wali) aur period (2 π ya general 2 L ) se cross hoti hain. Neeche ke aath examples har reachable branch ko kam se kam ek baar cover karte hain; last column mein woh branch likhi hai.
Cell
Symmetry
Period 2 L
Continuity
Twist
Example
A
Odd
2 π (L = π )
Jump
classic odd+jump
Ex 1: square wave
B
Even
2 π (L = π )
Continuous (corner)
corner, no jump
Ex 2: ∣ x ∣
C
Neither
2 π (L = π )
Jump
dono a n aur b n nonzero
Ex 3: f ( x ) = x + 1 on ( − π , π )
D
Even
General 2 L = 2 π
Continuous
L nπ x rakho
Ex 4: x 2 on ( − 2 , 2 )
E
Constant
koi bhi
Continuous
degenerate: no waves
Ex 5: f = 3
F
Odd
2 π
Jump
jump par value (naye coeffs nahi)
Ex 6: square-wave series at x = 0
G
Half-range
sirf L diya
extension choose karo
exam twist
Ex 7: f = 1 ki sine vs cosine series
H
Physics
2 π
Jump (block)
word problem, units
Ex 8: initial heat profile
Intuition Cells A aur F ek hi symmetry/continuity type kyun share karte hain
Cells A aur F dono "odd + jump" hain, lekin woh alag alag lessons dete hain. Cell A (Ex 1) poochhta hai "coefficients kya hain?" — ek computation. Cell F (Ex 6) poochhta hai "finished series discontinuity par kya output deti hai?" — ek convergence question jiska jawab coefficients seedha nahi dete. Same function, do genuinely different skills, isliye unhe alag examples milte hain.
Note karo ki matrix poori decision-tree span karta hai: even+jump Ex 8 mein hai (ek block jo symmetric hai par discontinuous), neither+jump Ex 3 mein hai, aur even+continuous on a general period Ex 4 mein. Koi bhi reachable branch skip nahi hua.
Worked example Square wave
f ( x ) = { − 1 + 1 − π < x < 0 0 < x < π , period 2 π
Forecast: Padhne se pehle — a n , b n mein se kaun zero hoga, aur answer mein sab n honge ya sirf kuch?
Step 1. f ( − x ) = − f ( x ) , toh f odd hai. Hence a n = 0 for all n .
Yeh step kyun? Odd × cosine (even) odd hota hai, aur symmetric interval [ − π , π ] par odd function ka integral 0 hota hai. Isse shuru hone se pehle hi saare cosine integrals delete ho jaate hain.
Step 2. b n = π 1 ∫ − π π f ( x ) sin n x d x = π 2 ∫ 0 π sin n x d x .
Yeh step kyun? f sin n x hai odd× odd = even, 0 ke baare mein symmetric, toh [ − π , π ] par integral, [ 0 , π ] par integral ka do guna hai, jahan f = + 1 hai.
Step 3. ∫ 0 π sin n x d x = [ n − c o s n x ] 0 π = n 1 − cos nπ = n 1 − ( − 1 ) n .
Yeh step kyun? cos nπ = ( − 1 ) n : even n ke liye + 1 hai (zero deta hai) aur odd n ke liye − 1 hai (2/ n deta hai). Sirf odd harmonics bachte hain.
Step 4. Toh b n = π 2 ⋅ n 1 − ( − 1 ) n , yaani b 1 = π 4 , b 3 = 3 π 4 , …
f ( x ) = π 4 ( sin x + 3 s i n 3 x + 5 s i n 5 x + ⋯ )
Verify: b 1 = π 4 ≈ 1.2732 , b 2 = 0 , b 3 = 3 π 4 ≈ 0.4244 . Yeh exactly parent ka result hai — consistent.
The figure below carries the geometry of this convergence.
Intuition Upar ki figure mein kya observe karna hai
Navy step target square wave hai. Orange curve sirf π 4 sin x hai (ek term) — ek single hump jo pehle se hi up-then-down shape capture kar leta hai. Violet curve n = 3 , 5 add karta hai aur magenta curve aath odd harmonics use karta hai: har extra term plateaus ko ± 1 ki taraf flatten karta hai aur vertical edge ko steep karta hai. Do cheezein notice karo: (1) x = 0 par har sine zero hai, toh saare partial sums midpoint value 0 se guzarte hain — pink/violet curves origin se cross karte hain, wahan ± 1 tak kabhi nahi pahunchte; (2) jump ke bilkul paas woh chota spike jo overshoot karta hai, terms badhane par nahi sihrta — woh Gibbs Phenomenon hai, ek permanent ≈ 9% overshoot.
Worked example Triangle-like
f ( x ) = ∣ x ∣ on ( − π , π ) , period 2 π
Forecast: ∣ x ∣ even hai aur iska koi jump nahi hai (bas 0 par ek sharp corner). Guess: kaun se coefficients vanish honge, aur series fast ya slow converge karegi?
Step 1. f ( − x ) = ∣ − x ∣ = ∣ x ∣ = f ( x ) — even . Toh b n = 0 .
Yeh step kyun? Even× sine (odd) odd hota hai ⇒ integrate karke 0 deta hai. Sirf cosines bachte hain.
Step 2. a 0 = π 1 ∫ − π π ∣ x ∣ d x = π 2 ∫ 0 π x d x = π 2 ⋅ 2 π 2 = π .
Yeh step kyun? a 0 /2 average height hai. ∣ x ∣ ka average π /2 hai, aur wakai a 0 /2 = π /2 hai. Ek built-in sanity check.
Step 3. a n = π 2 ∫ 0 π x cos n x d x . Integration by parts karo:
∫ 0 π x cos n x d x = [ n x s i n n x ] 0 π − ∫ 0 π n s i n n x d x = 0 + [ n 2 c o s n x ] 0 π = n 2 ( − 1 ) n − 1 .
Yeh step kyun? Integration by parts sahi tool hai kyunki integrand polynomial × trig hai — x differentiate karne se polynomial khatam ho jaata hai, aur hum pure trig integral ki taraf aate hain. (Boundary term n x s i n n x dono ends par vanish karta hai kyunki sin nπ = 0 aur lower limit par x = 0 hai.)
Step 4. a n = π 2 ⋅ n 2 ( − 1 ) n − 1 : even n ke liye zero, aur odd n ke liye − π n 2 4 .
∣ x ∣ = 2 π − π 4 ( cos x + 9 c o s 3 x + 25 c o s 5 x + ⋯ )
Verify: a 1 = − π 4 ≈ − 1.2732 , a 2 = 0 , a 3 = − 9 π 4 ≈ − 0.1415 . Kyunki coefficients 1/ n 2 ki tarah shrink karte hain (jump square wave ke 1/ n ke comparison mein), yeh faster converge karta hai — corner, jump se milder hai. x = 0 plug karo: 0 = 2 π − π 4 ∑ odd n 2 1 , jisse ∑ odd n 2 1 = 8 π 2 milta hai, jo ek known identity hai.
f ( x ) = x + 1 on ( − π , π ) , period 2 π
Forecast: Constant 1 add karne se odd symmetry toot jaati hai: f ( − x ) = − x + 1 = ± f ( x ) . Toh yeh neither even hai na odd — aur ± π par jump bhi hai. Guess: ab kaun se coefficients bachenge?
Step 1. f = x + 1 ko split karo. Constant 1 sirf a 0 mein contribute karta hai; linear x odd hai aur sirf b n mein contribute karta hai.
Yeh step kyun? Fourier coefficients f mein linear hain, toh har piece compute karo aur add karo. Yeh x + 1 ko poora attack karne se kaafi aasan hai.
Step 2. Constant se: a 0 = π 1 ∫ − π π 1 d x = 2 , aur baaki saare a n = π 1 ∫ − π π 1 ⋅ cos n x d x = 0 .
Yeh step kyun? ∫ − π π cos n x d x = n s i n nπ ⋅ 2 -type terms vanish karte hain kyunki sin nπ = 0 ; ek constant ko koi non-constant cosine nahi chahiye.
Step 3. Linear part se (sawtooth work reuse karte hue): b n = π 1 ∫ − π π x sin n x d x = n 2 ( − 1 ) n + 1 .
Yeh step kyun? x odd hai, sin n x odd hai, product even hai; by parts ∫ − π π x sin n x d x = 2 ∫ 0 π x sin n x d x = n − 2 π ( − 1 ) n , aur π se divide karne par n 2 ( − 1 ) n + 1 milta hai.
Step 4 (the twist). Dono ek nonzero average (a 0 /2 = 1 ) aur sines ka poora set appear karta hai:
x + 1 = 1 + 2 ( sin x − 2 s i n 2 x + 3 s i n 3 x − ⋯ ) .
b n mein har n present hai, aur ek standalone constant bhi hai — "neither" function ki pehchaan.
Verify: a 0 = 2 (toh average = 1 ✓, 2 π 1 ∫ − π π ( x + 1 ) d x = 1 se match karta hai), b 1 = 2 , b 2 = − 1 , b 3 = 3 2 . x = π /2 par: 2 π + 1 = 1 + 2 ( 1 − 3 1 + 5 1 − ⋯ ) = 1 + 2 ⋅ 4 π ✓ (phir Leibniz). Jump x = π par series midpoint 2 1 [ f ( π − ) + f ( − π + )] = 2 1 [( π + 1 ) + ( − π + 1 )] = 1 par converge karti hai.
f ( x ) = x 2 on ( − 2 , 2 ) , period 4 (toh L = 2 )
Forecast: Period ab 4 hai, 2 π nahi. Teen formulas mein se kaun si break ho jaati hai agar tum bhool jaate ho aur 2 nπ x ki jagah n x likh dete ho?
Step 1. x 2 even hai ⇒ b n = 0 . L nπ x = 2 nπ x har jagah use karo.
Yeh step kyun? Wave ko period 2 L = 4 mein poore cycles complete karne chahiye; argument L nπ x yeh guarantee karta hai. n x likhne par galat period use hogi aur result bakwaas aayega.
Step 2. a 0 = L 1 ∫ − L L x 2 d x = 2 1 ∫ − 2 2 x 2 d x = 2 1 ⋅ 3 2 ⋅ 2 3 = 2 1 ⋅ 3 16 = 3 8 .
Toh average a 0 /2 = 3 4 hai. ([ − 2 , 2 ] par x 2 ka average 4 1 ∫ − 2 2 x 2 = 3 4 hai ✓.)
Step 3a (parity reduction). a n = L 1 ∫ − L L x 2 cos L nπ x d x . Integrand even× even = even hai, x = 0 ke baare mein symmetric:
a n = L 1 ∫ − L L x 2 cos L nπ x d x = L 2 ∫ 0 L x 2 cos L nπ x d x .
Yeh step kyun? Symmetric integral ko aadha fold karne se algebra aadhi ho jaati hai aur limits clean rehti hain (0 se shuru).
Step 3b (pehla integration by parts). k = L nπ lo. u = x 2 , d v = cos ( k x ) d x ke saath:
∫ 0 L x 2 cos ( k x ) d x = [ k x 2 s i n k x ] 0 L − ∫ 0 L k 2 x s i n k x d x .
Boundary term k L 2 s i n ( k L ) = k L 2 s i n ( nπ ) = 0 hai (kyunki k L = nπ aur sin nπ = 0 ).
Yeh step kyun? x 2 differentiate karne par 2 x banta hai — ek power khatam, ek aur baaki.
Step 3c (doosra integration by parts). Baaki ∫ 0 L 2 x sin ( k x ) d x ke liye, u = 2 x , d v = sin ( k x ) d x lo:
∫ 0 L 2 x sin ( k x ) d x = [ k − 2 x c o s k x ] 0 L + ∫ 0 L k 2 c o s k x d x = k − 2 L c o s ( nπ ) + k 2 2 s i n ( nπ ) = k − 2 L ( − 1 ) n .
Step 3b aur 3c ko saath rakhte hain (3b ka minus sign dhyan mein rakho):
∫ 0 L x 2 cos ( k x ) d x = 0 − k 1 ( k − 2 L ( − 1 ) n ) = k 2 2 L ( − 1 ) n = n 2 π 2 2 L ( − 1 ) n L 2 = n 2 π 2 2 L 3 ( − 1 ) n .
Yeh step kyun? 2 x differentiate ho kar constant ban jaata hai, toh doosra part-step kaam khatam karta hai.
Step 3d. Therefore
a n = L 2 ⋅ n 2 π 2 2 L 3 ( − 1 ) n = n 2 π 2 4 L 2 ( − 1 ) n = n 2 π 2 16 ( − 1 ) n ( L = 2 ) .
Step 4. x 2 = 3 4 + π 2 16 ∑ n = 1 ∞ n 2 ( − 1 ) n cos 2 nπ x .
Verify: a 1 = π 2 − 16 ≈ − 1.6211 , a 2 = 4 π 2 16 = π 2 4 ≈ 0.4053 . x = 2 par (jahan cos 2 nπ ⋅ 2 = cos nπ = ( − 1 ) n ) series deti hai 3 4 + π 2 16 ∑ n 2 ( − 1 ) n ( − 1 ) n = 3 4 + π 2 16 ⋅ 6 π 2 = 3 4 + 3 8 = 4 = 2 2 ✓ (∑ 1/ n 2 = π 2 /6 use kiya).
f ( x ) = 3 on ( − π , π )
Forecast: Ek flat line — "no wave" banane mein kitni waves lagti hain?
Step 1. a 0 = π 1 ∫ − π π 3 d x = π 1 ⋅ 6 π = 6 .
Yeh step kyun? a 0 /2 average value 3 ke barabar hona chahiye, aur 6/2 = 3 ✓.
Step 2. a n = π 1 ∫ − π π 3 cos n x d x = 0 (cosine ke poore cycles integrate ho kar 0 dete hain), aur b n = 0 bhi.
Yeh step kyun? Constant mein koi oscillation nahi, toh use har non-constant wave ki zero zaroorat hai.
Result: series = a 0 /2 = 3 . Verify: trivially 3 = 3 for all x . Yeh poori theory ka limiting/degenerate corner hai: akela "0 -frequency" term.
Worked example Square-wave series
x = 0 par kahan converge karti hai?
Forecast: x = 0 par square wave undefined hai (woh − 1 se + 1 ki taraf jump karta hai). π 4 ( sin x + ⋯ ) wahan actually kaunsa number output karta hai?
Step 1. π 4 ( sin x + 3 s i n 3 x + ⋯ ) mein x = 0 plug karo. Har sin ( n ⋅ 0 ) = 0 , toh sum 0 hai.
Yeh step kyun? Hum danger point par series ko numerically test karte hain.
Step 2. Convergence theorem kehta hai ki jump par series midpoint ke barabar hoti hai: 2 1 [ f ( 0 + ) + f ( 0 − )] = 2 1 [( + 1 ) + ( − 1 )] = 0 .
Yeh step kyun? Yeh koi coincidence nahi — series hamesha jump par difference split karti hai (jump ke paas overshoot ke liye Gibbs Phenomenon dekho).
Verify: theorem se predicted 0 = sum se computed 0 ✓. Lesson: Fourier series mein "= " discontinuities par midpoint equality hai, literal equality nahi.
Intuition Half-range expansion kya hoti hai, aur
π 2 kahan se aata hai
Kabhi kabhi f sirf ( 0 , π ) par diya jaata hai — poore 2 π -period ka aadha. Fourier machinery use karne ke liye pehle ( − π , 0 ) par missing half invent karni padti hai. Hum isko choose karne ke liye free hain, aur do choices special hain: oddly extend karo (mirror-and-flip, sirf sine series force karta hai) ya evenly (seedha mirror, sirf cosine series force karta hai). Dono taraf se period 2 π ban jaata hai, toh L = π .
Coefficient formula mein π 2 kyun aata hai: general b n = L 1 ∫ − L L f sin L nπ x d x se L = π ke saath shuru karo, milta hai b n = π 1 ∫ − π π f sin n x d x . Odd extension ke liye integrand f sin n x even hai (odd× odd), 0 ke baare mein symmetric, toh [ − π , π ] integral [ 0 , π ] integral ka do guna hai — aur π 1 ban jaata hai π 2 :
b n = π 2 ∫ 0 π f ( x ) sin n x d x .
Bilkul waise hi, even extension ke liye f cos n x even hai, toh a n = π 2 ∫ 0 π f ( x ) cos n x d x . π 2 koi nayi baat nahi — yeh standard formula hai jo symmetry se aadha fold hua hai.
f ( x ) = 1 on ( 0 , π ) ko (a) sine series, (b) cosine series mein expand karo
Forecast: Same function, same interval — phir bhi tum do alag valid Fourier series nikaal sakte ho. Kaise? Kyunki tum choose karte ho ki f ko ( − π , 0 ) par kaise extend karo.
Step 1 (sine series — odd extension). f ko oddly reflect karo: f ( − x ) = − 1 on ( − π , 0 ) . Toh a n = 0 aur, upar ke half-range formula se,
b n = π 2 ∫ 0 π 1 ⋅ sin n x d x = π 2 [ n − c o s n x ] 0 π = π 2 ⋅ n 1 − c o s nπ = π 2 ⋅ n 1 − ( − 1 ) n .
Yeh step kyun? Odd extension sine-only series force karti hai; yeh exactly wahi hai jo PDE demands karta hai jab f = 0 dono ends par ho (Dirichlet boundary).
Step 2 (sine coefficients evaluate karo). b n even n ke liye 0 hai aur odd n ke liye π n 4 hai, toh
1 = π 4 ( sin x + 3 s i n 3 x + 5 s i n 5 x + ⋯ ) ( 0 < x < π ) .
Yeh step kyun? Constant 1 ki odd extension hi ± 1 square wave hai, toh coefficients exactly Ex 1 se match karte hain — ek satisfying consistency check.
Step 3 (cosine series — even extension). f ko evenly reflect karo: f ( − x ) = + 1 . Toh b n = 0 , aur
a 0 = π 2 ∫ 0 π 1 d x = 2 , a n = π 2 ∫ 0 π cos n x d x = π 2 ⋅ n s i n nπ = 0 ( n ≥ 1 ) .
Yeh step kyun? Constant ka even extension sirf constant hi hota hai — flat line ko average ke alaawa koi cosine nahi chahiye.
Step 4 (cosine series likho).
1 = 2 a 0 = 1 ( 0 < x < π ) .
Dono series bilkul alag lagti hain (ek infinite sine sum vs. ek single constant) phir bhi ( 0 , π ) par same f represent karti hain — yahi exam twist hai.
Verify: sine version at x = π /2 : π 4 ( 1 − 3 1 + 5 1 − ⋯ ) = π 4 ⋅ 4 π = 1 ✓. Cosine version = 1 trivially ✓. Dono given half-interval par 1 ke barabar hain.
π m lambi ek rod ka temperature u ( x , 0 ) = 100 °C middle third ( π /3 , 2 π /3 ) mein aur baaki jagah 0 °C se shuru hota hai. Initial profile ke sine coefficients nikalo.
Forecast: Hot block middle mein symmetrically baithta hai, toh profile rod ke centre ke baare mein symmetric hai phir bhi discontinuous hai (block edges par jump karti hai) — "even+jump" branch. Heat Equation par Separation of Variables modes sin ( n x ) deta hai; humein yeh blocky initial temperature unme expand karni hai. Kaun se coefficients dominate karenge?
Step 1. Model: u ( x , 0 ) = { 100 0 π /3 < x < 2 π /3 else on ( 0 , π ) . Sine (odd, half-range) expansion use karo kyunki rod ke ends 0 °C par held hain, yaani u ( 0 ) = u ( π ) = 0 .
Sine kyun? Boundary conditions cosines ko kill karte hain; sirf sin ( n x ) satisfy karta hai u ( 0 ) = u ( π ) = 0 . Ex 7 se half-range formula reuse karte hain: b n = π 2 ∫ 0 π u sin n x d x .
Step 2. b n = π 2 ∫ π /3 2 π /3 100 sin n x d x = π 200 [ n − cos n x ] π /3 2 π /3 = nπ 200 ( cos 3 nπ − cos 3 2 nπ ) .
Yeh step kyun? Sirf wahan integrate karo jahan u = 100 hai; baaki jagah integrand 0 hai.
Step 3. Amplitudes compute karo (°C): n = 1 ke saath, cos 3 π − cos 3 2 π = 2 1 − ( − 2 1 ) = 1 , toh b 1 = π 200 ≈ 63.66 . n = 2 ke saath, cos 3 2 π − cos 3 4 π = − 2 1 − ( − 2 1 ) = 0 , toh b 2 = 0 . n = 3 ke saath, cos π − cos 2 π = − 1 − 1 = − 2 , toh b 3 = 3 π 200 ⋅ ( − 2 ) = 3 π − 400 ≈ − 42.44 .
Units check kyun? b n same unit carry karta hai jitna u , yaani °C — yeh us temperature mode ki amplitude hai.
Step 4. Initial profile hai u ( x , 0 ) = ∑ n = 1 ∞ b n sin n x , aur jaise time chalti hai Heat Equation mode n ko e − n 2 t se multiply karta hai, toh high-n modes sabse jaldi decay karte hain — sharp block jaldi smooth ho jaata hai. (Wave Equation se compare karo, jahan har mode time mein oscillate karta hai.)
Verify: b 1 ≈ 63.66 °C, b 2 = 0 , b 3 ≈ − 42.44 °C. Dominant b 1 hump sabse slow-decaying hai, toh late-time temperature ek single sin x arch jaisi lagti hai — physically sensible.
Recall Kisi bhi Fourier problem se pehle scenario checklist
Symmetry? ::: Even ⇒ sirf a n ; odd ⇒ sirf b n ; neither ⇒ dono.
Period 2 L = ? ::: Hamesha L nπ x likho; sirf tab n x simplify karo jab L = π ho.
Jump present hai? ::: Series wahan midpoint par converge karti hai (Gibbs overshoot paas mein).
Half-range diya hai? ::: Tum even (cosine) ya odd (sine) extension choose karo boundary conditions ke hisaab se; π 1 symmetry se π 2 ban jaata hai.
Parent recipe & derivation: Fourier series motivation
Coefficients isolate kyun hote hain: Orthogonality of Functions
Symmetry shortcuts: Even and Odd Functions
Yeh series kahan se aati hain: Separation of Variables , Heat Equation aur Wave Equation par apply ki gayi
Jumps par overshoot: Gibbs Phenomenon
Continuous-spectrum sibling: Fourier Transform